adaptivegridrestart.hh

Go to the documentation of this file.

// -*- mode: C++; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 4 -*-
// vi: set et ts=4 sw=4 sts=4:
//
// SPDX-FileCopyrightInfo: Copyright © DuMux Project contributors, see AUTHORS.md in root folder
// SPDX-License-Identifier: GPL-3.0-or-later
//
#ifndef DUMUX_ADAPTIVEGRIDRESTART_HH
#define DUMUX_ADAPTIVEGRIDRESTART_HH
 
 
#if HAVE_DUNE_ALUGRID
#include <dune/alugrid/grid.hh>
#endif
 
#include <dune/grid/common/backuprestore.hh>
 
#include <dumux/common/exceptions.hh>
#include <dumux/common/parameters.hh>
 
namespace Dumux {
 
template<class Grid>
struct GridRestartCheck
{
    static const bool allowRestart = false;
};
 
// the specializations for grid managers that support restart
#if HAVE_DUNE_ALUGRID
template<int dim, int dimworld, Dune::ALUGridElementType elType, Dune::ALUGridRefinementType refinementType>
struct GridRestartCheck<Dune::ALUGrid<dim, dimworld, elType, refinementType> >
{
    static const bool allowRestart = true;
};
#endif
 
 
template <class Grid, bool allowGridRestart = GridRestartCheck<Grid>::allowRestart >
class AdaptiveGridRestart
{
public:
    template<class Problem>
    static void serializeGrid(Problem& problem)
    {
        DUNE_THROW(Dune::NotImplemented,
                   "Adaptive restart functionality currently only works for dune-ALUGrid.");
    }
 
    template<class Problem>
    static void restartGrid(Problem& problem)
    {
        DUNE_THROW(Dune::NotImplemented,
                   "Adaptive restart functionality currently only works for dune-ALUGrid.");
    }
};
 
template <class Grid>
class AdaptiveGridRestart<Grid, true>
{
public:
    template<class Problem>
    static void serializeGrid(Problem& problem)
    {
        std::string gridName = restartGridFileName_(problem);
#if HAVE_DUNE_ALUGRID
        Dune::BackupRestoreFacility<Grid>::backup(problem.grid(), gridName);
#else
        problem.grid().template writeGrid(gridName, problem.timeManager().time()
                                                    + problem.timeManager().timeStepSize());
#endif
    }
 
    template<class Problem>
    static void restartGrid(Problem& problem)
    {}
 
private:
    template<class Problem>
    static const std::string restartGridFileName_(Problem& problem)
    {
        int rank = problem.gridView().comm().rank();
        std::ostringstream oss;
        try {
            std::string name = getParam<std::string>("Problem.Name");
            oss << name;
        }
        catch (ParameterException &e)
        {
            std::cerr << e.what() << std::endl;
            std::cerr << "Taking name from problem.name(): " << problem.name() << std::endl;
            std::cerr << "Be sure to provide a parameter Problem.Name if you want to restart." << std::endl;
            oss << problem.name();
        }
        oss << "_time=" << problem.timeManager().time() + problem.timeManager().timeStepSize()
            << "_rank=" << rank << ".grs";
        return oss.str();
    }
};
 
 
} // namespace Dumux
#endif