ENH: Stream & density field openmp parallelization

This commit adds basic threading of the multi-stream and
density field calculation using openmp. This yields a
~2x speedup for this step, but performance can be improved
further with a bit more work.