17 #ifndef __deal2__mapping_q1_h
18 #define __deal2__mapping_q1_h
21 #include <deal.II/base/derivative_form.h>
22 #include <deal.II/base/config.h>
23 #include <deal.II/base/table.h>
24 #include <deal.II/base/qprojector.h>
25 #include <deal.II/grid/tria_iterator.h>
26 #include <deal.II/dofs/dof_accessor.h>
27 #include <deal.II/fe/mapping.h>
56 template <
int dim,
int spacedim=dim>
143 template <
int rank >
152 template <
int rank >
161 template <
int rank >
202 double shape (
const unsigned int qpoint,
203 const unsigned int shape_nr)
const;
209 double &
shape (
const unsigned int qpoint,
210 const unsigned int shape_nr);
218 const unsigned int shape_nr)
const;
226 const unsigned int shape_nr);
234 const unsigned int shape_nr)
const;
242 const unsigned int shape_nr);
299 std::vector<DerivativeForm<1,dim, spacedim > >
covariant;
341 std::vector<std::vector<Tensor<1,spacedim> > >
aux;
402 std::vector<double> &JxW_values,
407 CellSimilarity::Similarity &cell_similarity)
const;
415 const unsigned int face_no,
419 std::vector<double> &JxW_values,
429 const unsigned int face_no,
430 const unsigned int sub_no,
434 std::vector<double> &JxW_values,
463 const unsigned int n_orig_q_points,
480 const unsigned int n_orig_q_points,
488 const unsigned int npts,
490 const CellSimilarity::Similarity cell_similarity,
499 const unsigned int face_no,
500 const unsigned int subface_no,
501 const unsigned int npts,
503 const std::vector<double> &weights,
506 std::vector<double> &JxW_values,
593 transform_real_to_unit_cell_internal_codim1
773 InternalData &mdata)
const;
782 InternalData &mdata)
const;
791 InternalData &mdata)
const;
800 template <
int dim,
int spacedim=dim>
813 template<
int dim,
int spacedim>
817 const unsigned int shape_nr)
const
827 template<
int dim,
int spacedim>
831 const unsigned int shape_nr)
835 shape_values.size()));
840 template<
int dim,
int spacedim>
844 const unsigned int shape_nr)
const
848 shape_derivatives.size()));
854 template<
int dim,
int spacedim>
858 const unsigned int shape_nr)
862 shape_derivatives.size()));
867 template <
int dim,
int spacedim>
871 const unsigned int shape_nr)
const
875 shape_second_derivatives.size()));
881 template <
int dim,
int spacedim>
885 const unsigned int shape_nr)
889 shape_second_derivatives.size()));
895 template <
int dim,
int spacedim>
908 DEAL_II_NAMESPACE_CLOSE
void compute_data(const UpdateFlags flags, const Quadrature< dim > &quadrature, const unsigned int n_orig_q_points, InternalData &data) const
std::vector< DerivativeForm< 1, dim, spacedim > > contravariant
void compute_fill(const typename Triangulation< dim, spacedim >::cell_iterator &cell, const unsigned int npts, const DataSetDescriptor data_set, const CellSimilarity::Similarity cell_similarity, InternalData &data, std::vector< Point< spacedim > > &quadrature_points) const
virtual void compute_mapping_support_points(const typename Triangulation< dim, spacedim >::cell_iterator &cell, std::vector< Point< spacedim > > &a) const
Tensor< 2, dim > second_derivative(const unsigned int qpoint, const unsigned int shape_nr) const
virtual void fill_fe_subface_values(const typename Triangulation< dim, spacedim >::cell_iterator &cell, const unsigned int face_no, const unsigned int sub_no, const Quadrature< dim-1 > &quadrature, typename Mapping< dim, spacedim >::InternalDataBase &mapping_data, typename std::vector< Point< spacedim > > &quadrature_points, std::vector< double > &JxW_values, typename std::vector< Tensor< 1, spacedim > > &boundary_form, typename std::vector< Point< spacedim > > &normal_vectors) const
QProjector< dim >::DataSetDescriptor DataSetDescriptor
Point< dim > transform_real_to_unit_cell_initial_guess(const std::vector< Point< spacedim > > &vertex, const Point< spacedim > &p) const
virtual Mapping< dim, spacedim > * clone() const
std::vector< std::vector< Tensor< 1, dim > > > unit_tangentials
virtual void fill_fe_values(const typename Triangulation< dim, spacedim >::cell_iterator &cell, const Quadrature< dim > &quadrature, typename Mapping< dim, spacedim >::InternalDataBase &mapping_data, typename std::vector< Point< spacedim > > &quadrature_points, std::vector< double > &JxW_values, std::vector< DerivativeForm< 1, dim, spacedim > > &jacobians, std::vector< DerivativeForm< 2, dim, spacedim > > &jacobian_grads, std::vector< DerivativeForm< 1, spacedim, dim > > &inverse_jacobians, std::vector< Point< spacedim > > &cell_normal_vectors, CellSimilarity::Similarity &cell_similarity) const
virtual UpdateFlags update_once(const UpdateFlags flags) const
virtual bool preserves_vertex_locations() const
std::vector< Point< spacedim > > mapping_support_points
virtual Point< spacedim > transform_unit_to_real_cell(const typename Triangulation< dim, spacedim >::cell_iterator &cell, const Point< dim > &p) const
std::vector< std::vector< Tensor< 1, spacedim > > > aux
InternalData(const unsigned int n_shape_functions)
#define Assert(cond, exc)
std::vector< double > shape_values
virtual std::size_t memory_consumption() const
virtual Mapping< dim, spacedim >::InternalDataBase * get_data(const UpdateFlags, const Quadrature< dim > &quadrature) const
virtual void fill_fe_face_values(const typename Triangulation< dim, spacedim >::cell_iterator &cell, const unsigned int face_no, const Quadrature< dim-1 > &quadrature, typename Mapping< dim, spacedim >::InternalDataBase &mapping_data, typename std::vector< Point< spacedim > > &quadrature_points, std::vector< double > &JxW_values, typename std::vector< Tensor< 1, spacedim > > &boundary_form, typename std::vector< Point< spacedim > > &normal_vectors) const
Triangulation< dim, spacedim >::cell_iterator cell_of_current_support_points
::ExceptionBase & ExcIndexRange(int arg1, int arg2, int arg3)
std::vector< DerivativeForm< 1, dim, spacedim > > covariant
void transform_gradients(const VectorSlice< const std::vector< Tensor< rank, dim > > > input, VectorSlice< std::vector< Tensor< rank, spacedim > > > output, const typename Mapping< dim, spacedim >::InternalDataBase &internal, const MappingType type) const
Point< dim > transform_real_to_unit_cell_internal(const typename Triangulation< dim, spacedim >::cell_iterator &cell, const Point< spacedim > &p, const Point< dim > &initial_p_unit, InternalData &mdata) const
virtual Mapping< dim, spacedim >::InternalDataBase * get_face_data(const UpdateFlags flags, const Quadrature< dim-1 > &quadrature) const
Tensor< 1, dim > derivative(const unsigned int qpoint, const unsigned int shape_nr) const
std::vector< Tensor< 1, dim > > shape_derivatives
std::vector< Tensor< 2, dim > > shape_second_derivatives
void compute_face_data(const UpdateFlags flags, const Quadrature< dim > &quadrature, const unsigned int n_orig_q_points, InternalData &data) const
static const unsigned int n_shape_functions
void compute_fill_face(const typename Triangulation< dim, spacedim >::cell_iterator &cell, const unsigned int face_no, const unsigned int subface_no, const unsigned int npts, const DataSetDescriptor data_set, const std::vector< double > &weights, InternalData &mapping_data, std::vector< Point< spacedim > > &quadrature_points, std::vector< double > &JxW_values, std::vector< Tensor< 1, spacedim > > &boundary_form, std::vector< Point< spacedim > > &normal_vectors) const
void transform_differential_forms(const VectorSlice< const std::vector< DerivativeForm< rank, dim, spacedim > > > input, VectorSlice< std::vector< DerivativeForm< rank, spacedim, spacedim > > > output, const typename Mapping< dim, spacedim >::InternalDataBase &internal, const MappingType type) const
virtual UpdateFlags update_each(const UpdateFlags flags) const
double shape(const unsigned int qpoint, const unsigned int shape_nr) const
void transform_fields(const VectorSlice< const std::vector< Tensor< rank, dim > > > input, VectorSlice< std::vector< Tensor< rank, spacedim > > > output, const typename Mapping< dim, spacedim >::InternalDataBase &internal, const MappingType type) const
unsigned int n_shape_functions
virtual Point< dim > transform_real_to_unit_cell(const typename Triangulation< dim, spacedim >::cell_iterator &cell, const Point< spacedim > &p) const
void compute_shapes(const std::vector< Point< dim > > &unit_points, InternalData &data) const
Point< spacedim > transform_unit_to_real_cell_internal(const InternalData &mdata) const
virtual void compute_shapes_virtual(const std::vector< Point< dim > > &unit_points, InternalData &data) const
virtual Mapping< dim, spacedim >::InternalDataBase * get_subface_data(const UpdateFlags flags, const Quadrature< dim-1 > &quadrature) const