95 #include "param_elliptic.h" 106 assert(source.
get_etat() != ETATNONDEF) ;
117 if (sourva.get_etat() == ETATZERO) {
123 assert(sourva.get_etat() == ETATQCQ) ;
125 Valeur rho(sourva.get_mg()) ;
127 rho = *(sourva.c_cf) ;
140 Mtbl_cf resu = elliptic_solver (ope_var, *(rho.c_cf)) ;
164 double fact_dir,
double fact_neu )
const {
166 assert(source.
get_etat() != ETATNONDEF) ;
177 if (sourva.get_etat() == ETATZERO) {
183 assert(sourva.get_etat() == ETATQCQ) ;
185 Valeur rho(sourva.get_mg()) ;
187 rho = *(sourva.c_cf) ;
200 Mtbl_cf resu = elliptic_solver_boundary (ope_var, *(rho.c_cf), bound,
201 fact_dir, fact_neu) ;
227 double fact_dir,
double fact_neu )
const {
229 assert(source.
get_etat() != ETATNONDEF) ;
240 if (sourva.get_etat() == ETATZERO) {
246 assert(sourva.get_etat() == ETATQCQ) ;
248 Valeur rho(sourva.get_mg()) ;
250 rho = *(sourva.c_cf) ;
279 int nr = gri2d.
get_nr(0) ;
280 int nt = gri2d.
get_nt(0) ;
281 int np = gri2d.
get_np(0) ;
283 for(
int k=0; k<np+2; k++)
284 for (
int j=0; j<=nt-1; j++)
285 for(
int xi=0; xi<= nr-1; xi++)
291 Mtbl_cf resu = elliptic_solver_boundary (ope_var, *(rho.c_cf), bound2,
292 fact_dir, fact_neu) ;
316 Scalar& pot,
double val)
const {
318 assert(source.
get_etat() != ETATNONDEF) ;
329 if (sourva.get_etat() == ETATZERO) {
335 assert(sourva.get_etat() == ETATQCQ) ;
337 Valeur rho(sourva.get_mg()) ;
339 rho = *(sourva.c_cf) ;
352 Mtbl_cf resu = elliptic_solver_no_zec (ope_var, *(rho.c_cf), val) ;
373 Scalar& pot,
double val)
const {
375 assert(source.
get_etat() != ETATNONDEF) ;
386 if (sourva.get_etat() == ETATZERO) {
392 assert(sourva.get_etat() == ETATQCQ) ;
394 Valeur rho(sourva.get_mg()) ;
396 rho = *(sourva.c_cf) ;
409 Mtbl_cf resu = elliptic_solver_only_zec (ope_var, *(rho.c_cf), val) ;
430 const Scalar& source,
Scalar& pot,
double* amplis,
double* phases)
const {
432 assert(source.
get_etat() != ETATNONDEF) ;
443 if (sourva.get_etat() == ETATZERO) {
449 assert(sourva.get_etat() == ETATQCQ) ;
451 Valeur rho(sourva.get_mg()) ;
453 rho = *(sourva.c_cf) ;
466 Mtbl_cf resu = elliptic_solver_sin_zec (ope_var, *(rho.c_cf), amplis, phases) ;
491 assert(source.
get_etat() != ETATNONDEF) ;
502 if (sourva.get_etat() == ETATZERO) {
508 assert(sourva.get_etat() == ETATQCQ) ;
510 Valeur rho(sourva.get_mg()) ;
512 rho = *(sourva.c_cf) ;
526 Mtbl_cf resu = elliptic_solver_fixe_der_zero (valeur, ope_var, *(rho.c_cf)) ;
const Base_val & get_spectral_base() const
Returns the spectral bases of the Valeur va.
Mtbl_cf * c_cf
Coefficients of the spectral expansion of the function.
void ylm_i()
Inverse of ylm()
int get_np(int l) const
Returns the number of points in the azimuthal direction ( ) in domain no. l.
void coef() const
Computes the coeffcients of *this.
virtual void set_etat_zero()
Sets the logical state to ETATZERO (zero).
double * alpha
Array (size: mg->nzone ) of the values of in each domain.
void sol_elliptic_only_zec(Param_elliptic ¶ms, const Scalar &so, Scalar &uu, double val) const
General elliptic solver.
void ylm()
Computes the coefficients of *this.
const Mg3d * get_mg() const
Gives the Mg3d on which the mapping is defined.
Tensor field of valence 0 (or component of a tensorial field).
Base class for coordinate mappings.
Values and coefficients of a (real-value) function.
int get_etat() const
Returns the logical state ETATNONDEF (undefined), ETATZERO (null) or ETATQCQ (ordinary).
virtual void set_etat_qcq()
Sets the logical state to ETATQCQ (ordinary state).
void set_dzpuis(int)
Modifies the dzpuis flag.
Scalar var_F
Additive variable change function.
void sol_elliptic_no_zec(Param_elliptic ¶ms, const Scalar &so, Scalar &uu, double val) const
General elliptic solver.
Scalar var_G
Multiplicative variable change that must be sphericaly symetric !
This class contains the parameters needed to call the general elliptic solver.
void sol_elliptic(Param_elliptic ¶ms, const Scalar &so, Scalar &uu) const
General elliptic solver.
int get_nr(int l) const
Returns the number of points in the radial direction ( ) in domain no. l.
const Mg3d * mg
Pointer on the multi-grid Mgd3 on which this is defined.
void annule_hard()
Sets the Mtbl_cf to zero in a hard way.
void sol_elliptic_fixe_der_zero(double val, Param_elliptic ¶ms, const Scalar &so, Scalar &uu) const
General elliptic solver fixing the derivative at the origin and relaxing the continuity of the first ...
Coefficients storage for the multi-domain spectral method.
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.
int get_nt(int l) const
Returns the number of points in the co-latitude direction ( ) in domain no. l.
Valeur & set_spectral_va()
Returns va (read/write version)
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...
const Map & get_mp() const
Returns the mapping.
const Mg3d * get_angu_1dom() const
Returns the pointer on the associated mono-domain angular grid.
void sol_elliptic_sin_zec(Param_elliptic ¶ms, const Scalar &so, Scalar &uu, double *coefs, double *) const
General elliptic solver.
const Valeur & get_spectral_va() const
Returns va (read only version)