Benjamin Kirk wrote:
> 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
>
That's it. With /active_local iterators /I get second order convergence
in L2 norm.
Thank's for that.
But, it seems that Example 3 has the same problem since (current
version) uses
*MeshBase*::const_element_iterator el = mesh.elements_begin();
*const* MeshBase::const_element_iterator end_el = mesh.elements_end()
in order to iterate through the elements.
Mladen
-------------------------------------------------------------------------
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
