61 #include "param_elliptic.h" 68 assert(dynamic_cast<const Base_vect_spher*>(
triad) != 0x0) ;
70 assert(dynamic_cast<const Map_af*>(
mp) != 0x0) ;
71 assert( (dzfin == 0) || (dzfin == 2) ) ;
81 if (source_rr.
get_etat() != ETATZERO) {
88 for (
int lz=0; lz<nz; lz++)
97 resu.
set(1,1) = h_rr ;
115 Scalar source_eta = -resu(1,1).dsdr() ;
117 source_eta -= 3.*resu(1,1) ;
148 Scalar tautst = resu(1,2).dsdr() ;
158 tautst += 3 * resu(1,2) - resu(1,1).dsdt() ;
172 Scalar taupst = - resu(1,3).dsdr() ;
183 taupst -= 3 * resu(1,3) ;
226 *rcmp[
position(idx)] = - resu(1,1) - resu(3,3) ;
Sym_tensor_tt poisson(int dzfin=2) const
Computes the solution of a tensorial TT Poisson equation with *this as a source: ...
const Scalar & mu(Param *par=0x0) const
Gives the field (see member p_mu ).
int & set(int i)
Read/write of a particular element (index i ) (1D case)
void ylm_i()
Inverse of ylm()
virtual void set_etat_zero()
Sets the logical state to ETATZERO (zero).
const Metric *const met_div
Metric with respect to which the divergence and the trace are defined.
const Scalar & dsdt() const
Returns 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).
Scalar poisson() const
Solves the scalar Poisson equation with *this as a source.
void mult_sint()
Multiplication by .
Basic integer array class.
const Base_vect * triad
Vectorial basis (triad) with respect to which the tensor components are defined.
int get_etat() const
Returns the logical state ETATNONDEF (undefined), ETATZERO (null) or ETATQCQ (ordinary).
Scalar poisson_angu(double lambda=0) const
Solves the (generalized) angular Poisson equation with *this as source.
virtual void inc_dzpuis(int inc=1)
Increases by inc units the value of dzpuis and changes accordingly the values of the Scalar in the co...
Scalar ** cmp
Array of size n_comp of pointers onto the components.
Scalar sol_elliptic(Param_elliptic ¶ms) const
Resolution of a general elliptic equation, putting zero at infinity.
int get_nzone() const
Returns the number of domains.
void set_poisson_tens_rr(int zone)
Sets the operator to in all domains, for only.
const Scalar & stdsdp() const
Returns of *this .
void mult_cost()
Multiplication by .
This class contains the parameters needed to call the general elliptic solver.
Scalar & set(const Itbl &ind)
Returns the value of a component (read/write version).
virtual int position(const Itbl &ind) const
Returns the position in the array cmp of a component given by its indices.
const Scalar & operator()(const Itbl &ind) const
Returns the value of a component (read-only version).
void div_sint()
Division by .
Valeur & set_spectral_va()
Returns va (read/write version)
const Map *const mp
Mapping on which the numerical values at the grid points are defined.
void div_r_dzpuis(int ced_mult_r)
Division by r everywhere but with the output flag dzpuis set to ced_mult_r .
void mult_r_dzpuis(int ced_mult_r)
Multiplication by r everywhere but with the output flag dzpuis set to ced_mult_r ...
Transverse and traceless symmetric tensors of rank 2.