Questions tagged [convergence]

Questions related to whether the sequence of iterates generated by an iterative method has one or more limit points, and if those limit points have the correct properties.

Filter by
Sorted by
Tagged with
2
votes
0answers
29 views

Convergent Finite Difference Scheme for Parabolic Equation

Consider the PDE $$u_t = b_{11}u_{xx} + 2b_{12}u_{xy} + b_{22}u_{yy},$$ where $b_{11}, b_{22} > 0$, and $b_{12}^2 < b_1b_2$. In Strikwerda's book, the ADI scehme \begin{align*} \left( 1 - \frac{...
5
votes
1answer
91 views

Non-conforming bi-linear finite element

The four-noded bi-linear rectangle element, which sometimes goes under the name Melosh element, is non-conforming unless the element sides are aligned. Out of curiosity I have implemented this element ...
2
votes
1answer
72 views

Best way to check if SOR solution has converged for 2d matrix

I have written a SOR algorithm to solve the Laplace equation on a 2d grid. The outside of the grid is fixed at 0 and the central square is fixed at 10. I can obtain the fully converged solution for ...
3
votes
1answer
91 views

Determine stability of an algorithm?

This is related to a question I answered on Stack Overflow regarding calculating the square root of a number. I was thinking about it and realized that the formula is just the first in a family of ...
2
votes
4answers
128 views

Testing the time dependent Schrodinger Equation with an analytical solution?

I am numerically solving the Schrodinger Equation in 1D first and in higher dimension later, but I want to know the convergence rate of my numerical solver in different grid size and numerical methods....
1
vote
1answer
68 views

Unexpectedly Slow Convergence Implicit Euler

I'm solving the coupled ODE $$ \left[\begin{array}{c}x^\prime(z)\\p_x^\prime(z)\end{array}\right] = C(z)\cdot\left[\begin{array}{c}x(z)\\p_x(z)\end{array}\right] = \left[\begin{array}{cc}0& A(z)\\...
1
vote
1answer
58 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?
0
votes
0answers
61 views

A name for a numerical phenomena when using numerical methods

I have a nonlinear solver for equation $g= c_1f(x_1,y_1)+c_2f(x_2,y_2)$. Note that $c_1$ is much bigger than $c_2$. So after using Levenberg–Marquardt algorithm, I could only get $x_1$, $y_1$ and $...
2
votes
2answers
83 views

Implementing Gelfand’s formula for the spectral radius in Python - lack of convergence

For context: Gelfand's formula for the spectral radius is $\lim_{k\rightarrow \infty}|A^k|^{1/k}$ where $|\cdot|$ is any well-defined operator norm. I naively coded a function to calculate the $k$th ...
3
votes
0answers
51 views

Why GA convergence curves continue as two parallel lines?

I'm working on a optimization problem and using GA algorithm (in MATLAB, ga function). As you know MATLAB plots GA result with two curves, one for the best values and other to show the mean values ...
11
votes
1answer
182 views

Implications of thermodynamic inconsistency in CFD calculations

During my PhD work, I had to use tabulated values of thermodynamic properties of gases in some Computational Fluid Dynamics (CFD in short) simulations. My tables are discretized in temperature and ...
2
votes
3answers
121 views

How the gmres method iteration behaves for this **enfant terrible** matrix?

Recently, I have been studied my lessons about gmres iteration, probably the most popular iteration method for general large sparse linear system of equations Ax=b. And the convergence is obtained ...
0
votes
0answers
38 views

Divergence issues when using intrinsic cohesive elements approach

