Skip to content

Add a fisher forecast mode for pv #26

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
corentinravoux merged 13 commits into from
Jun 25, 2024
Merged

Add a fisher forecast mode for pv #26

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
13 commits
Select commit Hold shift + click to select a range
51c3eef
Skeleton for fisher matrix calculation
corentinravoux May 2, 2024
ef054cd
Merge branch 'main' into add-a-fisher-forecast-mode-for-pv
corentinravoux May 16, 2024
47e9a4b
moving and improvement of fisher class
corentinravoux May 16, 2024
087f334
large modification of vector loading for fisher, include scatter mode
corentinravoux May 16, 2024
7c8174d
working implementation of fisher matrix calculation
corentinravoux May 16, 2024
f86e967
add a script
corentinravoux May 16, 2024
f53dc00
adding symbolic generation of fisher terms
corentinravoux May 17, 2024
8cdb2ea
Creation of fisher terms for all models
corentinravoux May 17, 2024
6b1af8e
Merge pull request #28 from corentinravoux/main
corentinravoux Jun 3, 2024
0705a0c
correct bug in velocity error from saltfit
DamianoRosselli Jun 10, 2024
a5f6d92
fix on the script
corentinravoux Jun 18, 2024
e53779d
important fix for the fisher forecast
corentinravoux Jun 25, 2024
5d1a6df
return values dictionnary at the end of the minuit run
corentinravoux Jun 25, 2024
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
3 changes: 2 additions & 1 deletion flip/__init__.py
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -1,7 +1,8 @@
"""Init file of the flip package."""

import os

from . import covariance, fitter, gridding, likelihood, utils, power_spectra
from . import covariance, fisher, fitter, gridding, likelihood, power_spectra, utils

__version__ = "1.0.0"
__flip_dir_path__ = os.path.dirname(__file__)
299 changes: 299 additions & 0 deletions flip/covariance/adamsblake17plane/fisher_terms.py
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,299 @@
import numpy as np


def get_partial_derivative_coefficients(
model_type,
parameter_values_dict,
variant=None,
):
if model_type == "density":
return get_partial_derivative_coefficients_density(
parameter_values_dict,
variant=variant,
)
elif model_type == "velocity":
return get_partial_derivative_coefficients_velocity(
parameter_values_dict,
variant=variant,
)
elif model_type == "density_velocity":
return get_partial_derivative_coefficients_density_velocity(
parameter_values_dict,
variant=variant,
)
elif model_type == "full":
return get_partial_derivative_coefficients_full(
parameter_values_dict,
variant=variant,
)


def get_partial_derivative_coefficients_velocity(
parameter_values_dict,
variant=None,
):
if variant == "growth_index":
partial_coefficients_dict = {
"Omegam": {
"vv": [
2
* parameter_values_dict["Omegam"]
** (2 * parameter_values_dict["gamma"])
* parameter_values_dict["gamma"]
* parameter_values_dict["s8"] ** 2
/ parameter_values_dict["Omegam"],
],
},
"gamma": {
"vv": [
2
* parameter_values_dict["Omegam"]
** (2 * parameter_values_dict["gamma"])
* parameter_values_dict["s8"] ** 2
* np.log(parameter_values_dict["Omegam"]),
],
},
"s8": {
"vv": [
2
* parameter_values_dict["Omegam"]
** (2 * parameter_values_dict["gamma"])
* parameter_values_dict["s8"],
],
},
"bs8": {
"vv": [
0,
],
},
}
else:
partial_coefficients_dict = {
"fs8": {
"vv": [
2 * parameter_values_dict["fs8"],
],
},
"bs8": {
"vv": [
0,
],
},
}
return partial_coefficients_dict


