IR Tools

a MATLAB package of iterative regularization methods and large-scale test problems

Silvia Gazzola, Per Christian Hansen, James G. Nagy

Research output: Contribution to journalArticle

10 Citations (Scopus)
11 Downloads (Pure)

Abstract

This paper describes a new MATLAB software package of iterative regularization methods and test problems for large-scale linear inverse problems. The software package, called IR TOOLS, serves two related purposes: we provide implementations of a range of iterative solvers, including several recently proposed methods that are not available elsewhere, and we provide a set of large-scale test problems in the form of discretizations of 2D linear inverse problems. The solvers include iterative regularization methods where the regularization is due to the semi-convergence of the iterations, Tikhonov-type formulations where the regularization is explicitly formulated in the form of a regularization term, and methods that can impose bound constraints on the computed solutions. All the iterative methods are implemented in a very flexible fashion that allows the problem’s coefficient matrix to be available as a (sparse) matrix, a function handle, or an object. The most basic call to all of the various iterative methods requires only this matrix and the right hand side vector; if the method uses any special stopping criteria, regularization parameters, etc., then default values are set automatically by the code. Moreover, through the use of an optional input structure, the user can also have full control of any of the algorithm parameters. The test problems represent realistic large-scale problems found in image reconstruction and several other applications. Numerical examples illustrate the various algorithms and test problems available in this package.

Original languageEnglish
Pages (from-to)773-811
Number of pages39
JournalNumerical Algorithms
Volume81
Issue number3
Early online date3 Aug 2018
DOIs
Publication statusPublished - 1 Jul 2019

Keywords

  • Iterative regularization methods
  • Linear inverse problems
  • MATLAB
  • Semi-convergence
  • Test problems

ASJC Scopus subject areas

  • Applied Mathematics

Cite this

IR Tools : a MATLAB package of iterative regularization methods and large-scale test problems. / Gazzola, Silvia; Hansen, Per Christian; Nagy, James G.

In: Numerical Algorithms, Vol. 81, No. 3, 01.07.2019, p. 773-811.

Research output: Contribution to journalArticle

Gazzola, Silvia ; Hansen, Per Christian ; Nagy, James G. / IR Tools : a MATLAB package of iterative regularization methods and large-scale test problems. In: Numerical Algorithms. 2019 ; Vol. 81, No. 3. pp. 773-811.
@article{fe14784690974257b399092d829cbed1,
title = "IR Tools: a MATLAB package of iterative regularization methods and large-scale test problems",
abstract = "This paper describes a new MATLAB software package of iterative regularization methods and test problems for large-scale linear inverse problems. The software package, called IR TOOLS, serves two related purposes: we provide implementations of a range of iterative solvers, including several recently proposed methods that are not available elsewhere, and we provide a set of large-scale test problems in the form of discretizations of 2D linear inverse problems. The solvers include iterative regularization methods where the regularization is due to the semi-convergence of the iterations, Tikhonov-type formulations where the regularization is explicitly formulated in the form of a regularization term, and methods that can impose bound constraints on the computed solutions. All the iterative methods are implemented in a very flexible fashion that allows the problem’s coefficient matrix to be available as a (sparse) matrix, a function handle, or an object. The most basic call to all of the various iterative methods requires only this matrix and the right hand side vector; if the method uses any special stopping criteria, regularization parameters, etc., then default values are set automatically by the code. Moreover, through the use of an optional input structure, the user can also have full control of any of the algorithm parameters. The test problems represent realistic large-scale problems found in image reconstruction and several other applications. Numerical examples illustrate the various algorithms and test problems available in this package.",
keywords = "Iterative regularization methods, Linear inverse problems, MATLAB, Semi-convergence, Test problems",
author = "Silvia Gazzola and Hansen, {Per Christian} and Nagy, {James G.}",
year = "2019",
month = "7",
day = "1",
doi = "10.1007/s11075-018-0570-7",
language = "English",
volume = "81",
pages = "773--811",
journal = "Numerical Algorithms",
issn = "1017-1398",
publisher = "Springer Netherlands",
number = "3",

}

TY - JOUR

T1 - IR Tools

T2 - a MATLAB package of iterative regularization methods and large-scale test problems

AU - Gazzola, Silvia

AU - Hansen, Per Christian

AU - Nagy, James G.

PY - 2019/7/1

Y1 - 2019/7/1

N2 - This paper describes a new MATLAB software package of iterative regularization methods and test problems for large-scale linear inverse problems. The software package, called IR TOOLS, serves two related purposes: we provide implementations of a range of iterative solvers, including several recently proposed methods that are not available elsewhere, and we provide a set of large-scale test problems in the form of discretizations of 2D linear inverse problems. The solvers include iterative regularization methods where the regularization is due to the semi-convergence of the iterations, Tikhonov-type formulations where the regularization is explicitly formulated in the form of a regularization term, and methods that can impose bound constraints on the computed solutions. All the iterative methods are implemented in a very flexible fashion that allows the problem’s coefficient matrix to be available as a (sparse) matrix, a function handle, or an object. The most basic call to all of the various iterative methods requires only this matrix and the right hand side vector; if the method uses any special stopping criteria, regularization parameters, etc., then default values are set automatically by the code. Moreover, through the use of an optional input structure, the user can also have full control of any of the algorithm parameters. The test problems represent realistic large-scale problems found in image reconstruction and several other applications. Numerical examples illustrate the various algorithms and test problems available in this package.

AB - This paper describes a new MATLAB software package of iterative regularization methods and test problems for large-scale linear inverse problems. The software package, called IR TOOLS, serves two related purposes: we provide implementations of a range of iterative solvers, including several recently proposed methods that are not available elsewhere, and we provide a set of large-scale test problems in the form of discretizations of 2D linear inverse problems. The solvers include iterative regularization methods where the regularization is due to the semi-convergence of the iterations, Tikhonov-type formulations where the regularization is explicitly formulated in the form of a regularization term, and methods that can impose bound constraints on the computed solutions. All the iterative methods are implemented in a very flexible fashion that allows the problem’s coefficient matrix to be available as a (sparse) matrix, a function handle, or an object. The most basic call to all of the various iterative methods requires only this matrix and the right hand side vector; if the method uses any special stopping criteria, regularization parameters, etc., then default values are set automatically by the code. Moreover, through the use of an optional input structure, the user can also have full control of any of the algorithm parameters. The test problems represent realistic large-scale problems found in image reconstruction and several other applications. Numerical examples illustrate the various algorithms and test problems available in this package.

KW - Iterative regularization methods

KW - Linear inverse problems

KW - MATLAB

KW - Semi-convergence

KW - Test problems

UR - http://www.scopus.com/inward/record.url?scp=85051649136&partnerID=8YFLogxK

U2 - 10.1007/s11075-018-0570-7

DO - 10.1007/s11075-018-0570-7

M3 - Article

VL - 81

SP - 773

EP - 811

JO - Numerical Algorithms

JF - Numerical Algorithms

SN - 1017-1398

IS - 3

ER -