Changed gsm_noise_covar() so that it exploits matrix symmetry for speed.

Added a new file mpi_wrapper.py that is a template for using an mpi master/slave set up to compute the entries of a matrix.