%{
  extern double timestep (const face vector u);
%}