28char star_rot_dirac_diff_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Star/star_rot_dirac_diff.C,v 1.4 2014/10/13 08:53:39 j_novak Exp $" ;
41#include "star_rot_dirac_diff.h"
43#include "utilitaires.h"
79 primfrot(star.primfrot),
80 par_frot(star.par_frot),
81 omega_field(star.omega_field),
82 omega_min(star.omega_min),
83 omega_max(star.omega_max),
84 prim_field(star.prim_field)
176 ost <<
"Differentially rotating star in Dirac gauge" <<
'\n';
180 ost <<
"Differentially rotating star" <<
'\n';
181 ost <<
"-----------------------" <<
'\n';
183 ost <<
'\n'<<
"Parameters of F(Omega) : " <<
'\n';
187 <<
" Hz, " <<
omega_max / (2*
M_PI) * f_unit <<
" Hz" <<
'\n';
191 ost <<
"Central, equatorial value of Omega: "
197 ost <<
"Central, equatorial value of Omega/(2 Pi): "
204 ost <<
"Error on the virial identity GRV2 : " <<
'\n';
205 ost <<
"GRV2 = " <<
grv2() <<
'\n';
206 ost <<
"Error on the virial identity GRV3 : " <<
'\n';
207 ost <<
"GRV3 = " <<
grv3() <<
'\n';
209 ost <<
"Angular momentum J : "
210 <<
angu_mom()/( qpig / (4*
M_PI) *msol*msol) <<
" G M_sol^2 / c"
212 ost <<
"c J / (G M^2) : "
223 ost <<
"Ratio T/W : " <<
tsw() <<
'\n';
227 ost <<
"Circumferential equatorial radius R_circ : "
228 <<
r_circ()/km <<
" km" <<
'\n';
229 ost <<
"Circumferential polar radius Rp_circ : "
231 ost <<
"Coordinate equatorial radius r_eq : " <<
ray_eq()/km <<
" km"
234 ost <<
"Ellipticity sqrt(1-(Rp_circ/R_circ)^2) : " <<
ellipt() <<
endl ;
236 ost <<
"Compaction parameter M_g / R_circ : " <<
compact <<
'\n';
Equation of state base class.
Time evolution with partial storage (*** under development ***).
Base class for coordinate mappings.
const Mg3d * get_mg() const
Gives the Mg3d on which the mapping is defined.
Tensor field of valence 0 (or component of a tensorial field).
virtual void sauve(FILE *) const
Save in a file.
double val_grid_point(int l, int k, int j, int i) const
Returns the value of the field at a specified grid point.
Class for relativistic differentially rotating stars in Dirac gauge and maximal slicing.
virtual ~Star_rot_Dirac_diff()
Destructor.
Scalar omega_field
Field .
virtual ostream & operator>>(ostream &) const
Operator >> (virtual function called by the operator <<).
double funct_omega(double omeg) const
Evaluates , where F is the function defining the rotation profile.
virtual void sauve(FILE *) const
Save in a file.
double prim_funct_omega(double omeg) const
Evaluates the primitive of , where F is the function defining the rotation profile.
Star_rot_Dirac_diff(Map &mp_i, int nzet_i, const Eos &eos_i, double(*frot_i)(double, const Tbl &), double(*primfrot_i)(double, const Tbl &), const Tbl &par_frot_i)
Standard constructor.
void fait_prim_field()
Computes the member prim_field from omega_field .
double(* frot)(double, const Tbl &)
Function defining the rotation profile.
double omega_min
Minimum value of .
virtual double get_omega_c() const
Returns the central value of the rotation angular velocity ([f_unit] )
double(* primfrot)(double, const Tbl &)
Primitive of the function , which vanishes at the stellar center.
virtual double tsw() const
Ratio T/W.
Tbl par_frot
Parameters of the function .
void operator=(const Star_rot_Dirac_diff &)
Assignment to another Star_rot_Dirac_diff.
double omega_max
Maximum value of .
Class for relativistic rotating stars in Dirac gauge and maximal slicing.
virtual double mass_g() const
Gravitational mass.
virtual double ellipt() const
Ellipticity e.
virtual double grv3() const
Error on the virial identity GRV3.
virtual double angu_mom() const
Angular momentum.
virtual double r_circ() const
Circumferential equatorial radius.
double omega
Rotation angular velocity ([f_unit] )
virtual double rp_circ() const
Circumferential polar radius.
virtual double aplat() const
Flattening r_pole/r_eq.
void operator=(const Star_rot_Dirac &)
Assignment to another Star_rot_Dirac.
virtual double grv2() const
Error on the virial identity GRV2.
virtual ostream & operator>>(ostream &) const
Operator >> (virtual function called by the operator <<).
double ray_eq() const
Coordinate radius at , [r_unit].
virtual void sauve(FILE *) const
Save in a file.
Map & mp
Mapping associated with the star.
int nzet
Number of domains of *mp occupied by the star.
void sauve(FILE *) const
Save in a file.
Cmp pow(const Cmp &, int)
Power .
int fread_be(int *aa, int size, int nb, FILE *fich)
Reads integer(s) from a binary file according to the big endian convention.
int fwrite_be(const int *aa, int size, int nb, FILE *fich)
Writes integer(s) into a binary file according to the big endian convention.
Standard units of space, time and mass.