72 #include "graphique.h" 73 #include "utilitaires.h" 77 Cmp raccord_c1(
const Cmp& uu,
int l1) ;
107 for (
int i=1; i<=3; i++) {
108 v_orb.set(i) = www(i).val_grid_point(0, 0, 0, 0) ;
112 v_orb.std_spectral_base() ;
122 for (
int l=
nzet; l <= nzm1; l++) {
190 Cmp source_cmp (source) ;
191 Cmp zeta_h_cmp (zeta_h) ;
196 Cmp bb_cmp1 (bb(1)) ;
197 Cmp bb_cmp2 (bb(2)) ;
198 Cmp bb_cmp3 (bb(3)) ;
199 bb_cmp.set(0) = bb_cmp1 ;
200 bb_cmp.
set(1) = bb_cmp2 ;
201 bb_cmp.
set(2) = bb_cmp3 ;
203 source_cmp.
va.
ylm() ;
205 cout <<
"source" << endl <<
norme(source_cmp)<< endl ;
206 cout <<
"zeta_h " << endl <<
norme(zeta_h_cmp) << endl ;
207 cout <<
"bb(1)" << endl <<
norme(bb_cmp(0)) << endl ;
208 cout <<
"bb(2)" << endl <<
norme(bb_cmp(1)) << endl ;
209 cout <<
"bb(3)" << endl <<
norme(bb_cmp(2)) << endl ;
210 cout <<
"psiO" << endl <<
norme(psi0_cmp) << endl ;
216 cout <<
"psiO apres" << endl <<
norme(
psi0) << endl ;
228 double erreur =
diffrel(oper, source)(0) ;
230 cout <<
"Check of the resolution of the continuity equation : " 232 cout <<
" norme(source) : " <<
norme(source)(0)
233 <<
" diff oper/source : " << erreur << endl ;
243 for (
int i=1; i<=3; i++)
255 for (
int i=1; i<=3; i++) {
257 d_psi_i = raccord_c1(d_psi_i,
nzet) ;
virtual void set_etat_qcq()
Sets the logical state of all components to ETATQCQ (ordinary state).
Scalar psi0
Scalar potential of the non-translational part of fluid 4-velocity (in the irrotational case) ...
Component of a tensorial field *** DEPRECATED : use class Scalar instead ***.
Cmp exp(const Cmp &)
Exponential.
void add_int(const int &n, int position=0)
Adds the address of a new int to the list.
Map & mp
Mapping associated with the star.
void ylm_i()
Inverse of ylm()
Cmp der_nbar_ent(const Cmp &ent, int nzet, int l_min=0, Param *par=0x0) const
Computes the logarithmic derivative from the log-enthalpy and extra parameters.
virtual void annule(int l_min, int l_max)
Sets the Scalar to zero in several domains.
const Base_vect_spher & get_bvect_spher() const
Returns the orthonormal vectorial basis associated with the coordinates of the mapping.
Vector d_psi
Gradient of (in the irrotational case) (Spherical components with respect to the mapping of the star...
const Eos & eos
Equation of state of the stellar matter.
void ylm()
Computes the coefficients of *this.
const Mg3d * get_mg() const
Gives the Mg3d on which the mapping is defined.
const Scalar & laplacian(int ced_mult_r=4) const
Returns the Laplacian of *this.
Flat metric for tensor calculation.
Tensor field of valence 0 (or component of a tensorial field).
virtual void poisson_compact(const Cmp &source, const Cmp &aa, const Tenseur &bb, const Param &par, Cmp &psi) const =0
Resolution of the elliptic equation in the case where the stellar interior is covered by a single do...
virtual void std_spectral_base()
Sets the spectral bases of the Valeur va to the standard ones for a scalar field. ...
virtual void change_triad(const Base_vect &)
Sets a new vectorial basis (triad) of decomposition and modifies the components accordingly.
int get_etat() const
Returns the logical state ETATNONDEF (undefined), ETATZERO (null) or ETATQCQ (ordinary).
const Vector & derive_con(const Metric &gam) const
Returns the "contravariant" derivative of *this with respect to some metric , by raising the index of...
Tensor field of valence 1.
Tbl & set_domain(int l)
Read/write of the value in a given domain.
Tbl diffrel(const Cmp &a, const Cmp &b)
Relative difference between two Cmp (norme version).
int nzet
Number of domains of *mp occupied by the star.
Scalar gam_euler
Lorentz factor between the fluid and Eulerian observers.
Tbl norme(const Cmp &)
Sums of the absolute values of all the values of the Cmp in each domain.
Sym_tensor hij
Total deviation of the inverse conformal metric from the inverse flat metric.
const Base_vect * get_triad() const
Returns the vectorial basis (triad) on which the components are defined.
Metric_flat flat
Flat metric defined on the mapping (Spherical components with respect to the mapping of the star) ...
int get_nzone() const
Returns the number of domains.
Tenseur contract(const Tenseur &, int id1, int id2)
Self contraction of two indices of a Tenseur .
const Tensor & derive_cov(const Metric &gam) const
Returns the covariant derivative of this with respect to some metric .
virtual void inc_dzpuis(int inc=1)
Increases by inc units the value of dzpuis and changes accordingly the values in the compactified ext...
double velocity_potential(int mermax, double precis, double relax)
Computes the non-translational part of the velocity scalar potential by solving the continuity equat...
Tbl & set(int l)
Read/write of the value in a given domain.
const Base_vect_cart & get_bvect_cart() const
Returns the Cartesian basis associated with the coordinates (x,y,z) of the mapping, i.e.
void add_double(const double &x, int position=0)
Adds the the address of a new double to the list.
void set_triad(const Base_vect &new_triad)
Assigns a new vectorial basis (triad) of decomposition.
Valeur & set_spectral_va()
Returns va (read/write version)
Scalar & set(int)
Read/write access to a component.
void set_etat_qcq()
Sets the logical state to ETATQCQ (ordinary state).
const Vector & derive_cov(const Metric &gam) const
Returns the gradient (1-form = covariant vector) of *this
virtual void annule(int l_min, int l_max)
Sets the Tensor to zero in several domains.
Vector bsn
3-vector shift, divided by N, of the rotating coordinates, .
Valeur va
The numerical value of the Cmp.
Scalar psi4
Conformal factor .
const Metric_flat & flat_met_spher() const
Returns the flat metric associated with the spherical coordinates and with components expressed in th...
Tensor handling *** DEPRECATED : use class Tensor instead ***.
void add_int_mod(int &n, int position=0)
Adds the address of a new modifiable int to the list.
Tensor down(int ind, const Metric &gam) const
Computes a new tensor by lowering an index of *this.