Junchao,
>> The interpolated solution should already have the correct (old) boundary >> conditions. >> > Is it true? As you see, I call > constraints.distribute(interpolated_solution) to impose boundary > constraints, which are computed in setup_system() as follows > What I meant to say is that the interpolated vector has the old boundary values if you don't apply additional constraints. On cells that you refine, you can improve the approximation of the boundary values by applying a ConstraintMatrix (at time-dt) and hanging node constraints. For the current time step, you then initialize a different ConstraintMatrix containing both Dirichlet boundary values (at time instead of time-dt) and hanging node constraints and use this when assembling your system. Best, Daniel > -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to dealii+unsubscr...@googlegroups.com. For more options, visit https://groups.google.com/d/optout.
