Base cpu/gpu implementation

Naive n^2 neighbor search, viscosity, pressure, but no vorticity