23char op_mult_st_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Non_class_members/Operators/op_mult_st.C,v 1.8 2014/10/13 08:53:25 j_novak Exp $" ;
89void _mult_st_pas_prevu(Tbl * tb,
int& base) {
90 cout <<
"mult_st pas prevu..." << endl ;
91 cout <<
"Tbl: " << tb <<
" base: " << base << endl ;
100void _mult_st_t_cos(Tbl* tb,
int & b)
104 if (tb->get_etat() == ETATZERO) {
105 int base_r = b &
MSQ_R ;
106 int base_p = b &
MSQ_P ;
115 b = base_r | base_p |
T_SIN ;
122 assert(tb->get_etat() == ETATQCQ) ;
125 int nr = (tb->dim).dim[0] ;
126 int nt = (tb->dim).dim[1] ;
127 int np = (tb->dim).dim[2] ;
131 double* som =
new double [nr] ;
134 double* xo =
new double[(tb->dim).taille] ;
137 for (
int i=0; i<(tb->dim).taille; i++) {
154 for (
int i=0 ; i<nr ; i++) {
155 som[i] = -0.5*xci[i] ;
160 for (
int j=nt-2 ; j > 0 ; j--) {
165 for (
int i=0 ; i<nr ; i++ ) {
166 som[i] += 0.5*xci[i] ;
170 for (
int i=0; i<nr; i++) {
171 som[i] = -0.5*xci[i] ;
176 for (
int i=0 ; i<nr ; i++ ) {
177 xco[i] += 0.5*xci[i] ;
180 for (
int i=0; i<nr; i++) {
192 for (
int k=2 ; k<np+1 ; k++) {
200 for (
int i=0 ; i<nr ; i++) {
201 som[i] = -0.5*xci[i] ;
206 for (
int j=nt-2 ; j > 0 ; j--) {
211 for (
int i=0 ; i<nr ; i++ ) {
212 som[i] += 0.5*xci[i] ;
216 for (
int i=0; i<nr; i++) {
217 som[i] = -0.5*xci[i] ;
222 for (
int i = 0; i<nr; i++) {
223 xco[i] += 0.5*xci[i] ;
226 for (
int i=0; i<nr; i++) {
242 int base_r = b &
MSQ_R ;
243 int base_p = b &
MSQ_P ;
252 b = base_r | base_p |
T_SIN ;
261void _mult_st_t_sin(Tbl* tb,
int & b)
264 if (tb->get_etat() == ETATZERO) {
265 int base_r = b &
MSQ_R ;
266 int base_p = b &
MSQ_P ;
275 b = base_r | base_p |
T_COS ;
282 assert(tb->get_etat() == ETATQCQ) ;
285 int nr = (tb->dim).dim[0] ;
286 int nt = (tb->dim).dim[1] ;
287 int np = (tb->dim).dim[2] ;
291 double* som =
new double [nr] ;
294 double* xo =
new double[(tb->dim).taille] ;
297 for (
int i=0; i<(tb->dim).taille; i++) {
314 for (
int i=0 ; i<nr ; i++) {
315 som[i] = 0.5*xci[i] ;
320 for (
int j=nt-2 ; j > 0 ; j--) {
325 for (
int i=0 ; i<nr ; i++ ) {
326 som[i] -= 0.5*xci[i] ;
330 for (
int i=0; i<nr; i++) {
331 som[i] = 0.5*xci[i] ;
337 for (
int i =0; i<nr ; i++) {
349 for (
int k=2 ; k<np+1 ; k++) {
357 for (
int i=0 ; i<nr ; i++) {
358 som[i] = 0.5*xci[i] ;
363 for (
int j=nt-2 ; j > 0 ; j--) {
368 for (
int i=0 ; i<nr ; i++ ) {
369 som[i] -= 0.5*xci[i] ;
373 for (
int i=0 ; i<nr ; i++ ) {
374 som[i] = 0.5*xci[i] ;
380 for (
int i=0; i<nr; i++) {
396 int base_r = b &
MSQ_R ;
397 int base_p = b &
MSQ_P ;
406 b = base_r | base_p |
T_COS ;
414void _mult_st_t_cos_p(Tbl* tb,
int & b)
418 if (tb->get_etat() == ETATZERO) {
419 int base_r = b &
MSQ_R ;
420 int base_p = b &
MSQ_P ;
421 b = base_r | base_p |
T_SIN_I ;
426 assert(tb->get_etat() == ETATQCQ) ;
429 int nr = (tb->dim).dim[0] ;
430 int nt = (tb->dim).dim[1] ;
431 int np = (tb->dim).dim[2] ;
435 double* som =
new double [nr] ;
438 double* xo =
new double[(tb->dim).taille] ;
441 for (
int i=0; i<(tb->dim).taille; i++) {
458 for (
int i=0 ; i<nr ; i++) {
459 som[i] = -0.5*xci[i] ;
464 for (
int j=nt-2 ; j > 0 ; j--) {
469 for (
int i=0 ; i<nr ; i++ ) {
470 som[i] += 0.5*xci[i] ;
472 som[i] = -0.5*xci[i] ;
479 for (
int i=0 ; i<nr ; i++ ) {
480 xco[i] = som[i]+xci[i] ;
492 for (
int k=2 ; k<np+1 ; k++) {
500 for (
int i=0 ; i<nr ; i++) {
501 som[i] = -0.5*xci[i] ;
506 for (
int j=nt-2 ; j > 0 ; j--) {
511 for (
int i=0 ; i<nr ; i++ ) {
512 som[i] += 0.5*xci[i] ;
514 som[i] = -0.5*xci[i] ;
520 for (
int i = 0; i<nr; i++) {
521 xco[i] = xci[i] + som[i] ;
536 int base_r = b &
MSQ_R ;
537 int base_p = b &
MSQ_P ;
538 b = base_r | base_p |
T_SIN_I ;
545void _mult_st_t_sin_p(Tbl* tb,
int & b)
548 if (tb->get_etat() == ETATZERO) {
549 int base_r = b &
MSQ_R ;
550 int base_p = b &
MSQ_P ;
551 b = base_r | base_p |
T_COS_I ;
556 assert(tb->get_etat() == ETATQCQ) ;
559 int nr = (tb->dim).dim[0] ;
560 int nt = (tb->dim).dim[1] ;
561 int np = (tb->dim).dim[2] ;
565 double* som =
new double [nr] ;
568 double* xo =
new double[(tb->dim).taille] ;
571 for (
int i=0; i<(tb->dim).taille; i++) {
588 for (
int i=0 ; i<nr ; i++) {
589 som[i] = 0.5*xci[i] ;
594 for (
int j=nt-2 ; j > 0 ; j--) {
599 for (
int i=0 ; i<nr ; i++ ) {
600 som[i] -= 0.5*xci[i] ;
602 som[i] = 0.5*xci[i] ;
609 for (
int i =0; i<nr ; i++) {
622 for (
int k=2 ; k<np+1 ; k++) {
630 for (
int i=0 ; i<nr ; i++) {
631 som[i] = 0.5*xci[i] ;
636 for (
int j=nt-2 ; j > 0 ; j--) {
641 for (
int i=0 ; i<nr ; i++ ) {
642 som[i] -= 0.5*xci[i] ;
644 som[i] = 0.5*xci[i] ;
650 for (
int i=0; i<nr; i++) {
666 int base_r = b &
MSQ_R ;
667 int base_p = b &
MSQ_P ;
668 b = base_r | base_p |
T_COS_I ;
675void _mult_st_t_sin_i(Tbl* tb,
int & b)
678 if (tb->get_etat() == ETATZERO) {
679 int base_r = b &
MSQ_R ;
680 int base_p = b &
MSQ_P ;
681 b = base_r | base_p |
T_COS_P ;
686 assert(tb->get_etat() == ETATQCQ) ;
689 int nr = (tb->dim).dim[0] ;
690 int nt = (tb->dim).dim[1] ;
691 int np = (tb->dim).dim[2] ;
695 double* som =
new double [nr] ;
698 double* xo =
new double[(tb->dim).taille] ;
701 for (
int i=0; i<(tb->dim).taille; i++) {
718 for (
int i=0 ; i<nr ; i++) {
723 for (
int j=nt-1 ; j > 0 ; j--) {
727 for (
int i=0 ; i<nr ; i++ ) {
728 som[i] -= 0.5*xci[i] ;
730 som[i] = 0.5*xci[i] ;
735 for (
int i=0 ; i<nr ; i++) {
747 for (
int k=2 ; k<np+1 ; k++) {
755 for (
int i=0 ; i<nr ; i++) {
760 for (
int j=nt-1 ; j > 0 ; j--) {
764 for (
int i=0 ; i<nr ; i++ ) {
765 som[i] -= 0.5*xci[i] ;
767 som[i] = 0.5*xci[i] ;
772 for (
int i=0 ; i<nr ; i++) {
788 int base_r = b &
MSQ_R ;
789 int base_p = b &
MSQ_P ;
790 b = base_r | base_p |
T_COS_P ;
796void _mult_st_t_cos_i(Tbl* tb,
int & b)
799 if (tb->get_etat() == ETATZERO) {
800 int base_r = b &
MSQ_R ;
801 int base_p = b &
MSQ_P ;
802 b = base_r | base_p |
T_SIN_P ;
807 assert(tb->get_etat() == ETATQCQ) ;
810 int nr = (tb->dim).dim[0] ;
811 int nt = (tb->dim).dim[1] ;
812 int np = (tb->dim).dim[2] ;
816 double* som =
new double [nr] ;
819 double* xo =
new double[(tb->dim).taille] ;
822 for (
int i=0; i<(tb->dim).taille; i++) {
839 for (
int i=0 ; i<nr ; i++) {
844 for (
int j=nt-1 ; j > 0 ; j--) {
848 for (
int i=0 ; i<nr ; i++ ) {
849 som[i] += 0.5*xci[i] ;
851 som[i] = -0.5*xci[i] ;
855 for (
int i=0; i<nr; i++) {
867 for (
int k=2 ; k<np+1 ; k++) {
875 for (
int i=0 ; i<nr ; i++) {
880 for (
int j=nt-1 ; j > 0 ; j--) {
884 for (
int i=0 ; i<nr ; i++ ) {
885 som[i] += 0.5*xci[i] ;
887 som[i] = -0.5*xci[i] ;
891 for (
int i=0; i<nr; i++) {
908 int base_r = b &
MSQ_R ;
909 int base_p = b &
MSQ_P ;
910 b = base_r | base_p |
T_SIN_P ;
916void _mult_st_t_cossin_cp(Tbl* tb,
int & b)
919 if (tb->get_etat() == ETATZERO) {
920 int base_r = b &
MSQ_R ;
921 int base_p = b &
MSQ_P ;
927 assert(tb->get_etat() == ETATQCQ) ;
930 int nr = (tb->dim).dim[0] ;
931 int nt = (tb->dim).dim[1] ;
932 int np = (tb->dim).dim[2] ;
936 double* som =
new double [nr] ;
939 double* xo =
new double[(tb->dim).taille] ;
942 for (
int i=0; i<(tb->dim).taille; i++) {
959 for (
int i=0 ; i<nr ; i++) {
960 som[i] = -0.5*xci[i] ;
965 for (
int j=nt-2 ; j > 0 ; j--) {
970 for (
int i=0 ; i<nr ; i++ ) {
971 som[i] += 0.5*xci[i] ;
973 som[i] = -0.5*xci[i] ;
981 for (
int i = 0; i<nr; i++) {
982 xco[i] = xci[i] + som[i] ;
994 for (
int k=2 ; k<np+1 ; k++) {
1001 xci += nr * (nt-1) ;
1002 xco += nr * (nt-1) ;
1005 for (
int i=0 ; i<nr ; i++) {
1006 som[i] = -0.5*xci[i] ;
1011 for (
int j=nt-2 ; j > 0 ; j--) {
1016 for (
int i=0 ; i<nr ; i++ ) {
1017 som[i] += 0.5*xci[i] ;
1019 som[i] = -0.5*xci[i] ;
1025 for (
int i = 0; i<nr; i++) {
1026 xco[i] = xci[i] + som[i] ;
1036 xci += nr * (nt-1) ;
1037 xco += nr * (nt-1) ;
1040 for (
int i=0 ; i<nr ; i++) {
1045 for (
int j=nt-1 ; j > 0 ; j--) {
1049 for (
int i=0 ; i<nr ; i++ ) {
1050 som[i] -= 0.5*xci[i] ;
1052 som[i] = 0.5*xci[i] ;
1057 for (
int i=0 ; i<nr ; i++) {
1075 int base_r = b &
MSQ_R ;
1076 int base_p = b &
MSQ_P ;
1083void _mult_st_t_cossin_ci(Tbl* tb,
int & b)
1086 if (tb->get_etat() == ETATZERO) {
1087 int base_r = b &
MSQ_R ;
1088 int base_p = b &
MSQ_P ;
1094 assert(tb->get_etat() == ETATQCQ) ;
1097 int nr = (tb->dim).dim[0] ;
1098 int nt = (tb->dim).dim[1] ;
1099 int np = (tb->dim).dim[2] ;
1103 double* som =
new double [nr] ;
1106 double* xo =
new double[(tb->dim).taille] ;
1109 for (
int i=0; i<(tb->dim).taille; i++) {
1114 double* xi = tb->t ;
1122 xci += nr * (nt-1) ;
1123 xco += nr * (nt-1) ;
1126 for (
int i=0 ; i<nr ; i++) {
1131 for (
int j=nt-1 ; j > 0 ; j--) {
1135 for (
int i=0 ; i<nr ; i++ ) {
1136 som[i] += 0.5*xci[i] ;
1138 som[i] = -0.5*xci[i] ;
1142 for (
int i=0; i<nr; i++) {
1155 for (
int k=2 ; k<np+1 ; k++) {
1162 xci += nr * (nt-1) ;
1163 xco += nr * (nt-1) ;
1166 for (
int i=0 ; i<nr ; i++) {
1171 for (
int j=nt-1 ; j > 0 ; j--) {
1175 for (
int i=0 ; i<nr ; i++ ) {
1176 som[i] += 0.5*xci[i] ;
1178 som[i] = -0.5*xci[i] ;
1182 for (
int i=0; i<nr; i++) {
1194 xci += nr * (nt-1) ;
1195 xco += nr * (nt-1) ;
1198 for (
int i=0 ; i<nr ; i++) {
1199 som[i] = 0.5*xci[i] ;
1204 for (
int j=nt-2 ; j > 0 ; j--) {
1209 for (
int i=0 ; i<nr ; i++ ) {
1210 som[i] -= 0.5*xci[i] ;
1212 som[i] = 0.5*xci[i] ;
1218 for (
int i=0; i<nr; i++) {
1236 int base_r = b &
MSQ_R ;
1237 int base_p = b &
MSQ_P ;
1244void _mult_st_t_cossin_si(Tbl* tb,
int & b)
1247 if (tb->get_etat() == ETATZERO) {
1248 int base_r = b &
MSQ_R ;
1249 int base_p = b &
MSQ_P ;
1255 assert(tb->get_etat() == ETATQCQ) ;
1258 int nr = (tb->dim).dim[0] ;
1259 int nt = (tb->dim).dim[1] ;
1260 int np = (tb->dim).dim[2] ;
1264 double* som =
new double [nr] ;
1267 double* xo =
new double[(tb->dim).taille] ;
1270 for (
int i=0; i<(tb->dim).taille; i++) {
1275 double* xi = tb->t ;
1283 xci += nr * (nt-1) ;
1284 xco += nr * (nt-1) ;
1287 for (
int i=0 ; i<nr ; i++) {
1292 for (
int j=nt-1 ; j > 0 ; j--) {
1296 for (
int i=0 ; i<nr ; i++ ) {
1297 som[i] -= 0.5*xci[i] ;
1299 som[i] = 0.5*xci[i] ;
1304 for (
int i=0 ; i<nr ; i++) {
1316 for (
int k=2 ; k<np+1 ; k++) {
1323 xci += nr * (nt-1) ;
1324 xco += nr * (nt-1) ;
1327 for (
int i=0 ; i<nr ; i++) {
1332 for (
int j=nt-1 ; j > 0 ; j--) {
1336 for (
int i=0 ; i<nr ; i++ ) {
1337 som[i] -= 0.5*xci[i] ;
1339 som[i] = 0.5*xci[i] ;
1344 for (
int i=0 ; i<nr ; i++) {
1355 xci += nr * (nt-1) ;
1356 xco += nr * (nt-1) ;
1359 for (
int i=0 ; i<nr ; i++) {
1360 som[i] = -0.5*xci[i] ;
1365 for (
int j=nt-2 ; j > 0 ; j--) {
1370 for (
int i=0 ; i<nr ; i++ ) {
1371 som[i] += 0.5*xci[i] ;
1373 som[i] = -0.5*xci[i] ;
1379 for (
int i = 0; i<nr; i++) {
1380 xco[i] = xci[i] + som[i] ;
1397 int base_r = b &
MSQ_R ;
1398 int base_p = b &
MSQ_P ;
1406void _mult_st_t_cossin_sp(Tbl* tb,
int & b)
1409 if (tb->get_etat() == ETATZERO) {
1410 int base_r = b &
MSQ_R ;
1411 int base_p = b &
MSQ_P ;
1417 assert(tb->get_etat() == ETATQCQ) ;
1420 int nr = (tb->dim).dim[0] ;
1421 int nt = (tb->dim).dim[1] ;
1422 int np = (tb->dim).dim[2] ;
1426 double* som =
new double [nr] ;
1429 double* xo =
new double[(tb->dim).taille] ;
1432 for (
int i=0; i<(tb->dim).taille; i++) {
1437 double* xi = tb->t ;
1445 xci += nr * (nt-1) ;
1446 xco += nr * (nt-1) ;
1449 for (
int i=0 ; i<nr ; i++) {
1450 som[i] = 0.5*xci[i] ;
1455 for (
int j=nt-2 ; j > 0 ; j--) {
1460 for (
int i=0 ; i<nr ; i++ ) {
1461 som[i] -= 0.5*xci[i] ;
1463 som[i] = 0.5*xci[i] ;
1471 for (
int i=0; i<nr; i++) {
1483 for (
int k=2 ; k<np+1 ; k++) {
1490 xci += nr * (nt-1) ;
1491 xco += nr * (nt-1) ;
1494 for (
int i=0 ; i<nr ; i++) {
1499 for (
int j=nt-1 ; j > 0 ; j--) {
1503 for (
int i=0 ; i<nr ; i++ ) {
1504 som[i] += 0.5*xci[i] ;
1506 som[i] = -0.5*xci[i] ;
1510 for (
int i=0; i<nr; i++) {
1522 xci += nr * (nt-1) ;
1523 xco += nr * (nt-1) ;
1526 for (
int i=0 ; i<nr ; i++) {
1527 som[i] = 0.5*xci[i] ;
1532 for (
int j=nt-2 ; j > 0 ; j--) {
1537 for (
int i=0 ; i<nr ; i++ ) {
1538 som[i] -= 0.5*xci[i] ;
1540 som[i] = 0.5*xci[i] ;
1546 for (
int i=0; i<nr; i++) {
1564 int base_r = b &
MSQ_R ;
1565 int base_p = b &
MSQ_P ;
1573void _mult_st_t_cossin_c(Tbl* tb,
int & b)
1576 if (tb->get_etat() == ETATZERO) {
1577 int base_r = b &
MSQ_R ;
1578 int base_p = b &
MSQ_P ;
1594 assert(tb->get_etat() == ETATQCQ) ;
1597 int nr = (tb->dim).dim[0] ;
1598 int nt = (tb->dim).dim[1] ;
1599 int np = (tb->dim).dim[2] ;
1603 double* som =
new double [nr] ;
1606 double* xo =
new double[(tb->dim).taille] ;
1609 for (
int i=0; i<(tb->dim).taille; i++) {
1614 double* xi = tb->t ;
1622 xci += nr * (nt-1) ;
1623 xco += nr * (nt-1) ;
1626 for (
int i=0 ; i<nr ; i++) {
1627 som[i] = -0.5*xci[i] ;
1633 for (
int j=nt-2 ; j > 0 ; j--) {
1638 for (
int i=0 ; i<nr ; i++ ) {
1639 som[i] += 0.5*xci[i] ;
1643 for (
int i=0 ; i<nr ; i++ ) {
1644 som[i] = -0.5*xci[i] ;
1649 for (
int i=0; i<nr; i++) {
1650 xco[i] += 0.5*xci[i] ;
1653 for (
int i = 0; i<nr; i++) {
1665 for (
int k=2 ; k<np+1 ; k++) {
1672 xci += nr * (nt-1) ;
1673 xco += nr * (nt-1) ;
1676 for (
int i=0 ; i<nr ; i++) {
1677 som[i] = -0.5*xci[i] ;
1682 for (
int j=nt-2 ; j > 0 ; j--) {
1687 for (
int i=0 ; i<nr ; i++ ) {
1688 som[i] += 0.5*xci[i] ;
1692 for (
int i=0 ; i<nr ; i++ ) {
1693 som[i] = -0.5*xci[i] ;
1698 for (
int i=0; i<nr; i++) {
1699 xco[i] += 0.5*xci[i] ;
1702 for (
int i = 0; i<nr; i++) {
1713 xci += nr * (nt-1) ;
1714 xco += nr * (nt-1) ;
1717 for (
int i=0 ; i<nr ; i++) {
1718 som[i] = 0.5*xci[i] ;
1723 for (
int j=nt-2 ; j > 0 ; j--) {
1728 for (
int i=0 ; i<nr ; i++ ) {
1729 som[i] -= 0.5*xci[i] ;
1733 for (
int i=0 ; i<nr ; i++ ) {
1734 som[i] = 0.5*xci[i] ;
1740 for (
int i=0 ; i<nr ; i++) {
1758 int base_r = b &
MSQ_R ;
1759 int base_p = b &
MSQ_P ;
1776void _mult_st_t_cossin_s(Tbl* tb,
int & b)
1779 if (tb->get_etat() == ETATZERO) {
1780 int base_r = b &
MSQ_R ;
1781 int base_p = b &
MSQ_P ;
1797 assert(tb->get_etat() == ETATQCQ) ;
1800 int nr = (tb->dim).dim[0] ;
1801 int nt = (tb->dim).dim[1] ;
1802 int np = (tb->dim).dim[2] ;
1806 double* som =
new double [nr] ;
1809 double* xo =
new double[(tb->dim).taille] ;
1812 for (
int i=0; i<(tb->dim).taille; i++) {
1817 double* xi = tb->t ;
1825 xci += nr * (nt-1) ;
1826 xco += nr * (nt-1) ;
1829 for (
int i=0 ; i<nr ; i++) {
1830 som[i] = 0.5*xci[i] ;
1835 for (
int j=nt-2 ; j > 0 ; j--) {
1840 for (
int i=0 ; i<nr ; i++ ) {
1841 som[i] -= 0.5*xci[i] ;
1845 for (
int i=0 ; i<nr ; i++ ) {
1846 som[i] = 0.5*xci[i] ;
1852 for (
int i=0; i<nr; i++) {
1863 for (
int k=2 ; k<np+1 ; k++) {
1870 xci += nr * (nt-1) ;
1871 xco += nr * (nt-1) ;
1874 for (
int i=0 ; i<nr ; i++) {
1875 som[i] = -0.5*xci[i] ;
1880 for (
int j=nt-2 ; j > 0 ; j--) {
1885 for (
int i=0 ; i<nr ; i++ ) {
1886 som[i] += 0.5*xci[i] ;
1890 for (
int i=0 ; i<nr ; i++ ) {
1891 som[i] = -0.5*xci[i] ;
1895 for (
int i=0; i<nr; i++) {
1896 xco[i] += 0.5*xci[i] ;
1899 for (
int i=0; i<nr; i++) {
1911 xci += nr * (nt-1) ;
1912 xco += nr * (nt-1) ;
1915 for (
int i=0 ; i<nr ; i++) {
1916 som[i] = 0.5*xci[i] ;
1921 for (
int j=nt-2 ; j > 0 ; j--) {
1926 for (
int i=0 ; i<nr ; i++ ) {
1927 som[i] -= 0.5*xci[i] ;
1931 for (
int i=0 ; i<nr ; i++ ) {
1932 som[i] = 0.5*xci[i] ;
1938 for (
int i=0; i<nr; i++) {
1956 int base_r = b &
MSQ_R ;
1957 int base_p = b &
MSQ_P ;
#define T_COSSIN_SP
sin pair-cos impair alternes, sin pour m=0
#define T_SIN_P
dev. sin seulement, harmoniques paires
#define MSQ_R
Extraction de l'info sur R.
#define T_COSSIN_S
dev. cos-sin alternes, sin pour m=0
#define T_COSSIN_SI
sin impair-cos pair alternes, sin pour m=0
#define R_CHEBPI_I
Cheb. pair-impair suivant l impair pour l=0.
#define T_COS_P
dev. cos seulement, harmoniques paires
#define T_COSSIN_CI
cos impair-sin pair alternes, cos pour m=0
#define T_COSSIN_CP
cos pair-sin impair alternes, cos pour m=0
#define T_SIN_I
dev. sin seulement, harmoniques impaires
#define T_COS
dev. cos seulement
#define MSQ_P
Extraction de l'info sur Phi.
#define T_SIN
dev. sin seulement
#define T_COS_I
dev. cos seulement, harmoniques impaires
#define T_COSSIN_C
dev. cos-sin alternes, cos pour m=0
#define R_CHEBPI_P
Cheb. pair-impair suivant l pair pour l=0.