[Libmesh-users] Strange System::point_value behavior

[Andrew D. <and...@gm...>]

Hi all,

I'm observing some strange behavior using the
System::point_valuefunction.  I have a system of two variables (say
"u" and "v"), each added
using the function

system.add_variable(name, libMeshEnums::Order);

if they have the same libMesh::Enums::Order, there is no problem.  However,
if they have different orders I am noticing that system.point_value(varNum,
point) does not always work properly.

If the first one is lower order:

system.point_value(0, point) is correct
system.point_value(1, point) is wrong

if the first one is higher order:

system.point_value(0, point) is correct
system.point_value(1, point) results in a seg. fault

Does anyone know what I am doing wrong?

Thanks!
Andy