Constrained Optimisation of Rational Functions for Accelerating Subspace Iteration

10/21/2017
by   Konrad Kollnig, et al.
0

Earlier this decade, the so-called FEAST algorithm was released for computing the eigenvalues of a matrix in a given interval. Previously, rational filter functions have been examined as a parameter of FEAST. In this thesis, we expand on existing work with the following contributions: (i) Obtaining well-performing rational filter functions via standard minimisation algorithms, (ii) Obtaining constrained rational filter functions efficiently, and (iii) Improving existing rational filter functions algorithmically. Using our new rational filter functions, FEAST requires up to one quarter fewer iterations on average compared to state-of-art rational filter functions.

READ FULL TEXT
research
10/19/2020

Twice is enough for dangerous eigenvalues

We analyze the stability of a class of eigensolvers that target interior...
research
02/06/2023

Bounds on the eigenvalues of matrix rational functions

Upper and lower bounds on absolute values of the eigenvalues of a matrix...
research
03/08/2021

Fast randomized non-Hermitian eigensolver based on rational filtering and matrix partitioning

This paper describes a set of rational filtering algorithms to compute a...
research
04/11/2017

Non-Linear Least-Squares Optimization of Rational Filters for the Solution of Interior Eigenvalue Problems

Rational filter functions can be used to improve convergence of contour-...
research
10/22/2019

Bicubic partially blended rational quartic surface

This paper investigates some univariate and bivariate constrained interp...
research
04/26/2020

Designing a physically-feasible colour filter to make a camera more colorimetric

Previously, a method has been developed to find the best colour filter f...
research
10/28/2020

Simulating a coin with irrational bias using rational arithmetic

An algorithm is presented that, taking a sequence of unbiased coins as i...

Please sign up or login with your details

Forgot password? Click here to reset