Re: [Libmesh-users] Finding solution at a node

[Roy S. <roy...@ic...>]

On Mon, 5 Jun 2006, John Peterson wrote:

> David Knezevic writes:
> > var1 = system.current_solution(node->id());
> > var2 = system.current_solution(2*node->id());
> > var3 = system.current_solution(3*node->id());
> >
> > This is just a guess, maybe I'm way off...?
>
> Something like this would probably work for Lagrange elements.
>
> See src/mesh/gmv_io.C line 469....
>
> Although I can't see any obvious problems, I'd still double-check
> everything and be careful :)

I see one obvious problem - the vector "v" being handled in that
function is not a standard libMesh coefficient vector, it's
something that's been processed.  It looks like it's created in
EquationSystems::build_solution_vector(), which loops over elements
and uses FEInterface::nodal_soln() on each.

I'd forgotten nodal_soln() exists, but I think I recall upgrading some
of the non-Lagrange elements to support it correctly.  David, you
might want to give that a look.
---
Roy