faster disp cos/sin calculation.
(hoping even a dump compiler will figure the cos/sin pair, and use
vector operation for the multiply)
Put particles to mesh center.

Otherwise the top-right side (positive side) of the walls will
have more denisty than the bottom left side of the walls, for
the CM drifting of different masses are different.