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.
// here is the refinement
es.reinit(); // should this already do the projection of _other_vectors??
gives a runtime error message
PetscVector<T>::operator()(unsigned int) const [with T = Real]:
Assertion `((i >= this->first_local_index()) && (i <
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
Is this the way to go? Or should I use another (trivial) system for
the state variable (as Ben suggested for calculated output
Thanks for enligthenment!