30char binaire_global_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Binaire/binaire_global.C,v 1.7 2014/10/13 08:52:44 j_novak Exp $" ;
101 for (
int i=0; i<=1; i++) {
108 Cmp source =
pow(a2, 1.25) * ee
109 +
pow(a2, 0.25) * (ak2_auto + ak2_comp) / (4.*qpig) ;
150 for (
int i=0; i<=1; i++) {
165 dbe_auto + dbe_comp) ;
167 dnu_auto + dnu_comp) ;
169 Cmp source = lapse * ( a2 * (ee + se)
170 + (ak2_auto + ak2_comp)/qpig
171 - dndb()/qpig + dndn()/qpig ) ;
206 for (
int i=0; i<=1; i++) {
224 if (
et[i]->is_relativistic() ) {
228 rho =
pow(a2, 2.5) * (ee + pp) ;
240 Cmp source = rho * ( yy * vz - zz * vy ) ;
242 (source.
va).set_base( *(base[2]) ) ;
251 source = rho * ( zz * vx - xx * vz ) ;
253 (source.
va).set_base( *(base[2]) ) ;
262 source = rho * ( xx * vy - yy * vx ) ;
308 for (
int i=0; i<=1; i++) {
324 Cmp source = e_int + .5 * rho * vit2 + .5 * rho * nu ;
366 double vir_grav = 0 ;
368 for (
int i=0; i<=1; i++) {
383 Cmp source = 3*pp + rho * vit2 ;
387 source = .5 * rho * nu ;
393 *
p_virial = ( vir_mat + vir_grav ) / fabs(vir_grav) ;
423 double vir_pres = 0. ;
424 double vir_extr = 0. ;
425 double vir_grav = 0. ;
427 for (
int i=0; i<=1; i++) {
442 Cmp source = 2. * a2 * a1 * se ;
445 source = 1.5 * a1 * (ak2_auto + ak2_comp) / qpig ;
453 source = a1 * ( sprod1() - sprod2() - 2.*sprod3() )/qpig ;
496 double vir_pres = 0. ;
497 double vir_extr = 0. ;
498 double vir_grav = 0. ;
500 for (
int i=0; i<=1; i++) {
516 Cmp source = 2. * lapse * a2 * a1 * se ;
519 source = 1.5 * lapse * a1 * (ak2_auto + ak2_comp) / qpig ;
525 source = lapse * a1 * sprod() / qpig ;
Bases of the spectral expansions.
Tbl * p_angu_mom
Total angular momentum of the system.
double * p_virial_gb
Virial theorem error by E.Gourgoulhon and S.Bonazzola.
double mass_adm() const
Total ADM mass.
double virial_gb() const
Estimates the relative error on the virial theorem calculated by E.Gourgoulhon and S....
Etoile_bin * et[2]
Array of the two stars (to perform loops on the stars): {\tt et[0]} contains the address of {\tt star...
const Tbl & angu_mom() const
Total angular momentum.
double * p_mass_adm
Total ADM mass of the system.
double * p_total_ener
Total energy of the system.
double * p_virial
Virial theorem error.
Etoile_bin star2
Second star of the system.
double mass_kom() const
Total Komar mass.
double * p_mass_kom
Total Komar mass of the system.
double virial() const
Estimates the relative error on the virial theorem (for a relativistic one, it returns $|1 - M_{\rm K...
double * p_virial_fus
Virial theorem error by J.L.Friedman, K.Uryu, and M.Shibata.
Etoile_bin star1
First star of the system.
double total_ener() const
Total energy (excluding the rest mass energy).
double virial_fus() const
Estimates the relative error on the virial theorem calculated by J.L.Friedman, K.Uryu,...
Component of a tensorial field *** DEPRECATED : use class Scalar instead ***.
Valeur va
The numerical value of the Cmp
void std_base_scal()
Sets the spectral bases of the Valeur va to the standard ones for a scalar.
double integrale() const
Computes the integral over all space of *this .
virtual double mass_b() const
Baryon mass.
const Tenseur & get_d_logn_auto() const
Returns the gradient of logn_auto (Cartesian components with respect to ref_triad )
const Tenseur & get_logn_comp() const
Returns the part of the lapse logarithm (gravitational potential at the Newtonian limit) generated pr...
const Tenseur & get_d_beta_auto() const
Returns the gradient of beta_auto (Cartesian components with respect to ref_triad )
const Tenseur & get_akcar_comp() const
Returns the part of the scalar generated by shift_auto and shift_comp , i.e.
const Tenseur & get_akcar_auto() const
Returns the part of the scalar generated by shift_auto , i.e.
const Tenseur & get_d_logn_comp() const
Returns the gradient of logn_comp (Cartesian components with respect to ref_triad )
const Tenseur & get_d_beta_comp() const
Returns the gradient of beta_comp (Cartesian components with respect to ref_triad )
const Tenseur & get_nnn() const
Returns the total lapse function N.
const Map & get_mp() const
Returns the mapping.
bool is_relativistic() const
Returns true for a relativistic star, false for a Newtonian one.
const Tenseur & get_a_car() const
Returns the total conformal factor .
const Tenseur & get_logn_auto() const
Returns the logarithm of the part of the lapse N generated principaly by the star.
const Tenseur & get_u_euler() const
Returns the fluid 3-velocity with respect to the Eulerian observer.
const Tenseur & get_s_euler() const
Returns the trace of the stress tensor in the Eulerian frame.
const Tenseur & get_ener() const
Returns the proper total energy density.
const Tenseur & get_press() const
Returns the fluid pressure.
const Tenseur & get_ener_euler() const
Returns the total energy density with respect to the Eulerian observer.
const Tenseur & get_nbar() const
Returns the proper baryon density.
Base class for coordinate mappings.
Coord ya
Absolute y coordinate.
Coord za
Absolute z coordinate.
Coord xa
Absolute x coordinate.
void annule_hard()
Sets the Tbl to zero in a hard way.
double & set(int i)
Read/write of a particular element (index i) (1D case)
Tensor handling *** DEPRECATED : use class Tensor instead ***.
Cmp sqrt(const Cmp &)
Square root.
Cmp pow(const Cmp &, int)
Power .
Tenseur flat_scalar_prod(const Tenseur &t1, const Tenseur &t2)
Scalar product of two Tenseur when the metric is : performs the contraction of the last index of t1 w...
Standard units of space, time and mass.