jacobipointsgl.C
#include "tbl.h"
namespace Lorene {
double* jacobi(int , double) ;
double* pointsgausslobatto(int) ;
Tbl jacobipointsgl(int n) {
double* pointsgl = pointsgausslobatto(n) ;
Tbl jj(n+1,n+1) ;
jj.set_etat_qcq() ;
int i,k ;
for (i = 0 ; i < n+1 ; i++ ) {
double* yy = jacobi(n,pointsgl[i]) ;
for (k = 0 ; k < n+1 ; k++ ) {
jj.set(k,i) = yy[k] ;
}
delete [] yy ;
}
delete [] pointsgl ;
return jj ;
}
}
Lorene::Tbl
Basic array class.
Definition tbl.h:161
Lorene
Lorene prototypes.
Definition app_hor.h:64