20 #include "headcpp.hpp"
21 #include "spheric_symphi.hpp"
22 #include "array_math.hpp"
24 #include "tensor_impl.hpp"
30 const Param& par,
int type_exception,
const Val_domain& exception)
const {
38 assert ((type_exception>=1) && (type_exception <=4)) ;
72 cerr <<
"Unknow phi basis in Domain_shell_symphi::export_tau_val_domain_exception" << endl ;
77 for (
int k=kmin ; k<=kmax ; k++) {
101 cerr <<
"Unknow phi basis in Domain_shell_symphi::export_tau_val_domain_exception" << endl ;
109 if ((k==ktarget) && (j==jtarget)) {
110 switch (type_exception) {
115 sec.
set(pos_sec) = 0 ;
121 sec.
set(pos_sec) = 0 ;
124 cerr <<
"bad value for type_exception" << endl ;
134 pos_galerkin = pos_cf ;
135 pos_galerkin.
set(1) = 0 ;
136 if ((k==ktarget) && (j==jtarget)) {
137 switch (type_exception) {
142 sec.
set(pos_sec) = 0 ;
148 sec.
set(pos_sec) = 0 ;
151 cerr <<
"bad value for type_exception" << endl ;
165 if ((k==ktarget) && (j==jtarget)) {
166 switch (type_exception) {
171 sec.
set(pos_sec) = 0 ;
177 sec.
set(pos_sec) = 0 ;
180 cerr <<
"bad value for type_exception" << endl ;
190 pos_galerkin = pos_cf ;
191 pos_galerkin.
set(1) = 0 ;
192 if ((k==ktarget) && (j==jtarget)) {
193 switch (type_exception) {
198 sec.
set(pos_sec) = 0 ;
204 sec.
set(pos_sec) = 0 ;
207 cerr <<
"bad value for type_exception" << endl ;
220 if ((k==ktarget) && (j==jtarget)) {
221 switch (type_exception) {
226 sec.
set(pos_sec) = 0 ;
232 sec.
set(pos_sec) = 0 ;
235 cerr <<
"bad value for type_exception" << endl ;
245 pos_galerkin = pos_cf ;
246 pos_galerkin.
set(1) = 1 ;
247 if ((k==ktarget) && (j==jtarget)) {
248 switch (type_exception) {
253 sec.
set(pos_sec) = 0 ;
259 sec.
set(pos_sec) = 0 ;
262 cerr <<
"bad value for type_exception" << endl ;
276 if ((k==ktarget) && (j==jtarget)) {
277 switch (type_exception) {
282 sec.
set(pos_sec) = 0 ;
288 sec.
set(pos_sec) = 0 ;
291 cerr <<
"bad value for type_exception" << endl ;
301 pos_galerkin = pos_cf ;
302 pos_galerkin.
set(1) = 0 ;
303 if ((k==ktarget) && (j==jtarget)) {
304 switch (type_exception) {
309 sec.
set(pos_sec) = 0 ;
315 sec.
set(pos_sec) = 0 ;
318 cerr <<
"bad value for type_exception" << endl ;
331 cerr <<
"Unknow theta basis in Domain_shell_symphi::export_tau_val_domain_boundary_exception" << endl ;
340 const Param& par,
int type_exception,
const Tensor& exception,
343 if ((bound!=INNER_BC) && (bound!=OUTER_BC)) {
344 cerr <<
"Unknown boundary_exception in Domain_shell_symphi::export_tau_boundary_exception" << endl ;
354 cerr <<
"Valence " << val <<
" not implemented in Domain_shell_symphi::export_tau_boundary_exception" << endl ;
reference set(const Index &pos)
Read/write of an element.
Bases_container bases_1d
Arrays containing the various basis of decomposition.
void export_tau_val_domain_boundary_exception(const Val_domain &eq, int mlim, int bound, Array< double > &res, int &pos_res, int ncond, const Param ¶m, int type_exception, const Val_domain &exception) const
Exports all the residual equations corresponding to one tensorial one on a given boundary,...
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_exception(const Tensor &, int, int, Array< double > &, int &, const Array< int > &, const Param &, int, const Tensor &, int n_cmp=-1, Array< int > **p_cmp=0x0) const
Exports all the residual equations corresponding to one tensorial one on a given boundary,...
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.
const double & get_double(int position=0) const
Returns the reference of a double stored in the list.
const int & get_int(int position=0) const
Returns the reference of a int stored in the list.
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.