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-\alpha 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.