72 #include "utilitaires.h" 73 #include "graphique.h" 113 assert ((relax>0) && (relax<=1)) ;
115 cout <<
"-----------------------------------------------" << endl ;
116 cout <<
"Resolution PSI" << endl ;
133 lim_un.std_base_scal() ;
137 lim_deux.std_base_scal() ;
140 neumann_binaire (source_un, source_deux, lim_un, lim_deux,
Cmp log(const Cmp &)
Neperian logarithm.
void solve_phi(double precision, double relax)
Solve the equation for the logarithm of .
Component of a tensorial field *** DEPRECATED : use class Scalar instead ***.
Tenseur grad_n_tot
Total gradient of N .
void annule(int l)
Sets the Cmp to zero in a given domain.
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 flat_scalar_prod(const Tenseur &t1, const Tenseur &t2)
Scalar product of two Tenseur when the metric is : performs the contraction of the last index of t1 w...
Tenseur psi_auto
Part of generated by the hole.
Values and coefficients of a (real-value) function.
Tenseur n_comp
Part of N generated by the companion hole.
Tenseur_sym taij_tot
Total , which must be zero on the horizon of the regularisation on the shift has been done...
Tenseur_sym taij_auto
Part of generated by the hole.
Tenseur shift_auto
Part of generated by the hole.
Tenseur_sym taij_comp
Part of generated by the companion hole.
void set_etat_zero()
Sets the logical state to ETATZERO (zero).
Cmp & set()
Read/write for a scalar (see also operator=(const Cmp&) ).
Tenseur psi_comp
Part of generated by the companion hole.
void init_bhole_phi()
Initiates the black hole for a resolution with .
const Mg3d * get_angu() const
Returns the pointer on the associated angular grid.
Map_af & mp
Affine mapping.
Bhole hole1
Black hole one.
void init_phi()
Initiates the system for a resolution using the logarithm of .
double rayon
Radius of the horizon in LORENE's units.
Tenseur_sym tkij_auto
Auto .
void std_base_scal()
Sets the spectral bases of the Valeur va to the standard ones for a scalar.
Tenseur_sym tkij_tot
Total .
Tenseur grad_psi_tot
Total gradient of .
void set_omega(double ome)
Sets the orbital velocity to ome.
Cmp decouple
Function used to construct the part of generated by the hole from the total .
void raccord(int n)
Performs the matching of the nucleus with respect to the first shell.
void set_dzpuis(int)
Set a value to dzpuis.
Tenseur n_auto
Part of N generated by the hole.
void fait_psi_comp(const Bhole &comp)
Imports the part of due to the companion hole comp .
void set_etat_zero()
Sets the logical state to ETATZERO (zero state).
Tensor handling *** DEPRECATED : use class Tensor instead ***.
Tbl diffrelmax(const Cmp &a, const Cmp &b)
Relative difference between two Cmp (max version).
const Tenseur & gradient() const
Returns the gradient of *this (Cartesian coordinates)
Coord r
r coordinate centered on the grid
Bhole hole2
Black hole two.