void simulated_intrepid_integrate(VectorType& v) { if (v.rank() == 3) v(0,0,0) = 1.0; else v(0,0) = 1.0; }