def _simulate(spec):
import os
from buildingspy.simulate.Simulator import Simulator
if not spec["simulate"]:
return
wor_dir = create_working_directory()
out_dir = os.path.join(wor_dir, "simulations", spec["name"])
os.makedirs(out_dir)
# Update MODELICAPATH to get the right library version
os.environ["MODELICAPATH"] = ":".join([spec['lib_dir'], out_dir])
# Copy the models
# print("Copying models from {} to {}".format(CWD, wor_dir))
shutil.copytree(os.path.join(CWD, "VAVMultiZone"),
os.path.join(wor_dir, "VAVMultiZone"))
# Change the working directory so that the right checkout is loaded
os.chdir(os.path.join(wor_dir, "VAVMultiZone"))
# Write git information if the simulation is based on a github checkout
if 'git' in spec:
with open(os.path.join(out_dir, "version.txt"), "w+") as text_file:
text_file.write("branch={}\n".format(spec['git']['branch']))
text_file.write("commit={}\n".format(spec['git']['commit']))
s = Simulator(spec["model"], "dymola", outputDirectory=out_dir)
s.addPreProcessingStatement("OutputCPUtime:= true;")
s.addPreProcessingStatement("Advanced.ParallelizeCode = false;")
# s.addPreProcessingStatement("Advanced.EfficientMinorEvents = true;")
if not 'solver' in spec:
s.setSolver("Cvode")
if 'parameters' in spec:
s.addParameters(spec['parameters'])
s.setStartTime(spec["start_time"])
s.setStopTime(spec["stop_time"])
s.setTolerance(1E-5)
s.showGUI(False)
print("Starting simulation in {}".format(out_dir))
s.simulate()
# Copy results back
res_des = os.path.join(CWD, "simulations", spec["name"])
if os.path.isdir(res_des):
shutil.rmtree(res_des)
print("Copying results to {}".format(res_des))
shutil.move(out_dir, res_des)
# Delete the working directory
shutil.rmtree(wor_dir)
def test_addMethods(self):
'''
Tests the various add methods.
'''
import numpy as np
from buildingspy.io.outputfile import Reader
s = Simulator("MyModelicaLibrary.MyModel",
"dymola",
packagePath=self._packagePath)
s.addPreProcessingStatement("Advanced.StoreProtectedVariables:= true;")
s.addPostProcessingStatement(
"Advanced.StoreProtectedVariables:= false;")
s.addModelModifier(
"redeclare Modelica.Blocks.Sources.Step source(offset=-0.1, height=1.1, startTime=0.5)"
)
s.setStartTime(-1)
s.setStopTime(5)
s.setTimeOut(600)
s.setTolerance(1e-4)
s.setSolver("dassl")
s.setNumberOfIntervals(50)
s.setResultFile("myResults")
s.exitSimulator(True)
# s.deleteOutputFiles()
s.showGUI(False)
# s.printModelAndTime()
s.showProgressBar(False)
s.simulate()
# Read the result and test their validity
outDir = s.getOutputDirectory()
resultFile = os.path.abspath(os.path.join(outDir, "myResults.mat"))
r = Reader(resultFile, "dymola")
np.testing.assert_allclose(1.0, r.max('source.y'))
np.testing.assert_allclose(0.725, r.mean('source.y'))
np.testing.assert_allclose(0.725 * 6, r.integral('source.y'))
np.testing.assert_allclose(-0.1, r.min('source.y'))
# Delete output files
s.deleteOutputFiles()
s.deleteLogFiles()
def test_addMethods(self):
'''
Tests the various add methods.
'''
import os
import numpy as np
from buildingspy.io.outputfile import Reader
s = Simulator("MyModelicaLibrary.MyModel", "dymola")
s.addPreProcessingStatement("Advanced.StoreProtectedVariables:= true;")
s.addPostProcessingStatement("Advanced.StoreProtectedVariables:= false;")
s.addModelModifier('redeclare Modelica.Blocks.Sources.Step source(offset=-0.1, height=1.1, startTime=0.5)')
s.setStartTime(-1)
s.setStopTime(5)
s.setTimeOut(600)
s.setTolerance(1e-4)
s.setSolver("dassl")
s.setNumberOfIntervals(50)
s.setResultFile("myResults")
s.exitSimulator(True)
s.deleteOutputFiles()
s.showGUI(False)
# s.printModelAndTime()
s.showProgressBar(False)
s.simulate()
# Read the result and test their validity
outDir = s.getOutputDirectory()
resultFile = os.path.abspath(os.path.join(outDir, "myResults.mat"))
r=Reader(resultFile, "dymola")
np.testing.assert_allclose(1.0, r.max('source.y'))
np.testing.assert_allclose(0.725, r.mean('source.y'))
np.testing.assert_allclose(0.725*6, r.integral('source.y'))
np.testing.assert_allclose(-0.1, r.min('source.y'))
# Delete output files
s.deleteOutputFiles()
s.deleteLogFiles()
