23char prepa_helmholtz_plus_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Non_class_members/PDE/prepa_helmholtz_plus.C,v 1.5 2014/10/13 08:53:30 j_novak Exp $" ;
54#include "type_parite.h"
65Matrice _prepa_helmholtz_plus_nondege_pas_prevu(
const Matrice &so) {
67 cout <<
"Unknown case for prepa_helmholtz_plus_nondege" << endl ;
78Matrice _prepa_helmholtz_plus_nondege_r_cheb (
const Matrice &lap) {
83 Matrice res(n-non_dege, n-non_dege) ;
85 for (
int i=0 ; i<n-non_dege ; i++)
86 for (
int j=0 ; j<n-non_dege ; j++)
87 res.set(i, j) = lap(i, j+non_dege) ;
100Matrice _prepa_helmholtz_plus_nondege_r_chebp (
const Matrice &lap) {
105 Matrice res(n-non_dege, n-non_dege) ;
107 for (
int i=0 ; i<n-non_dege ; i++)
108 for (
int j=0 ; j<n-non_dege ; j++)
109 res.set(i, j) = lap(i, j+non_dege) ;
121Matrice prepa_helmholtz_plus_nondege(
const Matrice &ope,
int base_r) {
124 static Matrice (*prepa_helmholtz_plus_nondege[
MAX_BASE])
132 prepa_helmholtz_plus_nondege[i] =
133 _prepa_helmholtz_plus_nondege_pas_prevu ;
137 _prepa_helmholtz_plus_nondege_r_cheb ;
139 _prepa_helmholtz_plus_nondege_r_chebp ;
142 Matrice res(prepa_helmholtz_plus_nondege[base_r](ope)) ;
int get_dim(int i) const
Returns the dimension of the matrix.
#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)
#define R_CHEBP
base de Cheb. paire (rare) seulement