Bimetric 3+1 toolkit for spherical symmetry
 All Data Structures Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Defines
Defines | Functions
Functions of the prime state variables

Defines

#define gBet(m, n)   eq_gBet(m,n)
#define fBet(m, n)   eq_fBet(m,n)
#define gBet_r(m, n)   eq_gBet_r(m,n)
#define fBet_r(m, n)   eq_fBet_r(m,n)
#define eq_pr(m, n)   (p(m,n)/r(m,n))
#define eq_qr(m, n)   (q(m,n)/r(m,n))

Functions

Real BimetricEvolve::gA_r (Real m, Real n)
Real BimetricEvolve::fA_r (Real m, Real n)
Real BimetricEvolve::gB_r (Real m, Real n)
Real BimetricEvolve::fB_r (Real m, Real n)
Real BimetricEvolve::gK_r (Real m, Real n)
Real BimetricEvolve::fK_r (Real m, Real n)
Real BimetricEvolve::gKD_r (Real m, Real n)
Real BimetricEvolve::fKD_r (Real m, Real n)
Real BimetricEvolve::gDA_r (Real m, Real n)
Real BimetricEvolve::fDA_r (Real m, Real n)
Real BimetricEvolve::gDB_r (Real m, Real n)
Real BimetricEvolve::fDB_r (Real m, Real n)
Real BimetricEvolve::p_r (Real m, Real n)
Real BimetricEvolve::q_r (Real m, Real n)
Real BimetricEvolve::gAlp_r (Real m, Real n)
Real BimetricEvolve::fAlp_r (Real m, Real n)
Real BimetricEvolve::gDAlp_r (Real m, Real n)
Real BimetricEvolve::fDAlp_r (Real m, Real n)
Real BimetricEvolve::pfD_r (Real m, Real n)
Real BimetricEvolve::pfS_r (Real m, Real n)
Real BimetricEvolve::pftau_r (Real m, Real n)
Real BimetricEvolve::p_rr (Real m, Real n)
Real BimetricEvolve::q_rr (Real m, Real n)
Real BimetricEvolve::fAlp_rr (Real m, Real n)
Real BimetricEvolve::gK1 (Int m, Int n)
Real BimetricEvolve::gK2 (Int m, Int n)
Real BimetricEvolve::fK1 (Int m, Int n)
Real BimetricEvolve::fK2 (Int m, Int n)
Real BimetricEvolve::eq_pr_r (Real m, Real n)
Real BimetricEvolve::eq_qr_r (Real m, Real n)
Real BimetricEvolve::eq_BM_gAlp_t (Int m, Int n)
Real BimetricEvolve::eq_BM_gDAlp_t (Int m, Int n)
Real BimetricEvolve::eq_gA_t (Int m, Int n)
Real BimetricEvolve::eq_fA_t (Int m, Int n)
Real BimetricEvolve::eq_gB_t (Int m, Int n)
Real BimetricEvolve::eq_fB_t (Int m, Int n)
Real BimetricEvolve::eq_gDA_t (Int m, Int n)
Real BimetricEvolve::eq_fDA_t (Int m, Int n)
Real BimetricEvolve::eq_gDB_t (Int m, Int n)
Real BimetricEvolve::eq_fDB_t (Int m, Int n)
Real BimetricEvolve::eq_gSig_t (Int m, Int n)
Real BimetricEvolve::eq_fSig_t (Int m, Int n)
Real BimetricEvolve::eq_base_gK_t (Int m, Int n)
Real BimetricEvolve::eq_invr_gK_t (Int m, Int n)
Real BimetricEvolve::eq_gK_t (Int m, Int n)
Real BimetricEvolve::eq_base_gKD_t (Int m, Int n)
Real BimetricEvolve::eq_invr_gKD_t (Int m, Int n)
Real BimetricEvolve::eq_gKD_t (Int m, Int n)
Real BimetricEvolve::eq_base_fK_t (Int m, Int n)
Real BimetricEvolve::eq_invr_fK_t (Int m, Int n)
Real BimetricEvolve::eq_fK_t (Int m, Int n)
Real BimetricEvolve::eq_base_fKD_t (Int m, Int n)
Real BimetricEvolve::eq_invr_fKD_t (Int m, Int n)
Real BimetricEvolve::eq_fKD_t (Int m, Int n)
Real BimetricEvolve::eq_base_gK1_t (Int m, Int n)
Real BimetricEvolve::eq_invr_gK1_t (Int m, Int n)
Real BimetricEvolve::eq_base_gK2_t (Int m, Int n)
Real BimetricEvolve::eq_invr_gK2_t (Int m, Int n)
Real BimetricEvolve::eq_base_fK1_t (Int m, Int n)
Real BimetricEvolve::eq_invr_fK1_t (Int m, Int n)
Real BimetricEvolve::eq_base_fK2_t (Int m, Int n)
Real BimetricEvolve::eq_invr_fK2_t (Int m, Int n)
Real BimetricEvolve::eq_p_t (Int m, Int n)
Real BimetricEvolve::eq_p1_t (Int m, Int n)
Real BimetricEvolve::eq_base_p_g (Int m, Int n)
Real BimetricEvolve::eq_invr_p_g (Int m, Int n)
Real BimetricEvolve::eq_p_g (Int m, Int n)
Real BimetricEvolve::eq_base_p_f (Int m, Int n)
Real BimetricEvolve::eq_invr_p_f (Int m, Int n)
Real BimetricEvolve::eq_p_f (Int m, Int n)
Real BimetricEvolve::eq_gC_rho (Int m, Int n)
Real BimetricEvolve::eq_fC_rho (Int m, Int n)
Real BimetricEvolve::eq_gC_j (Int m, Int n)
Real BimetricEvolve::eq_fC_j (Int m, Int n)
Real BimetricEvolve::eq_C_bimConsLaw (Int m, Int n)
Real BimetricEvolve::eq_j_over_p (Int m, Int n)
Real BimetricEvolve::eq_C_1 (Int m, Int n)
Real BimetricEvolve::eq_C_2 (Int m, Int n)
Real BimetricEvolve::eq_C_3 (Int m, Int n)
Real BimetricEvolve::eq_C_4 (Int m, Int n)
Real BimetricEvolve::eq_g_rho (Int m, Int n)
Real BimetricEvolve::eq_f_rho (Int m, Int n)
Real BimetricEvolve::eq_g_j (Int m, Int n)
Real BimetricEvolve::eq_f_j (Int m, Int n)
Real BimetricEvolve::eq_g_JK1 (Int m, Int n)
Real BimetricEvolve::eq_f_JK1 (Int m, Int n)
Real BimetricEvolve::eq_g_JK2 (Int m, Int n)
Real BimetricEvolve::eq_f_JK2 (Int m, Int n)
Real BimetricEvolve::eq_pf_D (Int m, Int n)
Real BimetricEvolve::eq_pf_S (Int m, Int n)
Real BimetricEvolve::eq_pf_tau (Int m, Int n)
Real BimetricEvolve::eq_pf_v (Int m, Int n)
Real BimetricEvolve::eq_pf_v_r (Int m, Int n)
Real BimetricEvolve::eq_pf_D_t (Int m, Int n)
Real BimetricEvolve::eq_pf_S_t (Int m, Int n)
Real BimetricEvolve::eq_pf_tau_t (Int m, Int n)
Real BimetricEvolve::eq_pf_rho (Int m, Int n)
Real BimetricEvolve::eq_pf_j (Int m, Int n)
Real BimetricEvolve::eq_pf_J1 (Int m, Int n)
Real BimetricEvolve::eq_pf_J2 (Int m, Int n)
Real BimetricEvolve::eq_q (Int m, Int n)
Real BimetricEvolve::eq_gAlp (Int m, Int n)
Real BimetricEvolve::eq_gW (Int m, Int n)
Real BimetricEvolve::eq_fW (Int m, Int n)
Real BimetricEvolve::eq_cW (Int m, Int n)
Real BimetricEvolve::eq_fAlp (Int m, Int n)
Real BimetricEvolve::eq_gHorz (Int m, Int n)
Real BimetricEvolve::eq_fHorz (Int m, Int n)
Real BimetricEvolve::eq_gBet (Int m, Int n)
Real BimetricEvolve::eq_gBet_r (Int m, Int n)
Real BimetricEvolve::eq_gBetr_r (Int m, Int n)
Real BimetricEvolve::eq_fBet (Int m, Int n)
Real BimetricEvolve::eq_fBet_r (Int m, Int n)
Real BimetricEvolve::eq_fBetr_r (Int m, Int n)

Define Documentation

#define eq_pr (   m,
  n 
)    (p(m,n)/r(m,n))

Definition at line 339 of file bim-solver.cpp.

#define eq_qr (   m,
  n 
)    (q(m,n)/r(m,n))

Definition at line 340 of file bim-solver.cpp.

#define fBet (   m,
  n 
)    eq_fBet(m,n)

Definition at line 333 of file bim-solver.cpp.

#define fBet_r (   m,
  n 
)    eq_fBet_r(m,n)

Definition at line 335 of file bim-solver.cpp.

#define gBet (   m,
  n 
)    eq_gBet(m,n)
Todo:
gBet, fBet to be grid functions (after calc of gAlp, fAlp)

Definition at line 332 of file bim-solver.cpp.

#define gBet_r (   m,
  n 
)    eq_gBet_r(m,n)

Definition at line 334 of file bim-solver.cpp.

Function Documentation

Real BimetricEvolve::eq_base_fK1_t ( Int  m,
Int  n 
) [private]

Definition at line 135 of file eomEvolution.h.

{
    return -(k_f * fAlp(m,n) * f_JK1(m,n)) - (fAlp(m,n) * (2 * fKD_r(m,n) + fK_r(m,n))
    /*1*/  * p(m,n)) / (3. * fA(m,n) * Lt(m,n)) + fAlp(m,n) * ((fK(m,n) * (fK(m,n) + 2
    /*3*/  * fKD(m,n))) / 3. + (-fDAlp_r(m,n) + fDA(m,n) * (fDAlp(m,n) + 2 * fDB(m,n))
    /*2*/  - 2 * fDB_r(m,n) - pow2(fDAlp(m,n)) - 2 * pow2(fDB(m,n))) / pow2(fA(m,n)))
    /*0*/  + ((2 * fKD_r(m,n) + fK_r(m,n)) * q(m,n)) / 3.;
}
Real BimetricEvolve::eq_base_fK2_t ( Int  m,
Int  n 
) [private]

