A few simple unit tests for matrices
Particle constructor with omega
Remove p.r < r_max assert
Scalar multiplication for matrices