Ha ha. Dammit, what a silly mistake :-) Thanks Wolfgang! On Friday, August 26, 2016 at 11:08:17 PM UTC+2, Wolfgang Bangerth wrote: > > On 08/25/2016 02:34 AM, Jean-Paul Pelteret wrote: > > Hi Deepak, > > > > Ok, so I've created a minimal working example to demonstrate what you're > > seeing here. I'm not sufficiently familiar with the p-refinement part of > the > > hp functionality/implementation to tell whether this is truly the > expected > > behaviour or not. I see that the examples and description > > < > https://www.dealii.org/8.4.1/doxygen/deal.II/group__constraints.html#ga3eaa31a679484e80c193e74e8a967dc8> > > > > in the documentation mix h- and p-refinement, so are _perhaps_ not > > sufficiently clear to describe the expected outcome of this > configuration. Or > > maybe its a bug arising from the adjacent cells being on the same > refinement > > level (i.e. there are no actual hanging _nodes_). Perhaps someone else > could > > clarify this? > > It's not a bug but a feature :-) In your code, you call > > > // Build constraints > ConstraintMatrix constraints; > DoFTools::make_hanging_node_constraints (dof_handler, > constraints); > constraints.clear(); // ************* > > I suspect what you really wanted to call in the last line is > 'constraints.close()' -- because if you do so, the result is pretty much > what > I think you probably wanted ;-) > > Cheers > W. > > -- > ------------------------------------------------------------------------ > Wolfgang Bangerth email: bange...@colostate.edu > www: http://www.math.colostate.edu/~bangerth/ > >
-- 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.
