OATAO - Open Archive Toulouse Archive Ouverte Open Access Week

A conjugate gradient like method for p-norm minimization in functional spaces

Estatico, Claudio and Gratton, Serge and Lenti, Flavia and Titley-Peloquin, David A conjugate gradient like method for p-norm minimization in functional spaces. (2017) Numerische Mathematik, 137 (4). 895-922. ISSN 0029-599X

[img]
Preview
(Document in English)

PDF (Author's version) - Requires a PDF viewer such as GSview, Xpdf or Adobe Acrobat Reader
680kB

Official URL: https://doi.org/10.1007/s00211-017-0893-7

Abstract

We develop an iterative algorithm to recover the minimum p-norm solution of the functional linear equation Ax=b, where A:X⟶Y is a continuous linear operator between the two Banach spaces X=Lp, 1<p<2, and Y=Lr, r>1, with x∈X and b∈Y. The algorithm is conceived within the same framework of the Landweber method for functional linear equations in Banach spaces proposed by Schöpfer et al. (Inverse Probl 22:311–329, 2006). Indeed, the algorithm is based on using, at the n-th iteration, a linear combination of the steepest current “descent functional” A∗J(b−Axn) and the previous descent functional, where J denotes a duality map of the Banach space Y. In this regard, the algorithm can be viewed as a generalization of the classical conjugate gradient method on the normal equations in Hilbert spaces. We demonstrate that the proposed iterative algorithm converges strongly to the minimum p-norm solution of the functional linear equation Ax=b and that it is also a regularization method, by applying the discrepancy principle as stopping rule. According to the geometrical properties of Lp spaces, numerical experiments show that the method is fast, robust in terms of both restoration accuracy and stability, promotes sparsity and reduces the over-smoothness in reconstructing edges and abrupt intensity changes.

Item Type:Article
Audience (journal):International peer-reviewed journal
Uncontrolled Keywords:
Institution:French research institutions > Centre National de la Recherche Scientifique - CNRS (FRANCE)
Université de Toulouse > Institut National Polytechnique de Toulouse - Toulouse INP (FRANCE)
Université de Toulouse > Université Toulouse III - Paul Sabatier - UT3 (FRANCE)
Université de Toulouse > Université Toulouse - Jean Jaurès - UT2J (FRANCE)
Université de Toulouse > Université Toulouse 1 Capitole - UT1 (FRANCE)
Other partners > Università degli Studi di Genova - UNIGE (ITALY)
Other partners > Université McGill (CANADA)
Laboratory name:
Funders:
MIUR: Italian Ministry of Education, University and Research (Italie) - GNCS-INdAM : Gruppo Nazionale per il Calcolo Scientifico - Istituto Nazionale di Alta Matematica (Italie) - Fondation RTRA STAE : Réseau Thématique de Recherche Avancée - Sciences et Technologies pour l’Aéronautique et l’Espace (Toulouse, France)
Statistics:download
Deposited On:04 Jun 2020 13:24

Repository Staff Only: item control page