Definition at line 150 of file eomEvolution.h.

{
    return -(k_f * fAlp(m,n) * f_JK2(m,n)) + (fAlp(m,n) * (fKD_r(m,n) - fK_r(m,n)) 
    /*1*/  * p(m,n)) / (3. * fA(m,n) * Lt(m,n)) + fAlp(m,n) * ((fK(m,n) * (fK(m,n) 
    /*3*/  - fKD(m,n))) / 3. + (fDA(m,n) * fDB(m,n) - fDAlp(m,n) * fDB(m,n) - fDB_r(m,n)
    /*2*/  - 2 * pow2(fDB(m,n))) / pow2(fA(m,n))) + ((-fKD_r(m,n) + fK_r(m,n)) 
    /*1*/  * q(m,n)) / 3.;
}
Real BimetricEvolve::eq_base_fK_t ( Int  m,
Int  n 
) [private]

Definition at line 187 of file eomEvolution.h.

{
    return k_f * (-(fAlp(m,n) * f_JK1(m,n)) - 2 * fAlp(m,n) * f_JK2(m,n)) - (fAlp(m,n)
    /*1*/  * fK_r(m,n) * p(m,n)) / (fA(m,n) * Lt(m,n)) + fAlp(m,n) * ((-fDAlp_r(m,n) 
    /*2*/  - 2 * fDAlp(m,n) * fDB(m,n) + fDA(m,n) * (fDAlp(m,n) + 4 * fDB(m,n)) - 4 
    /*2*/  * fDB_r(m,n) - pow2(fDAlp(m,n)) - 6 * pow2(fDB(m,n))) / pow2(fA(m,n)) 
    /*1*/  + pow2(fK(m,n))) + fK_r(m,n) * q(m,n);
}
Real BimetricEvolve::eq_base_fKD_t ( Int  m,
Int  n 
) [private]

Definition at line 202 of file eomEvolution.h.

{
    return k_f * (-(fAlp(m,n) * f_JK1(m,n)) + fAlp(m,n) * f_JK2(m,n)) - (fAlp(m,n) 
    /*1*/  * fKD_r(m,n) * p(m,n)) / (fA(m,n) * Lt(m,n)) + fAlp(m,n) * (fK(m,n) 
    /*1*/  * fKD(m,n) + (-fDAlp_r(m,n) + fDAlp(m,n) * fDB(m,n) + fDA(m,n) * (fDAlp(m,n)
    /*3*/  + fDB(m,n)) - fDB_r(m,n) - pow2(fDAlp(m,n))) / pow2(fA(m,n))) + fKD_r(m,n)
    /*0*/  * q(m,n);
}
Real BimetricEvolve::eq_base_gK1_t ( Int  m,
Int  n 
) [private]

Definition at line 27 of file eomEvolution.h.

{
    return -(k_g * gAlp(m,n) * g_JK1(m,n)) + (gAlp(m,n) * (2 * gKD_r(m,n) + gK_r(m,n))
    /*1*/  * p(m,n)) / (3. * gA(m,n) * Lt(m,n)) + gAlp(m,n) * ((gK(m,n) * (gK(m,n) + 2
    /*3*/  * gKD(m,n))) / 3. + (-gDAlp_r(m,n) + gDA(m,n) * (gDAlp(m,n) + 2 * gDB(m,n))
    /*2*/  - 2 * gDB_r(m,n) - pow2(gDAlp(m,n)) - 2 * pow2(gDB(m,n))) / pow2(gA(m,n)))
    /*0*/  + ((2 * gKD_r(m,n) + gK_r(m,n)) * q(m,n)) / 3.;
}
Real BimetricEvolve::eq_base_gK2_t ( Int  m,
Int  n 
) [private]

Definition at line 42 of file eomEvolution.h.

{
    return -(k_g * gAlp(m,n) * g_JK2(m,n)) + (gAlp(m,n) * (-gKD_r(m,n) + gK_r(m,n)) 
    /*1*/  * p(m,n)) / (3. * gA(m,n) * Lt(m,n)) + gAlp(m,n) * ((gK(m,n) * (gK(m,n) 
    /*3*/  - gKD(m,n))) / 3. + (gDA(m,n) * gDB(m,n) - gDAlp(m,n) * gDB(m,n) - gDB_r(m,n)
    /*2*/  - 2 * pow2(gDB(m,n))) / pow2(gA(m,n))) + ((-gKD_r(m,n) + gK_r(m,n)) 
    /*1*/  * q(m,n)) / 3.;
}
Real BimetricEvolve::eq_base_gK_t ( Int  m,
Int  n 
) [private]

Definition at line 79 of file eomEvolution.h.

{
    return k_g * (-(gAlp(m,n) * g_JK1(m,n)) - 2 * gAlp(m,n) * g_JK2(m,n)) + (gAlp(m,n)
    /*1*/  * gK_r(m,n) * p(m,n)) / (gA(m,n) * Lt(m,n)) + gAlp(m,n) * ((-gDAlp_r(m,n) 
    /*2*/  - 2 * gDAlp(m,n) * gDB(m,n) + gDA(m,n) * (gDAlp(m,n) + 4 * gDB(m,n)) - 4 
    /*2*/  * gDB_r(m,n) - pow2(gDAlp(m,n)) - 6 * pow2(gDB(m,n))) / pow2(gA(m,n)) 
    /*1*/  + pow2(gK(m,n))) + gK_r(m,n) * q(m,n);
}
Real BimetricEvolve::eq_base_gKD_t ( Int  m,
Int  n 
) [private]

Definition at line 94 of file eomEvolution.h.

{
    return k_g * (-(gAlp(m,n) * g_JK1(m,n)) + gAlp(m,n) * g_JK2(m,n)) + (gAlp(m,n) 
    /*1*/  * gKD_r(m,n) * p(m,n)) / (gA(m,n) * Lt(m,n)) + gAlp(m,n) * (gK(m,n) 
    /*1*/  * gKD(m,n) + (-gDAlp_r(m,n) + gDAlp(m,n) * gDB(m,n) + gDA(m,n) * (gDAlp(m,n)
    /*3*/  + gDB(m,n)) - gDB_r(m,n) - pow2(gDAlp(m,n))) / pow2(gA(m,n))) + gKD_r(m,n)
    /*0*/  * q(m,n);
}
Real BimetricEvolve::eq_base_p_f ( Int  m,
Int  n 
) [private]

Definition at line 23 of file eomMiscEquations.h.

{
    return (2 * (3 * fDB(m,n) * fKD(m,n) + fKD_r(m,n) - fK_r(m,n)) * pow2(fB(m,n))) 
    /*0*/  / (3. * k_f * gA(m,n) * pow2(gB(m,n)) * P_2_1(R(m,n)));
}
Real BimetricEvolve::eq_base_p_g ( Int  m,
Int  n 
) [private]

Definition at line 10 of file eomMiscEquations.h.

{
    return (-2 * (3 * gDB(m,n) * gKD(m,n) + gKD_r(m,n) - gK_r(m,n))) / (3. * k_g 
    /*1*/  * fA(m,n) * P_2_1(R(m,n))) + pfS(m,n) / (fA(m,n) * gA(m,n) * pow2(gB(m,n)) 
    /*1*/  * P_2_1(R(m,n)));
}
Real BimetricEvolve::eq_BM_gAlp_t ( Int  m,
Int  n 
) [inline, private]

Definition at line 349 of file bim-solver.cpp.

                                             {
        return  -2 * gAlp(m,n) * gK(m,n);
    }
Real BimetricEvolve::eq_BM_gDAlp_t ( Int  m,
Int  n 
) [inline, private]

Definition at line 353 of file bim-solver.cpp.

                                              {
        return  -2 * gK_r(m,n);
    }
Real BimetricEvolve::eq_C_1 ( Int  m,
Int  n 
) [private]

Definition at line 47 of file eomConstraints.h.

{
    return -(k_g * g_rho(m,n)) + (2 * gDA(m,n) * gDB(m,n) - 2 * gDB_r(m,n) - 3 
    /*1*/  * pow2(gDB(m,n))) / pow2(gA(m,n)) + (pow2(gK(m,n)) - pow2(gKD(m,n))) / 3. 
    /*0*/  + (gSig(m,n) / (gA(m,n) * gB(m,n)) + (2 * gDA(m,n) - 6 * gDB(m,n) 
    /*2*/  + gSig(m,n)) / pow2(gA(m,n))) / r(m,n);
}
Real BimetricEvolve::eq_C_2 ( Int  m,
Int  n 
) [private]

Definition at line 55 of file eomConstraints.h.

{
    return -(k_f * f_rho(m,n)) + (2 * fDA(m,n) * fDB(m,n) - 2 * fDB_r(m,n) - 3 
    /*1*/  * pow2(fDB(m,n))) / pow2(fA(m,n)) + (pow2(fK(m,n)) - pow2(fKD(m,n))) / 3. 
    /*0*/  + (fSig(m,n) / (fA(m,n) * fB(m,n)) + (2 * fDA(m,n) - 6 * fDB(m,n) 
    /*2*/  + fSig(m,n)) / pow2(fA(m,n))) / r(m,n);
}
Real BimetricEvolve::eq_C_3 ( Int  m,
Int  n 
) [private]

Definition at line 63 of file eomConstraints.h.

{
    return k_f * ((gA(m,n) * (3 * gDB(m,n) * gKD(m,n) + gKD_r(m,n) - gK_r(m,n))) / 3.
    /*1*/  - (k_g * pfS(m,n)) / (2. * pow2(gB(m,n)))) + (k_g * fA(m,n) * (3 * fDB(m,n)
    /*2*/  * fKD(m,n) + fKD_r(m,n) - fK_r(m,n)) * pow2(R(m,n))) / 3. + (k_f * gA(m,n)
    /*1*/  * gKD(m,n) + k_g * fA(m,n) * fKD(m,n) * pow2(R(m,n))) / r(m,n);
}
Real BimetricEvolve::eq_C_4 ( Int  m,
Int  n 
) [private]

Definition at line 71 of file eomConstraints.h.

