21char pseudo_misner_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Bin_ns_bh/pseudo_misner.C,v 1.4 2014/10/13 08:52:43 j_novak Exp $" ;
49#include "et_bin_nsbh.h"
55#include "utilitaires.h"
59void Bin_ns_bh::pseudo_misner (
int& ite,
int itemax,
double relax,
60 double precis,
int bound_nn,
double lim_nn) {
67 source_n_prev.set_etat_zero() ;
70 par_poisson1.add_int(itemax, 0) ;
71 par_poisson1.add_double(relax, 0) ;
72 par_poisson1.add_double(precis, 1) ;
73 par_poisson1.add_int_mod(ite, 0) ;
74 par_poisson1.add_cmp_mod(source_n_prev) ;
79 source_psi_prev.set_etat_zero() ;
82 par_poisson2.add_int(itemax, 0) ;
83 par_poisson2.add_double(relax, 0) ;
84 par_poisson2.add_double(precis, 1) ;
85 par_poisson2.add_int_mod(ite, 0) ;
86 par_poisson2.add_cmp_mod(source_psi_prev) ;
109 source_n.set_std_base() ;
117 source_psi.set_std_base() ;
137 cout <<
"Step " << itere <<
" " << erreur << endl ;
139 if ((itere==itemax) || (erreur<precis))
Tenseur n_auto
Part of N generated by the hole.
void solve_psi_with_ns(double relax)
Solves the equation for ~:
void solve_lapse_with_ns(double relax, int bound_nn, double lim_nn)
Solves the equation for N ~:
double get_omega() const
Returns the orbital velocity.
Et_bin_nsbh star
The neutron star.
Bhole hole
The black hole.
double get_x_axe() const
Returns a constant reference to the black hole.
Tenseur confpsi_auto
Part of the conformal factor $\Psi$ generated principaly by the star.
Tenseur d_confpsi_comp
Gradient of {\tt confpsi_comp} (Cartesian components with respect to {\tt ref_triad})
void update_metric(const Bhole &comp)
Computes metric coefficients from known potentials, when the companion is a black hole.
Tenseur d_n_auto
Gradient of {\tt n_auto} (Cartesian components with respect to {\tt ref_triad})
virtual void kinematics(double omega, double x_axe)
Computes the quantities bsn and pot_centri .
Tenseur n_auto
Part of the lapse {\it N} generated principaly by the star.
void update_metric_der_comp(const Bhole &comp)
Computes the derivative of metric functions related to the companion black hole.
Tenseur confpsi
Total conformal factor $\Psi$.
Tenseur d_confpsi_auto
Gradient of {\tt confpsi_auto} (Cartesian components with respect to {\tt ref_triad})
virtual void hydro_euler()
Computes the hydrodynamical quantities relative to the Eulerian observer from those in the fluid fram...
void fait_d_psi()
Computes the gradient of the total velocity potential .
Tenseur nnn
Total lapse function.
const Map & get_mp() const
Returns the mapping.
virtual void equation_of_state()
Computes the proper baryon and energy density, as well as pressure from the enthalpy.
Tenseur ener_euler
Total energy density in the Eulerian frame.
Tenseur s_euler
Trace of the stress tensor in the Eulerian frame.
Cmp & set()
Read/write for a scalar (see also operator=(const Cmp&) ).
Tbl max(const Cmp &)
Maximum values of a Cmp in each domain.
Cmp pow(const Cmp &, int)
Power .
Tbl diffrelmax(const Cmp &a, const Cmp &b)
Relative difference between two Cmp (max version).
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.