25char som_tet_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Non_class_members/Operators/som_tet.C,v 1.8 2014/10/13 08:53:27 j_novak Exp $" ;
110void som_tet_pas_prevu
111 (
double* ,
const int,
const int,
const double,
double*) {
112 cout <<
"Mtbl_cf::val_point: theta basis not implemented yet ! "
122 (
double* ti,
const int nt,
const int np,
123 const double tet,
double* to) {
131 double* cosinus =
new double [nt] ;
132 for (j=0 ; j<nt ; j++) {
133 cosinus[j] =
cos(j * tet) ;
138 for (j=0 ; j<nt ; j++) {
139 *po += (*pi) * cosinus[j] ;
151 for (k=2 ; k<np+1 ; k++) {
153 for (j=0 ; j<nt ; j++) {
154 *po += (*pi) * cosinus[j] ;
172 (
double* ti,
const int nt,
const int np,
173 const double tet,
double* to) {
181 double* cosinus =
new double [nt] ;
182 for (j=0 ; j<nt ; j++) {
183 cosinus[j] =
cos(2*j * tet) ;
188 for (j=0 ; j<nt ; j++) {
189 *po += (*pi) * cosinus[j] ;
201 for (k=2 ; k<np+1 ; k++) {
203 for (j=0 ; j<nt ; j++) {
204 *po += (*pi) * cosinus[j] ;
221 (
double* ti,
const int nt,
const int np,
222 const double tet,
double* to) {
230 double* cosinus =
new double [nt] ;
231 for (j=0 ; j<nt-1 ; j++) {
232 cosinus[j] =
cos((2*j+1) * tet) ;
238 for (j=0 ; j<nt ; j++) {
239 *po += (*pi) * cosinus[j] ;
251 for (k=2 ; k<np+1 ; k++) {
253 for (j=0 ; j<nt ; j++) {
254 *po += (*pi) * cosinus[j] ;
275 (
double* ti,
const int nt,
const int np,
276 const double tet,
double* to) {
284 double* sinus =
new double [nt] ;
285 for (j=0 ; j<nt ; j++) {
286 sinus[j] =
sin(j * tet) ;
291 for (j=0 ; j<nt ; j++) {
292 *po += (*pi) * sinus[j] ;
304 for (k=2 ; k<np+1 ; k++) {
306 for (j=0 ; j<nt ; j++) {
307 *po += (*pi) * sinus[j] ;
324 (
double* ti,
const int nt,
const int np,
325 const double tet,
double* to) {
333 double* sinus =
new double [nt] ;
334 for (j=0 ; j<nt-1 ; j++) {
335 sinus[j] =
sin(2*j * tet) ;
341 for (j=0 ; j<nt ; j++) {
342 *po += (*pi) * sinus[j] ;
354 for (k=2 ; k<np+1 ; k++) {
356 for (j=0 ; j<nt ; j++) {
357 *po += (*pi) * sinus[j] ;
374 (
double* ti,
const int nt,
const int np,
375 const double tet,
double* to) {
383 double* sinus =
new double [nt] ;
384 for (j=0 ; j<nt-1 ; j++) {
385 sinus[j] =
sin( (2*j+1) * tet) ;
391 for (j=0 ; j<nt ; j++) {
392 *po += (*pi) * sinus[j] ;
404 for (k=2 ; k<np+1 ; k++) {
406 for (j=0 ; j<nt ; j++) {
407 *po += (*pi) * sinus[j] ;
424void som_tet_cossin_cp
425 (
double* ti,
const int nt,
const int np,
426 const double tet,
double* to) {
434 double* cossin =
new double [2*nt] ;
435 for (j=0 ; j<2*nt ; j +=2) {
436 cossin[j] =
cos(j * tet) ;
437 cossin[j+1] =
sin((j+1) * tet) ;
442 for (j=0 ; j<nt ; j++) {
443 *po += (*pi) * cossin[2*j] ;
455 for (k=2 ; k<np+1 ; k++) {
458 for (j=0 ; j<nt ; j++) {
459 *po += (*pi) * cossin[2*j + m] ;
476void som_tet_cossin_ci
477 (
double* ti,
const int nt,
const int np,
478 const double tet,
double* to) {
486 double* cossin =
new double [2*nt] ;
487 for (j=0 ; j<2*nt ; j +=2) {
488 cossin[j] =
cos((j+1) * tet) ;
489 cossin[j+1] =
sin(j * tet) ;
494 for (j=0 ; j<nt ; j++) {
495 *po += (*pi) * cossin[2*j] ;
507 for (k=2 ; k<np+1 ; k++) {
510 for (j=0 ; j<nt ; j++) {
511 *po += (*pi) * cossin[2*j + m] ;
528void som_tet_cossin_sp
529 (
double* ti,
const int nt,
const int np,
530 const double tet,
double* to) {
538 double* cossin =
new double [2*nt] ;
539 for (j=0 ; j<2*nt ; j +=2) {
540 cossin[j] =
sin(j * tet) ;
541 cossin[j+1] =
cos((j+1) * tet) ;
546 for (j=0 ; j<nt ; j++) {
547 *po += (*pi) * cossin[2*j] ;
559 for (k=2 ; k<np+1 ; k++) {
562 for (j=0 ; j<nt ; j++) {
563 *po += (*pi) * cossin[2*j + m] ;
580void som_tet_cossin_si
581 (
double* ti,
const int nt,
const int np,
582 const double tet,
double* to) {
590 double* cossin =
new double [2*nt] ;
591 for (j=0 ; j<2*nt ; j +=2) {
592 cossin[j] =
sin((j+1) * tet) ;
593 cossin[j+1] =
cos(j * tet) ;
598 for (j=0 ; j<nt ; j++) {
599 *po += (*pi) * cossin[2*j] ;
611 for (k=2 ; k<np+1 ; k++) {
614 for (j=0 ; j<nt ; j++) {
615 *po += (*pi) * cossin[2*j + m] ;
632 (
double* ti,
const int nt,
const int np,
633 const double tet,
double* to) {
641 double* cossin =
new double [2*nt] ;
642 for (j=0 ; j<2*nt ; j +=2) {
643 cossin[j] =
cos(j/2 * tet) ;
644 cossin[j+1] =
sin(j/2 * tet) ;
649 for (j=0 ; j<nt ; j++) {
650 *po += (*pi) * cossin[2*j] ;
662 for (k=2 ; k<np+1 ; k++) {
665 for (j=0 ; j<nt ; j++) {
666 *po += (*pi) * cossin[2*j + m] ;
683 (
double* ti,
const int nt,
const int np,
684 const double tet,
double* to) {
692 double* cossin =
new double [2*nt] ;
693 for (j=0 ; j<2*nt ; j +=2) {
694 cossin[j] =
sin(j/2 * tet) ;
695 cossin[j+1] =
cos(j/2 * tet) ;
700 for (j=0 ; j<nt ; j++) {
701 *po += (*pi) * cossin[2*j] ;
713 for (k=2 ; k<np+1 ; k++) {
716 for (j=0 ; j<nt ; j++) {
717 *po += (*pi) * cossin[2*j + m] ;
Cmp sin(const Cmp &)
Sine.
Cmp cos(const Cmp &)
Cosine.