32 const unsigned int order,
53std::vector<unsigned int>
56 std::vector<unsigned int> dpo(dim + 1, 0);
68 std::ostringstream namebuf;
69 namebuf <<
"FE_RannacherTurek"
78std::unique_ptr<FiniteElement<dim, dim>>
81 return std::make_unique<FE_RannacherTurek<dim>>(this->
order,
93 this->
weights = face_quadrature.get_weights();
95 for (
unsigned int q = 0; q < face_quadrature.size(); ++q)
114 std::vector<double> &nodal_values)
const
117 this->generalized_support_points.size());
120 const unsigned int q_points_per_face = this->
weights.size();
121 std::fill(nodal_values.begin(), nodal_values.end(), 0.0);
123 std::vector<Vector<double>>::const_iterator
value =
124 support_point_values.begin();
127 for (
unsigned int q = 0; q < q_points_per_face; ++q)
129 nodal_values[face] += (*value)[0] * this->
weights[q];
138#include "fe_rannacher_turek.inst"
FE_Poly(const ScalarPolynomialsBase< dim > &poly_space, const FiniteElementData< dim > &fe_data, const std::vector< bool > &restriction_is_additive_flags, const std::vector< ComponentMask > &nonzero_components)
virtual std::unique_ptr< FiniteElement< dim, dim > > clone() const override
std::vector< double > weights
const unsigned int n_face_support_points
FE_RannacherTurek(const unsigned int order=0, const unsigned int n_face_support_points=2)
virtual std::string get_name() const override
void initialize_support_points()
virtual void convert_generalized_support_point_values_to_dof_values(const std::vector< Vector< double > > &support_point_values, std::vector< double > &nodal_values) const override
std::vector< unsigned int > get_dpo_vector()
FiniteElementData(const std::vector< unsigned int > &dofs_per_object, const unsigned int n_components, const unsigned int degree, const Conformity conformity=unknown, const BlockIndices &block_indices=BlockIndices())
std::vector< Point< dim > > generalized_support_points
#define DEAL_II_NAMESPACE_OPEN
#define DEAL_II_NAMESPACE_CLOSE
static ::ExceptionBase & ExcNotImplemented()
#define Assert(cond, exc)
#define AssertDimension(dim1, dim2)
const bool IsBlockVector< VectorType >::value
static std_cxx20::ranges::iota_view< unsigned int, unsigned int > face_indices()