28char star_bhns_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Star_bhns/star_bhns.C,v 1.4 2014/10/13 08:53:40 j_novak Exp $" ;
65Cmp raccord_c1(
const Cmp& uu,
int l1) ;
78 d_psi(
mp_i, COV,
mp_i.get_bvect_cart()),
79 wit_w(
mp_i, CON,
mp_i.get_bvect_cart()),
81 bsn(
mp_i, CON,
mp_i.get_bvect_cart()),
90 d_lapconf_auto(
mp_i, COV,
mp_i.get_bvect_cart()),
91 d_lapconf_comp(
mp_i, COV,
mp_i.get_bvect_cart()),
92 shift_auto(
mp_i, CON,
mp_i.get_bvect_cart()),
93 shift_comp(
mp_i, CON,
mp_i.get_bvect_cart()),
94 shift_tot(
mp_i, CON,
mp_i.get_bvect_cart()),
95 d_shift_auto(
mp_i, 2, CON,
mp_i.get_bvect_cart()),
96 d_shift_comp(
mp_i, 2, CON,
mp_i.get_bvect_cart()),
100 d_confo_auto(
mp_i, COV,
mp_i.get_bvect_cart()),
101 d_confo_comp(
mp_i, COV,
mp_i.get_bvect_cart()),
103 taij_auto(
mp_i, CON,
mp_i.get_bvect_cart()),
104 taij_quad_auto(
mp_i),
109 ssjm1_wshift(
mp_i, CON,
mp_i.get_bvect_cart())
181 irrotational(star.irrotational),
189 pot_centri(star.pot_centri),
190 lapconf_auto(star.lapconf_auto),
191 lapconf_comp(star.lapconf_comp),
192 lapconf_tot(star.lapconf_tot),
193 lapse_auto(star.lapse_auto),
194 lapse_tot(star.lapse_tot),
195 d_lapconf_auto(star.d_lapconf_auto),
196 d_lapconf_comp(star.d_lapconf_comp),
197 shift_auto(star.shift_auto),
198 shift_comp(star.shift_comp),
199 shift_tot(star.shift_tot),
200 d_shift_auto(star.d_shift_auto),
201 d_shift_comp(star.d_shift_comp),
202 confo_auto(star.confo_auto),
203 confo_comp(star.confo_comp),
204 confo_tot(star.confo_tot),
205 d_confo_auto(star.d_confo_auto),
206 d_confo_comp(star.d_confo_comp),
208 taij_auto(star.taij_auto),
209 taij_quad_auto(star.taij_quad_auto),
211 ssjm1_lapconf(star.ssjm1_lapconf),
212 ssjm1_confo(star.ssjm1_confo),
213 ssjm1_khi(star.ssjm1_khi),
214 ssjm1_wshift(star.ssjm1_wshift)
225 d_psi(
mp_i, COV,
mp_i.get_bvect_cart()),
226 wit_w(
mp_i, CON,
mp_i.get_bvect_cart()),
228 bsn(
mp_i, CON,
mp_i.get_bvect_cart()),
237 d_lapconf_auto(
mp_i, COV,
mp_i.get_bvect_cart()),
238 d_lapconf_comp(
mp_i, COV,
mp_i.get_bvect_cart()),
240 shift_comp(
mp_i, CON,
mp_i.get_bvect_cart()),
241 shift_tot(
mp_i, CON,
mp_i.get_bvect_cart()),
242 d_shift_auto(
mp_i, 2, CON,
mp_i.get_bvect_cart()),
243 d_shift_comp(
mp_i, 2, CON,
mp_i.get_bvect_cart()),
247 d_confo_auto(
mp_i, COV,
mp_i.get_bvect_cart()),
248 d_confo_comp(
mp_i, COV,
mp_i.get_bvect_cart()),
250 taij_auto(
mp_i, CON,
mp_i.get_bvect_cart()),
251 taij_quad_auto(
mp_i),
502 ost <<
"Neutron star in a BHNS binary" <<
endl ;
503 ost <<
"-----------------------------" <<
endl ;
505 ost <<
"Coordinate radius R_eq_tow : "
507 ost <<
"Coordinate radius R_eq_opp : "
509 ost <<
"Coordinate radius R_eq_orb : "
511 ost <<
"Coordinate radius R_pole : "
513 ost <<
"Central enthalph H_ent : "
515 ost <<
"Lapse function at the center of NS : "
517 ost <<
"Conformal factor at the center of NS : "
519 ost <<
"shift(1) at the center of NS : "
521 ost <<
"shift(2) at the center of NS : "
523 ost <<
"shift(3) at the center of NS : "
557 for (
int i=1;
i<=3;
i++) {
558 v_orb.set(
i) = www(
i).val_grid_point(0,0,0,0) ;
566 for (
int i=1;
i<=3;
i++)
569 d_psi0.std_spectral_base() ;
573 for (
int i=1;
i<=3;
i++) {
574 if (
d_psi(
i).get_etat() == ETATZERO)
584 for (
int i=1;
i<=3;
i++) {
Component of a tensorial field *** DEPRECATED : use class Scalar instead ***.
Equation of state base class.
Time evolution with partial storage (*** under development ***).
Base class for coordinate mappings.
const Base_vect_cart & get_bvect_cart() const
Returns the Cartesian basis associated with the coordinates (x,y,z) of the mapping,...
const Mg3d * get_mg() const
Gives the Mg3d on which the mapping is defined.
Tensor field of valence 0 (or component of a tensorial field).
virtual void sauve(FILE *) const
Save in a file.
const Scalar & deriv(int i) const
Returns of *this , where .
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.
void annule_hard()
Sets the Scalar to zero in a hard way.
Class for stars in black hole-neutron star binaries.
Metric_flat flat
Flat metric defined on the mapping (Spherical components with respect to the mapping of the star ).
Scalar & set_confo_comp()
Read/write of the conformal factor generated by the companion black hole.
virtual void sauve(FILE *) const
Save in a file.
Vector shift_tot
Total shift vector.
Sym_tensor taij_auto
Part of the extrinsic curvature tensor generated by shift_auto , lapse_auto , and confo_auto .
Scalar lapconf_auto
Lapconf function generated by the star.
Scalar confo_comp
Conformal factor generated by the companion black hole.
Vector wit_w
Spatial projection of the fluid 3-velocity with respect to the co-orbiting observer.
Tensor d_shift_auto
Derivative of the shift vector generated by the star .
Vector bsn
3-vector shift, divided by N , of the rotating coordinates, .
Scalar & set_confo_auto()
Read/write of the conformal factor generated by the neutron star.
Vector & set_shift_auto()
Read/write of the shift vector generated by the neutron star.
Scalar & set_lapconf_auto()
Read/write of the lapconf function generated by the neutron star.
Scalar & set_pot_centri()
Read/write the centrifugal potential.
void fait_d_psi_bhns()
Computes the gradient of the total velocity potential .
Scalar lapse_tot
Total lapse function.
Vector ssjm1_wshift
Effective source at the previous step for the resolution of the vector Poisson equation for by means...
virtual void del_deriv() const
Deletes all the derived quantities.
Vector d_confo_comp
Derivative of the conformal factor generated by the companion black hole.
void operator=(const Star_bhns &)
Assignment to another Star_bhns.
Scalar confo_tot
Total conformal factor.
Scalar psi4
Fourth power of the total conformal factor.
Vector d_lapconf_auto
Derivative of the lapconf function generated by the star .
Vector d_psi
Gradient of (in the irrotational case) (Spherical components with respect to the mapping of the star...
Vector d_confo_auto
Derivative of the conformal factor generated by the star .
Scalar & set_lapconf_comp()
Read/write of the lapconf function generated by the companion black hole.
Scalar gam0
Lorentz factor between the co-orbiting observer and the Eulerian one.
Vector shift_comp
Shift vector generated by the companion black hole.
Scalar lapconf_tot
Total lapconf function.
Star_bhns(Map &mp_i, int nzet_i, const Eos &eos_i, bool irrot_i)
Standard constructor.
Vector d_lapconf_comp
Derivative of the lapconf function generated by the companion black hole.
Vector & set_shift_comp()
Read/write of the shift vector generated by the companion black hole.
Scalar gam
Lorentz factor between the fluid and the co-orbiting observer.
Scalar ssjm1_lapconf
Effective source at the previous step for the resolution of the Poisson equation for lapconf_auto .
bool irrotational
true for an irrotational star, false for a corotating one
Scalar pot_centri
Centrifugal potential.
void set_der_0x0() const
Sets to 0x0 all the pointers on derived quantities.
double * p_mass_b_bhns
Baryon mass.
Scalar lapse_auto
Lapse function generated by the "star".
void relax_bhns(const Star_bhns &star_prev, double relax_ent, double relax_met, int mer, int fmer_met)
Performs a relaxation on ent , lapse_auto , shift_auto , confo_auto .
Scalar psi0
Scalar potential of the non-translational part of fluid 4-velocity (in the irrotational case)
Scalar lapconf_comp
Lapconf function generated by the companion black hole.
Scalar taij_quad_auto
Part of the scalar generated by .
virtual ostream & operator>>(ostream &) const
Operator >> (virtual function called by the operator <<).
Vector shift_auto
Shift vector generated by the star.
Scalar ssjm1_khi
Effective source at the previous step for the resolution of the Poisson equation for the scalar by m...
Scalar confo_auto
Conformal factor generated by the star.
Scalar ssjm1_confo
Effective source at the previous step for the resolution of the Poisson equation for confo_auto .
Map_af mp_aff
Affine mapping for solving poisson's equations of metric quantities.
Tensor d_shift_comp
Derivative of the shift vector generated by the companion black hole.
Scalar loggam
Logarithm of the Lorentz factor between the fluid and the co-orbiting observer.
virtual ~Star_bhns()
Destructor.
double * p_mass_g_bhns
Gravitational mass.
void equation_of_state()
Computes the proper baryon and energy density, as well as pressure from the enthalpy.
virtual void del_deriv() const
Deletes all the derived quantities.
Scalar gam_euler
Lorentz factor between the fluid and Eulerian observers.
virtual void set_der_0x0() const
Sets to 0x0 all the pointers on derived quantities.
double ray_eq() const
Coordinate radius at , [r_unit].
virtual void sauve(FILE *) const
Save in a file.
double ray_eq_pis2() const
Coordinate radius at , [r_unit].
double ray_eq_pi() 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.
void operator=(const Star &)
Assignment to another Star.
int nzet
Number of domains of *mp occupied by the star.
double ray_pole() const
Coordinate radius at [r_unit].
Tensor field of valence 1.
virtual void change_triad(const Base_vect &)
Sets a new vectorial basis (triad) of decomposition and modifies the components accordingly.
Scalar & set(int)
Read/write access to a component.
Cmp exp(const Cmp &)
Exponential.
virtual void sauve(FILE *) const
Save in a binary file.
virtual void set_etat_nondef()
Sets the logical state of all components to ETATNONDEF (undefined state).
virtual void annule(int l_min, int l_max)
Sets the Tensor to zero in several domains.
virtual void set_etat_zero()
Sets the logical state of all components to ETATZERO (zero state).
Standard units of space, time and mass.