{
    return (fA(m,n) * (2 * (fKD(m,n) * P_1_1(R(m,n)) + (-gK(m,n) + gKD(m,n)) * Lt(m,n)
    /*3*/  * P_1_2(R(m,n))) + fK(m,n) * (-2 * P_1_1(R(m,n)) + 3 * P_2_1(R(m,n))))) 
    /*0*/  / 3. + (P_2_1(R(m,n)) * p_r(m,n)) / Lt(m,n) + (p(m,n) * ((2 * fA(m,n) 
    /*3*/  * P_1_2(R(m,n))) / gA(m,n) + (2 * gA(m,n) * P_1_1(R(m,n)) * R(m,n)) 
    /*2*/  / fA(m,n))) / r(m,n) + p(m,n) * ((2 * fA(m,n) * gDB(m,n) * P_1_2(R(m,n))) 
    /*1*/  / gA(m,n) + (2 * fDB(m,n) * gA(m,n) * P_1_1(R(m,n)) * R(m,n)) / fA(m,n)) 
    /*0*/  + gA(m,n) * (-((gK(m,n) + 2 * gKD(m,n)) * P_2_1(R(m,n))) / 3. - (2 
    /*2*/  * P_1_1(R(m,n)) * (gK(m,n) - gKD(m,n) + (-fK(m,n) + fKD(m,n)) * Lt(m,n) 
    /*3*/  * R(m,n))) / 3.);
}
Real BimetricEvolve::eq_C_bimConsLaw ( Int  m,
Int  n 
) [private]

Definition at line 34 of file eomConstraints.h.

{
    return (fA(m,n) * (2 * (fKD(m,n) * P_1_1(R(m,n)) + (-gK(m,n) + gKD(m,n)) * Lt(m,n)
    /*3*/  * P_1_2(R(m,n))) + fK(m,n) * (-2 * P_1_1(R(m,n)) + 3 * P_2_1(R(m,n))))) 
    /*0*/  / 3. + (P_2_1(R(m,n)) * p_r(m,n)) / Lt(m,n) + (p(m,n) * ((2 * fA(m,n) 
    /*3*/  * P_1_2(R(m,n))) / gA(m,n) + (2 * gA(m,n) * P_1_1(R(m,n)) * R(m,n)) 
    /*2*/  / fA(m,n))) / r(m,n) + p(m,n) * ((2 * fA(m,n) * gDB(m,n) * P_1_2(R(m,n))) 
    /*1*/  / gA(m,n) + (2 * fDB(m,n) * gA(m,n) * P_1_1(R(m,n)) * R(m,n)) / fA(m,n)) 
    /*0*/  + gA(m,n) * (-((gK(m,n) + 2 * gKD(m,n)) * P_2_1(R(m,n))) / 3. - (2 
    /*2*/  * P_1_1(R(m,n)) * (gK(m,n) - gKD(m,n) + (-fK(m,n) + fKD(m,n)) * Lt(m,n) 
    /*3*/  * R(m,n))) / 3.);
}
Real BimetricEvolve::eq_cW ( Int  m,
Int  n 
) [private]

Definition at line 144 of file eomLapseRatios.h.

{
    return 2 * pow2(fA(m,n)) * pow2(gA(m,n)) * pow2(gB(m,n)) * pow2(r(m,n)) 
    /*0*/  * pow2(R(m,n)) * pow3(Lt(m,n)) * P_2_1(R(m,n));
}
Real BimetricEvolve::eq_f_j ( Int  m,
Int  n 
) [private]

Definition at line 39 of file eomSources.h.

{
    return (gA(m,n) * p(m,n) * pow2(gB(m,n)) * P_2_1(R(m,n))) / pow2(fB(m,n));
}
Real BimetricEvolve::eq_f_JK1 ( Int  m,
Int  n 
) [private]

Definition at line 44 of file eomSources.h.

{
    return -(P_1_3(R(m,n)) / R(m,n)) + (gA(m,n) * ((1 - 2 * pow2(Lt(m,n))) 
    /*2*/  * P_2_1(R(m,n)) - 2 * P_1_2(R(m,n)) * R(m,n))) / (2. * fA(m,n) * Lt(m,n) 
    /*1*/  * pow2(R(m,n))) + gAlp(m,n) * (-((gA(m,n) * P_1_1(R(m,n))) / (fA(m,n) 
    /*3*/  * fAlp(m,n) * R(m,n))) + (P_2_1(R(m,n)) - 2 * P_1_2(R(m,n)) * R(m,n)) / (2. 
    /*2*/  * fAlp(m,n) * Lt(m,n) * pow2(R(m,n))));
}
Real BimetricEvolve::eq_f_JK2 ( Int  m,
Int  n 
) [private]

Definition at line 53 of file eomSources.h.

{
    return -(gA(m,n) * P_2_1(R(m,n))) / (2. * fA(m,n) * Lt(m,n) * pow2(R(m,n))) 
    /*0*/  - (gAlp(m,n) * P_2_1(R(m,n))) / (2. * fAlp(m,n) * Lt(m,n) * pow2(R(m,n))) 
    /*0*/  - P_2_2(R(m,n)) / pow2(R(m,n));
}
Real BimetricEvolve::eq_f_rho ( Int  m,
Int  n 
) [private]

Definition at line 33 of file eomSources.h.

{
    return -((gA(m,n) * Lt(m,n) * P_2_1(R(m,n))) / (fA(m,n) * pow2(R(m,n)))) 
    /*0*/  - P_2_2(R(m,n)) / pow2(R(m,n));
}
Real BimetricEvolve::eq_fA_t ( Int  m,
Int  n 
) [private]

Definition at line 114 of file eomEvolution.h.

{
    return -((fAlp(m,n) * fDAlp(m,n) * p(m,n)) / Lt(m,n)) + (fAlp(m,n) * pow2(p(m,n))
    /*1*/  * p_r(m,n)) / pow3(Lt(m,n)) + fAlp(m,n) * ((fA(m,n) * (-fK(m,n) - 2 
    /*3*/  * fKD(m,n))) / 3. - p_r(m,n) / Lt(m,n)) + fA(m,n) * fDA(m,n) * q(m,n) 
    /*0*/  + fA(m,n) * q_r(m,n);
}
Real BimetricEvolve::eq_fAlp ( Int  m,
Int  n 
) [private]

Definition at line 16 of file eomGauge.h.

{
    return -((gAlp(m,n) * gW(m,n)) / fW(m,n));
}
Real BimetricEvolve::eq_fB_t ( Int  m,
Int  n 
) [private]

Definition at line 122 of file eomEvolution.h.

{
    return eq_qr(m,n) * fB(m,n) + (fAlp(m,n) * fB(m,n) * (-fK(m,n) + fKD(m,n))) / 3. 
    /*0*/  - (eq_pr(m,n) * fAlp(m,n) * fB(m,n)) / (fA(m,n) * Lt(m,n)) - (fAlp(m,n) 
    /*1*/  * fB(m,n) * fDB(m,n) * p(m,n)) / (fA(m,n) * Lt(m,n)) + fB(m,n) * fDB(m,n) 
    /*0*/  * q(m,n);
}
Real BimetricEvolve::eq_fBet ( Int  m,
Int  n 
) [private]

Definition at line 53 of file eomMiscEquations.h.

{
    return -((fAlp(m,n) * p(m,n)) / (fA(m,n) * Lt(m,n))) + q(m,n);
}
Real BimetricEvolve::eq_fBet_r ( Int  m,
Int  n 
) [private]

Definition at line 64 of file eomMiscEquations.h.

{
    return (fAlp(m,n) * (fDA(m,n) - fDAlp(m,n)) * p(m,n)) / (fA(m,n) * Lt(m,n)) 
    /*0*/  - (fAlp(m,n) * p_r(m,n)) / (fA(m,n) * pow3(Lt(m,n))) + q_r(m,n);
}
Real BimetricEvolve::eq_fBetr_r ( Int  m,
Int  n 
) [private]

Definition at line 78 of file eomMiscEquations.h.

{
    return eq_qr_r(m,n) - (eq_pr_r(m,n) * fAlp(m,n)) / (fA(m,n) * pow3(Lt(m,n))) 
    /*0*/  + (fAlp(m,n) * pow3(p(m,n))) / (fA(m,n) * pow2(r(m,n)) * pow3(Lt(m,n))) 
    /*0*/  + (fAlp(m,n) * (fDA(m,n) - fDAlp(m,n)) * p(m,n)) / (fA(m,n) * Lt(m,n) 
    /*1*/  * r(m,n));
}
Real BimetricEvolve::eq_fC_j ( Int  m,
Int  n 
) [private]

Definition at line 28 of file eomConstraints.h.

{
    return (2 * (3 * fDB(m,n) * fKD(m,n) + fKD_r(m,n) - fK_r(m,n))) / 3. + (2 
    /*1*/  * fKD(m,n)) / r(m,n);
}
Real BimetricEvolve::eq_fC_rho ( Int  m,
Int  n 
) [private]

Definition at line 14 of file eomConstraints.h.

{
    return (2 * fDA(m,n) * fDB(m,n) - 2 * fDB_r(m,n) - 3 * pow2(fDB(m,n))) 
    /*0*/  / pow2(fA(m,n)) + (pow2(fK(m,n)) - pow2(fKD(m,n))) / 3. + (fSig(m,n) 
    /*1*/  / (fA(m,n) * fB(m,n)) + (2 * fDA(m,n) - 6 * fDB(m,n) + fSig(m,n)) 
    /*1*/  / pow2(fA(m,n))) / r(m,n);
}
Real BimetricEvolve::eq_fDA_t ( Int  m,
Int  n 
) [private]

Definition at line 159 of file eomEvolution.h.

{
    return (fAlp(m,n) * p(m,n) * (fDA(m,n) * fDAlp(m,n) * Power(Lt(m,n),4) 
    /*2*/  - fDAlp_r(m,n) * Power(Lt(m,n),4) - Power(Lt(m,n),4) * pow2(fDAlp(m,n)) + 3 
    /*2*/  * pow2(p_r(m,n)))) / (fA(m,n) * Power(Lt(m,n),5)) + fAlp(m,n) 
    /*0*/  * ((-(fDAlp(m,n) * (fK(m,n) + 2 * fKD(m,n))) - 2 * fKD_r(m,n) - fK_r(m,n)) 
    /*1*/  / 3. + (fDA(m,n) * p_r(m,n) - 2 * fDAlp(m,n) * p_r(m,n) - p_rr(m,n)) 
    /*1*/  / (fA(m,n) * pow3(Lt(m,n)))) + fDA_r(m,n) * q(m,n) + fDA(m,n) * q_r(m,n) 
    /*0*/  + q_rr(m,n);
}
Real BimetricEvolve::eq_fDB_t ( Int  m,
Int  n 
) [private]

