61 #include "hole_bhns.h" 110 cout <<
"!!!!! WARNING: Not yet available !!!!!" << endl ;
127 cout <<
"!!!!! WARNING: Not yet available !!!!!" << endl ;
146 cout <<
"!!!!! WARNING: Not yet available !!!!!" << endl ;
151 cout <<
"!!!!! WARNING: Not yet prepared !!!!!" << endl ;
161 cout <<
"!!!!! WARNING: Not yet available !!!!!" << endl ;
165 cout <<
"!!!!! WARNING: Not yet prepared !!!!!" << endl ;
175 for (
int j=0; j<nt; j++) {
176 for (
int k=0; k<np; k++) {
211 double mass = ggrav *
mass_bh ;
220 cout <<
"!!!!! WARNING: Not yet available !!!!!" << endl ;
235 cc = 2. * (
sqrt(13.) - 1.) / 3. ;
247 cout <<
"!!!!! WARNING: Not yet prepared !!!!!" << endl ;
253 cout <<
"!!!!! WARNING: Not yet prepared !!!!!" << endl ;
271 + (ome_orb -
omega_spin) * yy + ome_orb * (ori_y_bh - y_rot) ;
279 for (
int j=0; j<nt; j++) {
280 for (
int k=0; k<np; k++) {
285 bc.
base = *bases[0] ;
287 for (
int i=0; i<3; i++)
321 double mass = ggrav *
mass_bh ;
330 cout <<
"!!!!! WARNING: Not yet available !!!!!" << endl ;
345 cc = 2. * (
sqrt(13.) - 1.) / 3. ;
357 cout <<
"!!!!! WARNING: Not yet prepared !!!!!" << endl ;
363 cout <<
"!!!!! WARNING: Not yet prepared !!!!!" << endl ;
381 - (ome_orb -
omega_spin) * xx - ome_orb * (ori_x_bh - x_rot) ;
389 for (
int j=0; j<nt; j++) {
390 for (
int k=0; k<np; k++) {
395 bc.
base = *bases[1] ;
397 for (
int i=0; i<3; i++)
425 double mass = ggrav *
mass_bh ;
433 cout <<
"!!!!! WARNING: Not yet available !!!!!" << endl ;
448 cc = 2. * (
sqrt(13.) - 1.) / 3. ;
460 cout <<
"!!!!! WARNING: Not yet prepared !!!!!" << endl ;
466 cout <<
"!!!!! WARNING: Not yet prepared !!!!!" << endl ;
489 for (
int j=0; j<nt; j++) {
490 for (
int k=0; k<np; k++) {
495 bc.
base = *bases[2] ;
497 for (
int i=0; i<3; i++)
507 double y_rot)
const {
535 ll.set(1) = st % cp ;
536 ll.
set(2) = st % sp ;
567 Scalar tmp3 = -3.*lldllsh ;
574 double mass = ggrav *
mass_bh ;
578 cout <<
"!!!!! WARNING: Not yet available !!!!!" << endl ;
593 cc = 2. * (
sqrt(13.) - 1.) / 3. ;
605 cout <<
"!!!!! WARNING: Not yet prepared !!!!!" << endl ;
611 cout <<
"!!!!! WARNING: Not yet prepared !!!!!" << endl ;
634 -
sqrt(r_are)) / rr ;
647 for (
int j=0; j<nt; j++) {
648 for (
int k=0; k<np; k++) {
Map & mp
Mapping associated with the black hole.
Base_val ** std_base_vect_cart() const
Returns the standard spectral bases for the Cartesian components of a vector.
virtual void set_etat_qcq()
Sets the logical state of all components to ETATQCQ (ordinary state).
const Valeur bc_shift_z() const
Boundary condition on the apparent horizon of the black hole for the shift vector of the direction: ...
int get_np(int l) const
Returns the number of points in the azimuthal direction ( ) in domain no. l.
Cmp sqrt(const Cmp &)
Square root.
double get_ori_y() const
Returns the y coordinate of the origin.
double mass_bh
Gravitational mass of BH.
Standard units of space, time and mass.
const Mg3d * get_mg() const
Gives the Mg3d on which the mapping is defined.
Tensor field of valence 0 (or component of a tensorial field).
bool kerrschild
true for a Kerr-Schild background, false for a conformally flat background
double get_ori_x() const
Returns the x coordinate of the origin.
virtual void std_spectral_base()
Sets the spectral bases of the Valeur va to the standard ones for a scalar field. ...
Values and coefficients of a (real-value) function.
Tensor d_shift_comp
Derivative of the shift vector generated by the companion star.
Tensor field of valence 1.
void std_base_scal()
Sets the bases for spectral expansions (member base ) to the standard ones for a scalar.
double val_grid_point(int l, int k, int j, int i) const
Returns the value of the field at a specified grid point.
Scalar lapconf_comp
Lapconf function generated by the companion star.
const Valeur bc_shift_x(double ome_orb, double y_rot) const
Boundary condition on the apparent horizon of the black hole for the shift vector of the direction: ...
const Valeur bc_shift_y(double ome_orb, double x_rot) const
Boundary condition on the apparent horizon of the black hole for the shift vector of the direction: ...
Scalar confo_auto_rs
Part of the conformal factor from the numerical computation.
Vector d_lapconf_comp
Derivative of the lapconf function generated by the companion star.
virtual void inc_dzpuis(int inc=1)
Increases by inc units the value of dzpuis and changes accordingly the values of the Scalar in the co...
const Mg3d * get_angu() const
Returns the pointer on the associated angular grid.
Scalar confo_comp
Conformal factor generated by the companion star.
Base_val base
Bases on which the spectral expansion is performed.
Vector shift_auto_rs
Part of the shift vector from the numerical computation.
Vector shift_comp
Shift vector generated by the companion star.
Cmp pow(const Cmp &, int)
Power .
Scalar confo_tot
Total conformal factor.
Bases of the spectral expansions.
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).
Coord y
y coordinate centered on the grid
double omega_spin
Spin angular velocity of the black hole.
Vector d_confo_comp
Derivative of the conformal factor generated by the companion star.
Scalar lapconf_auto_rs
Part of the lapconf function from the numerical computation.
const Scalar & dsdr() const
Returns of *this .
Coord x
x coordinate centered on the grid
bool bc_lapconf_fs
true for the first type BC for the lapconf function, false for the second type BH ...
bool bc_lapconf_nd
true for the Neumann type BC for the lapconf function, false for the Dirichlet type BH ...
int get_nt(int l) const
Returns the number of points in the co-latitude direction ( ) in domain no. l.
const Scalar r_coord(bool neumann, bool first) const
Expresses the areal radial coordinate by that in spatially isotropic coordinates. ...
virtual void dec_dzpuis(int dec=1)
Decreases by dec units the value of dzpuis and changes accordingly the values of the Scalar in the co...
const Valeur bc_lapconf() const
Boundary condition on the apparent horizon of the black hole for the lapconf function: 2-D Valeur...
const Valeur bc_confo() const
Boundary condition on the apparent horizon of the black hole for the conformal factor: 2-D Valeur...
Scalar lapconf_tot
Total lapconf function.
Coord r
r coordinate centered on the grid
Tbl & set(int l)
Read/write of the value in a given domain (configuration space).