17 #include <deal.II/base/polynomials_p.h> 18 #include <deal.II/base/qprojector.h> 19 #include <deal.II/base/quadrature_lib.h> 20 #include <deal.II/base/std_cxx14/memory.h> 21 #include <deal.II/base/table.h> 23 #include <deal.II/dofs/dof_accessor.h> 25 #include <deal.II/fe/fe.h> 26 #include <deal.II/fe/fe_bdm.h> 27 #include <deal.II/fe/fe_tools.h> 28 #include <deal.II/fe/fe_values.h> 29 #include <deal.II/fe/mapping.h> 31 #include <deal.II/grid/tria.h> 32 #include <deal.II/grid/tria_iterator.h> 38 DEAL_II_NAMESPACE_OPEN
50 std::vector<bool>(dim, true)))
56 "Lowest order BDM element are degree 1, but you asked for degree 0"));
85 for (
unsigned int i = 0; i < GeometryInfo<dim>::max_children_per_face; ++i)
90 unsigned int target_row = 0;
91 for (
unsigned int d = 0; d < GeometryInfo<dim>::max_children_per_face; ++d)
92 for (
unsigned int i = 0; i < face_embeddings[d].
m(); ++i)
94 for (
unsigned int j = 0; j < face_embeddings[d].
n(); ++j)
117 std::ostringstream namebuf;
118 namebuf <<
"FE_BDM<" << dim <<
">(" << this->
degree - 1 <<
")";
120 return namebuf.str();
125 std::unique_ptr<FiniteElement<dim, dim>>
128 return std_cxx14::make_unique<FE_BDM<dim>>(*this);
137 std::vector<double> & nodal_values)
const 150 unsigned int dbase = 0;
152 unsigned int pbase = 0;
153 for (
unsigned int f = 0; f < GeometryInfo<dim>::faces_per_cell; ++f)
160 nodal_values[dbase + i] =
161 support_point_values[pbase + i]
175 nodal_values[dbase + i] = s;
205 nodal_values[dbase + i] = s;
215 std::vector<unsigned int>
228 std::vector<unsigned int> dpo(dim + 1, 0u);
245 return std::vector<bool>();
260 std::vector<bool> ret_val(dofs_per_cell,
false);
283 initialize_test_values(std::vector<std::vector<double>> &test_values,
285 const unsigned int deg)
288 std::vector<Tensor<1, dim>> dummy1;
289 std::vector<Tensor<2, dim>> dummy2;
290 std::vector<Tensor<3, dim>> dummy3;
291 std::vector<Tensor<4, dim>> dummy4;
293 test_values.resize(quadrature.
size());
295 for (
unsigned int k = 0; k < quadrature.
size(); ++k)
297 test_values[k].resize(poly.
n());
304 for (
unsigned int i = 0; i < poly.
n(); ++i)
306 test_values[k][i] *= quadrature.
weight(k);
316 initialize_test_values(std::vector<std::vector<double>> &,
334 QGauss<dim - 1> face_points(deg + 1);
338 for (
unsigned int k = 0; k < face_points.size(); ++k)
348 const unsigned int npoints =
354 for (
unsigned int k = 0;
359 0,
true,
false,
false, this->dofs_per_face));
373 const unsigned int ibase =
375 for (
unsigned int k = 0; k < cell_points.
size(); ++k)
391 #include "fe_bdm.inst" 393 DEAL_II_NAMESPACE_CLOSE
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< std::vector< double > > test_values_cell
void initialize_support_points(const unsigned int bdm_degree)
#define AssertDimension(dim1, dim2)
std::vector< Point< dim > > generalized_support_points
FullMatrix< double > interface_constraints
const unsigned int degree
std::vector< Point< dim-1 > > generalized_face_support_points
const Point< dim > & point(const unsigned int i) const
static std::vector< bool > get_ria_vector(const unsigned int degree)
FE_BDM(const unsigned int p)
unsigned int size() const
static::ExceptionBase & ExcMessage(std::string arg1)
static::ExceptionBase & ExcImpossibleInDim(int arg1)
std::vector< std::vector< FullMatrix< double > > > prolongation
void invert(const FullMatrix< number2 > &M)
void reinit(const TableIndices< N > &new_size, const bool omit_default_initialization=false)
#define Assert(cond, exc)
static::ExceptionBase & ExcDimensionMismatch(std::size_t arg1, std::size_t arg2)
double weight(const unsigned int i) const
virtual std::string get_name() const override
const unsigned int dofs_per_cell
static unsigned int compute_n_pols(const unsigned int n)
static Quadrature< dim > project_to_all_faces(const SubQuadrature &quadrature)
virtual std::unique_ptr< FiniteElement< dim, dim > > clone() const override
std::vector< std::vector< double > > test_values_face
static std::vector< unsigned int > get_dpo_vector(const unsigned int degree)
static unsigned int compute_n_pols(unsigned int degree)
const unsigned int dofs_per_face
void reinit_restriction_and_prolongation_matrices(const bool isotropic_restriction_only=false, const bool isotropic_prolongation_only=false)
void compute(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
FullMatrix< double > inverse_node_matrix
static::ExceptionBase & ExcInternalError()