Param_elliptic Class Reference
[General PDE solvers (under development)]

This class contains the parameters needed to call the general elliptic solver. More...

#include <param_elliptic.h>

Public Member Functions
	Param_elliptic (const Scalar &)
	Standard constructor from a `Scalar`.
	~Param_elliptic ()
	Destructor.
const Map_radial &	get_mp () const
	Returns the mapping.
double	get_alpha (int) const
double	get_beta (int) const
int	get_type (int) const
void	set_helmholtz_minus (int zone, double mas, Scalar &so)
	Set the operator to $\left(\Delta - m^2\right)$ in one domain (not in the nucleus).
void	set_poisson_pseudo_1d (Scalar &so)
	Set the operator to $\left(\Delta - 2 \partial_r / r\right)$ everywhere but in the compactified domain.
void	set_helmholtz_minus_pseudo_1d (int zone, double mas, Scalar &so)
	Set the operator to $\left(\Delta - 2 \partial_r / r - m^2\right)$ in one domain.
void	set_helmholtz_plus (int zone, double mas, Scalar &so)
	Set the operator to $\left(\Delta + m^2\right)$ in one domain (only in the shells).
void	set_poisson_2d (const Scalar &, bool indic=false)
	Set everything to do a 2d-Poisson, with or without l=0 (not put by default.
void	set_helmholtz_minus_2d (int zone, double mas, const Scalar &)
	Set the 2D Helmholtz operator (with minus sign).
void	set_sec_order_r2 (int zone, double a, double b, double c)
	Set the operator to $a r^2 \partial^2/\partial r^2 + b r \partial /\partial r + c$ in one domain (only in the shells).
void	set_sec_order (int zone, double a, double b, double c)
	Set the operator to $a \partial^2/\partial r^2 + b \partial /\partial r + c$ in one domain (only in the shells).
void	set_ope_vorton (int zone, Scalar &so)
	Set the operator to $\Delta - 2\partial /\partial r$ in one domain (not implemented in the nucleus).
void	set_poisson_vect_r (int zone, bool only_l_zero=false)
	Sets the operator to $\Delta + \frac{2}{r} \frac{\partial}{\partial r} + \frac{2 - l(l+1)}{r^2}$ in all domains, for $l \not= 0$ ; and to $\frac{\partial^2}{\partial r^2} + \frac{2}{r} \frac{\partial}{\partial r} - \frac{2}{r^2}$ in all domains otherwise.
void	set_poisson_vect_eta (int zone)
	Sets the operator to be a regular elliptic operator to solve for the $\eta$ component of the vector Poisson equation.
void	set_poisson_tens_rr (int zone)
	Sets the operator to $\Delta + \frac{4}{r} \frac{\partial}{\partial r} + \frac{6 - l(l+1)}{r^2}$ in all domains, for $l \geq 2$ only.
void	inc_l_quant (int zone)
	Increases the quantum number l in the domain `zone` .
void	set_variable_F (const Scalar &)
	Changes the variable function F.
void	set_variable_G (const Scalar &)
	Changes the variable function G.
double	F_plus (int zone, int k, int j) const
	Returns the value of F, for a given angular point, at the outer boundary of the domain `zone` ;.
double	F_minus (int zone, int k, int j) const
	Returns the value of F, for a given angular point, at the inner boundary of the domain `zone` ;.
double	dF_plus (int zone, int k, int j) const
	Returns the value of the radial derivative of F, for a given angular point, at the outer boundary of the domain `zone` ;.
double	dF_minus (int zone, int k, int j) const
	Returns the value of the radial derivative of F, for a given angular point, at the inner boundary of the domain `zone` ;.
double	G_plus (int zone) const
	Returns the value of G, for a given angular point, at the outer boundary of the domain `zone` ;.
double	G_minus (int zone) const
	Returns the value of G, for a given angular point, at the inner boundary of the domain `zone` ;.
double	dG_plus (int zone) const
	Returns the value of the radial derivative of G, for a given angular point, at the outer boundary of the domain `zone` ;.
double	dG_minus (int zone) const
	Returns the value of the radial derivative of G, for a given angular point, at the inner boundary of the domain `zone` ;.
Protected Attributes
int	type_map
	Type of mapping either MAP_AFF or MAP_LOG.
const Map_af *	mp_af
	The mapping, if affine.
const Map_log *	mp_log
	The mapping if log type.
Ope_elementary **	operateurs
	Array on the elementary operators.
Scalar	var_F
	Additive variable change function.
Scalar	var_G
	Multiplicative variable change that must be sphericaly symetric !
Itbl	done_F
	Stores what has been computed for `F`.
Itbl	done_G
	Stores what has been computed for `G`.
Tbl	val_F_plus
	Values of F at the outer boundaries of the various domains.
Tbl	val_F_minus
	Values of F at the inner boundaries of the various domains.
Tbl	val_dF_plus
	Values of the derivative of F at the outer boundaries of the various domains.
Tbl	val_dF_minus
	Values of the derivative of F at the inner boundaries of the various domains.
Tbl	val_G_plus
	Values of G at the outer boundaries of the various domains.
Tbl	val_G_minus
	Values of G at the inner boundaries of the various domains.
Tbl	val_dG_plus
	Values of the derivative of G at the outer boundaries of the various domains.
Tbl	val_dG_minus
	Values of the derivative of G at the inner boundaries of the various domains.
Private Member Functions
void	compute_val_F (int, int, int) const
	Computes the various values of `F`.
void	compute_val_G (int) const
	Computes the various values of `G`.
Friends
Mtbl_cf	elliptic_solver (const Param_elliptic &, const Mtbl_cf &)
Mtbl_cf	elliptic_solver_boundary (const Param_elliptic &, const Mtbl_cf &, const Mtbl_cf &bound, double fact_dir, double fact_neu)
Mtbl_cf	elliptic_solver_no_zec (const Param_elliptic &, const Mtbl_cf &, double)
Mtbl_cf	elliptic_solver_only_zec (const Param_elliptic &, const Mtbl_cf &, double)
Mtbl_cf	elliptic_solver_sin_zec (const Param_elliptic &, const Mtbl_cf &, double , double )
Mtbl_cf	elliptic_solver_fixe_der_zero (double, const Param_elliptic &, const Mtbl_cf &)
void	Map_af::sol_elliptic (Param_elliptic &, const Scalar &, Scalar &) const
void	Map_af::sol_elliptic_boundary (Param_elliptic &, const Scalar &, Scalar &, const Mtbl_cf &, double, double) const
void	Map_af::sol_elliptic_boundary (Param_elliptic &, const Scalar &, Scalar &, const Scalar &, double, double) const
void	Map_af::sol_elliptic_no_zec (Param_elliptic &, const Scalar &, Scalar &, double) const
void	Map_af::sol_elliptic_only_zec (Param_elliptic &, const Scalar &, Scalar &, double) const
void	Map_af::sol_elliptic_sin_zec (Param_elliptic &, const Scalar &, Scalar &, double , double ) const
void	Map_af::sol_elliptic_fixe_der_zero (double, Param_elliptic &, const Scalar &, Scalar &) const
void	Map_af::sol_elliptic_2d (Param_elliptic &, const Scalar &, Scalar &) const
void	Map_af::sol_elliptic_pseudo_1d (Param_elliptic &, const Scalar &, Scalar &) const
void	Map_log::sol_elliptic (Param_elliptic &, const Scalar &, Scalar &) const
void	Map_log::sol_elliptic_boundary (Param_elliptic &, const Scalar &, Scalar &, const Mtbl_cf &, double, double) const
void	Map_log::sol_elliptic_boundary (Param_elliptic &, const Scalar &, Scalar &, const Scalar &, double, double) const
void	Map_log::sol_elliptic_no_zec (Param_elliptic &, const Scalar &, Scalar &, double) const

Detailed Description

This class contains the parameters needed to call the general elliptic solver.

For every domain and every spherical harmonics, it contains the appropriate operator of type Ope_elementary and the appropriate variable given by a Change_var . ()

This class is only defined on an affine mapping Map_af or a logarithmic one Map_log

Definition at line 131 of file param_elliptic.h.

Constructor & Destructor Documentation

Param_elliptic::Param_elliptic ( const Scalar & so )

Standard constructor from a Scalar.

Parameters:

so [parameter] type of the source of the elliptic equation. The actual values are not used but *this will be constructed using the same number of points, domains and symetry than so .

This constructor initializes everything to solve a Poisson equation with non variable changes from domains to another.

Definition at line 179 of file param_elliptic.C.

References Itbl::annule_hard(), Scalar::annule_hard(), Valeur::base, Valeur::coef(), done_F, done_G, Scalar::dz_nonzero(), Valeur::get_base(), Scalar::get_dzpuis(), Map::get_mg(), get_mp(), Tensor::get_mp(), Mg3d::get_np(), Mg3d::get_nr(), Mg3d::get_nt(), Mg3d::get_nzone(), Scalar::get_spectral_va(), Map_log::get_type(), Base_val::give_quant_numbers(), mp_af, mp_log, operateurs, Valeur::set_base(), Tbl::set_etat_qcq(), Scalar::set_etat_qcq(), Scalar::set_spectral_va(), Scalar::std_spectral_base(), type_map, val_dF_minus, val_dF_plus, val_dG_minus, val_dG_plus, val_F_minus, val_F_plus, val_G_minus, val_G_plus, var_F, var_G, and Valeur::ylm().

Param_elliptic::~Param_elliptic ( )

Destructor.

Definition at line 306 of file param_elliptic.C.

References Map::get_mg(), get_mp(), Mg3d::get_np(), Mg3d::get_nt(), Mg3d::get_nzone(), and operateurs.

Member Function Documentation

void Param_elliptic::compute_val_F	(	int	zone,
		int	k,
		int	j
	)			const `[private]`

Computes the various values of F.

Definition at line 117 of file param_elliptic_val_lim.C.

References Tbl::annule_hard(), Valeur::base, Valeur::c_cf, done_F, Mtbl_cf::get_etat(), Map::get_mg(), get_mp(), Mg3d::get_nr(), Scalar::get_spectral_va(), Base_val::give_quant_numbers(), Itbl::set(), Tbl::set(), Tbl::set_etat_qcq(), val_dF_minus, val_dF_plus, val_F_minus, val_F_plus, and var_F.

void Param_elliptic::compute_val_G ( int zone ) const [private]

Computes the various values of G.

Definition at line 156 of file param_elliptic_val_lim.C.

References Tbl::annule_hard(), Valeur::base, Valeur::c_cf, done_G, Mtbl_cf::get_etat(), Map::get_mg(), get_mp(), Mg3d::get_nr(), Scalar::get_spectral_va(), Base_val::give_quant_numbers(), Itbl::set(), Tbl::set(), Tbl::set_etat_qcq(), val_dG_minus, val_dG_plus, val_G_minus, val_G_plus, and var_G.

double Param_elliptic::dF_minus	(	int	zone,
		int	k,
		int	j
	)			const

Returns the value of the radial derivative of F, for a given angular point, at the inner boundary of the domain zone ;.

Definition at line 75 of file param_elliptic_val_lim.C.

References compute_val_F(), done_F, and val_dF_minus.

double Param_elliptic::dF_plus	(	int	zone,
		int	k,
		int	j
	)			const

Returns the value of the radial derivative of F, for a given angular point, at the outer boundary of the domain zone ;.

Definition at line 67 of file param_elliptic_val_lim.C.

References compute_val_F(), done_F, and val_dF_plus.

double Param_elliptic::dG_minus ( int zone ) const

Returns the value of the radial derivative of G, for a given angular point, at the inner boundary of the domain zone ;.

Definition at line 108 of file param_elliptic_val_lim.C.

References compute_val_G(), done_G, and val_dG_minus.

double Param_elliptic::dG_plus ( int zone ) const

Returns the value of the radial derivative of G, for a given angular point, at the outer boundary of the domain zone ;.

Definition at line 100 of file param_elliptic_val_lim.C.

References compute_val_G(), done_G, and val_dG_plus.

double Param_elliptic::F_minus	(	int	zone,
		int	k,
		int	j
	)			const

Returns the value of F, for a given angular point, at the inner boundary of the domain zone ;.

Definition at line 59 of file param_elliptic_val_lim.C.

References compute_val_F(), done_F, and val_F_minus.

double Param_elliptic::F_plus	(	int	zone,
		int	k,
		int	j
	)			const

Returns the value of F, for a given angular point, at the outer boundary of the domain zone ;.

Definition at line 51 of file param_elliptic_val_lim.C.

References compute_val_F(), done_F, and val_F_plus.

double Param_elliptic::G_minus ( int zone ) const

Returns the value of G, for a given angular point, at the inner boundary of the domain zone ;.

Definition at line 92 of file param_elliptic_val_lim.C.

References compute_val_G(), done_G, and val_G_minus.

double Param_elliptic::G_plus ( int zone ) const

Returns the value of G, for a given angular point, at the outer boundary of the domain zone ;.

Definition at line 84 of file param_elliptic_val_lim.C.

References compute_val_G(), done_G, and val_G_plus.

const Map_radial & Param_elliptic::get_mp ( ) const

Returns the mapping.

Definition at line 114 of file param_elliptic.C.

References mp_af, mp_log, and type_map.

void Param_elliptic::inc_l_quant ( int zone )

Increases the quantum number l in the domain zone .

Definition at line 319 of file param_elliptic.C.

References Map::get_mg(), get_mp(), Mg3d::get_np(), Mg3d::get_nt(), Mg3d::get_nzone(), Ope_elementary::inc_l_quant(), and operateurs.

void Param_elliptic::set_helmholtz_minus	(	int	zone,
		double	mas,
		Scalar &	so
	)

Set the operator to $\left(\Delta - m^2\right)$ in one domain (not in the nucleus).

Parameters:

	zone	[input] : the domain.
	mas	[input] : the masse .
	so	[input] : the source (used only to get the right basis).

Definition at line 339 of file param_elliptic.C.

References Valeur::base, Valeur::coef(), Map::get_mg(), get_mp(), Mg3d::get_np(), Mg3d::get_nr(), Mg3d::get_nt(), Mg3d::get_nzone(), Scalar::get_spectral_va(), Map_log::get_type(), Base_val::give_quant_numbers(), mp_log, operateurs, Scalar::set_spectral_va(), type_map, and Valeur::ylm().

void Param_elliptic::set_helmholtz_minus_2d	(	int	zone,
		double	mas,
		const Scalar &	source
	)

Set the 2D Helmholtz operator (with minus sign).

Parameters:

	zone	[input] : the domain.
	mas	[input] : the masse parameter.

Definition at line 96 of file param_elliptic_2d.C.

References Valeur::base, Scalar::check_dzpuis(), Scalar::get_dzpuis(), Map::get_mg(), get_mp(), Mg3d::get_np(), Mg3d::get_nr(), Mg3d::get_nt(), Mg3d::get_nzone(), Scalar::get_spectral_va(), Base_val::give_quant_numbers(), operateurs, and type_map.

void Param_elliptic::set_helmholtz_minus_pseudo_1d	(	int	zone,
		double	mas,
		Scalar &	so
	)

Set the operator to $\left(\Delta - 2 \partial_r / r - m^2\right)$ in one domain.

Parameters:

	zone	[input] : the domain.
	mas	[input] : the masse .
	so	[input] : the source (used only to get the right basis).

Definition at line 91 of file param_elliptic_pseudo_1d.C.

References Valeur::base, Scalar::check_dzpuis(), Scalar::get_dzpuis(), Map::get_mg(), get_mp(), Mg3d::get_np(), Mg3d::get_nr(), Mg3d::get_nt(), Mg3d::get_nzone(), Scalar::get_spectral_va(), Base_val::give_quant_numbers(), operateurs, and type_map.

void Param_elliptic::set_helmholtz_plus	(	int	zone,
		double	mas,
		Scalar &	so
	)

Set the operator to $\left(\Delta + m^2\right)$ in one domain (only in the shells).

Parameters:

	zone	[input] : the domain.
	mas	[input] : the masse .
	so	[input] : the source (used only to get the right basis).

Definition at line 407 of file param_elliptic.C.

References Valeur::base, Valeur::coef(), Ope_elementary::get_base_r(), Map::get_mg(), get_mp(), Mg3d::get_np(), Mg3d::get_nr(), Mg3d::get_nt(), Mg3d::get_nzone(), Scalar::get_spectral_va(), Map_log::get_type(), Base_val::give_quant_numbers(), mp_log, operateurs, R_CHEBI, Scalar::set_spectral_va(), type_map, and Valeur::ylm().

void Param_elliptic::set_ope_vorton	(	int	zone,
		Scalar &	so
	)

Set the operator to $\Delta - 2\partial /\partial r$ in one domain (not implemented in the nucleus).

Parameters:

	zone	[input] : the domain.
	so	[input] : the source (used only to get the right basis).

Definition at line 701 of file param_elliptic.C.

References Valeur::base, Valeur::coef(), Scalar::get_dzpuis(), Map::get_mg(), get_mp(), Mg3d::get_np(), Mg3d::get_nr(), Mg3d::get_nt(), Mg3d::get_nzone(), Scalar::get_spectral_va(), Map_log::get_type(), Base_val::give_quant_numbers(), mp_log, operateurs, Scalar::set_spectral_va(), type_map, and Valeur::ylm().

void Param_elliptic::set_poisson_2d	(	const Scalar &	source,
		bool	indic = `false`
	)

Set everything to do a 2d-Poisson, with or without l=0 (not put by default.

..)

Definition at line 55 of file param_elliptic_2d.C.

References Valeur::base, Scalar::get_dzpuis(), Map::get_mg(), get_mp(), Mg3d::get_np(), Mg3d::get_nr(), Mg3d::get_nt(), Mg3d::get_nzone(), Scalar::get_spectral_va(), Base_val::give_quant_numbers(), operateurs, and type_map.

void Param_elliptic::set_poisson_pseudo_1d ( Scalar & so )

Set the operator to $\left(\Delta - 2 \partial_r / r\right)$ everywhere but in the compactified domain.

Parameters:

so

[input] : the source (used only to get the right basis).

Definition at line 55 of file param_elliptic_pseudo_1d.C.

References Valeur::base, Map::get_mg(), get_mp(), Mg3d::get_np(), Mg3d::get_nr(), Mg3d::get_nt(), Mg3d::get_nzone(), Scalar::get_spectral_va(), Base_val::give_quant_numbers(), operateurs, and type_map.

void Param_elliptic::set_poisson_tens_rr ( int zone )

Sets the operator to $\Delta + \frac{4}{r} \frac{\partial}{\partial r} + \frac{6 - l(l+1)}{r^2}$ in all domains, for $l \geq 2$ only.

Parameters:

zone

[input] : the domain.

Definition at line 556 of file param_elliptic.C.

References Ope_elementary::get_base_r(), Ope_poisson::get_dzpuis(), Ope_poisson::get_lquant(), Map::get_mg(), get_mp(), Mg3d::get_np(), Mg3d::get_nr(), Mg3d::get_nt(), Mg3d::get_nzone(), operateurs, and type_map.

void Param_elliptic::set_poisson_vect_eta ( int zone )

Sets the operator to be a regular elliptic operator to solve for the $\eta$ component of the vector Poisson equation.

The operator is $\frac{\partial^2}{\partial r^2} + \frac{2}{r} \frac{\partial}{\partial r} - \frac{l(l-1)}{r^2}$ (Poisson with the decrease of l by one unit) in all domains but the CED, for $l \not= 0$ ; it is not defined for l = 0. In the CED, the operator is also the Laplace one, but with l increased by one unit: $\frac{\partial^2}{\partial r^2} + \frac{2}{r} \frac{\partial}{\partial r} - \frac{(l+1)(l+2)}{r^2}$ . This is intended to solve the equation for $\eta$ arising in the decomposition of the vector Poisson equation.

Parameters:

zone

[input] : the domain.

Definition at line 529 of file param_elliptic.C.

References Ope_poisson::dec_l_quant(), Ope_poisson::get_lquant(), Map::get_mg(), get_mp(), Mg3d::get_np(), Mg3d::get_nt(), Mg3d::get_nzone(), Mg3d::get_type_r(), Ope_poisson::inc_l_quant(), and operateurs.

void Param_elliptic::set_poisson_vect_r	(	int	zone,
		bool	only_l_zero = `false`
	)

Sets the operator to $\Delta + \frac{2}{r} \frac{\partial}{\partial r} + \frac{2 - l(l+1)}{r^2}$ in all domains, for $l \not= 0$ ; and to $\frac{\partial^2}{\partial r^2} + \frac{2}{r} \frac{\partial}{\partial r} - \frac{2}{r^2}$ in all domains otherwise.

Parameters: