[Libmesh-users] State variables

[<an...@tn...>]

Hi

How should state variables (e.g. temperature) be defined? I tried to
use a quantity added with system.add_vector("temperature"), but I
cannot get the variable projected on a refined mesh. 

system.add_vector("T");

// here is the refinement

es.reinit();  // should this already do the projection of _other_vectors??
system.project_vector(system.get_vector("T"));

gives a runtime error message

/numeric/libmesh/include/numerics/petsc_vector.h:604: T
PetscVector<T>::operator()(unsigned int) const [with T = Real]:
Assertion `((i >= this->first_local_index()) && (i <
this->last_local_index()))' failed.

make: *** [run] Aborted

Even if this would work, I do not see how I can calculate an updated
state variable (based on the current solution) and store it in the
vector. 

Is this the way to go? Or should I use another (trivial) system for
the state variable (as Ben suggested for calculated output
quantities)? 

Thanks for enligthenment!

Martin