12#ifndef DUMUX_RESTART_HH 
   13#define DUMUX_RESTART_HH 
   14#warning "This header is deprecated and will be removed after release 2.10." 
   16#include <dune/common/exceptions.hh> 
   17#include <dune/common/fvector.hh> 
   18#include <dune/istl/bvector.hh> 
   19#include <dune/grid/common/rangegenerators.hh> 
   36    template <
class Gr
idView>
 
   37    static const std::string magicRestartCookie_(
const GridView &gridView)
 
   39        const std::string gridName = 
"blubb"; 
 
   40        const int dim = GridView::dimension;
 
   42        int numVertices = gridView.size(dim);
 
   43        int numElements = gridView.size(0);
 
   44        int numEdges = gridView.size(dim-1);
 
   45        int numCPUs = gridView.comm().size();
 
   46        int rank = gridView.comm().rank();
 
   48        std::ostringstream oss;
 
   49        oss << 
"DuMux restart file: " 
   50            << 
"gridName='"<<gridName<<
"' " 
   51            << 
"numCPUs="<<numCPUs<<
" " 
   52            << 
"myRank="<<rank<<
" " 
   53            << 
"numElements="<<numElements<<
" " 
   54            << 
"numEdges="<<numEdges<<
" " 
   55            << 
"numVertices="<<numVertices;
 
   60    template <
class Gr
idView>
 
   61    static const std::string restartFileName_(
const GridView &gridView,
 
   62                                              const std::string &simName,
 
   65        int rank = gridView.comm().rank();
 
   66        std::ostringstream oss;
 
   67        oss << simName<<
"_time="<<t<<
"_rank="<<rank<<
".drs";
 
   82    template <
class Problem>
 
   85        const std::string magicCookie = magicRestartCookie_(problem.gridView());
 
   86        fileName_ = restartFileName_(problem.gridView(),
 
   88                                     problem.timeManager().time()+problem.timeManager().timeStepSize());
 
   91        outStream_.open(fileName_.c_str());
 
   92        outStream_.precision(20);
 
 
  111        outStream_ << cookie << 
"\n";
 
 
  118    { outStream_ << 
"\n"; }
 
 
  125    template <
int codim, 
class Serializer, 
class Gr
idView>
 
  127                           const GridView &gridView)
 
  129        std::ostringstream oss;
 
  130        oss << 
"Entities: Codim " << codim;
 
  131        std::string cookie = oss.str();
 
  135        for (
const auto& entity : entities(gridView, Dune::Codim<codim>())) {
 
  136            serializer.serializeEntity(outStream_, entity);
 
 
  156    template <
class Problem>
 
  159        fileName_ = restartFileName_(problem.gridView(),
 
  165        inStream_.open(fileName_.c_str());
 
  166        if (!inStream_.good()) {
 
  167            DUNE_THROW(Dune::IOError,
 
  170                       << 
"' could not be opened properly");
 
  174        inStream_.seekg(0, std::ios::end);
 
  175        int pos = inStream_.tellg();
 
  177            DUNE_THROW(Dune::IOError,
 
  182        inStream_.seekg(0, std::ios::beg);
 
  184        const std::string magicCookie =
 
  185            magicRestartCookie_(problem.gridView());
 
 
  205        if (!inStream_.good())
 
  206            DUNE_THROW(Dune::IOError,
 
  207                       "Encountered unexpected EOF in restart file.");
 
  209        std::getline(inStream_, buf);
 
  211            DUNE_THROW(Dune::IOError,
 
  212                       "Could not start section '" << cookie << 
"'");
 
 
  221        std::getline(inStream_, dummy);
 
  222        for (
unsigned int i = 0; i < dummy.length(); ++i) {
 
  223            if (!std::isspace(dummy[i])) {
 
  224                DUNE_THROW(Dune::InvalidStateException,
 
  225                           "Encountered unread values while deserializing");
 
 
  236    template <
int codim, 
class Deserializer, 
class Gr
idView>
 
  238                             const GridView &gridView)
 
  240        std::ostringstream oss;
 
  241        oss << 
"Entities: Codim " << codim;
 
  242        std::string cookie = oss.str();
 
  248        for (
const auto& entity : entities(gridView, Dune::Codim<codim>())) {
 
  249            if (!inStream_.good()) {
 
  250                DUNE_THROW(Dune::IOError,
 
  251                           "Restart file is corrupted");
 
  254            std::getline(inStream_, curLine);
 
  255            std::istringstream curLineStream(curLine);
 
  256            deserializer.deserializeEntity(curLineStream, entity);
 
 
  274                                const std::string directory=
".")
 
  276        DUNE_THROW(Dune::NotImplemented,
 
  277                   "Restart::restartFileList()");
 
 
  282    std::string fileName_;
 
  283    std::ifstream inStream_;
 
  284    std::ofstream outStream_;
 
 
Load or save a state of a model to/from the harddisk.
Definition restart.hh:33
void deserializeBegin(Problem &problem, double t)
Start reading a restart file at a certain simulated time.
Definition restart.hh:157
void serializeSectionBegin(const std::string &cookie)
Start a new section in the serialized output.
Definition restart.hh:109
void deserializeEnd()
Stop reading the restart file.
Definition restart.hh:265
void deserializeEntities(Deserializer &deserializer, const GridView &gridView)
Deserialize all leaf entities of a codim in a grid.
Definition restart.hh:237
void deserializeSectionEnd()
End of a section in the serialized output.
Definition restart.hh:218
std::ostream & serializeStream()
The output stream to write the serialized data.
Definition restart.hh:101
const std::string & fileName() const
Returns the name of the file which is (de-)serialized.
Definition restart.hh:76
void deserializeSectionBegin(const std::string &cookie)
Start reading a new section of the restart file.
Definition restart.hh:203
void serializeEnd()
Finish the restart file.
Definition restart.hh:146
std::istream & deserializeStream()
The input stream to read the data which ought to be deserialized.
Definition restart.hh:195
void serializeSectionEnd()
End of a section in the serialized output.
Definition restart.hh:117
static void restartFileList(std::list< std::string > &fileList, const std::string directory=".")
Returns the list of restart files in the current directory.
Definition restart.hh:273
void serializeBegin(Problem &problem)
Write the current state of the model to disk.
Definition restart.hh:83
void serializeEntities(Serializer &serializer, const GridView &gridView)
Serialize all leaf entities of a codim in a gridView.
Definition restart.hh:126