32 template <
unsigned int dim>
39 std::shared_ptr<BulkIntegral> bulk_integral = std::make_shared<BulkIntegral>(shared_from_this(), dim, i_subset);
45 std::shared_ptr<BulkIntegral> EvalPoints::add_bulk<0>(
const Quadrature &quad)
49 std::shared_ptr<BulkIntegral> bulk_integral = std::make_shared<BulkIntegral>(shared_from_this(), 0, i_subset);
54 template <
unsigned int dim>
59 for (
unsigned int i=0; i<dim+1; ++i) {
64 std::shared_ptr<EdgeIntegral> edge_integral = std::make_shared<EdgeIntegral>(shared_from_this(), dim, i_subset);
70 template <
unsigned int dim>
74 std::shared_ptr<BulkIntegral> bulk_integral = this->
add_bulk<dim-1>(quad);
75 DebugOut() <<
"coupling bulk subset" << bulk_integral->get_subset_idx();
76 std::shared_ptr<EdgeIntegral> edge_integral = this->add_edge<dim>(quad);
77 DebugOut() <<
"coupling edge subset" << edge_integral->get_subset_idx();
78 return std::make_shared<CouplingIntegral>(edge_integral, bulk_integral);
81 template <
unsigned int dim>
85 std::shared_ptr<BulkIntegral> bulk_integral = this->
add_bulk<dim-1>(quad);
86 DebugOut() <<
"boundary bulk subset: " << bulk_integral->get_subset_idx()
87 <<
"begin: " <<
subset_begin(dim-1, bulk_integral->get_subset_idx());
88 std::shared_ptr<EdgeIntegral> edge_integral = this->add_edge<dim>(quad);
89 DebugOut() <<
"boundary edge subset" << edge_integral->get_subset_idx();
90 return std::make_shared<BoundaryIntegral>(edge_integral, bulk_integral);
94 : local_points_(dim), n_subsets_(0), dim_(dim)
100 template <
unsigned int dim>
102 ASSERT_GT(dim, 0).error(
"Dimension 0 not supported!\n");
103 unsigned int local_points_old_size = local_points_.size();
104 local_points_.resize(quad_points.
size() + local_points_old_size);
105 for (
unsigned int i=0; i<quad_points.
size(); ++i) {
108 local_points_.set(i+local_points_old_size) = quad_points.
vec<dim>(i);
117 subset_starts_[n_subsets_] = this->
size();
118 return n_subsets_ - 1;
122 template std::shared_ptr<BulkIntegral> EvalPoints::add_bulk<0>(
const Quadrature &);
123 template std::shared_ptr<BulkIntegral> EvalPoints::add_bulk<1>(
const Quadrature &);
124 template std::shared_ptr<BulkIntegral> EvalPoints::add_bulk<2>(
const Quadrature &);
125 template std::shared_ptr<BulkIntegral> EvalPoints::add_bulk<3>(
const Quadrature &);
126 template std::shared_ptr<EdgeIntegral> EvalPoints::add_edge<1>(
const Quadrature &);
127 template std::shared_ptr<EdgeIntegral> EvalPoints::add_edge<2>(
const Quadrature &);
128 template std::shared_ptr<EdgeIntegral> EvalPoints::add_edge<3>(
const Quadrature &);
129 template std::shared_ptr<CouplingIntegral> EvalPoints::add_coupling<2>(
const Quadrature &);
130 template std::shared_ptr<CouplingIntegral> EvalPoints::add_coupling<3>(
const Quadrature &);
131 template std::shared_ptr<BoundaryIntegral> EvalPoints::add_boundary<1>(
const Quadrature &);
132 template std::shared_ptr<BoundaryIntegral> EvalPoints::add_boundary<2>(
const Quadrature &);
133 template std::shared_ptr<BoundaryIntegral> EvalPoints::add_boundary<3>(
const Quadrature &);
#define ASSERT_LT(a, b)
Definition of comparative assert macro (Less Than) only for debug mode.
#define ASSERT_EQ(a, b)
Definition of comparative assert macro (EQual) 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)
std::array< int, EvalPoints::max_subsets+1 > subset_starts_
Indices of subsets data in local_points_ vector, used size is n_subsets_ + 1.
DimEvalPoints(unsigned int dim)
Constructor.
uint add_subset()
Adds new subset and its end size to subset_starts_ array.
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_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 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.
static const unsigned int max_subset_points
Maximal average number of points hold in 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.
Base class for quadrature rules on simplices in arbitrary dimensions.
Quadrature make_from_side(unsigned int sid) const
const Armor::Array< double > & get_points() const
Return a reference to the whole array of quadrature points.
#define DebugOut()
Macro defining 'debug' record of log.
Basic definitions of numerical quadrature rules.
Class RefElement defines numbering of vertices, sides, calculation of normal vectors etc.