From: Roy Stogner <roystgnr@ic...> - 2006-07-19 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
You can find this sort of information in the class documentation at
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.