53 #include "et_bin_bhns_extr.h" 87 r_bh.
set() =
pow( (xx+sepa)*(xx+sepa) + yy*yy + zz*zz, 0.5) ;
92 xx_cov.
set(0) = xx + sepa ;
98 xsr_cov = xx_cov / r_bh ;
102 msr = ggrav * mass / r_bh ;
110 for (
int i=0; i<3; i++)
111 tmp1.
set() += xsr_cov(i) %
bsn(i) ;
115 Tenseur tmp2 = 2.*msr % tmp1 % tmp1 ;
118 for (
int i=0; i<3; i++)
121 tmp2 =
a_car % tmp2 ;
137 xx_con.
set(0) = xx + sepa ;
143 xsr_con = xx_con / r_bh ;
151 for (
int i=0; i<3; i++)
152 tmp3.
set() += xsr_con(i) %
d_psi(i) ;
156 Tenseur tmp4 = -2.*msr % tmp3 % tmp3 / (1.+2.*msr) ;
159 for (
int i=0; i<3; i++)
162 tmp4 = tmp4 /
a_car ;
171 Tenseur vtmp2 = -2.* msr % tmp3 % xsr_con / (1.+2.*msr)
181 assert(vtmp1.
get_etat() == ETATQCQ) ;
183 for (
int i=0; i<3; i++) {
218 for (
int i=0; i<3; i++)
223 Tenseur stmp2 = 2.*msr % stmp1 % stmp1 ;
226 for (
int i=0; i<3; i++)
229 stmp2 =
a_car % stmp2 ;
239 Tenseur gtmp = 2.*msr % tmp1 % stmp1 ;
242 for (
int i=0; i<3; i++)
245 gtmp =
a_car % gtmp ;
274 cout <<
"BH-NS binary systems should be relativistic !!!" << endl ;
350 assert(vtmp.
get_etat() == ETATQCQ) ;
352 for (
int i=0; i<3; i++) {
419 cout <<
"BH-NS binary systems should be relativistic !!!" << endl ;
Cmp log(const Cmp &)
Neperian logarithm.
void annule(int l)
Sets the Tenseur to zero in a given domain.
const Base_vect & ref_triad
Reference triad ("absolute frame"), with respect to which the components of all the member Tenseur 's...
Cmp exp(const Cmp &)
Exponential.
void set_triad(const Base_vect &new_triad)
Assigns a new vectorial basis (triad) of decomposition.
Cmp sqrt(const Cmp &)
Square root.
void set_std_base()
Set the standard spectal basis of decomposition for each component.
Standard units of space, time and mass.
const Mg3d * get_mg() const
Gives the Mg3d on which the mapping is defined.
Tenseur s_euler
Trace of the stress tensor in the Eulerian frame.
Tenseur wit_w
Spatial projection of the fluid 3-velocity with respect to the co-orbiting observer.
double unsurc2
: unsurc2=1 for a relativistic star, 0 for a Newtonian one.
Tenseur press
Fluid pressure.
bool kerrschild
Indicator of the background metric: true for the Kerr-Shild metric, false for the conformally flat on...
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...
bool irrotational
true for an irrotational star, false for a corotating one
Cmp & set()
Read/write for a scalar (see also operator=(const Cmp&) ).
Tenseur u_euler
Fluid 3-velocity with respect to the Eulerian observer.
virtual void del_deriv() const
Deletes all the derived quantities.
Tenseur bsn
3-vector shift, divided by N , of the rotating coordinates, .
const Base_vect * get_triad() const
Returns the vectorial basis (triad) on which the components are defined.
Tenseur gam_euler
Lorentz factor between the fluid and Eulerian observers.
Map & mp
Mapping associated with the star.
int get_nzone() const
Returns the number of domains.
int get_etat() const
Returns the logical state.
Cmp pow(const Cmp &, int)
Power .
Active physical coordinates and mapping derivatives.
Tenseur a_car
Total conformal factor .
bool relativistic
Indicator of relativity: true for a relativistic star, false for a Newtonian one. ...
Tenseur ener
Total energy density in the fluid frame.
Tenseur loggam
Logarithm of the Lorentz factor between the fluid and the co-orbiting observer.
Tenseur ent
Log-enthalpy (relativistic case) or specific enthalpy (Newtonian case)
Coord y
y coordinate centered on the grid
Coord x
x coordinate centered on the grid
void set_etat_qcq()
Sets the logical state to ETATQCQ (ordinary state).
Tenseur ener_euler
Total energy density in the Eulerian frame.
void set_etat_zero()
Sets the logical state to ETATZERO (zero state).
Coord z
z coordinate centered on the grid
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...
Tensor handling *** DEPRECATED : use class Tensor instead ***.
Tenseur d_psi
Gradient of (in the irrotational case) (Cartesian components with respect to ref_triad ) ...