Noisy Accelerated Power Method for Eigenproblems with Applications

03/20/2019
by   Vien V. Mai, et al.
0

This paper introduces an efficient algorithm for finding the dominant generalized eigenvectors of a pair of symmetric matrices. Combining tools from approximation theory and convex optimization, we develop a simple scalable algorithm with strong theoretical performance guarantees. More precisely, the algorithm retains the simplicity of the well-known power method but enjoys the asymptotic iteration complexity of the powerful Lanczos method. Unlike these classic techniques, our algorithm is designed to decompose the overall problem into a series of subproblems that only need to be solved approximately. The combination of good initializations, fast iterative solvers, and appropriate error control in solving the subproblems lead to a linear running time in the input sizes compared to the superlinear time for the traditional methods. The improved running time immediately offers acceleration for several applications. As an example, we demonstrate how the proposed algorithm can be used to accelerate canonical correlation analysis, which is a fundamental statistical tool for learning of a low-dimensional representation of high-dimensional objects. Numerical experiments on real-world data sets confirm that our approach yields significant improvements over the current state-of-the-art.

READ FULL TEXT

page 1

page 2

page 3

page 4

research
04/13/2016

Efficient Algorithms for Large-scale Generalized Eigenvector Computation and Canonical Correlation Analysis

This paper considers the problem of canonical-correlation analysis (CCA)...
research
04/11/2022

Optimizing a low-dimensional convex function over a high-dimensional cube

For a matrix W ∈ℤ^m × n, m ≤ n, and a convex function g: ℝ^m →ℝ, we are ...
research
10/14/2020

Polar Deconvolution of Mixed Signals

The signal demixing problem seeks to separate the superposition of multi...
research
11/06/2020

A dynamic programming approach for generalized nearly isotonic optimization

Shape restricted statistical estimation problems have been extensively s...
research
12/12/2018

Massively scalable Sinkhorn distances via the Nyström method

The Sinkhorn distance, a variant of the Wasserstein distance with entrop...
research
06/26/2023

An orthogonalization-free implementation of the LOBPCG method in solving Kohn-Sham equation

In the classic implementation of the LOBPCG method, orthogonalization an...
research
03/12/2009

Adaptive Mesh Approach for Predicting Algorithm Behavior with Application to Visibility Culling in Computer Graphics

We propose a concise approximate description, and a method for efficient...

Please sign up or login with your details

Forgot password? Click here to reset