From: Roy Stogner <roystgnr@ic...>  20070910 17:40:55

On Mon, 10 Sep 2007, Yujie wrote: > I know the degree of freedom numbering is done elementbyelement. I also > know if one uses high order shape functions (p>1), he will obtain more > solution number than the nodes in the mesh. However, now I use Firstorder > Lagrange shape functions, In one tetrahedron element, to my knowledge, one > varialbe is relevant to one vertex of tetrahedron. Is it right? In this case, you can get the solution corresponding to a single Node with system.current_solution(node.dof_number(system.number(),var_num,0)) Unlike the two methods which I previously suggested to you, this method will fail on nonvertex nodes on nonLagrange elements, on interior nodes on quadratic Lagrange elements, and on nonvertex nodes on nonisoparametric Lagrange elements. If that's what you want, go nuts, but I hope you're beginning to see why I must again turn down your emailed offers to let me debug your application codes for free.  Roy 