21 char ope_poisson_pseudo_1d_non_dege_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Ope_elementary/Ope_poisson_pseudo_1d/ope_poisson_pseudo_1d_non_dege.C,v 1.3 2014/10/13 08:53:35 j_novak Exp $" ;
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 ;