Questions tagged [multigrid]

An approach to solving systems of equations by projecting the problem from a fine scale representation onto a coarser one. A coarse representation generally has fewer unknowns, making it faster to solve than the original problem. The coarse solution can then be projected back onto the finer problem as an initial guess of the solution to the finer problem.

Filter by
Sorted by
Tagged with
0 votes
0 answers
38 views

Multigrid/Two-grid method restriction and prolongation of residual

Starting from the problem $Au=f$, I'm not sure that I understand why a coarse grid solution is implemented on the coarse grid residual $r_c^{(k)}=P^Tr^{(k)}$, with $r^{(k)}=f-Au^{(k)}$ and $P^T$ ...
user avatar
0 votes
1 answer
45 views

Interpolation and Restriction operators in Multigrid

I saw in several places that interpolation operator ($P$) and restriction operator ($P^T$) are usually transposes of each other (up to multiplication by a constant). As I understood it related to ...
user avatar
0 votes
2 answers
104 views

Eigenvectors of Laplacian

I am studying introduction to Multigrid methods. In all tutorials, authors write that eigenvectors of Laplacian (1D, finite difference) are given as $w_k(x_i) = \sin(k \pi x_i),$ where $x_i$ is a ...
user avatar
1 vote
1 answer
66 views

Algebraic Multigrid on fine mesh vs. just starting with coarse mesh?

So if I've understood it correctly, the Algebraic Multigrid Method (AMG) basically takes a fine mesh, coarsens it, solves the coarse mesh and projects the solution back on the fine mesh. Wouldn't it ...
user avatar
  • 167
0 votes
0 answers
50 views

Geometric multigrid transfer operators constraints

I have been reading up on constraints used to construct prolongation and restriction operators for geometric multigrid. One such constraint is that the prolongation matrix is the transpose of the ...
user avatar
1 vote
1 answer
95 views

Prolongation and restriction operators in multigrid for high order PDEs

If I have the Poisson equation $\Delta u = f$ a standard transfer operator (for a regular grid) is the full weighting/bilinear interpolation scheme: $$K = \frac{1}{4}\begin{bmatrix}\frac{1}{4} & \...
user avatar
1 vote
1 answer
53 views

Restriction in (geometric) multigrid for vectors of non-even length

