15#ifndef dealii_manifold_lib_h
16#define dealii_manifold_lib_h
26#include <boost/signals2/connection.hpp>
34 namespace MappingQImplementation
37 class InverseQuadraticApproximation;
85template <
int dim,
int spacedim = dim>
100 virtual std::unique_ptr<Manifold<dim, spacedim>>
101 clone()
const override;
152 "Access the center with get_center() instead.")
167 static
Tensor<1, spacedim>
264template <
int dim,
int spacedim = dim>
279 virtual std::unique_ptr<Manifold<dim, spacedim>>
280 clone()
const override;
292 const double w)
const override;
358 "Access the center with get_center() instead.")
376 guess_new_point(const
ArrayView<const
Tensor<1, spacedim>> &directions,
377 const
ArrayView<const
double> &distances,
378 const
ArrayView<const
double> &weights) const;
398 do_get_new_points(const
ArrayView<const
Point<spacedim>> &surrounding_points,
421template <
int dim,
int spacedim = dim>
432 const double tolerance = 1e-10);
443 const double tolerance = 1e-10);
448 virtual std::unique_ptr<Manifold<dim, spacedim>>
449 clone()
const override;
490 get_direction()
const;
499 get_point_on_axis()
const;
505 get_tolerance()
const;
565template <
int dim,
int spacedim = dim>
582 const double eccentricity);
584 virtual std::unique_ptr<Manifold<dim, spacedim>>
585 clone()
const override;
603 push_forward_gradient(
const Point<spacedim> &chart_point)
const override;
609 get_major_axis_direction()
const;
621 get_eccentricity()
const;
670template <
int dim,
int spacedim = dim,
int chartdim = dim>
691 const double tolerance = 1e-10);
708 const double tolerance = 1e-10);
725 const std::string push_forward_expression,
726 const std::string pull_back_expression,
730 const std::string chart_vars =
732 const std::string space_vars =
734 const double tolerance = 1e-10,
735 const double h = 1e-8);
745 virtual std::unique_ptr<Manifold<dim, spacedim>>
746 clone()
const override;
777 push_forward_gradient(
const Point<chartdim> &chart_point)
const override;
867 static const int chartdim = 3;
868 static const int spacedim = 3;
875 TorusManifold(
const double centerline_radius,
const double inner_radius);
880 virtual std::unique_ptr<Manifold<dim, 3>>
881 clone()
const override;
887 pull_back(
const Point<3> &p)
const override;
893 push_forward(
const Point<3> &chart_point)
const override;
899 push_forward_gradient(
const Point<3> &chart_point)
const override;
905 get_centerline_radius()
const;
911 get_inner_radius()
const;
1058template <
int dim,
int spacedim = dim>
1075 virtual std::unique_ptr<Manifold<dim, spacedim>>
1076 clone()
const override;
1148 std::array<unsigned int, 20>
1149 get_possible_cells_around_points(
1160 compute_chart_points(
1212 push_forward_gradient(
1245 std::vector<internal::MappingQImplementation::
1246 InverseQuadraticApproximation<dim, spacedim>>
1258template <
int dim,
int spacedim>
1267template <
int dim,
int spacedim>
1276template <
int dim,
int spacedim>
1285template <
int dim,
int spacedim>
1289 return point_on_axis;
1294template <
int dim,
int spacedim>
1303template <
int dim,
int spacedim>
1312template <
int dim,
int spacedim>
1321template <
int dim,
int spacedim>
1325 return eccentricity;
1334 return centerline_radius;
1343 return inner_radius;
virtual void get_new_points(const ArrayView< const Point< spacedim > > &surrounding_points, const Table< 2, double > &weights, ArrayView< Point< spacedim > > new_points) const override
virtual Tensor< 1, spacedim > get_tangent_vector(const Point< spacedim > &x1, const Point< spacedim > &x2) const override
virtual Point< spacedim > get_intermediate_point(const Point< spacedim > &p1, const Point< spacedim > &p2, const double w) const override
virtual Point< spacedim > get_new_point(const ArrayView< const Point< spacedim > > &surrounding_points, const ArrayView< const double > &weights) const override
const Tensor< 1, spacedim > dxn
const Tensor< 1, spacedim > normal_direction
const Point< spacedim > point_on_axis
const Tensor< 1, spacedim > direction
double get_tolerance() const
const Point< spacedim > & get_point_on_axis() const
const Tensor< 1, spacedim > & get_direction() const
const Tensor< 1, spacedim > & get_major_axis_direction() const
const Point< spacedim > & get_center() const
double get_eccentricity() const
const Point< spacedim > center
const double eccentricity
const Tensor< 1, spacedim > direction
const FunctionParser< spacedim >::ConstMap const_map
const std::string chart_vars
const std::string pull_back_expression
const std::string space_vars
const double finite_difference_step
const std::string push_forward_expression
ObserverPointer< const Function< chartdim >, FunctionManifold< dim, spacedim, chartdim > > push_forward_function
ObserverPointer< const Function< spacedim >, FunctionManifold< dim, spacedim, chartdim > > pull_back_function
std::map< std::string, double > ConstMap
virtual void get_normals_at_vertices(const typename Triangulation< dim, spacedim >::face_iterator &face, FaceVertexNormals &face_vertex_normals) const
std::array< Tensor< 1, spacedim >, GeometryInfo< dim >::vertices_per_face > FaceVertexNormals
static Tensor< 1, spacedim > get_periodicity()
virtual std::unique_ptr< Manifold< dim, spacedim > > clone() const override
virtual DerivativeForm< 1, spacedim, spacedim > push_forward_gradient(const Point< spacedim > &chart_point) const override
const Point< spacedim > center
const Point< spacedim > p_center
virtual Tensor< 1, spacedim > normal_vector(const typename Triangulation< dim, spacedim >::face_iterator &face, const Point< spacedim > &p) const override
virtual Point< spacedim > push_forward(const Point< spacedim > &chart_point) const override
virtual Point< spacedim > pull_back(const Point< spacedim > &space_point) const override
const Point< spacedim > & get_center() const
const Point< spacedim > & get_center() const
double get_centerline_radius() const
double get_inner_radius() const
const Triangulation< dim, spacedim > * triangulation
boost::signals2::connection clear_signal
std::vector< internal::MappingQImplementation::InverseQuadraticApproximation< dim, spacedim > > quadratic_approximation
FlatManifold< dim > chart_manifold
std::vector< bool > coarse_cell_is_flat
#define DEAL_II_NAMESPACE_OPEN
#define DEAL_II_DEPRECATED_WITH_COMMENT(comment)
#define DEAL_II_NAMESPACE_CLOSE
const ::parallel::distributed::Triangulation< dim, spacedim > * triangulation