21 char ope_sec_order_mat_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Ope_elementary/Ope_sec_order/ope_sec_order_mat.C,v 1.3 2014/10/13 08:53:35 j_novak Exp $" ;
32 #include "ope_elementary.h"
40 Matrice _sec_order_mat_pas_prevu(
int,
double,
double,
double,
double) {
41 cout <<
"Sec_order : base not implemented..." << endl ;
54 Matrice _sec_order_mat_r_cheb (
int n,
double alpha,
55 double a,
double b,
double c) {
57 double* vect =
new double[n] ;
68 for (
int i=0 ; i<n ; i++) {
69 for (
int j=0 ; j<n ; j++)
73 d2sdx2_1d (n, &vect,
R_CHEB) ;
75 for (
int j=0 ; j<n ; j++)
76 dd.set(j,i) = vect[j] ;
80 for (
int i=0 ; i<n ; i++) {
81 for (
int j=0 ; j<n ; j++)
85 sxdsdx_1d (n, &vect,
R_CHEB) ;
87 for (
int j=0 ; j<n ; j++)
88 df.set(j,i) = vect[j] ;
92 for (
int i=0 ; i<n ; i++) {
93 for (
int j=0 ; j<n ; j++)
100 return a*dd/alpha/alpha+b*df/alpha+c*ff ;
109 double, double, double);
116 sec_order_mat[i] = _sec_order_mat_pas_prevu ;
119 sec_order_mat[
R_CHEB >>
TRA_R] = _sec_order_mat_r_cheb ;