pybird.greenfunction module

class pybird.greenfunction.GreenFunction(Omega0_m, w=None, quintessence=False, Omega0_k=0.0, vectorize=False)[source]

Bases: object

A class to compute cosmological Green’s functions and growth factors.

The GreenFunction class handles calculations related to linear perturbation theory, including growth factors, growth rates, and Green’s functions for different cosmological models (ΛCDM, wCDM, quintessence). It provides the time-dependent functions needed for perturbation theory calculations.

Omega0_m

Present-day matter density parameter.

Type:

float

w

Dark energy equation of state parameter.

Type:

float

quintessence

Whether to include clustering quintessence.

Type:

bool

Omega0_k

Present-day curvature density parameter.

Type:

float

vectorize

Whether to vectorize calculations for arrays of parameters.

Type:

bool

OmegaL_by_Omega_m

Ratio of dark energy to matter density.

Type:

float

wcdm

Whether using wCDM cosmology.

Type:

bool

epsrel

Relative tolerance for numerical integrations.

Type:

float

C()[source]

Compute the time-dependent function C(a) for quintessence models.

H()[source]

Compute the conformal Hubble parameter H(a).

H3()[source]

Compute H(a)^-3.

Omega_m()[source]

Compute the time-dependent matter density parameter Ω_m(a).

D()[source]

Compute the linear growth factor D(a).

DD()[source]

Compute the derivative of the growth factor D’(a).

fplus()[source]

Compute the growth rate f+(a) = aD’(a)/D(a).

Dminus()[source]

Compute the decay factor D-(a).

DDminus()[source]

Compute the derivative of the decay factor D-‘(a).

fminus()[source]

Compute the decay rate f-(a) = aD-‘(a)/D-(a).

W()[source]

Compute the Wronskian W(a) = D’(a)D-(a) - D(a)D-‘(a).

G1d(), G2d(), G1t(), G2t()[source]

Compute Green’s functions for perturbation theory.

I1d(), I2d(), I1t(), I2t()[source]

Compute second order coefficients.

mG1d(), mG2d(), mG1t(), mG2t()[source]

Compute second order time integrals.

G()[source]

Compute quintessence time function.

IU1d(), IU2d(), IU1t(), IU2t()[source]

Compute third order coefficients.

IV11d(), IV12d(), IV21d(), IV22d()[source]

Compute third order coefficients.

IV11t(), IV12t(), IV21t(), IV22t()[source]

Compute third order coefficients.

mU1d(), mU2d(), mU1t(), mU2t()[source]

Compute third order time integrals.

mV11d(), mV12d(), mV21d(), mV22d()[source]

Compute third order time integrals.

mV11t(), mV12t(), mV21t(), mV22t()[source]

Compute third order time integrals.

Y()[source]

Compute third order time function.

C(a)[source]
H(a)[source]

Conformal Hubble

H3(a)[source]
Omega_m(a)[source]
D(a)[source]

Growth factor

DD(a)[source]

Derivative of growth factor

fplus(a)[source]

Growth rate

Dminus(a)[source]

Decay factor

DDminus(a)[source]

Derivative of decay factor

fminus(a)[source]

Decay rate

W(a)[source]

Wronskian

G1d(a, ai)[source]
G2d(a, ai)[source]
G1t(a, ai)[source]
G2t(a, ai)[source]
I1d(ai, a)[source]
I2d(ai, a)[source]
I1t(ai, a)[source]
I2t(ai, a)[source]
mG1d(a)[source]
mG2d(a)[source]
mG1t(a)[source]
mG2t(a)[source]
G(a)[source]
IU1d(ai, a)[source]
IU2d(ai, a)[source]
IU1t(ai, a)[source]
IU2t(ai, a)[source]
IV11d(ai, a)[source]
IV12d(ai, a)[source]
IV21d(ai, a)[source]
IV22d(ai, a)[source]
IV11t(ai, a)[source]
IV12t(ai, a)[source]
IV21t(ai, a)[source]
IV22t(ai, a)[source]
mU1d(a)[source]
mU2d(a)[source]
mU1t(a)[source]
mU2t(a)[source]
mV11d(a)[source]
mV12d(a)[source]
mV21d(a)[source]
mV22d(a)[source]
mV11t(a)[source]
mV12t(a)[source]
mV21t(a)[source]
mV22t(a)[source]
Y(a)[source]