29char scalar_import_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Tensor/Scalar/scalar_import.C,v 1.5 2014/10/13 08:53:46 j_novak Exp $" ;
70 int nz =
mp->get_mg()->get_nzone() ;
96 * (
mp_d->get_bvect_cart()).get_align() ;
116 cout <<
"Scalar::import : unexpected value of align_rel : "
146 if (
cm_d.get_etat() == ETATZERO) {
151 if (
cm_d.get_etat() == ETATUN) {
161 if (
cm_d.get_dzpuis() != 0) {
163 "Scalar::import : the dzpuis of the Scalar to be imported must be zero !"
207 if ( (
mp->r).c == 0x0 ) (
mp->r).fait() ;
208 if ( (
mp->tet).c == 0x0 ) (
mp->tet).fait() ;
209 if ( (
mp->phi).c == 0x0 ) (
mp->phi).fait() ;
210 if ( (
mp->xa).c == 0x0 ) (
mp->xa).fait() ;
211 if ( (
mp->ya).c == 0x0 ) (
mp->ya).fait() ;
212 if ( (
mp->za).c == 0x0 ) (
mp->za).fait() ;
224 double precis = 1
e-15 ;
233 for (
int l=0;
l < nzet;
l++) {
235 int nr =
mg_a->get_nr(
l) ;
236 int nt =
mg_a->get_nt(
l) ;
237 int np =
mg_a->get_np(
l) ;
254 for (
int k=0 ;
k<np ;
k++) {
255 for (
int j=0 ;
j<nt ;
j++) {
256 for (
int i=0 ;
i<nr ;
i++) {
340 if (
cm_d.get_etat() == ETATZERO) {
345 if (
cm_d.get_etat() == ETATUN) {
354 int align = (
mp->get_bvect_cart()).get_align() ;
356 assert( align * (
mp_d->get_bvect_cart()).get_align() == -1 ) ;
360 if (
cm_d.get_dzpuis() != 0) {
362 "Scalar::import : the dzpuis of the Scalar to be imported must be zero !"
405 if ( (
mp->r).c == 0x0 ) (
mp->r).fait() ;
406 if ( (
mp->tet).c == 0x0 ) (
mp->tet).fait() ;
407 if ( (
mp->phi).c == 0x0 ) (
mp->phi).fait() ;
408 if ( (
mp->x).c == 0x0 ) (
mp->x).fait() ;
409 if ( (
mp->y).c == 0x0 ) (
mp->y).fait() ;
410 if ( (
mp->z).c == 0x0 ) (
mp->z).fait() ;
422 double precis = 1
e-15 ;
431 for (
int l=0;
l < nzet;
l++) {
433 int nr =
mg_a->get_nr(
l) ;
434 int nt =
mg_a->get_nt(
l) ;
435 int np =
mg_a->get_np(
l) ;
452 for (
int k=0 ;
k<np ;
k++) {
453 for (
int j=0 ;
j<nt ;
j++) {
454 for (
int i=0 ;
i<nr ;
i++) {
537 if (
cm_d.get_etat() == ETATZERO) {
541 if (
cm_d.get_etat() == ETATUN) {
550 int align = (
mp->get_bvect_cart()).get_align() ;
552 assert( align * (
mp_d->get_bvect_cart()).get_align() == 1 ) ;
556 if (
cm_d.get_dzpuis() != 0) {
558 "Scalar::import : the dzpuis of the Scalar to be imported must be zero !"
601 if ( (
mp->r).c == 0x0 ) (
mp->r).fait() ;
602 if ( (
mp->tet).c == 0x0 ) (
mp->tet).fait() ;
603 if ( (
mp->phi).c == 0x0 ) (
mp->phi).fait() ;
604 if ( (
mp->x).c == 0x0 ) (
mp->x).fait() ;
605 if ( (
mp->y).c == 0x0 ) (
mp->y).fait() ;
606 if ( (
mp->z).c == 0x0 ) (
mp->z).fait() ;
618 double precis = 1
e-15 ;
627 for (
int l=0;
l < nzet;
l++) {
629 int nr =
mg_a->get_nr(
l) ;
630 int nt =
mg_a->get_nt(
l) ;
631 int np =
mg_a->get_np(
l) ;
648 for (
int k=0 ;
k<np ;
k++) {
649 for (
int j=0 ;
j<nt ;
j++) {
650 for (
int i=0 ;
i<nr ;
i++) {
Time evolution with partial storage (*** under development ***).
Base class for coordinate mappings.
Tbl ** t
Array (size nzone ) of pointers on the Tbl 's.
void set_etat_qcq()
Sets the logical state to ETATQCQ (ordinary state).
Tensor field of valence 0 (or component of a tensorial field).
void set_etat_one()
Sets the logical state to ETATUN (one).
void import_align(int nzet, const Scalar &ci)
Assignment to another Scalar defined on a different mapping, when the two mappings have aligned Carte...
virtual void set_etat_qcq()
Sets the logical state to ETATQCQ (ordinary state).
void import_gal(int nzet, const Scalar &ci)
Assignment to another Scalar defined on a different mapping, when the two mappings do not have a part...
virtual void set_etat_zero()
Sets the logical state to ETATZERO (zero).
virtual void annule(int l_min, int l_max)
Sets the Scalar to zero in several domains.
void del_t()
Logical destructor.
void set_dzpuis(int)
Modifies the dzpuis flag.
friend Scalar sqrt(const Scalar &)
Square root.
Valeur va
The numerical value of the Scalar
void import_anti(int nzet, const Scalar &ci)
Assignment to another Scalar defined on a different mapping, when the two mappings have anti-aligned ...
void import(const Scalar &ci)
Assignment to another Scalar defined on a different mapping.
Values and coefficients of a (real-value) function.
void set_etat_c_qcq()
Sets the logical state to ETATQCQ (ordinary state) for values in the configuration space (Mtbl c ).
Mtbl * c
Values of the function at the points of the multi-grid
const Map *const mp
Mapping on which the numerical values at the grid points are defined.