tria_o1.h

#ifndef INCGb84bb391fa744cdd9159293bfcb5b311
#define INCGb84bb391fa744cdd9159293bfcb5b311
 
#include "geometry_interface.h"
 
namespace lf::geometry {
 
bool assertNonDegenerateTriangle(
    const Eigen::Matrix<double, Eigen::Dynamic, 3>& coords,
    double tol = 1.0E-8);
 
class TriaO1 : public Geometry {
 public:
  explicit TriaO1(Eigen::Matrix<double, Eigen::Dynamic, 3> coords);
 
  [[nodiscard]] dim_t DimLocal() const override { return 2; }
  [[nodiscard]] dim_t DimGlobal() const override { return coords_.rows(); }
  [[nodiscard]] base::RefEl RefEl() const override {
    return base::RefEl::kTria();
  }
  [[nodiscard]] Eigen::MatrixXd Global(
      const Eigen::MatrixXd& local) const override;
 
  [[nodiscard]] Eigen::MatrixXd Jacobian(
      const Eigen::MatrixXd& local) const override {
    return jacobian_.replicate(1, local.cols());
  }
  [[nodiscard]] Eigen::MatrixXd JacobianInverseGramian(
      const Eigen::MatrixXd& local) const override {
    return jacobian_inverse_gramian_.replicate(1, local.cols());
  }
  [[nodiscard]] Eigen::VectorXd IntegrationElement(
      const Eigen::MatrixXd& local) const override {
    return Eigen::VectorXd::Constant(local.cols(), integrationElement_);
  }
  [[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, lf::base::dim_t codim) const override;
 
 private:
  Eigen::Matrix<double, Eigen::Dynamic, 3> coords_;
  Eigen::Matrix<double, Eigen::Dynamic, 2> jacobian_;
  Eigen::Matrix<double, Eigen::Dynamic, 2> jacobian_inverse_gramian_;
  double integrationElement_{0};
};
}  // namespace lf::geometry
 
#endif  // INCGb84bb391fa744cdd9159293bfcb5b311