21char param_elliptic_pseudo_1d_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Param_elliptic/param_elliptic_pseudo_1d.C,v 1.3 2014/10/13 08:53:37 j_novak Exp $" ;
50#include "param_elliptic.h"
53#include "ope_elementary.h"
54#include "change_var.h"
62 cout <<
"set_poisson_pseudo_1d only defined for an affine mapping..." <<
endl ;
67 int nz =
get_mp().get_mg()->get_nzone() ;
74 for (
int l=0 ;
l<nz ;
l++) {
76 nr =
get_mp().get_mg()->get_nr(
l) ;
77 alpha = get_alpha (
l) ;
80 for (
int k=0 ;
k<
get_mp().get_mg()->get_np(
l)+1 ;
k++)
81 for (
int j=0 ;
j<
get_mp().get_mg()->get_nt(
l) ;
j++) {
85 if ((
k!=1) && (
l!=nz-1))
97 int dzpuis =
source.get_dzpuis() ;
101 cout <<
"set_helmholtz_minus_pseudo_1d only defined for an affine mapping..." <<
endl ;
106 int nz =
get_mp().get_mg()->get_nzone() ;
114 for (
int l=0 ;
l<nz ;
l++) {
116 nr =
get_mp().get_mg()->get_nr(
l) ;
117 alpha = get_alpha (
l) ;
118 beta = get_beta (
l) ;
120 for (
int k=0 ;
k<
get_mp().get_mg()->get_np(
l)+1 ;
k++)
121 for (
int j=0 ;
j<
get_mp().get_mg()->get_nt(
l) ;
j++) {
125 source.get_spectral_va().base.give_quant_numbers
128 alpha, beta, l_quant, masse, dzpuis) ;
Time evolution with partial storage (*** under development ***).
Class for the operator of the modified Helmholtz equation in pseudo-1d.
Class for the operator of the Poisson equation in pseudo 1d.
const Map_radial & get_mp() const
Returns the mapping.
void set_poisson_pseudo_1d(Scalar &so)
Set the operator to everywhere but in the compactified domain.
void set_helmholtz_minus_pseudo_1d(int zone, double mas, Scalar &so)
Set the operator to in one domain.
Ope_elementary ** operateurs
Array on the elementary operators.
int type_map
Type of mapping either MAP_AFF or MAP_LOG.
Tensor field of valence 0 (or component of a tensorial field).