Aquila-Consortium/csiborgtools
0
0
Fork 0
mirror of https://github.com/Richard-Sti/csiborgtools.git synced 2025-04-17 21:10:55 +00:00
Code Issues Projects Releases Packages Wiki Activity Actions

Add spacing

Browse source
This commit is contained in:
rstiskalek 2024-09-19 13:43:50 +01:00
parent c2bc5a6398
commit 140d5ea40e
1 changed files with 17 additions and 5 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

22
csiborgtools/flow/flow_model.py
View file

@ -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):
"""