Sparse Cholesky Factorization for Solving Nonlinear PDEs via Gaussian Processes