32 #include "ope_elementary.h" 40 Matrice _poisson_pseudo_1d_non_dege_pas_prevu(
const Matrice &lap,
int) {
41 cout <<
"Construction non degeneree pas prevue..." << endl ;
53 Matrice _poisson_pseudo_1d_non_dege_r_cheb (
const Matrice &lap,
int) {
58 Matrice res(n-2, n-2) ;
60 for (
int i=0 ; i<n-2 ; i++)
61 for (
int j=0 ; j<n-2 ; j++)
62 res.set(i, j) = lap(i, j+2) ;
77 Matrice _poisson_pseudo_1d_non_dege_r_chebp (
const Matrice &lap,
int l) {
80 assert (div(l, 2).rem == 0) ;
84 Matrice res(n-1, n-1) ;
86 for (
int i=0 ; i<n-1 ; i++)
87 for (
int j=0 ; j<n-1 ; j++)
88 res.set(i, j) = lap(i, j+1) ;
95 Matrice res(n-2, n-2) ;
97 for (
int i=0 ;i<n-2 ; i++)
98 for (
int j=0 ; j<n-2 ; j++)
99 res.set(i, j) = lap(i, j+2) ;
113 Matrice _poisson_pseudo_1d_non_dege_r_chebi (
const Matrice &lap,
int l) {
117 assert (div(l, 2).rem == 1) ;
120 Matrice res(n-1, n-1) ;
122 for (
int i=0 ; i<n-1 ; i++)
123 for (
int j=0 ; j<n-1 ; j++)
124 res.set(i, j) = lap(i, j+1) ;
131 Matrice res(n-2, n-2) ;
133 for (
int i=0 ;i<n-2 ; i++)
134 for (
int j=0 ; j<n-2 ; j++)
135 res.set(i, j) = lap(i, j+2) ;
160 poisson_pseudo_1d_non_dege[i] = _poisson_pseudo_1d_non_dege_pas_prevu ;
164 _poisson_pseudo_1d_non_dege_r_cheb ;
166 _poisson_pseudo_1d_non_dege_r_chebp ;
168 _poisson_pseudo_1d_non_dege_r_chebi ;
virtual void do_non_dege() const
Computes the non-degenerated matrix of the operator.
Matrice * ope_cl
Pointer on the banded-matrix of the operator.
virtual void do_ope_cl() const
Computes the banded-matrix of the operator.
int base_r
Radial basis of decomposition.
#define TRA_R
Translation en R, used for a bitwise shift (in hex)
#define R_CHEBI
base de Cheb. impaire (rare) seulement
#define R_CHEBP
base de Cheb. paire (rare) seulement
int get_dim(int i) const
Returns the dimension of the matrix.
#define MAX_BASE
Nombre max. de bases differentes.
Matrice * non_dege
Pointer on the non-degenerated matrix of the operator.
int l_quant
quantum number
#define R_CHEB
base de Chebychev ordinaire (fin)