Naive restriction operators in geometric multigrid that I have seen are typically implemented as a convolution and a subsequent averaging of every two entries in a vector $v^h$. For example: $$\tilde{...
user avatar
0 votes
0 answers
50 views

Differences between Two-Grid Correction and V-Cycle Scheme

I found in A Multigrid Tutorial: Second Edition by William L. Briggs‏, Van Emden Henson‏, Steve F. (here) the following Schemes: Two-Grid Correction Scheme, p. 37: And V-Cycle Scheme, p. 40 I want ...
user avatar
2 votes
0 answers
66 views

Multigrid method: linear solver and modified residual

I am trying to better understand the FAS multigrid algorithm for Euler equation in FV discretization. The usage of the modified residual (the residual with forcing) inside the different cases: ...
user avatar
  • 382
1 vote
0 answers
45 views

Multigrid: Gauss–Seidel eigenvalues and eigenvectors

I am trying to figure out the questions from Multigrid Tutorial by Briggs. However I am stuck on these two questions. https://www.researchgate.net/publication/...
user avatar
1 vote
1 answer
82 views

Algebraic multigrid as solver and as preconditioner

My question is around the efficiency of AMG. In which case AMG can perform better,as solver or as a preconditioner(for example a Krylov space method as CG)? Assume the case of elliptic pdes.
user avatar
  • 471
4 votes
0 answers
288 views

Construction of Prolongation and Restriction Operator for Geometric Multigrid (2D-FEM): Resulting in a Decreasing Solution

Consider the following problem, $$ -\Delta u(x) = f(x), \qquad x \in \Omega \\ u(x) = 0,\qquad x \in \partial \Omega$$ with $\Omega = [0,1]\times [0,1]$ being the domain and $\partial \Omega$ being ...
user avatar
  • 41
1 vote
0 answers
52 views

Multigrid Reduction In Time Convergence

I am trying to solve a 2D dynamic linear elasticity model parallel in time using Xbraid. The spatial domain is [0,1]x[0,1] and time domain [0,1]. For time integration I am using a backward Euler ...
user avatar
  • 471
2 votes
0 answers
151 views

Parallel In Time with Multigrid

I am trying to solve the linear finite element equation $M\ddot{u}+Ku=F(t)$, where $M$ is the mass matrix ,$K$ the stiffness matrix and $F(t)$ the external load vector, parallel in time using XBraid ...
user avatar
  • 471
1 vote
1 answer
162 views

coarsening coefficient matrixes (A2h, A4h...) for geometric multigrid method in 2-D/3-D

I am learning about multigrid methods from the textbook section 6.3 Multigrid Methods, which shows a geometric multigrid algorithm for 1-D examples in detail, including how to build restriction/...
user avatar
4 votes
1 answer
142 views

Why does smoothed aggregation multigrid method used as preconditioner in conjugate gradient slows down the solution time?

I'm solving a system of linear equations obtained from the FEM discretization of a simple linear elasticity problem on a cube with zero displacements at one plane and a load on the opposite one. The ...
user avatar
2 votes
1 answer
72 views

Calculating coarse grid matrix in geometric multigird

The coarse grid matrix is calculated via RAP where R,P are the restriction and interpolation matrix,respectively.By checking a typical MG algorithm I want to ask how to calculate efficiently coarse ...
user avatar
  • 471
1 vote
1 answer
111 views

Full approximation scheme - smoothers - literature recomendation

I would like to use full approximation scheme(FAS) for solving nonlinear PDE. I am looking into practical implementation of FAS in parallel(MPI) settings. I noticed the most common/effective smoother ...
user avatar
2 votes
0 answers
204 views

Residual of Poisson equation with periodic boundaries

I am trying to write a multigrid solver for Poisson's equation, $-\Delta u=f$, on the unit square, $\Omega=(0,1)^2$ with periodic boundaries. My primary source has been Multigrid by Trottenberg, ...
user avatar
2 votes
1 answer
103 views

How the number of pre-smoothing and post-smoothing steps affect the asymtotic convergence rate of geometrical Multigrid?

Does the convergence rate of multigrid depend on the total number of smoothing steps or on the number of pre and post smoothing steps seperately?
user avatar
  • 471
1 vote
1 answer
55 views

How coarsening rate affects MG convergence?

If we use an aggresive coarsening how does this affect the convergence of a MG cycle?Is the convergence slower?Do we need more cycles in this case?
user avatar
  • 471
2 votes
1 answer
70 views

Why post-smoothing in MG is needed?

In pre-smoothing we eliminate high frequencies,but what about post-smoothing?Is it used for the remain high frequencies that pre-smoothing didn't eliminate?
user avatar
  • 471
2 votes
0 answers
60 views

Avoid matrix multiplication in algebraic multigrid method

Currently when I try to solve a linear algebra system of the form of $A x =b$ I use the algebraic multigrid method. The algebraic multigrid method uses a Galerkin product to form the coarse grid ...
user avatar
1 vote
0 answers
35 views

What kind of problem or matrices are suitable for multigrid method?

For Poisson or Convection-diffusion equation as follows: $$ -\Delta u=f,\qquad u|_\Omega = g. $$ or $$ -\Delta u +\vec{w}.\nabla u=f,\qquad u|_\Omega = g. $$ using FDM or FEM discretization, we can ...
user avatar
  • 961
1 vote
1 answer
85 views

How to implement Geometric Multigrid in non-rectangular grids?

It is quite easy to implement multigrid on a rectangular grid but what about an non-rectangular?How to coarse a non-rectangular grid and apply multigrid(assume an easy non-rectangular grid capital ...
user avatar
  • 471
2 votes
0 answers
68 views

Multilevel minimization - boundary conditions

I am interested in minimizing $$min_{x \in R^{n^l}} f^l(x),$$ where $f^l(x)$ is nonlinear objective function arising from discretization of PDE. I would like to use nonlinear multilevel minimization ...
user avatar
1 vote
1 answer
194 views

How suitable is multigrid method for time-dependent PDEs?

For elliptic PDEs (Poisson-type), the multigrid method is very sufficient, but how about time-dependent problems (i.e parabolic or hyperbolic PDEs)? Is it efficient to solve such problems using a ...
user avatar
  • 471
1 vote
0 answers
60 views

Can the standard multigrid performance be used for time-dependent PDEs?

Consider a time dependent pde(i.e u(x,t)).I know when only space-coarsening is used the standard multigrid performance can be applied but what if instead we use only time-coarsening?Can we apply the ...
user avatar
  • 471
2 votes
2 answers
120 views

Can a direct method like Thomas be used in a multigrid method as a smoother?

As far as I know, multigrid uses stationary iterative methods as smoothers (i.e GS), but can we use a direct method also? For example, in case we have a tridiagonal system (for example 1D heat ...
user avatar
  • 471
2 votes
1 answer
83 views

Is there any method to incorporate minor changes into solved meshes to speed convergence in particle-in-cell solvers?

I'm trying to perform a $10^6$ timestep electrostatic particle-in-cell simulation on a rather large mesh, with very limited computational resources (a single GPU). Because of the large number of ...
user avatar
  • 123
3 votes
1 answer
100 views

How do multigrid approaches deal with Gibbs phenomenon?

I know (from https://scicomp.stackexchange.com/a/31339/20545, among others) that I need a certain mesh density in FEM, else I might get non-physical oscillations in my solution. How do multigrid ...
user avatar
  • 543
3 votes
1 answer
225 views

Null space for smoothed aggregation algebraic multigrid

I do not really get the point of null space usage for creating the prolongation operator for smoothed aggregation algebraic multigrid. I know what the null space is per definition and I know that the ...
user avatar
7 votes
2 answers
136 views

Algebraic multigrid for complex valued matrices

Assume one uses the classical AMG with Ruge-Stuben coarsening and direct interpolation for solving real valued problems. How can this approach be recycled to also solve complex valued problems like ...
user avatar
3 votes
1 answer
388 views

Multigrid preconditioner for conjugate gradient methods

When solving $A*x=b$ using preconditioned conjugate gradient methods one has to solve $z=K^{-1}*r$ for the preconditioning where $K$ is the preconditioner of $A$ and $r$ is the residual vector. ...
user avatar
0 votes
1 answer
44 views

Converting mass density to point mass approximation on a grid

In an nbody gravity simulation, instead of doing exact(all-pair brute force) solution, I added masses of each body into cells of a 3D grid(each cell is just a float value having a mass value). Then ...
user avatar
3 votes
0 answers
568 views

Neumann-Neumann boundary intersection

I am trying to solve a Vertex Centered, Finite Difference, Poisson equation $-\nabla^{2}u=f$ with Dirichlet boundaries on the left and bottom and Neumann boundaries on the top and right using ...
user avatar
1 vote
1 answer
581 views

Algebraic multigrid in PETSc

Consider a potential/poisson equation on a very large, complicated geometry. Currently, an self-written FEM and linear solvers from NumPy are used. Performance is, of course, not good enough for ...
user avatar
  • 1,443
4 votes
3 answers
334 views

For which problems Krylov subspace methods are preferred over multigrid methods?

As multigrid methods are known to have grid independent convergence rates with $O(N)$ computational cost, then why would one be interested in using Krylov subspace methods at all, for which ...
user avatar
  • 181
1 vote
2 answers
334 views

V-cycle Multigrid for 2D transient heat transfer on a square plate using finite difference

I'm currently developing a program to solve 2D transient state heat conduction on a square plate using the V-cycle multigrid. Althought my program is able to reach the steady state solution, it's ...
user avatar
1 vote
1 answer
368 views

What is the default smoother for the "PCMG" preconditioner in PETSc?

For a large parallel sparse matrix (mpiaij type matrix) in my code, I was experimenting with various preconditioners to see which one would do best with GMRES/BiCGSTAB. I tried the PCMG ...
user avatar
  • 158
3 votes
1 answer
134 views

Can I use multigrid to solve linear algebra problems that do not arise from a differential equation?

Can the Multigrid method be applied to solve a linear system of equations in the form $Ax=b$, that is not necessarily related to differential equations?
user avatar
5 votes
0 answers
171 views

Geometric Multigrid for Conform and Non–″ Elements: Restriction Operators

First of all, let me set up some notations. Suppose we have a hierarchy of meshes $\mathcal{T}_0 \subset \mathcal{T}_1 \subset \dots \subset \mathcal{T}_n$. For simplicity I restrict myself here to $...
user avatar
  • 861
1 vote
0 answers
85 views

What does it take to prove that a multigrid algorithm scales linearly with system size?

It is my understanding that the multigrid solution techniques are generally the preferable method to solve large Poisson problems. Now assume I have written a multigrid solver that is tailored to my ...
user avatar
3 votes
0 answers
202 views

Multigrid for Robin boundary conditions

I have a question regarding the treatment of Robin boundary conditions in a multigrid solver. I am solving the Poisson equation in $\Omega=(0,1)^2$ with Robin boundary conditions on the boundary, $$- \...
user avatar
  • 101
1 vote
0 answers
244 views

Algorithm for Adaptive Mesh Refinement

I am trying to implement Adaptive Mesh Refinement. I am not a Mathematics/Computational Science person so I will try to write the algorithm in a simpler way. I will be grateful if experts can comment ...
user avatar
3 votes
2 answers
735 views

Python environments for AMG and Gauss Seidel as solvers instead of preconditioners

I am working on block preconditioning and seemingly it is common to write customised Krylov solvers for them. Within each solver, the individual block linear system with preconditioners are ...
user avatar
  • 133
5 votes
1 answer
123 views

Does algebraic multigrid reuse its coarse grids?

Admittedly, I'm new to the subject, so this is probably a really simple question. Let's assume I want to solve the (large sparse) linear system Ax = b multiple times with algebraic multigrid and b ...
user avatar
7 votes
1 answer
280 views

Robust smoothers for geometric multigrid

I'm searching for robust smoothers for geometric multigrids. By robust I mean: Effective for high order approximations (say spectral element, spectral Discontinuous Galerkin), Parallel (suitable for ...
user avatar
2 votes
1 answer
490 views

How to define residual in multigrid approach?

I wish to solve the two-dimensional Navier Stokes equations using multigrid method on a collocated grid using the Predictor-Corrector method mentioned below. But first, let me elaborate on what I had ...
user avatar
12 votes
2 answers
519 views

Algebraic Multigrid: Why does the product of interpolation and restriction not result in something with norm 1?

I'm currently working with "A Multigrid Tutorial" by Briggs et al, Chapter 8. The construction of the interpolation operator is given as: Then construction of restriction operator and fine grid ...
user avatar
  • 1,443