We present an adaptive multigrid solver for application to the non-Hermitian Wilson-Dirac system of QCD. The key components leading to the success of our proposed algorithm are the use of an adaptive projection onto coarse grids that preserves the near null space of the system matrix together with a simplified form of the correction based on the so-called γ5-Hermitian symmetry of the Dirac operator. We demonstrate that the algorithm nearly eliminates critical slowing down in the chiral limit and that it has weak dependence on the lattice volume.
Download full-text PDF |
Source |
|---|---|
| http://dx.doi.org/10.1103/PhysRevLett.105.201602 | DOI Listing |
Publication Analysis
Top Keywords
Similar Publications
AdaCS: Adaptive Compressive Sensing With Restricted Isometry Property-Based Error-Clamping.
IEEE Trans Pattern Anal Mach Intell
July 2024
Scene-dependent adaptive compressive sensing (CS) has been a long pursuing goal that has huge potential to significantly improve the performance of CS. However, with no access to the ground truth, how to design the scene-dependent adaptive strategy is still an open problem. In this paper, a restricted isometry property (RIP) condition-based error-clamping is proposed, which could directly predict the reconstruction error, i.
View Article and Find Full Text PDFAdaptively Placed Multi-Grid Scene Representation Networks for Large-Scale Data Visualization.
IEEE Trans Vis Comput Graph
January 2024
Scene representation networks (SRNs) have been recently proposed for compression and visualization of scientific data. However, state-of-the-art SRNs do not adapt the allocation of available network parameters to the complex features found in scientific data, leading to a loss in reconstruction quality. We address this shortcoming with an adaptively placed multi-grid SRN (APMGSRN) and propose a domain decomposition training and inference technique for accelerated parallel training on multi-GPU systems.
View Article and Find Full Text PDFGoal-oriented adaptive finite element methods with optimal computational complexity.
Numer Math (Heidelb)
November 2022
TU Wien, Institute of Analysis and Scientific Computing, Wiedner Hauptstr. 8-10/E101/4, Vienna, 1040 Austria.
We consider a linear symmetric and elliptic PDE and a linear goal functional. We design and analyze a goal-oriented adaptive finite element method, which steers the adaptive mesh-refinement as well as the approximate solution of the arising linear systems by means of a contractive iterative solver like the optimally preconditioned conjugate gradient method or geometric multigrid. We prove linear convergence of the proposed adaptive algorithm with optimal algebraic rates.
View Article and Find Full Text PDFContinuum limits of discrete isoperimetric problems and Wulff shapes in lattices and quasicrystal tilings.
Calc Var Partial Differ Equ
October 2022
Facultad de Matemáticas, Pontificia Universidad Católica de Chile, Avda. Vicuña Mackenna 4860, 6904441 Macul, Santiago Chile.
We prove discrete-to-continuum convergence of interaction energies defined on lattices in the Euclidean space (with interactions beyond nearest neighbours) to a crystalline perimeter, and we discuss the possible Wulff shapes obtainable in this way. Exploiting the "multigrid construction" of quasiperiodic tilings (which is an extension of De Bruijn's "pentagrid" construction of Penrose tilings) we adapt the same techniques to also find the macroscopical homogenized perimeter when we microscopically rescale a given quasiperiodic tiling.
View Article and Find Full Text PDFSimulation of Multispecies Desmoplastic Cancer Growth via a Fully Adaptive Non-linear Full Multigrid Algorithm.
Front Physiol
July 2018
Department of Bioengineering, University of Louisville, Louisville, KY, United States.
A fully adaptive non-linear full multigrid (FMG) algorithm is implemented to computationally simulate a model of multispecies desmoplastic tumor growth in three spatial dimensions. The algorithm solves a thermodynamic mixture model employing a diffuse interface approach with Cahn-Hilliard-type fourth-order equations that are coupled, non-linear, and numerically stiff. The tumor model includes extracellular matrix (ECM) as a major component with elastic energy contribution in its chemical potential term.
View Article and Find Full Text PDFWant AI Summaries of new PubMed Abstracts delivered to your In-box?
Enter search terms and have AI summaries delivered each week - change queries or unsubscribe any time!