29char 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 $" ;
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++) {
const Scalar r_coord(bool neumann, bool first) const
Expresses the areal radial coordinate by that in spatially isotropic coordinates.
Map & mp
Mapping associated with the black hole.
Metric_flat flat
Flat metric defined on the mapping (Spherical components with respect to the mapping of the black hol...
bool kerrschild
true for a Kerr-Schild background, false for a conformally flat background
double mass_bh
Gravitational mass of BH.
Class for black holes in black hole-neutron star binaries.
Scalar confo_auto
Conformal factor generated by the black hole.
Scalar lapconf_auto
Lapconf function generated by the black hole.
Vector shift_tot
Total shift vector ;.
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 d_lapconf_auto_rs
Derivative of the part of the lapconf function from the numerical computation.
Scalar confo_auto_bh
Part of the conformal factor from the analytic background.
Tensor d_shift_auto_rs
Derivative of the part of the shift vector from the numerical computation.
Vector shift_auto_rs
Part of the shift vector from the numerical computation.
virtual void del_deriv() const
Deletes all the derived quantities.
Tensor d_shift_auto_bh
Derivative of the part of the shift vector from the analytic background.
bool bc_lapconf_fs
true for the first type BC for the lapconf function, false for the second type BH
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...
Vector d_lapconf_auto_bh
Derivative of the part of the lapconf function from the analytic background.
Vector d_confo_auto_rs
Derivative of the part of the conformal factor from the numerical computation.
Scalar confo_auto_rs
Part of the conformal factor from the numerical computation.
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.
Scalar confo_comp
Conformal factor generated by the companion star.
Scalar lapconf_comp
Lapconf function generated by the companion star.
Scalar lapconf_auto_bh
Part of the lapconf function from the analytic background.
Vector d_confo_auto
Derivative of the conformal factor generated by the black hole.
Vector shift_auto
Shift vector generated by the black hole.
Scalar lapconf_auto_rs
Part of the lapconf function from the numerical computation.
bool bc_lapconf_nd
true for the Neumann type BC for the lapconf function, false for the Dirichlet type BH
Scalar lapconf_tot
Total lapconf function.
Vector d_confo_auto_bh
Derivative of the part of the conformal factor from the analytic background.
Vector shift_comp
Shift vector generated by the companion star.
Scalar confo_tot
Total conformal factor.
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,...
Coord r
r coordinate centered on the grid
virtual const Sym_tensor & con() const
Read-only access to the contravariant representation.
Tensor field of valence 0 (or component of a tensorial field).
void raccord(int n)
Performs the matching of the nucleus with respect to the first shell.
const Scalar & deriv(int i) const
Returns of *this , where .
virtual void set_etat_qcq()
Sets the logical state to ETATQCQ (ordinary state).
virtual void std_spectral_base()
Sets the spectral bases of the Valeur va to the standard ones for a scalar field.
virtual void set_etat_zero()
Sets the logical state to ETATZERO (zero).
void import(const Scalar &ci)
Assignment to another Scalar defined on a different mapping.
Class for stars in black hole-neutron star binaries.
const Scalar & get_confo_auto() const
Returns the part of the conformal factor generated by the star.
const Vector & get_shift_auto() const
Returns the part of the shift vector generated by the star.
const Scalar & get_lapconf_auto() const
Returns the part of the lapconf function generated by the star.
const Map & get_mp() const
Returns the mapping.
Tensor field of valence 1.
virtual void std_spectral_base()
Sets the standard spectal bases of decomposition for each component.
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 sqrt(const Cmp &)
Square root.
Cmp pow(const Cmp &, int)
Power .
void set_triad(const Base_vect &new_triad)
Assigns a new vectorial basis (triad) of decomposition.
virtual void inc_dzpuis(int inc=1)
Increases by inc units the value of dzpuis and changes accordingly the values in the compactified ext...
const Base_vect * get_triad() const
Returns the vectorial basis (triad) on which the components are defined.
void annule_domain(int l)
Sets the Tensor to zero in a given domain.
virtual void set_etat_qcq()
Sets the logical state of all components to ETATQCQ (ordinary state).
virtual void set_etat_zero()
Sets the logical state of all components to ETATZERO (zero state).
Scalar & set(const Itbl &ind)
Returns the value of a component (read/write version).
virtual void std_spectral_base()
Sets the standard spectal bases of decomposition for each component.
Standard units of space, time and mass.