TNPACK—a truncated Newton minimization package for large-scale problems: II. Implementation examples

We have recently presented a FORTRAN package for solving unconstrained optimization problems by a truncated Newton algorithm. TNPACK is intended to solve problems for which some separability and sparsity-structure information of the Hessian is available. The Newton equations are solved approximately at each step by a Preconditioned Conjugate Gradient method, with adaptations to indefinite systems; the linear system involving the preconditioned is solved by a sparse modified Cholesky factorization based on the Yale Sparse Matrix Package. In this paper we describe implemental ion examples on two standard optimization problams and two real-life applications. Our intent is to aid users in their own applications, to highlight key options and parameters that may require tailoring to the problem and to describe application areas for which TNPACK is most !suited. These examples will illustrate various strategies for formulating preconditioners, applying reordering to them in order to minimize fill-in, enforcing truncation, and dealing with indefinite regions.

[1]  A. H. Sherman,et al.  Efficient implementation of sparse nonsymmetric Gaussian elimination without pivoting , 1975, SIGNUM.

[2]  J. Nocedal Updating Quasi-Newton Matrices With Limited Storage , 1980 .

[3]  Jorge J. Moré,et al.  Algorithm 566: FORTRAN Subroutines for Testing Unconstrained Optimization Software [C5], [E4] , 1981, TOMS.

[4]  Stanley C. Eisenstat,et al.  Algorithms and Data Structures for Sparse Symmetric Gaussian Elimination , 1981 .

[5]  Jorge J. Moré,et al.  Testing Unconstrained Optimization Software , 1981, TOMS.

[6]  Stanley C. Eisenstat,et al.  Yale sparse matrix package I: The symmetric codes , 1982 .

[7]  S. Lifson Potential Energy Functions for Structural Molecular Biology , 1983 .

[8]  Trond Steihaug,et al.  Truncated-newtono algorithms for large-scale unconstrained optimization , 1983, Math. Program..

[9]  P. Gill,et al.  Computing Forward-Difference Intervals for Numerical Optimization , 1983 .

[10]  M. Hartmann Molecular mechanics. Von ULRICH BURKERT und NORMAN L. ALLINGER. ACS Monograph 177. Washington: American Chemical Society 1982. 430 S., US $ 77.95 , 1984 .

[11]  A. Fogelson A MATHEMATICAL MODEL AND NUMERICAL METHOD FOR STUDYING PLATELET ADHESION AND AGGREGATION DURING BLOOD CLOTTING , 1984 .

[12]  S. Nash Preconditioning of Truncated-Newton Methods , 1985 .

[13]  J. Mccammon,et al.  Dynamics of Proteins and Nucleic Acids , 2018 .

[14]  Thomas C. Oppe,et al.  NSPCG (Nonsymmetric Preconditioned Conjugate Gradient) user's guide: Version 1. 0: A package for solving large sparse linear systems by various iterative methods , 1988 .

[15]  Jorge Nocedal,et al.  On the limited memory BFGS method for large scale optimization , 1989, Math. Program..

[16]  Claude Lemaréchal,et al.  Some numerical experiments with variable-storage quasi-Newton algorithms , 1989, Math. Program..

[17]  S. Nash,et al.  Assessing a search direction within a truncated-newton method , 1990 .

[18]  M. Mezei,et al.  A molecular dynamics simulation of a water droplet by the implicit-Euler/Langevin scheme , 1991 .

[19]  Tamar Schlick,et al.  New Approaches to Potential Energy Minimization and Molecular Dynamics Algorithms , 1991, Comput. Chem..

[20]  Tamar Schlick,et al.  TNPACK—A truncated Newton minimization package for large-scale problems: I. Algorithm and usage , 1992, TOMS.

[21]  A. Fogelson,et al.  Truncated newton methods and the modeling of complex immersed elastic structures , 1993 .

[22]  François-Xavier Le Dimet,et al.  Numerical Experience with Limited-Memory Quasi-Newton and Truncated Newton Methods , 1993, SIAM J. Optim..

[23]  Tamar Schlick,et al.  Modified Cholesky Factorizations for Sparse Preconditioners , 1993, SIAM J. Sci. Comput..