Related papers: A nested Krylov subspace method for the overlap op…
We propose a novel Krylov subspace method for estimating the finite impulse response (FIR) of a one-dimensional linear time-invariant systems. The method approximates the system's FIR using a kernel-based formulation combined with…
Consider a Dirac operator on an oriented compact surface endowed with a Riemannian metric and spin structure. Provided the area and the conformal class are fixed, how small can the $k$-th positive Dirac eigenvalue be? This problem mirrors…
An algorithm for unconstrained non-convex optimization is described, which does not evaluate the objective function and in which minimization is carried out, at each iteration, within a randomly selected subspace. It is shown that this…
It is shown in this paper that, almost all current prevalent iterative \mbox{methods} for solving linear system of equations can be classified as what we called extended Krylov subspace methods. In this paper a new type of iterative methods…
In a system where chiral symmetry is spontaneously broken, the low energy eigenmodes of the continuum Dirac operator are extended. On the lattice, due to discretization effects, the Dirac operator can have localized eigenmodes that affect…
Performing Bayesian inference on large spatio-temporal models requires extracting inverse elements of large sparse precision matrices for marginal variances, as well as estimating model hyperparameters. Although direct matrix factorizations…
We revisit quenched reduction with fermions and explain how some old problems can be avoided using the overlap Dirac operator.
We present a theoretical foundation for the Index theorem in naive and minimally doubled lattice fermions by studying the spectral flow of a Hermitean version of Dirac operators. We utilize the point splitting method to implement flavored…
We consider algebraic iterative reconstruction methods with applications in image reconstruction. In particular, we are concerned with methods based on an unmatched projector/backprojector pair; i.e., the backprojector is not the exact…
The Fr\'echet derivative $L_f(A,E)$ of the matrix function $f(A)$ plays an important role in many different applications, including condition number estimation and network analysis. We present several different Krylov subspace methods for…
In this proceeding we propose a new procedure to impose the Schroedinger functional Dirichlet boundary condition on the overlap Dirac operator and the domain-wall fermion using an orbifolding projection. With this procedure the zero mode…
We use the overlap formalism to define a topological index on the lattice. We study the spectral flow of the hermitian Wilson-Dirac operator and identify zero crossings with topological objects. We determine the topological susceptibility…
Block Krylov subspace methods (KSMs) comprise building blocks in many state-of-the-art solvers for large-scale matrix equations as they arise, e.g., from the discretization of partial differential equations. While extended and rational…
We present relaxation and preconditioning techniques which accelerate the inversion of the overlap operator by a factor of four on small lattices, with larger gains as the lattice size increases. These improvements can be used in both…
We present a new algorithmic paradigm for the decentralized solution of graph-structured optimization problems that arise in the estimation and control of network systems. A key and novel design concept of the proposed approach is that it…
We propose an adaptive randomized truncation estimator for Krylov subspace methods that optimizes the trade-off between the solution variance and the computational cost, while remaining unbiased. The estimator solves a constrained…
I apply a recently developed algorithm for reweighting simulations of lattice QCD from one quark mass to another to simulations performed with overlap fermions in the epsilon regime. I test it by computing the condensate from distributions…
We present simulation results for the 2-flavour Schwinger model with dynamical Ginsparg-Wilson fermions. Our Dirac operator is constructed by inserting an approximately chiral hypercube operator into the overlap formula, which yields the…
We study structure-preserving Krylov subspace methods for approximating the matrix-vector products f(H)b, where H is a large Hamiltonian matrix and f denotes either the matrix exponential or the related phi-function. Such computations are…
This is the first of two papers to describe a matrix sparsification algorithm that takes a general real or complex matrix as input and produces a sparse output matrix of the same size. The non-zero entries in the output are chosen to…