23char map_log_elliptic_C[] =
"$Header $" ;
68#include "param_elliptic.h"
78 assert(source.
get_etat() != ETATNONDEF) ;
95 assert(sourva.
get_etat() == ETATQCQ) ;
99 rho = *(sourva.
c_cf) ;
112 Mtbl_cf resu = elliptic_solver (ope_var, *(rho.
c_cf)) ;
134 double fact_dir,
double fact_neu)
const {
136 assert(source.
get_etat() != ETATNONDEF) ;
147 if (sourva.
get_etat() == ETATZERO) {
153 assert(sourva.
get_etat() == ETATQCQ) ;
157 rho = *(sourva.
c_cf) ;
170 Mtbl_cf resu = elliptic_solver_boundary (ope_var, *(rho.
c_cf), bound,
171 fact_dir, fact_neu) ;
193 double fact_dir,
double fact_neu)
const {
195 assert(source.
get_etat() != ETATNONDEF) ;
206 if (sourva.
get_etat() == ETATZERO) {
212 assert(sourva.
get_etat() == ETATQCQ) ;
216 rho = *(sourva.
c_cf) ;
242 int nr = gri2d.
get_nr(0) ;
243 int nt = gri2d.
get_nt(0) ;
244 int np = gri2d.
get_np(0) ;
246 for(
int k=0; k<np+2; k++)
247 for (
int j=0; j<=nt-1; j++)
248 for(
int xi=0; xi<= nr-1; xi++)
255 Mtbl_cf resu = elliptic_solver_boundary (ope_var, *(rho.
c_cf), bound2,
256 fact_dir, fact_neu) ;
276 Scalar& pot,
double val)
const {
278 assert(source.
get_etat() != ETATNONDEF) ;
289 if (sourva.
get_etat() == ETATZERO) {
295 assert(sourva.
get_etat() == ETATQCQ) ;
299 rho = *(sourva.
c_cf) ;
312 Mtbl_cf resu = elliptic_solver_no_zec (ope_var, *(rho.
c_cf), val) ;
void sol_elliptic_boundary(Param_elliptic ¶ms, const Scalar &so, Scalar &uu, const Mtbl_cf &bound, double fact_dir, double fact_neu) const
General elliptic solver including inner boundary conditions.
void sol_elliptic(Param_elliptic ¶ms, const Scalar &so, Scalar &uu) const
General elliptic solver.
void sol_elliptic_no_zec(Param_elliptic ¶ms, const Scalar &so, Scalar &uu, double) const
General elliptic solver.
Base class for coordinate mappings.
const Mg3d * get_mg() const
Gives the Mg3d on which the mapping is defined.
const Mg3d * mg
Pointer on the multi-grid Mgd3 on which this is defined
const Mg3d * get_angu_1dom() const
Returns the pointer on the associated mono-domain angular grid.
int get_np(int l) const
Returns the number of points in the azimuthal direction ( ) in domain no. l.
int get_nt(int l) const
Returns the number of points in the co-latitude direction ( ) in domain no. l.
int get_nr(int l) const
Returns the number of points in the radial direction ( ) in domain no. l.
Coefficients storage for the multi-domain spectral method.
Tbl & set(int l)
Read/write of the Tbl containing the coefficients in a given domain.
void annule_hard()
Sets the Mtbl_cf to zero in a hard way.
This class contains the parameters needed to call the general elliptic solver.
Scalar var_G
Multiplicative variable change that must be sphericaly symetric !
Scalar var_F
Additive variable change function.
Tensor field of valence 0 (or component of a tensorial field).
virtual void set_etat_qcq()
Sets the logical state to ETATQCQ (ordinary state).
bool check_dzpuis(int dzi) const
Returns false if the last domain is compactified and *this is not zero in this domain and dzpuis is n...
virtual void set_etat_zero()
Sets the logical state to ETATZERO (zero).
Valeur & set_spectral_va()
Returns va (read/write version)
const Valeur & get_spectral_va() const
Returns va (read only version)
int get_etat() const
Returns the logical state ETATNONDEF (undefined), ETATZERO (null) or ETATQCQ (ordinary).
const Base_val & get_spectral_base() const
Returns the spectral bases of the Valeur va
void set_dzpuis(int)
Modifies the dzpuis flag.
Values and coefficients of a (real-value) function.
int get_etat() const
Returns the logical state.
void ylm()
Computes the coefficients of *this.
Mtbl_cf * c_cf
Coefficients of the spectral expansion of the function.
void coef() const
Computes the coeffcients of *this.
const Mg3d * get_mg() const
Returns the Mg3d on which the this is defined.
void ylm_i()
Inverse of ylm()
const Map & get_mp() const
Returns the mapping.