28char altBH_QI_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Compobj/altBH_QI.C,v 1.5 2014/10/13 08:52:49 j_novak Exp $" ;
75 ifstream file(file_name) ;
77 cerr <<
"Problem in opening the file " << file_name << endl ;
92 file.ignore(1000,
'\n') ;
93 file.ignore(1000,
'\n') ;
95 cout <<
"nz : " << nz << endl ;
101 for (
int l=1; l<nz; l++) {
102 cout <<
"l = " << l << endl ;
106 double* r_iso =
new double[nr] ;
107 double* r_areal =
new double[nr] ;
108 double* psi4 =
new double[nr] ;
109 double* alpha =
new double[nr] ;
110 double* Krphi =
new double[nr] ;
111 double* beta_phi =
new double[nr] ;
115 r_areal[nr-1] = __infinity ;
116 r_iso[nr-1] = __infinity ;
122 for (
int i=0; i<nr_max; i++) {
130 file >> beta_phi[i] ;
131 cout <<
"r_iso, psi4, beta_phi : " << r_iso[i] <<
" " << psi4[i] <<
" " << beta_phi[i] << endl ;
134 if (l==1) r_inner = r_iso[0] ;
136 for (
int k=0; k<np; k++) {
137 for (
int j=0; j<nt; j++) {
138 for (
int i=0; i<nr; i++) {
254 ost << endl <<
"Alternative black hole spacetime in quasi-isotropic coordinates (class AltBH_QI) " << endl ;
273 if ( (
mp.
get_mg())->get_np(0) == 1) {
Alternative black hole spacetime in Quasi-Isotropic coordinates (under development).
char description2[256]
String describing the model.
void operator=(const AltBH_QI &)
Assignment to another AltBH_QI.
AltBH_QI(Map &mp_i, const char *file_name, double a_spin_i)
Standard constructor.
virtual ~AltBH_QI()
Destructor.
char description1[256]
String describing the model.
virtual void del_deriv() const
Deletes all the derived quantities.
Scalar krphi
K_{(r)(phi)} read in the file.
virtual ostream & operator>>(ostream &) const
Operator >> (virtual function called by the operator <<).
double a_spin
Spin parameter of the model.
virtual void set_der_0x0() const
Sets to 0x0 all the pointers on derived quantities.
virtual void extrinsic_curvature()
Computation of the extrinsic curvature.
virtual void sauve(FILE *) const
Save in a file.
Base class for axisymmetric stationary compact objects in Quasi-Isotropic coordinates (under developm...
void operator=(const Compobj_QI &)
Assignment to another Compobj_QI.
Scalar nphi
Metric coefficient .
virtual void del_deriv() const
Deletes all the derived quantities.
virtual ostream & operator>>(ostream &) const
Operator >> (virtual function called by the operator <<).
Scalar b_car
Square of the metric factor B.
Scalar bbb
Metric factor B.
Scalar a_car
Square of the metric factor A.
Sym_tensor kk
Extrinsic curvature tensor
Scalar nn
Lapse function N .
virtual void extrinsic_curvature()
Computation of the extrinsic curvature.
Map & mp
Mapping describing the coordinate system (r,theta,phi)
Base class for pure radial mappings.
Base class for coordinate mappings.
virtual void homothetie(double lambda)=0
Sets a new radial scale.
virtual double val_r(int l, double xi, double theta, double pphi) const =0
Returns the value of the radial coordinate r for a given in a given domain.
const Mg3d * get_mg() const
Gives the Mg3d on which the mapping is defined.
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.
int get_nzone() const
Returns the number of domains.
int get_nr(int l) const
Returns the number of points in the radial direction ( ) in domain no. l.
Tensor field of valence 0 (or component of a tensorial field).
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 annule(int l_min, int l_max)
Sets the Scalar to zero in several domains.
void mult_sint()
Multiplication by .
virtual void allocate_all()
Sets the logical state to ETATQCQ (ordinary state) and performs the memory allocation of all the elem...
double & set_grid_point(int l, int k, int j, int i)
Setting the value of the field at a given grid point.
Cmp sqrt(const Cmp &)
Square root.
Scalar & set(const Itbl &ind)
Returns the value of a component (read/write version).
Standard units of space, time and mass.