Tuning Multigrid Methods with Robust Optimization

by   Jed Brown, et al.

Local Fourier analysis is a useful tool for predicting and analyzing the performance of many efficient algorithms for the solution of discretized PDEs, such as multigrid and domain decomposition methods. The crucial aspect of local Fourier analysis is that it can be used to minimize an estimate of the spectral radius of a stationary iteration, or the condition number of a preconditioned system, in terms of a symbol representation of the algorithm. In practice, this is a "minimax" problem, minimizing with respect to solver parameters the appropriate measure of work, which involves maximizing over the Fourier frequency. Often, several algorithmic parameters may be determined by local Fourier analysis in order to obtain efficient algorithms. Analytical solutions to minimax problems are rarely possible beyond simple problems; the status quo in local Fourier analysis involves grid sampling, which is prohibitively expensive in high dimensions. In this paper, we propose and explore optimization algorithms to solve these problems efficiently. Several examples, with known and unknown analytical solutions, are presented to show the effectiveness of these approaches.



There are no comments yet.


page 15


A local Fourier analysis of additive Vanka relaxation for the Stokes equations

Multigrid methods are popular solution algorithms for many discretized P...

Local Fourier analysis of Balancing Domain Decomposition by Constraints algorithms

Local Fourier analysis is a commonly used tool for the analysis of multi...

Fourier spectral methods for nonlocal models

Efficient and accurate spectral solvers for nonlocal models in any spati...

Fourier-Domain Optimization for Image Processing

Image optimization problems encompass many applications such as spectral...

Trigonometric splines in spectral problems

Some questions of application of trigonometric splines in problems of sp...

Methods to Deal with Unknown Populational Minima during Parameter Inference

There is a myriad of phenomena that are better modelled with semi-infini...

IsoGeometric Approximations for Electromagnetic Problems in Axisymmetric Domains

We propose a numerical method for the solution of electromagnetic proble...
This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.