I am trying to preserve the computed solution on a coarser grid and project it to a refined mesh after performing AMR based on a specified criteria, after performing the solve in a cycle. I see that all of the examples (say step-6), the solution is reset to zero when performing solution.reinit() but I do not want this to happen. Is there a way to use VectorTools::interpolate or VectorTools::project to consistently perform an L2 projection on the locally refined mesh ? My initial tries with these functions did not yield much because I do not store a copy of the old dof_handler. Would this be the way to go ?
I was hoping this would happen internally since it doesn't make sense for the linear solver to start the iteration cycle with a zero guess when one could re-use the computed solution on a coarser mesh to get an estimate of the solution on the finer grid. Any help or pointers would be much appreciated. Also, the examples should probably show this method too because it can be much more efficient than starting with zero guesses. Thanks, Vijay _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
