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.
74
questions
2
votes
0answers
35 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 ...
1
vote
0answers
29 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 ...
0
votes
1answer
51 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 ...
2
votes
0answers
54 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 ...
0
votes
0answers
25 views
When semicoarsening is needed in multigrid method?
I am trying to understand multigrid in deep and all the coarsening techiques.
So,assuming a 2D grid when semicoarsening is prefered instead of standard coarsening?What is the error's behaviour when a ...
1
vote
1answer
117 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 ...
1
vote
0answers
56 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 ...
2
votes
2answers
65 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 ...
1
vote
1answer
44 views
Is there any method to incorporate minor changes into solved meshes to speed convergence in particle-in-cell solvers?
Apologies for the terrible title.
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). ...
3
votes
1answer
85 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 ...
2
votes
1answer
86 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 ...
0
votes
0answers
73 views
Algebraic multigrid for coupled equations
As far as I understand is algebraic multigrid(AMG) a method that was intentionally developed to solve linear systems where every grid point or node has a single DOF. When AMG should now be used for ...
7
votes
2answers
91 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 ...
1
vote
1answer
138 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.
...
0
votes
1answer
35 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 ...
3
votes
0answers
371 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 ...
1
vote
1answer
309 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 ...
3
votes
3answers
250 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 ...
1
vote
2answers
243 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 ...
1
vote
1answer
210 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 ...
3
votes
1answer
129 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?
4
votes
0answers
136 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 $...
1
vote
0answers
84 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 ...
3
votes
0answers
162 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,
$$- \...
1
vote
0answers
176 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 ...
3
votes
2answers
535 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 ...
5
votes
1answer
106 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 ...
5
votes
1answer
235 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 ...
2
votes
1answer
421 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 ...
12
votes
2answers
416 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 ...
2
votes
1answer
73 views
Solving new linear system that comes from an $p$ enrichment
Let's say I am solving a simple Poisson problem using a Mixed (DG) finite element method. If we use orthogonal polynomials as basis functions we can write the finite-dimensional linear system as
$$
...
0
votes
1answer
850 views
Quadtree type Grid
I would like to code for a quadtree type meshing but don't know how to do. If anyone can help or can share any starting code?
2
votes
1answer
131 views
Solve FEM matrix from coupled system
I'm developing an FEM solver for a coupled system. I have diffusion and potential equations which result in positive definite matrices for each equation, but the coupling makes the overall system ...
8
votes
1answer
338 views
Increasing V-cycles for constant Coarsest Grid Size and increasing Fine Grid size
Problem statement
I implemented geometric multigrid for $-\nabla^{2}=f$ where $f=\frac{3\pi^{2}}{4}sin \frac{\pi x}{2} sin \frac{\pi y}{2} sin \frac{\pi z}{2}$ on $\Omega \in [0,1]$ on a unit cube. ...
12
votes
3answers
269 views
Is it usual to have no convergence checking in Multigrid?
I just read Chapter 3 in "A Multigrid Tutorial" by Briggs/Henson/McCormick, link.
The text is about Multigrid cycles such as V-cycle, mu-cycle, FMG. What caught my eye: In most iterative procedures ...
9
votes
2answers
234 views
Multigrid on “not perfectly rectangular” grid
Multigrid introductions normally use a rectangular grid. Interpolation of values is then straight forward: Just interpolate linearly on the edge between two adjacent nodes of the coarse grid to find ...
1
vote
1answer
216 views
FENICS subdomains - restriction/ prolongation operators
I am trying to implement my own multigrid method in fenics. Is there any "smart/ fenics" way how to assemble subdomains and obtain restriction/ prolongation operators ?
Thanks!
3
votes
1answer
151 views
stabilizing advection-diffusion with multi-grid?
If one chooses to discetize the advection-diffusion (AD) equation using the standard Galerkin finite element method, stability issues may arise in cases of high Peclet number (i.e., high advection to ...
3
votes
1answer
153 views
Challenges in implementing Algebraic Multigrid on millions of processors
I just implemented an Algebraic Multigrid solver for a Mixed Dirichlet-Neumann Boundary Value problem and was surprised to see the speed-up as compared to a simple iterative solver for a large problem ...
1
vote
1answer
112 views
A doubt in Multigrid V-cycle
Assume I have 3 levels of grids. Finest Grid = level 2, Coarser Grid = level 1, Coarsest Grid = level 0.
Relax $u$ on $Au = b$ at level 2 for 3 times.
Find residual $r2$ at level 2, then restrict to ...
1
vote
2answers
162 views
Specific questions for 2-D Multigrid
I am simulating $\nabla^{2}u=0$ with mixed Dirichlet-Neumann boundary conditions on 2-D using 2-Grid method. Dirichlet ...
13
votes
1answer
852 views
How is Krylov-accelerated Multigrid (using MG as a preconditioner) motivated?
Multigrid (MG) may be used to solve a linear system $Ax=b$ by constructing an initial guess $x_0$ and repeating the following for $i=0,1..$ until convergence:
Compute the residual $r_i = b-Ax_i$
...
1
vote
0answers
50 views
Optimization of nonlocal stencil-like operator on subset of regular grid
I am trying to optimize the execution time for this particular piece of fortran code.
Details:
i_gc is a (ngpts, 3) array of containing (i,j,k) indices for each grid point. This is a subset of the ...
3
votes
2answers
223 views
Can F-cycle substitue FMG for update of existent solution?
I have a nicely working multigrid solver, which I use for solving the Poisson equation from an electrostatic problem. I solve this equation first without any charges, and then many times with a slowly ...
2
votes
2answers
101 views
Infinite Function Value on Dirichlet Boundary
I have been working on a multigrid solution to a non-homogeneous Dirichlet boundary value problem. However, the function goes to infinity on the boundary. This causes numerical overflow errors to be ...
0
votes
2answers
108 views
How to determine the number of c points in algebraic multi grid
I am trying to write an algebraic multi-grid solver (in c++). At a given level I determine which nodes are c-points and which nodes are f-points (where the total number of c and f points equals the ...
13
votes
3answers
3k views
Is the Thomas algorithm the fastest way to solve a symmetric diagonally dominant sparse tridiagonal linear system
I am wondering if the Thomas algorithm is the fastest way (provably?) to solve a symmetric diagonally dominate sparse tridiagonal system in terms of algorithmic complexity (not looking for ...
4
votes
1answer
247 views
Multigrid stops converging when more grid levels are used
I'm having a problem with multigrid code I wrote. If I solve Laplace's equation in 2D and use more than 5 grid levels, the V-cycles stop converging after a few cycles (see below, convergence factor > ...
0
votes
1answer
198 views
What libraries provide an implementation of multigrid?
I am working on numerical method of Multigrid. What's the available implementation(solver) (actually used in scientific computation) of multigrid method?
3
votes
1answer
353 views
Linear Algebra / Numerical Solution Of Matrix With Nullspace
I have a question relating to linear algebra.
We have a fluid solver that solves the poisson equation for pressures. When there are areas of the domain that are entirely enclosed by Neumann ...
