geometry_interface.h

#ifndef INCG7ed6b0d4d9244155819c464fc4eb9bbb
#define INCG7ed6b0d4d9244155819c464fc4eb9bbb
 
#include <lf/base/base.h>
 
#include <memory>
 
#include "refinement_pattern.h"
 
namespace lf::geometry {
 
class Geometry {
 protected:
  Geometry() = default;
  Geometry(const Geometry&) = default;
  Geometry(Geometry&&) = default;
  Geometry& operator=(const Geometry&) = default;
  Geometry& operator=(Geometry&&) = default;
 
 public:
  using dim_t = base::RefEl::dim_t;
 
  [[nodiscard]] virtual dim_t DimLocal() const = 0;
 
  [[nodiscard]] virtual dim_t DimGlobal() const = 0;
 
  [[nodiscard]] virtual base::RefEl RefEl() const = 0;
 
  [[nodiscard]] virtual Eigen::MatrixXd Global(
      const Eigen::MatrixXd& local) const = 0;
 
  [[nodiscard]] virtual Eigen::MatrixXd Jacobian(
      const Eigen::MatrixXd& local) const = 0;
 
  [[nodiscard]] virtual Eigen::MatrixXd JacobianInverseGramian(
      const Eigen::MatrixXd& local) const = 0;
 
  [[nodiscard]] virtual Eigen::VectorXd IntegrationElement(
      const Eigen::MatrixXd& local) const = 0;
 
  [[nodiscard]] virtual std::unique_ptr<Geometry> SubGeometry(
      dim_t codim, dim_t i) const = 0;
 
  [[nodiscard]] virtual std::vector<std::unique_ptr<Geometry>> ChildGeometry(
      const RefinementPattern& ref_pat, lf::base::dim_t codim) const = 0;
 
  [[nodiscard]] virtual bool isAffine() const { return true; }
 
  virtual ~Geometry() = default;
 
};  // class Geometry
 
double Volume(const Geometry& geo);
 
inline Eigen::MatrixXd Corners(const Geometry& geo) {
  return geo.Global(geo.RefEl().NodeCoords());
}
 
}  // namespace lf::geometry
 
#endif  // INCG7ed6b0d4d9244155819c464fc4eb9bbb