CoDoSol is a Matlab code for medium scale constrained nonlinear systems of equations F(x)=0, l< =x < = u, where F: R^n --> R^n, l and u are vectors of dimension n. Non-existent lower and upper bounds, i.e. entries of l and u equal to minus o plus infinity, are allowed. The code is based on an affine scaling trust region algorithm. The algorithm combines Newton method and trust region procedures where the merit function used is the norm of the nonlinear residual. The trust region and the corresponding scaled gradient are defined by suitable diagonal scaling avoiding the problem of running directly into a bound. The trust region problem is approximately solved by a constrained dogleg method. Only strictly feasible iterates are generated. A great flexibility in choosing the scaling matrix is allowed for application dependent purposes. If the problem has sparse Jacobians and a relatively big size, the user can choose to work with sparse memory storage. Then, the Newton step is computed via the built-in Matlab function LU with the syntax for calling the UMFPACK package when Matlab 6.5 or later versions are used. Various input/output options are provided, and we refer to the code itself for further documentation. CodoSol code has non-commercial purposes. We assume no responsibility for any errors that may exist in the provided routines. We welcome questions, comments, recommended changes and bug reports. If you use this code, the authors would appreciate your acknowledging having done so in the reports, publications, theses, etc., resulting from their use citing the accompayning paper: Bellavia S., Macconi M., Pieraccini S. (2012), Constrained dogleg methods for nonlinear systems with simple bounds. In: COMPUTATIONAL OPTIMIZATION AND APPLICATIONS, vol. 53 n. 3, pp. 771-794. - ISSN 0926-6003

CODOSOL: a bound-constrained nonlinear equations solver / S. BELLAVIA; S. PIERACCINI. - ELETTRONICO. - (2012).

CODOSOL: a bound-constrained nonlinear equations solver

BELLAVIA, STEFANIA;
2012

Abstract

CoDoSol is a Matlab code for medium scale constrained nonlinear systems of equations F(x)=0, l< =x < = u, where F: R^n --> R^n, l and u are vectors of dimension n. Non-existent lower and upper bounds, i.e. entries of l and u equal to minus o plus infinity, are allowed. The code is based on an affine scaling trust region algorithm. The algorithm combines Newton method and trust region procedures where the merit function used is the norm of the nonlinear residual. The trust region and the corresponding scaled gradient are defined by suitable diagonal scaling avoiding the problem of running directly into a bound. The trust region problem is approximately solved by a constrained dogleg method. Only strictly feasible iterates are generated. A great flexibility in choosing the scaling matrix is allowed for application dependent purposes. If the problem has sparse Jacobians and a relatively big size, the user can choose to work with sparse memory storage. Then, the Newton step is computed via the built-in Matlab function LU with the syntax for calling the UMFPACK package when Matlab 6.5 or later versions are used. Various input/output options are provided, and we refer to the code itself for further documentation. CodoSol code has non-commercial purposes. We assume no responsibility for any errors that may exist in the provided routines. We welcome questions, comments, recommended changes and bug reports. If you use this code, the authors would appreciate your acknowledging having done so in the reports, publications, theses, etc., resulting from their use citing the accompayning paper: Bellavia S., Macconi M., Pieraccini S. (2012), Constrained dogleg methods for nonlinear systems with simple bounds. In: COMPUTATIONAL OPTIMIZATION AND APPLICATIONS, vol. 53 n. 3, pp. 771-794. - ISSN 0926-6003
2012
Goal 9: Industry, Innovation, and Infrastructure
S. BELLAVIA; S. PIERACCINI
File in questo prodotto:
File Dimensione Formato  
codosol.zip

accesso aperto

Descrizione: Matlab Code
Tipologia: Altro
Licenza: Tutti i diritti riservati
Dimensione 10.07 kB
Formato Unknown
10.07 kB Unknown
codosol.pdf

accesso aperto

Descrizione: Papers describing the code
Tipologia: Pdf editoriale (Version of record)
Licenza: Tutti i diritti riservati
Dimensione 772.16 kB
Formato Adobe PDF
772.16 kB Adobe PDF

I documenti in FLORE sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificatore per citare o creare un link a questa risorsa: https://hdl.handle.net/2158/829293
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus ND
  • ???jsp.display-item.citation.isi??? ND
social impact