58 "Lowest order BDM element are degree 1, but you asked for degree 0"));
87 const unsigned int face_no = 0;
90 for (
unsigned int i = 0; i < GeometryInfo<dim>::max_children_per_face; ++i)
97 unsigned int target_row = 0;
98 for (
unsigned int d = 0; d < GeometryInfo<dim>::max_children_per_face; ++d)
99 for (
unsigned int i = 0; i < face_embeddings[d].
m(); ++i)
101 for (
unsigned int j = 0; j < face_embeddings[d].
n(); ++j)
140 std::ostringstream namebuf;
141 namebuf <<
"FE_BDM<" << dim <<
">(" << this->
degree - 1 <<
")";
143 return namebuf.str();
148std::unique_ptr<FiniteElement<dim, dim>>
151 return std::make_unique<FE_BDM<dim>>(*this);
160 std::vector<double> &nodal_values)
const
162 Assert(support_point_values.size() == this->generalized_support_points.size(),
164 this->generalized_support_points.size()));
166 Assert(nodal_values.size() == this->n_dofs_per_cell(),
173 unsigned int dbase = 0;
175 unsigned int pbase = 0;
183 nodal_values[dbase + i] =
184 support_point_values[pbase + i]
186 pbase += this->n_dofs_per_face(f);
198 nodal_values[dbase + i] = s;
206 const unsigned int face_no = 0;
233 nodal_values[dbase + i] = s;
243std::vector<unsigned int>
256 std::vector<unsigned int> dpo(dim + 1, 0u);
273 return std::vector<bool>();
311 initialize_test_values(std::vector<std::vector<double>> &test_values,
313 const unsigned int deg)
316 std::vector<Tensor<1, dim>> dummy1;
317 std::vector<Tensor<2, dim>> dummy2;
318 std::vector<Tensor<3, dim>> dummy3;
319 std::vector<Tensor<4, dim>> dummy4;
321 test_values.resize(quadrature.
size());
323 for (
unsigned int k = 0; k < quadrature.
size(); ++k)
325 test_values[k].resize(poly.
n());
332 for (
unsigned int i = 0; i < poly.
n(); ++i)
334 test_values[k][i] *= quadrature.
weight(k);
344 initialize_test_values(std::vector<std::vector<double>> &,
362 const QGauss<dim - 1> face_points(deg + 1);
367 const unsigned int face_no = 0;
371 for (
unsigned int k = 0; k < face_points.size(); ++k)
381 const unsigned int npoints =
389 for (
unsigned int f = 0; f < GeometryInfo<dim>::faces_per_cell; ++f)
396 for (
unsigned int k = 0; k < face_points.size(); ++k)
398 faces.
point(offset + k);
413 const unsigned int ibase =
415 for (
unsigned int k = 0; k < cell_points.size(); ++k)
431#include "fe_bdm.inst"
std::vector< std::vector< double > > test_values_face
void initialize_support_points(const unsigned int bdm_degree)
FE_BDM(const unsigned int p)
std::vector< std::vector< double > > test_values_cell
static std::vector< bool > get_ria_vector(const unsigned int degree)
virtual std::unique_ptr< FiniteElement< dim, dim > > clone() const override
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
void initialize_quad_dof_index_permutation_and_sign_change()
virtual std::string get_name() const override
static std::vector< unsigned int > get_dpo_vector(const unsigned int degree)
FullMatrix< double > inverse_node_matrix
std::vector< MappingKind > mapping_kind
FE_PolyTensor(const TensorPolynomialsBase< dim > &polynomials, const FiniteElementData< dim > &fe_data, const std::vector< bool > &restriction_is_additive_flags, const std::vector< ComponentMask > &nonzero_components)
const unsigned int degree
unsigned int n_dofs_per_cell() const
const unsigned int dofs_per_face
unsigned int n_dofs_per_face(unsigned int face_no=0, unsigned int child=0) const
ReferenceCell reference_cell() const
unsigned int n_unique_faces() const
const unsigned int dofs_per_cell
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())
void reinit_restriction_and_prolongation_matrices(const bool isotropic_restriction_only=false, const bool isotropic_prolongation_only=false)
std::vector< std::vector< Point< dim - 1 > > > generalized_face_support_points
FullMatrix< double > interface_constraints
std::vector< Point< dim > > generalized_support_points
std::vector< std::vector< FullMatrix< double > > > prolongation
static unsigned int n_polynomials(const unsigned int n)
void evaluate(const Point< dim > &unit_point, std::vector< double > &values, std::vector< Tensor< 1, dim > > &grads, std::vector< Tensor< 2, dim > > &grad_grads, std::vector< Tensor< 3, dim > > &third_derivatives, std::vector< Tensor< 4, dim > > &fourth_derivatives) const override
static unsigned int n_polynomials(const unsigned int degree)
static DataSetDescriptor face(const ReferenceCell &reference_cell, const unsigned int face_no, const bool face_orientation, const bool face_flip, const bool face_rotation, const unsigned int n_quadrature_points)
static Quadrature< dim > project_to_all_faces(const ReferenceCell &reference_cell, const hp::QCollection< dim - 1 > &quadrature)
const Point< dim > & point(const unsigned int i) const
double weight(const unsigned int i) const
unsigned int size() const
static constexpr unsigned char default_combined_face_orientation()
#define DEAL_II_NAMESPACE_OPEN
#define DEAL_II_NAMESPACE_CLOSE
#define Assert(cond, exc)
static ::ExceptionBase & ExcImpossibleInDim(int arg1)
#define AssertDimension(dim1, dim2)
static ::ExceptionBase & ExcInternalError()
static ::ExceptionBase & ExcDimensionMismatch(std::size_t arg1, std::size_t arg2)
static ::ExceptionBase & ExcMessage(std::string arg1)
static constexpr unsigned int max_children_per_face
static constexpr unsigned int faces_per_cell
static std_cxx20::ranges::iota_view< unsigned int, unsigned int > face_indices()
static constexpr std::array< unsigned int, faces_per_cell > unit_normal_direction