Researcher profile

K. Kahl

K. Kahl contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 17 - UnverifiedVerification L1Unclaimed author
4works
0followers
2topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

4 published item(s)

preprint2014arXiv

Algebraic distance for anisotropic diffusion problems: multilevel results

In this paper we motivate, discuss the implementation and present the resulting numerics for a new definition of strength of connection which is based on the notion of algebraic distance. This algebraic distance measure, combined with compatible relaxation, is used to choose suitable coarse grids and accurate interpolation operators for algebraic multigrid algorithms. The main tool of the proposed measure is the least squares functional defined using a set of relaxed test vectors. The motivating application is the anisotropic diffusion problem, in particular problems with non-grid aligned anisotropy. We demonstrate numerically that the measure yields a robust technique for determining strength of connectivity among variables, for both two-grid and multigrid solvers. %We illustrate the use of the measure to construct, in addition, an adaptive aggregation form of interpolation for the targeted anisotropic problems. %Our approach is not a two-level approach -- we provide preliminary results that show its extendability to multigrid. The proposed algebraic distance measure can also be used in any other coarsening procedure, assuming a rich enough set of the near-kernel components of the matrix for the targeted system is known or computed.

preprint2013arXiv

An adaptive aggregation based domain decomposition multilevel method for the lattice wilson dirac operator: multilevel results

In lattice QCD computations a substantial amount of work is spent in solving linear systems arising in Wilson's discretization of the Dirac equations. We show first numerical results of the extension of the two-level DD-αAMG method to a true multilevel method based on our parallel MPI-C implementation. Using additional levels pays off, allowing to cut down the core minutes spent on one system solve by a factor of approximately 700 compared to standard Krylov subspace methods and yielding another speed-up of a factor of 1.7 over the two-level approach.

preprint2012arXiv

2-norm error bounds and estimates for Lanczos approximations to linear systems and rational matrix functions

The Lanczos process constructs a sequence of orthonormal vectors v_m spanning a nested sequence of Krylov subspaces generated by a hermitian matrix A and some starting vector b. In this paper we show how to cheaply recover a secondary Lanczos process starting at an arbitrary Lanczos vector v_m. This secondary process is then used to efficiently obtain computable error estimates and error bounds for the Lanczos approximations to the action of a rational matrix function on a vector. This includes, as a special case, the Lanczos approximation to the solution of a linear system Ax = b. Our approach uses the relation between the Lanczos process and quadrature as developed by Golub and Meurant. It is different from methods known so far because of its use of the secondary Lanczos process. With our approach, it is now in particular possible to efficiently obtain {\em upper bounds} for the error in the {\em 2-norm}, provided a lower bound on the smallest eigenvalue of $A$ is known. This holds in particular for a large class of rational matrix functions including best rational approximations to the inverse square root and the sign function. We compare our approach to other existing error estimates and bounds known from the literature and include results of several numerical experiments.

preprint2010arXiv

A Bootstrap Algebraic Multilevel method for Markov Chains

This work concerns the development of an Algebraic Multilevel method for computing stationary vectors of Markov chains. We present an efficient Bootstrap Algebraic Multilevel method for this task. In our proposed approach, we employ a multilevel eigensolver, with interpolation built using ideas based on compatible relaxation, algebraic distances, and least squares fitting of test vectors. Our adaptive variational strategy for computation of the state vector of a given Markov chain is then a combination of this multilevel eigensolver and associated multilevel preconditioned GMRES iterations. We show that the Bootstrap AMG eigensolver by itself can efficiently compute accurate approximations to the state vector. An additional benefit of the Bootstrap approach is that it yields an accurate interpolation operator for many other eigenmodes. This in turn allows for the use of the resulting AMG hierarchy to accelerate the MLE steps using standard multigrid correction steps. The proposed approach is applied to a range of test problems, involving non-symmetric stochastic M-matrices, showing promising results for all problems considered.