On 09/27/2016 08:52 AM, Deepak Gupta wrote:
I created a simple example but it worked correctly and is attached (I modified
a test code given earlier my Jean-Paul). The difference between this simple
code and my code stated earlier is that between the two checks of hanging
nodes, the interpolate_boundary_values function is not called in hp_check.cpp.
However, in the piece of code sent earlier, there is this line:
VectorTools::interpolate_boundary_values(dof_handler,
42,
boundary_v,
boundary_values);
This is somehow causing the problem. If I comment this one, the hanging
constraints are correctly. I am trying to put this as well in my simplified
example. In the meantime, may you be you have an idea why this is messing up.
The online documentation is down, so could not check details of the function.
Well, VectorTools::interpolate_boundary_values() doesn't touch the
'constraints' variable, so it really shouldn't have any influence.
When trying to find issues such as this, I usually try to start with something
that doesn't work, and remove stuff as far as I can. So, start with your
non-working program and make it smaller and smaller until the issue goes away.
The part that you removed last was then the offending statement. Put it back
in, and that is your minimal testcase.
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.
