modified source of exp_stationary_sample_and_diagonal_approximation to remove dependence on make_stat_cov_matrix and to factor out the square root of the matrix