29char hole_bhns_bc_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Hole_bhns/hole_bhns_bc.C,v 1.4 2014/10/13 08:53:00 j_novak Exp $" ;
107 cout <<
"!!!!! WARNING: Not yet available !!!!!" << endl ;
124 cout <<
"!!!!! WARNING: Not yet available !!!!!" << endl ;
143 cout <<
"!!!!! WARNING: Not yet available !!!!!" << endl ;
148 cout <<
"!!!!! WARNING: Not yet prepared !!!!!" << endl ;
158 cout <<
"!!!!! WARNING: Not yet available !!!!!" << endl ;
162 cout <<
"!!!!! WARNING: Not yet prepared !!!!!" << endl ;
172 for (
int j=0; j<nt; j++) {
173 for (
int k=0; k<np; k++) {
208 double mass = ggrav *
mass_bh ;
217 cout <<
"!!!!! WARNING: Not yet available !!!!!" << endl ;
232 cc = 2. * (
sqrt(13.) - 1.) / 3. ;
244 cout <<
"!!!!! WARNING: Not yet prepared !!!!!" << endl ;
250 cout <<
"!!!!! WARNING: Not yet prepared !!!!!" << endl ;
268 + (ome_orb -
omega_spin) * yy + ome_orb * (ori_y_bh - y_rot) ;
276 for (
int j=0; j<nt; j++) {
277 for (
int k=0; k<np; k++) {
282 bc.
base = *bases[0] ;
284 for (
int i=0; i<3; i++)
318 double mass = ggrav *
mass_bh ;
327 cout <<
"!!!!! WARNING: Not yet available !!!!!" << endl ;
342 cc = 2. * (
sqrt(13.) - 1.) / 3. ;
354 cout <<
"!!!!! WARNING: Not yet prepared !!!!!" << endl ;
360 cout <<
"!!!!! WARNING: Not yet prepared !!!!!" << endl ;
378 - (ome_orb -
omega_spin) * xx - ome_orb * (ori_x_bh - x_rot) ;
386 for (
int j=0; j<nt; j++) {
387 for (
int k=0; k<np; k++) {
392 bc.
base = *bases[1] ;
394 for (
int i=0; i<3; i++)
422 double mass = ggrav *
mass_bh ;
430 cout <<
"!!!!! WARNING: Not yet available !!!!!" << endl ;
445 cc = 2. * (
sqrt(13.) - 1.) / 3. ;
457 cout <<
"!!!!! WARNING: Not yet prepared !!!!!" << endl ;
463 cout <<
"!!!!! WARNING: Not yet prepared !!!!!" << endl ;
486 for (
int j=0; j<nt; j++) {
487 for (
int k=0; k<np; k++) {
492 bc.
base = *bases[2] ;
494 for (
int i=0; i<3; i++)
504 double y_rot)
const {
532 ll.
set(1) = st % cp ;
533 ll.
set(2) = st % sp ;
564 Scalar tmp3 = -3.*lldllsh ;
571 double mass = ggrav *
mass_bh ;
575 cout <<
"!!!!! WARNING: Not yet available !!!!!" << endl ;
590 cc = 2. * (
sqrt(13.) - 1.) / 3. ;
602 cout <<
"!!!!! WARNING: Not yet prepared !!!!!" << endl ;
608 cout <<
"!!!!! WARNING: Not yet prepared !!!!!" << endl ;
631 -
sqrt(r_are)) / rr ;
644 for (
int j=0; j<nt; j++) {
645 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 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.
bool kerrschild
true for a Kerr-Schild background, false for a conformally flat background
double mass_bh
Gravitational mass of BH.
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_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 \fx\f direct...
Vector shift_auto_rs
Part of the shift vector from the numerical computation.
double omega_spin
Spin angular velocity of the black hole.
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 \fy\f direct...
bool bc_lapconf_fs
true for the first type BC for the lapconf function, false for the second type BH
Scalar confo_auto_rs
Part of the conformal factor from the numerical computation.
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...
Vector d_lapconf_comp
Derivative of the lapconf function generated by the companion star.
Scalar confo_comp
Conformal factor generated by the companion star.
Scalar lapconf_comp
Lapconf function generated by the companion star.
Scalar lapconf_auto_rs
Part of the lapconf function from the numerical computation.
Vector d_confo_comp
Derivative of the conformal factor generated by the companion star.
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.
Tensor d_shift_comp
Derivative of the shift vector generated by the companion star.
Vector shift_comp
Shift vector generated by the companion star.
Scalar confo_tot
Total conformal factor.
const Base_vect_cart & get_bvect_cart() const
Returns the Cartesian basis associated with the coordinates (x,y,z) of the mapping,...
Coord y
y coordinate centered on the grid
Coord r
r coordinate centered on the grid
double get_ori_y() const
Returns the y coordinate of the origin.
double get_ori_x() const
Returns the x coordinate of the origin.
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).
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.
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 Scalar & dsdr() const
Returns of *this .
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...
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.
Tensor field of valence 1.
virtual void std_spectral_base()
Sets the standard spectal bases of decomposition for each component.
Scalar & set(int)
Read/write access to a component.
Cmp sqrt(const Cmp &)
Square root.
Cmp pow(const Cmp &, int)
Power .
virtual void set_etat_qcq()
Sets the logical state of all components to ETATQCQ (ordinary state).
Standard units of space, time and mass.