28char strot_dirac_diff_equil_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Star/strot_dirac_diff_equil.C,v 1.3 2014/10/13 08:53:39 j_novak Exp $" ;
41#include "star_rot_dirac_diff.h"
43#include "utilitaires.h"
92 cout <<
"Star_rot_Dirac_diff::equilibrium: mer_change_omega < mer_rot !"
99 cout <<
"Star_rot_Dirac_diff::equilibrium: mer_fix_omega < mer_change_omega !"
114 diff.set_etat_qcq() ;
147 fichconv <<
"# diff_ent GRV2 max_triax vit_triax" <<
'\n' ;
154 "# |dH/dr_eq/dH/dr_pole| r_pole/r_eq ent_c"
168 cout <<
"-----------------------------------------------" <<
'\n' ;
169 cout <<
"step: " <<
mer <<
'\n' ;
290 cout <<
"U > c for l, i : " <<
l <<
" " <<
i
291 <<
" U = " <<
sqrt(
u2 ) <<
'\n' ;
296 cout <<
"**** VELOCITY OF LIGHT REACHED ****" <<
'\n' ;
298 cout <<
"New rotation frequency : "
299 <<
omega/(2.*
M_PI) * f_unit <<
" Hz" <<
'\n' ;
371 cout <<
"ent < 0 for l, i : " <<
l <<
" " <<
i
378 cout <<
"**** KEPLERIAN VELOCITY REACHED ****" <<
'\n' ;
380 cout <<
"New central rotation frequency : "
411 q_new.std_spectral_base() ;
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.
void fait_omega_field(double omeg_min, double omeg_max, double precis, int nitermax)
Computes (member omega_field ).
virtual 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.
Scalar omega_field
Field .
virtual void hydro_euler()
Computes the hydrodynamical quantities relative to the Eulerian observer from those in the fluid fram...
Tbl par_frot
Parameters of the function .
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.
void solve_logn_q(Scalar &ln_q_new) const
Solution of the two scalar Poisson equations for rotating stars in Dirac gauge.
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.
double ray_eq() const
Coordinate radius at , [r_unit].
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.
double & set(int i)
Read/write of a particular element (index i) (1D case)
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 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.