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 ;
double alpha
Parameter of the associated mapping.
Matrice * ope_mat
Pointer on the matrix representation of the operator.
int base_r
Radial basis of decomposition.
#define TRA_R
Translation en R, used for a bitwise shift (in hex)
double b_param
The parameter .
double c_param
The parameter .
double a_param
The parameter .
virtual void do_ope_mat() const
Computes the matrix of the operator.
int nr
Number of radial points.
#define MAX_BASE
Nombre max. de bases differentes.
#define R_CHEB
base de Chebychev ordinaire (fin)