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. Accompanying 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 / Bellavia, S.; Pieraccini, Sandra. - ELETTRONICO. - (2012).

CODOSOL: a bound-constrained nonlinear equations solver

PIERACCINI, SANDRA
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. Accompanying 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
File in questo prodotto:
File Dimensione Formato  
CoDoSol.gz

accesso aperto

Descrizione: Open software scritto dagli autori e pubblicato sotto licenza 2-clause BSD. Sostituisce il precedente allegato solo per il cambiamento di licenza.
Tipologia: Altro materiale allegato
Licenza: PUBBLICO - Tutti i diritti riservati
Dimensione 9.05 kB
Formato Unknown
9.05 kB Unknown Visualizza/Apri
Pubblicazioni consigliate

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

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/11583/2506460
 Attenzione

Attenzione! I dati visualizzati non sono stati sottoposti a validazione da parte dell'ateneo