64 #include "blackhole.h" 86 cout <<
"!!!!! WARNING: Not yet available !!!!!" << endl ;
157 for (
int j=0; j<nt; j++) {
158 for (
int k=0; k<np; k++) {
171 cout <<
"!!!!! WARNING: Not yet prepared !!!!!" << endl ;
178 cout <<
"!!!!! WARNING: Not yet prepared !!!!!" << endl ;
243 cout <<
"!!!!! WARNING: Not yet available !!!!!" << endl ;
265 for (
int j=0; j<nt; j++) {
266 for (
int k=0; k<np; k++) {
271 bc.
base = *bases[0] ;
274 for (
int i=0; i<3; i++)
313 cout <<
"!!!!! WARNING: Not yet available !!!!!" << endl ;
334 for (
int j=0; j<nt; j++) {
335 for (
int k=0; k<np; k++) {
340 bc.
base = *bases[1] ;
343 for (
int i=0; i<3; i++)
376 cout <<
"!!!!! WARNING: Not yet available !!!!!" << endl ;
394 for (
int j=0; j<nt; j++) {
395 for (
int k=0; k<np; k++) {
400 bc.
base = *bases[2] ;
403 for (
int i=0; i<3; i++)
423 double mass = ggrav *
mass_bh ;
448 cout <<
"!!!!! WARNING: Not yet available !!!!!" << endl ;
485 divshift =
shift(1).deriv(1) +
shift(2).deriv(2)
486 +
shift(3).deriv(3) ;
497 Scalar tmp2 = -3.*lldllsh ;
509 for (
int j=0; j<nt; j++) {
510 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.
const Valeur bc_lapconf(bool neumann, bool first) const
Boundary condition on the apparent horizon of the black hole for the lapse function: 2-D Valeur...
int get_np(int l) const
Returns the number of points in the azimuthal direction ( ) in domain no. l.
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
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.
void std_base_scal()
Sets the bases for spectral expansions (member base ) to the standard ones for a scalar.
void set_dzpuis(int)
Modifies the dzpuis flag.
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 confo
Conformal factor generated by the black hole.
const Valeur bc_shift_x(double omega_r) const
Boundary condition on the apparent horizon of the black hole for the shift vector of the direction: ...
const Mg3d * get_angu() const
Returns the pointer on the associated angular grid.
Base_val base
Bases on which the spectral expansion is performed.
int get_dzpuis() const
Returns dzpuis.
const Valeur bc_shift_y(double omega_r) const
Boundary condition on the apparent horizon of the black hole for the shift vector of the direction: ...
Cmp pow(const Cmp &, int)
Power .
const Valeur bc_shift_z() const
Boundary condition on the apparent horizon of the black hole for the shift vector of the direction: ...
Bases of the spectral expansions.
Coord y
y coordinate centered on the grid
Vector shift
Shift vector generated by the black hole.
const Scalar & dsdr() const
Returns of *this .
Coord x
x coordinate centered on the grid
int get_nt(int l) const
Returns the number of points in the co-latitude direction ( ) in domain no. l.
Scalar lapconf
A function (lapse function * conformal factor) lapconf generated by the black hole.
const Valeur bc_confo() const
Boundary condition on the apparent horizon of the black hole for the conformal factor: 2-D Valeur...
Coord r
r coordinate centered on the grid
Tbl & set(int l)
Read/write of the value in a given domain (configuration space).