Skip to main content

Questions tagged [iterative-method]

A method which produces a sequence of numerical approximations which converges (provided technical conditions are satisfied) to the solution of a problem, generally through repeated applications of some procedure. Examples include Newton's method for root finding, and Jacobi iteration for matrix-vector solves.

7 votes
2 answers
490 views

I have a fairly general question, which I don't think depend a lot on my particular application. Consider a linear system $Ax = b$ that I would like to solve with preconditioned MinRes. Call the ...
FEGirl's user avatar
  • 575
1 vote
1 answer
180 views

I was numerically solving ordinary differential equations and I encountered some very interesting properties. The forward Euler method for a simple harmonic oscillator diverges, it has energy drift ...
User198's user avatar
  • 173
0 votes
0 answers
88 views

I've been working on a simple multigrid solver that uses Jacobi iterations to solve the Poisson equation as a little exercise. What I'm finding, however is that my solver doesn't seem to converge, or ...
Alberto Cattaneo's user avatar
3 votes
3 answers
221 views

For large positive definite matrices $A$, conjugate gradients is the method of choice for solving linear systems $$Ax=b.$$ Convergence of conjugate gradients heavily relies on having a good ...
Daniel Adams's user avatar
1 vote
0 answers
57 views

In Constantine A. Balanis' book about antennas, he introduced the method of moments for current distribution over a finite dipole. I found that the method of moments works very bad for a half-...
kile's user avatar
  • 121
0 votes
0 answers
60 views

I am working on evaluating the integral: $ I(\vec{r}) = \int f(k) e^{i\vec{k} \cdot \vec{r}} d\vec{k}, $ for a system with a tight-binding dispersion relation given by: $ \epsilon_k = -2t \left[\cos(...
Lawerence's user avatar
1 vote
0 answers
124 views

I have a question regarding the computational complexity of the ILU preconditioner in Python. I am trying to implement an ILU(0) preconditioner using the following code: ILUfact = sla.spilu(...
Furk's user avatar
  • 21
3 votes
1 answer
243 views

The simplest possible matrix I can think of to use an arnoldi method is the identity matrix. In this case the Krylov sequence is just $\{v, v, v, v, \cdots\}$ for any $v$. Thus the span of the krylov ...
Makogan's user avatar
  • 515
1 vote
2 answers
190 views

I am interested in analyzing convergence of the Jacobi method to solve the linear system $Ax=b$, $$\begin{pmatrix} 2 & 4 \\ 1 & 1 \end{pmatrix} \begin{pmatrix} x_1 \\ x_2 \end{pmatrix} = \...
Diffusion's user avatar
  • 141
2 votes
1 answer
181 views

In my case, I am solving $AX=B$ repeatedly, but the solution usually doesn't change much. So it'd probably be faster than me when I start from the previous solution and iterative, rather than solving ...
Taylor Fang's user avatar
1 vote
1 answer
378 views

How can I define preconditioners (SPILU, SPAI, etc.) for sparse iterative methods (TFQMR, GMRES, CGS, etc.) for the matrix-free left-hand side? I defined $Ax=b$ using matrix-free $A$ (with ...
Furk's user avatar
  • 21
4 votes
0 answers
109 views

I have been using recently AMG as preconditioner for CG with several meshes for simple elliptic problems discretised with linear elements on "complicated" three dimensional geometries and I ...
FEGirl's user avatar
  • 575
1 vote
0 answers
197 views

I am looking for an Iterative Numerical PDE solver for 1D Burgers equation. I need to have access to the intermediate solutions of the Numerical Solver. By iterative methods, I mean techniques which ...
rajoy99's user avatar
  • 11
0 votes
0 answers
164 views

Problem: I want to solve the eigenvalue problem $$x=Ax$$ to the eigenvalue $1$ for a large matrix (roughly $N^3\times N^3$ and $N$ ranges from 10 to 100) where $A$ is stochastic (i.e. all entries are ...
Diplodokus's user avatar
0 votes
0 answers
56 views

If we want to use variable preconditioning in Conjugate Gradient, we can replace the Fletcher–Reeves by the Polak–Ribière formula (https://en.wikipedia.org/wiki/Conjugate_gradient_method#...
GS101's user avatar
  • 21

15 30 50 per page
1
2 3 4 5
21