slides - Friedrich-Schiller-Universität Jena

translation to array style
add types, declarations
recursive to iterative
detect loop parallelism
array x0(0,n), x1(0,n);
for int m = 1 to c {
for_parallel int i = 0 to n {
x1(i) = if (i==0) 0.0
else if (i==n) 0.0
else 1.0 – x0(i) + 0.5 * ( x0(i-1)+x0(i+1) );
}
swap(x1, x0);
}