20 #include "headcpp.hpp"
21 #include "bispheric.hpp"
22 #include "array_math.hpp"
24 #include "val_domain.hpp"
34 other_bound = CHI_ONE_BC ;
39 other_bound = ETA_MINUS_BC ;
43 other_bound = ETA_PLUS_BC ;
47 cerr <<
"Unknown boundary case in Domain_bispheric_rect::find_other_dom" << endl ;
60 Index copie_pos (pos_cf) ;
64 copie_pos.
set(0) = i ;
73 copie_pos.
set(0) = i ;
79 copie_pos.
set(1) = j ;
85 copie_pos.
set(1) = j ;
87 copie_pos.
set(0) = i ;
93 cerr <<
"Unknown boundary type in Domain_bispheric_rect::val_boundary" << endl ;
114 cerr <<
"Domain_bispheric_rect::nbr_points_boundary not yet implemented for boundary " << bound << endl ;
122 int pos_which = start ;
126 int mink = (basep==COS) ? 0 : 1 ;
135 which_coef[pos_which] =
new Index(pos) ;
138 for (k=mink ; k<maxk ; k++) {
142 which_coef[pos_which] =
new Index(pos) ;
150 for (k=mink ; k<maxk ; k++) {
152 which_coef[pos_which] =
new Index(pos) ;
157 cerr <<
"Unknown boundary in Domain_bispheric_rect::do_which_points_inside" << endl ;
Class for storing the basis of decompositions of a field.
Bases_container bases_1d
Arrays containing the various basis of decomposition.
double eta_minus
associated with .
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 int nbr_points_boundary(int, const Base_spectral &) const
Computes the number of relevant collocation points on a boundary.
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.
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