Definition at line 170 of file eomEvolution.h.

{
    return eq_qr_r(m,n) + (eq_pr(m,n) * fAlp(m,n) * (fDA(m,n) - fDAlp(m,n))) 
    /*0*/  / (fA(m,n) * Lt(m,n)) - (fAlp(m,n) * (-(fDA(m,n) * fDB(m,n)) + fDAlp(m,n) 
    /*2*/  * fDB(m,n) + fDB_r(m,n)) * p(m,n)) / (fA(m,n) * Lt(m,n)) + (fAlp(m,n) 
    /*1*/  * p(m,n) * pow2(eq_pr(m,n))) / (fA(m,n) * pow3(Lt(m,n))) + fAlp(m,n) 
    /*0*/  * ((fDAlp(m,n) * (-fK(m,n) + fKD(m,n)) + fKD_r(m,n) - fK_r(m,n)) / 3. 
    /*1*/  - (eq_pr_r(m,n) + fDB(m,n) * p_r(m,n)) / (fA(m,n) * pow3(Lt(m,n)))) 
    /*0*/  + fDB_r(m,n) * q(m,n) + fDB(m,n) * q_r(m,n);
}
Real BimetricEvolve::eq_fHorz ( Int  m,
Int  n 
) [private]

Definition at line 91 of file eomMiscEquations.h.

{
    return fDB(m,n) + (fA(m,n) * (-fK(m,n) + fKD(m,n))) / 3. + 1 / r(m,n);
}
Real BimetricEvolve::eq_fK_t ( Int  m,
Int  n 
) [inline, private]

Definition at line 388 of file bim-solver.cpp.

                                 {
        return eq_base_fK_t(m,n) + eq_invr_fK_t(m,n) / r(m,n);
    }
Real BimetricEvolve::eq_fKD_t ( Int  m,
Int  n 
) [inline, private]

Definition at line 395 of file bim-solver.cpp.

                                  {
        return eq_base_fKD_t(m,n) + eq_invr_fKD_t(m,n) / r(m,n);
    }
Real BimetricEvolve::eq_fSig_t ( Int  m,
Int  n 
) [private]

Definition at line 211 of file eomEvolution.h.

{
    return (eq_qr_r(m,n) * fA(m,n)) / fB(m,n) + (eq_qr(m,n) * fA(m,n) * (fDA(m,n) 
    /*2*/  - fDB(m,n))) / fB(m,n) + (eq_pr(m,n) * fAlp(m,n) * (-fDAlp(m,n) + fDB(m,n)))
    /*0*/  / (fB(m,n) * Lt(m,n)) + fAlp(m,n) * ((fA(m,n) * (3 * fDB(m,n) * fKD(m,n) 
    /*3*/  + fKD_r(m,n) - fK_r(m,n))) / (3. * fB(m,n)) - eq_pr_r(m,n) / (fB(m,n) 
    /*2*/  * pow3(Lt(m,n)))) + (fAlp(m,n) * p(m,n) * pow2(eq_pr(m,n))) / (fB(m,n) 
    /*1*/  * pow3(Lt(m,n))) - (k_f * fA(m,n) * fAlp(m,n) * gA(m,n) * p(m,n) 
    /*1*/  * P_2_1(R(m,n))) / (2. * fB(m,n) * pow2(R(m,n)));
}
Real BimetricEvolve::eq_fW ( Int  m,
Int  n 
) [private]

Definition at line 75 of file eomLapseRatios.h.

{
    return (4 * fA(m,n) * gA(m,n) * Lt(m,n) * p(m,n) * pow2(gB(m,n)) * pow2(R(m,n)) 
    /*1*/  * r(m,n) * (-(gA(m,n) * (fA(m,n) * Lt(m,n) * ((fK(m,n) - fKD(m,n)) * gDB(m,n)
    /*5*/  * (4 * P_1_1(R(m,n)) * P_1_2(R(m,n)) - P_0_2(R(m,n)) * P_2_1(R(m,n))) 
    /*5*/  + fDB(m,n) * (2 * fKD(m,n) * P_1_1(R(m,n)) * P_1_2(R(m,n)) + fK(m,n) 
    /*6*/  * P_1_2(R(m,n)) * (-2 * P_1_1(R(m,n)) + 3 * P_2_1(R(m,n))) + (gK(m,n) 
    /*7*/  - gKD(m,n)) * Lt(m,n) * (2 * pow2(P_1_2(R(m,n))) - P_0_3(R(m,n)) 
    /*7*/  * P_2_1(R(m,n))))) + 3 * fDB(m,n) * P_1_2(R(m,n)) * P_2_1(R(m,n)) * p_r(m,n))
    /*3*/  * R(m,n)) + fDB(m,n) * Lt(m,n) * pow2(gA(m,n)) * R(m,n) * (((gK(m,n) 
    /*5*/  - gKD(m,n)) * P_0_2(R(m,n)) + 2 * (fK(m,n) + 2 * fKD(m,n)) * Lt(m,n) 
    /*4*/  * P_1_1(R(m,n)) + (gK(m,n) + 2 * gKD(m,n)) * P_1_2(R(m,n))) * P_2_1(R(m,n)) 
    /*3*/  - 2 * (fK(m,n) - fKD(m,n)) * Lt(m,n) * (pow2(P_1_1(R(m,n))) - P_1_1(R(m,n)) 
    /*4*/  * P_1_2(R(m,n)) * R(m,n) + P_0_2(R(m,n)) * P_2_1(R(m,n)) * R(m,n))) 
    /*2*/  - gDB(m,n) * Lt(m,n) * pow2(fA(m,n)) * (-2 * (gK(m,n) - gKD(m,n)) 
    /*3*/  * pow2(P_1_2(R(m,n))) + (fK(m,n) - fKD(m,n)) * Lt(m,n) * (2 * P_1_2(R(m,n)) 
    /*4*/  * (P_1_1(R(m,n)) - P_2_1(R(m,n))) + P_0_3(R(m,n)) * P_2_1(R(m,n)) 
    /*4*/  * R(m,n))))) / 3. + 2 * pow2(Lt(m,n)) * (Lt(m,n) * pow3(fA(m,n)) * (-4 
    /*2*/  * pow2(gB(m,n)) * pow2(gDB(m,n)) * pow2(P_1_2(R(m,n))) * pow2(R(m,n)) 
    /*2*/  + pow2(gA(m,n)) * P_2_1(R(m,n)) * (-P_1_1(R(m,n)) + P_1_2(R(m,n)) * R(m,n)))
    /*1*/  + pow2(fDB(m,n)) * pow2(gB(m,n)) * pow3(gA(m,n)) * pow3(R(m,n)) * (-2 
    /*2*/  * pow2(p(m,n)) * P_0_2(R(m,n)) * P_2_1(R(m,n)) * R(m,n) - P_1_1(R(m,n)) 
    /*2*/  * (P_2_1(R(m,n)) - 4 * pow2(Lt(m,n)) * P_1_2(R(m,n)) * R(m,n))) + fA(m,n) 
    /*1*/  * fDB(m,n) * Lt(m,n) * pow2(gA(m,n)) * pow2(gB(m,n)) * pow2(R(m,n)) * (-2 
    /*2*/  * gDB(m,n) * R(m,n) * (2 * P_1_1(R(m,n)) * P_1_2(R(m,n)) + P_1_2(R(m,n)) 
    /*3*/  * P_2_1(R(m,n)) - 2 * pow2(P_1_2(R(m,n))) * R(m,n) + P_0_3(R(m,n)) 
    /*3*/  * P_2_1(R(m,n)) * R(m,n)) + fDB(m,n) * (P_1_1(R(m,n)) * P_2_1(R(m,n)) 
    /*3*/  - R(m,n) * (P_1_2(R(m,n)) * P_2_1(R(m,n)) + 4 * pow2(P_1_2(R(m,n))) * R(m,n)
    /*4*/  - 2 * P_0_3(R(m,n)) * P_2_1(R(m,n)) * R(m,n)))) + gA(m,n) * pow2(fA(m,n)) 
    /*1*/  * R(m,n) * (pow2(gA(m,n)) * P_2_1(R(m,n)) * (P_1_1(R(m,n)) - P_1_2(R(m,n)) 
    /*3*/  * R(m,n)) + gDB(m,n) * pow2(gB(m,n)) * R(m,n) * (2 * fDB(m,n) * (2 * (2 
    /*5*/  + pow2(p(m,n))) * pow2(P_1_2(R(m,n))) - pow2(Lt(m,n)) * P_0_3(R(m,n)) 
    /*4*/  * P_2_1(R(m,n))) * R(m,n) + gDB(m,n) * (3 * P_1_2(R(m,n)) * P_2_1(R(m,n)) - 4
    /*4*/  * pow2(P_1_2(R(m,n))) * R(m,n) + 2 * P_0_3(R(m,n)) * P_2_1(R(m,n)) 
    /*4*/  * R(m,n))))) + pow2(r(m,n)) * (k_g * pow2(fA(m,n)) * pow2(gA(m,n)) 
    /*1*/  * pow2(Lt(m,n)) * pow2(R(m,n)) * P_2_1(R(m,n)) * (2 * P_1_2(R(m,n)) 
    /*2*/  * (pfD(m,n) + pftau(m,n) - 3 * fA(m,n) * Lt(m,n) * pow2(gB(m,n)) 
    /*3*/  * P_2_1(R(m,n))) + gA(m,n) * pow2(gB(m,n)) * (-2 * P_1_2(R(m,n)) 
    /*3*/  * P_2_0(R(m,n)) + P_2_1(R(m,n)) * (-4 * P_1_1(R(m,n)) + P_2_1(R(m,n))))) 
    /*1*/  + k_f * fA(m,n) * pow2(gA(m,n)) * pow2(gB(m,n)) * pow2(Lt(m,n)) 
    /*1*/  * P_2_1(R(m,n)) * (fA(m,n) * gA(m,n) * (-pow2(P_2_1(R(m,n))) + 4 
    /*3*/  * P_1_2(R(m,n)) * P_2_1(R(m,n)) * R(m,n) + 2 * P_1_1(R(m,n)) * P_2_2(R(m,n))
    /*3*/  * R(m,n)) + 2 * Lt(m,n) * pow2(fA(m,n)) * (-(P_1_1(R(m,n)) * P_2_2(R(m,n)))
    /*3*/  + (P_1_3(R(m,n)) * P_2_1(R(m,n)) + P_1_2(R(m,n)) * P_2_2(R(m,n))) * R(m,n))
    /*2*/  + 2 * gA(m,n) * pow2(p(m,n)) * P_2_1(R(m,n)) * (-2 * gA(m,n) * Lt(m,n) 
    /*3*/  * P_1_1(R(m,n)) * R(m,n) + fA(m,n) * (P_1_1(R(m,n)) - P_2_1(R(m,n)) 
    /*4*/  + P_1_2(R(m,n)) * R(m,n)))) + (2 * fA(m,n) * Lt(m,n) * pow2(gA(m,n)) 
    /*2*/  * pow2(gB(m,n)) * pow2(R(m,n)) * (2 * (gK(m,n) - gKD(m,n)) * Lt(m,n) 
    /*3*/  * pow2(fA(m,n)) * pow2(p(m,n)) * P_1_2(R(m,n)) * (-2 * fKD(m,n) 
    /*4*/  * P_1_1(R(m,n)) + fK(m,n) * (2 * P_1_1(R(m,n)) - 3 * P_2_1(R(m,n)))) - 12 
    /*3*/  * (-fKD_r(m,n) + fK_r(m,n)) * gA(m,n) * p(m,n) * pow2(Lt(m,n)) 
    /*3*/  * P_1_1(R(m,n)) * P_2_1(R(m,n)) * R(m,n) + 4 * fA(m,n) * gA(m,n) * Lt(m,n) 
    /*3*/  * pow2(fK(m,n) - fKD(m,n)) * pow2(p(m,n)) * P_1_1(R(m,n)) * (-P_1_1(R(m,n)) 
    /*4*/  + P_2_1(R(m,n))) * R(m,n) - fA(m,n) * gA(m,n) * Lt(m,n) * P_2_1(R(m,n)) * (3
    /*4*/  * pow2(gK(m,n)) * P_1_2(R(m,n)) - 3 * pow2(gKD(m,n)) * P_1_2(R(m,n)) + 2 
    /*4*/  * (fK(m,n) - fKD(m,n)) * gKD(m,n) * Lt(m,n) * (P_0_2(R(m,n)) - 2 
    /*5*/  * P_1_2(R(m,n))) * R(m,n) - 2 * (fK(m,n) - fKD(m,n)) * gK(m,n) * Lt(m,n) 
    /*4*/  * (P_0_2(R(m,n)) + P_1_2(R(m,n))) * R(m,n) + pow2(fK(m,n) - fKD(m,n)) 
    /*4*/  * R(m,n) * (-P_1_1(R(m,n)) + 2 * pow2(Lt(m,n)) * P_0_2(R(m,n)) * R(m,n))) 
    /*3*/  + fA(m,n) * P_2_1(R(m,n)) * (-6 * P_1_2(R(m,n)) * p_r(m,n) * (-gK(m,n) 
    /*5*/  + gKD(m,n) + (fK(m,n) - fKD(m,n)) * Lt(m,n) * R(m,n)) + fA(m,n) 
    /*4*/  * pow2(Lt(m,n)) * (-((fK(m,n) - 7 * fKD(m,n)) * (fK(m,n) - fKD(m,n)) 
    /*6*/  * P_1_1(R(m,n))) + 2 * (gK(m,n) - gKD(m,n)) * Lt(m,n) * ((fK(m,n) + 2 
    /*7*/  * fKD(m,n)) * P_1_2(R(m,n)) + (fK(m,n) - fKD(m,n)) * P_0_3(R(m,n)) * R(m,n))
    /*5*/  + (fK(m,n) - fKD(m,n)) * (fK(m,n) * (2 * P_0_2(R(m,n)) - 3 * P_1_2(R(m,n)))
    /*6*/  * R(m,n) + fKD(m,n) * (-6 * P_2_1(R(m,n)) + (-2 * P_0_2(R(m,n)) + 3 
    /*8*/  * P_1_2(R(m,n))) * R(m,n))))))) / 9.);
}
Real BimetricEvolve::eq_g_j ( Int  m,
Int  n 
) [private]

