28char star_rot_dirac_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Star/star_rot_dirac.C,v 1.10 2014/10/13 08:53:39 j_novak Exp $" ;
81#include "star_rot_dirac.h"
83#include "utilitaires.h"
100 j_euler(
mpi, CON,
mpi.get_bvect_spher()),
102 flat(
mpi.flat_met_spher()),
104 aa(
mpi, CON,
mpi.get_bvect_spher()),
105 taa(
mpi, COV,
mpi.get_bvect_spher()),
107 hh(
mpi,
mpi.get_bvect_spher(), flat)
139 spectral_filter(star.spectral_filter),
144 j_euler(star.j_euler),
150 aa_quad(star.aa_quad),
170 j_euler(
mpi, CON,
mpi.get_bvect_spher()),
172 flat(
mpi.flat_met_spher()),
174 aa(
mpi, CON,
mpi.get_bvect_spher()),
175 taa(
mpi, COV,
mpi.get_bvect_spher()),
321 ost <<
"Rotating star in Dirac gauge" <<
endl ;
325 ost <<
"Uniformly rotating star" <<
endl ;
326 ost <<
"-----------------------" <<
endl ;
328 ost <<
"hydro sources of equations are filtered\n"
333 <<
" rad/s f : " <<
freq * f_unit <<
" Hz" <<
endl ;
334 ost <<
"Rotation period : " << 1000. / (
freq * f_unit) <<
" ms"
337 ost <<
"Error on the virial identity GRV2 : " <<
endl ;
339 ost <<
"Error on the virial identity GRV3 : " <<
endl ;
342 ost <<
"Angular momentum J : "
343 <<
angu_mom()/( qpig / (4*
M_PI) *msol*msol) <<
" G M_sol^2 / c"
345 ost <<
"c J / (G M^2) : "
357 ost <<
"Circumferential equatorial radius R_circ : "
360 ost <<
"Circumferential polar radius Rp_circ : "
362 ost <<
"Coordinate equatorial radius r_eq : " <<
ray_eq()/km <<
" km"
366 ost <<
"Ellipticity sqrt(1-(Rp_circ/R_circ)^2) : " <<
ellipt() <<
endl ;
369 ost <<
"Compaction parameter M_g / R_circ : " <<
compact <<
endl ;
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.
virtual void sauve(FILE *) const
Save in a file.
virtual void set_etat_nondef()
Sets the logical state to ETATNONDEF (undefined).
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.
Star_rot_Dirac(Map &mp_i, int nzet_i, const Eos &eos_i, int filter=0)
Standard constructor.
Sym_tensor_trans hh
is defined by .
virtual double angu_mom() const
Angular momentum.
virtual ostream & operator>>(ostream &) const
Operator >> (virtual function called by the operator <<).
virtual void del_deriv() const
Deletes all the derived quantities.
double * p_grv3
Error on the virial identity GRV3.
virtual double r_circ() const
Circumferential equatorial radius.
int spectral_filter
Spectral exponential filtering order.
double omega
Rotation angular velocity ([f_unit] )
double * p_r_circ
Circumferential equatorial radius.
virtual double tsw() const
Ratio T/W.
Vector j_euler
Momentum density 3-vector with respect to the Eulerian observer.
const Metric_flat & flat
flat metric (spherical components)
virtual void del_hydro_euler()
Sets to ETATNONDEF (undefined state) the hydrodynamical quantities relative to the Eulerian observer.
virtual void sauve(FILE *) const
Save in a file.
void set_der_0x0() const
Sets to 0x0 all the pointers on derived quantities.
virtual double rp_circ() const
Circumferential polar radius.
void update_metric()
Computes metric quantities from known potentials.
Scalar psi4
Conformal factor .
virtual void hydro_euler()
Computes the hydrodynamical quantities relative to the Eulerian observer from those in the fluid fram...
double * p_grv2
Error on the virial identity GRV2.
virtual ~Star_rot_Dirac()
Destructor.
double * p_angu_mom
Angular momentum.
double * p_rp_circ
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 void del_hydro_euler()
Sets to ETATNONDEF (undefined state) the hydrodynamical quantities relative to the Eulerian observer.
void equation_of_state()
Computes the proper baryon and energy density, as well as pressure from the enthalpy.
virtual void del_deriv() const
Deletes all the derived quantities.
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.
void operator=(const Star &)
Assignment to another Star.
Tensor field of valence 1.
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.
virtual void sauve(FILE *) const
Save in a binary file.
virtual void sauve(FILE *) const
Save in a binary file.
virtual void set_etat_nondef()
Sets the logical state of all components to ETATNONDEF (undefined state).
virtual void set_etat_zero()
Sets the logical state of all components to ETATZERO (zero state).
Standard units of space, time and mass.