29 char hole_bhns_upmetr_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Hole_bhns/hole_bhns_upmetr.C,v 1.3 2014/10/13 08:53:00 j_novak Exp $" ;
55 #include "hole_bhns.h"
56 #include "star_bhns.h"
57 #include "utilitaires.h"
131 double relax_jm1 = 1. - relax ;
145 rr.std_spectral_base() ;
161 ll.
set(1) = st % cp ;
162 ll.
set(2) = st % sp ;
190 for (
int i=1; i<=3; i++)
196 for (
int i=1; i<=3; i++) {
208 for (
int i=1; i<=3; i++) {
209 for (
int j=1; j<=3; j++) {
217 for (
int i=1; i<=3; i++) {
218 for (
int j=1; j<=3; j++) {
235 for (
int i=1; i<=3; i++)
241 for (
int i=1; i<=3; i++)
260 cc = 2. * (
sqrt(13.) - 1.) / 3. ;
272 cout <<
"!!!!! WARNING: Not yet prepared !!!!!" << endl ;
278 cout <<
"!!!!! WARNING: Not yet prepared !!!!!" << endl ;
293 + cc*cc*
pow(mass/r_are/rr, 4.)) *
sqrt(r_are) ;
306 for (
int i=1; i<=3; i++)
314 for (
int i=1; i<=3; i++)
320 for (
int i=1; i<=3; i++) {
322 * (mass/r_are/rr - 2.*cc*cc*
pow(mass/r_are/rr,4.))
325 * (
sqrt(1. - 2.*mass/r_are/rr + cc*cc*
pow(mass/r_are/rr,4.))-1.)
326 *
sqrt(1. - 2.*mass/r_are/rr + cc*cc*
pow(mass/r_are/rr,4.))
336 for (
int i=1; i<=3; i++)
340 for (
int i=1; i<=3; i++) {
341 for (
int j=1; j<=3; j++) {
349 for (
int i=1; i<=3; i++) {
350 for (
int j=1; j<=3; j++) {
353 -3.*
sqrt(1. - 2.*mass/r_are/rr
354 +cc*cc*
pow(mass/r_are/rr,4.))
366 for (
int i=1; i<=3; i++) {
367 for (
int j=1; j<=3; j++) {
373 for (
int i=1; i<=3; i++)
379 for (
int i=1; i<=3; i++) {
381 *(
sqrt(1. - 2.*mass/r_are/rr +cc*cc*
pow(mass/r_are/rr,4.)) - 1.)
391 for (
int i=1; i<=3; i++)
413 for (
int i=1; i<=3; i++) {
420 for (
int i=1; i<=3; i++) {
Scalar lapconf_auto_bh
Part of the lapconf function from the analytic background.
void raccord(int n)
Performs the matching of the nucleus with respect to the first shell.
Scalar confo_auto_bh
Part of the conformal factor from the analytic background.
virtual const Sym_tensor & con() const
Read-only access to the contravariant representation.
Vector d_lapconf_auto_bh
Derivative of the part of the lapconf function from the analytic background.
Vector d_confo_auto
Derivative of the conformal factor generated by the black hole.
virtual void del_deriv() const
Deletes all the derived quantities.
Metric_flat flat
Flat metric defined on the mapping (Spherical components with respect to the mapping of the black hol...
virtual void set_etat_zero()
Sets the logical state of all components to ETATZERO (zero state).
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.
const Scalar & deriv(int i) const
Returns of *this , where .
Vector shift_tot
Total shift vector ;.
Scalar confo_auto
Conformal factor generated by the black hole.
Scalar confo_auto_rs
Part of the conformal factor from the numerical computation.
double mass_bh
Gravitational mass of BH.
Coord r
r coordinate centered on the grid
Class for black holes in black hole-neutron star binaries.
const Scalar & get_confo_auto() const
Returns the part of the conformal factor generated by the star.
Vector shift_auto
Shift vector generated by the black hole.
Map & mp
Mapping associated with the black hole.
void import(const Scalar &ci)
Assignment to another Scalar defined on a different mapping.
Tensor field of valence 0 (or component of a tensorial field).
Cmp pow(const Cmp &, int)
Power .
void annule_domain(int l)
Sets the Tensor to zero in a given domain.
virtual void change_triad(const Base_vect &)
Sets a new vectorial basis (triad) of decomposition and modifies the components accordingly.
Scalar lapse_tot
Total lapse function.
Scalar lapse_auto
Lapse function of the "black hole" part.
virtual void set_etat_qcq()
Sets the logical state of all components to ETATQCQ (ordinary state).
Standard units of space, time and mass.
Scalar & set(const Itbl &ind)
Returns the value of a component (read/write version).
Scalar lapconf_comp
Lapconf function generated by the companion star.
Scalar confo_comp
Conformal factor generated by the companion star.
Tensor d_shift_auto_bh
Derivative of the part of the shift vector from the analytic background.
Scalar & set(int)
Read/write access to a component.
Class for stars in black hole-neutron star binaries.
Vector shift_auto_rs
Part of the shift vector from the numerical computation.
Vector shift_auto_bh
Part of the shift vector from the analytic background.
virtual void std_spectral_base()
Sets the standard spectal bases of decomposition for each component.
Vector shift_comp
Shift vector generated by the companion star.
const Vector & get_shift_auto() const
Returns the part of the shift vector generated by the star.
Vector d_confo_auto_bh
Derivative of the part of the conformal factor from the analytic background.
const Scalar r_coord(bool neumann, bool first) const
Expresses the areal radial coordinate by that in spatially isotropic coordinates.
virtual void set_etat_qcq()
Sets the logical state to ETATQCQ (ordinary state).
Tensor d_shift_auto_rs
Derivative of the part of the shift vector from the numerical computation.
const Base_vect_cart & get_bvect_cart() const
Returns the Cartesian basis associated with the coordinates (x,y,z) of the mapping,...
Cmp sqrt(const Cmp &)
Square root.
virtual void set_etat_zero()
Sets the logical state to ETATZERO (zero).
Tensor field of valence 1.
bool bc_lapconf_fs
true for the first type BC for the lapconf function, false for the second type BH
const Base_vect * get_triad() const
Returns the vectorial basis (triad) on which the components are defined.
bool bc_lapconf_nd
true for the Neumann type BC for the lapconf function, false for the Dirichlet type BH
Vector d_lapconf_auto_rs
Derivative of the part of the lapconf function from the numerical computation.
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 lapconf_auto
Lapconf function generated by the black hole.
virtual void std_spectral_base()
Sets the spectral bases of the Valeur va to the standard ones for a scalar field.
Scalar confo_tot
Total conformal factor.
const Map & get_mp() const
Returns the mapping.
Scalar lapconf_tot
Total lapconf function.
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...
void set_triad(const Base_vect &new_triad)
Assigns a new vectorial basis (triad) of decomposition.
const Scalar & get_lapconf_auto() const
Returns the part of the lapconf function generated by the star.
Base class for coordinate mappings.
Scalar lapconf_auto_rs
Part of the lapconf function from the numerical computation.
bool kerrschild
true for a Kerr-Schild background, false for a conformally flat background
Vector d_confo_auto_rs
Derivative of the part of the conformal factor from the numerical computation.
virtual void std_spectral_base()
Sets the standard spectal bases of decomposition for each component.