geometry_reader.h

Go to the documentation of this file.

#pragma once
 
#include <compile_time_options.h>
 
#include "geometry_common_includes.h"
 
#include <deal.II/grid/grid_in.h>
 
#include <fstream>
 
namespace ryujin
{
  namespace Geometries
  {
    template <int dim>
    class Reader : public Geometry<dim>
    {
    public:
      Reader(const std::string &subsection)
          : Geometry<dim>("reader", subsection)
      {
        filename_ = "ryujin.msh";
        this->add_parameter("filename",
                            filename_,
                            "The mesh file to read in via dealii::GridIn. This "
                            "class supports, among others, reading in Gmsh "
                            "*.msh files, and the *.ucd file format.");
      }
 
      void create_coarse_triangulation(
          dealii::Triangulation<dim> &triangulation) const final
      {
        dealii::GridIn<dim> gridin;
        gridin.attach_triangulation(triangulation);
        gridin.read(filename_);
      }
 
    private:
      std::string filename_;
    };
  } // namespace Geometries
} // namespace ryujin