1
const Float pi = acos(Float(-1.0));
2
struct u : std::unary_function<point,Float> {
3
Float operator() (const point& x) const {
6
case 1: return sin(pi*x[0]);
7
case 2: return sin(pi*x[0])*sin(pi*x[1]);
8
default: return sin(pi*x[0])*sin(pi*x[1])*sin(pi*x[2]);
11
u(size_t d1) : d(d1) {}