12#ifndef DUMUX_IO_GRID_MANAGER_ALU_HH
13#define DUMUX_IO_GRID_MANAGER_ALU_HH
17#include <dune/alugrid/grid.hh>
18#include <dune/alugrid/dgf.hh>
21#ifndef DUMUX_IO_GRID_MANAGER_BASE_HH
49template<
int dim,
int dimworld, Dune::ALUGr
idElementType elType, Dune::ALUGr
idRefinementType refinementType>
50class GridManager<
Dune::ALUGrid<dim, dimworld, elType, refinementType>>
51:
public GridManagerBase<Dune::ALUGrid<dim, dimworld, elType, refinementType>>
54 using Grid = Dune::ALUGrid<dim, dimworld, elType, refinementType>;
55 using ParentType = GridManagerBase<Grid>;
60 void init(
const std::string& modelParamGroup =
"",
bool adaptiveRestart =
false)
64 if (adaptiveRestart ||
hasParam(
"Restart") ||
hasParam(
"TimeManager.Restart"))
71 std::cerr <<
"Warning: You are using a deprecated restart mechanism. The usage will change in the future.\n";
74 const int rank = Dune::MPIHelper::getCommunication().rank();
76 std::ostringstream oss;
77 oss << name <<
"_time=" << restartTime <<
"_rank=" << rank <<
".grs";
78 std::cout <<
"Restoring an ALUGrid from " << oss.str() << std::endl;
79 ParentType::gridPtr() = std::shared_ptr<Grid>(Dune::BackupRestoreFacility<Grid>::restore(oss.str()));
80 ParentType::loadBalance();
88 ParentType::maybeRefineGrid(modelParamGroup);
89 ParentType::loadBalance();
96 if (elType == Dune::cube)
98 else if (elType == Dune::simplex)
101 DUNE_THROW(Dune::IOError,
"ALUGrid only supports Dune::cube or Dune::simplex as cell type!");
103 ParentType::maybeRefineGrid(modelParamGroup);
104 ParentType::loadBalance();
110 const auto prefix = modelParamGroup.empty() ? modelParamGroup : modelParamGroup +
".";
111 DUNE_THROW(ParameterException,
"Please supply one of the parameters "
112 << prefix +
"Grid.UpperRight"
113 <<
", or a grid file in " << prefix +
"Grid.File");
122 const std::string& modelParamGroup)
125 const std::string extension = ParentType::getFileExtension(fileName);
126 if(extension !=
"dgf" && extension !=
"msh")
127 DUNE_THROW(Dune::IOError,
"Grid type " << Dune::className<Grid>() <<
" only supports DGF (*.dgf) and Gmsh (*.msh) grid files but the specified filename has extension: *."<< extension);
130 if (extension ==
"dgf")
132 ParentType::enableDgfGridPointer_ =
true;
133 ParentType::dgfGridPtr() = Dune::GridPtr<Grid>(fileName.c_str(), Dune::MPIHelper::getCommunicator());
134 ParentType::gridData_ = std::make_shared<typename ParentType::GridData>(ParentType::dgfGridPtr());
136 if (extension ==
"msh")
144 ParentType::enableGmshDomainMarkers_ =
true;
149 std::vector<int> boundaryMarkersInsertionIndex, boundaryMarkers, faceMarkers, elementMarkers;
150 auto gridFactory = std::make_unique<Dune::GridFactory<Grid>>();
152 Dune::GmshReader<Grid>::read(*gridFactory, fileName, boundaryMarkersInsertionIndex, elementMarkers, verbose, boundarySegments);
153 ParentType::gridPtr() = std::shared_ptr<Grid>(gridFactory->createGrid());
156 boundaryMarkers.resize(ParentType::gridPtr()->numBoundarySegments(), 0);
157 faceMarkers.resize(ParentType::gridPtr()->leafGridView().size(1), 0);
158 const auto& indexSet = ParentType::gridPtr()->leafGridView().indexSet();
159 for (
const auto& element : elements(ParentType::
gridPtr()->leafGridView()))
161 for (
const auto& intersection : intersections(ParentType::
gridPtr()->leafGridView(),
element))
163 if (intersection.boundary() && gridFactory->wasInserted(intersection))
165 auto marker = boundaryMarkersInsertionIndex[gridFactory->insertionIndex(intersection)];
166 boundaryMarkers[intersection.boundarySegmentIndex()] = marker;
167 faceMarkers[indexSet.index(
element.template subEntity<1>(intersection.indexInInside()))] = marker;
172 ParentType::gridData_ = std::make_shared<typename ParentType::GridData>(ParentType::gridPtr(), std::move(gridFactory),
173 std::move(elementMarkers), std::move(boundaryMarkers), std::move(faceMarkers));
177 auto gridFactory = std::make_unique<Dune::GridFactory<Grid>>();
179 if (Dune::MPIHelper::getCommunication().rank() == 0)
180 Dune::GmshReader<Grid>::read(*gridFactory, fileName, verbose, boundarySegments);
182 ParentType::gridPtr() = std::shared_ptr<Grid>(gridFactory->createGrid());
190 template <
int dimension,
int dimensionworld, std::enable_if_t<dimension != dimensionworld,
int> = 0>
192 const std::string& modelParamGroup)
194 DUNE_THROW(Dune::IOError,
"ALUGrid currently only supports the creation of structured grids with dimension == dimensionworld. Consider reading in a grid file instead.");
200 template <
int dimension,
int dimensionworld, std::enable_if_t<dimension == dimensionworld,
int> = 0>
202 const std::string& modelParamGroup)
205 if (elType == Dune::cube)
207 else if (elType == Dune::simplex)
210 DUNE_THROW(Dune::IOError,
"ALUGrid only supports Dune::cube or Dune::simplex as cell type!");
219template<
int dim,
int dimworld, Dune::ALUGr
idElementType elType, Dune::ALUGr
idRefinementType refinementType>
220class BoundaryFlag<
Dune::ALUGrid<dim, dimworld, elType, refinementType>>
223 BoundaryFlag() : flag_(-1) {}
225 template<
class Intersection>
226 BoundaryFlag(
const Intersection& i) : flag_(-1)
229 flag_ = i.impl().boundaryId();
240namespace Grid::Capabilities {
244template<
int dim,
int dimworld, Dune::ALUGr
idElementType elType, Dune::ALUGr
idRefinementType refinementType>
245struct MultithreadingSupported<
Dune::ALUGrid<dim, dimworld, elType, refinementType>>
248 static bool eval(
const GV&)
Boundary flag to store e.g. in sub control volume faces.
std::size_t value_type
Definition boundaryflag.hh:39
value_type get() const
Definition boundaryflag.hh:41
Grid Grid
Definition gridmanager_base.hh:57
void makeGridFromFile(const std::string &fileName, const std::string &modelParamGroup)
Definition gridmanager_base.hh:187
void makeStructuredGrid(CellType cellType, const std::string &modelParamGroup)
Definition gridmanager_base.hh:270
std::shared_ptr< Grid > & gridPtr()
Definition gridmanager_base.hh:145
void init(const std::string &modelParamGroup="")
Definition gridmanager_base.hh:63
dune-grid capabilities compatibility layer
Provides a grid manager for all supported grid managers with input file interfaces....
T getParamFromGroup(Args &&... args)
A free function to get a parameter from the parameter tree singleton with a model group.
Definition parameters.hh:149
bool hasParamInGroup(const std::string ¶mGroup, const std::string ¶m)
Check whether a key exists in the parameter tree with a model group prefix.
Definition parameters.hh:165
bool hasParam(const std::string ¶m)
Check whether a key exists in the parameter tree.
Definition parameters.hh:157
Definition common/pdesolver.hh:24
static bool eval(const GV &)
Definition gridcapabilities.hh:69