Inference of Nonlinear Partial Differential Equations via Constrained Gaussian Processes