On Tue, 24 Mar 2009, Andrea Hawkins wrote: > SerialMesh::point has a check libmesh_assert (i < this->n_nodes()). > > I believe n_nodes() returns the number of nodes local to the processor.
No, n_nodes() is the number of global nodes. I'm actually not sure why we've removed that assert from ParallelMesh, in fact... it's true that ParallelMesh doesn't always know what n_nodes() even *is* (e.g. when refinement has just happened on some remote partition, adding new nodes there) but even when that data's out of sync it should still be good enough to test a semilocal index against. Perhaps I was planning to add a more specific replacement test in mapvector... > But, I want to pass in the global node number which can be larger > than that. That's fine; that function takes a global node number. > In briefly looking at the parallel mesh version, it looks like it is > returning a value from a vector that holds all the global nodes. Did > I not look close enough? Check out mapvector.h, and gaze in awe. Or mild interest, whatever. Anyway, that "vector" is just a shim over std::map, which made it easier to fit a sparse underlying container into the "everything's in vectors" assumptions of the existing Mesh code. --- Roy ------------------------------------------------------------------------------ Apps built with the Adobe(R) Flex(R) framework and Flex Builder(TM) are powering Web 2.0 with engaging, cross-platform capabilities. Quickly and easily build your RIAs with Flex Builder, the Eclipse(TM)based development software that enables intelligent coding and step-through debugging. Download the free 60 day trial. http://p.sf.net/sfu/www-adobe-com _______________________________________________ Libmesh-users mailing list Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users
