def runQsub(self, data, cluster_info):
logging.error("*" * 80)
logging.error("parametersweep.runQsub() modelType={0}".format(
data['modelType']))
logging.error("*" * 80)
modelDb = StochKitModelWrapper.get_by_id(data["modelID"])
path = os.path.abspath(os.path.dirname(__file__))
basedir = path + '/../'
dataDir = tempfile.mkdtemp(dir=basedir + 'output')
job = ParameterSweepJobWrapper()
job.user_id = self.user.user_id()
job.startTime = time.strftime("%Y-%m-%d-%H-%M-%S")
job.name = data["jobName"]
job.inData = json.dumps(data)
job.modelName = modelDb.name
job.outData = dataDir
job.status = "Pending"
job.output_stored = False
try:
templateData = {
"name": modelDb.name,
"modelType": modelDb.type,
"species": modelDb.species,
"parameters": modelDb.parameters,
"reactions": modelDb.reactions,
"speciesSelect": data['speciesSelect'],
"maxTime": data['maxTime'],
"increment": data['increment'],
"trajectories": data['trajectories'],
"seed": data['seed'],
"parameterA": data['parameterA'],
"minValueA": data['minValueA'],
"maxValueA": data['maxValueA'],
"stepsA": data['stepsA'],
"logA": data['logA'],
"parameterB": data['parameterB'],
"minValueB": data['minValueB'],
"maxValueB": data['maxValueB'],
"stepsB": data['stepsB'],
"logB": data['logB'],
"variableCount": data['variableCount'],
"isSpatial": modelDb.isSpatial,
"isLocal": True
}
if modelDb.isSpatial:
try:
meshWrapperDb = mesheditor.MeshWrapper.get_by_id(
modelDb.spatial["mesh_wrapper_id"])
except Exception as e:
logging.exception(e)
logging.error(
"No Mesh file set. Choose one in the Mesh tab of the Model Editor"
)
raise Exception(
"No Mesh file set. Choose one in the Mesh tab of the Model Editor"
)
try:
meshFileObj = fileserver.FileManager.getFile(
self, meshWrapperDb.meshFileId, noFile=False)
templateData["mesh"] = meshFileObj["data"]
except IOError as e:
logging.exception(e)
logging.error("Mesh file inaccessible. Try another mesh")
raise Exception("Mesh file inaccessible. Try another mesh")
templateData[
'reaction_subdomain_assignments'] = modelDb.spatial[
"reactions_subdomain_assignments"]
templateData[
'species_subdomain_assignments'] = modelDb.spatial[
"species_subdomain_assignments"]
templateData[
'species_diffusion_coefficients'] = modelDb.spatial[
"species_diffusion_coefficients"]
templateData['initial_conditions'] = modelDb.spatial[
"initial_conditions"]
templateData['subdomains'] = meshWrapperDb.subdomains
if data['modelType'] == "stochastic":
job.qsubHandle = pickle.dumps(
parametersweep_qsub.stochastic(templateData, cluster_info))
elif data['modelType'] == "deterministic":
job.qsubHandle = pickle.dumps(
parametersweep_qsub.deterministic(templateData,
cluster_info))
elif data['modelType'] == "spatial":
job.qsubHandle = pickle.dumps(
parametersweep_qsub.spatial(templateData, cluster_info))
else:
raise Exception(
"Trying to runQsub on unsupported modelType {0}".format(
data['modelType']))
job.resource = "qsub"
job.put()
except Exception as e:
logging.exception(e)
job.status = 'Failed'
job.delete(self)
raise
return job
def runQsub(self, data, cluster_info):
logging.error("*"*80)
logging.error("parametersweep.runQsub() modelType={0}".format(data['modelType']))
logging.error("*"*80)
modelDb = StochKitModelWrapper.get_by_id(data["modelID"])
path = os.path.abspath(os.path.dirname(__file__))
basedir = path + '/../'
dataDir = tempfile.mkdtemp(dir = basedir + 'output')
job = ParameterSweepJobWrapper()
job.user_id = self.user.user_id()
job.startTime = time.strftime("%Y-%m-%d-%H-%M-%S")
job.name = data["jobName"]
job.inData = json.dumps(data)
job.modelName = modelDb.name
job.outData = dataDir
job.status = "Pending"
job.output_stored = False
try:
templateData = {
"name" : modelDb.name,
"modelType" : modelDb.type,
"species" : modelDb.species,
"parameters" : modelDb.parameters,
"reactions" : modelDb.reactions,
"speciesSelect" : data['speciesSelect'],
"maxTime" : data['maxTime'],
"increment" : data['increment'],
"trajectories" : data['trajectories'],
"seed" : data['seed'],
"parameterA" : data['parameterA'],
"minValueA" : data['minValueA'],
"maxValueA" : data['maxValueA'],
"stepsA" : data['stepsA'],
"logA" : data['logA'],
"parameterB" : data['parameterB'],
"minValueB" : data['minValueB'],
"maxValueB" : data['maxValueB'],
"stepsB" : data['stepsB'],
"logB" : data['logB'],
"variableCount" : data['variableCount'],
"isSpatial" : modelDb.isSpatial,
"isLocal" : True
}
if modelDb.isSpatial:
try:
meshWrapperDb = mesheditor.MeshWrapper.get_by_id(modelDb.spatial["mesh_wrapper_id"])
except Exception as e:
logging.exception(e)
logging.error("No Mesh file set. Choose one in the Mesh tab of the Model Editor")
raise Exception("No Mesh file set. Choose one in the Mesh tab of the Model Editor")
try:
meshFileObj = fileserver.FileManager.getFile(self, meshWrapperDb.meshFileId, noFile = False)
templateData["mesh"] = meshFileObj["data"]
except IOError as e:
logging.exception(e)
logging.error("Mesh file inaccessible. Try another mesh")
raise Exception("Mesh file inaccessible. Try another mesh")
templateData['reaction_subdomain_assignments'] = modelDb.spatial["reactions_subdomain_assignments"]
templateData['species_subdomain_assignments'] = modelDb.spatial["species_subdomain_assignments"]
templateData['species_diffusion_coefficients'] = modelDb.spatial["species_diffusion_coefficients"]
templateData['initial_conditions'] = modelDb.spatial["initial_conditions"]
templateData['subdomains'] = meshWrapperDb.subdomains
if data['modelType'] == "stochastic":
job.qsubHandle = pickle.dumps(parametersweep_qsub.stochastic(templateData, cluster_info))
elif data['modelType'] == "deterministic":
job.qsubHandle = pickle.dumps(parametersweep_qsub.deterministic(templateData, cluster_info))
elif data['modelType'] == "spatial":
job.qsubHandle = pickle.dumps(parametersweep_qsub.spatial(templateData, cluster_info))
else:
raise Exception("Trying to runQsub on unsupported modelType {0}".format(data['modelType']))
job.resource = "qsub"
job.put()
except Exception as e:
logging.exception(e)
job.status='Failed'
job.delete(self)
raise
return job
