petsc-3.7.1 2016-05-15
Report Typos and Errors

MATMFFD_DS

the code for compute the "h" used in the finite difference matrix-free matrix vector product. This code implements the strategy in Dennis and Schnabel, "Numerical Methods for Unconstrained Optimization and Nonlinear Equations".

Options Database Keys

-mat_mffd_umin <umin> see MatMFFDDSSetUmin() - Many br Many br

Notes: Requires 2 norms and 1 inner product, but they are computed together Many brso only one parallel collective operation is needed. See MATMFFD_WP for a method Many br(with GMRES) that requires NO collective operations. Many br

Formula used

F'(u)*a = [F(u+h*a) - F(u)]/h where Many brh = error_rel*u'a/||a||^2 if |u'a| > umin*||a||_{1} Many br= error_rel*umin*sign(u'a)*||a||_{1}/||a||^2 otherwise Many brwhere Many brerror_rel = square root of relative error in function evaluation Many brumin = minimum iterate parameter Many br

See Also

MATMFFD, MatCreateMFFD(), MatCreateSNESMF(), MATMFFD_WP, MatMFFDDSSetUmin()

Level:intermediate
Location:
src/mat/impls/mffd/mffddef.c
Index of all Mat routines
Table of Contents for all manual pages
Index of all manual pages