Parallelize virial radius calculation and spherical filters. Add unit test for underlying code.