[dune-pdelab] Very large .vtu files
conf86 at web.de
conf86 at web.de
Fri Aug 21 23:52:29 CEST 2015
Sorry for my late response.
Thank you, I did it just like that.
If someone is interested, this is how I solved it:
// <<<9>>> time loop
U unew(gfs,0.0); // solution to be computed
while (time<tend-1e-8) {
// do time step
bctype.setTime(time+dt); // compute constraints
cc.clear(); // for this time step
Dune::PDELab::constraints(bctype,gfs,cc);
osm.apply(time,dt,uold,g,unew); // do one time step
if (fmod(time,10*dt)<=1e-3*dt){ // write every 10th step into a graphic file
// graphics
Dune::SubsamplingVTKWriter<GV> vtkwriter(gv,3);
Dune::PDELab::addSolutionToVTKWriter(vtkwriter,gfs,unew);
vtkwriter.write(fn.getName(),Dune::VTK::appendedraw);
fn.increment();
}
uold = unew; // advance time step
time += dt;
}
More information about the dune-pdelab
mailing list