n
it
ω
This GPU based script draws 10 cross-sections 2it weighted Jacobi iterations.
f.
_{i,j} = sin(πi/n) sin(πj/n) |

On square

Weighted Jacobi iterations are

Transformation (*) has eigenvectors (we put

with eigenvalues

λ

From

At last for small

so relaxation is very slow for long wave harmonics. As since convergence is faster for small

Simulations on GPU