29char blackhole_bc_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Black_hole/blackhole_bc.C,v 1.5 2014/10/13 08:52:45 j_novak Exp $" ;
83 cout <<
"!!!!! WARNING: Not yet available !!!!!" << endl ;
154 for (
int j=0; j<nt; j++) {
155 for (
int k=0; k<np; k++) {
168 cout <<
"!!!!! WARNING: Not yet prepared !!!!!" << endl ;
175 cout <<
"!!!!! WARNING: Not yet prepared !!!!!" << endl ;
240 cout <<
"!!!!! WARNING: Not yet available !!!!!" << endl ;
262 for (
int j=0; j<nt; j++) {
263 for (
int k=0; k<np; k++) {
268 bc.
base = *bases[0] ;
271 for (
int i=0; i<3; i++)
310 cout <<
"!!!!! WARNING: Not yet available !!!!!" << endl ;
331 for (
int j=0; j<nt; j++) {
332 for (
int k=0; k<np; k++) {
337 bc.
base = *bases[1] ;
340 for (
int i=0; i<3; i++)
373 cout <<
"!!!!! WARNING: Not yet available !!!!!" << endl ;
391 for (
int j=0; j<nt; j++) {
392 for (
int k=0; k<np; k++) {
397 bc.
base = *bases[2] ;
400 for (
int i=0; i<3; i++)
420 double mass = ggrav *
mass_bh ;
445 cout <<
"!!!!! WARNING: Not yet available !!!!!" << endl ;
482 divshift =
shift(1).deriv(1) +
shift(2).deriv(2)
483 +
shift(3).deriv(3) ;
494 Scalar tmp2 = -3.*lldllsh ;
506 for (
int j=0; j<nt; j++) {
507 for (
int k=0; k<np; k++) {
Bases of the spectral expansions.
const Valeur bc_confo() const
Boundary condition on the apparent horizon of the black hole for the conformal factor: 2-D Valeur.
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 \fx\f direct...
Map & mp
Mapping associated with the black hole.
Scalar lapconf
A function (lapse function * conformal factor) lapconf generated by the black hole.
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.
Vector shift
Shift vector generated by the black hole.
const Valeur bc_shift_z() const
Boundary condition on the apparent horizon of the black hole for the shift vector of the \fz\f direct...
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 \fy\f direct...
bool kerrschild
true for a Kerr-Schild background, false for a conformally flat background
Scalar confo
Conformal factor generated by the black hole.
double mass_bh
Gravitational mass of BH.
Coord y
y coordinate centered on the grid
Coord r
r coordinate centered on the grid
Coord x
x coordinate centered on the grid
const Mg3d * get_mg() const
Gives the Mg3d on which the mapping is defined.
const Mg3d * get_angu() const
Returns the pointer on the associated angular grid.
int get_np(int l) const
Returns the number of points in the azimuthal direction ( ) in domain no. l.
int get_nt(int l) const
Returns the number of points in the co-latitude direction ( ) in domain no. l.
Base_val ** std_base_vect_cart() const
Returns the standard spectral bases for the Cartesian components of a vector.
Tensor field of valence 0 (or component of a tensorial field).
int get_dzpuis() const
Returns dzpuis.
virtual void std_spectral_base()
Sets the spectral bases of the Valeur va to the standard ones for a scalar field.
double val_grid_point(int l, int k, int j, int i) const
Returns the value of the field at a specified grid point.
const Scalar & dsdr() const
Returns of *this .
void set_dzpuis(int)
Modifies the dzpuis flag.
Values and coefficients of a (real-value) function.
Tbl & set(int l)
Read/write of the value in a given domain (configuration space).
Base_val base
Bases on which the spectral expansion is performed.
void std_base_scal()
Sets the bases for spectral expansions (member base ) to the standard ones for a scalar.
Cmp pow(const Cmp &, int)
Power .
Standard units of space, time and mass.