PoissonOp
double u(?,?)
double f(?,?)
double h()
double ut(?,?) -> output