Skip to main content
The 2024 Developer Survey results are live! See the results

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.

-1 votes
1 answer
69 views

optimal gradient algorithm to determine best $α_k$

Let's consider an optimal-step gradient algorithm and assume that: $g(α) := f(X_k - α∇f(X_k)) = 2α^2-4α+17$, how can we determine the optimal $α_k$? Here is my simple solution: $g(α) = 2α^2-4α+17$ $g'(...
V_head's user avatar
  • 15
0 votes
1 answer
45 views

step-fixed algorithm first iterates

let us have the fixed-step gradient algorithm, with $p = 2$ and we assume that for $X = (x, y)$, $∇ f(X) = \begin{pmatrix} x -1\\ y -2 \end{pmatrix}$ Let me assume we intialize with $X_0 = (0,0)$ what ...
V_head's user avatar
  • 15
0 votes
1 answer
66 views

step-fixed algorithm to minimize f, which step to ensure convergence?

If we want to apply the fixed-step gradient algorithm to the minimization of $f(x) = \frac{1}{2}(Ax, x)$ where $A$ is a symmetric 2x2 matrix with eigenvalues $\lambda_1 > \lambda_2 > 0$, for ...
V_head's user avatar
  • 15
0 votes
1 answer
39 views

Estimating the rate of convergence of Projected Gradient Descent on constrained polynomial objectives

I am estimating the order of convergence of Projected Gradient Descent (GD) on quadratic polynomials with random coefficients independently drawn from Uniform(-1,1) and bounded by a unit hypercube. I'...
ufghd34's user avatar
  • 23
1 vote
1 answer
127 views

How to refine $h$ and $\Delta t$ for convergence tests on evolution PDE

Setting I am solving for $u(x,y,t)$ the wave equation $\partial_{tt} u - \partial_{xx} u = f$ on $(x,y)\in\Omega=[0,1]\times \mathbb{R}$ by splitting it into an equivalent first order system: $$\...
user1313292's user avatar
0 votes
0 answers
110 views

Is there a fast matrix-free inverse power iteration?

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
51 views

Calculating a 2D Ewald sum for a multipolar expansion

I am attempting to calculate the potential of a particle at the center of an infinite two-dimensional lattice as per the following reference: Reference: Lambin, PH & Senet, P. Ewald Summation of ...
JasonC's user avatar
  • 43
3 votes
1 answer
187 views

Stability of Euler forward method

I am trying to solve a linear system of ODEs of the form: $$ \frac{du}{dt} = A u, \quad u(0)=k$$ where $A$ is a 2x2 matrix and $u(t)$ is a 2x1 column vector. I want to solve this numerically, using ...
rainbow's user avatar
  • 31
3 votes
1 answer
150 views

Role of rotation's pivot point in optimization?

In this paper, the authors describe how to use locally rigid transformations (sampled on nodes in space) to deform mesh vertices. In the paper, rotations are relative to the pivot point, which ...
jordi's user avatar
  • 31