def assemble_results(self, point):
"""
Assemble geodetic data for given point in solution space.
Parameters
----------
point : :func:`pymc3.Point`
Dictionary with model parameters
Returns
-------
List with :class:`heart.GeodeticResult`
"""
logger.debug('Assembling geodetic data ...')
processed_synts = self.get_synthetics(point, outmode='stacked_arrays')
results = []
for i, data in enumerate(self.datasets):
res = data.displacement - processed_synts[i]
results.append(
heart.GeodeticResult(point=point,
processed_obs=data.displacement,
processed_syn=processed_synts[i],
processed_res=res))
return results
def assemble_geodetic_results(self, point):
"""
Assemble geodetic data for given point in solution space.
Parameters
----------
point : :func:`pymc3.Point`
Dictionary with model parameters
Returns
-------
List with :class:`heart.GeodeticResult`
"""
assert self._geodetic_flag
logger.debug('Assembling geodetic data ...')
if self._seismic_flag:
self._seismic_flag = False
reset_flag = True
else:
reset_flag = False
processed_synts = self.get_synthetics(point)['geodetic']
results = []
for i, target in enumerate(self.gtargets):
res = target.displacement - processed_synts[i]
results.append(heart.GeodeticResult(
processed_obs=target.displacement,
processed_syn=processed_synts[i],
processed_res=res))
if reset_flag:
self._seismic_flag = True
return results