def _ref_misfits(self, **outer_misfit_config):
misfits_by_s, misfits_by_sr = seismosizer.make_global_misfits(
self.ref_misfits_by_src,
self.ref_norms_by_src,
receiver_mask=self.receiver_mask,
**outer_misfit_config)
return misfits_by_s[0], misfits_by_sr[0, :]
def _best_source(self, return_misfits_by_r=False, **outer_misfit_config):
misfits_by_s, misfits_by_sr = seismosizer.make_global_misfits(
self.misfits_by_src, self.norms_by_src,
receiver_mask=self.receiver_mask,
**outer_misfit_config )
ibest = num.nanargmin(misfits_by_s)
if not num.isfinite(ibest):
ibest = 0
if not return_misfits_by_r:
return self.sources[ibest], misfits_by_s
else:
# misfit variability by receiver
misfits_varia_by_r = num.std(misfits_by_sr,0)
return self.sources[ibest], misfits_by_s, misfits_by_sr[ibest,:], misfits_varia_by_r
def _best_source(self, return_misfits_by_r=False, **outer_misfit_config):
misfits_by_s, misfits_by_sr = seismosizer.make_global_misfits(
self.misfits_by_src,
self.norms_by_src,
receiver_mask=self.receiver_mask,
**outer_misfit_config)
ibest = num.nanargmin(misfits_by_s)
if not num.isfinite(ibest):
ibest = 0
if not return_misfits_by_r:
return self.sources[ibest], misfits_by_s
else:
# misfit variability by receiver
misfits_varia_by_r = num.std(misfits_by_sr, 0)
return self.sources[ibest], misfits_by_s, misfits_by_sr[
ibest, :], misfits_varia_by_r
def _ref_misfits(self, **outer_misfit_config):
misfits_by_s, misfits_by_sr = seismosizer.make_global_misfits(
self.ref_misfits_by_src, self.ref_norms_by_src,
receiver_mask=self.receiver_mask, **outer_misfit_config)
return misfits_by_s[0], misfits_by_sr[0,:]
