slides - Friedrich-Schiller-Universität Jena

abstract parallelization
array x0(0,n), x1(0,n);
for int m = 1 to c {
loop tiling
for_parallel int i0 = 0 to n step n0 {
array y0(i0-1,i0+n0), y1(i0,i0+n0-1);
y0 = x0;
for_parallel int i = i0 to i0+n0-1 {
y1(i) = if (i==0) 0.0
else if (i==n) 0.0
else 1.0 – y0(i) + 0.5 * ( y0(i-1)+y0(i+1) );
}
x1 = y1;
}
swap(x1, x0);
}
local array y0, y1
data transfer at
y0=x0, x1=y1