Jacobi Iteration, The Loop

config var n : integer = 512;

eps : float = 0.00001;

region R = [1..n, 1..n];

var A, Temp : [R] float;

err : float;

direction N = [-1, 0]; S = [ 1, 0];

E = [ 0, 1]; W = [ 0,-1];

procedure Jacobi();

[R] begin

A := 0.0;

[N of R] A := 0.0; [W of R] A := 0.0;

[E of R] A := 0.0; [S of R] A := 1.0;

repeat

Temp := (A@N + A@E + A@W + A@S)/4.0;

err := max<< abs(Temp - A);

A := Temp;

until err < eps;

end;

end;

Previous slide Next slide Back to first slide View graphic version