annulus_triag_mesh_builder.h

#ifndef THESIS_MESH_ANNULUS_TRIAG_MESH_BUILDER_H
#define THESIS_MESH_ANNULUS_TRIAG_MESH_BUILDER_H
 
#define _USE_MATH_DEFINES
 
#include <lf/mesh/mesh.h>
#include <lf/mesh/utils/structured_mesh_builder.h>
 
namespace projects::ipdg_stokes {
 
namespace mesh {
 
class AnnulusTriagMeshBuilder {
 public:
  explicit AnnulusTriagMeshBuilder(
      std::unique_ptr<lf::mesh::MeshFactory> mesh_factory)
      : inner_radius_(0.25),
        outer_radius_(1),
        num_radial_cells_(1),
        num_angular_cells_(4),
        center_point_(Eigen::Vector2d::Zero()),
        mesh_factory_(std::move(mesh_factory)) {
    // Nothing to do here
  }
 
  void setInnerRadius(double r) { inner_radius_ = r; }
 
  void setOuterRadius(double r) { outer_radius_ = r; }
 
  void setCenterPoint(double x, double y) { center_point_ << x, y; }
 
  void setNumRadialCells(lf::base::size_type n) { num_radial_cells_ = n; }
 
  void setNumAngularCells(lf::base::size_type n) { num_angular_cells_ = n; }
 
  std::shared_ptr<lf::mesh::Mesh> Build();
 
 private:
  double inner_radius_;
  double outer_radius_;
  lf::base::size_type num_radial_cells_;
  lf::base::size_type num_angular_cells_;
  Eigen::Vector2d center_point_;
  std::unique_ptr<lf::mesh::MeshFactory> mesh_factory_;
};
 
}  // end namespace mesh
 
}  // end namespace projects::ipdg_stokes
 
#endif  // THESIS_MESH_ANNULUS_TRIAG_MESH_BUILDER_H