93
edits
Changes
→Assignment 3
completion.
</pre>
In our kernel, our main jacobi calculations look like this (when using global indexing):
<source>
a[ij] = (input + dxxinv*(b[ij - 1] + b[ij + 1])
+ dyyinv*(b[ij - m] + b[ij + m]))*dcent;
</source>
* b[ij - 1] is a cell's left neighbour
* b[ij + 1] is a cell's right neighbour
* b[ij - m] is a cell's top neighbour (we represent a 2D array as a 1D array, so subtract row length to go one cell up)
* b[ij + m] is a cell's bottom neighbour (add row to go one cell down)
[[File:Dps_915_gpusquad_needneighbour.png]]
[[File:Dps_915_gpusquad_ghostcellexample.png]]
