# Questions tagged [preconditioning]

For questions regarding design and implementation of preconditioners for solving linear systems.

13 questions
Filter by
Sorted by
Tagged with
11k views

### Why is my iterative linear solver not converging?

What can go wrong when using preconditoned Krylov methods from KSP (PETSc's linear solver package) to solve a sparse linear system such as those obtained by discretizing and linearizing partial ...
365 views

### How to directly compute the inverse of an ill-conditioned dense matrix

I know that it is generally a bad idea to compute the inverse matrix directly. However, if it is necessary to compute the inverse of an ill-conditioned invertible dense matrix, then what can I try? ...
478 views

### role of initial guess for iterative linear solver

Suppose we use a preconditioned iterative solver for a linear system. If the initial state for the solver can be chosen very close to the exact solution - does this reduce requirements for the ...
145 views

Suppose that I have a large (on the order of 10^6 unknowns) 3D scalar Poisson system which I would like to precondition. The boundary conditions have been treated so that the system is SPD. I.e., \... 5answers 3k views ### What is the advantage of multigrid over domain decomposition preconditioners, and vice versa? This is mostly aimed for elliptic PDEs over convex domains, so that I can get a good overview of the two methods. 2answers 2k views ### Which preconditioners (and solver) in PETSc for indefinite symmetric systems should I use? My system is a symmetric FE problem with lagrange multipliers (e.g. incompressible Stokes' flow): \begin{pmatrix}A & B^T \\ B & C\end{pmatrix} where C = 0 is the typical case (I have even ... 2answers 281 views ### Why does conjugate gradient work with this nonsymmetric preconditioner? In this previous thread the following multiplicative way to combine symmetric preconditioners P_1 and P_2 for the symmetric system Ax=b was suggested: \begin{align} P_\text{combo}^{-1} :=& ... 2answers 175 views ### Why iterative method: AMG preconditioned PCG is slower than Matlab direct method 'A\b'? Recently, I have met a question that a saying goes that for large linear system: iterative methods are required because of memory problem of direct methods. But when I implement some experiments ... 2answers 422 views ### Is there any way to do “double preconditioning” Question: Suppose that you have two different (factored) preconditioners for a symmetric positive definite matrix A:A \approx B^TB$$and$$A \approx C^TC,$$where the inverses of the factors B, ... 2answers 212 views ### Preconditioner for the GMRES method in the Uzawa algorithm I'm trying to solve \begin{equation}\left\{ \begin{split} \frac{\partial u}{\partial t}+(u\cdot\nabla)u-\nu\Delta u+\frac1\rho\nabla p&=f\;\;\;\text{in }\Lambda\\ u&=0\;\;\;\text{on }\partial\... 1answer 494 views ### Which preconditioning for large linear elasticity problem? The problem I want to solve is the displacement formulation of the linear elasticity :$$ \nabla \cdot \sigma = 0 \quad \text{in} \quad \Omega \\ \sigma = \lambda ( \nabla \cdot u ) I + \mu (\nabla \...
I am a physicist with limited numerical methods knowledge and I am trying to speed up the inversion of a very ill-conditioned problem ($rcond>10^{30}$). The same sparse square matrix is used ...