30char star_rot_lambda_grv2_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Star/star_rot_lambda_grv2.C,v 1.4 2014/10/13 08:53:39 j_novak Exp $" ;
67 cout <<
"Star_rot::lambda_grv2: the mapping of sou_m does not"
68 <<
endl <<
" belong to the class Map_radial !" <<
endl ;
74 sou_q.check_dzpuis(4) ;
87 for (
int l=0 ;
l<nz ;
l++) {
115 cout <<
"Star_rot::lambda_grv2: unknown type_r ! " <<
endl ;
132 for (
int l=0;
l<nz;
l++) {
167 cout <<
"Star_rot::lambda_grv2: unknown type_r ! " <<
endl ;
181 if (
sou_m.get_etat() == ETATZERO ) {
185 assert(
sou_m.get_spectral_va().get_etat() == ETATQCQ) ;
186 sou_m.get_spectral_va().coef_i() ;
187 s_m = *(
sou_m.get_spectral_va().c) ;
191 if (
sou_q.get_etat() == ETATZERO ) {
195 assert(
sou_q.get_spectral_va().get_etat() == ETATQCQ) ;
196 sou_q.get_spectral_va().coef_i() ;
197 s_q = *(
sou_q.get_spectral_va().c) ;
214 int*
ndl =
new int[nz+4] ;
216 for (
int l=0;
l<nz;
l++) {
226 for (
int l=0;
l<nz ;
l++) {
236 double*
erre =
new double [nz*
ndr] ;
238 for (
int l=0;
l<nz;
l++) {
241 for (
int i=0;
i<
ndl[1+
l];
i++) {
252 int taille =
ndrtp*nz ;
254 double*
tsou_m =
new double[ taille ] ;
255 double*
tsou_q =
new double[ taille ] ;
258 for (
int i=0;
i<taille;
i++) {
266 for (
int l=0;
l<nz;
l++) {
267 for (
int k=0;
k<
np1;
k++) {
268 for (
int j=0;
j<nt;
j++) {
282 for (
int l=0;
l<nz;
l++) {
283 for (
int k=0;
k<
np1;
k++) {
284 for (
int j=0;
j<nt;
j++) {
315 if (
int_q !=
double(0) ) {
Time evolution with partial storage (*** under development ***).
Base class for pure radial mappings.
const Grille3d * get_grille3d(int l) const
Returns a pointer on the 3D mono-grid for domain no. l.
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.
int get_nzone() const
Returns the number of domains.
int get_nr(int l) const
Returns the number of points in the radial direction ( ) in domain no. l.
int get_type_r(int l) const
Returns the type of sampling in the radial direction in domain no.
Tensor field of valence 0 (or component of a tensorial field).
static double lambda_grv2(const Scalar &sou_m, const Scalar &sou_q)
Computes the coefficient which ensures that the GRV2 virial identity is satisfied.