Definition at line 12 of file eomSources.h.

{
    return pfS(m,n) / (gA(m,n) * pow2(gB(m,n))) - fA(m,n) * p(m,n) * P_2_1(R(m,n));
}
Real BimetricEvolve::eq_g_JK1 ( Int  m,
Int  n 
) [private]

Definition at line 17 of file eomSources.h.

{
    return (pfD(m,n) + pfS(m,n) * pfv(m,n) + pftau(m,n)) / (2. * gA(m,n) 
    /*1*/  * pow2(gB(m,n))) - P_1_0(R(m,n)) + (fA(m,n) * (-2 * P_1_1(R(m,n)) + (1 - 2 
    /*3*/  * pow2(Lt(m,n))) * P_2_1(R(m,n)))) / (2. * gA(m,n) * Lt(m,n)) + (-((fA(m,n) 
    /*3*/  * fAlp(m,n) * P_1_2(R(m,n))) / gA(m,n)) + (fAlp(m,n) * (-2 * P_1_1(R(m,n)) 
    /*3*/  + P_2_1(R(m,n)))) / (2. * Lt(m,n))) / gAlp(m,n);
}
Real BimetricEvolve::eq_g_JK2 ( Int  m,
Int  n 
) [private]

Definition at line 26 of file eomSources.h.

{
    return (pfD(m,n) - pfS(m,n) * pfv(m,n) + pftau(m,n)) / (2. * gA(m,n) 
    /*1*/  * pow2(gB(m,n))) - P_2_0(R(m,n)) - (fA(m,n) * P_2_1(R(m,n))) / (2. * gA(m,n)
    /*1*/  * Lt(m,n)) - (fAlp(m,n) * P_2_1(R(m,n))) / (2. * gAlp(m,n) * Lt(m,n));
}
Real BimetricEvolve::eq_g_rho ( Int  m,
Int  n 
) [private]

Definition at line 6 of file eomSources.h.

{
    return (pfD(m,n) + pftau(m,n)) / (gA(m,n) * pow2(gB(m,n))) - P_2_0(R(m,n)) 
    /*0*/  - (fA(m,n) * Lt(m,n) * P_2_1(R(m,n))) / gA(m,n);
}
Real BimetricEvolve::eq_gA_t ( Int  m,
Int  n 
) [private]

Definition at line 6 of file eomEvolution.h.

{
    return (gAlp(m,n) * gDAlp(m,n) * p(m,n)) / Lt(m,n) - (gAlp(m,n) * pow2(p(m,n)) 
    /*1*/  * p_r(m,n)) / pow3(Lt(m,n)) + gAlp(m,n) * ((gA(m,n) * (-gK(m,n) - 2 
    /*3*/  * gKD(m,n))) / 3. + p_r(m,n) / Lt(m,n)) + gA(m,n) * gDA(m,n) * q(m,n) 
    /*0*/  + gA(m,n) * q_r(m,n);
}
Real BimetricEvolve::eq_gAlp ( Int  m,
Int  n 
) [private]

Definition at line 11 of file eomGauge.h.

{
    return 1;
}
Real BimetricEvolve::eq_gB_t ( Int  m,
Int  n 
) [private]

Definition at line 14 of file eomEvolution.h.

{
    return eq_qr(m,n) * gB(m,n) + (gAlp(m,n) * gB(m,n) * (-gK(m,n) + gKD(m,n))) / 3. 
    /*0*/  + (eq_pr(m,n) * gAlp(m,n) * gB(m,n)) / (gA(m,n) * Lt(m,n)) + (gAlp(m,n) 
    /*1*/  * gB(m,n) * gDB(m,n) * p(m,n)) / (gA(m,n) * Lt(m,n)) + gB(m,n) * gDB(m,n) 
    /*0*/  * q(m,n);
}
Real BimetricEvolve::eq_gBet ( Int  m,
Int  n 
) [private]

Definition at line 48 of file eomMiscEquations.h.

{
    return (gAlp(m,n) * p(m,n)) / (gA(m,n) * Lt(m,n)) + q(m,n);
}
Real BimetricEvolve::eq_gBet_r ( Int  m,
Int  n 
) [private]

Definition at line 58 of file eomMiscEquations.h.

{
    return (gAlp(m,n) * (-gDA(m,n) + gDAlp(m,n)) * p(m,n)) / (gA(m,n) * Lt(m,n)) 
    /*0*/  + (gAlp(m,n) * p_r(m,n)) / (gA(m,n) * pow3(Lt(m,n))) + q_r(m,n);
}
Real BimetricEvolve::eq_gBetr_r ( Int  m,
Int  n 
) [private]

Definition at line 70 of file eomMiscEquations.h.

{
    return eq_qr_r(m,n) + (eq_pr_r(m,n) * gAlp(m,n)) / (gA(m,n) * pow3(Lt(m,n))) 
    /*0*/  - (gAlp(m,n) * pow3(p(m,n))) / (gA(m,n) * pow2(r(m,n)) * pow3(Lt(m,n))) 
    /*0*/  + (gAlp(m,n) * (-gDA(m,n) + gDAlp(m,n)) * p(m,n)) / (gA(m,n) * Lt(m,n) 
    /*1*/  * r(m,n));
}
Real BimetricEvolve::eq_gC_j ( Int  m,
Int  n 
) [private]

Definition at line 22 of file eomConstraints.h.

{
    return (2 * (3 * gDB(m,n) * gKD(m,n) + gKD_r(m,n) - gK_r(m,n))) / 3. + (2 
    /*1*/  * gKD(m,n)) / r(m,n);
}
Real BimetricEvolve::eq_gC_rho ( Int  m,
Int  n 
) [private]

Definition at line 6 of file eomConstraints.h.