def get_partial_derivative_coefficients_density(
parameter_values_dict,
variant=None,
):
if variant == "growth_index":
partial_coefficients_dict = {
"Omegam": {
"gg": [
0,
],
},
"gamma": {
"gg": [
0,
],
},
"s8": {
"gg": [
0,
],
},
"bs8": {
"gg": [
2 * parameter_values_dict["bs8"],
],
},
}
else:
partial_coefficients_dict = {
"fs8": {
"gg": [
0,
],
},
"bs8": {
"gg": [
2 * parameter_values_dict["bs8"],
],
},
}
return partial_coefficients_dict


def get_partial_derivative_coefficients_density_velocity(
parameter_values_dict,
variant=None,
):
if variant == "growth_index":
partial_coefficients_dict = {
"Omegam": {
"gg": [
0,
],
"vv": [
2
* parameter_values_dict["Omegam"]
** (2 * parameter_values_dict["gamma"])
* parameter_values_dict["gamma"]
* parameter_values_dict["s8"] ** 2
/ parameter_values_dict["Omegam"],
],
},
"gamma": {
"gg": [
0,
],
"vv": [
2
* parameter_values_dict["Omegam"]
** (2 * parameter_values_dict["gamma"])
* parameter_values_dict["s8"] ** 2
* np.log(parameter_values_dict["Omegam"]),
],
},
"s8": {
"gg": [
0,
],
"vv": [
2
* parameter_values_dict["Omegam"]
** (2 * parameter_values_dict["gamma"])
* parameter_values_dict["s8"],
],
},
"bs8": {
"gg": [
2 * parameter_values_dict["bs8"],
],
"vv": [
0,
],
},
}
else:
partial_coefficients_dict = {
"fs8": {
"gg": [
0,
],
"vv": [
2 * parameter_values_dict["fs8"],
],
},
"bs8": {
"gg": [
2 * parameter_values_dict["bs8"],
],
"vv": [
0,
],
},
}
return partial_coefficients_dict


def get_partial_derivative_coefficients_full(
parameter_values_dict,
variant=None,
):
if variant == "growth_index":
partial_coefficients_dict = {
"Omegam": {
"gg": [
0,
],
"gv": [
parameter_values_dict["Omegam"] ** parameter_values_dict["gamma"]
* parameter_values_dict["bs8"]
* parameter_values_dict["gamma"]
* parameter_values_dict["s8"]
/ parameter_values_dict["Omegam"],
],
"vv": [
2
* parameter_values_dict["Omegam"]
** (2 * parameter_values_dict["gamma"])
* parameter_values_dict["gamma"]
* parameter_values_dict["s8"] ** 2
/ parameter_values_dict["Omegam"],
],
},
"gamma": {
"gg": [
0,
],
"gv": [
parameter_values_dict["Omegam"] ** parameter_values_dict["gamma"]
* parameter_values_dict["bs8"]
* parameter_values_dict["s8"]
* np.log(parameter_values_dict["Omegam"]),
],
"vv": [
2
* parameter_values_dict["Omegam"]
** (2 * parameter_values_dict["gamma"])
* parameter_values_dict["s8"] ** 2
* np.log(parameter_values_dict["Omegam"]),
],
},
"s8": {
"gg": [
0,
],
"gv": [
parameter_values_dict["Omegam"] ** parameter_values_dict["gamma"]
* parameter_values_dict["bs8"],
],
"vv": [
2
* parameter_values_dict["Omegam"]
** (2 * parameter_values_dict["gamma"])
* parameter_values_dict["s8"],
],
},
"bs8": {
"gg": [
2 * parameter_values_dict["bs8"],
],
"gv": [
parameter_values_dict["Omegam"] ** parameter_values_dict["gamma"]
* parameter_values_dict["s8"],
],
"vv": [
0,
],
},
}
else:
partial_coefficients_dict = {
"fs8": {
"gg": [
0,
],
"gv": [
parameter_values_dict["bs8"],
],
"vv": [
2 * parameter_values_dict["fs8"],
],
},
"bs8": {
"gg": [
2 * parameter_values_dict["bs8"],
],
"gv": [
parameter_values_dict["fs8"],
],
"vv": [
0,
],
},
}
return partial_coefficients_dict
Loading