scholarly journals Stability properties of a projector-splitting scheme for dynamical low rank approximation of random parabolic equations

2021 ◽  
Author(s):  
Yoshihito Kazashi ◽  
Fabio Nobile ◽  
Eva Vidličková
Keyword(s):  
Parabolic Equations ◽  
Parabolic Type ◽  
Implicit Scheme ◽  
Low Rank ◽  
Splitting Scheme ◽  
Low Rank Approximation ◽  
Unconditionally Stable ◽  
Fully Discrete ◽  
Rank Approximation ◽  
Smallest Singular Value

AbstractWe consider the Dynamical Low Rank (DLR) approximation of random parabolic equations and propose a class of fully discrete numerical schemes. Similarly to the continuous DLR approximation, our schemes are shown to satisfy a discrete variational formulation. By exploiting this property, we establish stability of our schemes: we show that our explicit and semi-implicit versions are conditionally stable under a “parabolic” type CFL condition which does not depend on the smallest singular value of the DLR solution; whereas our implicit scheme is unconditionally stable. Moreover, we show that, in certain cases, the semi-implicit scheme can be unconditionally stable if the randomness in the system is sufficiently small. Furthermore, we show that these schemes can be interpreted as projector-splitting integrators and are strongly related to the scheme proposed in [29, 30], to which our stability analysis applies as well. The analysis is supported by numerical results showing the sharpness of the obtained stability conditions.

scholarly journals Low-rank approximation of linear parabolic equations by space-time tensor Galerkin methods

2019 ◽  
Vol 53 (2) ◽  
pp. 635-658
Author(s):  
Thomas Boiveau ◽  
Virginie Ehrlacher ◽  
Alexandre Ern ◽  
Anthony Nouy
Keyword(s):  
Evolution Equations ◽  
Differential Operators ◽  
Space Time ◽  
Low Rank ◽  
Low Rank Approximation ◽  
Dual Norm ◽  
Fully Discrete ◽  
Linear Parabolic Equations ◽  
Parabolic Evolution Equations ◽  
Rank Approximation

We devise a space-time tensor method for the low-rank approximation of linear parabolic evolution equations. The proposed method is a Galerkin method, uniformly stable in the discretization parameters, based on a Minimal Residual formulation of the evolution problem in Hilbert–Bochner spaces. The discrete solution is sought in a linear trial space composed of tensors of discrete functions in space and in time and is characterized as the unique minimizer of a discrete functional where the dual norm of the residual is evaluated in a space semi-discrete test space. The resulting global space-time linear system is solved iteratively by a greedy algorithm. Numerical results are presented to illustrate the performance of the proposed method on test cases including non-selfadjoint and time-dependent differential operators in space. The results are also compared to those obtained using a fully discrete Petrov–Galerkin setting to evaluate the dual residual norm.

Improved MR image denoising via low- rank approximation and Laplacian-of-Gaussian edge detector

2020 ◽  
Vol 14 (12) ◽  
pp. 2791-2798
Author(s):  
Xiaoqun Qiu ◽  
Zhen Chen ◽  
Saifullah Adnan ◽  
Hongwei He
Keyword(s):  
Image Denoising ◽  
Low Rank ◽  
Edge Detector ◽  
Low Rank Approximation ◽  
Mr Image ◽  
Rank Approximation

scholarly journals Fast Joint Multiband Reconstruction From Wideband Images Based on Low-Rank Approximation

2020 ◽  
Vol 6 ◽  
pp. 922-933
Author(s):  
M. Amine Hadj-Youcef ◽  
Francois Orieux ◽  
Alain Abergel ◽  
Aurelia Fraysse
Keyword(s):  
Low Rank ◽  
Low Rank Approximation ◽  
Rank Approximation

scholarly journals Low-Rank Approximation of Difference between Correlation Matrices Using Inner Product

2021 ◽  
Vol 11 (10) ◽  
pp. 4582
Author(s):  
Kensuke Tanioka ◽  
Satoru Hiwa
Keyword(s):  
Dimensional Reduction ◽  
Low Rank ◽  
Inner Product ◽  
Low Rank Approximation ◽  
Correlation Matrices ◽  
Rank Matrix ◽  
Clustering Approach ◽  
Rank Approximation ◽  
The Difference ◽  
Low Rank Matrix

In the domain of functional magnetic resonance imaging (fMRI) data analysis, given two correlation matrices between regions of interest (ROIs) for the same subject, it is important to reveal relatively large differences to ensure accurate interpretation. However, clustering results based only on differences tend to be unsatisfactory and interpreting the features tends to be difficult because the differences likely suffer from noise. Therefore, to overcome these problems, we propose a new approach for dimensional reduction clustering. Methods: Our proposed dimensional reduction clustering approach consists of low-rank approximation and a clustering algorithm. The low-rank matrix, which reflects the difference, is estimated from the inner product of the difference matrix, not only from the difference. In addition, the low-rank matrix is calculated based on the majorize–minimization (MM) algorithm such that the difference is bounded within the range −1 to 1. For the clustering process, ordinal k-means is applied to the estimated low-rank matrix, which emphasizes the clustering structure. Results: Numerical simulations show that, compared with other approaches that are based only on differences, the proposed method provides superior performance in recovering the true clustering structure. Moreover, as demonstrated through a real-data example of brain activity measured via fMRI during the performance of a working memory task, the proposed method can visually provide interpretable community structures consisting of well-known brain functional networks, which can be associated with the human working memory system. Conclusions: The proposed dimensional reduction clustering approach is a very useful tool for revealing and interpreting the differences between correlation matrices, even when the true differences tend to be relatively small.

hIPPYlib

2021 ◽  
Vol 47 (2) ◽  
pp. 1-34
Author(s):  
Umberto Villa ◽  
Noemi Petra ◽  
Omar Ghattas
Keyword(s):  
Inverse Problems ◽  
Large Scale ◽  
Low Rank ◽  
Scalable Algorithms ◽  
Low Rank Approximation ◽  
Dimensional Parameter ◽  
Infinite Dimensional ◽  
Bayesian Inverse Problems ◽  
Rank Approximation ◽  
Extensible Software

We present an extensible software framework, hIPPYlib, for solution of large-scale deterministic and Bayesian inverse problems governed by partial differential equations (PDEs) with (possibly) infinite-dimensional parameter fields (which are high-dimensional after discretization). hIPPYlib overcomes the prohibitively expensive nature of Bayesian inversion for this class of problems by implementing state-of-the-art scalable algorithms for PDE-based inverse problems that exploit the structure of the underlying operators, notably the Hessian of the log-posterior. The key property of the algorithms implemented in hIPPYlib is that the solution of the inverse problem is computed at a cost, measured in linearized forward PDE solves, that is independent of the parameter dimension. The mean of the posterior is approximated by the MAP point, which is found by minimizing the negative log-posterior with an inexact matrix-free Newton-CG method. The posterior covariance is approximated by the inverse of the Hessian of the negative log posterior evaluated at the MAP point. The construction of the posterior covariance is made tractable by invoking a low-rank approximation of the Hessian of the log-likelihood. Scalable tools for sample generation are also discussed. hIPPYlib makes all of these advanced algorithms easily accessible to domain scientists and provides an environment that expedites the development of new algorithms.

Sign in / Sign up

Export Citation Format

Share Document