-
Convergence analysis of a primal-dual optimization-by-continuation algorithm
Authors:
Ignace Loris,
Simone Rebegoldi
Abstract:
We present a numerical iterative optimization algorithm for the minimization of a cost function consisting of a linear combination of three convex terms, one of which is differentiable, a second one is prox-simple and the third one is the composition of a linear map and a prox-simple function. The algorithm's special feature lies in its ability to approximate, in a single iteration run, the minimi…
▽ More
We present a numerical iterative optimization algorithm for the minimization of a cost function consisting of a linear combination of three convex terms, one of which is differentiable, a second one is prox-simple and the third one is the composition of a linear map and a prox-simple function. The algorithm's special feature lies in its ability to approximate, in a single iteration run, the minimizers of the cost function for many different values of the parameters determining the relative weight of the three terms in the cost function. A proof of convergence of the algorithm, based on an inexact variable metric approach, is also provided. As a special case, one recovers a generalization of the primal-dual algorithm of Chambolle and Pock, and also of the proximal-gradient algorithm. Finally, we show how it is related to a primal-dual iterative algorithm based on inexact proximal evaluations of the non-smooth terms of the cost function.
△ Less
Submitted 15 November, 2023;
originally announced November 2023.
-
On a fixed-point continuation method for a convex optimization problem
Authors:
Jean-Baptiste Fest,
Tommi Heikkilä,
Ignace Loris,
Ségolène Martin,
Luca Ratti,
Simone Rebegoldi,
Gesa Sarnighausen
Abstract:
We consider a variation of the classical proximal-gradient algorithm for the iterative minimization of a cost function consisting of a sum of two terms, one smooth and the other prox-simple, and whose relative weight is determined by a penalty parameter. This so-called fixed-point continuation method allows one to approximate the problem's trade-off curve, i.e. to compute the minimizers of the cos…
▽ More
We consider a variation of the classical proximal-gradient algorithm for the iterative minimization of a cost function consisting of a sum of two terms, one smooth and the other prox-simple, and whose relative weight is determined by a penalty parameter. This so-called fixed-point continuation method allows one to approximate the problem's trade-off curve, i.e. to compute the minimizers of the cost function for a whole range of values of the penalty parameter at once. The algorithm is shown to converge, and a rate of convergence of the cost function is also derived. Furthermore, it is shown that this method is related to iterative algorithms constructed on the basis of the $ε$-subdifferential of the prox-simple term. Some numerical examples are provided.
△ Less
Submitted 23 December, 2022;
originally announced December 2022.
-
Primal-dual splitting scheme with backtracking for handling with epigraphic constraint and sparse analysis regularization
Authors:
Laurence Denneulin,
Nelly Pustelnik,
Maud Langlois,
Ignace Loris,
Éric Thiébaut
Abstract:
The convergence of many proximal algorithms involving a gradient descent relies on its Lipschitz constant. To avoid computing it, backtracking rules can be used. While such a rule has already been designed for the forward-backward algorithm (FBwB), this scheme is not flexible enough when a non-differentiable penalization with a linear operator is added to a constraint. In this work we propose a ba…
▽ More
The convergence of many proximal algorithms involving a gradient descent relies on its Lipschitz constant. To avoid computing it, backtracking rules can be used. While such a rule has already been designed for the forward-backward algorithm (FBwB), this scheme is not flexible enough when a non-differentiable penalization with a linear operator is added to a constraint. In this work we propose a backtracking rule for the primal-dual scheme (PDwB), and evaluate its performance for the epigraphical constrained high dynamical reconstruction in high contrast polarimetric imaging, under TV penalization.
△ Less
Submitted 8 September, 2020;
originally announced September 2020.
-
On starting and stopping criteria for nested primal-dual iterations
Authors:
Jixin Chen,
Ignace Loris
Abstract:
The importance of an adequate inner loop starting point (as opposed to a sufficient inner loop stopping rule) is discussed in the context of a numerical optimization algorithm consisting of nested primal-dual proximal-gradient iterations. While the number of inner iterations is fixed in advance, convergence of the whole algorithm is still guaranteed by virtue of a warm-start strategy for the inner…
▽ More
The importance of an adequate inner loop starting point (as opposed to a sufficient inner loop stopping rule) is discussed in the context of a numerical optimization algorithm consisting of nested primal-dual proximal-gradient iterations. While the number of inner iterations is fixed in advance, convergence of the whole algorithm is still guaranteed by virtue of a warm-start strategy for the inner loop, showing that inner loop "starting rules" can be just as effective as "stopping rules" for guaranteeing convergence. The algorithm itself is applicable to the numerical solution of convex optimization problems defined by the sum of a differentiable term and two possibly non-differentiable terms. One of the latter terms should take the form of the composition of a linear map and a proximable function, while the differentiable term needs an accessible gradient. The algorithm reduces to the classical proximal gradient algorithm in certain special cases and it also generalizes other existing algorithms. In addition, under some conditions of strong convexity, we show a linear rate of convergence.
△ Less
Submitted 20 June, 2018;
originally announced June 2018.
-
On the convergence of a linesearch based proximal-gradient method for nonconvex optimization
Authors:
Silvia Bonettini,
Ignace Loris,
Federica Porta,
Marco Prato,
Simone Rebegoldi
Abstract:
We consider a variable metric linesearch based proximal gradient method for the minimization of the sum of a smooth, possibly nonconvex function plus a convex, possibly nonsmooth term. We prove convergence of this iterative algorithm to a critical point if the objective function satisfies the Kurdyka-Lojasiewicz property at each point of its domain, under the assumption that a limit point exists.…
▽ More
We consider a variable metric linesearch based proximal gradient method for the minimization of the sum of a smooth, possibly nonconvex function plus a convex, possibly nonsmooth term. We prove convergence of this iterative algorithm to a critical point if the objective function satisfies the Kurdyka-Lojasiewicz property at each point of its domain, under the assumption that a limit point exists. The proposed method is applied to a wide collection of image processing problems and our numerical tests show that our algorithm results to be flexible, robust and competitive when compared to recently proposed approaches able to address the optimization problems arising in the considered applications.
△ Less
Submitted 8 April, 2017; v1 submitted 12 May, 2016;
originally announced May 2016.
-
Variable metric inexact line-search based methods for nonsmooth optimization
Authors:
Silvia Bonettini,
Ignace Loris,
Federica Porta,
Marco Prato
Abstract:
We develop a new proximal-gradient method for minimizing the sum of a differentiable, possibly nonconvex, function plus a convex, possibly non differentiable, function. The key features of the proposed method are the definition of a suitable descent direction, based on the proximal operator associated to the convex part of the objective function, and an Armijo-like rule to determine the step size…
▽ More
We develop a new proximal-gradient method for minimizing the sum of a differentiable, possibly nonconvex, function plus a convex, possibly non differentiable, function. The key features of the proposed method are the definition of a suitable descent direction, based on the proximal operator associated to the convex part of the objective function, and an Armijo-like rule to determine the step size along this direction ensuring the sufficient decrease of the objective function. In this frame, we especially address the possibility of adopting a metric which may change at each iteration and an inexact computation of the proximal point defining the descent direction. For the more general nonconvex case, we prove that all limit points of the iterates sequence are stationary, while for convex objective functions we prove the convergence of the whole sequence to a minimizer, under the assumption that a minimizer exists. In the latter case, assuming also that the gradient of the smooth part of the objective function is Lipschitz, we also give a convergence rate estimate, showing the O(1/k) complexity with respect to the function values. We also discuss verifiable sufficient conditions for the inexact proximal point and we present the results of a numerical experience on a convex total variation based image restoration problem, showing that the proposed approach is competitive with another state-of-the-art method.
△ Less
Submitted 1 June, 2015;
originally announced June 2015.
-
Proceedings of the second "international Traveling Workshop on Interactions between Sparse models and Technology" (iTWIST'14)
Authors:
L. Jacques,
C. De Vleeschouwer,
Y. Boursier,
P. Sudhakar,
C. De Mol,
A. Pizurica,
S. Anthoine,
P. Vandergheynst,
P. Frossard,
C. Bilen,
S. Kitic,
N. Bertin,
R. Gribonval,
N. Boumal,
B. Mishra,
P. -A. Absil,
R. Sepulchre,
S. Bundervoet,
C. Schretter,
A. Dooms,
P. Schelkens,
O. Chabiron,
F. Malgouyres,
J. -Y. Tourneret,
N. Dobigeon
, et al. (42 additional authors not shown)
Abstract:
The implicit objective of the biennial "international - Traveling Workshop on Interactions between Sparse models and Technology" (iTWIST) is to foster collaboration between international scientific teams by disseminating ideas through both specific oral/poster presentations and free discussions. For its second edition, the iTWIST workshop took place in the medieval and picturesque town of Namur in…
▽ More
The implicit objective of the biennial "international - Traveling Workshop on Interactions between Sparse models and Technology" (iTWIST) is to foster collaboration between international scientific teams by disseminating ideas through both specific oral/poster presentations and free discussions. For its second edition, the iTWIST workshop took place in the medieval and picturesque town of Namur in Belgium, from Wednesday August 27th till Friday August 29th, 2014. The workshop was conveniently located in "The Arsenal" building within walking distance of both hotels and town center. iTWIST'14 has gathered about 70 international participants and has featured 9 invited talks, 10 oral presentations, and 14 posters on the following themes, all related to the theory, application and generalization of the "sparsity paradigm": Sparsity-driven data sensing and processing; Union of low dimensional subspaces; Beyond linear and convex inverse problem; Matrix/manifold/graph sensing/processing; Blind inverse problems and dictionary learning; Sparsity and computational neuroscience; Information theory, geometry and randomness; Complexity/accuracy tradeoffs in numerical methods; Sparsity? What's next?; Sparse machine learning and inference.
△ Less
Submitted 9 October, 2014; v1 submitted 2 October, 2014;
originally announced October 2014.
-
An efficient algorithm for structured sparse quantile regression
Authors:
Vahid Nassiri,
Ignace Loris
Abstract:
Quantile regression is studied in combination with a penalty which promotes structured (or group) sparsity. A mixed $\ell_{1,\infty}$-norm on the parameter vector is used to impose structured sparsity on the traditional quantile regression problem. An algorithm is derived to calculate the piece-wise linear solution path of the corresponding minimization problem. A Matlab implementation of the prop…
▽ More
Quantile regression is studied in combination with a penalty which promotes structured (or group) sparsity. A mixed $\ell_{1,\infty}$-norm on the parameter vector is used to impose structured sparsity on the traditional quantile regression problem. An algorithm is derived to calculate the piece-wise linear solution path of the corresponding minimization problem. A Matlab implementation of the proposed algorithm is provided and some applications of the methods are also studied.
△ Less
Submitted 25 February, 2013;
originally announced February 2013.
-
Iterative algorithms for total variation-like reconstructions in seismic tomography
Authors:
Ignace Loris,
Caroline Verhoeven
Abstract:
A qualitative comparison of total variation like penalties (total variation, Huber variant of total variation, total generalized variation, ...) is made in the context of global seismic tomography. Both penalized and constrained formulations of seismic recovery problems are treated. A number of simple iterative recovery algorithms applicable to these problems are described. The convergence speed o…
▽ More
A qualitative comparison of total variation like penalties (total variation, Huber variant of total variation, total generalized variation, ...) is made in the context of global seismic tomography. Both penalized and constrained formulations of seismic recovery problems are treated. A number of simple iterative recovery algorithms applicable to these problems are described. The convergence speed of these algorithms is compared numerically in this setting. For the constrained formulation a new algorithm is proposed and its convergence is proven.
△ Less
Submitted 6 April, 2012; v1 submitted 20 March, 2012;
originally announced March 2012.
-
An iterative algorithm for sparse and constrained recovery with applications to divergence-free current reconstructions in magneto-encephalography
Authors:
Ignace Loris,
Caroline Verhoeven
Abstract:
We propose an iterative algorithm for the minimization of a $\ell_1$-norm penalized least squares functional, under additional linear constraints. The algorithm is fully explicit: it uses only matrix multiplications with the three matrices present in the problem (in the linear constraint, in the data misfit part and in penalty term of the functional). None of the three matrices must be invertible.…
▽ More
We propose an iterative algorithm for the minimization of a $\ell_1$-norm penalized least squares functional, under additional linear constraints. The algorithm is fully explicit: it uses only matrix multiplications with the three matrices present in the problem (in the linear constraint, in the data misfit part and in penalty term of the functional). None of the three matrices must be invertible. Convergence is proven in a finite-dimensional setting. We apply the algorithm to a synthetic problem in magneto-encephalography where it is used for the reconstruction of divergence-free current densities subject to a sparsity promoting penalty on the wavelet coefficients of the current densities. We discuss the effects of imposing zero divergence and of imposing joint sparsity (of the vector components of the current density) on the current density reconstruction.
△ Less
Submitted 15 February, 2012;
originally announced February 2012.
-
Wavelets and wavelet-like transforms on the sphere and their application to geophysical data inversion
Authors:
Frederik J. Simons,
Ignace Loris,
Eugene Brevdo,
Ingrid C. Daubechies
Abstract:
Many flexible parameterizations exist to represent data on the sphere. In addition to the venerable spherical harmonics, we have the Slepian basis, harmonic splines, wavelets and wavelet-like Slepian frames. In this paper we focus on the latter two: spherical wavelets developed for geophysical applications on the cubed sphere, and the Slepian "tree", a new construction that combines a quadratic co…
▽ More
Many flexible parameterizations exist to represent data on the sphere. In addition to the venerable spherical harmonics, we have the Slepian basis, harmonic splines, wavelets and wavelet-like Slepian frames. In this paper we focus on the latter two: spherical wavelets developed for geophysical applications on the cubed sphere, and the Slepian "tree", a new construction that combines a quadratic concentration measure with wavelet-like multiresolution. We discuss the basic features of these mathematical tools, and illustrate their applicability in parameterizing large-scale global geophysical (inverse) problems.
△ Less
Submitted 8 September, 2011;
originally announced September 2011.
-
Solving or resolving global tomographic models with spherical wavelets, and the scale and sparsity of seismic heterogeneity
Authors:
Frederik J. Simons,
Ignace Loris,
Guust Nolet,
Ingrid C. Daubechies,
S. Voronin,
J. S. Judd,
P. A. Vetter,
J. Charlety,
C. Vonesch
Abstract:
We propose a class of spherical wavelet bases for the analysis of geophysical models and forthe tomographic inversion of global seismic data. Its multiresolution character allows for modeling with an effective spatial resolution that varies with position within the Earth. Our procedure is numerically efficient and can be implemented with parallel computing. We discuss two possible types of discret…
▽ More
We propose a class of spherical wavelet bases for the analysis of geophysical models and forthe tomographic inversion of global seismic data. Its multiresolution character allows for modeling with an effective spatial resolution that varies with position within the Earth. Our procedure is numerically efficient and can be implemented with parallel computing. We discuss two possible types of discrete wavelet transforms in the angular dimension of the cubed sphere. We discuss benefits and drawbacks of these constructions and apply them to analyze the information present in two published seismic wavespeed models of the mantle, for the statistics and power of wavelet coefficients across scales. The localization and sparsity properties of wavelet bases allow finding a sparse solution to inverse problems by iterative minimization of a combination of the $\ell_2$ norm of data fit and the $\ell_1$ norm on the wavelet coefficients. By validation with realistic synthetic experiments we illustrate the likely gains of our new approach in future inversions of finite-frequency seismic data and show its readiness for global seismic tomography.
△ Less
Submitted 14 April, 2011;
originally announced April 2011.
-
On a generalization of the iterative soft-thresholding algorithm for the case of non-separable penalty
Authors:
Ignace Loris,
Caroline Verhoeven
Abstract:
An explicit algorithm for the minimization of an $\ell_1$ penalized least squares functional, with non-separable $\ell_1$ term, is proposed. Each step in the iterative algorithm requires four matrix vector multiplications and a single simple projection on a convex set (or equivalently thresholding). Convergence is proven and a 1/N convergence rate is derived for the functional. In the special case…
▽ More
An explicit algorithm for the minimization of an $\ell_1$ penalized least squares functional, with non-separable $\ell_1$ term, is proposed. Each step in the iterative algorithm requires four matrix vector multiplications and a single simple projection on a convex set (or equivalently thresholding). Convergence is proven and a 1/N convergence rate is derived for the functional. In the special case where the matrix in the $\ell_1$ term is the identity (or orthogonal), the algorithm reduces to the traditional iterative soft-thresholding algorithm. In the special case where the matrix in the quadratic term is the identity (or orthogonal), the algorithm reduces to a gradient projection algorithm for the dual problem. By replacing the projection with a simple proximity operator, other convex non-separable penalties than those based on an $\ell_1$-norm can be handled as well.
△ Less
Submitted 4 August, 2011; v1 submitted 6 April, 2011;
originally announced April 2011.
-
Practical error estimates for sparse recovery in linear inverse problems
Authors:
Ignace Loris,
Caroline Verhoeven
Abstract:
The effectiveness of using model sparsity as a priori information when solving linear inverse problems is studied. We investigate the reconstruction quality of such a method in the non-idealized case and compute some typical recovery errors (depending on the sparsity of the desired solution, the number of data, the noise level on the data, and various properties of the measurement matrix); they ar…
▽ More
The effectiveness of using model sparsity as a priori information when solving linear inverse problems is studied. We investigate the reconstruction quality of such a method in the non-idealized case and compute some typical recovery errors (depending on the sparsity of the desired solution, the number of data, the noise level on the data, and various properties of the measurement matrix); they are compared to known theoretical bounds and illustrated on a magnetic tomography example.
△ Less
Submitted 16 July, 2010; v1 submitted 25 August, 2009;
originally announced August 2009.
-
Accelerating gradient projection methods for $\ell_1$-constrained signal recovery by steplength selection rules
Authors:
I. Loris,
M. Bertero,
C. De Mol,
R. Zanella,
L. Zanni
Abstract:
We propose a new gradient projection algorithm that compares favorably with the fastest algorithms available to date for $\ell_1$-constrained sparse recovery from noisy data, both in the compressed sensing and inverse problem frameworks. The method exploits a line-search along the feasible direction and an adaptive steplength selection based on recent strategies for the alternation of the well-k…
▽ More
We propose a new gradient projection algorithm that compares favorably with the fastest algorithms available to date for $\ell_1$-constrained sparse recovery from noisy data, both in the compressed sensing and inverse problem frameworks. The method exploits a line-search along the feasible direction and an adaptive steplength selection based on recent strategies for the alternation of the well-known Barzilai-Borwein rules. The convergence of the proposed approach is discussed and a computational study on both well-conditioned and ill-conditioned problems is carried out for performance evaluations in comparison with five other algorithms proposed in the literature.
△ Less
Submitted 25 February, 2009;
originally announced February 2009.
-
Nonlinear regularization techniques for seismic tomography
Authors:
I. Loris,
H. Douma,
G. Nolet,
I. Daubechies,
C. Regone
Abstract:
The effects of several nonlinear regularization techniques are discussed in the framework of 3D seismic tomography. Traditional, linear, $\ell_2$ penalties are compared to so-called sparsity promoting $\ell_1$ and $\ell_0$ penalties, and a total variation penalty. Which of these algorithms is judged optimal depends on the specific requirements of the scientific experiment. If the correct reproduct…
▽ More
The effects of several nonlinear regularization techniques are discussed in the framework of 3D seismic tomography. Traditional, linear, $\ell_2$ penalties are compared to so-called sparsity promoting $\ell_1$ and $\ell_0$ penalties, and a total variation penalty. Which of these algorithms is judged optimal depends on the specific requirements of the scientific experiment. If the correct reproduction of model amplitudes is important, classical damping towards a smooth model using an $\ell_2$ norm works almost as well as minimizing the total variation but is much more efficient. If gradients (edges of anomalies) should be resolved with a minimum of distortion, we prefer $\ell_1$ damping of Daubechies-4 wavelet coefficients. It has the additional advantage of yielding a noiseless reconstruction, contrary to simple $\ell_2$ minimization (`Tikhonov regularization') which should be avoided. In some of our examples, the $\ell_0$ method produced notable artifacts. In addition we show how nonlinear $\ell_1$ methods for finding sparse models can be competitive in speed with the widely used $\ell_2$ methods, certainly under noisy conditions, so that there is no need to shun $\ell_1$ penalizations.
△ Less
Submitted 18 August, 2010; v1 submitted 26 August, 2008;
originally announced August 2008.
-
On the performance of algorithms for the minimization of $\ell_1$-penalized functionals
Authors:
Ignace Loris
Abstract:
The problem of assessing the performance of algorithms used for the minimization of an $\ell_1$-penalized least-squares functional, for a range of penalty parameters, is investigated. A criterion that uses the idea of `approximation isochrones' is introduced. Five different iterative minimization algorithms are tested and compared, as well as two warm-start strategies. Both well-conditioned and…
▽ More
The problem of assessing the performance of algorithms used for the minimization of an $\ell_1$-penalized least-squares functional, for a range of penalty parameters, is investigated. A criterion that uses the idea of `approximation isochrones' is introduced. Five different iterative minimization algorithms are tested and compared, as well as two warm-start strategies. Both well-conditioned and ill-conditioned problems are used in the comparison, and the contrast between these two categories is highlighted.
△ Less
Submitted 8 January, 2009; v1 submitted 22 October, 2007;
originally announced October 2007.
-
L1Packv2: A Mathematica package for minimizing an $\ell_1$-penalized functional
Authors:
Ignace Loris
Abstract:
L1Packv2 is a Mathematica package that contains a number of algorithms that can be used for the minimization of an $\ell_1$-penalized least squares functional. The algorithms can handle a mix of penalized and unpenalized variables. Several instructive examples are given. Also, an implementation that yields an exact output whenever exact data are given is provided.
L1Packv2 is a Mathematica package that contains a number of algorithms that can be used for the minimization of an $\ell_1$-penalized least squares functional. The algorithms can handle a mix of penalized and unpenalized variables. Several instructive examples are given. Also, an implementation that yields an exact output whenever exact data are given is provided.
△ Less
Submitted 20 August, 2008; v1 submitted 19 October, 2007;
originally announced October 2007.
-
Sparse and stable Markowitz portfolios
Authors:
Joshua Brodie,
Ingrid Daubechies,
Christine De Mol,
Domenico Giannone,
Ignace Loris
Abstract:
We consider the problem of portfolio selection within the classical Markowitz mean-variance framework, reformulated as a constrained least-squares regression problem. We propose to add to the objective function a penalty proportional to the sum of the absolute values of the portfolio weights. This penalty regularizes (stabilizes) the optimization problem, encourages sparse portfolios (i.e. portf…
▽ More
We consider the problem of portfolio selection within the classical Markowitz mean-variance framework, reformulated as a constrained least-squares regression problem. We propose to add to the objective function a penalty proportional to the sum of the absolute values of the portfolio weights. This penalty regularizes (stabilizes) the optimization problem, encourages sparse portfolios (i.e. portfolios with only few active positions), and allows to account for transaction costs. Our approach recovers as special cases the no-short-positions portfolios, but does allow for short positions in limited number. We implement this methodology on two benchmark data sets constructed by Fama and French. Using only a modest amount of training data, we construct portfolios whose out-of-sample performance, as measured by Sharpe ratio, is consistently and significantly better than that of the naive evenly-weighted portfolio which constitutes, as shown in recent literature, a very tough benchmark.
△ Less
Submitted 29 May, 2008; v1 submitted 31 July, 2007;
originally announced August 2007.
-
Accelerated Projected Gradient Method for Linear Inverse Problems with Sparsity Constraints
Authors:
I. Daubechies,
M. Fornasier,
I. Loris
Abstract:
Regularization of ill-posed linear inverse problems via $\ell_1$ penalization has been proposed for cases where the solution is known to be (almost) sparse. One way to obtain the minimizer of such an $\ell_1$ penalized functional is via an iterative soft-thresholding algorithm. We propose an alternative implementation to $\ell_1$-constraints, using a gradient method, with projection on $\ell_1$-…
▽ More
Regularization of ill-posed linear inverse problems via $\ell_1$ penalization has been proposed for cases where the solution is known to be (almost) sparse. One way to obtain the minimizer of such an $\ell_1$ penalized functional is via an iterative soft-thresholding algorithm. We propose an alternative implementation to $\ell_1$-constraints, using a gradient method, with projection on $\ell_1$-balls. The corresponding algorithm uses again iterative soft-thresholding, now with a variable thresholding parameter. We also propose accelerated versions of this iterative method, using ingredients of the (linear) steepest descent method. We prove convergence in norm for one of these projected gradient methods, without and with acceleration.
△ Less
Submitted 27 August, 2008; v1 submitted 28 June, 2007;
originally announced June 2007.
-
Tomographic inversion using $\ell_1$-norm regularization of wavelet coefficients
Authors:
Ignace Loris,
Guust Nolet,
Ingrid Daubechies,
F. A. Dahlen
Abstract:
We propose the use of $\ell_1$ regularization in a wavelet basis for the solution of linearized seismic tomography problems $Am=d$, allowing for the possibility of sharp discontinuities superimposed on a smoothly varying background. An iterative method is used to find a sparse solution $m$ that contains no more fine-scale structure than is necessary to fit the data $d$ to within its assigned err…
▽ More
We propose the use of $\ell_1$ regularization in a wavelet basis for the solution of linearized seismic tomography problems $Am=d$, allowing for the possibility of sharp discontinuities superimposed on a smoothly varying background. An iterative method is used to find a sparse solution $m$ that contains no more fine-scale structure than is necessary to fit the data $d$ to within its assigned errors.
△ Less
Submitted 8 August, 2006;
originally announced August 2006.
-
Affine Toda-Sutherland Systems
Authors:
Avinash Khare,
I. Loris,
R. Sasaki
Abstract:
A cross between two well-known integrable multi-particle dynamics, an affine Toda molecule and a Sutherland system, is introduced for any affine root system. Though it is not completely integrable but partially integrable, or quasi exactly solvable, it inherits many remarkable properties from the parents. The equilibrium position is algebraic, i.e. proportional to the Weyl vector. The frequencie…
▽ More
A cross between two well-known integrable multi-particle dynamics, an affine Toda molecule and a Sutherland system, is introduced for any affine root system. Though it is not completely integrable but partially integrable, or quasi exactly solvable, it inherits many remarkable properties from the parents. The equilibrium position is algebraic, i.e. proportional to the Weyl vector. The frequencies of small oscillations near equilibrium are proportional to the affine Toda masses, which are essential ingredients of the exact factorisable S-matrices of affine Toda field theories. Some lower lying frequencies are integer times a coupling constant for which the corresponding exact quantum eigenvalues and eigenfunctions are obtained. An affine Toda-Calogero system, with a corresponding rational potential, is also discussed.
△ Less
Submitted 7 September, 2003;
originally announced September 2003.
-
Classical vs Quantum Mechanics: role of elementary excitations
Authors:
I. Loris,
R. Sasaki
Abstract:
Simple theorems relating a quantum mechanical system to the corresponding classical one at equilibrium and connecting the quantum eigenvalues to the frequencies of normal modes oscillations are presented. Corresponding to each quantum eigenfunction, a ` classical eigenfunction' is associated. Those belonging to `elementary excitations' play an important role.
Simple theorems relating a quantum mechanical system to the corresponding classical one at equilibrium and connecting the quantum eigenvalues to the frequencies of normal modes oscillations are presented. Corresponding to each quantum eigenfunction, a ` classical eigenfunction' is associated. Those belonging to `elementary excitations' play an important role.
△ Less
Submitted 16 October, 2003; v1 submitted 7 August, 2003;
originally announced August 2003.
-
Quantum & Classical Eigenfunctions in Calogero & Sutherland Systems
Authors:
I. Loris,
R. Sasaki
Abstract:
An interesting observation was reported by Corrigan-Sasaki that all the frequencies of small oscillations around equilibrium are " quantised" for Calogero and Sutherland (C-S) systems, typical integrable multi-particle dynamics. We present an analytic proof by applying recent results of Loris-Sasaki. Explicit forms of `classical' and quantum eigenfunctions are presented for C-S systems based on…
▽ More
An interesting observation was reported by Corrigan-Sasaki that all the frequencies of small oscillations around equilibrium are " quantised" for Calogero and Sutherland (C-S) systems, typical integrable multi-particle dynamics. We present an analytic proof by applying recent results of Loris-Sasaki. Explicit forms of `classical' and quantum eigenfunctions are presented for C-S systems based on any root systems.
△ Less
Submitted 9 September, 2003; v1 submitted 7 August, 2003;
originally announced August 2003.
