Concurrent multiprocessing methods for calculating nullspace and range space bases for multibody simulation