Although the focus of the paper is on #SmoothedParticleHydrodynamics, and the simulation of viscous fluids, the key finding is of more general interest: the numerical stability of #BiCGSTAB can be improved by rewriting the standard formulation https://en.wikipedia.org/wiki/Biconjugate_gradient_stabilized_method to avoid catastrophic cancellations in the computation of some important coefficients. Anyone interested in solving large linear systems would benefit from adopting the proposed alternative form of the method.

I'm still preparing my post series on #GPGPU, but in the mean time you can read some about in our most recent published paper:
“A numerically robust, parallel-friendly variant of #BiCGSTAB for the semi-implicit integration of the viscous term in #SmoothedParticleHydrodynamics”, freely accessible for the next 50 days from this link:
https://authors.elsevier.com/a/1fN8C508HsZiG