{
    return (2 * gDA(m,n) * gDB(m,n) - 2 * gDB_r(m,n) - 3 * pow2(gDB(m,n))) 
    /*0*/  / pow2(gA(m,n)) + (pow2(gK(m,n)) - pow2(gKD(m,n))) / 3. + (gSig(m,n) 
    /*1*/  / (gA(m,n) * gB(m,n)) + (2 * gDA(m,n) - 6 * gDB(m,n) + gSig(m,n)) 
    /*1*/  / pow2(gA(m,n))) / r(m,n);
}
Real BimetricEvolve::eq_gDA_t ( Int  m,
Int  n 
) [private]

Definition at line 51 of file eomEvolution.h.

{
    return (gAlp(m,n) * p(m,n) * (-(gDA(m,n) * gDAlp(m,n) * Power(Lt(m,n),4)) 
    /*2*/  + gDAlp_r(m,n) * Power(Lt(m,n),4) + Power(Lt(m,n),4) * pow2(gDAlp(m,n)) - 3 
    /*2*/  * pow2(p_r(m,n)))) / (gA(m,n) * Power(Lt(m,n),5)) + gAlp(m,n) 
    /*0*/  * ((-(gDAlp(m,n) * (gK(m,n) + 2 * gKD(m,n))) - 2 * gKD_r(m,n) - gK_r(m,n)) 
    /*1*/  / 3. + (-(gDA(m,n) * p_r(m,n)) + 2 * gDAlp(m,n) * p_r(m,n) + p_rr(m,n)) 
    /*1*/  / (gA(m,n) * pow3(Lt(m,n)))) + gDA_r(m,n) * q(m,n) + gDA(m,n) * q_r(m,n) 
    /*0*/  + q_rr(m,n);
}
Real BimetricEvolve::eq_gDB_t ( Int  m,
Int  n 
) [private]

Definition at line 62 of file eomEvolution.h.

{
    return eq_qr_r(m,n) + (eq_pr(m,n) * gAlp(m,n) * (-gDA(m,n) + gDAlp(m,n))) 
    /*0*/  / (gA(m,n) * Lt(m,n)) + (gAlp(m,n) * (-(gDA(m,n) * gDB(m,n)) + gDAlp(m,n) 
    /*2*/  * gDB(m,n) + gDB_r(m,n)) * p(m,n)) / (gA(m,n) * Lt(m,n)) - (gAlp(m,n) 
    /*1*/  * p(m,n) * pow2(eq_pr(m,n))) / (gA(m,n) * pow3(Lt(m,n))) + gAlp(m,n) 
    /*0*/  * ((gDAlp(m,n) * (-gK(m,n) + gKD(m,n)) + gKD_r(m,n) - gK_r(m,n)) / 3. 
    /*1*/  + (eq_pr_r(m,n) + gDB(m,n) * p_r(m,n)) / (gA(m,n) * pow3(Lt(m,n)))) 
    /*0*/  + gDB_r(m,n) * q(m,n) + gDB(m,n) * q_r(m,n);
}
Real BimetricEvolve::eq_gHorz ( Int  m,
Int  n 
) [private]

Definition at line 86 of file eomMiscEquations.h.

{
    return gDB(m,n) + (gA(m,n) * (-gK(m,n) + gKD(m,n))) / 3. + 1 / r(m,n);
}
Real BimetricEvolve::eq_gK_t ( Int  m,
Int  n 
) [inline, private]

Definition at line 374 of file bim-solver.cpp.

                                 {
        return eq_base_gK_t(m,n) + eq_invr_gK_t(m,n) / r(m,n);
    }
Real BimetricEvolve::eq_gKD_t ( Int  m,
Int  n 
) [inline, private]

Definition at line 381 of file bim-solver.cpp.

                                  {
        return eq_base_gKD_t(m,n) + eq_invr_gKD_t(m,n) / r(m,n);
    }
Real BimetricEvolve::eq_gSig_t ( Int  m,
Int  n 
) [private]

Definition at line 103 of file eomEvolution.h.

{
    return (eq_qr_r(m,n) * gA(m,n)) / gB(m,n) + (eq_qr(m,n) * gA(m,n) * (gDA(m,n) 
    /*2*/  - gDB(m,n))) / gB(m,n) + (eq_pr(m,n) * gAlp(m,n) * (gDAlp(m,n) - gDB(m,n))) 
    /*0*/  / (gB(m,n) * Lt(m,n)) - (k_g * gAlp(m,n) * pfS(m,n)) / (2. * pow3(gB(m,n))) 
    /*0*/  + gAlp(m,n) * ((gA(m,n) * (3 * gDB(m,n) * gKD(m,n) + gKD_r(m,n) - gK_r(m,n)))
    /*1*/  / (3. * gB(m,n)) + eq_pr_r(m,n) / (gB(m,n) * pow3(Lt(m,n)))) - (gAlp(m,n) 
    /*1*/  * p(m,n) * pow2(eq_pr(m,n))) / (gB(m,n) * pow3(Lt(m,n))) + (k_g * fA(m,n) 
    /*1*/  * gA(m,n) * gAlp(m,n) * p(m,n) * P_2_1(R(m,n))) / (2. * gB(m,n));
}
Real BimetricEvolve::eq_gW ( Int  m,
Int  n 
) [private]

Definition at line 5 of file eomLapseRatios.h.

