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.
