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

Within shared memory
Thread i − 1 Thread i Thread i + 1
Initial u i−2 , u i−1 , u i u i−1 , u i , u i+1 u i , u i+1 , u i+2
Reconstruct u L i−1 , u i−1, u R i−1 u L i , u i, u R i u L i+1 , u i+1, u R i+1
Slope limit u L∗
i−1 , u i−1, u R∗
i−1 u L∗
i , u i , u R∗
i u L∗
i+1 , u i+1, u R∗
i+1
Advance 1 2 ∆t uL i−1 , u i−1, u R i−1 u L i , u i, u R i u L i+1 , u i+1, u R i+1
syncthreads()
Solve R.P. f(u R i−2 , uL i−1 ) f(uR i−1 , uL i ) f(uR i , uL i+1 )
<strong>Finite</strong> <strong>Volume</strong> Methods
Laboratory for Scientific
Computing
10 / 22