TY - JOUR
T1 - A fast linear complementarity problem solver for fluid animation using high level algebra interfaces for GPU libraries
AU - Andersen, Michael
AU - Abel, Sarah Maria Niebe
AU - Erleben, Kenny
PY - 2017/12
Y1 - 2017/12
N2 - We address the task of computing solutions for a separating solid wall boundary condition model. We present a parallel, easy to implement, fluid linear complementarity problem solver. All that is needed is the implementation of linear operators, using an existing high-level sparse algebra GPU library. No low-level GPU programming is necessary. This means we can rely on the efficiency of a tried-and-tested library, requiring minimal debugging compared to writing more low level GPU kernels. The solver exploits matrix-vector products as computational building blocks. We block the matrix-vector products in a way that allows us to evaluate the products, without having to assemble the full systems. Our work shows speedup factors ranging up to two orders of magnitudes for larger grid resolutions.
AB - We address the task of computing solutions for a separating solid wall boundary condition model. We present a parallel, easy to implement, fluid linear complementarity problem solver. All that is needed is the implementation of linear operators, using an existing high-level sparse algebra GPU library. No low-level GPU programming is necessary. This means we can rely on the efficiency of a tried-and-tested library, requiring minimal debugging compared to writing more low level GPU kernels. The solver exploits matrix-vector products as computational building blocks. We block the matrix-vector products in a way that allows us to evaluate the products, without having to assemble the full systems. Our work shows speedup factors ranging up to two orders of magnitudes for larger grid resolutions.
U2 - 10.1016/j.cag.2017.09.006
DO - 10.1016/j.cag.2017.09.006
M3 - Journal article
SN - 0097-8493
VL - 69
SP - 36
EP - 48
JO - Computers and Graphics
JF - Computers and Graphics
ER -