24char op_d2dtet2_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Non_class_members/Operators/op_d2sdtet2.C,v 1.5 2014/10/13 08:53:24 j_novak Exp $" ;
82void _d2sdtet2_pas_prevu(Tbl* ,
int & b) {
83 cout <<
"Unknown theta basis in Mtbl_cf::d2sdt2() !" << endl ;
84 cout <<
" basis: " << hex << b << endl ;
90void _d2sdtet2_t_cos(Tbl* tb,
int &)
94 if (tb->get_etat() == ETATZERO) {
99 assert(tb->get_etat() == ETATQCQ) ;
102 int nr = (tb->dim).dim[0] ;
103 int nt = (tb->dim).dim[1] ;
104 int np = (tb->dim).dim[2] ;
108 static double* cx = 0 ;
109 static int nt_pre =0 ;
114 cx =
reinterpret_cast<double*
>(realloc(cx, nt *
sizeof(
double))) ;
115 for (
int i=0 ; i<nt ; i++) {
116 cx[i] = - double(i * i) ;
121 double* xo =
new double[(tb->dim).taille] ;
124 for (
int i=0; i<(tb->dim).taille; i++) {
134 for (
int j=0 ; j<nt ; j++) {
135 for (
int i=0 ; i<nr ; i++ ) {
136 *xco = cx[j] * (*xci) ;
147 int borne_phi = np + 1 ;
148 if (np == 1) borne_phi = 1 ;
150 for (
int k=2 ; k<borne_phi ; k++) {
151 for (
int j=0 ; j<nt ; j++) {
152 for (
int i=0 ; i<nr ; i++ ) {
153 *xco = cx[j] * (*xci) ;
170void _d2sdtet2_t_sin(Tbl* tb,
int &)
174 if (tb->get_etat() == ETATZERO) {
179 assert(tb->get_etat() == ETATQCQ) ;
182 int nr = (tb->dim).dim[0] ;
183 int nt = (tb->dim).dim[1] ;
184 int np = (tb->dim).dim[2] ;
188 static double* cx = 0 ;
189 static int nt_pre =0 ;
194 cx =
reinterpret_cast<double*
>(realloc(cx, nt *
sizeof(
double))) ;
195 for (
int i=0 ; i<nt ; i++) {
196 cx[i] = - double(i * i) ;
201 double* xo =
new double[(tb->dim).taille] ;
204 for (
int i=0; i<(tb->dim).taille; i++) {
213 int borne_phi = np + 1 ;
214 if (np == 1) borne_phi = 1 ;
216 for (
int k=0 ; k< borne_phi ; k++) {
217 for (
int j=0 ; j<nt ; j++) {
218 for (
int i=0 ; i<nr ; i++ ) {
219 *xco = cx[j] * (*xci) ;
236void _d2sdtet2_t_cos_p(Tbl* tb,
int &)
240 if (tb->get_etat() == ETATZERO) {
245 assert(tb->get_etat() == ETATQCQ) ;
248 int nr = (tb->dim).dim[0] ;
249 int nt = (tb->dim).dim[1] ;
250 int np = (tb->dim).dim[2] ;
254 static double* cx = 0 ;
255 static int nt_pre =0 ;
260 cx =
reinterpret_cast<double*
>(realloc(cx, nt *
sizeof(
double))) ;
261 for (
int i=0 ; i<nt ; i++) {
262 cx[i] = - (2*i) * (2*i) ;
267 double* xo =
new double[(tb->dim).taille] ;
270 for (
int i=0; i<(tb->dim).taille; i++) {
280 for (
int j=0 ; j<nt ; j++) {
281 for (
int i=0 ; i<nr ; i++ ) {
282 *xco = cx[j] * (*xci) ;
293 int borne_phi = np + 1 ;
294 if (np == 1) borne_phi = 1 ;
296 for (
int k=2 ; k<borne_phi ; k++) {
297 for (
int j=0 ; j<nt ; j++) {
298 for (
int i=0 ; i<nr ; i++ ) {
299 *xco = cx[j] * (*xci) ;
316void _d2sdtet2_t_sin_p(Tbl* tb,
int &)
320 if (tb->get_etat() == ETATZERO) {
325 assert(tb->get_etat() == ETATQCQ) ;
328 int nr = (tb->dim).dim[0] ;
329 int nt = (tb->dim).dim[1] ;
330 int np = (tb->dim).dim[2] ;
334 static double* cx = 0 ;
335 static int nt_pre =0 ;
340 cx =
reinterpret_cast<double*
>(realloc(cx, nt *
sizeof(
double))) ;
341 for (
int i=0 ; i<nt ; i++) {
342 cx[i] = - (2*i) * (2*i) ;
347 double* xo =
new double[(tb->dim).taille] ;
350 for (
int i=0; i<(tb->dim).taille; i++) {
359 int borne_phi = np + 1 ;
360 if (np == 1) borne_phi = 1 ;
362 for (
int k=0 ; k< borne_phi ; k++) {
363 for (
int j=0 ; j<nt ; j++) {
364 for (
int i=0 ; i<nr ; i++ ) {
365 *xco = cx[j] * (*xci) ;
382void _d2sdtet2_t_sin_i(Tbl* tb,
int &)
386 if (tb->get_etat() == ETATZERO) {
391 assert(tb->get_etat() == ETATQCQ) ;
394 int nr = (tb->dim).dim[0] ;
395 int nt = (tb->dim).dim[1] ;
396 int np = (tb->dim).dim[2] ;
400 static double* cx = 0 ;
401 static int nt_pre =0 ;
406 cx =
reinterpret_cast<double*
>(realloc(cx, nt *
sizeof(
double))) ;
407 for (
int i=0 ; i<nt ; i++) {
408 cx[i] = - (2*i+1) * (2*i+1) ;
413 double* xo =
new double[(tb->dim).taille] ;
416 for (
int i=0; i<(tb->dim).taille; i++) {
425 int borne_phi = np + 1 ;
426 if (np == 1) borne_phi = 1 ;
428 for (
int k=0 ; k< borne_phi ; k++) {
429 for (
int j=0 ; j<nt ; j++) {
430 for (
int i=0 ; i<nr ; i++ ) {
431 *xco = cx[j] * (*xci) ;
448void _d2sdtet2_t_cos_i(Tbl* tb,
int &)
452 if (tb->get_etat() == ETATZERO) {
457 assert(tb->get_etat() == ETATQCQ) ;
460 int nr = (tb->dim).dim[0] ;
461 int nt = (tb->dim).dim[1] ;
462 int np = (tb->dim).dim[2] ;
466 static double* cx = 0 ;
467 static int nt_pre =0 ;
472 cx =
reinterpret_cast<double*
>(realloc(cx, nt *
sizeof(
double))) ;
473 for (
int i=0 ; i<nt ; i++) {
474 cx[i] = - (2*i+1) * (2*i+1) ;
479 double* xo =
new double[(tb->dim).taille] ;
482 for (
int i=0; i<(tb->dim).taille; i++) {
491 int borne_phi = np + 1 ;
492 if (np == 1) borne_phi = 1 ;
494 for (
int k=0 ; k< borne_phi ; k++) {
495 for (
int j=0 ; j<nt ; j++) {
496 for (
int i=0 ; i<nr ; i++ ) {
497 *xco = cx[j] * (*xci) ;
514void _d2sdtet2_t_cossin_cp(Tbl* tb,
int &)
518 if (tb->get_etat() == ETATZERO) {
523 assert(tb->get_etat() == ETATQCQ) ;
526 int nr = (tb->dim).dim[0] ;
527 int nt = (tb->dim).dim[1] ;
528 int np = (tb->dim).dim[2] ;
532 static double* cxp = 0 ;
533 static double* cxi = 0 ;
534 static int nt_pre =0 ;
539 cxp =
reinterpret_cast<double*
>(realloc(cxp, nt *
sizeof(
double))) ;
540 cxi =
reinterpret_cast<double*
>(realloc(cxi, nt *
sizeof(
double))) ;
541 for (
int i=0 ; i<nt ; i++) {
542 cxp[i] = - (2*i) * (2*i) ;
543 cxi[i] = - (2*i+1) * (2*i+1) ;
548 double* xo =
new double[(tb->dim).taille] ;
551 for (
int i=0; i<(tb->dim).taille; i++) {
562 for (k=0 ; k<np+1 ; k += 4) {
563 for (
int m=0 ; m<2 ; m++) {
564 for (
int j=0 ; j<nt ; j++) {
565 for (
int i=0 ; i<nr ; i++ ) {
566 *xco = cxp[j] * (*xci) ;
579 for (k=2 ; k<np+1 ; k += 4) {
580 for (
int m=0 ; m<2 ; m++) {
581 for (
int j=0 ; j<nt ; j++) {
582 for (
int i=0 ; i<nr ; i++ ) {
583 *xco = cxi[j] * (*xci) ;
603void _d2sdtet2_t_cossin_sp(Tbl* tb,
int &)
607 if (tb->get_etat() == ETATZERO) {
612 assert(tb->get_etat() == ETATQCQ) ;
615 int nr = (tb->dim).dim[0] ;
616 int nt = (tb->dim).dim[1] ;
617 int np = (tb->dim).dim[2] ;
621 static double* cxp = 0 ;
622 static double* cxi = 0 ;
623 static int nt_pre =0 ;
628 cxp =
reinterpret_cast<double*
>(realloc(cxp, nt *
sizeof(
double))) ;
629 cxi =
reinterpret_cast<double*
>(realloc(cxi, nt *
sizeof(
double))) ;
630 for (
int i=0 ; i<nt ; i++) {
631 cxp[i] = - (2*i) * (2*i) ;
632 cxi[i] = - (2*i+1) * (2*i+1) ;
637 double* xo =
new double[(tb->dim).taille] ;
640 for (
int i=0; i<(tb->dim).taille; i++) {
651 for (k=0 ; k<np+1 ; k += 4) {
652 for (
int m=0 ; m<2 ; m++) {
653 for (
int j=0 ; j<nt ; j++) {
654 for (
int i=0 ; i<nr ; i++ ) {
655 *xco = cxp[j] * (*xci) ;
668 for (k=2 ; k<np+1 ; k += 4) {
669 for (
int m=0 ; m<2 ; m++) {
670 for (
int j=0 ; j<nt ; j++) {
671 for (
int i=0 ; i<nr ; i++ ) {
672 *xco = cxi[j] * (*xci) ;
692void _d2sdtet2_t_cossin_si(Tbl* tb,
int &)
696 if (tb->get_etat() == ETATZERO) {
701 assert(tb->get_etat() == ETATQCQ) ;
704 int nr = (tb->dim).dim[0] ;
705 int nt = (tb->dim).dim[1] ;
706 int np = (tb->dim).dim[2] ;
710 static double* cxp = 0 ;
711 static double* cxi = 0 ;
712 static int nt_pre =0 ;
717 cxp =
reinterpret_cast<double*
>(realloc(cxp, nt *
sizeof(
double))) ;
718 cxi =
reinterpret_cast<double*
>(realloc(cxi, nt *
sizeof(
double))) ;
719 for (
int i=0 ; i<nt ; i++) {
720 cxp[i] = - (2*i) * (2*i) ;
721 cxi[i] = - (2*i+1) * (2*i+1) ;
726 double* xo =
new double[(tb->dim).taille] ;
729 for (
int i=0; i<(tb->dim).taille; i++) {
740 for (k=0 ; k<np+1 ; k += 4) {
741 for (
int m=0 ; m<2 ; m++) {
742 for (
int j=0 ; j<nt ; j++) {
743 for (
int i=0 ; i<nr ; i++ ) {
744 *xco = cxi[j] * (*xci) ;
757 for (k=2 ; k<np+1 ; k += 4) {
758 for (
int m=0 ; m<2 ; m++) {
759 for (
int j=0 ; j<nt ; j++) {
760 for (
int i=0 ; i<nr ; i++ ) {
761 *xco = cxp[j] * (*xci) ;
781void _d2sdtet2_t_cossin_c(Tbl* tb,
int &)
785 if (tb->get_etat() == ETATZERO) {
790 assert(tb->get_etat() == ETATQCQ) ;
793 int nr = (tb->dim).dim[0] ;
794 int nt = (tb->dim).dim[1] ;
795 int np = (tb->dim).dim[2] ;
799 static double* cxp = 0 ;
800 static double* cxi = 0 ;
801 static int nt_pre =0 ;
806 cxp =
reinterpret_cast<double*
>(realloc(cxp, nt *
sizeof(
double))) ;
807 cxi =
reinterpret_cast<double*
>(realloc(cxi, nt *
sizeof(
double))) ;
808 for (
int i=0 ; i<nt ; i++) {
815 double* xo =
new double[(tb->dim).taille] ;
818 for (
int i=0; i<(tb->dim).taille; i++) {
829 for (k=0 ; k<np+1 ; k += 4) {
830 for (
int m=0 ; m<2 ; m++) {
831 for (
int j=0 ; j<nt ; j++) {
832 for (
int i=0 ; i<nr ; i++ ) {
833 *xco = cxp[j] * (*xci) ;
846 for (k=2 ; k<np+1 ; k += 4) {
847 for (
int m=0 ; m<2 ; m++) {
848 for (
int j=0 ; j<nt ; j++) {
849 for (
int i=0 ; i<nr ; i++ ) {
850 *xco = cxi[j] * (*xci) ;
870void _d2sdtet2_t_cossin_s(Tbl* tb,
int &)
874 if (tb->get_etat() == ETATZERO) {
879 assert(tb->get_etat() == ETATQCQ) ;
882 int nr = (tb->dim).dim[0] ;
883 int nt = (tb->dim).dim[1] ;
884 int np = (tb->dim).dim[2] ;
888 static double* cxp = 0 ;
889 static double* cxi = 0 ;
890 static int nt_pre =0 ;
895 cxp =
reinterpret_cast<double*
>(realloc(cxp, nt *
sizeof(
double))) ;
896 cxi =
reinterpret_cast<double*
>(realloc(cxi, nt *
sizeof(
double))) ;
897 for (
int i=0 ; i<nt ; i++) {
904 double* xo =
new double[(tb->dim).taille] ;
907 for (
int i=0; i<(tb->dim).taille; i++) {
918 for (k=0 ; k<np+1 ; k += 4) {
919 for (
int m=0 ; m<2 ; m++) {
920 for (
int j=0 ; j<nt ; j++) {
921 for (
int i=0 ; i<nr ; i++ ) {
922 *xco = cxp[j] * (*xci) ;
935 for (k=2 ; k<np+1 ; k += 4) {
936 for (
int m=0 ; m<2 ; m++) {
937 for (
int j=0 ; j<nt ; j++) {
938 for (
int i=0 ; i<nr ; i++ ) {
939 *xco = cxi[j] * (*xci) ;