Mihai, > solution_transfer.prepare_for_pure_refinement(); > [.] > solution_transfer.interpolate(q, optimality_solution);
As indicated in the documentation, the function corresponding to prepare_for_pure_refinement() is refine_interpolate(), not interpolate(). For the latter to work you need to call prepare_for_coarsening_and_refinement() prior to refinement, irrespective of the fact that you do not coarsen at all. The internal data structures used by the functions do not match. Best, Tobias _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
