58 #include "hole_bhns.h" 59 #include "star_bhns.h" 60 #include "utilitaires.h" 134 double relax_jm1 = 1. - relax ;
148 rr.std_spectral_base() ;
164 ll.set(1) = st % cp ;
165 ll.
set(2) = st % sp ;
193 for (
int i=1; i<=3; i++)
199 for (
int i=1; i<=3; i++) {
211 for (
int i=1; i<=3; i++) {
212 for (
int j=1; j<=3; j++) {
220 for (
int i=1; i<=3; i++) {
221 for (
int j=1; j<=3; j++) {
238 for (
int i=1; i<=3; i++)
244 for (
int i=1; i<=3; i++)
263 cc = 2. * (
sqrt(13.) - 1.) / 3. ;
275 cout <<
"!!!!! WARNING: Not yet prepared !!!!!" << endl ;
281 cout <<
"!!!!! WARNING: Not yet prepared !!!!!" << endl ;
296 + cc*cc*
pow(mass/r_are/rr, 4.)) *
sqrt(r_are) ;
309 for (
int i=1; i<=3; i++)
317 for (
int i=1; i<=3; i++)
323 for (
int i=1; i<=3; i++) {
325 * (mass/r_are/rr - 2.*cc*cc*
pow(mass/r_are/rr,4.))
328 * (
sqrt(1. - 2.*mass/r_are/rr + cc*cc*
pow(mass/r_are/rr,4.))-1.)
329 *
sqrt(1. - 2.*mass/r_are/rr + cc*cc*
pow(mass/r_are/rr,4.))
339 for (
int i=1; i<=3; i++)
343 for (
int i=1; i<=3; i++) {
344 for (
int j=1; j<=3; j++) {
352 for (
int i=1; i<=3; i++) {
353 for (
int j=1; j<=3; j++) {
356 -3.*
sqrt(1. - 2.*mass/r_are/rr
357 +cc*cc*
pow(mass/r_are/rr,4.))
369 for (
int i=1; i<=3; i++) {
370 for (
int j=1; j<=3; j++) {
376 for (
int i=1; i<=3; i++)
382 for (
int i=1; i<=3; i++) {
384 *(
sqrt(1. - 2.*mass/r_are/rr +cc*cc*
pow(mass/r_are/rr,4.)) - 1.)
394 for (
int i=1; i<=3; i++)
416 for (
int i=1; i<=3; i++) {
423 for (
int i=1; i<=3; i++) {
Map & mp
Mapping associated with the black hole.
void annule_domain(int l)
Sets the Tensor to zero in a given domain.
Scalar lapconf_auto_bh
Part of the lapconf function from the analytic background.
Class for stars in black hole-neutron star binaries.
virtual void set_etat_qcq()
Sets the logical state of all components to ETATQCQ (ordinary state).
Vector d_confo_auto
Derivative of the conformal factor generated by the black hole.
Cmp sqrt(const Cmp &)
Square root.
virtual void set_etat_zero()
Sets the logical state to ETATZERO (zero).
Scalar confo_auto
Conformal factor generated by the black hole.
double mass_bh
Gravitational mass of BH.
Standard units of space, time and mass.
Tensor field of valence 0 (or component of a tensorial field).
Vector d_lapconf_auto_bh
Derivative of the part of the lapconf function from the analytic background.
bool kerrschild
true for a Kerr-Schild background, false for a conformally flat background
Base class for coordinate mappings.
Vector shift_auto
Shift vector generated by the black hole.
Vector d_lapconf_auto
Derivative of the lapconf function generated by the black hole.
Tensor d_shift_auto
Derivative of the shift vector generated by the black hole.
Vector shift_tot
Total shift vector ;.
const Scalar & get_lapconf_auto() const
Returns the part of the lapconf function generated by the star.
virtual void std_spectral_base()
Sets the spectral bases of the Valeur va to the standard ones for a scalar field. ...
const Map & get_mp() const
Returns the mapping.
virtual void change_triad(const Base_vect &)
Sets a new vectorial basis (triad) of decomposition and modifies the components accordingly.
virtual void set_etat_qcq()
Sets the logical state to ETATQCQ (ordinary state).
Tensor field of valence 1.
Scalar lapconf_comp
Lapconf function generated by the companion star.
virtual void std_spectral_base()
Sets the standard spectal bases of decomposition for each component.
void raccord(int n)
Performs the matching of the nucleus with respect to the first shell.
Tensor d_shift_auto_bh
Derivative of the part of the shift vector from the analytic background.
Scalar confo_auto_rs
Part of the conformal factor from the numerical computation.
virtual const Sym_tensor & con() const
Read-only access to the contravariant representation.
const Vector & get_shift_auto() const
Returns the part of the shift vector generated by the star.
Scalar lapse_tot
Total lapse function.
Scalar lapse_auto
Lapse function of the "black hole" part.
Vector shift_auto_bh
Part of the shift vector from the analytic background.
const Base_vect * get_triad() const
Returns the vectorial basis (triad) on which the components are defined.
Scalar confo_comp
Conformal factor generated by the companion star.
Vector shift_auto_rs
Part of the shift vector from the numerical computation.
Vector shift_comp
Shift vector generated by the companion star.
const Scalar & deriv(int i) const
Returns of *this , where .
Cmp pow(const Cmp &, int)
Power .
const Scalar & get_confo_auto() const
Returns the part of the conformal factor generated by the star.
Vector d_lapconf_auto_rs
Derivative of the part of the lapconf function from the numerical computation.
Tensor d_shift_auto_rs
Derivative of the part of the shift vector from the numerical computation.
Vector d_confo_auto_bh
Derivative of the part of the conformal factor from the analytic background.
Class for black holes in black hole-neutron star binaries.
Scalar lapconf_auto
Lapconf function generated by the black hole.
virtual void inc_dzpuis(int inc=1)
Increases by inc units the value of dzpuis and changes accordingly the values in the compactified ext...
Scalar confo_tot
Total conformal factor.
void import(const Scalar &ci)
Assignment to another Scalar defined on a different mapping.
const Base_vect_cart & get_bvect_cart() const
Returns the Cartesian basis associated with the coordinates (x,y,z) of the mapping, i.e.
Scalar & set(const Itbl &ind)
Returns the value of a component (read/write version).
Metric_flat flat
Flat metric defined on the mapping (Spherical components with respect to the mapping of the black hol...
Scalar lapconf_auto_rs
Part of the lapconf function from the numerical computation.
bool bc_lapconf_fs
true for the first type BC for the lapconf function, false for the second type BH ...
void set_triad(const Base_vect &new_triad)
Assigns a new vectorial basis (triad) of decomposition.
Vector d_confo_auto_rs
Derivative of the part of the conformal factor from the numerical computation.
virtual void set_etat_zero()
Sets the logical state of all components to ETATZERO (zero state).
bool bc_lapconf_nd
true for the Neumann type BC for the lapconf function, false for the Dirichlet type BH ...
Scalar & set(int)
Read/write access to a component.
const Scalar r_coord(bool neumann, bool first) const
Expresses the areal radial coordinate by that in spatially isotropic coordinates. ...
void update_metric_bhns(const Star_bhns &star, const Hole_bhns &hole_prev, double relax)
Computes metric coefficients from known potentials with relaxation when the companion is a black hole...
virtual void std_spectral_base()
Sets the standard spectal bases of decomposition for each component.
virtual void del_deriv() const
Deletes all the derived quantities.
Scalar lapconf_tot
Total lapconf function.
Coord r
r coordinate centered on the grid
Scalar confo_auto_bh
Part of the conformal factor from the analytic background.