30 char et_bin_bhns_extr_hydro_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Etoile/et_bin_bhns_extr_hydro.C,v 1.3 2014/10/13 08:52:55 j_novak Exp $" ;
50 #include "et_bin_bhns_extr.h"
84 r_bh.
set() =
pow( (xx+sepa)*(xx+sepa) + yy*yy + zz*zz, 0.5) ;
89 xx_cov.
set(0) = xx + sepa ;
95 xsr_cov = xx_cov / r_bh ;
99 msr = ggrav * mass / r_bh ;
107 for (
int i=0; i<3; i++)
108 tmp1.
set() += xsr_cov(i) %
bsn(i) ;
112 Tenseur tmp2 = 2.*msr % tmp1 % tmp1 ;
115 for (
int i=0; i<3; i++)
118 tmp2 =
a_car % tmp2 ;
134 xx_con.
set(0) = xx + sepa ;
140 xsr_con = xx_con / r_bh ;
148 for (
int i=0; i<3; i++)
149 tmp3.
set() += xsr_con(i) %
d_psi(i) ;
153 Tenseur tmp4 = -2.*msr % tmp3 % tmp3 / (1.+2.*msr) ;
156 for (
int i=0; i<3; i++)
159 tmp4 = tmp4 /
a_car ;
168 Tenseur vtmp2 = -2.* msr % tmp3 % xsr_con / (1.+2.*msr)
178 assert(vtmp1.
get_etat() == ETATQCQ) ;
180 for (
int i=0; i<3; i++) {
215 for (
int i=0; i<3; i++)
220 Tenseur stmp2 = 2.*msr % stmp1 % stmp1 ;
223 for (
int i=0; i<3; i++)
226 stmp2 =
a_car % stmp2 ;
236 Tenseur gtmp = 2.*msr % tmp1 % stmp1 ;
239 for (
int i=0; i<3; i++)
242 gtmp =
a_car % gtmp ;
271 cout <<
"BH-NS binary systems should be relativistic !!!" << endl ;
347 assert(vtmp.
get_etat() == ETATQCQ) ;
349 for (
int i=0; i<3; i++) {
416 cout <<
"BH-NS binary systems should be relativistic !!!" << endl ;
Tenseur press
Fluid pressure.
Coord z
z coordinate centered on the grid
Tenseur a_car
Total conformal factor .
bool relativistic
Indicator of relativity: true for a relativistic star, false for a Newtonian one.
Tenseur u_euler
Fluid 3-velocity with respect to the Eulerian observer.
void set_etat_zero()
Sets the logical state to ETATZERO (zero state).
Tensor handling *** DEPRECATED : use class Tensor instead ***.
Tenseur ent
Log-enthalpy (relativistic case) or specific enthalpy (Newtonian case)
Cmp & set()
Read/write for a scalar (see also operator=(const Cmp&) ).
Tenseur bsn
3-vector shift, divided by N , of the rotating coordinates, .
Cmp log(const Cmp &)
Neperian logarithm.
Cmp exp(const Cmp &)
Exponential.
const Mg3d * get_mg() const
Gives the Mg3d on which the mapping is defined.
Tenseur d_psi
Gradient of (in the irrotational case) (Cartesian components with respect to ref_triad )
Map & mp
Mapping associated with the star.
const Base_vect * get_triad() const
Returns the vectorial basis (triad) on which the components are defined.
void set_triad(const Base_vect &new_triad)
Assigns a new vectorial basis (triad) of decomposition.
Cmp pow(const Cmp &, int)
Power .
const Base_vect & ref_triad
Reference triad ("absolute frame"), with respect to which the components of all the member Tenseur 's...
int get_etat() const
Returns the logical state.
Standard units of space, time and mass.
Tenseur ener
Total energy density in the fluid frame.
Tenseur loggam
Logarithm of the Lorentz factor between the fluid and the co-orbiting observer.
virtual Tenseur sprod(const Tenseur &t1, const Tenseur &t2) const
Performs the scalar product of two tensors by contracting the last index of t1 with the first index o...
Active physical coordinates and mapping derivatives.
bool kerrschild
Indicator of the background metric: true for the Kerr-Shild metric, false for the conformally flat on...
int get_nzone() const
Returns the number of domains.
Tenseur ener_euler
Total energy density in the Eulerian frame.
Cmp sqrt(const Cmp &)
Square root.
void hydro_euler_extr(const double &mass, const double &sepa)
Computes the hydrodynamical quantities relative to the Eulerian observer from those in the fluid fram...
Coord x
x coordinate centered on the grid
void set_etat_qcq()
Sets the logical state to ETATQCQ (ordinary state).
Tenseur wit_w
Spatial projection of the fluid 3-velocity with respect to the co-orbiting observer.
void annule(int l)
Sets the Tenseur to zero in a given domain.
Tenseur gam_euler
Lorentz factor between the fluid and Eulerian observers.
bool irrotational
true for an irrotational star, false for a corotating one
Coord y
y coordinate centered on the grid
virtual void del_deriv() const
Deletes all the derived quantities.
Tenseur s_euler
Trace of the stress tensor in the Eulerian frame.
double unsurc2
: unsurc2=1 for a relativistic star, 0 for a Newtonian one.
void set_std_base()
Set the standard spectal basis of decomposition for each component.