def test_to_honeybee_multiple_models():
"""Test to_honeybee with multiple honeybee models."""
pts_1 = (Point3D(0, 0, 3), Point3D(10, 0, 3), Point3D(10, 10, 3), Point3D(0, 10, 3))
pts_2 = (Point3D(10, 0, 3), Point3D(20, 0, 3), Point3D(20, 10, 3), Point3D(10, 10, 3))
pts_3 = (Point3D(0, 20, 3), Point3D(20, 20, 3), Point3D(20, 30, 3), Point3D(0, 30, 3))
room2d_1 = Room2D('Office1', Face3D(pts_1), 3)
room2d_2 = Room2D('Office2', Face3D(pts_2), 3)
room2d_3 = Room2D('Office3', Face3D(pts_3), 3)
story_big = Story('OfficeFloorBig', [room2d_3])
story = Story('OfficeFloor', [room2d_1, room2d_2])
story.solve_room_2d_adjacency(0.01)
story.set_outdoor_window_parameters(SimpleWindowRatio(0.4))
story.multiplier = 4
building = Building('OfficeBuilding', [story])
story_big.set_outdoor_window_parameters(SimpleWindowRatio(0.4))
story_big.multiplier = 4
building_big = Building('OfficeBuildingBig', [story_big])
tree_canopy_geo1 = Face3D.from_regular_polygon(6, 6, Plane(o=Point3D(5, -10, 6)))
tree_canopy_geo2 = Face3D.from_regular_polygon(6, 2, Plane(o=Point3D(-5, -10, 3)))
tree_canopy = ContextShade('TreeCanopy', [tree_canopy_geo1, tree_canopy_geo2])
model = Model('NewDevelopment', [building, building_big], [tree_canopy])
hb_models = model.to_honeybee('Building', 10, False, tolerance=0.01)
assert len(hb_models) == 2
assert isinstance(hb_models[0], hb_model.Model)
assert isinstance(hb_models[-1], hb_model.Model)
assert len(hb_models[-1].rooms) == 4
assert len(hb_models[-1].rooms[-1]) == 6
assert hb_models[-1].rooms[-1].volume == 600
assert hb_models[-1].rooms[-1].floor_area == 200
assert hb_models[-1].rooms[-1].exterior_wall_area == 180
assert len(hb_models[0].orphaned_shades) == 6
assert len(hb_models[-1].orphaned_shades) == 6
def test_properties():
"""Test various properties on the model."""
pts_1 = (Point3D(0, 0, 3), Point3D(10, 0, 3), Point3D(10, 10, 3), Point3D(0, 10, 3))
pts_2 = (Point3D(10, 0, 3), Point3D(20, 0, 3), Point3D(20, 10, 3), Point3D(10, 10, 3))
pts_3 = (Point3D(0, 20, 3), Point3D(20, 20, 3), Point3D(20, 30, 3), Point3D(0, 30, 3))
room2d_1 = Room2D('Office1', Face3D(pts_1), 3)
room2d_2 = Room2D('Office2', Face3D(pts_2), 3)
room2d_3 = Room2D('Office3', Face3D(pts_3), 3)
story_big = Story('OfficeFloorBig', [room2d_3])
story = Story('OfficeFloor', [room2d_1, room2d_2])
story.solve_room_2d_adjacency(0.01)
story.set_outdoor_window_parameters(SimpleWindowRatio(0.4))
story.multiplier = 4
building = Building('OfficeBuilding', [story])
story_big.set_outdoor_window_parameters(SimpleWindowRatio(0.4))
story_big.multiplier = 2
building_big = Building('OfficeBuildingBig', [story_big])
model = Model('NewDevelopment', [building, building_big])
assert model.average_story_count == 3
assert model.average_story_count_above_ground == 3
assert model.average_height == 12
assert model.average_height_above_ground == 9
assert model.footprint_area == 100 * 4
assert model.floor_area == (100 * 2 * 4) + (200 * 2)
assert model.exterior_wall_area == (90 * 2 * 4) + (180 * 2)
assert model.exterior_aperture_area == (90 * 2 * 4 * 0.4) + (180 * 2 * 0.4)
assert model.volume == (100 * 2 * 4 * 3) + (200 * 2 * 3)
def two_buildings():
pts_1 = (Point3D(0, 0, 3), Point3D(10, 0,
3), Point3D(10, 10,
3), Point3D(0, 10, 3))
pts_2 = (Point3D(10, 0, 3), Point3D(20, 0,
3), Point3D(20, 10,
3), Point3D(10, 10, 3))
pts_3 = (Point3D(0, 20, 3), Point3D(20, 20,
3), Point3D(20, 30,
3), Point3D(0, 30, 3))
room2d_1 = Room2D('Office1', Face3D(pts_1), 3)
room2d_2 = Room2D('Office2', Face3D(pts_2), 3)
room2d_3 = Room2D('Office3', Face3D(pts_3), 3)
story_big = Story('OfficeFloorBig', [room2d_3])
story = Story('OfficeFloor', [room2d_1, room2d_2])
story.solve_room_2d_adjacency(0.01)
story.set_outdoor_window_parameters(SimpleWindowRatio(0.3))
story.multiplier = 3
building = Building('OfficeBuilding', [story])
building.separate_top_bottom_floors()
story_big.set_outdoor_window_parameters(SimpleWindowRatio(0.6))
story_big.multiplier = 4
building_big = Building('OfficeBuildingBig', [story_big])
building_big.separate_top_bottom_floors()
tree_canopy_geo1 = Face3D.from_regular_polygon(6, 6,
Plane(o=Point3D(5, -10, 6)))
tree_canopy_geo2 = Face3D.from_regular_polygon(
6, 2, Plane(o=Point3D(-5, -10, 3)))
tree_canopy = ContextShade('TreeCanopy',
[tree_canopy_geo1, tree_canopy_geo2])
model = Model('NewDevelopment', [building, building_big], [tree_canopy])
return model
def test_to_geojson():
"""Test the Model to_geojson method."""
pts_1 = (Point3D(50, 50, 3), Point3D(60, 50, 3), Point3D(60, 60, 3), Point3D(50, 60, 3))
pts_2 = (Point3D(60, 50, 3), Point3D(70, 50, 3), Point3D(70, 60, 3), Point3D(60, 60, 3))
pts_3 = (Point3D(50, 70, 3), Point3D(70, 70, 3), Point3D(70, 80, 3), Point3D(50, 80, 3))
room2d_1 = Room2D('Residence1', Face3D(pts_1), 3)
room2d_2 = Room2D('Residence2', Face3D(pts_2), 3)
room2d_3 = Room2D('Retail', Face3D(pts_3), 3)
story_big = Story('RetailFloor', [room2d_3])
story = Story('ResidenceFloor', [room2d_1, room2d_2])
story.solve_room_2d_adjacency(0.01)
story.set_outdoor_window_parameters(SimpleWindowRatio(0.4))
story.multiplier = 3
building = Building('ResidenceBuilding', [story])
story_big.set_outdoor_window_parameters(SimpleWindowRatio(0.4))
story_big.multiplier = 1
building_big = Building('RetailBuildingBig', [story_big])
pts_1 = (Point3D(0, 0, 3), Point3D(0, 5, 3), Point3D(15, 5, 3), Point3D(15, 0, 3))
pts_2 = (Point3D(15, 0, 3), Point3D(15, 15, 3), Point3D(20, 15, 3), Point3D(20, 0, 3))
pts_3 = (Point3D(0, 5, 3), Point3D(0, 20, 3), Point3D(5, 20, 3), Point3D(5, 5, 3))
pts_4 = (Point3D(5, 15, 3), Point3D(5, 20, 3), Point3D(20, 20, 3), Point3D(20, 15, 3))
pts_5 = (Point3D(-5, -5, 3), Point3D(-10, -5, 3), Point3D(-10, -10, 3), Point3D(-5, -10, 3))
room2d_1 = Room2D('Office1', Face3D(pts_1), 3)
room2d_2 = Room2D('Office2', Face3D(pts_2), 3)
room2d_3 = Room2D('Office3', Face3D(pts_3), 3)
room2d_4 = Room2D('Office4', Face3D(pts_4), 3)
room2d_5 = Room2D('Office5', Face3D(pts_5), 3)
int_rms = Room2D.intersect_adjacency(
[room2d_1, room2d_2, room2d_3, room2d_4, room2d_5], 0.01)
story = Story('OfficeFloor', int_rms)
story.rotate_xy(5, Point3D(0, 0, 0))
story.solve_room_2d_adjacency(0.01)
story.set_outdoor_window_parameters(SimpleWindowRatio(0.4))
story.multiplier = 5
building_mult = Building('OfficeBuilding', [story])
tree_canopy_geo1 = Face3D.from_regular_polygon(6, 6, Plane(o=Point3D(5, -10, 6)))
tree_canopy_geo2 = Face3D.from_regular_polygon(6, 2, Plane(o=Point3D(-5, -10, 3)))
tree_canopy = ContextShade('TreeCanopy', [tree_canopy_geo1, tree_canopy_geo2])
model = Model('TestGeoJSON', [building, building_big, building_mult], [tree_canopy])
location = Location('Boston', 'MA', 'USA', 42.366151, -71.019357)
geojson_folder = './tests/geojson/'
model.to_geojson(location, folder=geojson_folder)
geo_fp = os.path.join(
geojson_folder, model.identifier, '{}.geojson'.format(model.identifier))
assert os.path.isfile(geo_fp)
nukedir(os.path.join(geojson_folder, model.identifier), True)
def test_buildings_to_honeybee():
"""Test the buildings_to_honeybee method."""
pts_1 = (Point3D(0, 0, 3), Point3D(10, 0,
3), Point3D(10, 10,
3), Point3D(0, 10, 3))
pts_2 = (Point3D(10, 0, 3), Point3D(20, 0,
3), Point3D(20, 10,
3), Point3D(10, 10, 3))
pts_3 = (Point3D(0, 20, 3), Point3D(20, 20,
3), Point3D(20, 30,
3), Point3D(0, 30, 3))
room2d_1 = Room2D('Office1', Face3D(pts_1), 3)
room2d_2 = Room2D('Office2', Face3D(pts_2), 3)
room2d_3 = Room2D('Office3', Face3D(pts_3), 3)
story_big = Story('Office_Floor_Big', [room2d_3])
story = Story('Office_Floor', [room2d_1, room2d_2])
story.solve_room_2d_adjacency(0.01)
story.set_outdoor_window_parameters(SimpleWindowRatio(0.4))
story.multiplier = 4
building = Building('Office_Building', [story])
story_big.set_outdoor_window_parameters(SimpleWindowRatio(0.4))
story_big.multiplier = 4
building_big = Building('Office_Building_Big', [story_big])
hb_models = Building.buildings_to_honeybee([building, building_big], None,
None, False)
assert len(hb_models) == 2
assert isinstance(hb_models[0], Model)
assert len(hb_models[0].orphaned_shades) == 4
assert isinstance(hb_models[-1], Model)
assert len(hb_models[-1].rooms) == 4
assert len(hb_models[-1].rooms[-1]) == 6
assert hb_models[-1].rooms[-1].volume == 600
assert hb_models[-1].rooms[-1].floor_area == 200
assert hb_models[-1].rooms[-1].exterior_wall_area == 180
hb_models = Building.buildings_to_honeybee([building, building_big], None,
5, True)
assert len(hb_models) == 2
assert isinstance(hb_models[0], Model)
assert len(hb_models[0].orphaned_shades) == 0
assert len(hb_models[0].rooms) == 2
for room in hb_models[0].rooms:
assert room.multiplier == 4
hb_models = Building.buildings_to_honeybee([building, building_big],
cap=True)
assert len(hb_models) == 2
assert isinstance(hb_models[0], Model)
assert len(hb_models[0].orphaned_shades) == 5
def test_to_dict():
"""Test the Building to_dict method."""
pts_1 = (Point3D(0, 0, 2), Point3D(10, 0,
2), Point3D(10, 10,
2), Point3D(0, 10, 2))
pts_2 = (Point3D(10, 0, 3), Point3D(20, 0,
3), Point3D(20, 10,
3), Point3D(10, 10, 3))
room2d_1 = Room2D('Office1', Face3D(pts_1), 5)
room2d_2 = Room2D('Office2', Face3D(pts_2), 3)
story = Story('Office_Floor', [room2d_1, room2d_2])
story.solve_room_2d_adjacency(0.01)
story.set_outdoor_window_parameters(SimpleWindowRatio(0.4))
story.set_outdoor_shading_parameters(Overhang(1))
story.multiplier = 4
building = Building('Office_Building', [story])
building.separate_top_bottom_floors()
bd = building.to_dict()
assert bd['type'] == 'Building'
assert bd['identifier'] == 'Office_Building'
assert bd['display_name'] == 'Office_Building'
assert 'unique_stories' in bd
assert len(bd['unique_stories']) == 3
assert 'properties' in bd
assert bd['properties']['type'] == 'BuildingProperties'
def test_reflect():
"""Test the Building reflect method."""
pts = (Point3D(1, 1, 2), Point3D(2, 1, 2), Point3D(2, 2,
2), Point3D(1, 2, 2))
plane = Plane(Vector3D(0, 0, 1), Point3D(0, 0, 2))
room = Room2D('Square_Shoebox', Face3D(pts, plane), 3)
story = Story('Office_Floor', [room])
story.multiplier = 4
building = Building('Office_Building', [story])
origin_1 = Point3D(1, 0, 2)
normal_1 = Vector3D(1, 0, 0)
plane_1 = Plane(normal_1, origin_1)
test_1 = building.duplicate()
test_1.reflect(plane_1)
assert test_1.unique_stories[0].room_2ds[0].floor_geometry[
-1].x == pytest.approx(1, rel=1e-3)
assert test_1.unique_stories[0].room_2ds[0].floor_geometry[
-1].y == pytest.approx(1, rel=1e-3)
assert test_1.unique_stories[0].room_2ds[0].floor_geometry[
-1].z == pytest.approx(2, rel=1e-3)
assert test_1.unique_stories[0].room_2ds[0].floor_geometry[
1].x == pytest.approx(0, rel=1e-3)
assert test_1.unique_stories[0].room_2ds[0].floor_geometry[
1].y == pytest.approx(2, rel=1e-3)
assert test_1.unique_stories[0].room_2ds[0].floor_geometry[
1].z == pytest.approx(2, rel=1e-3)
def test_building_shade_representation():
"""Test the Building shade_representation method."""
pts_1 = (Point3D(0, 0, 3), Point3D(0, 10,
3), Point3D(10, 10,
3), Point3D(10, 0, 3))
pts_2 = (Point3D(10, 0, 3), Point3D(10, 10,
3), Point3D(20, 10,
3), Point3D(20, 0, 3))
pts_3 = (Point3D(0, 10, 3), Point3D(0, 20,
3), Point3D(10, 20,
3), Point3D(10, 10, 3))
pts_4 = (Point3D(10, 10, 3), Point3D(10, 20,
3), Point3D(20, 20,
3), Point3D(20, 10, 3))
room2d_1 = Room2D('Office1', Face3D(pts_1), 3)
room2d_2 = Room2D('Office2', Face3D(pts_2), 3)
room2d_3 = Room2D('Office3', Face3D(pts_3), 3)
room2d_4 = Room2D('Office4', Face3D(pts_4), 3)
story = Story('Office_Floor', [room2d_1, room2d_2, room2d_3, room2d_4])
story.solve_room_2d_adjacency(0.01)
story.multiplier = 4
building = Building('Office_Building', [story])
shade_rep = building.shade_representation(tolerance=0.01)
assert len(shade_rep) == 8
shd_area = sum([shd.area for shd in shade_rep])
assert shd_area == building.exterior_wall_area
def test_building_footprint_courtyard():
"""Test the building footprint method with courtyard geometry."""
pts_1 = (Point3D(0, 0, 3), Point3D(0, 5, 3), Point3D(15, 5,
3), Point3D(15, 0, 3))
pts_2 = (Point3D(15, 0, 3), Point3D(15, 15,
3), Point3D(20, 15,
3), Point3D(20, 0, 3))
pts_3 = (Point3D(0, 5, 3), Point3D(0, 20, 3), Point3D(5, 20,
3), Point3D(5, 5, 3))
pts_4 = (Point3D(5, 15, 3), Point3D(5, 20,
3), Point3D(20, 20,
3), Point3D(20, 15, 3))
room2d_1 = Room2D('Office1', Face3D(pts_1), 3)
room2d_2 = Room2D('Office2', Face3D(pts_2), 3)
room2d_3 = Room2D('Office3', Face3D(pts_3), 3)
room2d_4 = Room2D('Office4', Face3D(pts_4), 3)
int_rms = Room2D.intersect_adjacency(
[room2d_1, room2d_2, room2d_3, room2d_4], 0.01)
story = Story('Office_Floor', int_rms)
story.rotate_xy(5, Point3D(0, 0, 0))
story.solve_room_2d_adjacency(0.01)
story.set_outdoor_window_parameters(SimpleWindowRatio(0.4))
story.multiplier = 4
building = Building('Office_Building', [story])
footprint = building.footprint(0.01)
assert len(footprint) == 1
assert isinstance(footprint[0], Face3D)
assert len(footprint[0].boundary) == 8
assert len(footprint[0].holes) == 1
assert len(footprint[0].holes[0]) == 4
def test_to_dict():
"""Test the Building to_dict method with energy properties."""
mass_set = ConstructionSet('Thermal Mass Construction Set')
pts_1 = (Point3D(0, 0, 3), Point3D(0, 10, 3), Point3D(10, 10, 3), Point3D(10, 0, 3))
pts_2 = (Point3D(10, 0, 3), Point3D(10, 10, 3), Point3D(20, 10, 3), Point3D(20, 0, 3))
pts_3 = (Point3D(0, 10, 3), Point3D(0, 20, 3), Point3D(10, 20, 3), Point3D(10, 10, 3))
pts_4 = (Point3D(10, 10, 3), Point3D(10, 20, 3), Point3D(20, 20, 3), Point3D(20, 10, 3))
room2d_1 = Room2D('Office1', Face3D(pts_1), 3)
room2d_2 = Room2D('Office2', Face3D(pts_2), 3)
room2d_3 = Room2D('Office3', Face3D(pts_3), 3)
room2d_4 = Room2D('Office4', Face3D(pts_4), 3)
story = Story('OfficeFloor', [room2d_1, room2d_2, room2d_3, room2d_4])
story.solve_room_2d_adjacency(0.01)
story.set_outdoor_window_parameters(SimpleWindowRatio(0.4))
story.multiplier = 4
building = Building('OfficeBuilding', [story])
bd = building.to_dict()
assert 'properties' in bd
assert bd['properties']['type'] == 'BuildingProperties'
assert 'energy' in bd['properties']
assert bd['properties']['energy']['type'] == 'BuildingEnergyProperties'
assert 'construction_set' not in bd['properties']['energy'] or \
bd['properties']['energy']['construction_set'] is None
building.properties.energy.construction_set = mass_set
bd = building.to_dict()
assert bd['properties']['energy']['construction_set'] is not None
def test_building_footprint_simple():
"""Test the building footprint method with simple geometry."""
pts_1 = (Point3D(0, 0, 3), Point3D(0, 10,
3), Point3D(10, 10,
3), Point3D(10, 0, 3))
pts_2 = (Point3D(10, 0, 3), Point3D(10, 10,
3), Point3D(20, 10,
3), Point3D(20, 0, 3))
pts_3 = (Point3D(0, 10, 3), Point3D(0, 20,
3), Point3D(10, 20,
3), Point3D(10, 10, 3))
pts_4 = (Point3D(10, 10, 3), Point3D(10, 20,
3), Point3D(20, 20,
3), Point3D(20, 10, 3))
room2d_1 = Room2D('Office1', Face3D(pts_1), 3)
room2d_2 = Room2D('Office2', Face3D(pts_2), 3)
room2d_3 = Room2D('Office3', Face3D(pts_3), 3)
room2d_4 = Room2D('Office4', Face3D(pts_4), 3)
story = Story('Office_Floor', [room2d_1, room2d_2, room2d_3, room2d_4])
story.solve_room_2d_adjacency(0.01)
story.set_outdoor_window_parameters(SimpleWindowRatio(0.4))
story.multiplier = 4
building = Building('Office_Building', [story])
footprint = building.footprint(0.01)
assert len(footprint) == 1
assert isinstance(footprint[0], Face3D)
assert footprint[0].holes is None
assert len(footprint[0].vertices) == 8
assert footprint[0].min == Point3D(0, 0, 3)
assert footprint[0].max == Point3D(20, 20, 3)
def test_reflect():
"""Test the Model reflect method."""
pts = (Point3D(1, 1, 2), Point3D(2, 1, 2), Point3D(2, 2, 2), Point3D(1, 2, 2))
plane = Plane(Vector3D(0, 0, 1), Point3D(0, 0, 2))
room = Room2D('SquareShoebox', Face3D(pts, plane), 3)
story = Story('OfficeFloor', [room])
story.multiplier = 4
building = Building('OfficeBuilding', [story])
pts = (Point3D(1, 1, 2), Point3D(2, 1, 2), Point3D(2, 2, 2), Point3D(1, 2, 2))
plane = Plane(Vector3D(0, 0, 1), Point3D(0, 0, 2))
awning_canopy = ContextShade('AwningCanopy', [Face3D(pts, plane)])
model = Model('NewDevelopment', [building], [awning_canopy])
origin_1 = Point3D(1, 0, 2)
normal_1 = Vector3D(1, 0, 0)
plane_1 = Plane(normal_1, origin_1)
test_1 = model.duplicate()
test_1.reflect(plane_1)
assert test_1.buildings[0].unique_stories[0].room_2ds[0].floor_geometry[-1].x == pytest.approx(1, rel=1e-3)
assert test_1.buildings[0].unique_stories[0].room_2ds[0].floor_geometry[-1].y == pytest.approx(1, rel=1e-3)
assert test_1.buildings[0].unique_stories[0].room_2ds[0].floor_geometry[-1].z == pytest.approx(2, rel=1e-3)
assert test_1.buildings[0].unique_stories[0].room_2ds[0].floor_geometry[1].x == pytest.approx(0, rel=1e-3)
assert test_1.buildings[0].unique_stories[0].room_2ds[0].floor_geometry[1].y == pytest.approx(2, rel=1e-3)
assert test_1.buildings[0].unique_stories[0].room_2ds[0].floor_geometry[1].z == pytest.approx(2, rel=1e-3)
assert test_1.context_shades[0][0][-1].x == pytest.approx(1, rel=1e-3)
assert test_1.context_shades[0][0][-1].y == pytest.approx(1, rel=1e-3)
assert test_1.context_shades[0][0][-1].z == pytest.approx(2, rel=1e-3)
assert test_1.context_shades[0][0][1].x == pytest.approx(0, rel=1e-3)
assert test_1.context_shades[0][0][1].y == pytest.approx(2, rel=1e-3)
assert test_1.context_shades[0][0][1].z == pytest.approx(2, rel=1e-3)
def test_set_all_room_2d_program_type():
"""Test the set_all_room_2d_program_type method on a Building."""
pts_1 = (Point3D(0, 0, 3), Point3D(0, 10, 3), Point3D(10, 10, 3), Point3D(10, 0, 3))
pts_2 = (Point3D(10, 0, 3), Point3D(10, 10, 3), Point3D(20, 10, 3), Point3D(20, 0, 3))
pts_3 = (Point3D(0, 10, 3), Point3D(0, 20, 3), Point3D(10, 20, 3), Point3D(10, 10, 3))
pts_4 = (Point3D(10, 10, 3), Point3D(10, 20, 3), Point3D(20, 20, 3), Point3D(20, 10, 3))
room2d_1 = Room2D('Office1', Face3D(pts_1), 3)
room2d_2 = Room2D('Office2', Face3D(pts_2), 3)
room2d_3 = Room2D('Office3', Face3D(pts_3), 3)
room2d_4 = Room2D('Office4', Face3D(pts_4), 3)
story = Story('OfficeFloor', [room2d_1, room2d_2, room2d_3, room2d_4])
story.solve_room_2d_adjacency(0.01)
story.set_outdoor_window_parameters(SimpleWindowRatio(0.4))
story.multiplier = 4
building = Building('OfficeBuilding', [story])
lab_program = office_program.duplicate()
lab_program.identifier = 'Bio Laboratory'
lab_program.electric_equipment.watts_per_area = 50
lab_program.ventilation.flow_per_person = 0
lab_program.ventilation.flow_per_area = 0
lab_program.ventilation.air_changes_per_hour = 6
building.properties.energy.set_all_room_2d_program_type(lab_program)
assert all(room_2d.properties.energy.program_type == lab_program
for room_2d in building.unique_room_2ds)
def test_to_from_dfpkl_methods():
"""Test the to/from dfpkl methods."""
pts_1 = (Point3D(0, 0, 3), Point3D(0, 10, 3), Point3D(10, 10, 3), Point3D(10, 0, 3))
pts_2 = (Point3D(10, 0, 3), Point3D(10, 10, 3), Point3D(20, 10, 3), Point3D(20, 0, 3))
pts_3 = (Point3D(0, 10, 3), Point3D(0, 20, 3), Point3D(10, 20, 3), Point3D(10, 10, 3))
pts_4 = (Point3D(10, 10, 3), Point3D(10, 20, 3), Point3D(20, 20, 3), Point3D(20, 10, 3))
room2d_1 = Room2D('Office1', Face3D(pts_1), 3)
room2d_2 = Room2D('Office2', Face3D(pts_2), 3)
room2d_3 = Room2D('Office3', Face3D(pts_3), 3)
room2d_4 = Room2D('Office4', Face3D(pts_4), 3)
story = Story('OfficeFloor', [room2d_1, room2d_2, room2d_3, room2d_4])
story.solve_room_2d_adjacency(0.01)
story.set_outdoor_window_parameters(SimpleWindowRatio(0.4))
story.multiplier = 4
building = Building('OfficeBuilding', [story])
tree_canopy_geo1 = Face3D.from_regular_polygon(6, 6, Plane(o=Point3D(5, -10, 6)))
tree_canopy_geo2 = Face3D.from_regular_polygon(6, 2, Plane(o=Point3D(-5, -10, 3)))
tree_canopy = ContextShade('TreeCanopy', [tree_canopy_geo1, tree_canopy_geo2])
model = Model('NewDevelopment', [building], [tree_canopy])
model_dfpkl = model.to_dfpkl('test')
assert os.path.isfile(model_dfpkl)
new_model = Model.from_dfpkl(model_dfpkl)
assert isinstance(new_model, Model)
os.remove(model_dfpkl)
def test_model_properties_setability():
"""Test the setting of properties on the Model."""
pts_1 = (Point3D(0, 0, 3), Point3D(0, 10, 3), Point3D(10, 10, 3), Point3D(10, 0, 3))
pts_2 = (Point3D(10, 0, 3), Point3D(10, 10, 3), Point3D(20, 10, 3), Point3D(20, 0, 3))
pts_3 = (Point3D(0, 10, 3), Point3D(0, 20, 3), Point3D(10, 20, 3), Point3D(10, 10, 3))
pts_4 = (Point3D(10, 10, 3), Point3D(10, 20, 3), Point3D(20, 20, 3), Point3D(20, 10, 3))
room2d_1 = Room2D('Office1', Face3D(pts_1), 3)
room2d_2 = Room2D('Office2', Face3D(pts_2), 3)
room2d_3 = Room2D('Office3', Face3D(pts_3), 3)
room2d_4 = Room2D('Office4', Face3D(pts_4), 3)
story = Story('OfficeFloor', [room2d_1, room2d_2, room2d_3, room2d_4])
story.solve_room_2d_adjacency(0.01)
story.set_outdoor_window_parameters(SimpleWindowRatio(0.4))
story.multiplier = 4
building = Building('OfficeBuilding', [story])
model = Model('NewDevelopment', [building])
model.to_rectangular_windows()
model.display_name = 'TestBuilding'
assert model.display_name == 'TestBuilding'
model.units = 'Feet'
assert model.units == 'Feet'
model.tolerance = 0.1
assert model.tolerance == 0.1
model.angle_tolerance = 0.01
assert model.angle_tolerance == 0.01
model.tolerance = 0.01
assert model.tolerance == 0.01
def test_set_all_room_2d_hvac():
"""Test the set_all_room_2d_hvac method on a Building."""
pts_1 = (Point3D(0, 0, 3), Point3D(0, 10, 3), Point3D(10, 10, 3), Point3D(10, 0, 3))
pts_2 = (Point3D(10, 0, 3), Point3D(10, 10, 3), Point3D(20, 10, 3), Point3D(20, 0, 3))
pts_3 = (Point3D(0, 10, 3), Point3D(0, 20, 3), Point3D(10, 20, 3), Point3D(10, 10, 3))
pts_4 = (Point3D(10, 10, 3), Point3D(10, 20, 3), Point3D(20, 20, 3), Point3D(20, 10, 3))
room2d_1 = Room2D('Office1', Face3D(pts_1), 3)
room2d_2 = Room2D('Office2', Face3D(pts_2), 3)
room2d_3 = Room2D('Office3', Face3D(pts_3), 3)
room2d_4 = Room2D('Office4', Face3D(pts_4), 3)
story = Story('OfficeFloor', [room2d_1, room2d_2, room2d_3, room2d_4])
story.solve_room_2d_adjacency(0.01)
story.set_outdoor_window_parameters(SimpleWindowRatio(0.4))
story.multiplier = 4
building = Building('OfficeBuilding', [story])
sensible = 0.8
latent = 0.7
ideal_air_sys = IdealAirSystem('Office Ideal Air', sensible_heat_recovery=sensible,
latent_heat_recovery=latent)
building.properties.energy.set_all_room_2d_hvac(ideal_air_sys, False)
assert all(isinstance(room.properties.energy.hvac, IdealAirSystem)
for room in building.unique_room_2ds)
assert all(room.properties.energy.hvac.sensible_heat_recovery == sensible
for room in building.unique_room_2ds)
assert all(room.properties.energy.hvac.latent_heat_recovery == latent
for room in building.unique_room_2ds)
def test_scale():
"""Test the Model scale method."""
pts_1 = (Point3D(0, 0, 3), Point3D(10, 0, 3), Point3D(10, 10, 3), Point3D(0, 10, 3))
pts_2 = (Point3D(10, 0, 3), Point3D(20, 0, 3), Point3D(20, 10, 3), Point3D(10, 10, 3))
room2d_1 = Room2D('Office1', Face3D(pts_1), 3)
room2d_2 = Room2D('Office2', Face3D(pts_2), 3)
story = Story('OfficeFloor', [room2d_1, room2d_2])
story.solve_room_2d_adjacency(0.01)
story.multiplier = 4
building = Building('OfficeBuilding', [story])
awning_geo1 = Face3D.from_rectangle(6, 6, Plane(o=Point3D(5, -10, 6)))
awning_geo2 = Face3D.from_rectangle(2, 2, Plane(o=Point3D(-5, -10, 3)))
awning_canopy_1 = ContextShade('AwningCanopy1', [awning_geo1])
awning_canopy_2 = ContextShade('AwningCanopy2', [awning_geo2])
model = Model('NewDevelopment', [building], [awning_canopy_1, awning_canopy_2])
new_m = model.duplicate()
new_m.scale(2)
assert new_m.buildings[0].unique_stories[0].room_2ds[0].floor_geometry[0] == Point3D(0, 0, 6)
assert new_m.buildings[0].unique_stories[0].room_2ds[0].floor_geometry[1] == Point3D(20, 0, 6)
assert new_m.buildings[0].unique_stories[0].room_2ds[0].floor_geometry[2] == Point3D(20, 20, 6)
assert new_m.buildings[0].unique_stories[0].room_2ds[0].floor_geometry[3] == Point3D(0, 20, 6)
assert new_m.buildings[0].unique_stories[0].room_2ds[1].floor_geometry[0] == Point3D(20, 0, 6)
assert new_m.buildings[0].unique_stories[0].room_2ds[1].floor_geometry[1] == Point3D(40, 0, 6)
assert new_m.buildings[0].unique_stories[0].room_2ds[1].floor_geometry[2] == Point3D(40, 20, 6)
assert new_m.buildings[0].unique_stories[0].room_2ds[1].floor_geometry[3] == Point3D(20, 20, 6)
assert new_m.buildings[0].floor_area == building.floor_area * 2 ** 2
assert new_m.context_shades[0][0][0] == Point3D(10, -20, 12)
assert new_m.context_shades[1][0][0] == Point3D(-10, -20, 6)
def test_energy_properties():
"""Test the existence of the Model energy properties."""
pts_1 = (Point3D(0, 0, 3), Point3D(0, 10,
3), Point3D(10, 10,
3), Point3D(10, 0, 3))
pts_2 = (Point3D(10, 0, 3), Point3D(10, 10,
3), Point3D(20, 10,
3), Point3D(20, 0, 3))
pts_3 = (Point3D(0, 10, 3), Point3D(0, 20,
3), Point3D(10, 20,
3), Point3D(10, 10, 3))
pts_4 = (Point3D(10, 10, 3), Point3D(10, 20,
3), Point3D(20, 20,
3), Point3D(20, 10, 3))
room2d_1 = Room2D('Office1', Face3D(pts_1), 3)
room2d_2 = Room2D('Office2', Face3D(pts_2), 3)
room2d_3 = Room2D('Office3', Face3D(pts_3), 3)
room2d_4 = Room2D('Office4', Face3D(pts_4), 3)
story = Story('OfficeFloor', [room2d_1, room2d_2, room2d_3, room2d_4])
story.solve_room_2d_adjacency(0.01)
story.set_outdoor_window_parameters(SimpleWindowRatio(0.4))
story.multiplier = 4
for room in story.room_2ds:
room.properties.energy.program_type = office_program
room.properties.energy.add_default_ideal_air()
building = Building('OfficeBuilding', [story])
tree_canopy_geo1 = Face3D.from_regular_polygon(6, 6,
Plane(o=Point3D(5, -10, 6)))
tree_canopy_geo2 = Face3D.from_regular_polygon(
6, 2, Plane(o=Point3D(-5, -10, 3)))
tree_canopy = ContextShade('TreeCanopy',
[tree_canopy_geo1, tree_canopy_geo2])
bright_leaves = ShadeConstruction('Bright Light Leaves', 0.5, 0.5, True)
tree_canopy.properties.energy.construction = bright_leaves
tree_trans = ScheduleRuleset.from_constant_value('Tree Transmittance', 0.5,
schedule_types.fractional)
tree_canopy.properties.energy.transmittance_schedule = tree_trans
model = Model('NewDevelopment', [building], [tree_canopy])
assert hasattr(model.properties, 'energy')
assert isinstance(model.properties.energy, ModelEnergyProperties)
assert isinstance(model.properties.host, Model)
assert len(model.properties.energy.materials) == 0
for mat in model.properties.energy.materials:
assert isinstance(mat, _EnergyMaterialBase)
assert len(model.properties.energy.constructions) == 1
for cnst in model.properties.energy.constructions:
assert isinstance(cnst, (WindowConstruction, OpaqueConstruction,
ShadeConstruction, AirBoundaryConstruction))
assert len(model.properties.energy.shade_constructions) == 1
assert len(model.properties.energy.construction_sets) == 0
assert len(model.properties.energy.schedule_type_limits) == 3
assert len(model.properties.energy.schedules) == 8
assert len(model.properties.energy.shade_schedules) == 1
assert len(model.properties.energy.program_types) == 1
def test_model_add_objects():
"""Test the addition of objects to a Model and getting objects by identifier."""
pts_1 = (Point3D(0, 0, 3), Point3D(0, 10, 3), Point3D(10, 10, 3), Point3D(10, 0, 3))
pts_2 = (Point3D(10, 0, 3), Point3D(10, 10, 3), Point3D(20, 10, 3), Point3D(20, 0, 3))
pts_3 = (Point3D(0, 20, 3), Point3D(0, 30, 3), Point3D(10, 30, 3), Point3D(10, 20, 3))
pts_4 = (Point3D(10, 20, 3), Point3D(10, 30, 3), Point3D(20, 30, 3), Point3D(20, 20, 3))
room2d_1 = Room2D('Office1', Face3D(pts_1), 3)
room2d_2 = Room2D('Office2', Face3D(pts_2), 3)
room2d_3 = Room2D('Office3', Face3D(pts_3), 3)
room2d_4 = Room2D('Office4', Face3D(pts_4), 3)
story_1 = Story('OfficeFloor1', [room2d_1, room2d_2])
story_2 = Story('OfficeFloor2', [room2d_3, room2d_4])
story_1.solve_room_2d_adjacency(0.01)
story_1.set_outdoor_window_parameters(SimpleWindowRatio(0.4))
story_1.multiplier = 4
story_2.solve_room_2d_adjacency(0.01)
story_2.set_outdoor_window_parameters(SimpleWindowRatio(0.4))
story_2.multiplier = 2
building_1 = Building('OfficeBuilding1', [story_1])
building_2 = Building('OfficeBuilding2', [story_2])
tree_canopy_geo1 = Face3D.from_regular_polygon(6, 6, Plane(o=Point3D(5, -10, 6)))
tree_canopy_geo2 = Face3D.from_regular_polygon(6, 2, Plane(o=Point3D(-5, -10, 3)))
tree_canopy_1 = ContextShade('TreeCanopy1', [tree_canopy_geo1])
tree_canopy_2 = ContextShade('TreeCanopy2', [tree_canopy_geo2])
model = Model('NewDevelopment', [building_1], [tree_canopy_1])
assert len(model.buildings) == 1
assert len(model.context_shades) == 1
with pytest.raises(AssertionError):
model.add_building(tree_canopy_2)
model.add_building(building_2)
assert len(model.buildings) == 2
with pytest.raises(AssertionError):
model.add_context_shade(building_2)
model.add_context_shade(tree_canopy_2)
assert len(model.context_shades) == 2
assert len(model.buildings_by_identifier(['OfficeBuilding1'])) == 1
with pytest.raises(ValueError):
model.buildings_by_identifier(['NotABuilding'])
assert len(model.context_shade_by_identifier(['TreeCanopy1'])) == 1
with pytest.raises(ValueError):
model.context_shade_by_identifier(['NotAShade'])
def two_buildings_with_terrain():
pts_1 = (Point3D(0, 0, 3), Point3D(10, 0,
3), Point3D(10, 10,
3), Point3D(0, 10, 3))
pts_2 = (Point3D(10, 0, 3), Point3D(20, 0,
3), Point3D(20, 10,
3), Point3D(10, 10, 3))
pts_3 = (Point3D(0, 20, 3), Point3D(20, 20,
3), Point3D(20, 30,
3), Point3D(0, 30, 3))
room2d_1 = Room2D('Office1', Face3D(pts_1), 3)
room2d_2 = Room2D('Office2', Face3D(pts_2), 3)
room2d_3 = Room2D('Office3', Face3D(pts_3), 3)
story_big = Story('OfficeFloorBig', [room2d_3])
story = Story('OfficeFloor', [room2d_1, room2d_2])
story.solve_room_2d_adjacency(0.01)
story.set_outdoor_window_parameters(SimpleWindowRatio(0.3))
story.multiplier = 3
building = Building('OfficeBuilding', [story])
building.separate_top_bottom_floors()
story_big.set_outdoor_window_parameters(SimpleWindowRatio(0.6))
story_big.multiplier = 4
building_big = Building('OfficeBuildingBig', [story_big])
building_big.separate_top_bottom_floors()
tree_canopy_geo1 = Face3D.from_regular_polygon(6, 6,
Plane(o=Point3D(5, -10, 6)))
tree_canopy_geo2 = Face3D.from_regular_polygon(
6, 2, Plane(o=Point3D(-5, -10, 3)))
tree_canopy = ContextShade('TreeCanopy',
[tree_canopy_geo1, tree_canopy_geo2])
model = Model('NewDevelopment', [building, building_big], [tree_canopy])
ter_geo = Face3D(
(Point3D(0, 0, 0), Point3D(30, 0, 0), Point3D(30, 30,
0), Point3D(0, 30, 3)))
terrain = Terrain([ter_geo])
model.properties.uwg.terrain = terrain
model.properties.uwg.traffic.watts_per_area = 7
model.properties.uwg.tree_coverage_fraction = 0.25
model.properties.uwg.grass_coverage_fraction = 0.4
return model
def test_model_add_model():
"""Test the addition of one Model to another."""
pts_1 = (Point3D(0, 0, 3), Point3D(0, 10, 3), Point3D(10, 10, 3), Point3D(10, 0, 3))
pts_2 = (Point3D(10, 0, 3), Point3D(10, 10, 3), Point3D(20, 10, 3), Point3D(20, 0, 3))
pts_3 = (Point3D(0, 20, 3), Point3D(0, 30, 3), Point3D(10, 30, 3), Point3D(10, 20, 3))
pts_4 = (Point3D(10, 20, 3), Point3D(10, 30, 3), Point3D(20, 30, 3), Point3D(20, 20, 3))
room2d_1 = Room2D('Office1', Face3D(pts_1), 3)
room2d_2 = Room2D('Office2', Face3D(pts_2), 3)
room2d_3 = Room2D('Office3', Face3D(pts_3), 3)
room2d_4 = Room2D('Office4', Face3D(pts_4), 3)
story_1 = Story('OfficeFloor1', [room2d_1, room2d_2])
story_2 = Story('OfficeFloor2', [room2d_3, room2d_4])
story_1.solve_room_2d_adjacency(0.01)
story_1.set_outdoor_window_parameters(SimpleWindowRatio(0.4))
story_1.multiplier = 4
story_2.solve_room_2d_adjacency(0.01)
story_2.set_outdoor_window_parameters(SimpleWindowRatio(0.4))
story_2.multiplier = 2
building_1 = Building('OfficeBuilding1', [story_1])
building_2 = Building('OfficeBuilding2', [story_2])
tree_canopy_geo1 = Face3D.from_regular_polygon(6, 6, Plane(o=Point3D(5, -10, 6)))
tree_canopy_geo2 = Face3D.from_regular_polygon(6, 2, Plane(o=Point3D(-5, -10, 3)))
tree_canopy_1 = ContextShade('TreeCanopy1', [tree_canopy_geo1])
tree_canopy_2 = ContextShade('TreeCanopy2', [tree_canopy_geo2])
model_1 = Model('NewDevelopment1', [building_1], [tree_canopy_1])
model_2 = Model('NewDevelopment2', [building_2], [tree_canopy_2])
assert len(model_1.buildings) == 1
assert len(model_1.context_shades) == 1
assert len(model_2.buildings) == 1
assert len(model_2.context_shades) == 1
combined_model = model_1 + model_2
assert len(combined_model.buildings) == 2
assert len(combined_model.context_shades) == 2
assert len(model_1.buildings) == 1
assert len(model_1.context_shades) == 1
model_1 += model_2
assert len(model_1.buildings) == 2
assert len(model_1.context_shades) == 2
def test_building_init():
"""Test the initialization of Building objects and basic properties."""
pts_1 = (Point3D(0, 0, 3), Point3D(0, 10,
3), Point3D(10, 10,
3), Point3D(10, 0, 3))
pts_2 = (Point3D(10, 0, 3), Point3D(10, 10,
3), Point3D(20, 10,
3), Point3D(20, 0, 3))
pts_3 = (Point3D(0, 10, 3), Point3D(0, 20,
3), Point3D(10, 20,
3), Point3D(10, 10, 3))
pts_4 = (Point3D(10, 10, 3), Point3D(10, 20,
3), Point3D(20, 20,
3), Point3D(20, 10, 3))
room2d_1 = Room2D('Office1', Face3D(pts_1), 3)
room2d_2 = Room2D('Office2', Face3D(pts_2), 3)
room2d_3 = Room2D('Office3', Face3D(pts_3), 3)
room2d_4 = Room2D('Office4', Face3D(pts_4), 3)
story = Story('Office_Floor', [room2d_1, room2d_2, room2d_3, room2d_4])
story.solve_room_2d_adjacency(0.01)
story.set_outdoor_window_parameters(SimpleWindowRatio(0.4))
story.multiplier = 4
building = Building('Office_Building_1234', [story])
building.display_name = 'Office Building'
str(building) # test the string representation
assert building.identifier == 'Office_Building_1234'
assert building.display_name == 'Office Building'
assert len(building.unique_stories) == len(
building.unique_stories_above_ground) == 1
assert len(building.all_stories()) == 4
assert len(building.unique_room_2ds) == 4
assert len(building.all_room_2ds()) == 16
for story in building.unique_stories:
assert isinstance(story, Story)
assert story.has_parent
for story in building.all_stories():
assert isinstance(story, Story)
assert story.has_parent
for room in building.unique_room_2ds:
assert isinstance(room, Room2D)
assert room.has_parent
for room in building.all_room_2ds():
assert isinstance(room, Room2D)
assert room.has_parent
assert building.height == 15
assert building.story_count == building.story_count_above_ground == 4
assert building.height_from_first_floor == building.height_above_ground == 12
assert building.footprint_area == 100 * 4
assert building.floor_area == 100 * 4 * 4
assert building.exterior_wall_area == 60 * 4 * 4
assert building.exterior_aperture_area == 60 * 4 * 4 * 0.4
assert building.volume == 100 * 3 * 4 * 4
def test_check_duplicate_identifiers():
"""Test check_duplicate_building_identifiers."""
pts_1 = (Point3D(0, 0, 3), Point3D(0, 10, 3), Point3D(10, 10, 3), Point3D(10, 0, 3))
pts_2 = (Point3D(10, 0, 3), Point3D(10, 10, 3), Point3D(20, 10, 3), Point3D(20, 0, 3))
pts_3 = (Point3D(0, 20, 3), Point3D(0, 30, 3), Point3D(10, 30, 3), Point3D(10, 20, 3))
pts_4 = (Point3D(10, 20, 3), Point3D(10, 30, 3), Point3D(20, 30, 3), Point3D(20, 20, 3))
room2d_1 = Room2D('Office1', Face3D(pts_1), 3)
room2d_2 = Room2D('Office2', Face3D(pts_2), 3)
room2d_3 = Room2D('Office3', Face3D(pts_3), 3)
room2d_4 = Room2D('Office4', Face3D(pts_4), 3)
story_1 = Story('OfficeFloor1', [room2d_1, room2d_2])
story_2 = Story('OfficeFloor2', [room2d_3, room2d_4])
story_1.solve_room_2d_adjacency(0.01)
story_1.set_outdoor_window_parameters(SimpleWindowRatio(0.4))
story_1.multiplier = 4
story_2.solve_room_2d_adjacency(0.01)
story_2.set_outdoor_window_parameters(SimpleWindowRatio(0.4))
story_2.multiplier = 2
building_1 = Building('OfficeBuilding', [story_1])
building_2 = Building('OfficeBuilding', [story_2])
tree_canopy_geo1 = Face3D.from_regular_polygon(6, 6, Plane(o=Point3D(5, -10, 6)))
tree_canopy_geo2 = Face3D.from_regular_polygon(6, 2, Plane(o=Point3D(-5, -10, 3)))
tree_canopy_1 = ContextShade('TreeCanopy', [tree_canopy_geo1])
tree_canopy_2 = ContextShade('TreeCanopy', [tree_canopy_geo2])
model_1 = Model('NewDevelopment1', [building_1], [tree_canopy_1])
model_2 = Model('NewDevelopment2', [building_2], [tree_canopy_2])
assert model_1.check_duplicate_building_identifiers(False) == ''
assert model_1.check_duplicate_context_shade_identifiers(False) == ''
model_1.add_model(model_2)
assert model_1.check_duplicate_building_identifiers(False) != ''
with pytest.raises(ValueError):
model_1.check_duplicate_building_identifiers(True)
assert model_1.check_duplicate_context_shade_identifiers(False) != ''
with pytest.raises(ValueError):
model_1.check_duplicate_context_shade_identifiers(True)
def test_check_duplicate_construction_set_identifiers():
"""Test the check_duplicate_construction_set_identifiers method."""
pts_1 = (Point3D(0, 0, 3), Point3D(0, 10,
3), Point3D(10, 10,
3), Point3D(10, 0, 3))
pts_2 = (Point3D(10, 0, 3), Point3D(10, 10,
3), Point3D(20, 10,
3), Point3D(20, 0, 3))
room2d_1 = Room2D('Office1', Face3D(pts_1), 3)
room2d_2 = Room2D('Office2', Face3D(pts_2), 3)
story = Story('OfficeFloor', [room2d_1, room2d_2])
story.solve_room_2d_adjacency(0.01)
story.set_outdoor_window_parameters(SimpleWindowRatio(0.4))
story.multiplier = 4
for room in story.room_2ds:
room.properties.energy.program_type = office_program
room.properties.energy.add_default_ideal_air()
building = Building('OfficeBuilding', [story])
building.separate_top_bottom_floors()
constr_set = ConstructionSet('Attic Construction Set')
polyiso = EnergyMaterial('PolyIso', 0.2, 0.03, 43, 1210, 'MediumRough')
roof_constr = OpaqueConstruction('Attic Roof Construction',
[roof_membrane, polyiso, wood])
floor_constr = OpaqueConstruction('Attic Floor Construction',
[wood, insulation, wood])
constr_set.floor_set.interior_construction = floor_constr
constr_set.roof_ceiling_set.exterior_construction = roof_constr
building.unique_room_2ds[
-1].properties.energy.construction_set = constr_set
building.unique_room_2ds[
-2].properties.energy.construction_set = constr_set
tree_canopy_geo1 = Face3D.from_regular_polygon(6, 6,
Plane(o=Point3D(5, -10, 6)))
tree_canopy_geo2 = Face3D.from_regular_polygon(
6, 2, Plane(o=Point3D(-5, -10, 3)))
tree_canopy = ContextShade('TreeCanopy',
[tree_canopy_geo1, tree_canopy_geo2])
model = Model('NewDevelopment', [building], [tree_canopy])
assert model.properties.energy.check_duplicate_construction_set_identifiers(
False)
constr_set2 = ConstructionSet('Attic Construction Set')
building.unique_room_2ds[
-2].properties.energy.construction_set = constr_set2
assert not model.properties.energy.check_duplicate_construction_set_identifiers(
False)
with pytest.raises(ValueError):
model.properties.energy.check_duplicate_construction_set_identifiers(
True)
def default_building():
pts_1 = (Point3D(0, 0, 3), Point3D(10, 0, 3), Point3D(10, 10, 3), Point3D(0, 10, 3))
pts_2 = (Point3D(10, 0, 3), Point3D(20, 0, 3), Point3D(20, 10, 3),
Point3D(10, 10, 3))
room2d_1 = Room2D('Office1', Face3D(pts_1), 3)
room2d_2 = Room2D('Office2', Face3D(pts_2), 3)
story = Story('OfficeFloor', [room2d_1, room2d_2])
story.solve_room_2d_adjacency(0.01)
story.set_outdoor_window_parameters(SimpleWindowRatio(0.3))
story.multiplier = 3
building = Building('OfficeBuilding', [story])
building.separate_top_bottom_floors()
return building
def test_to_honeybee():
"""Test the to_honeybee method."""
pts_1 = (Point3D(0, 0, 3), Point3D(10, 0,
3), Point3D(10, 10,
3), Point3D(0, 10, 3))
pts_2 = (Point3D(10, 0, 3), Point3D(20, 0,
3), Point3D(20, 10,
3), Point3D(10, 10, 3))
room2d_1 = Room2D('Office1', Face3D(pts_1), 3)
room2d_2 = Room2D('Office2', Face3D(pts_2), 3)
story = Story('Office_Floor', [room2d_1, room2d_2])
story.solve_room_2d_adjacency(0.01)
story.set_outdoor_window_parameters(SimpleWindowRatio(0.4))
story.multiplier = 4
building = Building('Office_Building', [story])
hb_model = building.to_honeybee(False, 0.01)
assert isinstance(hb_model, Model)
assert len(hb_model.rooms) == 8
assert len(hb_model.rooms[0]) == 6
assert hb_model.rooms[0].volume == 300
assert hb_model.rooms[0].floor_area == 100
assert hb_model.rooms[0].exterior_wall_area == 90
assert hb_model.rooms[0].exterior_aperture_area == pytest.approx(90 * 0.4,
rel=1e-3)
assert hb_model.rooms[0].average_floor_height == 3
assert hb_model.rooms[0].check_solid(0.01, 1) == ''
assert isinstance(hb_model.rooms[0][1].boundary_condition, Outdoors)
assert isinstance(hb_model.rooms[0][2].boundary_condition, Surface)
assert hb_model.rooms[0][2].boundary_condition.boundary_condition_object == \
hb_model.rooms[1][4].identifier
assert len(hb_model.rooms[0][1].apertures) == 1
assert len(hb_model.rooms[0][2].apertures) == 0
assert hb_model.check_duplicate_room_identifiers() == ''
assert hb_model.check_duplicate_face_identifiers() == ''
assert hb_model.check_duplicate_sub_face_identifiers() == ''
assert hb_model.check_missing_adjacencies() == ''
hb_model = building.to_honeybee(True, 0.01)
assert len(hb_model.rooms) == 2
for room in hb_model.rooms:
assert room.multiplier == 4
assert hb_model.check_duplicate_room_identifiers() == ''
assert hb_model.check_duplicate_face_identifiers() == ''
assert hb_model.check_duplicate_sub_face_identifiers() == ''
assert hb_model.check_missing_adjacencies() == ''
def test_check_duplicate_program_type_identifiers():
"""Test the check_duplicate_program_type_identifiers method."""
pts_1 = (Point3D(0, 0, 3), Point3D(0, 10,
3), Point3D(10, 10,
3), Point3D(10, 0, 3))
pts_2 = (Point3D(10, 0, 3), Point3D(10, 10,
3), Point3D(20, 10,
3), Point3D(20, 0, 3))
room2d_1 = Room2D('Office1', Face3D(pts_1), 3)
room2d_2 = Room2D('Office2', Face3D(pts_2), 3)
story = Story('OfficeFloor', [room2d_1, room2d_2])
story.solve_room_2d_adjacency(0.01)
story.set_outdoor_window_parameters(SimpleWindowRatio(0.4))
story.multiplier = 4
for room in story.room_2ds:
room.properties.energy.program_type = office_program
room.properties.energy.add_default_ideal_air()
building = Building('OfficeBuilding', [story])
building.separate_top_bottom_floors()
attic_program_type = plenum_program.duplicate()
attic_program_type.identifier = 'Attic Space'
schedule = ScheduleRuleset.from_constant_value('Always Dim', 1,
schedule_types.fractional)
lighting = Lighting('Attic Lighting', 3, schedule)
attic_program_type.lighting = lighting
building.unique_room_2ds[
-1].properties.energy.program_type = attic_program_type
building.unique_room_2ds[
-2].properties.energy.program_type = attic_program_type
tree_canopy_geo1 = Face3D.from_regular_polygon(6, 6,
Plane(o=Point3D(5, -10, 6)))
tree_canopy_geo2 = Face3D.from_regular_polygon(
6, 2, Plane(o=Point3D(-5, -10, 3)))
tree_canopy = ContextShade('TreeCanopy',
[tree_canopy_geo1, tree_canopy_geo2])
model = Model('NewDevelopment', [building], [tree_canopy])
assert model.properties.energy.check_duplicate_program_type_identifiers(
False)
attic_program_type.unlock()
attic_program_type.identifier = office_program.identifier
attic_program_type.lock()
assert not model.properties.energy.check_duplicate_program_type_identifiers(
False)
with pytest.raises(ValueError):
model.properties.energy.check_duplicate_program_type_identifiers(True)
def test_to_honeybee_missing_adjacency():
"""Test the to_honeybee method with a missing adjacency."""
pts_1 = (Point3D(0, 0, 3), Point3D(10, 0, 3), Point3D(10, 10, 3), Point3D(0, 10, 3))
pts_2 = (Point3D(10, 0, 3), Point3D(20, 0, 3), Point3D(20, 10, 3), Point3D(10, 10, 3))
pts_3 = (Point3D(0, 20, 3), Point3D(20, 20, 3), Point3D(20, 30, 3), Point3D(0, 30, 3))
room2d_1 = Room2D('Office1', Face3D(pts_1), 3)
room2d_2 = Room2D('Office2', Face3D(pts_2), 3)
room2d_3 = Room2D('Office3', Face3D(pts_3), 3)
story_big = Story('OfficeFloorBig', [room2d_3])
Room2D.solve_adjacency([room2d_1, room2d_2])
story = Story('OfficeFloor', [room2d_1])
story.set_outdoor_window_parameters(SimpleWindowRatio(0.4))
story.multiplier = 4
building = Building('OfficeBuilding', [story])
story_big.set_outdoor_window_parameters(SimpleWindowRatio(0.4))
story_big.multiplier = 4
building_big = Building('OfficeBuildingBig', [story_big])
model = Model('NewDevelopment', [building, building_big])
hb_models = model.to_honeybee('District', None, False, tolerance=0.01)
assert len(hb_models) == 1
assert isinstance(hb_models[0], hb_model.Model)
def test_rotate_xy():
"""Test the Building rotate_xy method."""
pts = (Point3D(1, 1, 2), Point3D(2, 1, 2), Point3D(2, 2,
2), Point3D(1, 2, 2))
plane = Plane(Vector3D(0, 0, 1), Point3D(0, 0, 2))
room = Room2D('Square_Shoebox', Face3D(pts, plane), 3)
story = Story('Office_Floor', [room])
story.multiplier = 4
building = Building('Office_Building', [story])
origin_1 = Point3D(1, 1, 0)
test_1 = building.duplicate()
test_1.rotate_xy(180, origin_1)
assert test_1.unique_stories[0].room_2ds[0].floor_geometry[
0].x == pytest.approx(1, rel=1e-3)
assert test_1.unique_stories[0].room_2ds[0].floor_geometry[
0].y == pytest.approx(1, rel=1e-3)
assert test_1.unique_stories[0].room_2ds[0].floor_geometry[
0].z == pytest.approx(2, rel=1e-3)
assert test_1.unique_stories[0].room_2ds[0].floor_geometry[
2].x == pytest.approx(0, rel=1e-3)
assert test_1.unique_stories[0].room_2ds[0].floor_geometry[
2].y == pytest.approx(0, rel=1e-3)
assert test_1.unique_stories[0].room_2ds[0].floor_geometry[
2].z == pytest.approx(2, rel=1e-3)
test_2 = building.duplicate()
test_2.rotate_xy(90, origin_1)
assert test_2.unique_stories[0].room_2ds[0].floor_geometry[
0].x == pytest.approx(1, rel=1e-3)
assert test_2.unique_stories[0].room_2ds[0].floor_geometry[
0].y == pytest.approx(1, rel=1e-3)
assert test_2.unique_stories[0].room_2ds[0].floor_geometry[
0].z == pytest.approx(2, rel=1e-3)
assert test_2.unique_stories[0].room_2ds[0].floor_geometry[
2].x == pytest.approx(0, rel=1e-3)
assert test_2.unique_stories[0].room_2ds[0].floor_geometry[
2].y == pytest.approx(2, rel=1e-3)
assert test_2.unique_stories[0].room_2ds[0].floor_geometry[
2].z == pytest.approx(2, rel=1e-3)
assert building.floor_area == test_1.floor_area
assert building.volume == test_1.volume
assert building.height_from_first_floor == test_1.height_from_first_floor
assert building.height == 14
def test_model_init():
"""Test the initialization of Model objects and basic properties."""
pts_1 = (Point3D(0, 0, 3), Point3D(0, 10, 3), Point3D(10, 10, 3), Point3D(10, 0, 3))
pts_2 = (Point3D(10, 0, 3), Point3D(10, 10, 3), Point3D(20, 10, 3), Point3D(20, 0, 3))
pts_3 = (Point3D(0, 10, 3), Point3D(0, 20, 3), Point3D(10, 20, 3), Point3D(10, 10, 3))
pts_4 = (Point3D(10, 10, 3), Point3D(10, 20, 3), Point3D(20, 20, 3), Point3D(20, 10, 3))
room2d_1 = Room2D('Office1', Face3D(pts_1), 3)
room2d_2 = Room2D('Office2', Face3D(pts_2), 3)
room2d_3 = Room2D('Office3', Face3D(pts_3), 3)
room2d_4 = Room2D('Office4', Face3D(pts_4), 3)
story = Story('Office_Floor', [room2d_1, room2d_2, room2d_3, room2d_4])
story.solve_room_2d_adjacency(0.01)
story.set_outdoor_window_parameters(SimpleWindowRatio(0.4))
story.multiplier = 4
building = Building('Office_Building', [story])
tree_canopy_geo1 = Face3D.from_regular_polygon(6, 6, Plane(o=Point3D(5, -10, 6)))
tree_canopy_geo2 = Face3D.from_regular_polygon(6, 2, Plane(o=Point3D(-5, -10, 3)))
tree_canopy = ContextShade('Tree_Canopy', [tree_canopy_geo1, tree_canopy_geo2])
model = Model('New_Development', [building], [tree_canopy])
str(model) # test the string representation of the object
assert model.identifier == 'New_Development'
assert model.display_name == 'New_Development'
assert model.units == 'Meters'
assert model.tolerance == 0.01
assert model.angle_tolerance == 1.0
assert len(model.buildings) == 1
assert isinstance(model.buildings[0], Building)
assert len(model.context_shades) == 1
assert isinstance(model.context_shades[0], ContextShade)
assert model.average_story_count == 4
assert model.average_story_count_above_ground == 4
assert model.average_height == 15
assert model.average_height_above_ground == 12
assert model.footprint_area == 100 * 4
assert model.floor_area == 100 * 4 * 4
assert model.exterior_wall_area == 60 * 4 * 4
assert model.exterior_aperture_area == 60 * 4 * 4 * 0.4
assert model.volume == 100 * 3 * 4 * 4
assert model.min.x == pytest.approx(-6.73, rel=1e-2)
assert model.min.y == pytest.approx(-16, rel=1e-2)
assert model.max == Point2D(20, 20)
