Ted,
> What is is that I'm doing incorrectly ?
>
> My code looks something like so:
>
> DataOutStack<dim> data_out_stack;
> data_out_stack.declare_data_vector("U", DataOutStack<dim>::dof_vector);
> const std::string filename = "solution-1d-stacked.vtk";
> std::ofstream output (filename.c_str());
> for (timestep_number=1, time=time_step; time<=500; time+=time_step,
> ++timestep_number) {
> ..... // Compute Something
>
> // Now for output
>
> data_out_stack.new_parameter_value(time, time_step);
> data_out_stack.attach_dof_handler(dof_handler);
> data_out_stack.add_data_vector(solution_u, "U");
> data_out_stack.build_patches(2);
> data_out_stack.finish_parameter_value();
> data_out_stack.write_vtk (output);
> }
This looks reasonable, though you appear to be writing to the same output file
in every time step, appending in the second time step what was already there
in the first. It may be that your VTK writer only reads the first data set,
which corresponds to the first time step. It will of course also take unduly
long to write the same data over and over again, and the file should become
massive this way...
So, does moving the write_vtk call out of the loop help?
W.
-------------------------------------------------------------------------
Wolfgang Bangerth email: [email protected]
www: http://www.math.tamu.edu/~bangerth/
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
