“…This important specialization facilitates use of an inexact Newton method by eliminating the need to identify and implement the Jacobian. Of course, efficiency of JFNK depends critically on preconditioning the inner Krylov subspace method, and the need to track changes in the Jacobian places a premium on preconditioners with low setup cost; see (Knoll and Keyes, 2004;Anderson et al, 1999;Gropp et al, 2000;Pernice and Tocci, 2001;Chacón and Knoll, 2003;Mousseau, 2004;Knoll et al, 2005;Reynolds et al, 2006;Yeckel et al, 2009). Implementations of inexact Newton methods are available in several high-performance software libraries Heroux et al, 2005;Hindmarsh et al, 2005).…”