### Top 10 Arxiv Papers Today in Numerical Analysis

##### #1. Accelerating Nonnegative Matrix Factorization Algorithms using Extrapolation
###### Andersen Man Shun Ang, Nicolas Gillis
In this paper, we propose a general framework to accelerate significantly the algorithms for nonnegative matrix factorization (NMF). This framework is inspired from the extrapolation scheme used to accelerate gradient methods in convex optimization and from the method of parallel tangents. However, the use of extrapolation in the context of the two-block exact coordinate descent algorithms tackling the non-convex NMF problems is novel. We illustrate the performance of this approach on two state-of-the-art NMF algorithms, namely, accelerated hierarchical alternating least squares (A-HALS) and alternating nonnegative least squares (ANLS), using synthetic, image and document data sets.
more | pdf | html
None.
None.
###### Other stats
Sample Sizes : None.
Authors: 2
Total Words: 8005
Unqiue Words: 2015

##### #2. Ideal Preconditioners for Saddle Point Systems with a Rank-Deficient Leading Block
We consider the iterative solution of symmetric saddle point systems with a rank-deficient leading block. We develop two preconditioners that, under certain assumptions on the rank structure of the system, yield a preconditioned matrix with a constant number of eigenvalues. We then derive some properties of the inverse of a particular class of saddle point system and exploit these to develop a third preconditioner, which remains ideal even when the earlier assumptions on rank structure are relaxed.
more | pdf | html
None.
None.
None.
###### Other stats
Sample Sizes : None.
Authors: 1
Total Words: 5952
Unqiue Words: 1184

##### #3. The trouble with tensor ring decompositions
###### Kim Batselier
The tensor train decomposition decomposes a tensor into a "train" of 3-way tensors that are interconnected through the summation of auxiliary indices. The decomposition is stable, has a well-defined notion of rank and enables the user to perform various linear algebra operations on vectors and matrices of exponential size in a computationally efficient manner. The tensor ring decomposition replaces the train by a ring through the introduction of one additional auxiliary variable. This article discusses a major issue with the tensor ring decomposition: its inability to compute an exact minimal-rank decomposition from a decomposition with sub-optimal ranks. Both the contraction operation and Hadamard product are motivated from applications and it is shown through simple examples how the tensor ring-rounding procedure fails to retrieve minimal-rank decompositions with these operations. These observations, together with the already known issue of not being able to find a best low-rank tensor ring approximation to a given tensor...
more | pdf | html
###### Tweets
ComputerPapers: The trouble with tensor ring decompositions. https://t.co/mEjCDV2S2W
None.
None.
###### Other stats
Sample Sizes : None.
Authors: 1
Total Words: 4373
Unqiue Words: 1390

##### #4. Solutions of Quadratic First-Order ODEs applied to Computer Vision Problems
###### David Casillas-Perez, Daniel Pizarro
The article proves the existence of a maximum of two possible solutions to the initial value problem composed by the planar-perspective equation and an initial condition. This initial value problem has a geometric interpretation. Solutions are curves than pass trough the initial condition which is a point of the plane.
more | pdf | html
None.
None.
###### Other stats
Sample Sizes : [0, 1, 2, 3, 4, 5, 6, 0, 1, 2, 3, 4, 6, 5, 2, 3, 4, 5, 6, 1, 0]
Authors: 2
Total Words: 9311
Unqiue Words: 1793

##### #5. Fast and Accurate Tensor Completion with Tensor Trains: A System Identification Approach
###### Ching-Yun Ko, Kim Batselier, Wenjian Yu, Ngai Wong
This paper proposes a new tensor completion method based on tensor trains and system identification. The to be completed tensor is modeled as a low-rank tensor train system and the coordinates and corresponding tensor entries are interpreted as system inputs and outputs, respectively. A novel tensor train initialization procedure is proposed specifically for image and video completion, which ensures faster convergence of the completion algorithm. The tensor train framework is also shown to easily accommodate Total Variation and Tikhonov regularization due to their low-rank tensor train representations. Image and video inpainting experiments verify the superiority of the proposed scheme in terms of both speed, accuracy, and scalability, where a speedup of up to 60X is observed compared to state-of-the-art tensor completion methods at a similar accuracy.
more | pdf | html
###### Tweets
arxiv_cscv: Fast and Accurate Tensor Completion with Total Variation Regularized Tensor Trains https://t.co/N07dJL95IF
arxiv_cscv: Fast and Accurate Tensor Completion with Total Variation Regularized Tensor Trains https://t.co/N07dJL95IF
###### Github

Fast and accurate tensor completion using tensor trains with a system identification approach

Repository: TTC
User: IRENEKO
Language: Matlab
Stargazers: 1
Subscribers: 0
Forks: 1
Open Issues: 0
views: 539575
likes: 1303
dislikes: 49
favorites: 0
###### Other stats
Sample Sizes : None.
Authors: 4
Total Words: 10303
Unqiue Words: 2708

##### #6. A Class of Multirate Infinitesimal GARK Methods
Differential equations arising in many practical applications are characterized by multiple time scales. Multirate time integration seeks to solve them efficiently by discretizing each scale with a different, appropriate time step, while ensuring the overall accuracy and stability of the numerical solution. In a seminal paper Knoth and Wolke (APNUM, 1998) proposed a hybrid solution approach: discretize the slow component with an explicit Runge-Kutta method, and advance the fast component via a modified fast differential equation. The idea led to the development of multirate infinitesimal step (MIS) methods by Wensch et al. (BIT, 2009.)G\"{u}nther and Sandu (BIT, 2016) explained MIS schemes as a particular case of multirate General-structure Additive Runge-Kutta (MR-GARK) methods. The hybrid approach offers extreme flexibility in the choice of the numerical solution process for the fast component. This work constructs a family of multirate infinitesimal GARK schemes (MRI-GARK) that extends the hybrid dynamics approachin multiple...
more | pdf | html
###### Tweets
MathPaper: A Class of Multirate Infinitesimal GARK Methods. https://t.co/jFVMniz6Mp
None.
None.
###### Other stats
Sample Sizes : None.
Authors: 1
Total Words: 13495
Unqiue Words: 2515

