This repo contains an implementation of Incomplete-Cholesky preconditioned conjugate gradient algorithm using c++ and cuBLAS/cuSPARSE (CUDA 11.0) which I used to make a 2D semi-lagrangain fluid simulatoin.
You can find the fluid simulation here.
See this tutorial written in Chinese for more implementation details.
The algorithm I used:
