PRIMME_SVDS: A High-Performance Preconditioned SVD Solver for Accurate Large-Scale Computations

07/05/2016
by   Lingfei Wu, et al.
0

The increasing number of applications requiring the solution of large scale singular value problems have rekindled interest in iterative methods for the SVD. Some promising recent ad- vances in large scale iterative methods are still plagued by slow convergence and accuracy limitations for computing smallest singular triplets. Furthermore, their current implementations in MATLAB cannot address the required large problems. Recently, we presented a preconditioned, two-stage method to effectively and accurately compute a small number of extreme singular triplets. In this research, we present a high-performance software, PRIMME SVDS, that implements our hybrid method based on the state-of-the-art eigensolver package PRIMME for both largest and smallest singular values. PRIMME SVDS fills a gap in production level software for computing the partial SVD, especially with preconditioning. The numerical experiments demonstrate its superior performance compared to other state-of-the-art software and its good parallel performance under strong and weak scaling.

READ FULL TEXT

page 17

page 21

research
04/03/2022

An iterative formula of smallest singular value of nonsingular matrices

We obtain an iterative formula that converges incrementally to the small...
research
12/15/2022

A skew-symmetric Lanczos bidiagonalization method for computing several largest eigenpairs of a large skew-symmetric matrix

The spectral decomposition of a real skew-symmetric matrix A can be math...
research
04/29/2021

High-Performance Partial Spectrum Computation for Symmetric eigenvalue problems and the SVD

Current dense symmetric eigenvalue (EIG) and singular value decompositio...
research
07/26/2019

Massively Scaling Seismic Processing on Sunway TaihuLight Supercomputer

Common Midpoint (CMP) and Common Reflection Surface (CRS) are widely use...
research
01/05/2023

A tensor bidiagonalization method for higher-order singular value decomposition with applications

The need to know a few singular triplets associated with the largest sin...
research
05/30/2018

Accelerating Large-Scale Data Analysis by Offloading to High-Performance Computing Libraries using Alchemist

Apache Spark is a popular system aimed at the analysis of large data set...
research
04/27/2023

Singularity Distance Computations for 3-RPR Manipulators using Extrinsic Metrics

It is well-known that parallel manipulators are prone to singularities. ...

Please sign up or login with your details

Forgot password? Click here to reset