Indeed Wolfgang is right. With the stated change, the example code of JP
works now. So I need to relate the error in mine.
Thanks Wolfgang and Jean-Paul for the help.
Best
Deepak
On 26-08-16 23:08, 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.
--
Deepak K. Gupta
PhD Candidate/Researcher
Structural Optimization & Mechanics
TU Delft / Precision & Microsystems Engineering
Faculty of Mechanical, Maritime & Materials Engineering (3mE)
Mekelweg 2, 2628 CD Delft, The Netherlands
T +31 (0)15 27 86818
d.k.gu...@tudelft.nl
Room 3mE.G.1.150
--
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.
