21char 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"
40Matrice _sec_order_mat_pas_prevu(
int,
double,
double,
double,
double) {
41 cout <<
"Sec_order : base not implemented..." << endl ;
54Matrice _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 ;
Time evolution with partial storage (*** under development ***).
Matrice * ope_mat
Pointer on the matrix representation of the operator.
double alpha
Parameter of the associated mapping.
int base_r
Radial basis of decomposition.
int nr
Number of radial points.
double b_param
The parameter .
virtual void do_ope_mat() const
Computes the matrix of the operator.
double a_param
The parameter .
double c_param
The parameter .
#define MAX_BASE
Nombre max. de bases differentes.
#define TRA_R
Translation en R, used for a bitwise shift (in hex)
#define R_CHEB
base de Chebychev ordinaire (fin)