MeshManager.h

Go to the documentation of this file.

#ifndef MESH_MANAGER_H
#define MESH_MANAGER_H
 
#include <deal.II/grid/tria.h>
#include <deal.II/distributed/tria.h>
#include <deal.II/distributed/fully_distributed_tria.h>
#include <deal.II/grid/grid_generator.h>
#include <deal.II/grid/grid_out.h>
#include <deal.II/grid/grid_in.h>
#include <deal.II/grid/grid_tools.h>
#include <deal.II/base/utilities.h>
#include <deal.II/base/conditional_ostream.h>
#include <deal.II/numerics/data_out.h>
 
#include <filesystem>
#include <vector>
#include <string>
#include <fstream>
#include <stdexcept>
 
using namespace dealii;
 
template <int dim, int spacedim = dim>
class MeshManager {
public:
    MeshManager(const MPI_Comm& comm);
 
    template<typename TriangulationType>
    void generate_mesh(TriangulationType& tria,
                      const std::string& grid_generator_function,
                      const std::string& grid_generator_arguments,
                      const unsigned int n_refinements,
                      const bool use_tetrahedral_mesh);
 
    void load_source_mesh(parallel::fullydistributed::Triangulation<dim, spacedim>& source_mesh,
                         const std::string& source_file = "");
 
    void load_target_mesh(Triangulation<dim, spacedim>& target_mesh);
 
    void load_mesh_at_level(parallel::fullydistributed::Triangulation<dim, spacedim>& source_mesh,
                           DoFHandler<dim, spacedim>& dof_handler_source,
                           const std::string& mesh_file);
 
    void save_meshes(const parallel::fullydistributed::Triangulation<dim, spacedim>& source_mesh,
                    const Triangulation<dim, spacedim>& target_mesh);
 
    template<typename TriangulationType>
    bool write_mesh(const TriangulationType& mesh,
                   const std::string& filepath,
                   const std::vector<std::string>& formats,
                   const std::vector<double>* cell_data = nullptr,
                   const std::string& data_name = "cell_data");
 
    static std::vector<std::string> get_mesh_hierarchy_files(
        const std::string& dir = "output/data_multilevel/source_multilevel");
 
private:
    MPI_Comm mpi_communicator;
    const unsigned int n_mpi_processes;
    const unsigned int this_mpi_process;
    ConditionalOStream pcout;
 
    // Constants
    const std::string mesh_directory = "output/data_mesh";
};
 
// Include template implementations
#include "MeshManager.templates.h"
 
#endif // MESH_MANAGER_H