test-chill/test-cases/examples/chill/jacobi1.c


1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19


#ifndef N
#define N 512
#endif

/*
<test name=jacobi define="{'N':512}">
procedure int jacobi(
    in out float[N][N] a = matrix [i,j] random(2,-2))
</test>
*/
int jacobi(float a[N][N]) {
    int t, i;
	for (t = 2; t <= 100; t++)
		for (i = 2; i <= N - 1; i++)
			a[t][i] = a[t - 1][i - 1] + a[t - 1][i] + a[t - 1][i + 1];

	return 0;
}