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
