Thanks John! Calling mesh.prepare_for_use() right after the triangulation step did the trick.
What does mesh.prepare_for_use do ? Thanks for your help. Regards, Xavier On Wed, Apr 20, 2016 at 1:09 PM, John Peterson <jwpeter...@gmail.com> wrote: > > > On Wed, Apr 20, 2016 at 1:54 PM, Praphat Fernandes < > fernandes.prap...@gmail.com> wrote: > >> Hello, >> >> I'm attempting to use a mesh(like the one generated in Misc Example 6, >> using the Triangle Interface) to solve a Poisson problem like Intro. >> Example 3. >> >> Here are the problems I encounter: >> >> 1. Within the assemble_poisson function, >> MeshBase::const_element_iterator el = >> mesh.active_local_elements_begin(); and >> MeshBase::const_element_iterator el = >> mesh.active_local_elements_end(); are the same. >> >> 2. On the other hand I can iterate through the elements by accessing them >> directly through mesh.elem(i) for index i within the correct range. >> > > Hmm... somehow *none* of the elements are active? You can try inserting a > mesh.prepare_for_use() after the triangulation step to see if that makes > any difference. > > If your example code is materially different from misc_ex6, can you send > it to me? I'll take a look... > > -- > John > ------------------------------------------------------------------------------ Find and fix application performance issues faster with Applications Manager Applications Manager provides deep performance insights into multiple tiers of your business applications. It resolves application problems quickly and reduces your MTTR. Get your free trial! https://ad.doubleclick.net/ddm/clk/302982198;130105516;z _______________________________________________ Libmesh-users mailing list Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users
