# Updating the symmetric indefinite factorization with applications in a modified Newton's method. [SYMUPD; also PIVIXI and PIV2X2 for pivoting]

## Abstract

In recent years the use of quasi-Newton methods in optimization algorithms has inspired much of the research in an area of numerical linear algebra called updating matrix factorizations. Previous research in this area has been concerned with updating the factorization of a symmetric positive definite matrix. Here, a numerical algorithm is presented for updating the Symmetric Indefinite Factorization of Bunch and Parlett. The algorithm requires only O(n/sup 2/) arithmetic operations to update the factorization of an n x n symmetric matrix when modified by a rank-one matrix. An error analysis of this algorithm is given. Computational results are presented that investigate the timing and accuracy of this algorithm. Another algorithm is presented for the unconstrained minimization of a nonlinear functional. The algorithm is a modification of Newton's method. At points where the Hessian is indefinite the search for the next iterate is conducted along a quadratic curve in the plane spanned by a direction of negative curvature and a gradient-related descent direction. The stopping criteria for this search take into account the second-order derivative information. The result is that the iterates are shown to converge globally to a critical point at which the Hessian is positive semidefinite. Computational results aremore »

- Authors:

- Publication Date:

- Research Org.:
- Argonne National Lab., IL (USA)

- Sponsoring Org.:
- US Energy Research and Development Administration (ERDA)

- OSTI Identifier:
- 7220580

- Report Number(s):
- ANL-77-49

- DOE Contract Number:
- W-31-109-ENG-38

- Resource Type:
- Technical Report

- Resource Relation:
- Other Information: Thesis

- Country of Publication:
- United States

- Language:
- English

- Subject:
- 99 GENERAL AND MISCELLANEOUS//MATHEMATICS, COMPUTING, AND INFORMATION SCIENCE; ALGORITHMS; COMPUTER CODES; P CODES; S CODES; MATRICES; FACTORIZATION; OPTIMIZATION; ALGEBRA; COMPUTER CALCULATIONS; FUNCTIONALS; NONLINEAR PROBLEMS; FUNCTIONS; MATHEMATICAL LOGIC; MATHEMATICS; 990200* - Mathematics & Computers

### Citation Formats

```
Sorensen, D C.
```*Updating the symmetric indefinite factorization with applications in a modified Newton's method. [SYMUPD; also PIVIXI and PIV2X2 for pivoting]*. United States: N. p., 1977.
Web. doi:10.2172/7220580.

```
Sorensen, D C.
```*Updating the symmetric indefinite factorization with applications in a modified Newton's method. [SYMUPD; also PIVIXI and PIV2X2 for pivoting]*. United States. doi:10.2172/7220580.

```
Sorensen, D C. Wed .
"Updating the symmetric indefinite factorization with applications in a modified Newton's method. [SYMUPD; also PIVIXI and PIV2X2 for pivoting]". United States. doi:10.2172/7220580. https://www.osti.gov/servlets/purl/7220580.
```

```
@article{osti_7220580,
```

title = {Updating the symmetric indefinite factorization with applications in a modified Newton's method. [SYMUPD; also PIVIXI and PIV2X2 for pivoting]},

author = {Sorensen, D C},

abstractNote = {In recent years the use of quasi-Newton methods in optimization algorithms has inspired much of the research in an area of numerical linear algebra called updating matrix factorizations. Previous research in this area has been concerned with updating the factorization of a symmetric positive definite matrix. Here, a numerical algorithm is presented for updating the Symmetric Indefinite Factorization of Bunch and Parlett. The algorithm requires only O(n/sup 2/) arithmetic operations to update the factorization of an n x n symmetric matrix when modified by a rank-one matrix. An error analysis of this algorithm is given. Computational results are presented that investigate the timing and accuracy of this algorithm. Another algorithm is presented for the unconstrained minimization of a nonlinear functional. The algorithm is a modification of Newton's method. At points where the Hessian is indefinite the search for the next iterate is conducted along a quadratic curve in the plane spanned by a direction of negative curvature and a gradient-related descent direction. The stopping criteria for this search take into account the second-order derivative information. The result is that the iterates are shown to converge globally to a critical point at which the Hessian is positive semidefinite. Computational results are presented which indicate that the method is promising. 6 figures, 9 tables.},

doi = {10.2172/7220580},

journal = {},

number = ,

volume = ,

place = {United States},

year = {1977},

month = {6}

}