19 #ifndef EVAL_POINTS_HH_
20 #define EVAL_POINTS_HH_
43 class EvalPoints :
public std::enable_shared_from_this<EvalPoints> {
58 inline unsigned int size(
unsigned int dim)
const {
63 template<
unsigned int dim>
64 inline arma::vec::fixed<dim>
local_point(
unsigned int local_point_idx)
const {
65 ASSERT_GT(dim, 0).error(
"Dimension 0 not supported!\n");
70 inline int subset_begin(
unsigned int dim,
unsigned int idx)
const {
75 inline int subset_end(
unsigned int dim,
unsigned int idx)
const {
80 inline int subset_size(
unsigned int dim,
unsigned int idx)
const {
85 inline unsigned int n_subsets(
unsigned int dim)
const {
93 template <
unsigned int dim>
97 template <
unsigned int dim>
101 template <
unsigned int dim>
105 template <
unsigned int dim>
114 for (
uint i=0; i<4; ++i)
126 inline unsigned int size()
const {
132 template<
unsigned int dim>
133 inline arma::vec::fixed<dim>
local_point(
unsigned int local_point_idx)
const {
162 template <
unsigned int dim>
#define ASSERT_LT(a, b)
Definition of comparative assert macro (Less Than) only for debug mode.
#define ASSERT_GT(a, b)
Definition of comparative assert macro (Greater Than) only for debug mode.
ArmaVec< Type, nr > vec(uint mat_index) const
unsigned int size() const
Subobject holds evaluation points data of one dimension (0,1,2,3)
unsigned int n_subsets_
Number of subset.
int subset_begin(unsigned int idx) const
Return begin index of appropriate subset data.
std::array< int, EvalPoints::max_subsets+1 > subset_starts_
Indices of subsets data in local_points_ vector, used size is n_subsets_ + 1.
unsigned int size() const
Return size of evaluation points object (number of points).
arma::vec::fixed< dim > local_point(unsigned int local_point_idx) const
Return local coordinates of given local point.
DimEvalPoints(unsigned int dim)
Constructor.
unsigned int n_subsets() const
Return number of subsets.
uint add_subset()
Adds new subset and its end size to subset_starts_ array.
unsigned int dim_
Dimension of local points.
int subset_end(unsigned int idx) const
Return end index of appropriate subset data.
Armor::Array< double > local_points_
Local coords of points vector.
void add_local_points(const Armor::Array< double > &quad_points)
Adds set of local point to local_points_ (bulk or side).
int subset_size(unsigned int idx) const
Return number of local points corresponding to subset.
Class holds local coordinations of evaluating points (bulk and sides) specified by element dimension.
int subset_begin(unsigned int dim, unsigned int idx) const
Return begin index of appropriate subset data.
std::shared_ptr< CouplingIntegral > add_coupling(const Quadrature &)
The same as add_bulk but for points between side points of element of dim and bulk points of element ...
static const unsigned int undefined_dim
Undefined dimension of new (empty) object.
unsigned int n_subsets(unsigned int dim) const
Return number of subsets.
unsigned int size(unsigned int dim) const
Return size of evaluation points object (number of points).
std::shared_ptr< BoundaryIntegral > add_boundary(const Quadrature &)
The same as add_bulk but for edge points on boundary sides.
unsigned int max_size() const
Return maximal size of evaluation points objects.
int subset_end(unsigned int dim, unsigned int idx) const
Return end index of appropriate subset data.
unsigned int max_size_
Maximal number of used EvalPoints.
static const unsigned int max_subset_points
Maximal average number of points hold in subset.
int subset_size(unsigned int dim, unsigned int idx) const
Return number of local points corresponding to subset.
std::shared_ptr< BulkIntegral > add_bulk(const Quadrature &)
std::shared_ptr< EdgeIntegral > add_edge(const Quadrature &)
The same as add_bulk but for edge points on sides.
std::array< DimEvalPoints, 4 > dim_eval_points_
Sub objects of dimensions 0,1,2,3.
static constexpr unsigned int max_subsets
Maximal number of hold subsets.
arma::vec::fixed< dim > local_point(unsigned int local_point_idx) const
Return local coordinates of given local point and appropriate dim.
Base class for quadrature rules on simplices in arbitrary dimensions.
Implementation of range helper class.