{
    return (4 * fA(m,n) * gA(m,n) * Lt(m,n) * p(m,n) * pow2(gB(m,n)) * pow2(R(m,n)) 
    /*1*/  * r(m,n) * (fDB(m,n) * Lt(m,n) * pow2(gA(m,n)) * R(m,n) * (2 * (fK(m,n) 
    /*4*/  - fKD(m,n)) * pow2(P_1_1(R(m,n))) * R(m,n) + (gK(m,n) - gKD(m,n)) * Lt(m,n) 
    /*3*/  * (P_1_1(R(m,n)) * (2 * P_1_1(R(m,n)) + P_2_1(R(m,n))) + P_0_2(R(m,n)) 
    /*4*/  * P_2_1(R(m,n)) * R(m,n))) + fA(m,n) * gDB(m,n) * (-3 * P_1_2(R(m,n)) 
    /*3*/  * P_2_1(R(m,n)) * p_r(m,n) * R(m,n) + fA(m,n) * Lt(m,n) * (2 * Lt(m,n) 
    /*4*/  * P_1_2(R(m,n)) * (2 * (gK(m,n) - gKD(m,n)) * P_1_1(R(m,n)) + (gK(m,n) + 2 
    /*6*/  * gKD(m,n)) * P_2_1(R(m,n))) + (-(fKD(m,n) * P_0_2(R(m,n))) + 2 * (gK(m,n) 
    /*6*/  - gKD(m,n)) * Lt(m,n) * P_0_3(R(m,n)) + fK(m,n) * (P_0_2(R(m,n)) - 3 
    /*6*/  * P_1_2(R(m,n)))) * P_2_1(R(m,n)) * R(m,n))) + fA(m,n) * gA(m,n) * Lt(m,n) 
    /*2*/  * R(m,n) * (fDB(m,n) * (gK(m,n) - gKD(m,n)) * (-4 * P_1_1(R(m,n)) 
    /*4*/  * P_1_2(R(m,n)) + P_0_2(R(m,n)) * P_2_1(R(m,n))) - gDB(m,n) 
    /*3*/  * (-(P_1_2(R(m,n)) * (2 * (gK(m,n) - gKD(m,n)) * P_1_1(R(m,n)) + (gK(m,n) + 2
    /*7*/  * gKD(m,n)) * P_2_1(R(m,n)))) + (fK(m,n) - fKD(m,n)) * Lt(m,n) 
    /*4*/  * (P_1_1(R(m,n)) * (2 * P_1_1(R(m,n)) + P_2_1(R(m,n))) + P_0_2(R(m,n)) 
    /*5*/  * P_2_1(R(m,n)) * R(m,n)))))) / 3. + pow2(r(m,n)) * (k_g * gA(m,n) 
    /*1*/  * pow2(fA(m,n)) * pow2(Lt(m,n)) * pow2(R(m,n)) * P_2_1(R(m,n)) * (2 * fA(m,n)
    /*2*/  * (-pftau(m,n) + pfD(m,n) * (-1 + pow2(Lt(m,n))) + (-(pfS(m,n) * pfv(m,n))
    /*4*/  + pftau(m,n)) * pow2(Lt(m,n))) * P_1_2(R(m,n)) + 4 * fA(m,n) * pow2(p(m,n))
    /*2*/  * P_1_2(R(m,n)) * (-pfD(m,n) - pftau(m,n) + fA(m,n) * Lt(m,n) 
    /*3*/  * pow2(gB(m,n)) * P_2_1(R(m,n))) + 2 * Lt(m,n) * pow2(gA(m,n)) 
    /*2*/  * pow2(gB(m,n)) * (-2 * P_1_1(R(m,n)) * P_2_0(R(m,n)) + (-P_1_0(R(m,n)) 
    /*4*/  + P_2_0(R(m,n))) * P_2_1(R(m,n))) + gA(m,n) * (Lt(m,n) * (pfD(m,n) - pfS(m,n)
    /*4*/  * pfv(m,n) + pftau(m,n)) * (2 * P_1_1(R(m,n)) - P_2_1(R(m,n))) + fA(m,n) 
    /*3*/  * pow2(gB(m,n)) * (-2 * P_1_2(R(m,n)) * P_2_0(R(m,n)) + P_2_1(R(m,n)) * (-4 
    /*5*/  * P_1_1(R(m,n)) + P_2_1(R(m,n)))))) + k_f * pow2(fA(m,n)) * pow2(gA(m,n)) 
    /*1*/  * pow2(gB(m,n)) * pow2(Lt(m,n)) * P_2_1(R(m,n)) * (6 * gA(m,n) * Lt(m,n) 
    /*2*/  * P_1_1(R(m,n)) * P_2_1(R(m,n)) * R(m,n) + fA(m,n) * (P_2_1(R(m,n)) 
    /*3*/  * (P_2_1(R(m,n)) + 2 * P_1_2(R(m,n)) * R(m,n)) - 2 * P_1_1(R(m,n)) 
    /*3*/  * (P_2_1(R(m,n)) - P_2_2(R(m,n)) * R(m,n)))) + (2 * gA(m,n) * Lt(m,n) 
    /*2*/  * pow2(fA(m,n)) * pow2(gB(m,n)) * pow2(R(m,n)) * (-12 * fA(m,n) 
    /*3*/  * (-gKD_r(m,n) + gK_r(m,n)) * p(m,n) * pow2(Lt(m,n)) * P_1_2(R(m,n)) 
    /*3*/  * P_2_1(R(m,n)) - 2 * fA(m,n) * gA(m,n) * Lt(m,n) * pow2(gK(m,n) - gKD(m,n))
    /*3*/  * pow2(p(m,n)) * P_1_2(R(m,n)) * (2 * P_1_1(R(m,n)) + P_2_1(R(m,n))) + 2 
    /*3*/  * (fK(m,n) - fKD(m,n)) * Lt(m,n) * pow2(gA(m,n)) * pow2(p(m,n)) 
    /*3*/  * P_1_1(R(m,n)) * (2 * (gK(m,n) - gKD(m,n)) * P_1_1(R(m,n)) + (gK(m,n) + 2 
    /*5*/  * gKD(m,n)) * P_2_1(R(m,n))) * R(m,n) + pow2(gA(m,n)) * pow2(Lt(m,n)) 
    /*3*/  * P_2_1(R(m,n)) * (-2 * (fK(m,n) - fKD(m,n)) * (gK(m,n) + 2 * gKD(m,n)) 
    /*4*/  * Lt(m,n) * P_1_1(R(m,n)) * R(m,n) + pow2(gK(m,n) - gKD(m,n)) * (-4 
    /*5*/  * P_1_1(R(m,n)) + P_2_1(R(m,n)) - 2 * P_0_2(R(m,n)) * R(m,n)) + 2 * (gK(m,n)
    /*5*/  - gKD(m,n)) * R(m,n) * (-((gK(m,n) + 2 * gKD(m,n)) * P_1_2(R(m,n))) 
    /*5*/  + (fK(m,n) - fKD(m,n)) * Lt(m,n) * (P_1_1(R(m,n)) + P_0_2(R(m,n)) * R(m,n))))
    /*3*/  + gA(m,n) * P_2_1(R(m,n)) * (6 * ((fK(m,n) - fKD(m,n)) * P_1_1(R(m,n)) 
    /*5*/  + (gK(m,n) - gKD(m,n)) * Lt(m,n) * P_1_2(R(m,n))) * p_r(m,n) * R(m,n) 
    /*4*/  + fA(m,n) * Lt(m,n) * (3 * (fK(m,n) - fKD(m,n)) * (fK(m,n) + fKD(m,n)) 
    /*5*/  * P_1_1(R(m,n)) * R(m,n) + 2 * (gK(m,n) - gKD(m,n)) * Lt(m,n) * (fKD(m,n) 
    /*6*/  * P_0_2(R(m,n)) - fK(m,n) * (P_0_2(R(m,n)) - 3 * P_1_2(R(m,n)))) * R(m,n) 
    /*5*/  - pow2(gK(m,n) - gKD(m,n)) * (3 * P_1_2(R(m,n)) + 2 * pow2(Lt(m,n)) 
    /*6*/  * P_0_3(R(m,n)) * R(m,n)))))) / 9.) + 2 * pow2(Lt(m,n)) * R(m,n) * (4 
    /*1*/  * Lt(m,n) * pow2(fDB(m,n)) * pow2(gB(m,n)) * pow2(P_1_1(R(m,n))) 
    /*1*/  * pow3(gA(m,n)) * pow3(R(m,n)) + pow3(fA(m,n)) * (pow2(gA(m,n)) 
    /*2*/  * P_2_1(R(m,n)) * (P_1_1(R(m,n)) - P_1_2(R(m,n)) * R(m,n)) + pow2(gB(m,n)) 
    /*2*/  * pow2(gDB(m,n)) * R(m,n) * (-4 * pow2(Lt(m,n)) * P_1_1(R(m,n)) 
    /*3*/  * P_1_2(R(m,n)) + P_2_1(R(m,n)) * ((1 - 2 * pow2(p(m,n))) * P_1_2(R(m,n)) - 2
    /*4*/  * pow2(p(m,n)) * P_0_3(R(m,n)) * R(m,n)))) - fA(m,n) * fDB(m,n) 
    /*1*/  * pow2(gA(m,n)) * pow2(gB(m,n)) * pow2(R(m,n)) * (fDB(m,n) * (-2 
    /*3*/  * P_0_2(R(m,n)) * P_2_1(R(m,n)) * R(m,n) + P_1_1(R(m,n)) * (P_2_1(R(m,n)) + 4
    /*4*/  * P_1_2(R(m,n)) * R(m,n))) + 2 * gDB(m,n) * (-2 * P_1_1(R(m,n)) 
    /*3*/  * P_1_2(R(m,n)) * R(m,n) + pow2(Lt(m,n)) * (2 * pow2(P_1_1(R(m,n))) 
    /*4*/  + P_1_1(R(m,n)) * P_2_1(R(m,n)) + P_0_2(R(m,n)) * P_2_1(R(m,n)) * R(m,n)))) 
    /*1*/  + gA(m,n) * Lt(m,n) * pow2(fA(m,n)) * R(m,n) * (pow2(gA(m,n)) * P_2_1(R(m,n))
    /*2*/  * (-2 * P_1_1(R(m,n)) + P_2_1(R(m,n))) + gDB(m,n) * pow2(gB(m,n)) * (2 
    /*3*/  * fDB(m,n) * (4 * P_1_1(R(m,n)) * P_1_2(R(m,n)) - P_0_2(R(m,n)) 
    /*4*/  * P_2_1(R(m,n))) * R(m,n) + gDB(m,n) * (-(P_2_1(R(m,n)) * (P_2_1(R(m,n)) - 2
    /*6*/  * P_0_2(R(m,n)) * R(m,n))) + 4 * P_1_1(R(m,n)) * (P_2_1(R(m,n)) 
    /*5*/  - P_1_2(R(m,n)) * R(m,n))))));
}
Real BimetricEvolve::eq_invr_fK1_t ( Int  m,
Int  n 
) [private]

Definition at line 130 of file eomEvolution.h.

{
    return (2 * fAlp(m,n) * (fDA(m,n) - 2 * fDB(m,n))) / pow2(fA(m,n));
}
Real BimetricEvolve::eq_invr_fK2_t ( Int  m,
Int  n 
) [private]

Definition at line 144 of file eomEvolution.h.

{
    return fAlp(m,n) * (fSig(m,n) / (fA(m,n) * fB(m,n)) + (fDA(m,n) - fDAlp(m,n) - 4 
    /*2*/  * fDB(m,n) + fSig(m,n)) / pow2(fA(m,n)));
}
Real BimetricEvolve::eq_invr_fK_t ( Int  m,
Int  n 
) [private]

Definition at line 181 of file eomEvolution.h.

{
    return fAlp(m,n) * ((2 * fSig(m,n)) / (fA(m,n) * fB(m,n)) + (2 * (2 * fDA(m,n) 
    /*3*/  - fDAlp(m,n) - 6 * fDB(m,n) + fSig(m,n))) / pow2(fA(m,n)));
}
Real BimetricEvolve::eq_invr_fKD_t ( Int  m,
Int  n 
) [private]

Definition at line 196 of file eomEvolution.h.

{
    return fAlp(m,n) * (-(fSig(m,n) / (fA(m,n) * fB(m,n))) + (fDA(m,n) + fDAlp(m,n) 
    /*2*/  - fSig(m,n)) / pow2(fA(m,n)));
}
Real BimetricEvolve::eq_invr_gK1_t ( Int  m,
Int  n 
) [private]

Definition at line 22 of file eomEvolution.h.

{
    return (2 * gAlp(m,n) * (gDA(m,n) - 2 * gDB(m,n))) / pow2(gA(m,n));
}
Real BimetricEvolve::eq_invr_gK2_t ( Int  m,
Int  n 
) [private]

Definition at line 36 of file eomEvolution.h.

{
    return gAlp(m,n) * (gSig(m,n) / (gA(m,n) * gB(m,n)) + (gDA(m,n) - gDAlp(m,n) - 4 
    /*2*/  * gDB(m,n) + gSig(m,n)) / pow2(gA(m,n)));
}
Real BimetricEvolve::eq_invr_gK_t ( Int  m,
Int  n 
) [private]

Definition at line 73 of file eomEvolution.h.

{
    return gAlp(m,n) * ((2 * gSig(m,n)) / (gA(m,n) * gB(m,n)) + (2 * (2 * gDA(m,n) 
    /*3*/  - gDAlp(m,n) - 6 * gDB(m,n) + gSig(m,n))) / pow2(gA(m,n)));
}
Real BimetricEvolve::eq_invr_gKD_t ( Int  m,
Int  n 
) [private]

Definition at line 88 of file eomEvolution.h.

{
    return gAlp(m,n) * (-(gSig(m,n) / (gA(m,n) * gB(m,n))) + (gDA(m,n) + gDAlp(m,n) 
    /*2*/  - gSig(m,n)) / pow2(gA(m,n)));
}
Real BimetricEvolve::eq_invr_p_f ( Int  m,
Int  n 
) [private]

Definition at line 17 of file eomMiscEquations.h.

{
    return (2 * fKD(m,n) * pow2(fB(m,n))) / (k_f * gA(m,n) * pow2(gB(m,n)) 
    /*1*/  * P_2_1(R(m,n)));
}
Real BimetricEvolve::eq_invr_p_g ( Int  m,
Int  n 
) [private]

Definition at line 5 of file eomMiscEquations.h.

{
    return (-2 * gKD(m,n)) / (k_g * fA(m,n) * P_2_1(R(m,n)));
}
Real BimetricEvolve::eq_j_over_p ( Int  m,
Int  n 
) [private]
Real BimetricEvolve::eq_p1_t ( Int  m,
Int  n 
) [private]
Real BimetricEvolve::eq_p_f ( Int  m,
Int  n 
) [inline, private]

Definition at line 427 of file bim-solver.cpp.

                                {
        return eq_base_p_f(m,n) + eq_invr_p_f(m,n) / r(m,n);
    }
Real BimetricEvolve::eq_p_g ( Int  m,
Int  n 
) [inline, private]

Definition at line 420 of file bim-solver.cpp.

                                {
        return eq_base_p_g(m,n) + eq_invr_p_g(m,n) / r(m,n);
    }
