88 void write (
double time, Dune::VTK::OutputType type = Dune::VTK::ascii)
91 unsigned int count = timesteps_.size();
92 timesteps_.push_back(time);
96 vtkWriter_->write(Dune::concatPaths(path_,seqName(count)),type);
98 vtkWriter_->pwrite(seqName(count), path_,extendpath_,type);
102 std::ofstream pvdFile;
103 pvdFile.exceptions(std::ios_base::badbit | std::ios_base::failbit |
104 std::ios_base::eofbit);
105 std::string pvdname = name_ +
".pvd";
106 pvdFile.open(pvdname.c_str());
107 pvdFile <<
"<?xml version=\"1.0\"?> \n"
108 <<
"<VTKFile type=\"Collection\" version=\"0.1\" byte_order=\"" << Dune::VTK::getEndiannessString() <<
"\"> \n"
109 <<
"<Collection> \n";
110 for (
unsigned int i=0; i<=count; i++)
113 std::string piecepath;
114 std::string fullname;
117 fullname = vtkWriter_->getSerialPieceName(seqName(i), piecepath);
120 piecepath = Dune::concatPaths(path_, extendpath_);
121 fullname = vtkWriter_->getParallelHeaderName(seqName(i), piecepath, size_);
123 pvdFile <<
"<DataSet timestep=\"" << timesteps_[i]
124 <<
"\" group=\"\" part=\"0\" name=\"\" file=\""
125 << fullname <<
"\"/> \n";
127 pvdFile <<
"</Collection> \n"
128 <<
"</VTKFile> \n" << std::flush;