Jacobian-free Newton-Krylov methods with GPU acceleration for computing nonlinear ship wave patterns