20 #include "headcpp.hpp"
22 #include "spheric.hpp"
24 #include "val_domain.hpp"
31 other_bound = OUTER_BC ;
35 other_bound = INNER_BC ;
38 cerr <<
"Unknown boundary case in Domain_shell::find_other_dom" << endl ;
51 Index copie_pos (pos_cf) ;
55 copie_pos.
set(0) = i ;
64 copie_pos.
set(0) = i ;
69 cerr <<
"Unknown boundary type in Domain_shell::val_boundary" << endl ;
78 if ((bound!=INNER_BC) && (bound!=OUTER_BC)) {
79 cerr <<
"Unknown boundary in Domain_shell::nbr_points_boundary" << endl ;
84 int first_base_theta = (*bb.
bases_1d[1]) (0) ;
85 int nbrj = ((first_base_theta==COS_EVEN) || (first_base_theta==SIN_ODD)) ?
nbr_points(1)-1 :
nbr_points(1)-2 ;
92 int pos_which = start ;
103 cerr <<
"Unknown boundary in Domain_shell::do_which_points_inside" << endl ;
109 int first_base_theta = (*bb.
bases_1d[1]) (0) ;
110 int maxj = ((first_base_theta==COS_EVEN) || (first_base_theta==SIN_ODD)) ?
nbr_points(1) :
nbr_points(1)-1 ;
114 for (
int j=0 ; j<maxj ; j++) {
116 if ((k==0) || (j!=0)) {
117 which_coef[pos_which] =
new Index(pos) ;
Class for storing the basis of decompositions of a field.
Bases_container bases_1d
Arrays containing the various basis of decomposition.
virtual void do_which_points_boundary(int, const Base_spectral &, Index **, int) const
Lists all the indices corresponding to true collocation points on a boundary.
virtual void find_other_dom(int, int, int &, int &) const
Gives the informations corresponding the a touching neighboring domain.
virtual double val_boundary(int, const Val_domain &, const Index &) const
Computes the value of a field at a boundary.
virtual int nbr_points_boundary(int, const Base_spectral &) const
Computes the number of relevant collocation points on a boundary.
Dim_array nbr_coefs
Number of coefficients.
Dim_array nbr_points
Number of colocation points.
Class that gives the position inside a multi-dimensional Array.
int & set(int i)
Read/write of the position in a given dimension.
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.
Array< double > get_coef() const