When I model the strain localisation of a microscopic sample (or say RVE ) with cohesive elements approach, the convergence performance looks very terrible. I have to use extremely time increments (...
3
votes
1answer
73 views

Conjugate Gradient for singular 2D poisson finite element with Neumann Boundary Conditions

Heavily edited question after I realised partly what the problem was I have programmed a simple 2D square finite element solution to the Poisson equation $-\Delta u = f$ The source function ...
2
votes
1answer
108 views

Number of GMRES iterations increase when stepping forward in time, using the Newton method

I am solving a system of nonlinear time-dependent equations using the Newton method in a finite-element-setting, i.e. first I create the jacobian matrix for the current time, and afterwards I try to ...
0
votes
0answers
23 views

Question about the visible and hidden neurons in neural networks methods

My problem is the following : I found the ground state energy (for the Ising model) with neural networks (more specifically RBM). I reproduced the same result but by increasing every time the ratio $=...
0
votes
0answers
32 views

Solutions in Ngspice do not converge to right value

I recently noticed while using Ngspice that the solution converges to a deviated value if the number of iterations per unit period for a square wave is less than 100, i.e. if I am using a square wave ...
6
votes
0answers
61 views

Why not use the preconditioned residual as termination criterion for preconditioned CG?

I have a Poisson equation with wildly varying material parameters (1 .. 1000), wildly varying element sizes (5 nm .. 100 um) and some quite anisotropic (tetrahedral) elements (100 nm x 100 um). I use (...
0
votes
0answers
49 views

Why does GMRES converge much slower for large Dirichlet boundary conditions?

I'm trying to numerically solve a simple Laplace equation in 2D, with a nonlinear source term: $\nabla^2 u = u^2$ with boundary conditions as $u=0$ everywhere except for $y=1$ where $u=u_0$. I'm ...
8
votes
1answer
185 views

Should we always expect FEM error plots to be straight lines?

The error estimates in FEM are usually of the form $$||u^h-u||\leq Ch.$$ Taking logarithm on both sides, we obtain $$\log ||u^h-u||\leq \log C + \log h.$$ This estimate implies that the error lies ...
2
votes
0answers
26 views

Difference between contraction rate and convergence rate

I am slightly confused about the concept of contraction rate. For me it sounds equivalent to a convergence rate. Could somebody clear up the difference for me? I have those two definitions of ...
1
vote
1answer
89 views

Abaqus, ANSYS, and FVM solver for thermal expansion problem converges to different values

Is it reasonable for a FEM and FVM code to converge to slightly different solutions for the same physical problem (identical BCs, geometry, properties, etc...), provided stability constraints are ...
2
votes
1answer
66 views

Blowup of error in Conjugate Gradient method with periodic Dirichlet Poisson matrix

My problem is that the L2-Norm of the residual for the periodic Poisson matrix $P$ is initially decreasing but starts to blow up after a certain number of iterations. The blowup happens earlier the ...
0
votes
0answers
29 views

Fit exponential convergence

I'm working with a numerical algorithm whose output $y$ asymptotically approaches a certain unknown value $a$. I expect an exponential convergence, i.e. the data $y$ given by my algorithm should be ...
2
votes
0answers
44 views

How to determine the order of convergence of the Euler-Maruyama method?

This question is originally posted in Quant.StackExchange but has been unanswered for some time so I ask in here. To make this simple let us consider the Geometric Brownian Motions (GBM). My ...
1
vote
1answer
50 views

Simplest way to precondition Uzawa iteration

I have a diffusion problem with an internal circular dirichlet constraint and a side condition which shall enforce a certain global volume integral. $\nabla(D \nabla u(x)) = 0$ outer boundary ...
7
votes
1answer
232 views

Lack of quadratic convergence in Newton's method

It is well-known that Newton's method can converge quadratically, if initial guess is close enough and if the arising linear systems are solved accurately. I am applying Newton's method to highly ill-...
0
votes
0answers
45 views

Correct order of convergence

I have a sequence of points which was obtained from an iterative algorithm, and I computed the order of convergence $p$ of the method using the formula $$ p \approx \frac{\log({\rm err}(k+2))-\log({\...
3
votes
3answers
204 views

$H^1$-convergence rate of finite element method for Poisson equation, depending on element order

I wanted to verify my FEM-program by applying the method of manufactured solutions, while solving the Poisson equation in two dimensions using the continuous Galerkin method $$-\nabla^2u=f$$ with $$u=...
4
votes
2answers
101 views

Analytical convergent sequence and numerical divergent sequence

Is it possible to construct a sequence that converges in theory but when computed numerically with a computer program is diverging. I feel that today our computer programs doesn't allow such ...
-1
votes
1answer
92 views

Finite element convergence rate and possion's ratio

I am running simulations of a cantilever beam where it is fixed on one end and negative force applied to the other end. The first simulation is with 4-node linear quadrilateral elements and the other ...
3
votes
2answers
173 views

Asymptotic error of forward Euler

I'm trying to understand the asymptotic error behaviour of forward Euler (finite difference method), as timesteps are decreased (refined), so I feel trust in the method of manufactured solutions (MMS)....
1
vote
1answer
113 views

Analytical testcase for 2D/3D anisotropic Diffusion (Heat Kernel)

I want to verify and compare different Discretizations of the anisotropic diffusion equation in 2D / 3D. In order to both test the timestepping and the spatial discretisations I had a look at using ...
3
votes
0answers
139 views

Unstable convergence of a Poisson equation

What could be the reason that the solution of a Poisson equation is smooth when obtained by an iterative solver, only if the maximum residual is set to a high value (e.g. 0.1)? When the maximum ...
4
votes
0answers
80 views

Preconditioned residual converges, but true residual doesn't

I'm using Albany w/ Trilinos to solve an elasticity problem with thermal expansion mismatch. I'm using block GMRES with MueLu preconditioning. It works for problem size of several million dofs, but ...
4
votes
1answer
157 views

Convergence rate Jacobi/Gauss-Seidel with mesh resolution

In the book A Multigrid Tutorial - Briggs, Henson. McCormick in the beginning of Chapter 3, it is mentioned that ...because the convergence factor behaves as 1-$O(h^{2})$, the coarse grid ...
0
votes
1answer
148 views

Can I convert CUDA core to CPU core and use it as cpu core while running any program?

I was using Metatrader5 and have designed a strategy for trading using MQL5 programming language. While I was running a Strategy Optimization process, I saw the it will need 10,00= iterations or ...
0
votes
1answer
48 views

Conflicting definition of limit point

This question was raised at a different place without sufficient answers. Definition 1: We say that a vector $x \in R^n$ is a limit point of a sequence $\{x_k\}$ in $R^n$ if there exists a ...
1
vote
0answers
51 views

Decrease in slope during convergence analysis

I am using the method of manufactured solutions to perform the order of accuracy testing. I am using a cube for the testing. The cube is size 1m on all sides. I used 5 refinements: $dx = dy = dz = ...
3
votes
2answers
675 views

Numerical Sensitivity in Density of States of Tight-binding model

I'm working with the tight-binding model, and I'm trying to learn the basics of how to compute the Density of States (DOS) $N(E)$ numerically. The DOS is given by $$N(E) = \frac{1}{N}\sum_k \delta(...
3
votes
1answer
180 views

Numerical Lax-Wendroff scheme order of convergence on Burgers equation

I was suggested to move that question here. The question to be as follows. Statement of the problem Is it possible to achieve the second order of convergence (OOC) of Lax-Wendroff (LxW) scheme ...
0
votes
1answer
75 views

How to report non-monotonic runtimes in convergence plots

Let's say I have an algorithm that can be tuned by a parameter $h>0$ and is expected to converge as $h\to 0$. I want to study the computational complexity of this algorithm, i.e., how the ...
1
vote
0answers
59 views

BFGS convergence problem

I would like first to state that it is beyond my capability to identify whether this is a BFGS issue or a R package problem. I've been doing some mixed logit regression using the R package "mlogit". ...
1
vote
0answers
75 views

Successive iteration method for solving eigenvalue ploblem

I have a question concerning the branch of successive iteration methods (Newton, Runge-Kutta). I definitely know (or can read in Wikipedia) the implementation of these methods. But I was wondering ...
3
votes
1answer
405 views

Scaling/nondimensionalization for numerical optimization

I have a numerical optimization problem that I am trying to scale appropriately, in order to allow for the solver to achieve faster and more accurate results. I found a paper here that had a short ...
2
votes
2answers
660 views

Convergence problem for Poisson equation with periodic BC

I have written Poisson solvers using two different methods: A classic Jacobi scheme and one using the multigrid solver Hypre. I made up a couple of test cases ensuring the validity of those solvers. ...
1
vote
2answers
112 views

First-order ODE scheme implementation giving less than first-order convergence?

I am solving the initial value problem $$ \frac{d}{dt} (E C_g) = -\delta, \quad E(0) = E_0, $$ for $E$, where $E$ and $C_g$ are functions of $t$, $C_g$ is completely known, and $\delta$ is a function ...
1
vote
0answers
169 views

Accuracy of finite difference method for heat equation on a disk

To study an approximation for the heat equation $$\frac{\partial^2 u}{\partial r^2}+\frac{1}{r}\frac{\partial u}{\partial r}+\frac{1}{r^2}\frac{\partial^2 u}{\partial\theta^2}=f(r,\theta)$$ on the ...
4
votes
0answers
396 views

Convergence rate of Picard iterations

Given a first order ODE $y'(x)=f(x,y)$ with the initial condition $y(x_0)=x_0$ such that it satisfies Picard thoerem of existence and uniquness, one can compute the solution by Picard iterations : $$ ...
0
votes
1answer
346 views

Line search bracketing for proximal gradient. Is it good idea?

Maybe my question is obvious but i cannot find any good source which answers it I trying to learn about proximal gradient. One thing which is not clear for me is particular algorithm for line search. ...