24char trigo_ini_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Non_class_members/Coef/FFT991/trigo_ini.C,v 1.4 2014/10/15 12:48:22 j_novak Exp $" ;
93double *trigo_ini(
int n )
98static double *table_trigo[NMAX] ;
100static int tbn[NMAX] ;
101static int trois = 3 ;
108 for ( i=0 ; i < nwork ; i++ ) {
109 if ( tbn[i] == n ) indice = i ;
114 if ( nwork >= NMAX ) {
115 cout <<
"trigo_ini : nwork >= NMAX !" << endl ;
118 indice = nwork ; nwork++ ; tbn[indice] = n ;
121 table_trigo[indice] = (
double *) malloc(
sizeof(
double) * (3*n/2 + 1) ) ;
122 if ( table_trigo[indice] == 0 ) {
123 cout <<
"trigo_ini : malloc error !" << endl ;
127 F77_fftrig( table_trigo[indice], &n, &trois ) ;
133 return table_trigo[indice] ;