Author | Title | Year | Journal/Proceedings | Reftype | DOI/URL |
---|---|---|---|---|---|
Vande Velde, S., Piron, A., Leo, O. and Loris, I. | SmartFACS: a computational deconvolution method for blood and tumor samples in humans and mice | 2024 | article | ||
Abstract: An important confounding factor in the analysis of expression data is the cellular composition of the samples. This is often overlooked and can lead to a loss of information from less abundant cell types. The distribution and abundance of immune cells play a central role in cancer immunology and hematology. For this reason, we developed a computational cell sorter called SmartFACS to estimate the percentages of over 20 immune cell types from bulk RNA-seq data. The algorithm at the heart of this program is used for the first time in the field of deconvolution. The reference profiles contained in SmartFACS for the detection of the cell types enable highly accurate analysis of tumor samples as well as blood samples from healthy patients or those suffering from autoimmune diseases. SmartFACS is also suitable for mouse data, which is very important for research applications. A comparison with current deconvolution methods has been carried out to demonstrate the powerful and unique function of SmartFACS in the estimation of tumor infiltrates and the analysis of blood samples. Finally, two concrete applications are presented: one in immuno-oncology and the other in vaccinology. | |||||
BibTeX:
@article{Velde2024, author = {Vande Velde, Sylvie and Anthony Piron and Oberdan Leo and Ignace Loris}, title = {SmartFACS: a computational deconvolution method for blood and tumor samples in humans and mice}, year = {2024}, note = {To be submitted.} } |
|||||
Fest, J.-B., Heikkilä, T., Loris, I., Martin, S., Ratti, L., Rebegoldi, S. and Sarnighausen, G. | On a fixed-point continuation method for a convex optimization problem | 2023 | Advanced Techniques in Optimization for Machine learning and Imaging | inproceedings | URL |
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 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. | |||||
BibTeX:
@inproceedings{Fest2022, author = {Jean-Baptiste Fest and Tommi Heikkilä and Ignace Loris and Ségolène Martin and Luca Ratti and Simone Rebegoldi and Gesa Sarnighausen}, title = {On a fixed-point continuation method for a convex optimization problem}, booktitle = {Advanced Techniques in Optimization for Machine learning and Imaging}, publisher = {Springer}, year = {2023}, note = {Accepted.}, url = {https://arxiv.org/abs/2212.12256} } |
|||||
Loris, I. and Rebegoldi, S. | Convergence analysis of a primal-dual optimization-by-continuation algorithm | 2023 | Journal of Computational and Applied Mathematics, pp. 22 | article | URL |
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 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. | |||||
BibTeX:
@article{Loris2023, author = {Ignace Loris and Simone Rebegoldi}, title = {Convergence analysis of a primal-dual optimization-by-continuation algorithm}, journal = {Journal of Computational and Applied Mathematics}, year = {2023}, pages = {22}, note = {Submitted.}, url = {https://arxiv.org/abs/2311.09123} } |
|||||
Denneulin, L., Pustelnik, N., Langlois, M., Loris, I. and Thiébaut, E. | Primal-dual splitting scheme with backtracking for handling with epigraphic constraint and sparse analysis regularization | 2020 | International Traveling Workshop on Interactions between low-complexity data models and Sensing Techniques (iTWIST 2020) | inproceedings | URL |
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 nondifferentiable 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. | |||||
BibTeX:
@inproceedings{Denneulin2020, author = {Laurence Denneulin and Nelly Pustelnik and Maud Langlois and Ignace Loris and Eric Thiébaut}, title = {Primal-dual splitting scheme with backtracking for handling with epigraphic constraint and sparse analysis regularization}, booktitle = {International Traveling Workshop on Interactions between low-complexity data models and Sensing Techniques (iTWIST 2020)}, year = {2020}, url = {https://arxiv.org/abs/2009.03576} } |
|||||
Chen, J. and Loris, I. | On starting and stopping criteria for nested primal-dual iterations | 2019 | Numerical algorithms Vol. 82(2), pp. 605-621 |
article | DOI URL |
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 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. | |||||
BibTeX:
@article{Chen.Loris2019, author = {Jixin Chen and Ignace Loris}, title = {On starting and stopping criteria for nested primal-dual iterations}, journal = {Numerical algorithms}, year = {2019}, volume = {82}, number = {2}, pages = {605--621}, url = {https://arxiv.org/abs/1806.07677}, doi = {https://doi.org/10.1007/s11075-018-0616-x} } |
|||||
Bonettini, S., Loris, I., Porta, F., Prato, M. and Rebegoldi, S. | On the convergence of a linesearch based proximal-gradient method for nonconvex optimization | 2017 | Inverse Problems Vol. 33(5), pp. 055005 |
article | DOI URL |
Abstract: We consider a variable metric line-search 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 minimum point if the objective function satisfies the Kurdyka-Łojasiewicz 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 if compared to recently proposed approaches able to address the optimization problems arising in the considered applications. |
|||||
BibTeX:
@article{Bonettini2016, author = {S. Bonettini and I. Loris and F. Porta and M. Prato and S. Rebegoldi}, title = {On the convergence of a linesearch based proximal-gradient method for nonconvex optimization}, journal = {Inverse Problems}, year = {2017}, volume = {33}, number = {5}, pages = {055005}, url = {https://arxiv.org/abs/1605.03791}, doi = {https://doi.org/10.1088/1361-6420/aa5bfd} } |
|||||
Bonettini, S., Loris, I., Porta, F. and Prato, M. | Variable metric inexact line-search based methods for nonsmooth optimization | 2016 | Siam Journal on Optimization Vol. 26(2), pp. 891-921 |
article | DOI URL |
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 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 𝒪(1k) 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. |
|||||
BibTeX:
@article{Bonettini2015, author = {S. Bonettini and I. Loris and F. Porta and M. Prato}, title = {Variable metric inexact line-search based methods for nonsmooth optimization}, journal = {Siam Journal on Optimization}, year = {2016}, volume = {26}, number = {2}, pages = {891--921}, url = {http://arxiv.org/abs/1506.00385}, doi = {https://doi.org/10.1137/15M1019325} } |
|||||
Prato, M., Bonettini, S., Loris, I., Porta, F. and Rebegoldi, S. | On the constrained minimization of smooth Kurdyka-Łojasiewicz functions with the scaled gradient projection method | 2016 | Journal of Physics: Conference Series Vol. 756(1), pp. 012001 |
article | DOI |
Abstract: The scaled gradient projection (SGP) method is a first-order optimization
method applicable to the constrained minimization of smooth functions and exploiting a scaling matrix multiplying the gradient and a variable steplength parameter to improve the convergence of the scheme. For a general nonconvex function, the limit points of the sequence generated by SGP have been proved to be stationary, while in the convex case and with some restrictions on the choice of the scaling matrix the sequence itself converges to a constrained minimum point. In this paper we extend these convergence results by showing that the SGP sequence converges to a limit point provided that the objective function satisfies the Kurdyka-Łojasiewicz property at each point of its domain and its gradient is Lipschitz continuous. |
|||||
BibTeX:
@article{Prato.ea2016, author = {Marco Prato and Silvia Bonettini and Ignace Loris and Federica Porta and Simone Rebegoldi}, title = {On the constrained minimization of smooth Kurdyka-Łojasiewicz functions with the scaled gradient projection method}, journal = {Journal of Physics: Conference Series}, year = {2016}, volume = {756}, number = {1}, pages = {012001}, note = {Proceedings of 6th International Workshop on New Computational Methods for Inverse Problems (Cachan, 20/05/2016).}, doi = {https://doi.org/10.1088/1742-6596/756/1/012001} } |
|||||
Porta, F. and Loris, I. | On some steplength approaches for proximal algorithms | 2015 | Applied Mathematics and Computation Vol. 253, pp. 345-362 |
article | DOI |
Abstract: We discuss a number of novel steplength selection schemes for proximal-based
convex optimization algorithms. In particular, we consider the problem where the Lipschitz constant of the gradient of the smooth part of the objective function is unknown. We generalize two optimization algorithms of Khobotov type and prove convergence. We also take into account possible inaccurate computation of the proximal operator of the non-smooth part of the objective function. Secondly, we show convergence of an iterative algorithm with Armijo-type steplength rule, and discuss its use with an approximate computation of the proximal operator. Numerical experiments show the efficiency of the methods in comparison to some existing schemes. |
|||||
BibTeX:
@article{Porta2015, author = {Federica Porta and Ignace Loris}, title = {On some steplength approaches for proximal algorithms}, journal = {Applied Mathematics and Computation}, year = {2015}, volume = {253}, pages = {345--362}, doi = {https://doi.org/10.1016/j.amc.2014.12.079} } |
|||||
Loris, I. | Numerical algorithms for non-smooth optimization applicable to seismic recovery | 2014 | Handbook of Geomathematics, pp. 1-33 | incollection | DOI |
Abstract: Inverse problems in seismic tomography are often cast in the form
of an optimization problem involving a cost function composed of a data misfit term and regularizing constraint or penalty. Depending on the noise model that is assumed to underlie the data acquisition, these optimization problems may be non-smooth. Another source of lack of smoothness (differentiability) of the cost function may arise from the regularization method chosen to handle the ill-posed nature of the inverse problem. A numerical algorithm that is well suited to handle minimization problems involving two non-smooth convex functions and two linear operators is studied. The emphasis lies on the use of some simple proximity operators that allow for the iterative solution of non-smooth convex optimization problems. Explicit formulas for several of these proximity operators are given and their application to seismic tomography is demonstrated. |
|||||
BibTeX:
@incollection{Loris2013, author = {Ignace Loris}, title = {Numerical algorithms for non-smooth optimization applicable to seismic recovery}, booktitle = {Handbook of Geomathematics}, publisher = {Springer}, year = {2014}, pages = {1--33}, edition = {Second edition.}, doi = {https://doi.org/10.1007/978-3-642-27793-1_65-3} } |
|||||
Nassiri, V. and Loris, I. | An efficient algorithm for structured sparse quantile regression | 2014 | Computational Statistics Vol. 29(5), pp. 1321-1343 |
article | DOI URL |
Abstract: An efficient algorithm is derived for solving the quantile regression
problem combined with a group sparsity promoting penalty. The group sparsity of the regression parameters is achieved by using a _1,∞-norm penalty (or constraint) on the regression parameters. The algorithm is efficient in the sense that it obtains the regression parameters for a wide range of penalty parameters, thus enabling easy application of a model selection criteria afterwards. A Matlab implementation of the proposed algorithm is provided and some applications of the methods are studied. |
|||||
BibTeX:
@article{Nassiri2014, author = {Vahid Nassiri and Ignace Loris}, title = {An efficient algorithm for structured sparse quantile regression}, journal = {Computational Statistics}, year = {2014}, volume = {29}, number = {5}, pages = {1321-1343}, url = {http://arxiv.org/abs/1302.6088}, doi = {https://doi.org/10.1007/s00180-014-0494-1} } |
|||||
Schretter, C., Loris, I., Dooms, A. and Schelkens, P. | Total Variation reconstruction from quasi-random samples | 2014 | Proceedings of the second ``international Traveling Workshop on Interactions between Sparse models and Technology'' (iTWIST 2014), pp. 57-58 | inproceedings | URL |
Abstract: Pseudo-random numbers are often used for generating incoherent uniformly
distributed sample distributions. However randomness is a sufficient -- not necessary -- condition to ensure incoherence. If one wants to reconstruct an image from few samples, choosing a globally optimized set of evenly distributed points could capture the visual content more efficiently. This work compares classical random sampling with a simple construction based on properties of the fractional Golden ratio sequence and the Hilbert space filling curve. Images are then reconstructed using a total variation prior. Results show improvements in terms of peak signal to noise ratio over pseudo-random sampling. |
|||||
BibTeX:
@inproceedings{Schretter2014, author = {Colas Schretter and Ignace Loris and Ann Dooms and Peter Schelkens}, title = {Total Variation reconstruction from quasi-random samples}, booktitle = {Proceedings of the second ``international Traveling Workshop on Interactions between Sparse models and Technology'' (iTWIST 2014)}, year = {2014}, pages = {57--58}, url = {http://arxiv.org/abs/1410.0719} } |
|||||
Charléty, J., Voronin, S., Nolet, G., Loris, I., Simons, F.J., Sigloch, K. and Daubechies, I.C. | Global seismic tomography with sparsity constraints: Comparison with smoothing and damping regularization | 2013 | Journal of Geophysical Research - Solid Earth Vol. 118, pp. 4887-4899 |
article | DOI |
Abstract: We present a realistic application of an inversion scheme for global
seismic tomography that uses as prior information the sparsity of a solution, defined as having few nonzero coefficients under the action of a linear transformation. In this paper, the sparsifying transform is a wavelet transform. We use an accelerated iterative soft-thresholding algorithm for a regularization strategy, which produces sparse models in the wavelet domain. The approach and scheme we present may be of use for preserving sharp edges in a tomographic reconstruction and minimizing the number of features in the solution warranted by the data. The method is tested on a data set of time delays for finite-frequency tomography using the USArray network, the first application in global seismic tomography to real data. The approach presented should also be suitable for other imaging problems. From a comparison with a more traditional inversion using damping and smoothing constraints, we show that (1) we generally retrieve similar features, (2) fewer nonzero coefficients under a properly chosen representation (such as wavelets) are needed to explain the data at the same level of root-mean-square misfit, (3) the model is sparse or compressible in the wavelet domain, and (4) we do not need to construct a heterogeneous mesh to capture the available resolution. |
|||||
BibTeX:
@article{Charlety2013, author = {Jean Charléty and Sergey Voronin and Guust Nolet and Ignace Loris and Frederik J. Simons and Karin Sigloch and Ingrid C. Daubechies}, title = {Global seismic tomography with sparsity constraints: Comparison with smoothing and damping regularization}, journal = {Journal of Geophysical Research - Solid Earth}, year = {2013}, volume = {118}, pages = {4887--4899}, doi = {https://doi.org/10.1002/jgrb.50326} } |
|||||
Loris, I. and Verhoeven, C. | An iterative algorithm for sparse and constrained recovery with applications to divergence-free current reconstructions in magneto-encephalography | 2013 | Computational Optimization and Applications Vol. 54(2), pp. 399-416 |
article | DOI URL |
Abstract: We propose an iterative algorithm for the minimization of a _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. |
|||||
BibTeX:
@article{Loris2011, author = {Ignace Loris and Caroline Verhoeven}, title = {An iterative algorithm for sparse and constrained recovery with applications to divergence-free current reconstructions in magneto-encephalography}, journal = {Computational Optimization and Applications}, year = {2013}, volume = {54}, number = {2}, pages = {399--416}, note = {Special issue on optimization methods for inverse problems in imaging.}, url = {http://arxiv.org/abs/1202.3362}, doi = {https://doi.org/10.1007/s10589-012-9482-y} } |
|||||
Nassiri, V. and Loris, I. | A generalized quantile regression model | 2013 | Journal of Applied Statistics Vol. 40(5), pp. 1090-1105 |
article | DOI |
Abstract: A new class of probability distributions, the so-called connected
double truncated gamma distribution, is introduced. We show that using this class as the error distribution of a linear model leads to a generalized quantile regression model that combines desirable properties of both least squares and quantile regression methods: robustness to outliers and differentiable loss function. |
|||||
BibTeX:
@article{Nassiri2013, author = {Vahid Nassiri and Ignace Loris}, title = {A generalized quantile regression model}, journal = {Journal of Applied Statistics}, year = {2013}, volume = {40}, number = {5}, pages = {1090--1105}, doi = {https://doi.org/10.1080/02664763.2013.780158} } |
|||||
Cloquet, C., Loris, I., Verhoeven, C. and Defrise, M. | GISTA reconstructs faster with a restart strategy and even faster with a FISTA-like reconstruction | 2012 | 2012 IEEE Nuclear Science Symposium and Medical Imaging Conference, pp. 2334-2338 | inproceedings | DOI |
Abstract: In X-Ray CT, the cone-beam geometry leads to specific artifacts. Moreover,
as the concern about dose-related health effects rises, the struggle to reduce the dose leads to increased image noise. To overcome these image quality issues, a popular recent trend models CT images by flat regions separated by sharp edges, and incorporate this knowledge into a reconstruction algorithm using a total variation (TV) penalty term. Recently, Loris and Verhoeven (2011) designed the Generalization of the Iterative Soft Thresholding Algorithm (GISTA), that has proven convergence and can handle conveniently a non smooth penalty term, like the TV penalty term. To our knowledge, GISTA has not yet been used to reconstruct medical image data. Therefore, we would like to introduce it to the reconstruction community, and present first results of phantom data acquired on a scanner consisting of a cone beam X-ray source and a flat panel detector. We also propose to accelerate the convergence at the initial iterations, using an innovative restart-strategy. |
|||||
BibTeX:
@inproceedings{Cloquet2012, author = {Christophe Cloquet and Ignace Loris and Caroline Verhoeven and Michel Defrise}, title = {GISTA reconstructs faster with a restart strategy and even faster with a FISTA-like reconstruction}, booktitle = {2012 IEEE Nuclear Science Symposium and Medical Imaging Conference}, year = {2012}, pages = {2334--2338}, doi = {https://doi.org/10.1109/NSSMIC.2012.6551530} } |
|||||
Loris, I. | A generalization of the iterative soft-thresholding algorithm for non-separable penalties | 2012 | Vol. 9(2)Oberwolfach Reports, pp. 1811-1814 |
inproceedings | DOI |
Abstract: The present research report focuses on the efficient use of sparse
representations for the regularization and solution of ill-posed problems. In particular, we discuss simple iterative algorithms for the minimization of certain convex functionals encountered in this area. |
|||||
BibTeX:
@inproceedings{Loris2012b, author = {Ignace Loris}, title = {A generalization of the iterative soft-thresholding algorithm for non-separable penalties}, booktitle = {Oberwolfach Reports}, publisher = {European Mathematical Society}, year = {2012}, volume = {9}, number = {2}, pages = {1811--1814}, doi = {https://doi.org/10.4171/OWR/2012/29} } |
|||||
Loris, I. | Wavelets: A Concise Guide. Amir-Homayoon Najmi. 270 pp. The Johns Hopkins University Press, Baltimore, 2012. Price: $45.00 (paper) ISBN 978-1-4214-0496-6 [BibTeX] |
2012 | American Journal of Physics Vol. 80(12), pp. 1113-1113 |
article | DOI |
BibTeX:
@article{Loris2012a, author = {Ignace Loris}, title = {Wavelets: A Concise Guide. Amir-Homayoon Najmi. 270 pp. The Johns Hopkins University Press, Baltimore, 2012. Price: $45.00 (paper) ISBN 978-1-4214-0496-6}, journal = {American Journal of Physics}, year = {2012}, volume = {80}, number = {12}, pages = {1113--1113}, doi = {https://doi.org/10.1119/1.4742757} } |
|||||
Loris, I. and Verhoeven, C. | Iterative algorithms for total variation-like reconstructions in seismic tomography | 2012 | International Journal on Geomathematics Vol. 3(2), pp. 179-208 |
article | DOI URL |
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 of these algorithms is compared numerically in this setting. For the constrained formulation a new algorithm is proposed and its convergence is proven. |
|||||
BibTeX:
@article{Loris2012, author = {Ignace Loris and Caroline Verhoeven}, title = {Iterative algorithms for total variation-like reconstructions in seismic tomography}, journal = {International Journal on Geomathematics}, year = {2012}, volume = {3}, number = {2}, pages = {179--208}, url = {http://arxiv.org/abs/1203.4451}, doi = {https://doi.org/10.1007/s13137-012-0036-3} } |
|||||
Nassiri, V. and Loris, I. | On log-concavity of skew-symmetric distributions and their applications in penalized linear models | 2012 | Proceedings of the 43rd Annual Iranian Mathematics Conference, pp. 1037-1039 | inproceedings | URL |
Abstract: Log-concavity of the skew-symmetric class of distributions is studied.
Also the possibility of using them as error distribution in a sparse linear model is investigated. A procedure to estimate the penalized model is discussed. |
|||||
BibTeX:
@inproceedings{Nassiri2012, author = {Vahid Nassiri and Ignace Loris}, title = {On log-concavity of skew-symmetric distributions and their applications in penalized linear models}, booktitle = {Proceedings of the 43rd Annual Iranian Mathematics Conference}, year = {2012}, pages = {1037--1039}, note = {The 43rd Annual Iranian Mathematics Conference (27--30/8/2012, Tabriz, Iran).}, url = {http://imc43.tabrizu.ac.ir/en/} } |
|||||
Loris, I. and Verhoeven, C. | On a generalization of the iterative soft-thresholding algorithm for the case of non-separable penalty | 2011 | Inverse Problems Vol. 27(12), pp. 125007 |
article | DOI URL |
Abstract: An explicit algorithm for the minimization of an _1 penalized least squares functional, with non-separable _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 _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 _1-norm can be handled as well. | |||||
BibTeX:
@article{Loris.Verhoeven2011, author = {Ignace Loris and Caroline Verhoeven}, title = {On a generalization of the iterative soft-thresholding algorithm for the case of non-separable penalty}, journal = {Inverse Problems}, year = {2011}, volume = {27}, number = {12}, pages = {125007}, url = {http://arxiv.org/abs/1104.1087}, doi = {https://doi.org/10.1088/0266-5611/27/12/125007} } |
|||||
Simons, F.J., Loris, I., Brevdo, E. and Daubechies, I.C. | Wavelets and wavelet-like transforms on the sphere and their application to geophysical data inversion | 2011 | Vol. 8138Wavelets and Sparsity XIV, pp. X1-X15 |
inproceedings | DOI URL |
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 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. |
|||||
BibTeX:
@inproceedings{Simons2011, author = {Frederik J. Simons and Ignace Loris and Eugene Brevdo and Ingrid C. Daubechies}, title = {Wavelets and wavelet-like transforms on the sphere and their application to geophysical data inversion}, booktitle = {Wavelets and Sparsity XIV}, publisher = {SPIE}, year = {2011}, volume = {8138}, pages = {X1--X15}, url = {http://arxiv.org/abs/1109.1718}, doi = {https://doi.org/10.1117/12.892285} } |
|||||
Simons, F.J., Loris, I., Nolet, G., Daubechies, I.C., Voronin, S., Judd, J.S., Vetter, P.A., Charléty, J. and Vonesch, C. | Solving or resolving global tomographic models with spherical wavelets, and the scale and sparsity of seismic heterogeneity | 2011 | Geophysical Journal International Vol. 187(2), pp. 969-988 |
article | DOI URL |
Abstract: We propose a class of spherical wavelet bases for the analysis of
geophysical models and for the 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 _2 norm of data fit and the _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 |
|||||
BibTeX:
@article{Simons.Loris.ea2011, author = {Frederik J. Simons and Ignace Loris and Guust Nolet and Ingrid C. Daubechies and S. Voronin and J. S. Judd and P. A. Vetter and J. Charléty and C. Vonesch}, title = {Solving or resolving global tomographic models with spherical wavelets, and the scale and sparsity of seismic heterogeneity}, journal = {Geophysical Journal International}, year = {2011}, volume = {187}, number = {2}, pages = {969--988}, url = {http://arxiv.org/abs/1104.3151}, doi = {https://doi.org/10.1111/j.1365-246X.2011.05190.x} } |
|||||
Loris, I., Douma, H., Nolet, G., Daubechies, I. and Regone, C. | Nonlinear regularization techniques for seismic tomography | 2010 | Journal of Computational Physics Vol. 229(3), pp. 890-905 |
article | DOI URL |
Abstract: The effects of several nonlinear regularization techniques are discussed
in the framework of 3D seismic tomography. Traditional, linear, _2 penalties are compared to so-called sparsity promoting _1 and _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 _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 _1 damping of Daubechies-4 wavelet coefficients. It has the additional advantage of yielding a noiseless reconstruction, contrary to simple _2 minimization (`Tikhonov regularization') which should be avoided. In some of our examples, the _0 method produced notable artifacts. In addition we show how nonlinear _1 methods for finding sparse models can be competitive in speed with the widely used _2 methods, certainly under noisy conditions, so that there is no need to shun _1 penalizations. |
|||||
BibTeX:
@article{Loris.Douma.ea2009, author = {Loris, I. and Douma, H. and Nolet, G. and Daubechies, I. and Regone, C.}, title = {Nonlinear regularization techniques for seismic tomography}, journal = {Journal of Computational Physics}, year = {2010}, volume = {229}, number = {3}, pages = {890--905}, url = {http://arxiv.org/abs/0808.3472}, doi = {https://doi.org/10.1016/j.jcp.2009.10.020} } |
|||||
Loris, I. and Verhoeven, C. | Practical error estimates for sparse recovery in linear inverse problems. | 2010 | School: Vrije Universiteit Brussel | techreport | URL |
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 are compared to known theoretical bounds and illustrated on a magnetic tomography example. |
|||||
BibTeX:
@techreport{Loris.Verhoeven2010, author = {Ignace Loris and Caroline Verhoeven}, title = {Practical error estimates for sparse recovery in linear inverse problems.}, school = {Vrije Universiteit Brussel}, year = {2010}, url = {http://arxiv.org/abs/0908.3636} } |
|||||
Brodie, J., Daubechies, I., De Mol, C., Giannone, D. and Loris, I. | Sparse and stable Markowitz portfolios | 2009 | Proceedings of the National Academy of Sciences of the USA Vol. 106(30), pp. 12267-12272 |
article | DOI URL |
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. 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. |
|||||
BibTeX:
@article{Brodie.Daubechies.ea2009, author = {Brodie, Joshua and Daubechies, Ingrid and De Mol, Christine and Giannone, Domenico and Loris, Ignace}, title = {Sparse and stable Markowitz portfolios}, journal = {Proceedings of the National Academy of Sciences of the USA}, year = {2009}, volume = {106}, number = {30}, pages = {12267--12272}, url = {http://arxiv.org/abs/0708.0046}, doi = {https://doi.org/10.1073/pnas.0904287106} } |
|||||
Loris, I. | On the performance of algorithms for the minimization of _1-penalized functionals | 2009 | Inverse Problems Vol. 25(3), pp. 035008 (16pp) |
article | DOI URL |
Abstract: The problem of assessing the performance of algorithms used for the
minimization of an _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. |
|||||
BibTeX:
@article{Loris2009, author = {Loris, Ignace}, title = {On the performance of algorithms for the minimization of _1-penalized functionals}, journal = {Inverse Problems}, year = {2009}, volume = {25}, number = {3}, pages = {035008 (16pp)}, url = {http://arxiv.org/abs/0710.4082}, doi = {https://doi.org/10.1088/0266-5611/25/3/035008} } |
|||||
Loris, I., Bertero, M., De Mol, C., Zanella, R. and Zanni, L. | Accelerating gradient projection methods for _1- constrained signal recovery by steplength selection rules | 2009 | Applied and Computational Harmonic Analysis Vol. 27(2), pp. 247-254 |
article | DOI URL |
Abstract: We propose a new gradient projection algorithm that compares favorably
with the fastest algorithms available to date for _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. |
|||||
BibTeX:
@article{Loris.Bertero.ea2009, author = {Loris, I. and Bertero, M. and De Mol, C. and Zanella, R. and Zanni, L.}, title = {Accelerating gradient projection methods for _1- constrained signal recovery by steplength selection rules}, journal = {Applied and Computational Harmonic Analysis}, year = {2009}, volume = {27}, number = {2}, pages = {247--254}, url = {http://arxiv.org/abs/0902.4424}, doi = {https://doi.org/10.1016/j.acha.2009.02.003} } |
|||||
Daubechies, I., Fornasier, M. and Loris, I. | Accelerated projected gradient method for linear inverse problems with sparsity constraints | 2008 | Journal of Fourier Analysis and Applications Vol. 14(5-6), pp. 764-792 |
article | DOI URL |
Abstract: Regularization of ill-posed linear inverse problems via _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 _1 penalized functional is via an iterative soft-thresholding algorithm. We propose an alternative implementation to _1-constraints, using a gradient method, with projection on _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. |
|||||
BibTeX:
@article{Daubechies.Fornasier.ea2008, author = {Ingrid Daubechies and Massimo Fornasier and Ignace Loris}, title = {Accelerated projected gradient method for linear inverse problems with sparsity constraints}, journal = {Journal of Fourier Analysis and Applications}, year = {2008}, volume = {14}, number = {5--6}, pages = {764--792}, url = {http://arxiv.org/abs/0706.4297}, doi = {https://doi.org/10.1007/s00041-008-9039-8} } |
|||||
Loris, I. | L1Packv2: A Mathematica package for minimizing an _1-penalized functional | 2008 | Computer Physics Communications Vol. 179(12), pp. 895-902 |
article | DOI URL |
Abstract: L1Packv2 is a Mathematica package that contains a number of algorithms
that can be used for the minimization of an _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. |
|||||
BibTeX:
@article{Loris2008, author = {Loris, Ignace}, title = {L1Packv2: A Mathematica package for minimizing an _1-penalized functional}, journal = {Computer Physics Communications}, year = {2008}, volume = {179}, number = {12}, pages = {895--902}, url = {http://arxiv.org/abs/0710.3728}, doi = {https://doi.org/10.1016/j.cpc.2008.07.010} } |
|||||
Loris, I., Nolet, G., Daubechies, I. and Dahlen, F.A. | Tomographic inversion using _1-norm regularization of wavelet coefficients | 2007 | Geophysical Journal International Vol. 170(1), pp. 359-370 |
article | DOI URL |
Abstract: We propose the use of _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. |
|||||
BibTeX:
@article{Loris.Nolet.ea2007, author = {Loris, Ignace and Nolet,Guust and Daubechies, Ingrid and Dahlen, F. A.}, title = {Tomographic inversion using _1-norm regularization of wavelet coefficients}, journal = {Geophysical Journal International}, year = {2007}, volume = {170}, number = {1}, pages = {359--370}, url = {http://arxiv.org/abs/physics/0608094}, doi = {https://doi.org/10.1111/j.1365-246X.2007.03409.x} } |
|||||
Khare, A., Loris, I. and Sasaki, R. | Affine Toda-Sutherland systems | 2004 | Journal of Physics A-Mathematical and General Vol. 37(5), pp. 1665-1679 |
article | DOI URL |
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 frequencies of small oscillations near equilibrium are proportional to the affine Toda masses, which are essential ingredients of the exact factorizable 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. |
|||||
BibTeX:
@article{Khare.Loris.ea2004, author = {Khare, A. and Loris, I. and Sasaki, R.}, title = {Affine Toda-Sutherland systems}, journal = {Journal of Physics A-Mathematical and General}, year = {2004}, volume = {37}, number = {5}, pages = {1665--1679}, url = {http://arxiv.org/abs/hep-th/0309077}, doi = {https://doi.org/10.1088/0305-4470/37/5/013} } |
|||||
Loris, I. and Sasaki, R. | Quantum and classical eigenfunctions in Calogero and Sutherland systems | 2004 | Journal of Physics A-Mathematical and General Vol. 37(1), pp. 211-237 |
article | DOI URL |
Abstract: An interesting observation was reported by Corrigan-Sasaki that all
the frequencies of small oscillations around equilibrium are 'quantized' for Calogero and Sutherland (CS) 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 CS systems based on any root system |
|||||
BibTeX:
@article{Loris.Sasaki2004, author = {Loris, I. and Sasaki, R.}, title = {Quantum and classical eigenfunctions in Calogero and Sutherland systems}, journal = {Journal of Physics A-Mathematical and General}, year = {2004}, volume = {37}, number = {1}, pages = {211--237}, url = {http://arxiv.org/abs/hep-th/0308052}, doi = {https://doi.org/10.1088/0305-4470/37/1/015} } |
|||||
Loris, I. and Sasaki, R. | Quantum vs classical mechanics: role of elementary excitations | 2004 | Physics Letters A Vol. 327(2-3), pp. 152-157 |
article | DOI URL |
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. |
|||||
BibTeX:
@article{Loris.Sasaki2004a, author = {Loris, I. and Sasaki, R.}, title = {Quantum vs classical mechanics: role of elementary excitations}, journal = {Physics Letters A}, year = {2004}, volume = {327}, number = {2-3}, pages = {152--157}, url = {http://arxiv.org/abs/quant-ph/0308040}, doi = {https://doi.org/10.1016/j.physleta.2004.05.015} } |
|||||
Loris, I. | Bilinear representations of integrable equations | 2002 | Theoretical and Mathematical Physics Vol. 133(2), pp. 1549-1556 |
article | DOI |
Abstract: We present a method for deriving recursion operators and canonical
Lax pairs directly from bilinear identities of the KP type. Examples include the KdV equation, the Boussinesq equation, and a real equivalent of the nonlinear Schrödinger equation. |
|||||
BibTeX:
@article{Loris2002, author = {Loris, I.}, title = {Bilinear representations of integrable equations}, journal = {Theoretical and Mathematical Physics}, year = {2002}, volume = {133}, number = {2}, pages = {1549--1556}, note = {Proceedings of the NEEDS'01 Conference.}, doi = {https://doi.org/10.1023/A:1021103012057} } |
|||||
Lambert, F., Loris, I. and Springael, J. | Classical Darboux transformations and the KP hierarchy | 2001 | Inverse Problems Vol. 17(4), pp. 1067-1074 |
article | DOI |
Abstract: Classical Darboux transformations together with partitional polynomials
are used as elementary tools for the construction of members of a basic hierarchy of integrable nonlinear partial differential equations (the KP hierarchy). |
|||||
BibTeX:
@article{Lambert.Loris.ea2001, author = {Lambert, F. and Loris, I. and Springael, J.}, title = {Classical Darboux transformations and the KP hierarchy}, journal = {Inverse Problems}, year = {2001}, volume = {17}, number = {4}, pages = {1067--1074}, doi = {https://doi.org/10.1088/0266-5611/17/4/333} } |
|||||
Lambert, F., Loris, I., Springael, J. and Willox, R. | On the Hirota representation of soliton equations with one tau-function | 2001 | Journal of the Physical Society of Japan Vol. 70(3), pp. 605-608 |
article | DOI |
Abstract: Alternative Hirota representations in terms of a single tau-function
are derived for a variaty of solition equations, including the sine-Gordon and Tzitzeica equations. The relevance of these representations with respect to known bilinear representations of integrable hierarchies is briefly discussed. The essentials of the derivations method are presented. |
|||||
BibTeX:
@article{Lambert.Loris.ea2001a, author = {Lambert, F. and Loris, I. and Springael, J. and Willox, R.}, title = {On the Hirota representation of soliton equations with one tau-function}, journal = {Journal of the Physical Society of Japan}, year = {2001}, volume = {70}, number = {3}, pages = {605--608}, doi = {https://doi.org/10.1143/JPSJ.70.605} } |
|||||
Loris, I. | Dimensional reductions of BKP and CKP hierarchies | 2001 | Journal of Physics A-Mathematical and General Vol. 34(16), pp. 3447-3459 |
article | DOI |
Abstract: A discussion of dimensional reductions, which are not classical symmetry
reductions, is made for the BKP and CKP hierarchies of integrable evolution equations. A novel direct method for testing Pfaffian solutions to bilinearidentities is presented and applied to these reductions. |
|||||
BibTeX:
@article{Loris2001, author = {Loris, I.}, title = {Dimensional reductions of BKP and CKP hierarchies}, journal = {Journal of Physics A-Mathematical and General}, year = {2001}, volume = {34}, number = {16}, pages = {3447--3459}, doi = {https://doi.org/10.1088/0305-4470/34/16/313} } |
|||||
Loris, I. | Solutions of coupled Korteweg-de Vries systems | 2001 | Journal of the Physical Society of Japan Vol. 70(3), pp. 662-665 |
article | DOI |
Abstract: A class of solutions in determinant form to a set of coupled KdV equations
is derived. The bilinear (Hirota) form of the coupled equations and links with the CKP system serve as a guide. |
|||||
BibTeX:
@article{Loris2001a, author = {Loris, I.}, title = {Solutions of coupled Korteweg-de Vries systems}, journal = {Journal of the Physical Society of Japan}, year = {2001}, volume = {70}, number = {3}, pages = {662--665}, doi = {https://doi.org/10.1143/JPSJ.70.662} } |
|||||
Loris, I. | Recursion operator for a constrained BKP system | 2000 | Proceedings of the workshop on Nonlinearity, integrability, and all that: 20 years after NEEDS '79, pp. 325-330 | inproceedings | |
Abstract: Symmetry reductions of the KP system of B-type are discussed. A recursion
operator (mapping symmetries into symmetries) is derived for the main example. |
|||||
BibTeX:
@inproceedings{Loris2000, author = {Loris, Ignace}, title = {Recursion operator for a constrained BKP system}, booktitle = {Proceedings of the workshop on Nonlinearity, integrability, and all that: 20 years after NEEDS '79}, publisher = {World Scientific, Singapore}, year = {2000}, pages = {325--330} } |
|||||
Loris, I. | On reduced CKP equations | 1999 | Inverse Problems Vol. 15(4), pp. 1099-1109 |
article | DOI |
Abstract: Symmetry reductions of the CKP hierarchy are discussed in a pseudodifferential
and tau-function context. Solutions of the resulting nonlinear partial differential equations are obtained via the methods of gauge transformations and of tau functions. Reductions of other 2+1-dimensional hierarchies related to KP are also briefly investigated. |
|||||
BibTeX:
@article{Loris1999, author = {Loris, I.}, title = {On reduced CKP equations}, journal = {Inverse Problems}, year = {1999}, volume = {15}, number = {4}, pages = {1099--1109}, doi = {https://doi.org/10.1088/0266-5611/15/4/317} } |
|||||
Loris, I. and Willox, R. | Symmetry reductions of the BKP hierarchy | 1999 | Journal of Mathematical Physics Vol. 40(3), pp. 1420-1431 |
article | DOI |
Abstract: A general symmetry of the bilinear BKP hierarchy is studied in terms
of tau functions. We use this symmetry to define reductions of the BKP hierarchy, among which new integrable systems can be found. The reductions are connected to constraints on the Lax operator as well as on the bilinear formulation. A class of solutions for the reduced equations is derived |
|||||
BibTeX:
@article{Loris.Willox1999, author = {Loris, I. and Willox, R.}, title = {Symmetry reductions of the BKP hierarchy}, journal = {Journal of Mathematical Physics}, year = {1999}, volume = {40}, number = {3}, pages = {1420--1431}, doi = {https://doi.org/10.1063/1.532812} } |
|||||
Willox, R. and Loris, I. | An algebraic description of generalized k-constraints | 1999 | Journal of Physics A-Mathematical and General Vol. 32(10), pp. 2027-2036 |
article | DOI |
Abstract: The generalized k-constrained KP hierarchy is shown to correspond
to a so-called pseudo-reduction of the two-dimensional Toda lattice hierarchy, described in a free-fermion approach which is adapted to the case of two singularities in the spectral parameter range. Wronskian solutions are discussed and, in particular, soliton solutions are recovered through a p^k + c/p = q^k + c/q reduction of the Toda solitons. |
|||||
BibTeX:
@article{Willox.Loris1999a, author = {Willox, R. and Loris, I.}, title = {An algebraic description of generalized k-constraints}, journal = {Journal of Physics A-Mathematical and General}, year = {1999}, volume = {32}, number = {10}, pages = {2027--2036}, doi = {https://doi.org/10.1088/0305-4470/32/10/018} } |
|||||
Willox, R. and Loris, I. | KP constraints from reduced multi-component hierarchies | 1999 | Journal of Mathematical Physics Vol. 40(12), pp. 6501-6525 |
article | DOI |
Abstract: The m-vector k-constrained Kadomtsev�Petviashvili KP hierarchy
is shown to be a ''pseudo''-reduction of the (m+1)-component KP hierarchy. To facilitate the implementation of this reduction on the level of the solutions, the typical multicomponent KP solutions are mapped onto solutions of a Toda molecule-type equation from which Wronskian and Grammian solutions for the constrained KP hierarchy follow. The reduction of the associated linear systems is discussed and its importance for the choice of bilinear representation of the reduced systems is explained. |
|||||
BibTeX:
@article{Willox.Loris1999, author = {Willox, R. and Loris, I.}, title = {KP constraints from reduced multi-component hierarchies}, journal = {Journal of Mathematical Physics}, year = {1999}, volume = {40}, number = {12}, pages = {6501--6525}, doi = {https://doi.org/10.1063/1.533104} } |
|||||
Loris, I. | Symmetry reductions in the tau-function approach to integrability | 1998 | School: Vrije Universiteit Brussel | phdthesis | |
Abstract: It is the aim of this thesis to investigate how different 1+1 dimensional
``soliton'' systems and the various properties characterizing their ``integrability'', may be derived from a single universal equation through reduction procedures. These reductions are implemented by imposing judiciously chosen constraints on the field in which these equations take their most fundamental form: the ``tau-function''. Soliton systems, such as the non-linear Schrödinger equation or the Korteweg-de Vries equation, feature among the set of integrable non-linear partial differential equations. The term `soliton' refers to the existence of particular solutions which describe the interaction of a (unlimited) number of localized pulses. The soliton phenomenon was intensively studied in the years following its identification. However, over the past decade, attention and research efforts have shifted towards the concept of ``integrability'' in systems with an infinite number of degrees of freedom. For these systems, various criteria and tests of integrability exist. Different criteria require the existence of an infinite number of conservation laws (or an infinite number of symmetries), or the absence of certain movable singularities in the general solution (Painlevé property), or the existence of an underlying linear Lax system and a corresponding inverse scattering problem (IST) … It is a generally accepted view that, despite the existence of partial results, the problem of understanding the notion of integrability is a difficult one. It is not our purpose to settle this question here. Yet, just as the different integrability criteria are inter-related, so do a number of different 1+1 dimensional soliton systems exhibit astoundingly similar features. It is our intention to show that these features spring from the characteristics of an underlying 2+1 dimensional integrable hierarchy: the Kadomtsev-Petviashvili (KP) hierarchy. The important point to be understood is how 1+1 dimensional soliton systems can be regarded as reductions of a higher dimensional integrable equation. More precisely, using the symmetries of this KP hierarchy, we shall present a new and unified approach to reductions of the KP hierarchy to 1+1 dimensional integrable systems. The appropriate technique for the implementation of theses reductions is the method of symmetry reductions. As the set of symmetries (of a given equation) exhibits a linear structure (i.e. the defining equations for symmetries are linear partial differential equations) one can impose the condition that a certain symmetry be identically equal to zero. This condition will not be in contradiction with the integrable structure of the hierarchy, and hence represents a sound reduction technique. In order to apply these reductions successfully, one must express them in the appropriate field variable. The fundamental field-variable in the KP hierarchy is the `tau-function'. This tau-function ((t_1,t_2,)) is the solution of Hirota-type partial differential equations, i.e. equation which are quadratically homogeneous in the field variable. The first topic in our investigation is the determination of a general symmetry for the bilinear equations. This symmetry can be written in terms of an eigenfunction potential (eigenfunctions are the solutions of the underlying linear problem). An important result is the fact that this potential can in its turn also be expressed entirely in terms of tau-functions. The `k-constrained' KP hierarchy is defined by imposing a symmetry constraint between the elementary symmetry _t_kτ and a eigenfunction potential symmetry. We show that this reduction can be expressed entirely in terms of tau-functions. A characterization of the corresponding Lax operator is given. Additional bilinear identities for the KP tau-function, which give a bilinear description of the constraint, are derived. Our description of the reduced hierarchies gives access to the determination of general classes of determinant-type solutions. The proposed technique can also be used to determine solutions to the unreduced KP hierarchy and sets itself apart from known determinant techniques by the fact that solutions are shown to exist for all the equations in the hierarchy at once. Finally, a generalized k-constraint is introduced which is especially adapted to the description of systems with non-zero boundary conditions. This extension is physically relevant in view of the fact that the non-linear Schrödinger equation (modeling light waves in non-linear media) also exhibits solutions with non-vanishing boundary conditions. A bilinear description of this reduction is given and soliton solutions for the reduced systems are derived. Furthermore a Bäcklund transformation and a Lax pair are derived as well. |
|||||
BibTeX:
@phdthesis{Loris1998, author = {Loris, Ignace}, title = {Symmetry reductions in the tau-function approach to integrability}, school = {Vrije Universiteit Brussel}, year = {1998} } |
|||||
Willox, R., Tokihiro, T., Loris, I. and Satsuma, J. | The fermionic approach to Darboux transformations | 1998 | Inverse Problems Vol. 14(3), pp. 745-762 |
article | DOI |
Abstract: Starting from the free fermion description of the one-component KP
hierarchy, we establish a connection between this approach and the theory of Darboux and binary Darboux transformations. Certain difference identities---allowing for the treatment of both continuous as well as discrete evolution equations---turn out to be crucial: first to show that any solution of the associated (adjoint) linear problems can always be expressed as a superposition of KP (adjoint) wavefunctions and then to interpret Darboux (and binary Darboux) transformations as Bäcklund transformations in the fermion language. |
|||||
BibTeX:
@article{Willox.Tokihiro.ea1998, author = {Willox, R. and Tokihiro, T. and Loris, I. and Satsuma, J.}, title = {The fermionic approach to Darboux transformations}, journal = {Inverse Problems}, year = {1998}, volume = {14}, number = {3}, pages = {745--762}, doi = {https://doi.org/10.1088/0266-5611/14/3/022} } |
|||||
Loris, I. and Willox, R. | Bilinear form and solutions of the k-constrained Kadomtsev-Petviashvili hierarchy | 1997 | Inverse Problems Vol. 13(2), pp. 411-420 |
article | DOI |
Abstract: We show how to derive an alternative bilinear formulation for the
k-constrained Kadomtsev�Petviashvili hierarchy. This Hirota form allows for the easy identification of a broad class of solutions to these equations |
|||||
BibTeX:
@article{Loris.Willox1997c, author = {Loris, I. and Willox, R.}, title = {Bilinear form and solutions of the k-constrained Kadomtsev-Petviashvili hierarchy}, journal = {Inverse Problems}, year = {1997}, volume = {13}, number = {2}, pages = {411--420}, doi = {https://doi.org/10.1088/0266-5611/13/2/014} } |
|||||
Loris, I. and Willox, R. | KP symmetry reductions and a generalized constraint | 1997 | Journal of Physics A-Mathematical and General Vol. 30(19), pp. 6925-6938 |
article | DOI |
Abstract: We study the link between symmetry reductions and constraints of the
Kadomtsev-Petviashvili equations in terms of the tau function. We propose a generalization---adapted to non-zero boundary conditions---of the standard constraints, and show a particular class of solutions (solitons). |
|||||
BibTeX:
@article{Loris.Willox1997a, author = {Loris, I. and Willox, R.}, title = {KP symmetry reductions and a generalized constraint}, journal = {Journal of Physics A-Mathematical and General}, year = {1997}, volume = {30}, number = {19}, pages = {6925--6938}, doi = {https://doi.org/10.1088/0305-4470/30/19/027} } |
|||||
Loris, I. and Willox, R. | On solutions of constrained Kadomtsev-Petviashvili equations: Grammians | 1997 | Journal of Mathematical Physics Vol. 38(10), pp. 5190-5197 |
article | DOI |
Abstract: We show the existence of Grammian-type solutions for the (vector)
k-constrained Kadomtsev-Petviashvili (KP) equations. To introduce the method we give a novel proof for the presence of Grammian solutions for the bilinear l-modified KP hierarchies |
|||||
BibTeX:
@article{Loris.Willox1997, author = {Loris, I. and Willox, R.}, title = {On solutions of constrained Kadomtsev-Petviashvili equations: Grammians}, journal = {Journal of Mathematical Physics}, year = {1997}, volume = {38}, number = {10}, pages = {5190--5197}, doi = {https://doi.org/10.1063/1.531937} } |
|||||
Loris, I. and Willox, R. | On solutions of constrained KP equations | 1997 | Journal of Mathematical Physics Vol. 38(1), pp. 283-291 |
article | DOI |
Abstract: We derive solutions of general Wronskian form for the (vector) constrained
KP hierarchies. As one explicit example we discuss rational solutions. In order to introduce our method, we give a direct, elementary proof of the existence of Wronskian solutions for the l-modified KP hierarchies (l:0,1,…). |
|||||
BibTeX:
@article{Loris.Willox1997b, author = {Loris, I. and Willox, R.}, title = {On solutions of constrained KP equations}, journal = {Journal of Mathematical Physics}, year = {1997}, volume = {38}, number = {1}, pages = {283--291}, doi = {https://doi.org/10.1063/1.531843} } |
|||||
Loris, I. and Willox, R. | Recent results on dimensional reductions of the Kadomtsev-Petviashvili equation | 1997 | Proceedings of the Fourth National Conference on Theoretical and Applied Mechanics, pp. 27-30 | inproceedings | |
Abstract: In the present paper, we study a dimensional reductions of the Kadomtsev-Petviashvili
(KP) equations, amonst which one finds the non-linear Schrödinger equation (with non-zero boundary conditions). In particular we use the Darboux transformation to obtain the N soliton solution (∀ N). |
|||||
BibTeX:
@inproceedings{Loris.Willox1997d, author = {Loris, Ignace and Willox, Ralph}, title = {Recent results on dimensional reductions of the Kadomtsev-Petviashvili equation}, booktitle = {Proceedings of the Fourth National Conference on Theoretical and Applied Mechanics}, year = {1997}, pages = {27--30} } |
|||||
Pelinovsky, D., Springael, J., Lambert, F. and Loris, I. | On modified NLS, Kaup and NLBq equations: differential transformations and bilinearization | 1997 | Journal of Physics A-Mathematical and General Vol. 30(24), pp. 8705-8717 |
article | DOI |
Abstract: New transformations between the nonlinear Schrödinger, Kaup and
non-local Boussinesq equations as well as their modified counterparts are found and analysed. The bilinear representations of these equations, including an alternative bilinear form of the Chen-Lee-Liu equation, are obtained by a direct method based on the Bell's exponential polynomials. Explicit Wronskian solutions to these equations are also presented. |
|||||
BibTeX:
@article{Pelinovsky.Springael.ea1997, author = {Pelinovsky, D. and Springael, J. and Lambert, F. and Loris, I.}, title = {On modified NLS, Kaup and NLBq equations: differential transformations and bilinearization}, journal = {Journal of Physics A-Mathematical and General}, year = {1997}, volume = {30}, number = {24}, pages = {8705--8717}, doi = {https://doi.org/10.1088/0305-4470/30/24/029} } |
|||||
Willox, R., Loris, I. and Gilson, C.R. | Binary Darboux transformations for constrained KP hierarchies | 1997 | Inverse Problems Vol. 13(3), pp. 849-865 |
article | DOI |
Abstract: We describe how Darboux transformations and binary Darboux transformations
can be constructed for (vector-) constrained KP hierarchies. These transformations are then used to obtain explicit classes of Wronskian and Grammian solutions for these hierarchies. The relationship between these two types of solutions is also discussed. |
|||||
BibTeX:
@article{Willox.Loris.ea1997, author = {Willox, R. and Loris, I. and Gilson, C. R.}, title = {Binary Darboux transformations for constrained KP hierarchies}, journal = {Inverse Problems}, year = {1997}, volume = {13}, number = {3}, pages = {849--865}, doi = {https://doi.org/10.1088/0266-5611/13/3/019} } |
|||||
Loris, I., Lambert, F. and Willox, R. | New ways of applying the Hirota method in soliton theory | 1996 | Journal of Technical Physics Vol. 37, pp. 519-522 |
article | |
Abstract: It is shown how a sech-squared soliton hierarchy may be obtained from
the classical Boussinesq hierarchy. A collective approach to the problem of its bilinearization is discussed. Solutions and relations to other integrable systems are investigated from a bilinear point of view. |
|||||
BibTeX:
@article{Loris.Lambert.ea1996, author = {Loris, I. and Lambert, F. and Willox, R.}, title = {New ways of applying the Hirota method in soliton theory}, journal = {Journal of Technical Physics}, year = {1996}, volume = {37}, pages = {519--522}, note = {Proceedings of the Conference on Nonlinear Dynamics, Chaotic and Complex Systems (NDCCS'95), Zakopane, Poland.} } |
|||||
Loris, I. and Willox, R. | Soliton solutions of Wronskian form to the nonlocal Boussinesq equation | 1996 | Journal of the Physical Society of Japan Vol. 65(2), pp. 383-388 |
article | DOI |
Abstract: We investigate how one may write the soliton solutions of a nonlocal
Boussinesq equation in Wronskian form, and subsequently prove the existence of N-soliton solutions making use of the bilinear form of this equation. This technique also allows us to construct a bilinear Bäcklund transformation for this equation, mapping N-soliton solutions on (N+1)-soliton solutions. Our analysis extends the results previously obtained by Hirota for the Classical Boussinesq system to actual (c≠ 0) ''pq=c''-reductions performed on Wronskians. |
|||||
BibTeX:
@article{Loris.Willox1996, author = {Loris, I. and Willox, R.}, title = {Soliton solutions of Wronskian form to the nonlocal Boussinesq equation}, journal = {Journal of the Physical Society of Japan}, year = {1996}, volume = {65}, number = {2}, pages = {383--388}, doi = {https://doi.org/10.1143/JPSJ.65.383} } |
|||||
Springael, J., Hu, X.B. and Loris, I. | Bilinear characterization of higher order Ito-equations | 1996 | Journal of the Physical Society of Japan Vol. 65(5), pp. 1222-1226 |
article | DOI |
Abstract: We construct an infinite family of soliton equations, the lowest member
of which corresponds to the Ito-equation: using the bilinear formalism, we obtain a generic bilinear form leading to a recursion operator for this family. A generic Lax-pair and bilinear Bäcklund-transformation are also reported. |
|||||
BibTeX:
@article{Springael.Hu.ea1996, author = {Springael, J. and Hu, X. B. and Loris, I.}, title = {Bilinear characterization of higher order Ito-equations}, journal = {Journal of the Physical Society of Japan}, year = {1996}, volume = {65}, number = {5}, pages = {1222--1226}, doi = {https://doi.org/10.1143/JPSJ.65.1222} } |
|||||
Willox, R. and Loris, I. | Symmetry constraints of the KP hierarchy and a nonlocal Boussinesq equation | 1996 | Vol. 2VII International Conference Symmetry Methods in Physics, pp. 603-609 |
inproceedings | |
Abstract: In this paper we use the Lax pair for the nonlocal Boussinesq equation
in order to show that this equation can be interpreted as a 1-constraint (or corresponding symmetry constraint) on the KP linear problem. A link with density constraints on the modified-KP hierarchy, leading to the classical Boussinesq eaution is also disclosed. |
|||||
BibTeX:
@inproceedings{Willox.Loris1996, author = {Willox, R. and Loris, I.}, title = {Symmetry constraints of the KP hierarchy and a nonlocal Boussinesq equation}, booktitle = {VII International Conference Symmetry Methods in Physics}, publisher = {Joint Institute for Nuclear Research}, year = {1996}, volume = {2}, pages = {603--609} } |
|||||
Willox, R., Loris, I. and Springael, J. | The nlBq-hierarchy as a pq=C reduction of the KP-hierarchy | 1996 | Proceedings of the Workshop``Non-linear Physics, Theory and Experiment", pp. 321-329 | inproceedings | |
Abstract: In this paper it is shown that the so-called nonlocal Boussinesq equation
can be regarded as pq=c reduction of the KP hierarchy. The reduction procedure can be carried out explicitly on the N-soliton solutions of the KP hierarchy once a suitable bilinear form of the nonlocal Boussinesq equation is obtained. It is however necessary to impose an extra constraint on the KP evolution in order to implement the reduction on the bilinear forms themselves. This special constraint allows for yet another class of solutions. |
|||||
BibTeX:
@inproceedings{Willox.Loris.ea1996, author = {Willox, R. and Loris, I. and Springael, J.}, title = {The nlBq-hierarchy as a pq=C reduction of the KP-hierarchy}, booktitle = {Proceedings of the Workshop``Non-linear Physics, Theory and Experiment"}, publisher = {World Scientific, Singapore}, year = {1996}, pages = {321--329} } |
|||||
Lambert, F., Loris, I., Springael, J. and Willox, R. | A direct bilinearization scheme based on the use of partition polynomials | 1995 | Proceedings of the NEEDS'94 workshop at Los Alamos NL, pp. 102-111 | inproceedings | |
Abstract: A systematic procedure for the bilinerization of classes of soliton
equations (and related ordinary differential equations) is presented. The method, based on the use of generalized Bell polynomials, is illustrated with several examples. |
|||||
BibTeX:
@inproceedings{Lambert.Loris.ea1995, author = {Lambert, F. and Loris, I. and Springael, J. and Willox, R.}, title = {A direct bilinearization scheme based on the use of partition polynomials}, booktitle = {Proceedings of the NEEDS'94 workshop at Los Alamos NL}, publisher = {World Scientific, Singapore}, year = {1995}, pages = {102--111} } |
|||||
Willox, R., Loris, I. and Springael, J. | Bilinearization of the nonlocal Boussinesq equation | 1995 | Journal of Physics A-Mathematical and General Vol. 28(20), pp. 5963-5972 |
article | DOI |
Abstract: A single-field bilinear system generating the so-called non-local
Boussinesq equation is constructed. From the bilinearization procedure it can be seen that the associated hierarchy of soliton systems which we construct shares part of the solution set of a hierarchy related to the Kadomtsev-Petviashvili equation. A bilinear representation of the recursion operator for the Kaup hierarchy is essential in the construction and a systematic way of obtaining such a representation from just two-soliton considerations is presented. |
|||||
BibTeX:
@article{Willox.Loris.ea1995, author = {Willox, R. and Loris, I. and Springael, J.}, title = {Bilinearization of the nonlocal Boussinesq equation}, journal = {Journal of Physics A-Mathematical and General}, year = {1995}, volume = {28}, number = {20}, pages = {5963--5972}, doi = {https://doi.org/10.1088/0305-4470/28/20/024} } |
|||||
Lambert, F., Loris, I., Springael, J. and Willox, R. | On a direct bilinearization method - Kaup's higher-order water-wave equation as a modified nonlocal Boussinesq equation | 1994 | Journal of Physics A-Mathematical and General Vol. 27(15), pp. 5325-5334 |
article | DOI |
Abstract: A systematic procedure for the bilinearization of classes of soliton
equations is developed with the help of a generalization of Bell's exponential polynomials. Application of this procedure to Kaup's higher-order wave equation discloses several links with other soliton systems. In particular, it is found that the Kaup equation is the modified version of a sech square soliton system which constitutes an alternative to the good Boussinesq equation. |
|||||
BibTeX:
@article{Lambert.Loris.ea1994, author = {Lambert, F. and Loris, I. and Springael, J. and Willox, R.}, title = {On a direct bilinearization method - Kaup's higher-order water-wave equation as a modified nonlocal Boussinesq equation}, journal = {Journal of Physics A-Mathematical and General}, year = {1994}, volume = {27}, number = {15}, pages = {5325--5334}, doi = {https://doi.org/10.1088/0305-4470/27/15/028} } |
|||||
Loris, I. | Rechtstreekse bilinearisatie van solitonsystemen: Een systematische aanpak met veralgemeende Bell-polynomen [BibTeX] |
1994 | School: Vrije Universiteit Brussel | mastersthesis | |
BibTeX:
@mastersthesis{Loris1994, author = {Loris, Ignace}, title = {Rechtstreekse bilinearisatie van solitonsystemen: Een systematische aanpak met veralgemeende Bell-polynomen}, school = {Vrije Universiteit Brussel}, year = {1994} } |