7,040 research outputs found
A Self-learning Algebraic Multigrid Method for Extremal Singular Triplets and Eigenpairs
A self-learning algebraic multigrid method for dominant and minimal singular
triplets and eigenpairs is described. The method consists of two multilevel
phases. In the first, multiplicative phase (setup phase), tentative singular
triplets are calculated along with a multigrid hierarchy of interpolation
operators that approximately fit the tentative singular vectors in a collective
and self-learning manner, using multiplicative update formulas. In the second,
additive phase (solve phase), the tentative singular triplets are improved up
to the desired accuracy by using an additive correction scheme with fixed
interpolation operators, combined with a Ritz update. A suitable generalization
of the singular value decomposition is formulated that applies to the coarse
levels of the multilevel cycles. The proposed algorithm combines and extends
two existing multigrid approaches for symmetric positive definite eigenvalue
problems to the case of dominant and minimal singular triplets. Numerical tests
on model problems from different areas show that the algorithm converges to
high accuracy in a modest number of iterations, and is flexible enough to deal
with a variety of problems due to its self-learning properties.Comment: 29 page
Algebraic Multigrid for Disordered Systems and Lattice Gauge Theories
The construction of multigrid operators for disordered linear lattice
operators, in particular the fermion matrix in lattice gauge theories, by means
of algebraic multigrid and block LU decomposition is discussed. In this
formalism, the effective coarse-grid operator is obtained as the Schur
complement of the original matrix. An optimal approximation to it is found by a
numerical optimization procedure akin to Monte Carlo renormalization, resulting
in a generalized (gauge-path dependent) stencil that is easily evaluated for a
given disorder field. Applications to preconditioning and relaxation methods
are investigated.Comment: 43 pages, 14 figures, revtex4 styl
Localization in Lattice Gauge Theory and a New Multigrid Method
We show numerically that the lowest eigenmodes of the 2-dimensional
Laplace-operator with SU(2) gauge couplings are strongly localized. A
connection is drawn to the Anderson-Localization problem. A new Multigrid
algorithm, capable to deal with these modes, shows no critical slowing down for
this problem.Comment: LATeX style, 11 pages (plus 4 figure pages). Figure pages are
available as uuencoded ps-file via anonymous ftp from x4u2.desy.de, get
pub/outgoing/baeker/heplat.uu. DESY-preprint 94-07
A multigrid continuation method for elliptic problems with folds
We introduce a new multigrid continuation method for computing solutions of nonlinear elliptic eigenvalue problems which contain limit points (also called turning points or folds). Our method combines the frozen tau technique of Brandt with pseudo-arc length continuation and correction of the parameter on the coarsest grid. This produces considerable storage savings over direct continuation methods,as well as better initial coarse grid approximations, and avoids complicated algorithms for determining the parameter on finer grids. We provide numerical results for second, fourth and sixth order approximations to the two-parameter, two-dimensional stationary reaction-diffusion problem: Δu+λ exp(u/(1+au)) = 0.
For the higher order interpolations we use bicubic and biquintic splines. The convergence rate is observed to be independent of the occurrence of limit points
Adaptive Aggregation Based Domain Decomposition Multigrid for the Lattice Wilson Dirac Operator
In lattice QCD computations a substantial amount of work is spent in solving
discretized versions of the Dirac equation. Conventional Krylov solvers show
critical slowing down for large system sizes and physically interesting
parameter regions. We present a domain decomposition adaptive algebraic
multigrid method used as a precondtioner to solve the "clover improved" Wilson
discretization of the Dirac equation. This approach combines and improves two
approaches, namely domain decomposition and adaptive algebraic multigrid, that
have been used seperately in lattice QCD before. We show in extensive numerical
test conducted with a parallel production code implementation that considerable
speed-up over conventional Krylov subspace methods, domain decomposition
methods and other hierarchical approaches for realistic system sizes can be
achieved.Comment: Additional comparison to method of arXiv:1011.2775 and to
mixed-precision odd-even preconditioned BiCGStab. Results of numerical
experiments changed slightly due to more systematic use of odd-even
preconditionin
- …