#include "taylor.icc"
typedef g u_exact;
struct p_exact {
  Float operator() (const point& x) const {
    return - Re*(cos(2*pi*x[0]) + cos(2*pi*x[1]))/4 
           - (!have_kinetic_energy ? 0 : Re*(norm2(u(x))/2 - 0.25));
  }
  p_exact(Float Re1=0, bool have_kinetic_energy1=false)
  : u(), pi(acos(Float(-1.0))), Re(Re1), have_kinetic_energy(have_kinetic_energy1) {}
  u_exact u; const Float pi, Re; bool have_kinetic_energy;
};