##### #7. Random Walk Laplacian and Network Centrality Measures
###### Daniel Boley, Alejandro Buendia, Golshan Golnari
Random walks over directed graphs are used to model activities in many domains, such as social networks, influence propagation, and Bayesian graphical models. They are often used to compute the importance or centrality of individual nodes according to a variety of different criteria. Here we show how the pseudoinverse of the "random walk" Laplacian can be used to quickly compute measures such as the average number of visits to a given node and various centrality and betweenness measures for individual nodes, both for the network in general and in the case a subset of nodes is to be avoided. We show that with a single matrix inversion it is possible to rapidly compute many such quantities.
more | pdf | html
None.
###### Tweets
ComputerPapers: Random Walk Laplacian and Network Centrality Measures. https://t.co/qoSDsM29Ei
SRoyLee: Random Walk Laplacian and Network Centrality Measures - https://t.co/9MSZ77hPRo
None.
None.
###### Other stats
Sample Sizes : None.
Authors: 3
Total Words: 5411
Unqiue Words: 1367

##### #8. Choosing the optimal multi-point iterative method for the Colebrook flow friction equation -- Numerical validation
###### Pavel Praks, Dejan Brkic
The Colebrook equation $\zeta$ is implicitly given in respect to the unknown flow friction factor $\lambda$; $\lambda=\zeta(Re,\epsilon^*,\lambda)$ which cannot be expressed explicitly in exact way without simplifications and use of approximate calculus. Common approach to solve it is through the Newton-Raphson iterative procedure or through the fixed-point iterative procedure. Both requires in some case even eight iterations. On the other hand numerous more powerful iterative methods such as three-or two-point methods, etc. are available. The purpose is to choose optimal iterative method in order to solve the implicit Colebrook equation for flow friction accurately using the least possible number of iterations. The methods are thoroughly tested and those which require the least possible number of iterations to reach the accurate solution are identified. The most powerful three-point methods require in worst case only two iterations to reach final solution. The recommended representatives are Sharma-Guha-Gupta, Sharma-Sharma,...
more | pdf | html
None.
###### Tweets
M157q_News_RSS: Choosing the optimal multi-point iterative method for the Colebrook flow friction equation -- Numerical validation. https://t.co/dzUJb98STs
ComputerPapers: Choosing the optimal multi-point iterative method for the Colebrook flow friction equation -- Numerical validation. https://t.co/YvSbRDGcJK
None.
None.
###### Other stats
Sample Sizes : None.
Authors: 2
Total Words: 5640
Unqiue Words: 1722

##### #9. Minimizing convex quadratic with variable precision Krylov methods
###### S. Gratton, E. Simon, Ph. L. Toint
Iterative algorithms for the solution of convex quadratic optimization problems are investigated, which exploit inaccurate matrix-vector products. Theoretical bounds on the performance of a Conjugate Gradients and a Full-Orthormalization methods are derived, the necessary quantities occurring in the theoretical bounds estimated and new practical algorithms derived. Numerical experiments suggest that the new methods have significant potential, including in the steadily more important context of multi-precision computations.
more | pdf | html
None.
###### Tweets
ComputerPapers: Minimizing convex quadratic with variable precision Krylov methods. https://t.co/WrwStZ4DK6
None.
None.
###### Other stats
Sample Sizes : None.
Authors: 3
Total Words: 12603
Unqiue Words: 3383

##### #10. Fast Matrix Inversion and Determinant Computation for Polarimetric Synthetic Aperture Radar
###### D. F. G. Coelho, R. J. Cintra, A. C. Frery, V. S. Dimitrov
This paper introduces a fast algorithm for simultaneous inversion and determinant computation of small sized matrices in the context of fully Polarimetric Synthetic Aperture Radar (PolSAR) image processing and analysis. The proposed fast algorithm is based on the computation of the adjoint matrix and the symmetry of the input matrix. The algorithm is implemented in a general purpose graphical processing unit (GPGPU) and compared to the usual approach based on Cholesky factorization. The assessment with simulated observations and data from an actual PolSAR sensor show a speedup factor of about two when compared to the usual Cholesky factorization. Moreover, the expressions provided here can be implemented in any platform.
more | pdf | html
None.
None.
None.
###### Other stats
Sample Sizes : None.
Authors: 4
Total Words: 5037
Unqiue Words: 1815

Assert is a website where the best academic papers on arXiv (computer science, math, physics), bioRxiv (biology), BITSS (reproducibility), EarthArXiv (earth science), engrXiv (engineering), LawArXiv (law), PsyArXiv (psychology), SocArXiv (social science), and SportRxiv (sport research) bubble to the top each day.

Papers are scored (in real-time) based on how verifiable they are (as determined by their Github repos) and how interesting they are (based on Twitter).

To see top papers, follow us on twitter @assertpub_ (arXiv), @assert_pub (bioRxiv), and @assertpub_dev (everything else).

To see beautiful figures extracted from papers, follow us on Instagram.

Tracking 57,756 papers.

###### Search
Sort results based on if they are interesting or reproducible.
Interesting
Reproducible
Online
###### Stats
Tracking 57,756 papers.