20 #include "headcpp.hpp"
21 #include "polar_periodic.hpp"
23 #include "array_math.hpp"
25 #include "tensor_impl.hpp"
49 cerr <<
"Unknown time basis in Domain_polar_periodic_shell_export_tau_val_domain_boundary" << endl ;
54 for (
int k=mink ; k<maxk ; k++) {
79 cerr <<
"Unknown theta basis in Domain_polar_periodic_shell_export_tau_val_domain" << endl ;
84 for (
int j=minj ; j<maxj ; j++) {
96 #ifndef REMOVE_ALL_CHECKS
98 if ((bound!=OUTER_BC) && (bound!=INNER_BC)) {
99 cerr <<
"Unknown boundary in Domain_polar_periodic_shell::export_tau_boundary" << endl ;
115 else for (
int i=0 ; i<n_cmp ; i++) {
116 if ((*p_cmp[i])(0)==1)
118 if ((*p_cmp[i])(0)==2)
120 if ((*p_cmp[i])(0)==3)
134 else for (
int i=0 ; i<n_cmp ; i++) {
135 if (((*p_cmp[i])(0)==1) && ((*p_cmp[i])(1)==1))
137 if (((*p_cmp[i])(0)==1) && ((*p_cmp[i])(1)==2))
139 if (((*p_cmp[i])(0)==1) && ((*p_cmp[i])(1)==3))
141 if (((*p_cmp[i])(0)==2) && ((*p_cmp[i])(1)==2))
143 if (((*p_cmp[i])(0)==2) && ((*p_cmp[i])(1)==3))
145 if (((*p_cmp[i])(0)==3) && ((*p_cmp[i])(1)==3))
163 else for (
int i=0 ; i<n_cmp ; i++) {
164 if (((*p_cmp[i])(0)==1) && ((*p_cmp[i])(1)==1))
166 if (((*p_cmp[i])(0)==1) && ((*p_cmp[i])(1)==2))
168 if (((*p_cmp[i])(0)==1) && ((*p_cmp[i])(1)==3))
170 if (((*p_cmp[i])(0)==2) && ((*p_cmp[i])(1)==1))
172 if (((*p_cmp[i])(0)==2) && ((*p_cmp[i])(1)==2))
174 if (((*p_cmp[i])(0)==2) && ((*p_cmp[i])(1)==3))
176 if (((*p_cmp[i])(0)==3) && ((*p_cmp[i])(1)==1))
178 if (((*p_cmp[i])(0)==3) && ((*p_cmp[i])(1)==2))
180 if (((*p_cmp[i])(0)==3) && ((*p_cmp[i])(1)==3))
186 cerr <<
"Valence " << val <<
" not implemented in Domain_polar_periodic_shell::export_tau_boundary" << endl ;
reference set(const Index &pos)
Read/write of an element.
Bases_container bases_1d
Arrays containing the various basis of decomposition.
virtual double val_boundary(int, const Val_domain &, const Index &) const
Computes the value of a field at a boundary.
virtual void export_tau_boundary(const Tensor &, int, int, Array< double > &, int &, const Array< int > &, int n_cmp=-1, Array< int > **p_cmp=0x0) const
Exports all the residual equations corresponding to a tensorial one on a given boundary It makes use ...
void export_tau_val_domain_boundary(const Val_domain &eq, int bound, Array< double > &res, int &pos_res, int ncond) const
Exports all the residual equations corresponding to a tensorial one on a given boundary It makes use ...
Dim_array nbr_coefs
Number of coefficients.
Class that gives the position inside a multi-dimensional Array.
int & set(int i)
Read/write of the position in a given dimension.
int get_n_comp() const
Returns the number of stored components.
int get_valence() const
Returns the valence.
Class for storing the basis of decompositions of a field and its values on both the configuration and...
bool check_if_zero() const
Check whether the logical state is zero or not.
void coef() const
Computes the coefficients.
const Base_spectral & get_base() const
Returns the basis of decomposition.