Skip to content
Snippets Groups Projects
Commit 73b38d5d authored by Andreas Popp's avatar Andreas Popp
Browse files

Added overview print statements

parent cc5e1580
No related branches found
No related tags found
1 merge request!55Jax port of essential resolve to resolve.re
Pipeline #262371 failed
...@@ -9,10 +9,8 @@ from .likelihood_models import * ...@@ -9,10 +9,8 @@ from .likelihood_models import *
# The classes from .likelihoods model are: # The classes from .likelihoods model are:
# - ModelCalibrationLikelihoodFixedCovariance # - ModelCalibrationLikelihoodFixedCovariance
# - ModelCalibrationLikelihoodVariableCovariance # - ModelCalibrationLikelihoodVariableCovariance
# - ModelImagingLikelihoodFixedCovarianceCalibrationField # - ModelImagingLikelihoodFixedCovariance
# - ModelImagingLikelihoodFixedCovarianceCalibrationOperator # - ModelImagingLikelihoodVariableCovariance
# - ModelImagingLikelihoodVariableCovarianceCalibrationField
# - ModelImagingLikelihoodVariableCovarianceCalibrationOperator
from ..util import _obj2list, _duplicate from ..util import _obj2list, _duplicate
from ..data.observation import Observation from ..data.observation import Observation
...@@ -81,11 +79,15 @@ def CalibrationLikelihood( ...@@ -81,11 +79,15 @@ def CalibrationLikelihood(
flagged_data = jnp.asarray(oo.vis.val)[mask] flagged_data = jnp.asarray(oo.vis.val)[mask]
if log_inv_cov is None: if log_inv_cov is None:
if label is None:
print(f"| Imaging Likelihood {ii} |")
else:
print(f"| Imaging Likelihood {ii} | {label} |")
model = ModelCalibrationLikelihoodFixedCovariance(cop,model_vis,mask) model = ModelCalibrationLikelihoodFixedCovariance(cop,model_vis,mask)
flagged_inv_cov = jnp.asarray(oo.weight.val)[mask] flagged_inv_cov = jnp.asarray(oo.weight.val)[mask]
lh = jft.Gaussian(data=flagged_data, noise_cov_inv=flagged_inv_cov) lh = jft.Gaussian(data=flagged_data, noise_cov_inv=flagged_inv_cov)
else: else:
model = ModelCalibrationLikelihoodVariableCovariance(cop,model_vis,log_inv_cov,mask) model = ModelCalibrationLikelihoodVariableCovariance(cop,model_vis,log_inv_cov,mask)
...@@ -202,12 +204,16 @@ def ImagingLikelihood( ...@@ -202,12 +204,16 @@ def ImagingLikelihood(
flagged_data = jnp.asarray(oo.vis.val)[mask] flagged_data = jnp.asarray(oo.vis.val)[mask]
if log_inv_cov is None: if log_inv_cov is None:
if label is None:
print(f"| Imaging Likelihood {ii} |")
else:
print(f"| Imaging Likelihood {ii} | {label} |")
model = ModelImagingLikelihoodFixedCovariance(R,sky_operator,mask,cop,cfld) model = ModelImagingLikelihoodFixedCovariance(R,sky_operator,mask,cop,cfld)
flagged_inv_cov = jnp.asarray(oo.weight.val)[mask] flagged_inv_cov = jnp.asarray(oo.weight.val)[mask]
lh = jft.Gaussian(data=flagged_data, noise_cov_inv=flagged_inv_cov) lh = jft.Gaussian(data=flagged_data, noise_cov_inv=flagged_inv_cov)
else: else:
model = ModelImagingLikelihoodVariableCovariance(R,sky_operator,log_inv_cov,mask,cop,cfld) model = ModelImagingLikelihoodVariableCovariance(R,sky_operator,log_inv_cov,mask,cop,cfld)
... ...
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please to comment