segment_o1.h

#ifndef INCG0ec30f6fc5554ec9bea4d9b83e83ea42
#define INCG0ec30f6fc5554ec9bea4d9b83e83ea42
 
#include "geometry_interface.h"
 
namespace lf::geometry {
 
class SegmentO1 : public Geometry {
 public:
  explicit SegmentO1(Eigen::Matrix<double, Eigen::Dynamic, 2> coords)
      : coords_(std::move(coords)) {}
 
  [[nodiscard]] dim_t DimLocal() const override { return 1; }
  [[nodiscard]] dim_t DimGlobal() const override { return coords_.rows(); }
  [[nodiscard]] base::RefEl RefEl() const override {
    return base::RefEl::kSegment();
  }
  [[nodiscard]] Eigen::MatrixXd Global(
      const Eigen::MatrixXd& local) const override;
  [[nodiscard]] Eigen::MatrixXd Jacobian(
      const Eigen::MatrixXd& local) const override;
  [[nodiscard]] Eigen::MatrixXd JacobianInverseGramian(
      const Eigen::MatrixXd& local) const override;
  [[nodiscard]] Eigen::VectorXd IntegrationElement(
      const Eigen::MatrixXd& local) const override;
  [[nodiscard]] std::unique_ptr<Geometry> SubGeometry(dim_t codim,
                                                      dim_t i) const override;
  [[nodiscard]] std::vector<std::unique_ptr<Geometry>> ChildGeometry(
      const RefinementPattern& ref_pat, base::dim_t codim) const override;
 
 private:
  Eigen::Matrix<double, Eigen::Dynamic, 2> coords_;
};
 
}  // namespace lf::geometry
 
#endif  // INCG0ec30f6fc5554ec9bea4d9b83e83ea42