Many questions in science and engineering give rise to linear ill-posed problems, whose solution is known to satisfy box constraints, such as nonnegativity. The solution of discretized versions of these problems is highly sensitive to perturbations in the data, discretization errors, and round-off errors introduced during the computations. It is therefore often beneficial to impose known constraints during the solution process. This paper describes a two-phase algorithm for the solution of large-scale box-constrained linear discrete ill-posed problems. The first phase applies a cascadic multilevel method and imposes the constraints on each level by orthogonal projection. The second phase improves the computed approximate solution on the finest level by an active set method. The latter allows several indices of the active set to be updated simultaneously. This reduces the computational effort significantly, when compared to standard active set methods that update one index at a time. Applications to image restoration are presented

A hybrid multilevel-active set method for large box-constrained discrete ill-posed problems

MORIGI, SERENA;REICHEL, LOTHAR;SGALLARI, FIORELLA
2011

Abstract

Many questions in science and engineering give rise to linear ill-posed problems, whose solution is known to satisfy box constraints, such as nonnegativity. The solution of discretized versions of these problems is highly sensitive to perturbations in the data, discretization errors, and round-off errors introduced during the computations. It is therefore often beneficial to impose known constraints during the solution process. This paper describes a two-phase algorithm for the solution of large-scale box-constrained linear discrete ill-posed problems. The first phase applies a cascadic multilevel method and imposes the constraints on each level by orthogonal projection. The second phase improves the computed approximate solution on the finest level by an active set method. The latter allows several indices of the active set to be updated simultaneously. This reduces the computational effort significantly, when compared to standard active set methods that update one index at a time. Applications to image restoration are presented
2011
S.Morigi; R.Plemmons; L.Reichel; F.Sgallari
File in questo prodotto:
Eventuali allegati, non sono esposti

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/11585/110216
 Attenzione

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

Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 8
  • ???jsp.display-item.citation.isi??? 7
social impact