24char map_af_fait_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Map/map_af_fait.C,v 1.13 2014/10/13 08:53:02 j_novak Exp $" ;
116Mtbl* map_af_fait_r(
const Map* cvi) {
119 const Map_af* cv =
static_cast<const Map_af*
>(cvi) ;
120 const Mg3d* mg = cv->
get_mg() ;
124 Mtbl* mti =
new Mtbl(mg) ;
125 mti->set_etat_qcq() ;
128 double* alpha = cv->alpha ;
129 double* beta = cv->beta ;
132 for (
int l=0 ; l<nz ; l++) {
133 int ir = mg->get_nr(l);
134 int it = mg->get_nt(l) ;
135 int ip = mg->get_np(l) ;
136 const Grille3d* g = mg->get_grille3d(l) ;
137 Tbl* tb = (mti->t)[l] ;
139 double* p_r = tb->t ;
141 switch(mg->get_type_r(l)) {
143 for (k=0 ; k<ip ; k++) {
144 for (j=0 ; j<it ; j++) {
145 for (i=0 ; i<ir ; i++) {
146 *p_r = alpha[l] * (g->x)[i] + beta[l] ;
154 for (k=0 ; k<ip ; k++) {
155 for (j=0 ; j<it ; j++) {
156 for (i=0 ; i<ir ; i++) {
157 *p_r = 1./(alpha[l] * (g->x)[i] + beta[l]) ;
165 cout <<
"Map_af_fait_r: unknown type_r !\n" ;
180Mtbl* map_af_fait_tet(
const Map* cvi) {
183 const Map_af* cv =
static_cast<const Map_af*
>(cvi) ;
184 const Mg3d* mg = cv->
get_mg() ;
188 Mtbl* mti =
new Mtbl(mg) ;
189 mti->set_etat_qcq() ;
192 for (
int l=0 ; l<nz ; l++) {
193 int ir = mg->get_nr(l);
194 int it = mg->get_nt(l);
195 int ip = mg->get_np(l);
196 const Grille3d* g = mg->get_grille3d(l) ;
197 Tbl* tb = (mti->t)[l] ;
199 double* p_r = tb->t ;
200 for (k=0 ; k<ip ; k++) {
201 for (j=0 ; j<it ; j++) {
202 for (i=0 ; i<ir ; i++) {
218Mtbl* map_af_fait_phi(
const Map* cvi) {
221 const Map_af* cv =
static_cast<const Map_af*
>(cvi) ;
222 const Mg3d* mg = cv->
get_mg() ;
226 Mtbl* mti =
new Mtbl(mg) ;
227 mti->set_etat_qcq() ;
230 for (
int l=0 ; l<nz ; l++) {
231 int ir = mg->get_nr(l);
232 int it = mg->get_nt(l);
233 int ip = mg->get_np(l);
234 const Grille3d* g = mg->get_grille3d(l) ;
235 Tbl* tb = (mti->t)[l] ;
237 double* p_r = tb->t ;
238 for (k=0 ; k<ip ; k++) {
239 for (j=0 ; j<it ; j++) {
240 for (i=0 ; i<ir ; i++) {
256Mtbl* map_af_fait_x(
const Map* cvi) {
259 const Mg3d* mg = cvi->
get_mg() ;
262 Mtbl* mti =
new Mtbl(mg) ;
264 *mti = (cvi->r) * (cvi->sint) * (cvi->cosp) ;
274Mtbl* map_af_fait_y(
const Map* cvi) {
277 const Mg3d* mg = cvi->
get_mg() ;
280 Mtbl* mti =
new Mtbl(mg) ;
282 *mti = (cvi->r) * (cvi->sint) * (cvi->sinp) ;
292Mtbl* map_af_fait_z(
const Map* cvi) {
295 const Mg3d* mg = cvi->
get_mg() ;
298 Mtbl* mti =
new Mtbl(mg) ;
300 *mti = (cvi->r) * (cvi->cost) ;
310Mtbl* map_af_fait_xa(
const Map* cvi) {
313 const Mg3d* mg = cvi->
get_mg() ;
316 Mtbl* mti =
new Mtbl(mg) ;
318 double r_phi = cvi->get_rot_phi() ;
319 double t_x = cvi->get_ori_x() ;
321 if ( fabs(r_phi) < 1.e-14 ) {
322 *mti = (cvi->x) + t_x ;
324 else if ( fabs(r_phi - M_PI) < 1.e-14 ) {
325 *mti = - (cvi->x) + t_x ;
328 Mtbl phi = cvi->phi + r_phi ;
329 *mti = (cvi->r) * (cvi->sint) *
cos(phi) + t_x ;
340Mtbl* map_af_fait_ya(
const Map* cvi) {
343 const Mg3d* mg = cvi->
get_mg() ;
346 Mtbl* mti =
new Mtbl(mg) ;
348 double r_phi = cvi->get_rot_phi() ;
349 double t_y = cvi->get_ori_y() ;
351 if ( fabs(r_phi) < 1.e-14 ) {
352 *mti = (cvi->y) + t_y ;
354 else if ( fabs(r_phi - M_PI) < 1.e-14 ) {
355 *mti = - (cvi->y) + t_y ;
358 Mtbl phi = cvi->phi + r_phi ;
359 *mti = (cvi->r) * (cvi->sint) *
sin(phi) + t_y ;
370Mtbl* map_af_fait_za(
const Map* cvi) {
373 const Mg3d* mg = cvi->
get_mg() ;
375 double t_z = cvi->get_ori_z() ;
378 Mtbl* mti =
new Mtbl(mg) ;
380 *mti = (cvi->r) * (cvi->cost) + t_z ;
390Mtbl* map_af_fait_sint(
const Map* cvi) {
393 const Map_af* cv =
static_cast<const Map_af*
>(cvi) ;
394 const Mg3d* mg = cv->
get_mg() ;
398 Mtbl* mti =
new Mtbl(mg) ;
399 mti->set_etat_qcq() ;
402 for (
int l=0 ; l<nz ; l++) {
403 int ir = mg->get_nr(l);
404 int it = mg->get_nt(l);
405 int ip = mg->get_np(l);
406 const Grille3d* g = mg->get_grille3d(l) ;
407 Tbl* tb = (mti->t)[l] ;
409 double* p_r = tb->t ;
410 for (k=0 ; k<ip ; k++) {
411 for (j=0 ; j<it ; j++) {
412 for (i=0 ; i<ir ; i++) {
413 *p_r =
sin(g->tet[j]) ;
424Mtbl* map_af_fait_cost(
const Map* cvi) {
427 const Map_af* cv =
static_cast<const Map_af*
>(cvi) ;
428 const Mg3d* mg = cv->
get_mg() ;
432 Mtbl* mti =
new Mtbl(mg) ;
433 mti->set_etat_qcq() ;
436 for (
int l=0 ; l<nz ; l++) {
437 int ir = mg->get_nr(l);
438 int it = mg->get_nt(l);
439 int ip = mg->get_np(l);
440 const Grille3d* g = mg->get_grille3d(l) ;
441 Tbl* tb = (mti->t)[l] ;
443 double* p_r = tb->t ;
444 for (k=0 ; k<ip ; k++) {
445 for (j=0 ; j<it ; j++) {
446 for (i=0 ; i<ir ; i++) {
447 *p_r =
cos(g->tet[j]) ;
458Mtbl* map_af_fait_sinp(
const Map* cvi) {
461 const Map_af* cv =
static_cast<const Map_af*
>(cvi) ;
462 const Mg3d* mg = cv->
get_mg() ;
466 Mtbl* mti =
new Mtbl(mg) ;
467 mti->set_etat_qcq() ;
470 for (
int l=0 ; l<nz ; l++) {
471 int ir = mg->get_nr(l);
472 int it = mg->get_nt(l);
473 int ip = mg->get_np(l);
474 const Grille3d* g = mg->get_grille3d(l) ;
475 Tbl* tb = (mti->t)[l] ;
477 double* p_r = tb->t ;
478 for (k=0 ; k<ip ; k++) {
479 for (j=0 ; j<it ; j++) {
480 for (i=0 ; i<ir ; i++) {
481 *p_r =
sin(g->phi[k]) ;
492Mtbl* map_af_fait_cosp(
const Map* cvi) {
495 const Map_af* cv =
static_cast<const Map_af*
>(cvi) ;
496 const Mg3d* mg = cv->
get_mg() ;
500 Mtbl* mti =
new Mtbl(mg) ;
501 mti->set_etat_qcq() ;
504 for (
int l=0 ; l<nz ; l++) {
505 int ir = mg->get_nr(l);
506 int it = mg->get_nt(l);
507 int ip = mg->get_np(l);
508 const Grille3d* g = mg->get_grille3d(l) ;
509 Tbl* tb = (mti->t)[l] ;
511 double* p_r = tb->t ;
512 for (k=0 ; k<ip ; k++) {
513 for (j=0 ; j<it ; j++) {
514 for (i=0 ; i<ir ; i++) {
515 *p_r =
cos(g->phi[k]) ;
532Mtbl* map_af_fait_xsr(
const Map* cvi) {
535 const Map_af* cv =
static_cast<const Map_af*
>(cvi) ;
536 const Mg3d* mg = cv->
get_mg() ;
540 Mtbl* mti =
new Mtbl(mg) ;
541 mti->set_etat_qcq() ;
544 double* alpha = cv->alpha ;
545 double* beta = cv->beta ;
548 for (
int l=0 ; l<nz ; l++) {
549 int ir = mg->get_nr(l);
550 int it = mg->get_nt(l) ;
551 int ip = mg->get_np(l) ;
552 const Grille3d* g = mg->get_grille3d(l) ;
553 Tbl* tb = (mti->t)[l] ;
555 double* p_r = tb->t ;
557 switch(mg->get_type_r(l)) {
561 for (k=0 ; k<ip ; k++) {
562 for (j=0 ; j<it ; j++) {
563 for (i=0 ; i<ir ; i++) {
564 *p_r = 1. / alpha[l] ;
572 for (k=0 ; k<ip ; k++) {
573 for (j=0 ; j<it ; j++) {
575 *p_r = 1. / beta[l] ;
579 for (i=0 ; i<ir ; i++) {
580 *p_r = 1. / ( alpha[l] * (g->x)[i] + beta[l] ) ;
588 assert(beta[l] == - alpha[l]) ;
589 for (k=0 ; k<ip ; k++) {
590 for (j=0 ; j<it ; j++) {
591 for (i=0 ; i<ir ; i++) {
592 *p_r = 1. / alpha[l] ;
600 cout <<
"map_af_fait_xsr: unknown type_r !" << endl ;
617Mtbl* map_af_fait_dxdr(
const Map* cvi) {
620 const Map_af* cv =
static_cast<const Map_af*
>(cvi) ;
621 const Mg3d* mg = cv->
get_mg() ;
625 Mtbl* mti =
new Mtbl(mg) ;
626 mti->set_etat_qcq() ;
629 double* alpha = cv->alpha ;
632 for (
int l=0 ; l<nz ; l++) {
633 int ir = mg->get_nr(l);
634 int it = mg->get_nt(l) ;
635 int ip = mg->get_np(l) ;
636 Tbl* tb = (mti->t)[l] ;
638 double* p_r = tb->t ;
640 switch(mg->get_type_r(l)) {
643 for (k=0 ; k<ip ; k++) {
644 for (j=0 ; j<it ; j++) {
645 for (i=0 ; i<ir ; i++) {
646 *p_r = 1. / alpha[l] ;
654 for (k=0 ; k<ip ; k++) {
655 for (j=0 ; j<it ; j++) {
656 for (i=0 ; i<ir ; i++) {
657 *p_r = - 1. / alpha[l] ;
666 cout <<
"map_af_fait_dxdr: unknown type_r !" << endl ;
681Mtbl* map_af_fait_drdt(
const Map* cvi) {
684 const Mg3d* mg = cvi->
get_mg() ;
687 Mtbl* mti =
new Mtbl(mg) ;
688 mti->set_etat_zero() ;
699Mtbl* map_af_fait_stdrdp(
const Map* cvi) {
702 const Mg3d* mg = cvi->
get_mg() ;
705 Mtbl* mti =
new Mtbl(mg) ;
706 mti->set_etat_zero() ;
717Mtbl* map_af_fait_srdrdt(
const Map* cvi) {
720 const Mg3d* mg = cvi->
get_mg() ;
723 Mtbl* mti =
new Mtbl(mg) ;
724 mti->set_etat_zero() ;
735Mtbl* map_af_fait_srstdrdp(
const Map* cvi) {
738 const Mg3d* mg = cvi->
get_mg() ;
741 Mtbl* mti =
new Mtbl(mg) ;
742 mti->set_etat_zero() ;
753Mtbl* map_af_fait_sr2drdt(
const Map* cvi) {
756 const Mg3d* mg = cvi->
get_mg() ;
759 Mtbl* mti =
new Mtbl(mg) ;
760 mti->set_etat_zero() ;
771Mtbl* map_af_fait_sr2stdrdp(
const Map* cvi) {
774 const Mg3d* mg = cvi->
get_mg() ;
777 Mtbl* mti =
new Mtbl(mg) ;
778 mti->set_etat_zero() ;
789Mtbl* map_af_fait_d2rdx2(
const Map* cvi) {
792 const Mg3d* mg = cvi->
get_mg() ;
795 Mtbl* mti =
new Mtbl(mg) ;
796 mti->set_etat_zero() ;
807Mtbl* map_af_fait_lapr_tp(
const Map* cvi) {
810 const Mg3d* mg = cvi->
get_mg() ;
813 Mtbl* mti =
new Mtbl(mg) ;
814 mti->set_etat_zero() ;
825Mtbl* map_af_fait_d2rdtdx(
const Map* cvi) {
828 const Mg3d* mg = cvi->
get_mg() ;
831 Mtbl* mti =
new Mtbl(mg) ;
832 mti->set_etat_zero() ;
843Mtbl* map_af_fait_sstd2rdpdx(
const Map* cvi) {
846 const Mg3d* mg = cvi->
get_mg() ;
849 Mtbl* mti =
new Mtbl(mg) ;
850 mti->set_etat_zero() ;
861Mtbl* map_af_fait_sr2d2rdt2(
const Map* cvi) {
864 const Mg3d* mg = cvi->
get_mg() ;
867 Mtbl* mti =
new Mtbl(mg) ;
868 mti->set_etat_zero() ;
int get_nzone() const
Returns the number of domains.
const Mg3d * get_mg() const
Gives the Mg3d on which the Mtbl is defined.
Cmp sin(const Cmp &)
Sine.
Cmp cos(const Cmp &)
Cosine.