BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

James Brannick

James Brannick appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.NAhep-latastro-ph.SRNumerical Analysis

Catalog footprint

What is connected

13works
4topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

13 published item(s)

preprint2023arXiv

A Bootstrap Multigrid Eigensolver

This paper introduces bootstrap multigrid methods for solving eigenvalue problems arising from the discretization of partial differential equations. Inspired by the full bootstrap algebraic multigrid (BAMG) setup algorithm that includes an AMG eigensolver, it is illustrated how the algorithm can be simplified for the case of a discretized partial differential equation (PDE), thereby developing a bootstrap geometric multigrid (BMG) approach. We illustrate numerically the efficacy of the BMG method for: (1) recovering eigenvalues having large multiplicity, (2) computing interior eigenvalues, and (3) approximating shifted indefinite eigenvalue problems. Numerical experiments are presented to illustrate the basic components and ideas behind the success of the overall bootstrap multigrid approach. For completeness, we present a simplified error analysis of a two-grid bootstrap algorithm for the Laplace-Beltrami eigenvalue problem.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Microphysics of Neutron Star Outer Envelopes in the Periodized, Magnetic Thomas-Fermi Model

Static and dynamic properties of low density outer envelopes of neutron stars are calculated within the nonlinear magnetic Thomas-Fermi model, assuming degenerate electrons. A novel domain decomposition enables proper description of lattice symmetry and may be seen as a prototype for the general class of problems involving nonlinear charge screening of periodic, quasi-low-dimensionality structures, e.g. liquid crystals. We describe a scalable implementation of the method using Hypre. Phase velocity of long wavelength transverse phonons is found to be a factor of 5-7 larger than in the corresponding Coulomb crystal model, which could have implications for low temperature phonon-mediated thermal conductivity. Other findings include $c'<0$ elastic instabilities for both bcc and fcc lattices, reminiscent of the situation in some light actinides, and suggestive of a symmetry-lowering transition to a tetragonal or orthorhombic lattice.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

An Adaptively Constructed Algebraic Multigrid Preconditioner for Irreducible Markov Chains

The computation of stationary distributions of Markov chains is an important task in the simulation of stochastic models. The linear systems arising in such applications involve non-symmetric M-matrices, making algebraic multigrid methods a natural choice for solving these systems. In this paper we investigate extensions and improvements of the bootstrap algebraic multigrid framework for solving these systems. This is achieved by reworking the bootstrap setup process to use singular vectors instead of eigenvectors in constructing interpolation and restriction. We formulate a result concerning the convergence speed of GMRES for singular systems and experimentally justify why rapid convergence of the proposed method can be expected. We demonstrate its fast convergence and the favorable scaling behavior for various test problems.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Bootstrap Algebraic Multigrid: status report, open problems, and outlook

This paper provides an overview of the main ideas driving the bootstrap algebraic multigrid methodology, including compatible relaxation and algebraic distances for defining effective coarsening strategies, the least squares method for computing accurate prolongation operators and the bootstrap cycles for computing the test vectors that are used in the least squares process. We review some recent research in the development, analysis and application of bootstrap algebraic multigrid and point to open problems in these areas. Results from our previous research as well as some new results for some model diffusion problems with highly oscillatory diffusion coefficient are presented to illustrate the basic components of the BAMG algorithm.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Local Fourier Analysis of Multigrid Methods with Polynomial Smoothers and Aggressive coarsening

We focus on the study of multigrid methods with aggressive coarsening and polynomial smoothers for the solution of the linear systems corresponding to finite difference/element discretizations of the Laplace equation. Using local Fourier analysis we determine automatically the optimal values for the parameters involved in defining the polynomial smoothers and achieve fast convergence of cycles with aggressive coarsening. We also present numerical tests supporting the theoretical results and the heuristic ideas. The methods we introduce are highly parallelizable and efficient multigrid algorithms on structured and semi-structured grids in two and three spatial dimensions.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Multigrid Preconditioning for the Overlap Operator in Lattice QCD

The overlap operator is a lattice discretization of the Dirac operator of quantum chromodynamics, the fundamental physical theory of the strong interaction between the quarks. As opposed to other discretizations it preserves the important physical property of chiral symmetry, at the expense of requiring much more effort when solving systems with this operator. We present a preconditioning technique based on another lattice discretization, the Wilson-Dirac operator. The mathematical analysis precisely describes the effect of this preconditioning in the case that the Wilson-Dirac operator is normal. Although this is not exactly the case in realistic settings, we show that current smearing techniques indeed drive the Wilson-Dirac operator towards normality, thus providing a motivation why our preconditioner works well in computational practice. Results of numerical experiments in physically relevant settings show that our preconditioning yields accelerations of up to one order of magnitude.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Aggregation-based aggressive coarsening with polynomial smoothing

