Related papers: Hierarchically deflated conjugate residual
While classic work in convex-concave min-max optimization relies on average-iterate convergence results, the emergence of nonconvex applications such as training Generative Adversarial Networks has led to renewed interest in last-iterate…
We consider fermionic (Dirac or Majorana) cold thermal relic dark-matter coupling to standard-model particles through the effective dimension-5 Higgs portal operators $\Lambda^{-1} \ \mathcal{O}_{\text{DM}} \cdot H^\dagger H$, where…
This work is about a new two-level solver for Helmholtz equations discretized by finite elements. The method is inspired by two-grid methods for finite-difference Helmholtz problems as well as by previous work on two-level…
We study the robust matrix completion problem for the low-rank Hankel matrix, which detects the sparse corruptions caused by extreme outliers while we try to recover the original Hankel matrix from the partial observation. In this paper, we…
The overlap Dirac operator at nonzero quark chemical potential involves the computation of the sign function of a non-Hermitian matrix. In this talk we present an iterative method, first proposed by us in Ref. [1], which allows for an…
This paper presents a general framework of high-order finite difference (HFD) schemes for the tempered fractional Laplacian (TFL) based on new generating functions obtained from the discrete symbols. Specifically, for sufficiently smooth…
The numerical computation of many hadronic correlation functions is exceedingly difficult due to the exponentially decreasing signal-to-noise ratio with the distance between source and sink. Multilevel integration methods, using independent…
It is shown that the effective 4D lattice Dirac operator of optimal lattice domain-wall fermions with finite N_s (in the fifth dimension) is exponentially local for sufficiently smooth gauge background.
Getting standard multigrid to work efficiently for the high-frequency Helmholtz equation has been an open problem in applied mathematics for years. Much effort has been dedicated to finding solution methods which can use multigrid…
We present numerical results for the 2-flavour Schwinger model with dynamical chiral lattice fermions. We insert an approximately chiral hypercube Dirac operator into the overlap formula to construct the overlap hypercube operator. This is…
According to the necessary requirements for a chirally symmetric Dirac operator, we present a systematic construction of such operators. We formulate a criterion for the hermitian operator which enters the construction such that the doubled…
Starting from the chiral Lagrangian for Wilson fermions at nonzero lattice spacing we have obtained compact expressions for all spectral correlation functions of the Hermitian Wilson Dirac operator in the $\epsilon$-domain of QCD with…
The efficient and accurate QR decomposition for matrices with hierarchical low-rank structures, such as HODLR and hierarchical matrices, has been challenging. Existing structure-exploiting algorithms are prone to numerical instability as…
Lattice simulations of Quantum Chromodynamics (QCD) enable one to calculate the low-energy properties of the strong interaction among quarks and gluons based on the first principle. The most time-consuming part of the numerical simulations…
We propose a discretisation scheme based on the Dirac-Kahler formalism (DK) in which the algebraic relations between continuum operators ${\wedge, d, \star}$ are captured by their discrete analogues, allowing the construction of the…
We introduce a multigrid multilevel Monte Carlo method for stochastic trace estimation in lattice QCD based on orthogonal projections. This formulation extends the previously proposed oblique decomposition and it is assessed on three…
We generalize the technique of [Solving Dirichlet boundary-value problems on curved domains by extensions from subdomains, SIAM J. Sci. Comput. 34, pp. A497--A519 (2012)] to elliptic problems with mixed boundary conditions and elliptic…
This paper develops a high-order selective discontinuous Galerkin (SDG) method for solving elliptic interface problems on interface-unfitted Cartesian meshes. This method applies the discontinuous Galerkin (DG) formulation on interface…
This paper investigates a category of constrained fractional optimization problems that emerge in various practical applications. The objective function for this category is characterized by the ratio of a numerator and denominator, both…
This paper presents a scalable multigrid preconditioner targeting large-scale systems arising from discontinuous Petrov-Galerkin (DPG) discretizations of high-frequency wave operators. This work is built on previously developed multigrid…