20 #include "headcpp.hpp"
22 #include "spheric.hpp"
24 #include "val_domain.hpp"
31 other_bound = OUTER_BC ;
34 cerr <<
"Unknown boundary case in Domain_compact::find_other_dom" << endl ;
46 Index copie_pos (pos_cf) ;
50 copie_pos.
set(0) = i ;
59 copie_pos.
set(0) = i ;
64 cerr <<
"Unknown boundary type in Domain_compact::val_boundary" << endl ;
73 if ((bound!=INNER_BC) && (bound!=OUTER_BC)) {
74 cerr <<
"Unknown boundary in Domain_compac::nbr_points_boundary" << endl ;
78 int first_base_theta = (*bb.
bases_1d[1]) (0) ;
79 int nbrj = ((first_base_theta==COS_EVEN) || (first_base_theta==SIN_ODD)) ?
nbr_points(1)-1 :
nbr_points(1)-2 ;
86 int pos_which = start ;
96 cerr <<
"Unknown boundary in Domain_compact::do_which_points_inside" << endl ;
100 int first_base_theta = (*bb.
bases_1d[1]) (0) ;
101 int maxj = ((first_base_theta==COS_EVEN) || (first_base_theta==SIN_ODD)) ?
nbr_points(1) :
nbr_points(1)-1 ;
105 for (
int j=0 ; j<maxj ; j++) {
107 if ((k==0) || (j!=0)) {
108 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 double val_boundary(int, const Val_domain &, const Index &) const
Computes the value of a field at a boundary.
virtual void find_other_dom(int, int, int &, int &) const
Gives the informations corresponding the a touching neighboring domain.
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 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