def test_Schedule(self):
""" Tests for Schedule (Schdef class) """
refDOE, refBEM, Schedule = uwg.readDOE(serialize_output=False)
schlst = [
'Elec', 'Light', 'Gas', 'Occ', 'Cool', 'Heat', 'SWH', 'Qelec',
'Qlight', 'Nocc', 'Qgas', 'Vent', 'Vswh'
]
for si in range(len(schlst)):
schid = schlst[si]
# Define file path
matlab_path = os.path.join(self.DIR_MATLAB_PATH,
"matlab_ref_sch_{}_.txt".format(schid))
# check file exists
if not os.path.exists(matlab_path):
raise Exception("Failed to open {}!".format(matlab_path))
matlab_file = open(matlab_path, 'r')
assert len(Schedule) == pytest.approx(16.0, abs=1e-3)
for bldType in range(1): # bldType
assert len(Schedule[bldType]) == pytest.approx(3.0, abs=1e-3)
for bldEra in range(1): # bltEra
assert len(Schedule[bldType][bldEra]) == pytest.approx(
16.0, abs=1e-3)
for climateZone in range(1): # ZoneType
if schid in schlst[:
7]: #3x24 matrix of schedule for SWH (WD,Sat,Sun)
matlab_ref_str = matlab_file.readline()
# replace [,] -> "" && split at " "
matlab_ref_str = matlab_ref_str.replace(
"[",
"").replace("]",
"").replace(";",
"_").replace(" ", "_")
matlab_ref_str = "".join(matlab_ref_str.split())
matlab_ref_str = matlab_ref_str.split("_")
matlab_ref_value = [
float(x) for x in matlab_ref_str
]
# can compare equality of flat lists [0.12, 4, 5,1.12, 5, 6] == pytest.approx([0.124, 4., 5, 1.125, 5, 6], abs=1e-2)
assert len(matlab_ref_value) == pytest.approx(
24.0 * 3., abs=1e-3)
# Calculate tolerances and set minimum one for all
tol = min([
self.CALCULATE_TOLERANCE(x, 15.0)
for x in matlab_ref_value
])
else:
matlab_ref_value = float(matlab_file.readline())
tol = self.CALCULATE_TOLERANCE(
matlab_ref_value, 15.0)
if schid == 'Elec': #3x24 matrix of schedule for SWH (WD,Sat,Sun)
flatten_sch = reduce(
lambda x, y: x + y,
Schedule[bldType][bldEra][climateZone].Elec)
assert flatten_sch == pytest.approx(matlab_ref_value, abs=tol), \
'btype={},era={},czone={}'.format(bldType+1, bldEra+1, climateZone+1)
elif schid == 'Light': #3x24 matrix of schedule for SWH (WD,Sat,Sun)
flatten_sch = reduce(
lambda x, y: x + y,
Schedule[bldType][bldEra][climateZone].Light)
assert flatten_sch == pytest.approx(matlab_ref_value, abs=tol), \
'btype={},era={},czone={}'.format(bldType+1, bldEra+1, climateZone+1)
elif schid == 'Gas': #3x24 matrix of schedule for SWH (WD,Sat,Sun)
flatten_sch = reduce(
lambda x, y: x + y,
Schedule[bldType][bldEra][climateZone].Gas)
assert flatten_sch == pytest.approx(matlab_ref_value, abs=tol), \
'btype={},era={},czone={}'.format(bldType+1, bldEra+1, climateZone+1)
elif schid == 'Occ': #3x24 matrix of schedule for SWH (WD,Sat,Sun)
flatten_sch = reduce(
lambda x, y: x + y,
Schedule[bldType][bldEra][climateZone].Occ)
assert flatten_sch == pytest.approx(matlab_ref_value, abs=tol), \
'btype={},era={},czone={}'.format(bldType+1, bldEra+1, climateZone+1)
elif schid == 'Cool': #3x24 matrix of schedule for SWH (WD,Sat,Sun)
flatten_sch = reduce(
lambda x, y: x + y,
Schedule[bldType][bldEra][climateZone].Cool)
assert flatten_sch == pytest.approx(matlab_ref_value, abs=tol), \
'btype={},era={},czone={}'.format(bldType+1, bldEra+1, climateZone+1)
elif schid == 'Heat': #3x24 matrix of schedule for SWH (WD,Sat,Sun)
flatten_sch = reduce(
lambda x, y: x + y,
Schedule[bldType][bldEra][climateZone].Heat)
assert flatten_sch == pytest.approx(matlab_ref_value, abs=tol), \
'btype={},era={},czone={}'.format(bldType+1, bldEra+1, climateZone+1)
elif schid == 'SWH': #3x24 matrix of schedule for SWH (WD,Sat,Sun)
flatten_sch = reduce(
lambda x, y: x + y,
Schedule[bldType][bldEra][climateZone].SWH)
assert flatten_sch == pytest.approx(matlab_ref_value, abs=tol), \
'btype={},era={},czone={}'.format(bldType+1, bldEra+1, climateZone+1)
elif schid == 'Qelec':
assert Schedule[bldType][bldEra][climateZone].Qelec == pytest.approx(matlab_ref_value, abs=tol), \
'btype={},era={},czone={}'.format(bldType+1, bldEra+1, climateZone+1)
elif schid == 'Qlight':
assert Schedule[bldType][bldEra][climateZone].Qlight == pytest.approx(matlab_ref_value, abs=tol), \
'btype={},era={},czone={}'.format(bldType+1, bldEra+1, climateZone+1)
elif schid == 'Nocc':
assert Schedule[bldType][bldEra][climateZone].Nocc == pytest.approx(matlab_ref_value, abs=tol), \
'btype={},era={},czone={}'.format(bldType+1, bldEra+1, climateZone+1)
elif schid == 'Qgas':
assert Schedule[bldType][bldEra][climateZone].Qgas == pytest.approx(matlab_ref_value, abs=tol), \
'btype={},era={},czone={}'.format(bldType+1, bldEra+1, climateZone+1)
elif schid == 'Vent':
assert Schedule[bldType][bldEra][climateZone].Vent == pytest.approx(matlab_ref_value, abs=tol), \
'btype={},era={},czone={}'.format(bldType+1, bldEra+1, climateZone+1)
elif schid == 'Vswh':
assert Schedule[bldType][bldEra][climateZone].Vswh == pytest.approx(matlab_ref_value, abs=tol), \
'btype={},era={},czone={}'.format(bldType+1, bldEra+1, climateZone+1)
matlab_file.close()
def test_refDOE(self):
""" Tests for refDOE (Building class)"""
# Run readDOE.py
refDOE, refBEM, Schedule = uwg.readDOE(serialize_output=False)
bldlst = [
'floorHeight', 'intHeat', 'intHeatNight', 'intHeatDay',
'intHeatFRad', 'intHeatFLat', 'infil', 'vent', 'glazingRatio',
'uValue', 'shgc', 'condType', 'cop', 'coolSetpointDay',
'coolSetpointNight', 'heatSetpointDay', 'heatSetpointNight',
'coolCap', 'heatEff', 'mSys', 'indoorTemp', 'indoorHum', 'heatCap',
'copAdj', 'fanMax'
]
for bi in range(len(bldlst)):
bldid = bldlst[bi]
matlab_path = os.path.join(self.DIR_MATLAB_PATH,
"matlab_ref_{}_.txt".format(bldid))
# check file exists
if not os.path.exists(matlab_path):
raise Exception("Failed to open {}!".format(matlab_path))
matlab_file = open(matlab_path, 'r')
assert len(refDOE) == pytest.approx(16.0, abs=1e-3)
for bldType in range(16): # bldType
assert len(refDOE[bldType]) == pytest.approx(3.0, abs=1e-3)
for bldEra in range(3): # bltEra
assert len(refDOE[bldType][bldEra]) == pytest.approx(
16.0, abs=1e-3)
for climateZone in range(16): # ZoneType
matlab_ref_value = matlab_file.readline()
if bldid != 'condType':
matlab_ref_value = float(matlab_ref_value)
tol = self.CALCULATE_TOLERANCE(
matlab_ref_value, 14.0)
else:
matlab_ref_value = 'AIR'
# run tests
if bldid == 'floorHeight':
assert refDOE[bldType][bldEra][climateZone].floorHeight == pytest.approx(matlab_ref_value, abs=tol), \
'btype={},era={},czone={}'.format(bldType+1, bldEra+1, climateZone+1)
elif bldid == 'intHeat':
assert refDOE[bldType][bldEra][climateZone].intHeat == pytest.approx(matlab_ref_value, abs=tol), \
'btype={},era={},czone={}'.format(bldType+1, bldEra+1, climateZone+1)
elif bldid == 'intHeatNight':
assert refDOE[bldType][bldEra][climateZone].intHeatNight == pytest.approx(matlab_ref_value, abs=tol), \
'btype={},era={},czone={}'.format(bldType+1, bldEra+1, climateZone+1)
elif bldid == 'intHeatDay':
assert refDOE[bldType][bldEra][climateZone].intHeatDay == pytest.approx(matlab_ref_value, abs=tol), \
'btype={},era={},czone={}'.format(bldType+1, bldEra+1, climateZone+1)
elif bldid == 'intHeatFRad':
assert refDOE[bldType][bldEra][climateZone].intHeatFRad == pytest.approx(matlab_ref_value, abs=tol), \
'btype={},era={},czone={}'.format(bldType+1, bldEra+1, climateZone+1)
elif bldid == 'intHeatFLat':
assert refDOE[bldType][bldEra][climateZone].intHeatFLat == pytest.approx(matlab_ref_value, abs=tol), \
'btype={},era={},czone={}'.format(bldType+1, bldEra+1, climateZone+1)
elif bldid == 'infil':
assert refDOE[bldType][bldEra][climateZone].infil == pytest.approx(matlab_ref_value, abs=tol), \
'btype={},era={},czone={}'.format(bldType+1, bldEra+1, climateZone+1)
elif bldid == 'vent':
assert refDOE[bldType][bldEra][climateZone].vent == pytest.approx(matlab_ref_value, abs=tol), \
'btype={},era={},czone={}'.format(bldType+1, bldEra+1, climateZone+1)
elif bldid == 'glazingRatio':
assert refDOE[bldType][bldEra][climateZone].glazingRatio == pytest.approx(matlab_ref_value, abs=tol), \
'btype={},era={},czone={}'.format(bldType+1, bldEra+1, climateZone+1)
elif bldid == 'uValue':
assert refDOE[bldType][bldEra][climateZone].uValue == pytest.approx(matlab_ref_value, abs=tol), \
'btype={},era={},czone={}'.format(bldType+1, bldEra+1, climateZone+1)
elif bldid == 'shgc':
assert refDOE[bldType][bldEra][climateZone].shgc == pytest.approx(matlab_ref_value, abs=tol), \
'btype={},era={},czone={}'.format(bldType+1, bldEra+1, climateZone+1)
elif bldid == 'condType':
assert refDOE[bldType][bldEra][climateZone].condType == matlab_ref_value, \
'btype={},era={},czone={}'.format(bldType+1, bldEra+1, climateZone+1)
elif bldid == 'cop':
assert refDOE[bldType][bldEra][climateZone].cop == pytest.approx(matlab_ref_value, abs=tol), \
'btype={},era={},czone={}'.format(bldType+1, bldEra+1, climateZone+1)
elif bldid == 'coolSetpointDay':
assert refDOE[bldType][bldEra][climateZone].coolSetpointDay == pytest.approx(matlab_ref_value, abs=tol), \
'btype={},era={},czone={}'.format(bldType+1, bldEra+1, climateZone+1)
elif bldid == 'coolSetpointNight':
assert refDOE[bldType][bldEra][climateZone].coolSetpointNight == pytest.approx(matlab_ref_value, abs=tol), \
'btype={},era={},czone={}'.format(bldType+1, bldEra+1, climateZone+1)
elif bldid == 'heatSetpointDay':
assert refDOE[bldType][bldEra][climateZone].heatSetpointDay == pytest.approx(matlab_ref_value, abs=tol), \
'btype={},era={},czone={}'.format(bldType+1, bldEra+1, climateZone+1)
elif bldid == 'heatSetpointNight':
assert refDOE[bldType][bldEra][climateZone].heatSetpointNight == pytest.approx(matlab_ref_value, abs=tol), \
'btype={},era={},czone={}'.format(bldType+1, bldEra+1, climateZone+1)
elif bldid == 'coolCap': # knocking the testing tolerance one order of magnitude down due to excel extraction problem
assert refDOE[bldType][bldEra][climateZone].coolCap == pytest.approx(matlab_ref_value, abs=10**(math.log10(tol)+1)), \
'btype={},era={},czone={}'.format(bldType+1, bldEra+1, climateZone+1)
elif bldid == 'heatEff':
assert refDOE[bldType][bldEra][climateZone].heatEff == pytest.approx(matlab_ref_value, abs=tol), \
'btype={},era={},czone={}'.format(bldType+1, bldEra+1, climateZone+1)
elif bldid == 'mSys':
assert refDOE[bldType][bldEra][climateZone].mSys == pytest.approx(matlab_ref_value, abs=tol), \
'btype={},era={},czone={}'.format(bldType+1, bldEra+1, climateZone+1)
elif bldid == 'indoorTemp':
assert refDOE[bldType][bldEra][climateZone].indoorTemp == pytest.approx(matlab_ref_value, abs=tol), \
'btype={},era={},czone={}'.format(bldType+1, bldEra+1, climateZone+1)
elif bldid == 'indoorHum':
assert refDOE[bldType][bldEra][climateZone].indoorHum == pytest.approx(matlab_ref_value, abs=tol), \
'btype={},era={},czone={}'.format(bldType+1, bldEra+1, climateZone+1)
elif bldid == 'heatCap': # knocking the testing tolerance one order of magnitude down due to excel extraction problem
assert refDOE[bldType][bldEra][climateZone].heatCap == pytest.approx(matlab_ref_value, abs=10**(math.log10(tol)+1)), \
'btype={},era={},czone={}'.format(bldType+1, bldEra+1, climateZone+1)
elif bldid == 'copAdj':
assert refDOE[bldType][bldEra][climateZone].copAdj == pytest.approx(matlab_ref_value, abs=tol), \
'btype={},era={},czone={}'.format(bldType+1, bldEra+1, climateZone+1)
elif bldid == 'canyon_fraction':
assert refDOE[bldType][bldEra][climateZone].canyon_fraction == pytest.approx(matlab_ref_value, abs=tol), \
'btype={},era={},czone={}'.format(bldType+1, bldEra+1, climateZone+1)
elif bldid == 'fanMax':
assert refBEM[bldType][bldEra][climateZone].building.FanMax == pytest.approx(matlab_ref_value, abs=tol), \
'btype={},era={},czone={}'.format(bldType+1, bldEra+1, climateZone+1)
matlab_file.close()
def test_refBEM(self):
""" Tests for refBEM (BEMDef class) """
refDOE, refBEM, Schedule = uwg.readDOE(serialize_output=False)
elementlst = [
'albedo', 'emissivity', 'layerThickness', 'layerThermalCond',
'layerVolHeat', 'vegCoverage', 'layerTemp', 'horizontal'
]
bemlst = [("mass", copy.deepcopy(elementlst)),
("wall", copy.deepcopy(elementlst)),
("roof", copy.deepcopy(elementlst))]
for bemi in range(len(bemlst)):
for ei in range(len(elementlst)):
bemid = bemlst[bemi][0] + "_" + bemlst[bemi][1][ei]
matlab_path = os.path.join(
self.DIR_MATLAB_PATH,
"matlab_ref_bemdef_{}.txt".format(bemid))
# check file exists
if not os.path.exists(matlab_path):
raise Exception("Failed to open {}!".format(matlab_path))
matlab_file = open(matlab_path, 'r')
assert len(refBEM) == pytest.approx(16.0, abs=1e-3)
for bldType in range(16): # bldType
assert len(refBEM[bldType]) == pytest.approx(3.0, abs=1e-3)
for bldEra in range(3): # bltEra
assert len(refBEM[bldType][bldEra]) == pytest.approx(
16.0, abs=1e-3)
for climateZone in range(16): # ZoneType
# next line
matlab_ref_value = float(matlab_file.readline())
tol = self.CALCULATE_TOLERANCE(
matlab_ref_value, 15.0)
# run tests
if bemid == 'mass_albedo':
assert refBEM[bldType][bldEra][climateZone].mass.albedo == pytest.approx(matlab_ref_value, abs=tol), \
'btype={},era={},czone={}'.format(bldType+1, bldEra+1, climateZone+1)
elif bemid == 'mass_emissivity':
assert refBEM[bldType][bldEra][climateZone].mass.emissivity == pytest.approx(matlab_ref_value, abs=tol), \
'btype={},era={},czone={}'.format(bldType+1, bldEra+1, climateZone+1)
elif bemid == 'mass_layerThickness':
assert refBEM[bldType][bldEra][climateZone].mass.layerThickness[0] == pytest.approx(matlab_ref_value, abs=tol), \
'btype={},era={},czone={}'.format(bldType+1, bldEra+1, climateZone+1)
elif bemid == 'mass_layerThermalCond':
assert refBEM[bldType][bldEra][climateZone].mass.layerThermalCond[0] == pytest.approx(matlab_ref_value, abs=tol), \
'btype={},era={},czone={}'.format(bldType+1, bldEra+1, climateZone+1)
elif bemid == 'mass_layerVolHeat':
assert refBEM[bldType][bldEra][climateZone].mass.layerVolHeat[0] == pytest.approx(matlab_ref_value, abs=tol), \
'btype={},era={},czone={}'.format(bldType+1, bldEra+1, climateZone+1)
elif bemid == 'mass_vegCoverage':
assert refBEM[bldType][bldEra][climateZone].mass.vegCoverage == pytest.approx(matlab_ref_value, abs=tol), \
'btype={},era={},czone={}'.format(bldType+1, bldEra+1, climateZone+1)
elif bemid == 'mass_layerTemp':
assert refBEM[bldType][bldEra][climateZone].mass.layerTemp[0] == pytest.approx(matlab_ref_value, abs=tol), \
'btype={},era={},czone={}'.format(bldType+1, bldEra+1, climateZone+1)
elif bemid == 'mass_horizontal':
assert refBEM[bldType][bldEra][climateZone].mass.horizontal == pytest.approx(matlab_ref_value, abs=tol), \
'btype={},era={},czone={}'.format(bldType+1, bldEra+1, climateZone+1)
elif bemid == 'wall_albedo':
assert refBEM[bldType][bldEra][climateZone].wall.albedo == pytest.approx(matlab_ref_value, abs=tol), \
'btype={},era={},czone={}'.format(bldType+1, bldEra+1, climateZone+1)
elif bemid == 'wall_emissivity':
assert refBEM[bldType][bldEra][climateZone].wall.emissivity == pytest.approx(matlab_ref_value, abs=tol), \
'btype={},era={},czone={}'.format(bldType+1, bldEra+1, climateZone+1)
elif bemid == 'wall_layerThickness':
assert refBEM[bldType][bldEra][climateZone].wall.layerThickness[0] == pytest.approx(matlab_ref_value, abs=tol), \
'btype={},era={},czone={}'.format(bldType+1, bldEra+1, climateZone+1)
elif bemid == 'wall_layerThermalCond':
assert refBEM[bldType][bldEra][climateZone].wall.layerThermalCond[0] == pytest.approx(matlab_ref_value, abs=tol), \
'btype={},era={},czone={}'.format(bldType+1, bldEra+1, climateZone+1)
elif bemid == 'wall_layerVolHeat':
assert refBEM[bldType][bldEra][climateZone].wall.layerVolHeat[0] == pytest.approx(matlab_ref_value, abs=tol), \
'btype={},era={},czone={}'.format(bldType+1, bldEra+1, climateZone+1)
elif bemid == 'wall_vegCoverage':
assert refBEM[bldType][bldEra][climateZone].wall.vegCoverage == pytest.approx(matlab_ref_value, abs=tol), \
'btype={},era={},czone={}'.format(bldType+1, bldEra+1, climateZone+1)
elif bemid == 'wall_layerTemp':
assert refBEM[bldType][bldEra][climateZone].wall.layerTemp[0] == pytest.approx(matlab_ref_value, abs=tol), \
'btype={},era={},czone={}'.format(bldType+1, bldEra+1, climateZone+1)
elif bemid == 'wall_horizontal':
assert refBEM[bldType][bldEra][climateZone].wall.horizontal == pytest.approx(matlab_ref_value, abs=tol), \
'btype={},era={},czone={}'.format(bldType+1, bldEra+1, climateZone+1)
elif bemid == 'roof_albedo':
assert refBEM[bldType][bldEra][climateZone].roof.albedo == pytest.approx(matlab_ref_value, abs=tol), \
'btype={},era={},czone={}'.format(bldType+1, bldEra+1, climateZone+1)
elif bemid == 'roof_emissivity':
assert refBEM[bldType][bldEra][climateZone].roof.emissivity == pytest.approx(matlab_ref_value, abs=tol), \
'btype={},era={},czone={}'.format(bldType+1, bldEra+1, climateZone+1)
elif bemid == 'roof_layerThickness':
assert refBEM[bldType][bldEra][climateZone].roof.layerThickness[0] == pytest.approx(matlab_ref_value, abs=tol), \
'btype={},era={},czone={}'.format(bldType+1, bldEra+1, climateZone+1)
elif bemid == 'roof_layerThermalCond':
assert refBEM[bldType][bldEra][climateZone].roof.layerThermalCond[0] == pytest.approx(matlab_ref_value, abs=tol), \
'btype={},era={},czone={}'.format(bldType+1, bldEra+1, climateZone+1)
elif bemid == 'roof_layerVolHeat':
assert refBEM[bldType][bldEra][climateZone].roof.layerVolHeat[0] == pytest.approx(matlab_ref_value, abs=tol), \
'btype={},era={},czone={}'.format(bldType+1, bldEra+1, climateZone+1)
elif bemid == 'roof_vegCoverage':
assert refBEM[bldType][bldEra][climateZone].roof.vegCoverage == pytest.approx(matlab_ref_value, abs=tol), \
'btype={},era={},czone={}'.format(bldType+1, bldEra+1, climateZone+1)
elif bemid == 'roof_layerTemp':
assert refBEM[bldType][bldEra][climateZone].roof.layerTemp[0] == pytest.approx(matlab_ref_value, abs=tol), \
'btype={},era={},czone={}'.format(bldType+1, bldEra+1, climateZone+1)
elif bemid == 'roof_horizontal':
assert refBEM[bldType][bldEra][climateZone].roof.horizontal == pytest.approx(matlab_ref_value, abs=tol), \
'btype={},era={},czone={}'.format(bldType+1, bldEra+1, climateZone+1)
elif bemid == 'roof_horizontal':
assert refBEM[bldType][bldEra][climateZone].building.FanMax == pytest.approx(matlab_ref_value, abs=tol), \
'btype={},era={},czone={}'.format(bldType+1, bldEra+1, climateZone+1)
matlab_file.close()