“…Generalized Newton methods have been extensively studied for solving piecewise linear systems, such as linear complementarity problems arising from the discretization of American options pricing problems [23] and obstacle problems [16], the discrete HJB equation [2,29] and piecewise linear systems arising in the numerical solution of the free-surface hydrodynamics models [3,4,28]. It has been verified that this type of methods possess a finite termination property, i.e., they are able to find a solution in a finite number of iterations under suitable conditions [2,3,10,12,27,28]. In addition, if the generalized Jacobi matrix is an M-matrix, these methods converge globally.…”