File: jacobipointsgl.C

package info (click to toggle)
lorene 0.0.0~cvs20161116%2Bdfsg-2
  • links: PTS, VCS
  • area: main
  • in suites: forky, sid, trixie
  • size: 26,472 kB
  • sloc: cpp: 212,946; fortran: 21,645; makefile: 1,750; sh: 4
file content (26 lines) | stat: -rw-r--r-- 429 bytes parent folder | download | duplicates (2) 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
 
#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 ;
}
}