GMGPolar v1.0.3

Kurzfassung

GMGPolar offers a taylored Geometric Multigrid solver for geometries described by curvilinear coordinates (in its simplest case, by polar coordinates). GMGPolar uses an implicit extrapolation scheme to raise the convergence order. For nonstandard finite elements, this result has been shown theoretically and for nine-point finite difference stars, we observed up to convergence order four in all numerical simulations. This extrapolation scheme can be used likewise for other domains described by tensor- or product-format meshes in 2D. For curvilinear domains, suited smoothing schemes have been developed and implemented. GMGPolar offers a matrix and matrix-free version for most of its operators. Its original application is on cross-sectional views of realistic Tokamak geometries for fusion plasma simulations.