Methods to Load Balance a GCR Pressure Solver Using a Stencil Framework on Multi- and Many-Core Architectures

Algorithm 4

 do 3 k=2,lp-1

 do 3 j=2,mp-1

 do 3 i=2,np-1

3 out_p(i,j,k)=cCoeff(

 ⋅ in_p(i+1,j,k)+in_p(i-1,j,k)+

 ⋅ in_p(i,j+1,k)+in_p(i,j-1,k)+

 ⋅ in_p(i,j,k+1)+in_p(i,j,k-1))