petsc-3.7.1 2016-05-15
MATMFFD_WP
Implements an alternative approach for computing the differencing parameter h used with the finite difference based matrix-free Jacobian. This code implements the strategy of M. Pernice and H. Walker: h = error_rel * sqrt(1 + ||U||) / ||a|| Many br
Notes
1) || U || does not change between linear iterations so is reused Many br2) In GMRES || a || == 1 and so does not need to ever be computed except at restart Many brwhen it is recomputed. Many br
Reference: M. Pernice and H. F. Walker, "NITSOL: A Newton Iterative Many brSolver for Nonlinear Systems", SIAM J. Sci. Stat. Comput.", 1998, Many brvol 19, pp. 302--318. Many br
Options Database Keys
-mat_mffd_compute_normu -Compute the norm of u everytime see MatMFFDWPSetComputeNormU() Many br
Many br
Notes: Requires no global collectives when used with GMRES Many br
Formula used
F'(u)*a = [F(u+h*a) - F(u)]/h where Many br
See Also
MATMFFD, MatCreateMFFD(), MatCreateSNESMF(), MATMFFD_DS
Level:intermediate
Location:src/mat/impls/mffd/wp.c
Index of all Mat routines
Table of Contents for all manual pages
Index of all manual pages