On 10/05/2016 09:03 AM, Deepak Gupta wrote:
Since it works for the combination of p=1 and 2, I do not think hanging nodes
are a problem. I looked at the hp-paper as well and do not see something that
might be going wrong related to the hanging nodes. I am still trying to narrow
down the problem and come up with a simplified example, however, I still
thought of raising this question. Based on the above scenario, can someone
figure out what might be going wrong?
I am now aware of anything that would cause this, but it certainly sounds like
a bug. A small testcase would definitely be useful!
In the past, I have debugged problems such as this by interpolating a known
and simple function onto the mesh, and then applying the constraints via
ConstraintMatrix::distribute() to the vector with the interpolant. If the
constraints are wrong, you will see this if you just output this as a field
(with sufficiently high n_subdivisions passed to DataOut::build_patches()).
Best
W.
--
------------------------------------------------------------------------
Wolfgang Bangerth email: [email protected]
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 [email protected].
For more options, visit https://groups.google.com/d/optout.
