On Friday, July 22, 2016 at 4:16:47 PM UTC-5, Daniel Arndt wrote: > > 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. >
OK, I see. I need to build two ConstraintMatrix. There is a little waste here since both have the same mesh and same hanging node constraints. > > 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.
