From 140d5ea40e738564da6093ee822247a07bfebe2d Mon Sep 17 00:00:00 2001 From: rstiskalek Date: Thu, 19 Sep 2024 13:43:50 +0100 Subject: [PATCH] Add spacing --- csiborgtools/flow/flow_model.py | 22 +++++++++++++++++----- 1 file changed, 17 insertions(+), 5 deletions(-) diff --git a/csiborgtools/flow/flow_model.py b/csiborgtools/flow/flow_model.py index 5000176..6494927 100644 --- a/csiborgtools/flow/flow_model.py +++ b/csiborgtools/flow/flow_model.py @@ -634,11 +634,6 @@ def sample_calibration(Vext_min, Vext_max, Vmono_min, Vmono_max, beta_min, } -############################################################################### -# PV calibration model # -############################################################################### - - def sample_gaussian_hyperprior(param, name, xmin, xmax): """Sample MNR Gaussian hyperprior mean and standard deviation.""" mean = sample(f"{param}_mean_{name}", Uniform(xmin, xmax)) @@ -646,6 +641,11 @@ def sample_gaussian_hyperprior(param, name, xmin, xmax): return mean, std +############################################################################### +# PV calibration model without absolute calibration # +############################################################################### + + class PV_LogLikelihood(BaseFlowValidationModel): """ Peculiar velocity validation model log-likelihood. @@ -953,6 +953,18 @@ class PV_LogLikelihood(BaseFlowValidationModel): return ll0 + jnp.sum(logsumexp(ll, axis=0)) + self.norm +############################################################################### +# PV calibration model with absolute calibration # +############################################################################### + + + + +############################################################################### +# Combining several catalogues # +############################################################################### + + def PV_validation_model(models, distmod_hyperparams_per_model, field_calibration_hyperparams, inference_method): """