WebGL-PIC-FLIP-Fluid
WebGL-PIC-FLIP-Fluid copied to clipboard
According to the original paper, the Incomplete Poisson preconditioner is slightly different
The Incomplete Poisson preconditioner according to the original paper is M = (I - LD^-1)(I - D^-1L^T) i.e. it uses the transpose of L in the second term.
It would be interesting to see whether that makes a performance difference compared to your implementation M = (I - LD^-1)(I - D^-1L)?