Hey!
I'm implementing something similar following Bridson's textbook as well! Sometimes I look at this codebase as a reference ;)
Wondering however whether or not this is a mistake on my end or not. It seems the A Matrix as constructed in your pressureSolve() function has coefficients that are not multiplied by the timestep length, as the textbook suggests they should be.
Maybe you just assume a timestep of 1 and that works out for you?
Thanks