17 #include <deal.II/base/qprojector.h> 18 #include <deal.II/base/tensor_product_polynomials.h> 19 #include <deal.II/base/tensor_product_polynomials_const.h> 20 #include <deal.II/base/tensor_product_polynomials_bubbles.h> 21 #include <deal.II/base/polynomials_p.h> 22 #include <deal.II/base/polynomials_piecewise.h> 23 #include <deal.II/base/polynomials_rannacher_turek.h> 24 #include <deal.II/fe/fe_poly.h> 25 #include <deal.II/fe/fe_values.h> 26 #include <deal.II/fe/fe_poly.templates.h> 28 DEAL_II_NAMESPACE_OPEN
39 const ::internal::FEValuesImplementation::MappingRelatedData<1,2> &mapping_data,
48 const InternalData &fe_data =
static_cast<const InternalData &
> (fe_internal);
54 for (
unsigned int k=0; k<this->dofs_per_cell; ++k)
62 for (
unsigned int k=0; k<this->dofs_per_cell; ++k)
68 for (
unsigned int k=0; k<this->dofs_per_cell; ++k)
69 for (
unsigned int i=0; i<quadrature.
size(); ++i)
70 for (
unsigned int j=0; j<2; ++j)
72 mapping_data.jacobian_pushed_forward_grads[i][j]
78 for (
unsigned int k=0; k<this->dofs_per_cell; ++k)
84 for (
unsigned int k=0; k<this->dofs_per_cell; ++k)
85 correct_third_derivatives(output_data, mapping_data, quadrature.
size(), k);
99 const ::internal::FEValuesImplementation::MappingRelatedData<2,3> &mapping_data,
107 const InternalData &fe_data =
static_cast<const InternalData &
> (fe_internal);
113 for (
unsigned int k=0; k<this->dofs_per_cell; ++k)
121 for (
unsigned int k=0; k<this->dofs_per_cell; ++k)
127 for (
unsigned int k=0; k<this->dofs_per_cell; ++k)
128 for (
unsigned int i=0; i<quadrature.
size(); ++i)
129 for (
unsigned int j=0; j<3; ++j)
131 mapping_data.jacobian_pushed_forward_grads[i][j]
137 for (
unsigned int k=0; k<this->dofs_per_cell; ++k)
143 for (
unsigned int k=0; k<this->dofs_per_cell; ++k)
144 correct_third_derivatives(output_data, mapping_data, quadrature.
size(), k);
157 const ::internal::FEValuesImplementation::MappingRelatedData<1,2> &mapping_data,
167 const InternalData &fe_data =
static_cast<const InternalData &
> (fe_internal);
173 for (
unsigned int k=0; k<this->dofs_per_cell; ++k)
181 for (
unsigned int k=0; k<this->dofs_per_cell; ++k)
187 for (
unsigned int k=0; k<this->dofs_per_cell; ++k)
188 for (
unsigned int i=0; i<quadrature.
size(); ++i)
189 for (
unsigned int j=0; j<2; ++j)
191 mapping_data.jacobian_pushed_forward_grads[i][j]
197 for (
unsigned int k=0; k<this->dofs_per_cell; ++k)
203 for (
unsigned int k=0; k<this->dofs_per_cell; ++k)
204 correct_third_derivatives(output_data, mapping_data, quadrature.
size(), k);
217 const ::internal::FEValuesImplementation::MappingRelatedData<2,3> &mapping_data,
222 const InternalData &fe_data =
static_cast<const InternalData &
> (fe_internal);
228 for (
unsigned int k=0; k<this->dofs_per_cell; ++k)
236 for (
unsigned int k=0; k<this->dofs_per_cell; ++k)
242 for (
unsigned int k=0; k<this->dofs_per_cell; ++k)
243 for (
unsigned int i=0; i<quadrature.
size(); ++i)
244 for (
unsigned int j=0; j<3; ++j)
246 mapping_data.jacobian_pushed_forward_grads[i][j]
252 for (
unsigned int k=0; k<this->dofs_per_cell; ++k)
258 for (
unsigned int k=0; k<this->dofs_per_cell; ++k)
259 correct_third_derivatives(output_data, mapping_data, quadrature.
size(), k);
264 #include "fe_poly.inst" 266 DEAL_II_NAMESPACE_CLOSE
virtual void transform(const ArrayView< const Tensor< 1, dim > > &input, const MappingType type, const typename Mapping< dim, spacedim >::InternalDataBase &internal, const ArrayView< Tensor< 1, spacedim > > &output) const =0
ArrayView< typename std::remove_reference< typename std::iterator_traits< Iterator >::reference >::type > make_array_view(const Iterator begin, const Iterator end)
Third derivatives of shape functions.
#define Assert(cond, exc)
Abstract base class for mapping classes.
Second derivatives of shape functions.
unsigned int size() const
Shape function gradients.
static ::ExceptionBase & ExcInternalError()