>> I don't seem to see the attachment... But I have to ask:
>>
>> If the assembly function was taken from example 3 (2 does no assembly(?)),
>> are you constraining the hanging degrees of freedom?
>
> That shouldn't be the problem, if he's just doing uniform refinement.
Roy's right, that's not the problem, but I think this is:
MeshBase::const_element_iterator el =
mesh.local_elements_begin();
const MeshBase::const_element_iterator end_el =
mesh.local_elements_end();
for ( ; el != end_el; ++el)
{
...
What happens if you change the loop to
MeshBase::const_element_iterator el =
mesh.active_local_elements_begin();
const MeshBase::const_element_iterator end_el =
mesh.active_local_elements_end();
Note that in the code you sent you are assembling the system for *all* the
elements, including the coarser mesh levels.
What does that do?
-Ben
-------------------------------------------------------------------------
This SF.Net email is sponsored by the Moblin Your Move Developer's challenge
Build the coolest Linux based applications with Moblin SDK & win great prizes
Grand prize is a trip for two to an Open Source event anywhere in the world
http://moblin-contest.org/redirect.php?banner_id=100&url=/
_______________________________________________
Libmesh-users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/libmesh-users
