fixing Jacobi Neumann
Add PDEOperator*ProductFun
lap^2