24char op_d2sdphi2_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Non_class_members/Operators/op_d2sdphi2.C,v 1.3 2014/10/13 08:53:24 j_novak Exp $" ;
76void _d2sdphi2_pas_prevu(Tbl* ,
int & b) {
77 cout <<
"Unknown phi basis in Mtbl_cf::d2sdp2() !" << endl ;
78 cout <<
" basis: " << hex << b << endl ;
84void _d2sdphi2_p_cossin(Tbl* tb,
int & )
88 if (tb->get_etat() == ETATZERO) {
93 assert(tb->get_etat() == ETATQCQ) ;
96 int nr = (tb->dim).dim[0] ;
97 int nt = (tb->dim).dim[1] ;
98 int np = (tb->dim).dim[2] ;
102 static double* cx = 0 ;
103 static int np_pre =0 ;
108 cx =
reinterpret_cast<double*
>(realloc(cx, (np+2) *
sizeof(
double))) ;
109 for (
int i=0 ; i<np+2 ; i++) {
110 cx[i] = - (i/2) * (i/2) ;
115 double* xo =
new double[(tb->dim).taille] ;
118 for (
int i=0; i<(tb->dim).taille; i++) {
128 for (
int j=0 ; j<nt ; j++) {
129 for (
int i=0 ; i<nr ; i++ ) {
130 *xco = cx[0] * (*xci) ;
141 for (
int k=2 ; k<np+1 ; k++) {
142 for (
int j=0 ; j<nt ; j++) {
143 for (
int i=0 ; i<nr ; i++ ) {
144 *xco = cx[k] * (*xci) ;