23char bhole_solve_phi_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Bhole_binaire/bhole_solve_phi.C,v 1.4 2014/10/13 08:52:40 j_novak Exp $" ;
69#include "utilitaires.h"
110 assert ((relax>0) && (relax<=1)) ;
112 cout <<
"-----------------------------------------------" << endl ;
113 cout <<
"Resolution PSI" << endl ;
137 neumann_binaire (source_un, source_deux, lim_un, lim_deux,
Bhole hole1
Black hole one.
void set_omega(double ome)
Sets the orbital velocity to ome.
void solve_phi(double precision, double relax)
Solve the equation for the logarithm of .
Bhole hole2
Black hole two.
void init_phi()
Initiates the system for a resolution using the logarithm of .
Tenseur psi_auto
Part of generated by the hole.
Tenseur shift_auto
Part of generated by the hole.
Tenseur psi_comp
Part of generated by the companion hole.
Tenseur_sym tkij_auto
Auto .
void init_bhole_phi()
Initiates the black hole for a resolution with .
Tenseur grad_psi_tot
Total gradient of .
Tenseur_sym taij_auto
Part of generated by the hole.
Tenseur grad_n_tot
Total gradient of N .
Tenseur n_auto
Part of N generated by the hole.
Tenseur_sym taij_comp
Part of generated by the companion hole.
double rayon
Radius of the horizon in LORENE's units.
Tenseur n_comp
Part of N generated by the companion hole.
Map_af & mp
Affine mapping.
Tenseur_sym taij_tot
Total , which must be zero on the horizon of the regularisation on the shift has been done.
void fait_psi_comp(const Bhole &comp)
Imports the part of due to the companion hole comp .
Tenseur_sym tkij_tot
Total .
Cmp decouple
Function used to construct the part of generated by the hole from the total .
Component of a tensorial field *** DEPRECATED : use class Scalar instead ***.
void std_base_scal()
Sets the spectral bases of the Valeur va to the standard ones for a scalar.
void annule(int l)
Sets the Cmp to zero in a given domain.
void set_etat_zero()
Sets the logical state to ETATZERO (zero).
void set_dzpuis(int)
Set a value to dzpuis.
void raccord(int n)
Performs the matching of the nucleus with respect to the first shell.
Coord r
r 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.
Tensor handling *** DEPRECATED : use class Tensor instead ***.
Cmp & set()
Read/write for a scalar (see also operator=(const Cmp&) ).
const Tenseur & gradient() const
Returns the gradient of *this (Cartesian coordinates)
void set_std_base()
Set the standard spectal basis of decomposition for each component.
void set_etat_zero()
Sets the logical state to ETATZERO (zero state).
Values and coefficients of a (real-value) function.
void std_base_scal()
Sets the bases for spectral expansions (member base ) to the standard ones for a scalar.
Tbl diffrelmax(const Cmp &a, const Cmp &b)
Relative difference between two Cmp (max version).
Cmp log(const Cmp &)
Neperian logarithm.
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...