Computational Methods for O+PDEs