From: Roy Stogner <roystgnr@ic...>  20060719 00:48:16

On Tue, 18 Jul 2006, Dong Xu wrote: > Can anyone tell me how to get the xyz physical locations of all the nodes in > the assembled system matrices? Node is a subclass of Point which is a subclass of TypeVector, so you can get the physical location of any node n by setting x = n(0), y = n(1), and z = n(2). You can loop over all nodes in the mesh by getting an iterator from Mesh::nodes_begin() and incrementing it until you reach Mesh::nodes_end(). You can find this sort of information in the class documentation at http://libmesh.sourceforge.net/doxygen/classes.php Unless you're solving a single scalar equation with isoparametric Lagrange elements, your system matrix will have multiple degrees of freedom (and/or zero degrees of freedom) on some nodes, and "all the nodes in the assembled system matrices" doesn't make sense.  Roy Stogner 