29char single_bound_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Isol_hor/single_bound.C,v 1.3 2014/10/13 08:53:01 j_novak Exp $" ;
57#include "time_slice.h"
63#include "utilitaires.h"
96 for (
int k=0 ; k<nnp ; k++)
97 for (
int j=0 ; j<nnt ; j++)
106const Valeur Single_hor::boundary_nn_Neu(
double cc)
const {
124 for (
int k=0 ; k<nnp ; k++)
125 for (
int j=0 ; j<nnt ; j++)
128 nn_bound.std_base_scal() ;
135const Valeur Single_hor::boundary_nn_Dir(
double cc)
const {
150 tmp = (tmp + rho *
nn)/(1 + rho) ;
161 for (
int k=0 ; k<nnp ; k++)
162 for (
int j=0 ; j<nnt ; j++)
165 nn_bound.std_base_scal() ;
174const Valeur Single_hor:: boundary_beta_x(
double om_orb,
175 double om_loc)
const {
179 assert ((orientation == 0) || (orientation == M_PI)) ;
180 int aligne = (orientation == 0) ? 1 : -1 ;
195 ya_mtbl.set_etat_qcq() ;
199 yy_mtbl.set_etat_qcq() ;
202 for (
int k=0 ; k<nnp ; k++)
203 for (
int j=0 ; j<nnt ; j++)
204 lim_x.
set(0, k, j, 0) = aligne * om_orb * ya_mtbl(1, k, j, 0)
205 + (om_loc-om_orb)* yy_mtbl(1, k, j, 0)
206 + tmp_vect(1).val_grid_point(1, k, j, 0) ;
217const Valeur Single_hor:: boundary_beta_y(
double om_orb,
218 double om_loc)
const {
222 assert ((orientation == 0) || (orientation == M_PI)) ;
223 int aligne = (orientation == 0) ? 1 : -1 ;
239 xa_mtbl.set_etat_qcq() ;
243 xx_mtbl.set_etat_qcq() ;
246 for (
int k=0 ; k<nnp ; k++)
247 for (
int j=0 ; j<nnt ; j++)
248 lim_y.
set(0, k, j, 0) = - aligne *om_orb * xa_mtbl(1, k, j, 0) -
249 (om_loc-om_orb)*xx_mtbl(1, k, j, 0)
250 + tmp_vect(2).val_grid_point(1, k, j, 0) ;
260const Valeur Single_hor:: boundary_beta_z()
const {
274 for (
int k=0 ; k<nnp ; k++)
275 for (
int j=0 ; j<nnt ; j++)
276 lim_z.set(0, k, j, 0) = tmp_vect(3).val_grid_point(1, k, j, 0) ;
void set(const Map *mp, Mtbl *(*construct)(const Map *))
Semi-constructor from a mapping and a method.
const Base_vect_cart & get_bvect_cart() const
Returns the Cartesian basis associated with the coordinates (x,y,z) of the mapping,...
Coord y
y coordinate centered on the grid
Coord ya
Absolute y coordinate.
const Base_vect_spher & get_bvect_spher() const
Returns the orthonormal vectorial basis associated with the coordinates of the mapping.
Coord x
x coordinate centered on the grid
double get_rot_phi() const
Returns the angle between the x –axis and X –axis.
Coord xa
Absolute x coordinate.
const Mg3d * get_mg() const
Gives the Mg3d on which the mapping is defined.
const Metric_flat & flat_met_spher() const
Returns the flat metric associated with the spherical coordinates and with components expressed in th...
Metric for tensor calculation.
virtual const Vector & radial_vect() const
Returns the radial vector normal to a spherical slicing and pointing toward spatial infinity.
const Mg3d * get_angu() const
Returns the pointer on the associated angular grid.
int get_np(int l) const
Returns the number of points in the azimuthal direction ( ) in domain no. l.
int get_nt(int l) const
Returns the number of points in the co-latitude direction ( ) in domain no. l.
Base_val ** std_base_vect_cart() const
Returns the standard spectral bases for the Cartesian components of a vector.
Tensor field of valence 0 (or component of a tensorial field).
const Vector & derive_cov(const Metric &gam) const
Returns the gradient (1-form = covariant vector) of *this
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 Metric & get_gam() const
metric
const Valeur boundary_psi_app_hor() const
Neumann boundary condition for.
const Sym_tensor & get_k_dd() const
k_dd
Metric_flat ff
3 metric flat
Scalar psi
Conformal factor .
Map_af & mp
Affine mapping.
Vector dn
Covariant derivative of the lapse with respect to the flat metric .
Scalar nn
Lapse function .
Scalar trK
Trace of the extrinsic curvature.
Metric tgam
3 metric tilde
Values and coefficients of a (real-value) function.
Tbl & set(int l)
Read/write of the value in a given domain (configuration space).
void std_base_scal()
Sets the bases for spectral expansions (member base ) to the standard ones for a scalar.
Tensor field of valence 1.
virtual void std_spectral_base()
Sets the standard spectal bases of decomposition for each component.
const Scalar & divergence(const Metric &) const
The divergence of this with respect to a Metric .
virtual void change_triad(const Base_vect &)
Sets a new vectorial basis (triad) of decomposition and modifies the components accordingly.
Scalar & set(int)
Read/write access to a component.
Tenseur contract(const Tenseur &, int id1, int id2)
Self contraction of two indices of a Tenseur .