52 #include "et_bin_nsbh.h" 55 #include "bin_ns_bh.h" 57 #include "graphique.h" 58 #include "utilitaires.h" 62 void Bin_ns_bh::pseudo_misner (
int& ite,
int itemax,
double relax,
63 double precis,
int bound_nn,
double lim_nn) {
70 source_n_prev.set_etat_zero() ;
73 par_poisson1.add_int(itemax, 0) ;
74 par_poisson1.add_double(relax, 0) ;
75 par_poisson1.add_double(precis, 1) ;
76 par_poisson1.add_int_mod(ite, 0) ;
77 par_poisson1.add_cmp_mod(source_n_prev) ;
82 source_psi_prev.set_etat_zero() ;
85 par_poisson2.add_int(itemax, 0) ;
86 par_poisson2.add_double(relax, 0) ;
87 par_poisson2.add_double(precis, 1) ;
88 par_poisson2.add_int_mod(ite, 0) ;
89 par_poisson2.add_cmp_mod(source_psi_prev) ;
112 source_n.set_std_base() ;
120 source_psi.set_std_base() ;
140 cout <<
"Step " << itere <<
" " << erreur << endl ;
142 if ((itere==itemax) || (erreur<precis))
const Map & get_mp() const
Returns the mapping.
Standard units of space, time and mass.
virtual void kinematics(double omega, double x_axe)
Computes the quantities bsn and pot_centri .
double get_omega() const
Returns the orbital velocity.
Tenseur nnn
Total lapse function.
Tenseur s_euler
Trace of the stress tensor in the Eulerian frame.
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...
void update_metric_der_comp(const Bhole &comp)
Computes the derivative of metric functions related to the companion black hole.
Tenseur n_auto
Part of the lapse { N} generated principaly by the star.
Tenseur confpsi
Total conformal factor $$.
Tenseur d_confpsi_comp
Gradient of { confpsi_comp} (Cartesian components with respect to { ref_triad})
void solve_psi_with_ns(double relax)
Solves the equation for ~: with the condition that on the horizon, where f is the value of on th...
Cmp & set()
Read/write for a scalar (see also operator=(const Cmp&) ).
Et_bin_nsbh star
The neutron star.
void fait_d_psi()
Computes the gradient of the total velocity potential .
virtual void equation_of_state()
Computes the proper baryon and energy density, as well as pressure from the enthalpy.
Tbl max(const Cmp &)
Maximum values of a Cmp in each domain.
Cmp pow(const Cmp &, int)
Power .
Tenseur d_confpsi_auto
Gradient of { confpsi_auto} (Cartesian components with respect to { ref_triad})
Tenseur d_n_auto
Gradient of { n_auto} (Cartesian components with respect to { ref_triad})
void update_metric(const Bhole &comp)
Computes metric coefficients from known potentials, when the companion is a black hole...
Tenseur confpsi_auto
Part of the conformal factor $$ generated principaly by the star.
double get_x_axe() const
Returns a constant reference to the black hole.
void solve_lapse_with_ns(double relax, int bound_nn, double lim_nn)
Solves the equation for N ~: with the condition that N =0 on the horizon.
Bhole hole
The black hole.
Tenseur n_auto
Part of N generated by the hole.
Tenseur ener_euler
Total energy density in the Eulerian frame.
Tbl diffrelmax(const Cmp &a, const Cmp &b)
Relative difference between two Cmp (max version).
virtual void hydro_euler()
Computes the hydrodynamical quantities relative to the Eulerian observer from those in the fluid fram...