Hi,
> I am trying to add mesh refinement to step-35, by using
> KellyErrorEstimator, However, after several hundreds time
> steps, I got following errors:
> [..]
> dof->first < n_dofs
Looking at your code, this might be the problem:
static std::map<unsigned int, double> bval;
if (reinit_prec)
VectorTools::interpolate_boundary_values (dof_handler_pressure, 3,
ZeroFunction<dim>(), bval);
you need to clear() bval in the if clause before calling
interpolate_..., otherwise you are just adding new boundary values and
an invalid with an index>n_dofs might appear.
--
Timo Heister
http://www.math.tamu.edu/~heister/
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
