 |
Bimetric 3+1 toolkit for spherical symmetry
|
|
Bimetric 3+1 toolkit for spherical symmetry
|
00001 /** 00002 * @file bimetricModel.h 00003 * @brief Implements the parameters which define bimetric models. 00004 * @author Mikica Kocic 00005 * @copyright GNU General Public License (GPLv3). 00006 */ 00007 00008 ///////////////////////////////////////////////////////////////////////////////////////// 00009 /** @defgroup g2 The bimetric model */ 00010 ///////////////////////////////////////////////////////////////////////////////////////// 00011 /*@{*/ 00012 #ifndef V_SIGN 00013 #define V_SIGN - //!< V_SIGN should be defined either as empty or as `-`. 00014 #endif 00015 00016 /** A structure that holds parameters which define a specific bimetric model. 00017 */ 00018 struct BimetricModel 00019 { 00020 std::string title; //!< Description 00021 00022 Real k_g; //!< kappa_g parameter 00023 Real k_f; //!< kappa_f parameter 00024 Real b_0; //!< beta_0 parameter 00025 Real b_1; //!< beta_1 parameter 00026 Real b_2; //!< beta_2 parameter 00027 Real b_3; //!< beta_3 parameter 00028 Real b_4; //!< beta_4 parameter 00029 00030 /** Constructs the bimetric model based on values from the parameter file. 00031 */ 00032 BimetricModel( Parameters& params ) 00033 { 00034 params.get( "model.title", title, "" ); 00035 params.get( "model.k_g", k_g, 1.0 ); 00036 params.get( "model.k_f", k_f, 1.0 ); 00037 params.get( "model.b_0", b_0, 0.0 ); 00038 params.get( "model.b_1", b_1, 0.0 ); 00039 params.get( "model.b_2", b_2, 0.0 ); 00040 params.get( "model.b_3", b_3, 0.0 ); 00041 params.get( "model.b_4", b_4, 0.0 ); 00042 00043 slog << "Bimetric Model: " << std::endl << std::endl 00044 << " title = " << title << std::endl 00045 << " k_g = " << k_g << ", k_f = " << k_f 00046 << ", b_0 = " << b_0 << ", b_1 = " << b_1 00047 << ", b_2 = " << b_2 << ", b_3 = " << b_3 00048 << ", b_4 = " << b_4 << std::endl << std::endl; 00049 } 00050 00051 /** Returns true if the bimetric model comprises two decoupled GR sectors. 00052 */ 00053 bool isGR () 00054 { 00055 return b_1 == 0 && b_2 == 0 && b_3 == 0; 00056 } 00057 00058 ///////////////////////////////////////////////////////////////////////////////////// 00059 /** @defgroup g3 Shifted elementary symmetric polynomials */ 00060 ///////////////////////////////////////////////////////////////////////////////////// 00061 /*@{*/ 00062 inline Real P_0_0( Real R ) { return V_SIGN b_0; } 00063 inline Real P_0_1( Real R ) { return V_SIGN b_1; } 00064 inline Real P_0_2( Real R ) { return V_SIGN b_2; } 00065 inline Real P_0_3( Real R ) { return V_SIGN b_3; } 00066 inline Real P_0_4( Real R ) { return V_SIGN b_4; } 00067 00068 inline Real P_1_0( Real R ) { return V_SIGN ( b_0 + b_1 * R ); } 00069 inline Real P_1_1( Real R ) { return V_SIGN ( b_1 + b_2 * R ); } 00070 inline Real P_1_2( Real R ) { return V_SIGN ( b_2 + b_3 * R ); } 00071 inline Real P_1_3( Real R ) { return V_SIGN ( b_3 + b_4 * R ); } 00072 00073 inline Real P_2_0( Real R ) { return V_SIGN ( b_0 + ( 2 * b_1 + b_2 * R ) * R ); } 00074 inline Real P_2_1( Real R ) { return V_SIGN ( b_1 + ( 2 * b_2 + b_3 * R ) * R ); } 00075 inline Real P_2_2( Real R ) { return V_SIGN ( b_2 + ( 2 * b_3 + b_4 * R ) * R ); } 00076 }; 00077 /*@}*/ 00078 /////////////////////////////////////////////////////////////////////////////////////////
1.8.0 
