read_field_factory.h

Go to the documentation of this file.

// SPDX-FileCopyrightText: © 2026 PRISMS Center at the University of Michigan
// SPDX-License-Identifier: GNU Lesser General Public Version 2.1
 
#pragma once
 
#include <deal.II/base/exceptions.h>
 
#include <prismspf/core/types.h>
 
#include <prismspf/field_input/read_binary.h>
#include <prismspf/field_input/read_field_base.h>
 
#include <prismspf/user_inputs/spatial_discretization.h>
 
#include <prismspf/utilities/utilities.h>
 
#ifdef PRISMS_PF_WITH_VTK
#  include <prismspf/field_input/read_vtk.h>
#endif
 
PRISMS_PF_BEGIN_NAMESPACE

 
enum class Type
{
  ReadUnstructuredVTK,
  ReadBinary
};
 
template <unsigned int dim, typename number>
std::shared_ptr<ReadFieldBase<dim, number>>
create_reader(const InitialConditionFile       &ic_file,
              const SpatialDiscretization<dim> &spatial_discretization)
{
  switch (ic_file.format)
    {
      case InitialConditionFile::DataFormatType::VTKUnstructuredGrid:
#ifdef PRISMS_PF_WITH_VTK
        return std::make_shared<ReadUnstructuredVTK<dim, number>>(ic_file,
                                                                  spatial_discretization);
#else
        AssertThrow(false,
                    dealii::ExcMessage(
                      "You are trying to read a VTK file as an input; however, PRISMS-PF "
                      "was not built with VTK. Please reconfig PRISMS-PF with VTK using "
                      "-D PRISMS_PF_WITH_VTK=ON"));
#endif
      case InitialConditionFile::DataFormatType::FlatBinary:
        return std::make_shared<ReadBinary<dim, number>>(ic_file, spatial_discretization);
      default:
        AssertThrow(false, UnreachableCode());
    }
}
 
PRISMS_PF_END_NAMESPACE