28char strot_dirac_equilibrium_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Star/strot_dirac_equilibrium.C,v 1.13 2014/10/13 08:53:40 j_novak Exp $" ;
89#include "star_rot_dirac.h"
91#include "utilitaires.h"
120 assert( ( type_t == SYM) || (type_t == NONSYM) ) ;
142 cout <<
"Star_rot_Dirac::equilibrium: mer_change_omega < mer_rot !"
149 cout <<
"Star_rot_Dirac::equilibrium: mer_fix_omega < mer_change_omega !"
206 fichconv <<
"# diff_ent GRV2 max_triax vit_triax" <<
endl ;
213 "# |dH/dr_eq/dH/dr_pole| r_pole/r_eq ent_c"
227 cout <<
"-----------------------------------------------" <<
endl ;
329 cout <<
"U > c for l, i : " <<
l <<
" " <<
i
335 cout <<
"**** VELOCITY OF LIGHT REACHED ****" <<
endl ;
337 cout <<
"New rotation frequency : "
407 cout <<
"ent < 0 for l, i : " <<
l <<
" " <<
i
414 cout <<
"**** KEPLERIAN VELOCITY REACHED ****" <<
endl ;
416 cout <<
"New rotation frequency : "
447 q_new.std_spectral_base() ;
507 cout <<
"Discrep. baryon mass <-> wanted bar. mass : " <<
xx
512 cout <<
"Discrep. grav. mass <-> wanted grav. mass : " <<
xx
518 double ax = .5 * ( 2. +
xx ) / (1. +
xx ) ;
520 cout <<
" xprog, xx, ax, fact : " <<
xprog <<
" " <<
Time evolution with partial storage (*** under development ***).
Basic integer array class.
virtual void homothetie(double lambda)=0
Sets a new radial scale.
const Base_vect_spher & get_bvect_spher() const
Returns the orthonormal vectorial basis associated with the coordinates of the mapping.
const Mg3d * get_mg() const
Gives the Mg3d on which the mapping is defined.
virtual const Sym_tensor & cov() const
Read-only access to the covariant representation.
int get_type_t() const
Returns the type of sampling in the direction: SYM : : symmetry with respect to the equatorial pl...
int get_nt(int l) const
Returns the number of points in the co-latitude direction ( ) in domain no. l.
int get_nzone() const
Returns the number of domains.
int get_nr(int l) const
Returns the number of points in the radial direction ( ) in domain no. l.
Tensor field of valence 0 (or component of a tensorial field).
virtual void std_spectral_base()
Sets the spectral bases of the Valeur va to the standard ones for a scalar field.
double val_grid_point(int l, int k, int j, int i) const
Returns the value of the field at a specified grid point.
virtual void set_etat_zero()
Sets the logical state to ETATZERO (zero).
virtual void annule(int l_min, int l_max)
Sets the Scalar to zero in several domains.
void mult_rsint()
Multiplication by everywhere; dzpuis is not changed.
virtual double mass_g() const
Gravitational mass.
Sym_tensor_trans hh
is defined by .
double omega
Rotation angular velocity ([f_unit] )
const Metric_flat & flat
flat metric (spherical components)
void solve_logn_f(Scalar &ln_f_new) const
Solution of the two scalar Poisson equations for rotating stars in Dirac gauge.
void solve_qqq(Scalar &q_new) const
Solution of the two scalar Poisson equations for rotating stars in Dirac gauge.
void update_metric()
Computes metric quantities from known potentials.
void solve_hij(Sym_tensor_trans &hij_new) const
Solution of the tensor Poisson equation for rotating stars in Dirac gauge.
virtual void hydro_euler()
Computes the hydrodynamical quantities relative to the Eulerian observer from those in the fluid fram...
void solve_logn_q(Scalar &ln_q_new) const
Solution of the two scalar Poisson equations for rotating stars in Dirac gauge.
virtual double mass_b() const
Baryonic mass.
void equilibrium(double ent_c, double omega0, double fact_omega, int nzadapt, const Tbl &ent_limit, const Itbl &icontrol, const Tbl &control, double mbar_wanted, double aexp_mass, Tbl &diff)
Computes an equilibrium configuration.
void solve_shift(Vector &shift_new) const
Solution of the shift equation for rotating stars in Dirac gauge.
virtual double grv2() const
Error on the virial identity GRV2.
Scalar logn
Logarithm of the lapse N .
Scalar nn
Lapse function N .
void equation_of_state()
Computes the proper baryon and energy density, as well as pressure from the enthalpy.
Scalar gam_euler
Lorentz factor between the fluid and Eulerian observers.
Vector u_euler
Fluid 3-velocity with respect to the Eulerian observer.
Map & mp
Mapping associated with the star.
int nzet
Number of domains of *mp occupied by the star.
Transverse symmetric tensors of rank 2.
Tensor field of valence 1.
Scalar & set(int)
Read/write access to a component.
Cmp sqrt(const Cmp &)
Square root.
Cmp log10(const Cmp &)
Basis 10 logarithm.
Tbl diffrel(const Cmp &a, const Cmp &b)
Relative difference between two Cmp (norme version).
Cmp pow(const Cmp &, int)
Power .
Cmp abs(const Cmp &)
Absolute value.
Cmp log(const Cmp &)
Neperian logarithm.
Tenseur contract(const Tenseur &, int id1, int id2)
Self contraction of two indices of a Tenseur .
Standard units of space, time and mass.