60 #include "time_slice.h" 63 #include "evolution.h" 65 #include "graphique.h" 66 #include "utilitaires.h" 99 for (
int k=0 ; k<nnp ; k++)
100 for (
int j=0 ; j<nnt ; j++)
103 psi_bound.std_base_scal() ;
109 const Valeur Single_hor::boundary_nn_Neu(
double cc)
const {
127 for (
int k=0 ; k<nnp ; k++)
128 for (
int j=0 ; j<nnt ; j++)
131 nn_bound.std_base_scal() ;
138 const Valeur Single_hor::boundary_nn_Dir(
double cc)
const {
153 tmp = (tmp + rho *
nn)/(1 + rho) ;
164 for (
int k=0 ; k<nnp ; k++)
165 for (
int j=0 ; j<nnt ; j++)
168 nn_bound.std_base_scal() ;
177 const Valeur Single_hor:: boundary_beta_x(
double om_orb,
178 double om_loc)
const {
182 assert ((orientation == 0) || (orientation == M_PI)) ;
183 int aligne = (orientation == 0) ? 1 : -1 ;
198 ya_mtbl.set_etat_qcq() ;
202 yy_mtbl.set_etat_qcq() ;
205 for (
int k=0 ; k<nnp ; k++)
206 for (
int j=0 ; j<nnt ; j++)
207 lim_x.
set(0, k, j, 0) = aligne * om_orb * ya_mtbl(1, k, j, 0)
208 + (om_loc-om_orb)* yy_mtbl(1, k, j, 0)
209 + tmp_vect(1).val_grid_point(1, k, j, 0) ;
220 const Valeur Single_hor:: boundary_beta_y(
double om_orb,
221 double om_loc)
const {
225 assert ((orientation == 0) || (orientation == M_PI)) ;
226 int aligne = (orientation == 0) ? 1 : -1 ;
242 xa_mtbl.set_etat_qcq() ;
246 xx_mtbl.set_etat_qcq() ;
249 for (
int k=0 ; k<nnp ; k++)
250 for (
int j=0 ; j<nnt ; j++)
251 lim_y.
set(0, k, j, 0) = - aligne *om_orb * xa_mtbl(1, k, j, 0) -
252 (om_loc-om_orb)*xx_mtbl(1, k, j, 0)
253 + tmp_vect(2).val_grid_point(1, k, j, 0) ;
263 const Valeur Single_hor:: boundary_beta_z()
const {
277 for (
int k=0 ; k<nnp ; k++)
278 for (
int j=0 ; j<nnt ; j++)
279 lim_z.set(0, k, j, 0) = tmp_vect(3).val_grid_point(1, k, j, 0) ;
Coord xa
Absolute x coordinate.
Base_val ** std_base_vect_cart() const
Returns the standard spectral bases for the Cartesian components of a vector.
Metric for tensor calculation.
const Metric & get_gam() const
metric
Vector dn
Covariant derivative of the lapse with respect to the flat metric .
Metric tgam
3 metric tilde
int get_np(int l) const
Returns the number of points in the azimuthal direction ( ) in domain no. l.
const Sym_tensor & get_k_dd() const
k_dd
const Base_vect_spher & get_bvect_spher() const
Returns the orthonormal vectorial basis associated with the coordinates of the mapping.
const Mg3d * get_mg() const
Gives the Mg3d on which the mapping is defined.
Tensor field of valence 0 (or component of a tensorial field).
virtual const Vector & radial_vect() const
Returns the radial vector normal to a spherical slicing and pointing toward spatial infinity...
virtual void std_spectral_base()
Sets the spectral bases of the Valeur va to the standard ones for a scalar field. ...
Values and coefficients of a (real-value) function.
virtual void change_triad(const Base_vect &)
Sets a new vectorial basis (triad) of decomposition and modifies the components accordingly.
double get_rot_phi() const
Returns the angle between the x –axis and X –axis.
void set(const Map *mp, Mtbl *(*construct)(const Map *))
Semi-constructor from a mapping and a method.
Tensor field of valence 1.
double val_grid_point(int l, int k, int j, int i) const
Returns the value of the field at a specified grid point.
const Valeur boundary_psi_app_hor() const
Neumann boundary condition for.
Scalar psi
Conformal factor .
Scalar trK
Trace of the extrinsic curvature.
const Mg3d * get_angu() const
Returns the pointer on the associated angular grid.
const Scalar & divergence(const Metric &) const
The divergence of this with respect to a Metric .
Tenseur contract(const Tenseur &, int id1, int id2)
Self contraction of two indices of a Tenseur .
Coord ya
Absolute y coordinate.
const Base_vect_cart & get_bvect_cart() const
Returns the Cartesian basis associated with the coordinates (x,y,z) of the mapping, i.e.
Map_af & mp
Affine mapping.
Scalar & set(const Itbl &ind)
Returns the value of a component (read/write version).
Coord y
y coordinate centered on the grid
Coord x
x coordinate centered on the grid
int get_nt(int l) const
Returns the number of points in the co-latitude direction ( ) in domain no. l.
Scalar & set(int)
Read/write access to a component.
Metric_flat ff
3 metric flat
const Vector & derive_cov(const Metric &gam) const
Returns the gradient (1-form = covariant vector) of *this
const Metric_flat & flat_met_spher() const
Returns the flat metric associated with the spherical coordinates and with components expressed in th...
Scalar nn
Lapse function .
Tbl & set(int l)
Read/write of the value in a given domain (configuration space).