PyCosmo 2.1.3__tar.gz → 2.2.1__tar.gz

{PyCosmo-2.1.3 → pycosmo-2.2.1}/AUTHORS.rst

@@ -19,4 +19,3 @@ Contact information
 
 If you have any suggestions or questions about PyCosmo feel free to email us
 at pycosmo@lists.phys.ethz.ch.
-

pycosmo-2.2.1/MANIFEST.in

@@ -0,0 +1 @@
+prune tests

pycosmo-2.2.1/PKG-INFO

@@ -0,0 +1,55 @@
+Metadata-Version: 2.4
+Name: PyCosmo
+Version: 2.2.1
+Summary: A multi-purpose cosmology calculation tool
+Author-email: Uwe Schmitt et al <uwe.schmitt@id.ethz.ch>
+License: GPLv3
+Project-URL: Source, https://cosmo-docs.phys.ethz.ch/PyCosmo
+Classifier: Development Status :: 4 - Beta
+Classifier: Intended Audience :: Developers
+Classifier: Intended Audience :: Science/Research
+Classifier: License :: OSI Approved :: GNU General Public License v3 (GPLv3)
+Classifier: Natural Language :: English
+Classifier: Operating System :: OS Independent
+Classifier: Programming Language :: C
+Classifier: Programming Language :: Python :: 3
+Classifier: Programming Language :: Python :: 3.12
+Classifier: Programming Language :: Python :: 3.11
+Classifier: Programming Language :: Python :: 3.10
+Classifier: Programming Language :: Python :: 3.9
+Classifier: Topic :: Scientific/Engineering :: Astronomy
+Requires-Python: >=3.10
+Description-Content-Type: text/markdown; charset=UTF-8
+License-File: LICENSE
+License-File: AUTHORS.rst
+Requires-Dist: numpy>=2
+Requires-Dist: recfast4py>=0.2.3
+Requires-Dist: dill
+Requires-Dist: matplotlib
+Requires-Dist: numba>=0.60
+Requires-Dist: sympy2c>=2.3.1
+Requires-Dist: numexpr>=2.10
+Requires-Dist: scipy
+Requires-Dist: setuptools>=75.6.0
+Requires-Dist: mpmath==1.3
+Dynamic: license-file
+
+# PyCosmo
+
+PyCosmo is a multi-purpose cosmology calculation tool in Python. It is
+designed to be interactive and user friendly.
+
+The code includes calculations of:
+  - Background
+  - Linear perturbations (incl. fitting functions and a Boltzmann solver)
+  - Non-linear perturbations
+  - Series of observables (e.g., weak lensing and galaxy clustering
+    power spectrum, ...)
+  - Several model extensions (massive neutrinos, wCDM, ...)
+
+
+The code documentation can be found at
+[https://cosmo-docs.phys.ethz.ch/PyCosmo](https://cosmo-docs.phys.ethz.ch/PyCosmo)
+
+Interactive tutorial notebooks are available on the
+[PyCosmoHUB](https://pycosmohub.com).

pycosmo-2.2.1/README.md

@@ -0,0 +1,19 @@
+# PyCosmo
+
+PyCosmo is a multi-purpose cosmology calculation tool in Python. It is
+designed to be interactive and user friendly.
+
+The code includes calculations of:
+  - Background
+  - Linear perturbations (incl. fitting functions and a Boltzmann solver)
+  - Non-linear perturbations
+  - Series of observables (e.g., weak lensing and galaxy clustering
+    power spectrum, ...)
+  - Several model extensions (massive neutrinos, wCDM, ...)
+
+
+The code documentation can be found at
+[https://cosmo-docs.phys.ethz.ch/PyCosmo](https://cosmo-docs.phys.ethz.ch/PyCosmo)
+
+Interactive tutorial notebooks are available on the
+[PyCosmoHUB](https://pycosmohub.com).

pycosmo-2.2.1/pyproject.toml

@@ -0,0 +1,113 @@
+[build-system]
+requires = ["setuptools >= 64", "numpy", "Cython"]
+build-backend = "setuptools.build_meta"
+
+
+[project]
+name = "PyCosmo"
+version = "2.2.1"
+license = {text = "GPLv3"}
+description = "A multi-purpose cosmology calculation tool"
+authors = [{name = "Uwe Schmitt et al", email = "uwe.schmitt@id.ethz.ch"}]
+classifiers = [
+    "Development Status :: 4 - Beta",
+    "Intended Audience :: Developers",
+    "Intended Audience :: Science/Research",
+    "License :: OSI Approved :: GNU General Public License v3 (GPLv3)",
+    "Natural Language :: English",
+    "Operating System :: OS Independent",
+    "Programming Language :: C",
+    "Programming Language :: Python :: 3",
+    "Programming Language :: Python :: 3.12",
+    "Programming Language :: Python :: 3.11",
+    "Programming Language :: Python :: 3.10",
+    "Programming Language :: Python :: 3.9",
+    "Topic :: Scientific/Engineering :: Astronomy",
+]
+
+urls = {Source = "https://cosmo-docs.phys.ethz.ch/PyCosmo"}
+requires-python = ">=3.10"
+dependencies = [
+    "numpy>=2",
+    "recfast4py>=0.2.3",
+    "dill",
+    "matplotlib",
+    "numba>=0.60",
+    "sympy2c>=2.3.1",
+    "numexpr>=2.10",
+    "scipy",
+    "setuptools>=75.6.0",
+    "mpmath==1.3",
+]
+
+[project.readme]
+file = "README.md"
+content-type = "text/markdown; charset=UTF-8"
+
+[dependency-groups]
+dev = [
+    "Cython",
+    "twine",
+    "pytest",
+    "pytest-cov",
+    "pytest-regtest>=2.2.0a2",
+    "ruff",
+    "isort",
+    "pre-commit",
+    "pdbp>=1.7.1",
+]
+
+docs = [
+    "sphinx < 7",
+    "sphinx-pynpoint-theme",
+    "myst-parser",
+]
+
+[project.scripts]
+recompile = "PyCosmo:Cosmo.recompile_from_cli"
+
+[tool.setuptools]
+zip-safe = false
+include-package-data = true
+
+[tool.setuptools.package-data]
+"*" = [
+    "*.pyx",
+    "*.c",
+    "*.h",
+    "config/*.ini",
+]
+
+
+[tool.ruff]
+line-length    = 88
+extend-exclude = ["benchmarks/", "*.ipynb"]
+
+    [tool.ruff.format]
+    quote-style           = "double"
+    line-ending           = "lf"
+    docstring-code-format = false
+
+    [tool.ruff.lint]
+    extend-select = [
+        "RUF100",             # unnecessary `noqa` directives
+    ]
+    ignore = [
+        "E501",
+        "E731",               # lambda assignment
+        "E741",               # ambiguous variable name
+    ]
+
+    [tool.ruff.lint.per-file-ignores]
+    '__init__.py' = [
+        "F401",               # imported but unused
+    ]
+
+
+[tool.pytest.ini_options]
+testpaths = ["tests"]
+addopts = ["--verbose", "--durations=0", "--durations-min=5"]
+markers = [
+    "slow: Marks tests as slow to run. Pass `--runslow` to include them."
+]
+console_output_style = "count"

pycosmo-2.2.1/setup.cfg

@@ -0,0 +1,4 @@
+[egg_info]
+tag_build = 
+tag_date = 0
+

pycosmo-2.2.1/setup.py

@@ -0,0 +1,41 @@
+#!/usr/bin/env python
+import numpy
+from Cython.Build import cythonize
+from setuptools import Extension, setup
+
+files = [
+    "const.c",
+    "main.c",
+    "polevl.c",
+    "sici.c",
+    "sicif.c",
+    "polevlf.c",
+    "logf.c",
+    "sinf.c",
+    "constf.c",
+    "mtherr.c",
+]
+ext_modules = []
+
+
+for stem in (
+    "halo_integral",
+    "halo_integral_1h",
+    "halo_integral_2h",
+    "halo_integral_HI_1h",
+    "halo_integral_HI_2h",
+    "rho_HI_integral",
+    "rho_m_integral",
+):
+    ext_modules.append(
+        Extension(
+            name=f"PyCosmo.cython.{stem}",
+            ext_modules=cythonize(f"src/PyCosmo/cython/{stem}.pyx"),
+            sources=["src/PyCosmo/cython/" + file for file in files + [stem + ".c"]],
+            include_dirs=[numpy.get_include()],
+        )
+    )
+
+setup(
+    ext_modules=ext_modules,
+)

{PyCosmo-2.1.3 → pycosmo-2.2.1}/src/PyCosmo/BoltzmannSolver/Fields.py

@@ -65,7 +65,6 @@ class Fields(object):
         return super().__dir__() + list(map(str, self._symbols))
 
     def __getattr__(self, symbol):
-
         # needed to support pickling:
         if "_offsets" not in vars(self):
             raise AttributeError()

{PyCosmo-2.1.3 → pycosmo-2.2.1}/src/PyCosmo/BoltzmannSolver/LinearPerturbationBoltzmann.py

@@ -16,7 +16,6 @@
 
 import itertools
 import warnings
-from datetime import datetime
 from functools import lru_cache
 from multiprocessing import current_process
 
@@ -26,8 +25,9 @@ from ..LinearPerturbationBase import LinearPerturbationBase
 from .Fields import Fields
 from .LsodaSolver import LsodaSolver
 
+
 # use decorator to print messages only once per model:
-@lru_cache(maxsize=1)  # noqa: E302
+@lru_cache(maxsize=1)
 def print_new_traces(cache_file):
     print()
     print("new traces detected! you might want to run")
@@ -38,36 +38,6 @@ def print_new_traces(cache_file):
 proc_name = current_process().name
 
 
-def trace_call(frame, event, arg):
-    if proc_name.endswith("-2"):
-        return trace_call
-    if event in ("line",):
-        return trace_call
-    fname = frame.f_code.co_filename
-    if any(
-        n in fname
-        for n in (
-            "pprint",
-            "dill",
-            "pickle",
-            "numba",
-            "logging",
-            "ini_handling",
-            "llvmlite",
-            "traceback",
-            "weakref",
-        )
-    ):
-        return trace_call
-    print(datetime.now(), end=" ")
-    fr = frame
-    print(proc_name, f"{event:8s}", end="")
-    while (fr := fr.f_back) is not None:
-        print("  ", end="")
-    print(fname, frame.f_code.co_name, frame.f_lineno, flush=True)
-    return trace_call
-
-
 class LinearPerturbationBoltzmann(LinearPerturbationBase):
     """
     Class for computing linear perturbations by solving the Einstein-Boltzmann ODE
@@ -196,9 +166,7 @@ class LinearPerturbationBoltzmann(LinearPerturbationBase):
                 pk = norm * growth**2 * k**self._params.n * T_k**2
             else:
                 # pk = norm * np.outer(growth.T ** 2, k ** self._params.n * T_k ** 2).T
-                pk = (
-                    norm * growth**2 * (k**self._params.n * T_k**2).reshape(-1, 1)
-                )
+                pk = norm * growth**2 * (k**self._params.n * T_k**2).reshape(-1, 1)
 
         if (self._params.pk_norm_type == "A_s") or (
             self._params.pk_norm_type == "sigma8"
@@ -363,9 +331,9 @@ class LinearPerturbationBoltzmann(LinearPerturbationBase):
 
         :return: D(a): growth factor normalised to 1 at a=1 [1]
         """
-        assert (
-            k is not None
-        ), "growth factor of Einstein-Boltzmann model requires a k value"
+        assert k is not None, (
+            "growth factor of Einstein-Boltzmann model requires a k value"
+        )
         # TODO: add option where D is normalised to a in matter dominated era
         grid = np.log(np.atleast_1d(a))
         if grid[-1] != 0.0:

{PyCosmo-2.1.3 → pycosmo-2.2.1}/src/PyCosmo/BoltzmannSolver/LsodaSolver.py

@@ -21,7 +21,6 @@ import numpy as np
 
 
 class LsodaSolver:
-
     SOLVER_ID = "linear_perturbation"
 
     _wrappers = {}
@@ -89,7 +88,6 @@ class LsodaSolver:
         enable_sparse_lu_solver=True,
         keep_lna0=False,
     ):
-
         self._setup_globals(k)
         self.wrapper.set_sec_factor(sec_factor)
 

{PyCosmo-2.1.3 → pycosmo-2.2.1}/src/PyCosmo/BoltzmannSolver/ScipySolver.py

@@ -24,7 +24,6 @@ from .Fields import Fields
 
 
 class ScipySolver:
-
     _wrappers = {}
 
     def __init__(
@@ -301,7 +300,7 @@ def merge(wrapper, t_0, t_switch, t_1, y_0, y_1):  # pragma: no cover
         + a
         * H_Mpc
         * wrapper.taudot_interp_ufunc(a)
-        / k**2
+        / k** 2
         * (
             u_b
             - k / (a * H_Mpc) * wrapper.c_s_interp_ufunc(a) ** 2 * delta_b

{PyCosmo-2.1.3 → pycosmo-2.2.1}/src/PyCosmo/CosmologyCore.py

@@ -16,6 +16,8 @@
 
 # pragma: no cover
 
+# ruff: noqa
+
 import sympy as sp
 from sympy2c import (
     ERROR,
@@ -141,9 +143,7 @@ omega_l_a = IfThenElse(flat_universe, (1 - omega_m - omega_r), omega_l)
 omega_k = 1 - omega_m - omega_r - omega_l_a
 
 # Friedmann equation
-H = H0 * (omega_r * a**-4 + omega_m * a**-3 + omega_k * a**-2 + omega_l_a) ** (
-    1 / 2
-)
+H = H0 * (omega_r * a**-4 + omega_m * a**-3 + omega_k * a**-2 + omega_l_a) ** (1 / 2)
 
 R = (3 / 4) * omega_b / omega_gamma * a
 
@@ -522,14 +522,10 @@ def rhs(H=H_Mpc, H0=H0_Mpc, R=R, l_max=l_max):
     dTheta_dlan[0] = -k_H * Theta[1] - dPhi_dlan  # 0th photon multipole
     dTheta_dlan[1] = k_H / 3 * (Theta[0] - 2 * Theta[2] + Psi) + taudot_interp(a) / (
         a * H
-    ) * (
-        Theta[1] - u_b / 3
-    )  # 1st photon multipole
+    ) * (Theta[1] - u_b / 3)  # 1st photon multipole
     dTheta_dlan[2] = k_H / 5 * (2 * Theta[1] - 3 * Theta[3]) + taudot_interp(a) / (
         a * H
-    ) * (
-        Theta[2] - Pi / 10
-    )  # 2nd photon multipole
+    ) * (Theta[2] - Pi / 10)  # 2nd photon multipole
 
     dTheta_P_dlan[0] = k_H / 1 * -Theta_P[1] + taudot_interp(a) / (a * H) * (
         Theta_P[0] - Pi / 2
@@ -540,9 +536,7 @@ def rhs(H=H_Mpc, H0=H0_Mpc, R=R, l_max=l_max):
     )  # 1st photon polarization multipole
     dTheta_P_dlan[2] = k_H / 5 * (2 * Theta_P[1] - 3 * Theta_P[3]) + taudot_interp(
         a
-    ) / (a * H) * (
-        Theta_P[2] - Pi / 10
-    )  # 2nd photon polarization multipole
+    ) / (a * H) * (Theta_P[2] - Pi / 10)  # 2nd photon polarization multipole
 
     for l in range(3, l_max):  # lth photons multipoles
         dTheta_dlan[l] = (

{PyCosmo-2.1.3 → pycosmo-2.2.1}/src/PyCosmo/CosmologyCore_massivenu.py

@@ -15,6 +15,9 @@
 # program.  If not, see <http://www.gnu.org/licenses/>.
 
 
+# pragma: no cover
+# ruff: noqa
+
 import numpy as np
 import sympy as sp
 from sympy2c import (
@@ -121,6 +124,7 @@ globals_ = Globals(
     T_mnu,
 )
 
+
 # check functions used to impose conditions
 # on the parameters
 def check_massive_nu_total_mass(params):
@@ -767,14 +771,10 @@ def rhs(H=H_Mpc, H0=H0_Mpc, amnu=amnu):
     dTheta_dlan[0] = -k_H * Theta[1] - dPhi_dlan  # 0th photon multipole
     dTheta_dlan[1] = k_H / 3 * (Theta[0] - 2 * Theta[2] + Psi) + taudot_interp(a) / (
         a * H
-    ) * (
-        Theta[1] - u_b / 3
-    )  # 1st photon multipole
+    ) * (Theta[1] - u_b / 3)  # 1st photon multipole
     dTheta_dlan[2] = k_H / 5 * (2 * Theta[1] - 3 * Theta[3]) + taudot_interp(a) / (
         a * H
-    ) * (
-        Theta[2] - Pi / 10
-    )  # 2nd photon multipole
+    ) * (Theta[2] - Pi / 10)  # 2nd photon multipole
 
     dTheta_P_dlan[0] = k_H / 1 * -Theta_P[1] + taudot_interp(a) / (a * H) * (
         Theta_P[0] - Pi / 2
@@ -785,9 +785,7 @@ def rhs(H=H_Mpc, H0=H0_Mpc, amnu=amnu):
     )  # 1st photon polarization multipole
     dTheta_P_dlan[2] = k_H / 5 * (2 * Theta_P[1] - 3 * Theta_P[3]) + taudot_interp(
         a
-    ) / (a * H) * (
-        Theta_P[2] - Pi / 10
-    )  # 2nd photon polarization multipole
+    ) / (a * H) * (Theta_P[2] - Pi / 10)  # 2nd photon polarization multipole
 
     for l in range(3, l_max):  # lth photons multipoles
         dTheta_dlan[l] = (
@@ -845,7 +843,6 @@ def rhs(H=H_Mpc, H0=H0_Mpc, amnu=amnu):
 
     for m in range(2, l_max_mnu):
         for j in range(nq):
-
             dM_dlan[m * nq + j] = (
                 qq[j]
                 * k
@@ -862,10 +859,7 @@ def rhs(H=H_Mpc, H0=H0_Mpc, amnu=amnu):
             1
             / (a * H)
             * (
-                qq[j]
-                * k
-                / sp.sqrt(amnu**2 + qq[j] ** 2)
-                * M[(l_max_mnu - 1) * nq + j]
+                qq[j] * k / sp.sqrt(amnu**2 + qq[j] ** 2) * M[(l_max_mnu - 1) * nq + j]
                 - (l_max_mnu + 1) / eta_interp(a) * M[l_max_mnu * nq + j]
             )
         )

{PyCosmo-2.1.3 → pycosmo-2.2.1}/src/PyCosmo/CosmologyCore_massivenu_rsa.py

@@ -14,6 +14,10 @@
 # You should have received a copy of the GNU General Public License along with this
 # program.  If not, see <http://www.gnu.org/licenses/>.
 
+
+# pragma: no cover
+# ruff: noqa
+
 from CosmologyCore_massivenu import *
 from sympy2c import OdeCombined
 
@@ -140,7 +144,6 @@ def rhs_rsa(H=H_Mpc, H0=H0_Mpc, R=R, amnu=amnu, l_max=2):
     # lth multipoles massive neutrinos
     for m in range(2, l_max_mnu):
         for j in range(nq):
-
             dM_dlan[m * nq + j] = (
                 qq[j]
                 * k
@@ -153,10 +156,7 @@ def rhs_rsa(H=H_Mpc, H0=H0_Mpc, R=R, amnu=amnu, l_max=2):
             1
             / (a * H)
             * (
-                qq[j]
-                * k
-                / sp.sqrt(amnu**2 + qq[j] ** 2)
-                * M[(l_max_mnu - 1) * nq + j]
+                qq[j] * k / sp.sqrt(amnu**2 + qq[j] ** 2) * M[(l_max_mnu - 1) * nq + j]
                 - (l_max_mnu + 1) / eta_interp(a) * M[l_max_mnu * nq + j]
             )
         )
@@ -232,13 +232,7 @@ def merge(wrapper, t_0, t_switch, t_1, y_0, y_1):  # pragma: no cover
     T_0 = Tcmb * T_mnu
     amnu = a * massive_nu_total_mass / (N_massive_nu * kb * T_0)
     rho_ev = (
-        3
-        * H0**2
-        * (hbar) ** 3
-        * c_ms**3
-        / (8 * np.pi * G)
-        / evc2
-        / (mpc / 1000) ** 2
+        3 * H0**2 * (hbar) ** 3 * c_ms**3 / (8 * np.pi * G) / evc2 / (mpc / 1000) ** 2
     ) / (kb * T_0) ** 4
     C_nu = N_massive_nu / (2 * np.pi) ** 2 / rho_ev
 

{PyCosmo-2.1.3 → pycosmo-2.2.1}/src/PyCosmo/CosmologyCore_massivenu_wcdm.py

@@ -14,6 +14,8 @@
 # You should have received a copy of the GNU General Public License along with this
 # program.  If not, see <http://www.gnu.org/licenses/>.
 
+# pragma: no cover
+# ruff: noqa
 
 import numpy as np
 import sympy as sp
@@ -129,6 +131,7 @@ globals_ = Globals(
     cs_de2,
 )
 
+
 # check functions used to impose conditions
 # on the parameters
 def check_massive_nu_total_mass(params):
@@ -816,14 +819,10 @@ def rhs(H=H_Mpc, H0=H0_Mpc, R=R, amnu=amnu):
     dTheta_dlan[0] = -k_H * Theta[1] - dPhi_dlan  # 0th photon multipole
     dTheta_dlan[1] = k_H / 3 * (Theta[0] - 2 * Theta[2] + Psi) + taudot_interp(a) / (
         a * H
-    ) * (
-        Theta[1] - u_b / 3
-    )  # 1st photon multipole
+    ) * (Theta[1] - u_b / 3)  # 1st photon multipole
     dTheta_dlan[2] = k_H / 5 * (2 * Theta[1] - 3 * Theta[3]) + taudot_interp(a) / (
         a * H
-    ) * (
-        Theta[2] - Pi / 10
-    )  # 2nd photon multipole
+    ) * (Theta[2] - Pi / 10)  # 2nd photon multipole
 
     dTheta_P_dlan[0] = k_H / 1 * -Theta_P[1] + taudot_interp(a) / (a * H) * (
         Theta_P[0] - Pi / 2
@@ -834,9 +833,7 @@ def rhs(H=H_Mpc, H0=H0_Mpc, R=R, amnu=amnu):
     )  # 1st photon polarization multipole
     dTheta_P_dlan[2] = k_H / 5 * (2 * Theta_P[1] - 3 * Theta_P[3]) + taudot_interp(
         a
-    ) / (a * H) * (
-        Theta_P[2] - Pi / 10
-    )  # 2nd photon polarization multipole
+    ) / (a * H) * (Theta_P[2] - Pi / 10)  # 2nd photon polarization multipole
 
     for l in range(3, l_max):  # lth photons multipoles
         dTheta_dlan[l] = (
@@ -894,7 +891,6 @@ def rhs(H=H_Mpc, H0=H0_Mpc, R=R, amnu=amnu):
 
     for m in range(2, l_max_mnu):
         for j in range(nq):
-
             dM_dlan[m * nq + j] = (
                 qq[j]
                 * k
@@ -911,10 +907,7 @@ def rhs(H=H_Mpc, H0=H0_Mpc, R=R, amnu=amnu):
             1
             / (a * H)
             * (
-                qq[j]
-                * k
-                / sp.sqrt(amnu**2 + qq[j] ** 2)
-                * M[(l_max_mnu - 1) * nq + j]
+                qq[j] * k / sp.sqrt(amnu**2 + qq[j] ** 2) * M[(l_max_mnu - 1) * nq + j]
                 - (l_max_mnu + 1) / eta_interp(a) * M[l_max_mnu * nq + j]
             )
         )

{PyCosmo-2.1.3 → pycosmo-2.2.1}/src/PyCosmo/CosmologyCore_massivenu_wcdm_rsa.py

@@ -14,6 +14,9 @@
 # You should have received a copy of the GNU General Public License along with this
 # program.  If not, see <http://www.gnu.org/licenses/>.
 
+# pragma: no cover
+# ruff: noqa
+
 from CosmologyCore_massivenu_wcdm import *
 from sympy2c import OdeCombined
 
@@ -157,7 +160,6 @@ def rhs_rsa(H=H_Mpc, H0=H0_Mpc, R=R, amnu=amnu, l_max=2):
     # lth multipoles massive neutrinos
     for m in range(2, l_max_mnu):
         for j in range(nq):
-
             dM_dlan[m * nq + j] = (
                 qq[j]
                 * k
@@ -170,10 +172,7 @@ def rhs_rsa(H=H_Mpc, H0=H0_Mpc, R=R, amnu=amnu, l_max=2):
             1
             / (a * H)
             * (
-                qq[j]
-                * k
-                / sp.sqrt(amnu**2 + qq[j] ** 2)
-                * M[(l_max_mnu - 1) * nq + j]
+                qq[j] * k / sp.sqrt(amnu**2 + qq[j] ** 2) * M[(l_max_mnu - 1) * nq + j]
                 - (l_max_mnu + 1) / eta_interp(a) * M[l_max_mnu * nq + j]
             )
         )
@@ -256,13 +255,7 @@ def merge(wrapper, t_0, t_switch, t_1, y_0, y_1):  # pragma: no cover
     T_0 = Tcmb * T_mnu
     amnu = a * massive_nu_total_mass / (N_massive_nu * kb * T_0)
     rho_ev = (
-        3
-        * H0**2
-        * (hbar) ** 3
-        * c_ms**3
-        / (8 * np.pi * G)
-        / evc2
-        / (mpc / 1000) ** 2
+        3 * H0**2 * (hbar) ** 3 * c_ms**3 / (8 * np.pi * G) / evc2 / (mpc / 1000) ** 2
     ) / (kb * T_0) ** 4
     C_nu = N_massive_nu / (2 * np.pi) ** 2 / rho_ev
 

{PyCosmo-2.1.3 → pycosmo-2.2.1}/src/PyCosmo/CosmologyCore_rsa.py

@@ -15,6 +15,7 @@
 # program.  If not, see <http://www.gnu.org/licenses/>.
 
 # pragma: no cover
+# ruff: noqa
 
 from CosmologyCore import *
 from sympy2c import OdeCombined