Mathematics – Numerical Analysis
Scientific paper
2009-05-07
Mathematics
Numerical Analysis
12 pages, 2 figures
Scientific paper
A new Jacobian approximation is developed for use in quasi-Newton methods for solving systems of nonlinear equations. The new hypersecant Jacobian approximation is intended for the special case where the evaluation of the functions whose roots are sought dominates the computation time, and additionally the Jacobian is sparse. One example of such a case is the solution of the discretized transport equation to calculate particle and energy fluxes in a fusion plasma. The hypersecant approximation of the Jacobian is calculated using function values from previous Newton iterations, similarly to the Broyden approximation. Unlike Broyden, the hypersecant Jacobian converges to the finite-difference approximation of the Jacobian. The calculation of the hypersecant Jacobian elements requires solving small, dense linear systems, where the coefficient matrices can be ill-conditioned or even exactly singular. Singular-value decomposition (SVD) is therefore used. Convergence comparisons of the hypersecant method, the standard Broyden method, and the colored finite differencing of the PETSc SNES solver are presented.
Carlsson Johan
Cary John R.
No associations
LandOfFree
The hypersecant Jacobian approximation for quasi-Newton solves of sparse nonlinear systems does not yet have a rating. At this time, there are no reviews or comments for this scientific paper.
If you have personal experience with The hypersecant Jacobian approximation for quasi-Newton solves of sparse nonlinear systems, we encourage you to share that experience with our LandOfFree.com community. Your opinion is very important and The hypersecant Jacobian approximation for quasi-Newton solves of sparse nonlinear systems will most certainly appreciate the feedback.
Profile ID: LFWR-SCP-O-455810