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