def post_execute(self, result):
"""
:param result:
a dictionary (trt_model_id, gsim) -> haz_curves or an empty
dictionary if hazard_curves_from_gmfs is false
"""
oq = self.oqparam
if not oq.hazard_curves_from_gmfs and not oq.ground_motion_fields:
return
if oq.hazard_curves_from_gmfs:
ClassicalCalculator.post_execute.__func__(self, result)
if oq.compare_with_classical: # compute classical curves
export_dir = os.path.join(oq.export_dir, 'cl')
if not os.path.exists(export_dir):
os.makedirs(export_dir)
oq.export_dir = export_dir
# use a different datastore
self.cl = ClassicalCalculator(oq, self.monitor)
self.cl.datastore.parent = self.datastore
result = self.cl.run(pre_execute=False, clean_up=False)
for imt in self.mean_curves.dtype.fields:
rdiff, index = max_rel_diff_index(self.cl.mean_curves[imt],
self.mean_curves[imt])
logging.warn(
'Relative difference with the classical '
'mean curves for IMT=%s: %d%% at site index %d', imt,
rdiff * 100, index)
