Hi everyone, i was trying to implement jos stam paper for fluid dynamics from 2D in 3D but i think that along the way I have messed things up. Well when the simulation runs if you hit the 'v' (toggle density and velocity) on your keyboard you will see a cube (3D) instead of a rectangle (2D) but the problems that arise are :

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.