LORENE
jacobipointsgl.C
1#include "tbl.h"
2
3namespace Lorene {
4double* jacobi(int , double) ;
5double* pointsgausslobatto(int) ;
6
7Tbl jacobipointsgl(int n) {
8
9 double* pointsgl = pointsgausslobatto(n) ;
10
11 Tbl jj(n+1,n+1) ;
12 jj.set_etat_qcq() ;
13
14 int i,k ;
15
16 for (i = 0 ; i < n+1 ; i++ ) {
17 double* yy = jacobi(n,pointsgl[i]) ;
18 for (k = 0 ; k < n+1 ; k++ ) {
19 jj.set(k,i) = yy[k] ;
20 }
21 delete [] yy ;
22 }
23 delete [] pointsgl ;
24 return jj ;
25}
26}
Lorene prototypes.
Definition app_hor.h:64