28char kerr_QI_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Compobj/kerr_QI.C,v 1.4 2014/10/13 08:52:49 j_novak Exp $" ;
72 double hh =
sqrt(mm2-aa2) ;
74 double r_hor = hh / double(2) ;
80 Mtbl cost2 = cost*cost ;
81 Mtbl sint2 = sint*sint ;
84 a_car = 1 + 2*
mm/r + (3*mm2 + aa2*(2*cost2-1))/(2*r2) + hh2*
mm/(2*r*r2) +
90 Mtbl rBL = r + hh2/(4*r) +
mm ;
92 Mtbl sigma = rBL2 + aa2*cost2 ;
93 Mtbl rBLovr = 1 +
mm/r + r_hor*r_hor/r2 ;
96 b_car = rBLovr * ( rBLovr + 2*aa2*
mm*sint2 / (r*sigma) ) + aa2 / r2 ;
103 nn = (1 - r_hor*r_hor / r2) /
bbb ;
108 nphi = 2*
aa*
mm / (sigma*(rBL+aa2/rBL) + 2*aa2*
mm*sint2) ;
209 ost << endl <<
"Kerr spacetime in quasi-isotropic coordinates (class Kerr_QI) " << endl ;
211 ost <<
"M = " <<
mm <<
" a = " <<
aa << endl ;
Base class for axisymmetric stationary compact objects in Quasi-Isotropic coordinates (under developm...
void operator=(const Compobj_QI &)
Assignment to another Compobj_QI.
Scalar nphi
Metric coefficient .
virtual void del_deriv() const
Deletes all the derived quantities.
virtual ostream & operator>>(ostream &) const
Operator >> (virtual function called by the operator <<).
Scalar b_car
Square of the metric factor B.
Scalar bbb
Metric factor B.
Scalar a_car
Square of the metric factor A.
Scalar nn
Lapse function N .
Map & mp
Mapping describing the coordinate system (r,theta,phi)
Active physical coordinates and mapping derivatives.
Kerr spacetime in Quasi-Isotropic coordinates (under development).
virtual ~Kerr_QI()
Destructor.
virtual void del_deriv() const
Deletes all the derived quantities.
Kerr_QI(Map &mp_i, double mass, double a_over_m)
Standard constructor.
void operator=(const Kerr_QI &)
Assignment to another Kerr_QI.
double aa
angular momentum parameter
virtual void set_der_0x0() const
Sets to 0x0 all the pointers on derived quantities.
virtual void sauve(FILE *) const
Save in a file.
virtual ostream & operator>>(ostream &) const
Operator >> (virtual function called by the operator <<).
Base class for coordinate mappings.
Coord r
r coordinate centered on the grid
virtual void std_spectral_base()
Sets the spectral bases of the Valeur va to the standard ones for a scalar field.
Tbl & set_domain(int l)
Read/write of the value in a given domain.
int get_etat() const
Returns the logical state ETATNONDEF (undefined), ETATZERO (null) or ETATQCQ (ordinary).
Cmp sqrt(const Cmp &)
Square root.
Standard units of space, time and mass.