Implementing Finite Volume algorithms on GPUs - many-core.group ...

Template over blockDim
Initially allocated shared memory dynamically
Drawback:
Offsets only calculated at run-time
Overflow of shared memory only found at run-time
So template all functions over blockDim.x and blockDim.y:
template
global void getSLICflux GPU(Grid u, Grid flux, float dt, int
coord)
{
shared float temp[4][NUM VARS][blockDim y][blockDim x];
}
Compiler can pre-compute memory offsets
Use of excess shared-memory can be detected at compile-time
Drawback: Must determine block-size at compile-time
Time now: 8.23s (109.7×)
<strong>Finite</strong> <strong>Volume</strong> Methods
Laboratory for Scientific
Computing
13 / 22