12#include "geometry_interface.h"
39 explicit TriaO2(Eigen::Matrix<double, Eigen::Dynamic, 6> coords);
47 [[nodiscard]] Eigen::MatrixXd
Global(
48 const Eigen::MatrixXd &local)
const override;
49 [[nodiscard]] Eigen::MatrixXd
Jacobian(
50 const Eigen::MatrixXd &local)
const override;
52 const Eigen::MatrixXd &local)
const override;
54 const Eigen::MatrixXd &local)
const override;
56 [[nodiscard]] std::unique_ptr<Geometry>
SubGeometry(dim_t codim,
57 dim_t i)
const override;
59 [[nodiscard]] std::vector<std::unique_ptr<Geometry>>
ChildGeometry(
66 Eigen::Matrix<double, Eigen::Dynamic, 6>
coords_;
72 Eigen::Matrix<double, Eigen::Dynamic, 1>
alpha_;
73 Eigen::Matrix<double, Eigen::Dynamic, 2>
beta_;
74 Eigen::Matrix<double, Eigen::Dynamic, 2>
gamma_;
75 Eigen::Matrix<double, Eigen::Dynamic, 1>
delta_;
Represents a reference element with all its properties.
static constexpr RefEl kTria()
Returns the reference triangle.
Interface class for shape information on a mesh cell in the spirit of parametric finite element metho...
Abstract interface class for encoding topological local refinement
A second-order triangle in the plane or in 3D space.
Eigen::MatrixXd Global(const Eigen::MatrixXd &local) const override
Map a number of points in local coordinates into the global coordinate system.
dim_t DimGlobal() const override
Dimension of the image of this mapping.
Eigen::Matrix< double, Eigen::Dynamic, 6 > coords_
Coordinates of the 6 vertices/midpoints, stored in matrix columns.
Eigen::Matrix< double, Eigen::Dynamic, 2 > beta_
Eigen::MatrixXd JacobianInverseGramian(const Eigen::MatrixXd &local) const override
Evaluate the Jacobian * Inverse Gramian ( ) simultaneously at numPoints.
std::unique_ptr< Geometry > SubGeometry(dim_t codim, dim_t i) const override
Construct a new Geometry() object that describes the geometry of the i-th sub-entity with codimension...
Eigen::Matrix< double, Eigen::Dynamic, 1 > alpha_
Eigen::Matrix< double, Eigen::Dynamic, 2 > gamma_
dim_t DimLocal() const override
Dimension of the domain of this mapping.
std::vector< std::unique_ptr< Geometry > > ChildGeometry(const RefinementPattern &ref_pat, lf::base::dim_t codim) const override
Generate geometry objects for child entities created in the course of refinement.
Eigen::Matrix< double, Eigen::Dynamic, 1 > delta_
Eigen::VectorXd IntegrationElement(const Eigen::MatrixXd &local) const override
The integration element (factor appearing in integral transformation formula, see below) at number of...
Eigen::Matrix< double, Eigen::Dynamic, 2 > gamma_x_2_
base::RefEl RefEl() const override
The Reference element that defines the domain of this mapping.
TriaO2(Eigen::Matrix< double, Eigen::Dynamic, 6 > coords)
Constructor building triangle from vertex/midpoint coordinates.
Eigen::MatrixXd Jacobian(const Eigen::MatrixXd &local) const override
Evaluate the jacobian of the mapping simultaneously at numPoints points.
unsigned int dim_t
type for dimensions and co-dimensions and numbers derived from them
Defines the Geometry interface and provides a number of classes that implement this interface + addit...