On Thu, Dec 12, 2013 at 11:15 AM, Andrew Davis <andrew....@gmail.com> wrote: > Hi, > > I suppose you are right - it would depend on the node number in the mesh. > So in that case, the set dofIndices contains the correct numbers. However, > if I look at the points (the line cout << pt << endl;) they only come from a > quarter of the domain. > > I am building a mesh via build_square (its a unit box). When the parameter > has linear shape functions the points that are printed come from the domain > [0, 1]x[0, 0.25].
In general, I don't think there's any reason that dof "12" should have any relationship to node "12" in the Mesh... If you want to check that, you can try looping over the Nodes and comparing node->id() and node->dof_number(sys_num, var_num, 0); (The 0 applies to LAGRANGE elements.) -- John ------------------------------------------------------------------------------ Rapidly troubleshoot problems before they affect your business. Most IT organizations don't have a clear picture of how application performance affects their revenue. With AppDynamics, you get 100% visibility into your Java,.NET, & PHP application. Start your 15-day FREE TRIAL of AppDynamics Pro! http://pubads.g.doubleclick.net/gampad/clk?id=84349831&iu=/4140/ostg.clktrk _______________________________________________ Libmesh-users mailing list Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users
