16 #ifndef dealii_fe_poly_face_h 17 #define dealii_fe_poly_face_h 20 #include <deal.II/base/qprojector.h> 21 #include <deal.II/base/std_cxx14/memory.h> 23 #include <deal.II/fe/fe.h> 26 DEAL_II_NAMESPACE_OPEN
56 template <
class PolynomialType,
57 int dim = PolynomialType::dimension + 1,
88 virtual std::unique_ptr<
98 return std_cxx14::make_unique<InternalData>();
101 std::unique_ptr<typename FiniteElement<dim, spacedim>::InternalDataBase>
112 auto data = std_cxx14::make_unique<InternalData>();
115 const unsigned int n_q_points = quadrature.
size();
118 std::vector<double> values(0);
119 std::vector<
Tensor<1, dim - 1>> grads(0);
120 std::vector<Tensor<2, dim - 1>> grad_grads(0);
121 std::vector<Tensor<3, dim - 1>> empty_vector_of_3rd_order_tensors;
122 std::vector<Tensor<4, dim - 1>> empty_vector_of_4th_order_tensors;
129 values.resize(poly_space.n());
130 data->shape_values.resize(poly_space.n(),
131 std::vector<double>(n_q_points));
132 for (
unsigned int i = 0; i < n_q_points; ++i)
134 poly_space.compute(quadrature.
point(i),
138 empty_vector_of_3rd_order_tensors,
139 empty_vector_of_4th_order_tensors);
141 for (
unsigned int k = 0; k < poly_space.n(); ++k)
142 data->shape_values[k][i] = values[k];
153 return std::move(data);
156 std::unique_ptr<typename FiniteElement<dim, spacedim>::InternalDataBase>
163 &output_data)
const override 179 const ::internal::FEValuesImplementation::MappingRelatedData<dim,
185 &output_data)
const override;
190 const unsigned int face_no,
194 const ::internal::FEValuesImplementation::MappingRelatedData<dim,
200 &output_data)
const override;
203 fill_fe_subface_values(
205 const unsigned int face_no,
206 const unsigned int sub_no,
210 const ::internal::FEValuesImplementation::MappingRelatedData<dim,
216 &output_data)
const override;
252 DEAL_II_NAMESPACE_CLOSE
virtual UpdateFlags requires_update_flags(const UpdateFlags update_flags) const override
std::vector< std::vector< double > > shape_values
const Point< dim > & point(const unsigned int i) const
std::unique_ptr< typename FiniteElement< dim, spacedim >::InternalDataBase > get_face_data(const UpdateFlags update_flags, const Mapping< dim, spacedim > &, const Quadrature< dim-1 > &quadrature,::internal::FEValuesImplementation::FiniteElementRelatedData< dim, spacedim > &) const override
unsigned int size() const
#define Assert(cond, exc)
virtual std::unique_ptr< typename FiniteElement< dim, spacedim >::InternalDataBase > get_data(const UpdateFlags, const Mapping< dim, spacedim > &, const Quadrature< dim > &,::internal::FEValuesImplementation::FiniteElementRelatedData< dim, spacedim > &) const override
Abstract base class for mapping classes.
Second derivatives of shape functions.
PolynomialType poly_space
std::unique_ptr< typename FiniteElement< dim, spacedim >::InternalDataBase > get_subface_data(const UpdateFlags update_flags, const Mapping< dim, spacedim > &mapping, const Quadrature< dim-1 > &quadrature,::internal::FEValuesImplementation::FiniteElementRelatedData< dim, spacedim > &output_data) const override
unsigned int get_degree() const
Shape function gradients.
static::ExceptionBase & ExcNotImplemented()
FE_PolyFace(const PolynomialType &poly_space, const FiniteElementData< dim > &fe_data, const std::vector< bool > &restriction_is_additive_flags)
const std::vector< bool > restriction_is_additive_flags