
Multigrid MethodsSince there exists a simple subdivision scheme for bsplines, websplines are ideally suited for multigrid techniques. A sequence of nested grids of width h, 2h, 4h, ... is used to reduce the residual error of the approximated solution and accelerate the iterative solution of the Galerkin system. As a model problem, we consider Poisson's equation div(grad(u))=1 with homogeneous Dirichlet boundary conditions on the following domain (left). The affiliated solution u is depicted on the right. Bsplines of degree 2 have been used on the following grids of width h=0.5, 1, and 2. Outer bsplines are marked red, extended bsplines green, and unextended bsplines blue. Grid transfer of correction vectors between the three webspaces of dimension 966, 260, and 78 can efficiently be performed with the aid of sparse projection matrices with the following sparsity pattern: Results of the dynamic multigrid solver (dmg) and a standard ssor preconditioned conjugate gradient solver (pcg) to achieve a residual norm <1e10 are compared in the following table:
Each row represents one complete test using all grids of width 2^{n},2^{n+1},...,2 in case of the dmg solver, respectively one grid of width 2^{n} in case of the pcg solver. The second column lists the number of basis functions for the given grid width. Where pcg iter represents the number of pcg iterations, dmg iter represents the weighted number of ssor smoothing steps. A smoothing step on the finest grid of width h is weighted by 1, on the grid of width 2h by 1/4 and so on. The rate is computed by the iterth root of the last residual norm. Since dmg smoothing steps are less expensive than pcg steps, the rates are not directly comparable. This is reflected by the ratio of solver times, which shows the immense speedup by using the dynamic multigrid solver. LiteratureK. Höllig, U. Reif, J. Wipper: Multigrid methods with websplines. Numerische Mathematik, Volume 91, Number 2, pp. 237256, 2002.See Publications for a complete list of WEBpublications.
Author: Joachim Wipper ; Last modification: 2006/08/31 09:38:06 UTC.
