23char op_mult_ct_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Non_class_members/Operators/op_mult_ct.C,v 1.8 2014/10/13 08:53:25 j_novak Exp $" ;
89void _mult_ct_pas_prevu(Tbl * tb,
int& base) {
90 cout <<
"mult_ct pas prevu..." << endl ;
91 cout <<
"Tbl: " << tb <<
" base: " << base << endl ;
100void _mult_ct_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_COS ;
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] ;
175 for (
int i=0 ; i<nr ; i++ ) {
176 xco[i] += 0.5*xci[i] ;
179 for (
int i=0; i<nr; i++)
191 for (
int k=2 ; k<np+1 ; k++) {
199 for (
int i=0 ; i<nr ; i++) {
200 som[i] = 0.5*xci[i] ;
205 for (
int j=nt-2 ; j > 0 ; j--) {
210 for (
int i=0 ; i<nr ; i++ ) {
211 som[i] += 0.5*xci[i] ;
215 for (
int i=0; i<nr; i++)
216 som[i] = 0.5*xci[i] ;
220 for (
int i = 0; i<nr; i++) {
221 xco[i] += 0.5*xci[i] ;
224 for (
int i=0; i<nr; i++)
239 int base_r = b &
MSQ_R ;
240 int base_p = b &
MSQ_P ;
249 b = base_r | base_p |
T_COS ;
258void _mult_ct_t_sin(Tbl* tb,
int & b)
261 if (tb->get_etat() == ETATZERO) {
262 int base_r = b &
MSQ_R ;
263 int base_p = b &
MSQ_P ;
272 b = base_r | base_p |
T_SIN ;
279 assert(tb->get_etat() == ETATQCQ) ;
282 int nr = (tb->dim).dim[0] ;
283 int nt = (tb->dim).dim[1] ;
284 int np = (tb->dim).dim[2] ;
288 double* som =
new double [nr] ;
291 double* xo =
new double[(tb->dim).taille] ;
294 for (
int i=0; i<(tb->dim).taille; i++) {
311 for (
int i=0 ; i<nr ; i++) {
312 som[i] = 0.5*xci[i] ;
317 for (
int j=nt-2 ; j > 0 ; j--) {
322 for (
int i=0 ; i<nr ; i++ ) {
323 som[i] += 0.5*xci[i] ;
327 for (
int i=0; i<nr; i++) {
328 som[i] = 0.5*xci[i] ;
334 for (
int i =0; i<nr ; i++) {
346 for (
int k=2 ; k<np+1 ; k++) {
354 for (
int i=0 ; i<nr ; i++) {
355 som[i] = 0.5*xci[i] ;
360 for (
int j=nt-2 ; j > 0 ; j--) {
365 for (
int i=0 ; i<nr ; i++ ) {
366 som[i] += 0.5*xci[i] ;
370 for (
int i=0; i<nr; i++) {
371 som[i] = 0.5*xci[i] ;
377 for (
int i=0; i<nr; i++) {
393 int base_r = b &
MSQ_R ;
394 int base_p = b &
MSQ_P ;
403 b = base_r | base_p |
T_SIN ;
411void _mult_ct_t_cos_p(Tbl* tb,
int & b)
415 if (tb->get_etat() == ETATZERO) {
416 int base_r = b &
MSQ_R ;
417 int base_p = b &
MSQ_P ;
418 b = base_r | base_p |
T_COS_I ;
423 assert(tb->get_etat() == ETATQCQ) ;
426 int nr = (tb->dim).dim[0] ;
427 int nt = (tb->dim).dim[1] ;
428 int np = (tb->dim).dim[2] ;
432 double* som =
new double [nr] ;
435 double* xo =
new double[(tb->dim).taille] ;
438 for (
int i=0; i<(tb->dim).taille; i++) {
455 for (
int i=0 ; i<nr ; i++) {
456 som[i] = 0.5*xci[i] ;
461 for (
int j=nt-2 ; j > 0 ; j--) {
466 for (
int i=0 ; i<nr ; i++ ) {
467 som[i] += 0.5*xci[i] ;
469 som[i] = 0.5*xci[i] ;
477 for (
int i=0 ; i<nr ; i++ ) {
478 xco[i] = som[i]+xci[i] ;
491 for (
int k=2 ; k<np+1 ; k++) {
499 for (
int i=0 ; i<nr ; i++) {
500 som[i] = 0.5*xci[i] ;
505 for (
int j=nt-2 ; j > 0 ; j--) {
510 for (
int i=0 ; i<nr ; i++ ) {
511 som[i] += 0.5*xci[i] ;
513 som[i] = 0.5*xci[i] ;
519 for (
int i = 0; i<nr; i++) {
520 xco[i] = xci[i] + som[i] ;
535 int base_r = b &
MSQ_R ;
536 int base_p = b &
MSQ_P ;
537 b = base_r | base_p |
T_COS_I ;
544void _mult_ct_t_sin_p(Tbl* tb,
int & b)
547 if (tb->get_etat() == ETATZERO) {
548 int base_r = b &
MSQ_R ;
549 int base_p = b &
MSQ_P ;
550 b = base_r | base_p |
T_SIN_I ;
555 assert(tb->get_etat() == ETATQCQ) ;
558 int nr = (tb->dim).dim[0] ;
559 int nt = (tb->dim).dim[1] ;
560 int np = (tb->dim).dim[2] ;
564 double* som =
new double [nr] ;
567 double* xo =
new double[(tb->dim).taille] ;
570 for (
int i=0; i<(tb->dim).taille; i++) {
587 for (
int i=0 ; i<nr ; i++) {
588 som[i] = 0.5*xci[i] ;
593 for (
int j=nt-2 ; j > 0 ; j--) {
598 for (
int i=0 ; i<nr ; i++ ) {
599 som[i] += 0.5*xci[i] ;
601 som[i] = 0.5*xci[i] ;
608 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_SIN_I ;
675void _mult_ct_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_SIN_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] ;
734 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] ;
771 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_SIN_P ;
796void _mult_ct_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_COS_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] ;
856 for (
int i=0 ; i<nr ; i++) {
868 for (
int k=2 ; k<np+1 ; k++) {
876 for (
int i=0 ; i<nr ; i++) {
881 for (
int j=nt-1 ; j > 0 ; j--) {
885 for (
int i=0 ; i<nr ; i++ ) {
886 som[i] += 0.5*xci[i] ;
888 som[i] = 0.5*xci[i] ;
893 for (
int i=0 ; i<nr ; i++) {
909 int base_r = b &
MSQ_R ;
910 int base_p = b &
MSQ_P ;
911 b = base_r | base_p |
T_COS_P ;
917void _mult_ct_t_cossin_cp(Tbl* tb,
int & b)
920 if (tb->get_etat() == ETATZERO) {
921 int base_r = b &
MSQ_R ;
922 int base_p = b &
MSQ_P ;
928 assert(tb->get_etat() == ETATQCQ) ;
931 int nr = (tb->dim).dim[0] ;
932 int nt = (tb->dim).dim[1] ;
933 int np = (tb->dim).dim[2] ;
937 double* som =
new double [nr] ;
940 double* xo =
new double[(tb->dim).taille] ;
943 for (
int i=0; i<(tb->dim).taille; i++) {
960 for (
int i=0 ; i<nr ; i++) {
961 som[i] = 0.5*xci[i] ;
966 for (
int j=nt-2 ; j > 0 ; j--) {
971 for (
int i=0 ; i<nr ; i++ ) {
972 som[i] += 0.5*xci[i] ;
974 som[i] = 0.5*xci[i] ;
982 for (
int i = 0; i<nr; i++) {
983 xco[i] = xci[i] + som[i] ;
995 for (
int k=2 ; k<np+1 ; k++) {
1002 xci += nr * (nt-1) ;
1003 xco += nr * (nt-1) ;
1006 for (
int i=0 ; i<nr ; i++) {
1007 som[i] = 0.5*xci[i] ;
1012 for (
int j=nt-2 ; j > 0 ; j--) {
1017 for (
int i=0 ; i<nr ; i++ ) {
1018 som[i] += 0.5*xci[i] ;
1020 som[i] = 0.5*xci[i] ;
1026 for (
int i = 0; i<nr; i++) {
1027 xco[i] = xci[i] + som[i] ;
1037 xci += nr * (nt-1) ;
1038 xco += nr * (nt-1) ;
1041 for (
int i=0 ; i<nr ; i++) {
1046 for (
int j=nt-1 ; j > 0 ; j--) {
1050 for (
int i=0 ; i<nr ; i++ ) {
1051 som[i] += 0.5*xci[i] ;
1053 som[i] = 0.5*xci[i] ;
1057 for (
int i=0; i<nr; i++) {
1076 int base_r = b &
MSQ_R ;
1077 int base_p = b &
MSQ_P ;
1084void _mult_ct_t_cossin_ci(Tbl* tb,
int & b)
1087 if (tb->get_etat() == ETATZERO) {
1088 int base_r = b &
MSQ_R ;
1089 int base_p = b &
MSQ_P ;
1095 assert(tb->get_etat() == ETATQCQ) ;
1098 int nr = (tb->dim).dim[0] ;
1099 int nt = (tb->dim).dim[1] ;
1100 int np = (tb->dim).dim[2] ;
1104 double* som =
new double [nr] ;
1107 double* xo =
new double[(tb->dim).taille] ;
1110 for (
int i=0; i<(tb->dim).taille; i++) {
1115 double* xi = tb->t ;
1123 xci += nr * (nt-1) ;
1124 xco += nr * (nt-1) ;
1127 for (
int i=0 ; i<nr ; i++) {
1132 for (
int j=nt-1 ; j > 0 ; j--) {
1136 for (
int i=0 ; i<nr ; i++ ) {
1137 som[i] += 0.5*xci[i] ;
1139 som[i] = 0.5*xci[i] ;
1144 for (
int i=0 ; i<nr ; i++) {
1156 for (
int k=2 ; k<np+1 ; k++) {
1163 xci += nr * (nt-1) ;
1164 xco += nr * (nt-1) ;
1167 for (
int i=0 ; i<nr ; i++) {
1172 for (
int j=nt-1 ; j > 0 ; j--) {
1176 for (
int i=0 ; i<nr ; i++ ) {
1177 som[i] += 0.5*xci[i] ;
1179 som[i] = 0.5*xci[i] ;
1184 for (
int i=0 ; i<nr ; i++) {
1195 xci += nr * (nt-1) ;
1196 xco += nr * (nt-1) ;
1199 for (
int i=0 ; i<nr ; i++) {
1200 som[i] = 0.5*xci[i] ;
1205 for (
int j=nt-2 ; j > 0 ; j--) {
1210 for (
int i=0 ; i<nr ; i++ ) {
1211 som[i] += 0.5*xci[i] ;
1213 som[i] = 0.5*xci[i] ;
1219 for (
int i=0; i<nr; i++) {
1237 int base_r = b &
MSQ_R ;
1238 int base_p = b &
MSQ_P ;
1245void _mult_ct_t_cossin_si(Tbl* tb,
int & b)
1248 if (tb->get_etat() == ETATZERO) {
1249 int base_r = b &
MSQ_R ;
1250 int base_p = b &
MSQ_P ;
1256 assert(tb->get_etat() == ETATQCQ) ;
1259 int nr = (tb->dim).dim[0] ;
1260 int nt = (tb->dim).dim[1] ;
1261 int np = (tb->dim).dim[2] ;
1265 double* som =
new double [nr] ;
1268 double* xo =
new double[(tb->dim).taille] ;
1271 for (
int i=0; i<(tb->dim).taille; i++) {
1276 double* xi = tb->t ;
1284 xci += nr * (nt-1) ;
1285 xco += nr * (nt-1) ;
1288 for (
int i=0 ; i<nr ; i++) {
1293 for (
int j=nt-1 ; j > 0 ; j--) {
1297 for (
int i=0 ; i<nr ; i++ ) {
1298 som[i] += 0.5*xci[i] ;
1300 som[i] = 0.5*xci[i] ;
1304 for (
int i=0; i<nr; i++) {
1317 for (
int k=2 ; k<np+1 ; k++) {
1324 xci += nr * (nt-1) ;
1325 xco += nr * (nt-1) ;
1328 for (
int i=0 ; i<nr ; i++) {
1333 for (
int j=nt-1 ; j > 0 ; j--) {
1337 for (
int i=0 ; i<nr ; i++ ) {
1338 som[i] += 0.5*xci[i] ;
1340 som[i] = 0.5*xci[i] ;
1344 for (
int i=0; i<nr; i++) {
1356 xci += nr * (nt-1) ;
1357 xco += nr * (nt-1) ;
1360 for (
int i=0 ; i<nr ; i++) {
1361 som[i] = 0.5*xci[i] ;
1366 for (
int j=nt-2 ; j > 0 ; j--) {
1371 for (
int i=0 ; i<nr ; i++ ) {
1372 som[i] += 0.5*xci[i] ;
1374 som[i] = 0.5*xci[i] ;
1380 for (
int i = 0; i<nr; i++) {
1381 xco[i] = xci[i] + som[i] ;
1398 int base_r = b &
MSQ_R ;
1399 int base_p = b &
MSQ_P ;
1407void _mult_ct_t_cossin_sp(Tbl* tb,
int & b)
1410 if (tb->get_etat() == ETATZERO) {
1411 int base_r = b &
MSQ_R ;
1412 int base_p = b &
MSQ_P ;
1418 assert(tb->get_etat() == ETATQCQ) ;
1421 int nr = (tb->dim).dim[0] ;
1422 int nt = (tb->dim).dim[1] ;
1423 int np = (tb->dim).dim[2] ;
1427 double* som =
new double [nr] ;
1430 double* xo =
new double[(tb->dim).taille] ;
1433 for (
int i=0; i<(tb->dim).taille; i++) {
1438 double* xi = tb->t ;
1446 xci += nr * (nt-1) ;
1447 xco += nr * (nt-1) ;
1450 for (
int i=0 ; i<nr ; i++) {
1451 som[i] = 0.5*xci[i] ;
1456 for (
int j=nt-2 ; j > 0 ; j--) {
1461 for (
int i=0 ; i<nr ; i++ ) {
1462 som[i] += 0.5*xci[i] ;
1464 som[i] = 0.5*xci[i] ;
1472 for (
int i=0; i<nr; i++) {
1484 for (
int k=2 ; k<np+1 ; k++) {
1491 xci += nr * (nt-1) ;
1492 xco += nr * (nt-1) ;
1495 for (
int i=0 ; i<nr ; i++) {
1500 for (
int j=nt-1 ; j > 0 ; j--) {
1504 for (
int i=0 ; i<nr ; i++ ) {
1505 som[i] += 0.5*xci[i] ;
1507 som[i] = 0.5*xci[i] ;
1512 for (
int i=0 ; i<nr ; i++) {
1523 xci += nr * (nt-1) ;
1524 xco += nr * (nt-1) ;
1527 for (
int i=0 ; i<nr ; i++) {
1528 som[i] = 0.5*xci[i] ;
1533 for (
int j=nt-2 ; j > 0 ; j--) {
1538 for (
int i=0 ; i<nr ; i++ ) {
1539 som[i] += 0.5*xci[i] ;
1541 som[i] = 0.5*xci[i] ;
1547 for (
int i=0; i<nr; i++) {
1565 int base_r = b &
MSQ_R ;
1566 int base_p = b &
MSQ_P ;
1574void _mult_ct_t_cossin_c(Tbl* tb,
int & b)
1577 if (tb->get_etat() == ETATZERO) {
1578 int base_r = b &
MSQ_R ;
1579 int base_p = b &
MSQ_P ;
1595 assert(tb->get_etat() == ETATQCQ) ;
1598 int nr = (tb->dim).dim[0] ;
1599 int nt = (tb->dim).dim[1] ;
1600 int np = (tb->dim).dim[2] ;
1604 double* som =
new double [nr] ;
1607 double* xo =
new double[(tb->dim).taille] ;
1610 for (
int i=0; i<(tb->dim).taille; i++) {
1615 double* xi = tb->t ;
1623 xci += nr * (nt-1) ;
1624 xco += nr * (nt-1) ;
1627 for (
int i=0 ; i<nr ; i++) {
1628 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++) {
1759 int base_r = b &
MSQ_R ;
1760 int base_p = b &
MSQ_P ;
1777void _mult_ct_t_cossin_s(Tbl* tb,
int & b)
1780 if (tb->get_etat() == ETATZERO) {
1781 int base_r = b &
MSQ_R ;
1782 int base_p = b &
MSQ_P ;
1798 assert(tb->get_etat() == ETATQCQ) ;
1801 int nr = (tb->dim).dim[0] ;
1802 int nt = (tb->dim).dim[1] ;
1803 int np = (tb->dim).dim[2] ;
1807 double* som =
new double [nr] ;
1810 double* xo =
new double[(tb->dim).taille] ;
1813 for (
int i=0; i<(tb->dim).taille; i++) {
1818 double* xi = tb->t ;
1826 xci += nr * (nt-1) ;
1827 xco += nr * (nt-1) ;
1830 for (
int i=0 ; i<nr ; i++) {
1831 som[i] = 0.5*xci[i] ;
1836 for (
int j=nt-2 ; j > 0 ; j--) {
1841 for (
int i=0 ; i<nr ; i++ ) {
1842 som[i] += 0.5*xci[i] ;
1846 for (
int i=0 ; i<nr ; i++ ) {
1847 som[i] = 0.5*xci[i] ;
1853 for (
int i=0; i<nr; i++) {
1864 for (
int k=2 ; k<np+1 ; k++) {
1871 xci += nr * (nt-1) ;
1872 xco += nr * (nt-1) ;
1875 for (
int i=0 ; i<nr ; i++) {
1876 som[i] = 0.5*xci[i] ;
1881 for (
int j=nt-2 ; j > 0 ; j--) {
1886 for (
int i=0 ; i<nr ; i++ ) {
1887 som[i] += 0.5*xci[i] ;
1891 for (
int i=0 ; i<nr ; i++ ) {
1892 som[i] = 0.5*xci[i] ;
1898 for (
int i=0; i<nr; i++) {
1899 xco[i] += 0.5*xci[i] ;
1902 for (
int i=0 ; i<nr ; i++) {
1913 xci += nr * (nt-1) ;
1914 xco += nr * (nt-1) ;
1917 for (
int i=0 ; i<nr ; i++) {
1918 som[i] = 0.5*xci[i] ;
1923 for (
int j=nt-2 ; j > 0 ; j--) {
1928 for (
int i=0 ; i<nr ; i++ ) {
1929 som[i] += 0.5*xci[i] ;
1933 for (
int i=0 ; i<nr ; i++ ) {
1934 som[i] = 0.5*xci[i] ;
1940 for (
int i=0; i<nr; i++) {
1958 int base_r = b &
MSQ_R ;
1959 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.