deltavpn buffer hopefully should be computed right now. The matrix is split into 3 vectors, so we can make use of the atomicAdd function. Also a bugfix which prevented writing into the velBn buffer. You can now see that the particles actually accelerate instead of falling at a constant speed. Next thing would be to implement collisions on gridnodes and particles.