Real BimetricEvolve::eq_p_t ( Int  m,
Int  n 
) [private]

Definition at line 29 of file eomMiscEquations.h.

{
    return eq_pr(m,n) * p(m,n) * ((-2 * fB(m,n) * gAlp(m,n) * P_1_1(R(m,n))) 
    /*1*/  / (fA(m,n) * gA(m,n) * (1 + Lt(m,n)) * P_2_1(R(m,n))) + (2 * fAlp(m,n) 
    /*2*/  * fB(m,n) * P_1_2(R(m,n))) / (fA(m,n) * gA(m,n) * (1 + Lt(m,n)) 
    /*2*/  * P_2_1(R(m,n)))) + (fAlp(m,n) * ((-2 * gDB(m,n) * P_1_2(R(m,n))) 
    /*2*/  / P_2_1(R(m,n)) - (2 * fB(m,n) * (fSig(m,n) - gSig(m,n) * Lt(m,n)) 
    /*3*/  * P_1_2(R(m,n))) / (fA(m,n) * P_2_1(R(m,n))))) / gA(m,n) + p(m,n) 
    /*0*/  * ((gAlp(m,n) * (gK(m,n) + 2 * gKD(m,n) + (2 * (gK(m,n) - gKD(m,n)) 
    /*4*/  * P_1_1(R(m,n))) / P_2_1(R(m,n)))) / 3. + (fAlp(m,n) * (2 * fKD(m,n) 
    /*3*/  * P_1_1(R(m,n)) + fK(m,n) * (-2 * P_1_1(R(m,n)) + 3 * P_2_1(R(m,n))))) / (3.
    /*2*/  * P_2_1(R(m,n)))) + p_r(m,n) * q(m,n) + gAlp(m,n) * ((Lt(m,n) 
    /*2*/  * (-gDAlp(m,n) - (2 * gDB(m,n) * P_1_1(R(m,n))) / P_2_1(R(m,n))) + (2 
    /*3*/  * fB(m,n) * (gSig(m,n) - fSig(m,n) * Lt(m,n)) * P_1_1(R(m,n))) / (fA(m,n) 
    /*3*/  * P_2_1(R(m,n)))) / gA(m,n) + (2 * fDB(m,n) * P_1_1(R(m,n)) * R(m,n)) 
    /*1*/  / (fA(m,n) * P_2_1(R(m,n)))) + (fAlp(m,n) * Lt(m,n) * (fDAlp(m,n) + (2 
    /*3*/  * fDB(m,n) * P_1_2(R(m,n)) * R(m,n)) / P_2_1(R(m,n)))) / fA(m,n);
}
Real BimetricEvolve::eq_pf_D ( Int  m,
Int  n 
) [private]
Real BimetricEvolve::eq_pf_D_t ( Int  m,
Int  n 
) [private]

Definition at line 6 of file eomMatter.h.

{
    return gBet_r(m,n) * pfD(m,n) + gBet(m,n) * pfD_r(m,n) + gAlp(m,n) * (-(gDAlp(m,n)
    /*2*/  * pfD(m,n) * pfv(m,n)) - pfD_r(m,n) * pfv(m,n) - pfD(m,n) * pfv_r(m,n)) 
    /*0*/  + (2 * gBet(m,n) * pfD(m,n) - 2 * gAlp(m,n) * pfD(m,n) * pfv(m,n)) / r(m,n);
}
Real BimetricEvolve::eq_pf_j ( Int  m,
Int  n 
) [private]

Definition at line 47 of file eomMatter.h.

{
    return pfS(m,n) / (gA(m,n) * pow2(gB(m,n)));
}
Real BimetricEvolve::eq_pf_J1 ( Int  m,
Int  n 
) [private]

Definition at line 52 of file eomMatter.h.

{
    return (pfS(m,n) * pfv(m,n)) / (gA(m,n) * pow2(gB(m,n)));
}
Real BimetricEvolve::eq_pf_J2 ( Int  m,
Int  n 
) [private]

Definition at line 57 of file eomMatter.h.

{
    return 0;
}
Real BimetricEvolve::eq_pf_rho ( Int  m,
Int  n 
) [private]

Definition at line 42 of file eomMatter.h.

{
    return (pfD(m,n) + pftau(m,n)) / (gA(m,n) * pow2(gB(m,n)));
}
Real BimetricEvolve::eq_pf_S ( Int  m,
Int  n 
) [private]
Real BimetricEvolve::eq_pf_S_t ( Int  m,
Int  n 
) [private]

Definition at line 13 of file eomMatter.h.

{
    return 2 * gBet_r(m,n) * pfS(m,n) + gBet(m,n) * pfS_r(m,n) + gAlp(m,n) * (gDA(m,n)
    /*1*/  * pfS(m,n) * pfv(m,n) - pfS_r(m,n) * pfv(m,n) - pfS(m,n) * pfv_r(m,n) 
    /*1*/  - gDAlp(m,n) * (pfD(m,n) + pfS(m,n) * pfv(m,n) + pftau(m,n))) + (2 
    /*1*/  * gBet(m,n) * pfS(m,n) - 2 * gAlp(m,n) * pfS(m,n) * pfv(m,n)) / r(m,n);
}
Real BimetricEvolve::eq_pf_tau ( Int  m,
Int  n 
) [private]
Real BimetricEvolve::eq_pf_tau_t ( Int  m,
Int  n 
) [private]

Definition at line 21 of file eomMatter.h.

{
    return gBet_r(m,n) * pftau(m,n) + gBet(m,n) * pftau_r(m,n) + gAlp(m,n) * (gK1(m,n)
    /*1*/  * pfS(m,n) * pfv(m,n) - gDAlp(m,n) * pfv(m,n) * pftau(m,n) - pfv_r(m,n) 
    /*1*/  * pftau(m,n) - pfv(m,n) * pftau_r(m,n) - (gDAlp(m,n) * pfS(m,n)) 
    /*1*/  / pow2(gA(m,n))) + (2 * gBet(m,n) * pftau(m,n) - 2 * gAlp(m,n) * pfv(m,n) 
    /*1*/  * pftau(m,n)) / r(m,n);
}
Real BimetricEvolve::eq_pf_v ( Int  m,
Int  n 
) [private]

Definition at line 30 of file eomMatter.h.

{
    return pfS(m,n) / ((TINY_Real + pfD(m,n) + pftau(m,n)) * pow2(gA(m,n)));
}
Real BimetricEvolve::eq_pf_v_r ( Int  m,
Int  n 
) [private]

Definition at line 35 of file eomMatter.h.

{
    return (-2 * gDA(m,n) * pfS(m,n) * (pfD(m,n) + pftau(m,n)) + pfS_r(m,n) 
    /*1*/  * (pfD(m,n) + pftau(m,n)) - pfS(m,n) * (pfD_r(m,n) + pftau_r(m,n))) 
    /*0*/  / (pow2(gA(m,n)) * pow2(TINY_Real + pfD(m,n) + pftau(m,n)));
}
Real BimetricEvolve::eq_pr_r ( Real  m,
Real  n 
) [private]
Real BimetricEvolve::eq_q ( Int  m,
Int  n 
) [private]

Definition at line 6 of file eomGauge.h.

{
    return 0;
}
Real BimetricEvolve::eq_qr_r ( Real  m,
Real  n 
) [private]
Real BimetricEvolve::fA_r ( Real  m,
Real  n 
) [private]
Real BimetricEvolve::fAlp_r ( Real  m,
Real  n 
) [private]
Real BimetricEvolve::fAlp_rr ( Real  m,
Real  n 
) [private]
Real BimetricEvolve::fB_r ( Real  m,
Real  n 
) [private]
Real BimetricEvolve::fDA_r ( Real  m,
Real  n 
) [private]
Real BimetricEvolve::fDAlp_r ( Real  m,
Real  n 
) [private]
Real BimetricEvolve::fDB_r ( Real  m,
Real  n 
) [private]
Real BimetricEvolve::fK1 ( Int  m,
Int  n 
) [inline, private]

Definition at line 320 of file bim-solver.cpp.

{ return ( fK(m,n) + 2 * fKD(m,n) ) / 3; }
Real BimetricEvolve::fK2 ( Int  m,
Int  n 
) [inline, private]

Definition at line 321 of file bim-solver.cpp.

{ return ( fK(m,n) - fKD(m,n) ) / 3; }
Real BimetricEvolve::fK_r ( Real  m,
Real  n 
) [private]
Real BimetricEvolve::fKD_r ( Real  m,
Real  n 
) [private]
Real BimetricEvolve::gA_r ( Real  m,
Real  n 
) [private]
Real BimetricEvolve::gAlp_r ( Real  m,
Real  n 
) [private]
Real BimetricEvolve::gB_r ( Real  m,
Real  n 
) [private]
Real BimetricEvolve::gDA_r ( Real  m,
Real  n 
) [private]
Real BimetricEvolve::gDAlp_r ( Real  m,
Real  n 
) [private]
Real BimetricEvolve::gDB_r ( Real  m,
Real  n 
) [private]
Real BimetricEvolve::gK1 ( Int  m,
Int  n 
) [inline, private]

Definition at line 318 of file bim-solver.cpp.

{ return ( gK(m,n) + 2 * gKD(m,n) ) / 3; }
Real BimetricEvolve::gK2 ( Int  m,
Int  n 
) [inline, private]

Definition at line 319 of file bim-solver.cpp.

{ return ( gK(m,n) - gKD(m,n) ) / 3; }
Real BimetricEvolve::gK_r ( Real  m,
Real  n 
) [private]
Real BimetricEvolve::gKD_r ( Real  m,
Real  n 
) [private]
Real BimetricEvolve::p_r ( Real  m,
Real  n 
) [private]
Real BimetricEvolve::p_rr ( Real  m,
Real  n 
) [private]
Real BimetricEvolve::pfD_r ( Real  m,
Real  n 
) [private]
Real BimetricEvolve::pfS_r ( Real  m,
Real  n 
) [private]
Real BimetricEvolve::pftau_r ( Real  m,
Real  n 
) [private]
Real BimetricEvolve::q_r ( Real  m,
Real  n 
) [private]
Real BimetricEvolve::q_rr ( Real  m,
Real  n 
) [private]