67 assert (
rayon ==
sqrt (masse*masse-moment*moment)/2.) ;
70 omega = moment/2/masse/(masse+
sqrt(masse*masse-moment*moment)) ;
74 grand_r =
mp.
r + (masse*masse-moment*moment)/4/
mp.
r + masse ;
78 sigma = moment*moment*
mp.
cost*
mp.
cost + grand_r*grand_r ;
82 grand_a = 1 + 2*masse/
mp.
r +
91 n_phi = (2*moment*masse*grand_r) / (sigma*(grand_r*grand_r+moment*moment)
99 carre = 1-(2*masse*grand_r)/sigma + (4*moment*moment*masse*masse
101 (sigma*sigma*(grand_r*grand_r+moment*moment)+2*moment*moment*sigma*masse*
133 for (
int i=0 ; i<2 ; i++) {
Component of a tensorial field *** DEPRECATED : use class Scalar instead ***.
void init_kerr(double masse, double moment)
Set the inital values to those of Kerr.
Cmp sqrt(const Cmp &)
Square root.
void annule(int l)
Sets the Cmp to zero in a given domain.
const Valeur & mult_sp() const
Returns applied to *this.
void set_std_base()
Set the standard spectal basis of decomposition for each component.
const Mg3d * get_mg() const
Gives the Mg3d on which the mapping is defined.
Tenseur psi_auto
Part of generated by the hole.
Values and coefficients of a (real-value) function.
void inc_dzpuis()
dzpuis += 1 ;
Tenseur shift_auto
Part of generated by the hole.
void mult_r()
Multiplication by r everywhere.
void set_etat_zero()
Sets the logical state to ETATZERO (zero).
Cmp & set()
Read/write for a scalar (see also operator=(const Cmp&) ).
Map_af & mp
Affine mapping.
void set_val_inf(double val)
Sets the value of the Cmp to val at infinity.
const Valeur & mult_st() const
Returns applied to *this.
double rayon
Radius of the horizon in LORENE's units.
Cmp pow(const Cmp &, int)
Power .
void std_base_scal()
Sets the spectral bases of the Valeur va to the standard ones for a scalar.
void set_val_hor(double val, int zone)
Sets the value of the Cmp to val on the inner boudary of the shell number zone .This is usefull for d...
const Valeur & mult_cp() const
Returns applied to *this.
void raccord(int n)
Performs the matching of the nucleus with respect to the first shell.
double omega
Angular velocity in LORENE's units.
void set_dzpuis(int)
Set a value to dzpuis.
Tenseur n_auto
Part of N generated by the hole.
void set_etat_qcq()
Sets the logical state to ETATQCQ (ordinary state).
Valeur va
The numerical value of the Cmp.
Coord r
r coordinate centered on the grid