_priv_/bim-ss/eom_sources_8h_source.html

00001 /** @file  eomSources.h
00002  *  @brief The stress-energy tensor components.
00003  *  @image html sources.png
00004  */
00005
00006 Real BimetricEvolve::eq_g_rho( Int m, Int n )
00007 {
00008     return (pfD(m,n) + pftau(m,n)) / (gA(m,n) * pow2(gB(m,n))) - P_2_0(R(m,n))
00009     /*0*/  - (fA(m,n) * Lt(m,n) * P_2_1(R(m,n))) / gA(m,n);
00010 }
00011
00012 Real BimetricEvolve::eq_g_j( Int m, Int n )
00013 {
00014     return pfS(m,n) / (gA(m,n) * pow2(gB(m,n))) - fA(m,n) * p(m,n) * P_2_1(R(m,n));
00015 }
00016
00017 Real BimetricEvolve::eq_g_JK1( Int m, Int n )
00018 {
00019     return (pfD(m,n) + pfS(m,n) * pfv(m,n) + pftau(m,n)) / (2. * gA(m,n)
00020     /*1*/  * pow2(gB(m,n))) - P_1_0(R(m,n)) + (fA(m,n) * (-2 * P_1_1(R(m,n)) + (1 - 2
00021     /*3*/  * pow2(Lt(m,n))) * P_2_1(R(m,n)))) / (2. * gA(m,n) * Lt(m,n)) + (-((fA(m,n)
00022     /*3*/  * fAlp(m,n) * P_1_2(R(m,n))) / gA(m,n)) + (fAlp(m,n) * (-2 * P_1_1(R(m,n))
00023     /*3*/  + P_2_1(R(m,n)))) / (2. * Lt(m,n))) / gAlp(m,n);
00024 }
00025
00026 Real BimetricEvolve::eq_g_JK2( Int m, Int n )
00027 {
00028     return (pfD(m,n) - pfS(m,n) * pfv(m,n) + pftau(m,n)) / (2. * gA(m,n)
00029     /*1*/  * pow2(gB(m,n))) - P_2_0(R(m,n)) - (fA(m,n) * P_2_1(R(m,n))) / (2. * gA(m,n)
00030     /*1*/  * Lt(m,n)) - (fAlp(m,n) * P_2_1(R(m,n))) / (2. * gAlp(m,n) * Lt(m,n));
00031 }
00032
00033 Real BimetricEvolve::eq_f_rho( Int m, Int n )
00034 {
00035     return -((gA(m,n) * Lt(m,n) * P_2_1(R(m,n))) / (fA(m,n) * pow2(R(m,n))))
00036     /*0*/  - P_2_2(R(m,n)) / pow2(R(m,n));
00037 }
00038
00039 Real BimetricEvolve::eq_f_j( Int m, Int n )
00040 {
00041     return (gA(m,n) * p(m,n) * pow2(gB(m,n)) * P_2_1(R(m,n))) / pow2(fB(m,n));
00042 }
00043
00044 Real BimetricEvolve::eq_f_JK1( Int m, Int n )
00045 {
00046     return -(P_1_3(R(m,n)) / R(m,n)) + (gA(m,n) * ((1 - 2 * pow2(Lt(m,n)))
00047     /*2*/  * P_2_1(R(m,n)) - 2 * P_1_2(R(m,n)) * R(m,n))) / (2. * fA(m,n) * Lt(m,n)
00048     /*1*/  * pow2(R(m,n))) + gAlp(m,n) * (-((gA(m,n) * P_1_1(R(m,n))) / (fA(m,n)
00049     /*3*/  * fAlp(m,n) * R(m,n))) + (P_2_1(R(m,n)) - 2 * P_1_2(R(m,n)) * R(m,n)) / (2.
00050     /*2*/  * fAlp(m,n) * Lt(m,n) * pow2(R(m,n))));
00051 }
00052
00053 Real BimetricEvolve::eq_f_JK2( Int m, Int n )
00054 {
00055     return -(gA(m,n) * P_2_1(R(m,n))) / (2. * fA(m,n) * Lt(m,n) * pow2(R(m,n)))
00056     /*0*/  - (gAlp(m,n) * P_2_1(R(m,n))) / (2. * fAlp(m,n) * Lt(m,n) * pow2(R(m,n)))
00057     /*0*/  - P_2_2(R(m,n)) / pow2(R(m,n));
00058 }