This paper develops an algebraic multigrid preconditioner for the graph Laplacian. The proposed approach uses aggressive coarsening based on the aggregation framework in the setup phase and a polynomial smoother with sufficiently large degree within a (nonlinear) Algebraic Multilevel Iteration as a preconditioner to the flexible Conjugate Gradient iteration in the solve phase. We show that by combining these techniques it is possible to design a simple and scalable algorithm. Results of the algorithm applied to graph Laplacian systems arising from the standard linear finite element discretization of the scalar Poisson problem are reported.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Bootstrap Algebraic Multigrid for the 2D Wilson Dirac system

We develop an algebraic multigrid method for solving the non-Hermitian Wilson discretization of the 2-dimensional Dirac equation. The proposed approach uses a bootstrap setup algorithm based on a multigrid eigensolver. It computes test vectors which define the least squares interpolation operators by working mainly on coarse grids, leading to an efficient and integrated self learning process for defining algebraic multigrid interpolation. The algorithm is motivated by the γ-symmetry of the Dirac equation, which carries over to the Wilson discretization. This discrete γ-symmetry is used to reduce a general Petrov Galerkin bootstrap setup algorithm to a Galerkin method for the Hermitian and indefinite formulation of the Wilson matrix. Kaczmarz relaxation is used as the multigrid smoothing scheme in both the setup and solve phases of the resulting Galerkin algorithm. The overall method is applied to the odd-even reduced Wilson matrix, which also fulfills the discrete γ-symmetry. Extensive numerical results are presented to motivate the design and demonstrate the effectiveness of the proposed approach.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Parallel Unsmoothed Aggregation Algebraic Multigrid Algorithms on GPUs

We design and implement a parallel algebraic multigrid method for isotropic graph Laplacian problems on multicore Graphical Processing Units (GPUs). The proposed AMG method is based on the aggregation framework. The setup phase of the algorithm uses a parallel maximal independent set algorithm in forming aggregates and the resulting coarse level hierarchy is then used in a K-cycle iteration solve phase with a $\ell^1$-Jacobi smoother. Numerical tests of a parallel implementation of the method for graphics processors are presented to demonstrate its effectiveness.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Algebraic multilevel preconditioners for the graph Laplacian based on matching in graphs

This paper presents estimates of the convergence rate and complexity of an algebraic multilevel preconditioner based on piecewise constant coarse vector spaces applied to the graph Laplacian. A bound is derived on the energy norm of the projection operator onto any piecewise constant vector space, which results in an estimate of the two-level convergence rate where the coarse level graph is obtained by matching. The two-level convergence of the method is then used to establish the convergence of an Algebraic Multilevel Iteration that uses the two-level scheme recursively. On structured grids, the method is proven to have convergence rate $\approx (1-1/\log n)$ and $O(n\log n)$ complexity for each cycle, where $n$ denotes the number of unknowns in the given problem. Numerical results of the algorithm applied to various graph Laplacians are reported. It is also shown that all the theoretical estimates derived for matching can be generalized to the case of aggregates containing more than two vertices.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

An algebraic distances measure of AMG strength of connection

Algebraic multigrid is an iterative method that is often optimal for solving the matrix equations that arise in a wide variety of applications, including discretized partial differential equations. It automatically constructs a sequence of increasingly smaller matrix problems that enable efficient resolution of all scales present in the solution. One of the main components of the method is an adequate choice of coarse grids. The current coarsening methodology is based on measuring how a so-called algebraically smooth error value at one point depends on the error values at its neighbors. Such a concept of strength of connection is well understood for operators whose principal part is an M-matrix; however, the strength concept for more general matrices is not yet clearly understood, and this lack of knowledge limits the scope of AMG applicability. The purpose of this paper is to motivate a general definition of strength of connection, based on the notion of algebraic distances, discuss its implementation, and present the results of initial numerical experiments. The algebraic distance measure, we propose, uses as its main tool a least squares functional, which is also applied to define interpolation.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The role of multigrid algorithms for LQCD

We report on the first successful QCD multigrid algorithm which demonstrates constant convergence rates independent of quark mass and lattice volume for the Wilson Dirac operator. The new ingredient is the adaptive method for constructing the near null space on which the coarse grid multigrid Dirac operator acts. In addition we speculate on future prospects for extending this algorithm to the Domain Wall and Staggered discretizations, its exceptional suitability for high performance GPU code and its potential impact on simulations at the physical pion mass.

0 citations0 reviewsNo code linkedOpen access