30char scalar_import_asymy_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Tensor/Scalar/scalar_import_asymy.C,v 1.5 2014/10/13 08:53:46 j_novak Exp $" ;
74 int nz =
mp->get_mg()->get_nzone() ;
100 * (
mp_d->get_bvect_cart()).get_align() ;
115 cout <<
"Scalar::import_asymy : unexpected value of align_rel : "
136 if (
cm_d.get_etat() == ETATZERO) {
140 if (
cm_d.get_etat() == ETATUN) {
149 int align = (
mp->get_bvect_cart()).get_align() ;
151 assert( align * (
mp_d->get_bvect_cart()).get_align() == -1 ) ;
155 if (
cm_d.get_dzpuis() != 0) {
157 "Scalar::import_anti_asymy : the dzpuis of the Scalar to be imported"
158 <<
" must be zero !" <<
endl ;
203 if ( (
mp->r).c == 0x0 ) (
mp->r).fait() ;
204 if ( (
mp->tet).c == 0x0 ) (
mp->tet).fait() ;
205 if ( (
mp->phi).c == 0x0 ) (
mp->phi).fait() ;
206 if ( (
mp->x).c == 0x0 ) (
mp->x).fait() ;
207 if ( (
mp->y).c == 0x0 ) (
mp->y).fait() ;
208 if ( (
mp->z).c == 0x0 ) (
mp->z).fait() ;
220 double precis = 1
e-15 ;
229 for (
int l=0;
l < nzet;
l++) {
231 int nr =
mg_a->get_nr(
l) ;
232 int nt =
mg_a->get_nt(
l) ;
233 int np =
mg_a->get_np(
l) ;
266 for (
int k=1 ;
k<np/2 ;
k++) {
267 for (
int j=0 ;
j<nt ;
j++) {
268 for (
int i=0 ;
i<nr ;
i++) {
339 for (
int k=np/2+1 ;
k<np ;
k++) {
346 for (
int j=0 ;
j<nt ;
j++) {
347 for (
int i=0 ;
i<nr ;
i++) {
348 *
ptx = - (*ptx_symy) ;
383 if (
cm_d.get_etat() == ETATZERO) {
387 if (
cm_d.get_etat() == ETATUN) {
396 int align = (
mp->get_bvect_cart()).get_align() ;
398 assert( align * (
mp_d->get_bvect_cart()).get_align() == 1 ) ;
402 if (
cm_d.get_dzpuis() != 0) {
404 "Scalar::import_align_asymy : the dzpuis of the Scalar to be imported"
405 <<
" must be zero !" <<
endl ;
449 if ( (
mp->r).c == 0x0 ) (
mp->r).fait() ;
450 if ( (
mp->tet).c == 0x0 ) (
mp->tet).fait() ;
451 if ( (
mp->phi).c == 0x0 ) (
mp->phi).fait() ;
452 if ( (
mp->x).c == 0x0 ) (
mp->x).fait() ;
453 if ( (
mp->y).c == 0x0 ) (
mp->y).fait() ;
454 if ( (
mp->z).c == 0x0 ) (
mp->z).fait() ;
466 double precis = 1
e-15 ;
475 for (
int l=0;
l < nzet;
l++) {
477 int nr =
mg_a->get_nr(
l) ;
478 int nt =
mg_a->get_nt(
l) ;
479 int np =
mg_a->get_np(
l) ;
513 for (
int k=1 ;
k<np/2 ;
k++) {
514 for (
int j=0 ;
j<nt ;
j++) {
515 for (
int i=0 ;
i<nr ;
i++) {
586 for (
int k=np/2+1 ;
k<np ;
k++) {
593 for (
int j=0 ;
j<nt ;
j++) {
594 for (
int i=0 ;
i<nr ;
i++) {
595 *
ptx = - (*ptx_symy) ;
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_asymy(const Scalar &ci)
Assignment to another Scalar defined on a different mapping.
void import_anti_asymy(int nzet, const Scalar &ci)
Assignment to another Scalar defined on a different mapping, when the two mappings have anti-aligned ...
virtual void set_etat_qcq()
Sets the logical state to ETATQCQ (ordinary state).
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_align_asymy(int nzet, const Scalar &ci)
Assignment to another Scalar defined on a different mapping, when the two mappings have aligned Carte...
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.