23char legendre_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Non_class_members/Coef/legendre.C,v 1.6 2014/10/13 08:53:13 j_novak Exp $" ;
95double* legendre(
int m,
int nt) {
103 double dt = M_PI / double(2*(nt-1)) ;
108 double* resu =
new double[(lmax-m+1)*nt] ;
111 double* cost =
new double[nt] ;
118 for (j=0; j<nt; j++) {
128 for (j=0; j<nt; j++) {
130 double s =
sin(j*dt) ;
131 for (i=1 ; i<2*m; i+=2) {
152 for (j=0; j<nt; j++) {
153 cost[j] =
cos(j*dt) ;
156 for (j=0; j<nt; j++) {
157 resu[nt+j] = cost[j] * (2.*m+1) * resu[j] ;
164 for (l=m+2; l < lmax+1 ; l++) {
166 int i_lm1 = nt*(l-1-m) ;
167 int i_lm2 = nt*(l-2-m) ;
172 for (j=0; j<nt; j++) {
173 resu[i_l+j] = ( cost[j] * a * resu[i_lm1+j]
174 - b * resu[i_lm2+j] ) / c ;
Cmp sin(const Cmp &)
Sine.
Cmp cos(const Cmp &)
Cosine.