discontinuous_scalar_reference_finite_element.h

#ifndef PROJECTS_DPG_DISCONTINUOUS_SCALAR_REFERENCE_FINITE_ELEMENT
#define PROJECTS_DPG_DISCONTINUOUS_SCALAR_REFERENCE_FINITE_ELEMENT
 
#include <lf/fe/fe.h>
#include <lf/uscalfe/uscalfe.h>
 
#include <iostream>
#include <typeinfo>
 
#include "dpg.h"
 
namespace projects::dpg {
 
template <typename SCALAR>
class DiscontinuousScalarReferenceFiniteElement
    : public lf::fe::ScalarReferenceFiniteElement<SCALAR> {
 public:
  DiscontinuousScalarReferenceFiniteElement() = default;
 
  DiscontinuousScalarReferenceFiniteElement(
      const DiscontinuousScalarReferenceFiniteElement&) = delete;
  DiscontinuousScalarReferenceFiniteElement(
      DiscontinuousScalarReferenceFiniteElement&&) noexcept = default;
  DiscontinuousScalarReferenceFiniteElement& operator=(
      const DiscontinuousScalarReferenceFiniteElement&) = delete;
  DiscontinuousScalarReferenceFiniteElement& operator=(
      DiscontinuousScalarReferenceFiniteElement&&) noexcept = default;
 
  explicit DiscontinuousScalarReferenceFiniteElement(
      std::shared_ptr<const lf::fe::ScalarReferenceFiniteElement<SCALAR>> cfe)
      : lf::fe::ScalarReferenceFiniteElement<SCALAR>(), cfe_(std::move(cfe)) {}
 
  [[nodiscard]] inline bool isInitialized() const { return (cfe_ != nullptr); }
 
  [[nodiscard]] lf::base::RefEl RefEl() const override {
    LF_ASSERT_MSG(isInitialized(), "Not initialized!");
    return cfe_->RefEl();
  }
 
  [[nodiscard]] unsigned Degree() const override {
    LF_ASSERT_MSG(isInitialized(), "Not Initialized");
    return cfe_->Degree();
  }
 
  [[nodiscard]] size_type NumRefShapeFunctions() const override {
    LF_ASSERT_MSG(isInitialized(), "Not initialized");
    return cfe_->NumRefShapeFunctions();
  }
 
  [[nodiscard]] size_type NumRefShapeFunctions(dim_t codim) const override {
    LF_ASSERT_MSG(isInitialized(), "Not initialized");
    return (codim == 0) ? cfe_->NumRefShapeFunctions() : 0;
  }
 
  [[nodiscard]] size_type NumRefShapeFunctions(
      dim_t codim, sub_idx_t /*subidx*/) const override {
    LF_ASSERT_MSG(isInitialized(), "Not initialized");
    return (codim == 0) ? cfe_->NumRefShapeFunctions() : 0;
  }
 
  [[nodiscard]] Eigen::Matrix<SCALAR, Eigen::Dynamic, Eigen::Dynamic>
  EvalReferenceShapeFunctions(const Eigen::MatrixXd& local) const override {
    LF_ASSERT_MSG(isInitialized(), "Not initialized");
    return cfe_->EvalReferenceShapeFunctions(local);
  }
 
  [[nodiscard]] Eigen::Matrix<SCALAR, Eigen::Dynamic, Eigen::Dynamic>
  GradientsReferenceShapeFunctions(
      const Eigen::MatrixXd& local) const override {
    LF_ASSERT_MSG(isInitialized(), "Not initialized");
    return cfe_->GradientsReferenceShapeFunctions(local);
  }
 
  [[nodiscard]] Eigen::MatrixXd EvaluationNodes() const override {
    LF_ASSERT_MSG(isInitialized(), "Not initialized");
    return cfe_->EvaluationNodes();
  }
 
  [[nodiscard]] size_type NumEvaluationNodes() const override {
    LF_ASSERT_MSG(isInitialized(), "Not initialized");
    return cfe_->NumEvaluationNodes();
  }
 
  [[nodiscard]] Eigen::Matrix<SCALAR, 1, Eigen::Dynamic> NodalValuesToDofs(
      const Eigen::Matrix<SCALAR, 1, Eigen::Dynamic>& nodvals) const override {
    LF_ASSERT_MSG(isInitialized(), "Not initialized");
    return cfe_->NodalValuesToDofs(nodvals);
  }
 
  ~DiscontinuousScalarReferenceFiniteElement() override = default;
 
 private:
  std::shared_ptr<const lf::fe::ScalarReferenceFiniteElement<SCALAR>> cfe_;
};
 
}  // namespace projects::dpg
 
#endif  // PROJECTS_DPG_DISCONTINUOUS_SCALAR_REFERENCE_FINITE_ELEMENT