23 char prepa_helmholtz_minus_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Non_class_members/PDE/prepa_helmholtz_minus.C,v 1.8 2014/10/13 08:53:30 j_novak Exp $" ;
68 #include "type_parite.h"
79 Matrice _prepa_helmholtz_minus_nondege_pas_prevu(
const Matrice &so) {
81 cout <<
"Unknown case for prepa_helmholtz_minus_nondege" << endl ;
91 Matrice _prepa_helmholtz_minus_nondege_r_cheb (
const Matrice &lap) {
96 Matrice res(n-non_dege, n-non_dege) ;
98 for (
int i=0 ; i<n-non_dege ; i++)
99 for (
int j=0 ; j<n-non_dege ; j++)
100 res.set(i, j) = lap(i, j+non_dege) ;
110 Matrice _prepa_helmholtz_minus_nondege_r_chebu (
const Matrice &lap) {
115 Matrice res(n-non_dege, n-non_dege) ;
117 for (
int i=0 ; i<n-non_dege ; i++)
118 for (
int j=0 ; j<n-non_dege ; j++)
119 res.set(i, j) = lap(i, j+non_dege) ;
128 Matrice _prepa_helmholtz_minus_nondege_r_chebp (
const Matrice &lap) {
133 Matrice res(n-non_dege, n-non_dege) ;
135 for (
int i=0 ; i<n-non_dege ; i++)
136 for (
int j=0 ; j<n-non_dege ; j++)
137 res.set(i, j) = lap(i, j+non_dege) ;
146 Matrice _prepa_helmholtz_minus_nondege_r_chebi (
const Matrice &lap) {
151 Matrice res(n-non_dege, n-non_dege) ;
153 for (
int i=0 ; i<n-non_dege ; i++)
154 for (
int j=0 ; j<n-non_dege ; j++)
155 res.set(i, j) = lap(i, j+non_dege) ;
166 Matrice prepa_helmholtz_minus_nondege(
const Matrice &ope,
int base_r) {
169 static Matrice (*prepa_helmholtz_minus_nondege[
MAX_BASE])
177 prepa_helmholtz_minus_nondege[i] =
178 _prepa_helmholtz_minus_nondege_pas_prevu ;
182 _prepa_helmholtz_minus_nondege_r_cheb ;
184 _prepa_helmholtz_minus_nondege_r_chebu ;
186 _prepa_helmholtz_minus_nondege_r_chebp ;
188 _prepa_helmholtz_minus_nondege_r_chebi ;
191 Matrice res(prepa_helmholtz_minus_nondege[base_r](ope)) ;