67 if ( (
mp.
get_mg())->get_np(0) == 1) {
79 assert( dynamic_cast<const Map_radial*>(&
mp) != 0x0 ) ;
82 dnpdr.
va = (dnpdr.
va).mult_st() ;
86 dnpdt.
va = (dnpdt.
va).mult_ct() ;
89 Cmp tmp = dnpdr + dnpdt ;
101 dnpdr =
nphi().dsdr() ;
102 dnpdt =
nphi().srdsdt() ;
105 dnpdr.
va = (dnpdr.
va).mult_ct() ;
108 if (dnpdt.get_etat() == ETATQCQ) {
109 dnpdt.va = (dnpdt.va).mult_st() ;
112 tmp = dnpdr - dnpdt ;
154 for (
int i=0; i<3; i++) {
155 for (
int j=i; j<3; j++) {
156 tkij.
set(i, j) = dn(i, j) + dn(j, i) ;
158 tkij.
set(i, i) -= double(2) /double(3) * divn() ;
165 assert( divn.
get_etat() == ETATZERO ) ;
181 for (
int i=0; i<3; i++) {
182 for (
int j=0; j<3; j++) {
Component of a tensorial field *** DEPRECATED : use class Scalar instead ***.
const Mg3d * get_mg() const
Gives the Mg3d on which the mapping is defined.
int get_etat() const
Returns the logical state.
Tenseur nnn
Total lapse function.
Tenseur nphi
Metric coefficient .
Tenseur b_car
Square of the metric factor B.
Tenseur shift
Total shift vector.
void extrinsic_curvature()
Computes tkij and ak_car from shift , nnn and b_car .
Cmp & set()
Read/write for a scalar (see also operator=(const Cmp&) ).
void mult_rsint()
Multiplication by .
Map & mp
Mapping associated with the star.
int get_etat() const
Returns the logical state.
Tenseur contract(const Tenseur &, int id1, int id2)
Self contraction of two indices of a Tenseur .
void set_etat_qcq()
Sets the logical state to ETATQCQ (ordinary state).
void set_etat_zero()
Sets the logical state to ETATZERO (zero state).
Valeur va
The numerical value of the Cmp.
Tenseur_sym tkij
Tensor related to the extrinsic curvature tensor by .
Tensor handling *** DEPRECATED : use class Tensor instead ***.
const Tenseur & gradient() const
Returns the gradient of *this (Cartesian coordinates)