refinement_pattern.h

 
#ifndef INCG3178e8d1e7bf4366bcb00cdb4ebbf5fb
#define INCG3178e8d1e7bf4366bcb00cdb4ebbf5fb
#include <utility>
 
#include "lf/base/base.h"
 
namespace lf::geometry {
 
class RefinementPattern {
 protected:
  RefinementPattern(const RefinementPattern&) = default;
  RefinementPattern(RefinementPattern&&) = default;
  RefinementPattern& operator=(const RefinementPattern&) = default;
  RefinementPattern& operator=(RefinementPattern&&) = default;
 
 public:
  explicit RefinementPattern(lf::base::RefEl ref_el)
      : ref_el_(std::move(ref_el)), lattice_const_(6) {}
 
  RefinementPattern(lf::base::RefEl ref_el, lf::base::size_type lattice_const)
      : ref_el_(std::move(ref_el)), lattice_const_(lattice_const) {
    // Lattice constant N should be a multiple of six in order to be able to
    // define both the barycenter (lattice coordinates [N/3,N/3,N/3], and
    // the midpoints of edges. (lattice coordinates e.g. [N/2,N/2,0]).
    // Of course all lattice coordinates must be integers and must add up to N.
    LF_VERIFY_MSG(lattice_const % 6 == 0,
                  "Lattice constant should be multiple of 6");
  }
 
  [[nodiscard]] lf::base::RefEl RefEl() const { return ref_el_; }
  [[nodiscard]] lf::base::size_type LatticeConst() const {
    return lattice_const_;
  }
  [[nodiscard]] virtual lf::base::size_type NumChildren(
      lf::base::dim_t codim) const = 0;
  [[nodiscard]] virtual std::vector<
      Eigen::Matrix<int, Eigen::Dynamic, Eigen::Dynamic>>
  ChildPolygons(lf::base::dim_t codim) const = 0;
 
  virtual ~RefinementPattern() = default;
 
 protected:
  // TODO(craffael): make this a pure virtual base class
  // NOLINTNEXTLINE
  lf::base::RefEl ref_el_; 
  // NOLINTNEXTLINE
  lf::base::size_type lattice_const_; 
};
 
bool isParallelogram(
    const Eigen::Matrix<int, Eigen::Dynamic, Eigen::Dynamic>& polygon);
 
}  // namespace lf::geometry
 
#endif  // INCG3178e8d1e7bf4366bcb00cdb4ebbf5fb