23char cirleg_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Non_class_members/Coef/cirleg.C,v 1.3 2014/10/13 08:53:11 j_novak Exp $" ;
57#include "utilitaires.h"
63void cirleg(
const int* deg,
const int* dimc,
double* cf,
const int* dimf,
81 cout <<
"cirleg: nr > n3c : nr = " << nr <<
" , n3c = "
87 cout <<
"cirleg: nr > n3f : nr = " << nr <<
" , n3f = "
93 cout <<
"cirleg: n1c > n1f : n1c = " << n1c <<
" , n1f = "
99 cout <<
"cirleg: n2c > n2f : n2c = " << n2c <<
" , n2f = "
107 int n2n3f = n2f * n3f ;
108 int n2n3c = n2c * n3c ;
116 int borne_phi = ( n1c > 1 ) ? n1c-1 : 1 ;
118 double* colloc =
new double[nr] ;
119 legendre_collocation_points(nr, colloc) ;
121 for (
int j=0; j< borne_phi; j++) {
125 for (
int k=0; k<n2c; k++) {
127 int i0 = n2n3c * j + n3c * k ;
128 double* cf0 = cf + i0 ;
130 i0 = n2n3f * j + n3f * k ;
131 double* ff0 = ff + i0 ;
133 for (
int i = 0; i<nr; i++) {
134 double x0 = colloc[i] ;
137 double som = cf0[0] + cf0[1]*x0 ;
138 for (
int h=2; h<nr; h++) {
139 double Pip2 = (2. - 1./double(h))*x0*Pip1
140 - (1. - 1./
double(h))*Pi ;
154void cirlegp(
const int* deg,
const int* dimc,
double* cf,
155 const int* dimf,
double* ff)
172 cout <<
"cirlegp: nr > n3c : nr = " << nr <<
" , n3c = "
178 cout <<
"cirlegp: nr > n3f : nr = " << nr <<
" , n3f = "
184 cout <<
"cirlegp: n1c > n1f : n1c = " << n1c <<
" , n1f = "
190 cout <<
"cirlegp: n2c > n2f : n2c = " << n2c <<
" , n2f = "
198 int dnm1 = 2*nr - 1 ;
202 int n2n3f = n2f * n3f ;
203 int n2n3c = n2c * n3c ;
211 int borne_phi = ( n1c > 1 ) ? n1c-1 : 1 ;
213 double* colloc =
new double[dnm1] ;
214 legendre_collocation_points(dnm1, colloc) ;
216 for (
int j=0; j< borne_phi; j++) {
220 for (
int k=0; k<n2c; k++) {
222 int i0 = n2n3c * j + n3c * k ;
223 double* cf0 = cf + i0 ;
225 i0 = n2n3f * j + n3f * k ;
226 double* ff0 = ff + i0 ;
228 for (
int i = 0; i<nr; i++) {
229 double x0 = colloc[nm1+i] ;
232 double som = cf0[0] ;
233 for (
int h=2; h<dnm1; h++) {
234 double Pip2 = (2. - 1./double(h))*x0*Pip1
235 - (1. - 1./
double(h))*Pi ;
236 if (h%2 == 0) som += cf0[h/2]*Pip2 ;
252void cirlegi(
const int* deg,
const int* dimc,
double* cf,
253 const int* dimf,
double* ff)
270 cout <<
"cirlegi: nr > n3c : nr = " << nr <<
" , n3c = "
276 cout <<
"cirlegi: nr > n3f : nr = " << nr <<
" , n3f = "
282 cout <<
"cirlegi: n1c > n1f : n1c = " << n1c <<
" , n1f = "
288 cout <<
"cirlegi: n2c > n2f : n2c = " << n2c <<
" , n2f = "
296 int dnm1 = 2*nr - 1 ;
300 int n2n3f = n2f * n3f ;
301 int n2n3c = n2c * n3c ;
309 int borne_phi = ( n1c > 1 ) ? n1c-1 : 1 ;
311 double* colloc =
new double[dnm1] ;
312 legendre_collocation_points(dnm1, colloc) ;
314 for (
int j=0; j< borne_phi; j++) {
318 for (
int k=0; k<n2c; k++) {
320 int i0 = n2n3c * j + n3c * k ;
321 double* cf0 = cf + i0 ;
323 i0 = n2n3f * j + n3f * k ;
324 double* ff0 = ff + i0 ;
326 for (
int i = 0; i<nr; i++) {
327 double x0 = colloc[nm1+i] ;
330 double som = cf0[0]*x0 ;
331 for (
int h=2; h<dnm1; h++) {
332 double Pip2 = (2. - 1./double(h))*x0*Pip1
333 - (1. - 1./
double(h))*Pi ;
334 if (h%2 == 1) som += cf0[h/2]*Pip2 ;