1) It takes a long time to display densities and velocities since you pressed the 'v' keyboard button
2) The default size of rectangular matrix in jos stam paper was 64, in order for my simulation to run i have reduce it to 16, as soon as i try to increase to N = 64 (default in jos stam example) the simulation blows up.
If anyone has any idea as to what i am doing wrong please let me know.
If you can redirect me to the correct solution i would be more than appreciated.
I have attached Jos Stam code and also my modifications on Jos Stam's code. If you would like please take a look to have a concrete understanding on their differences.