Changes

To avoid data-racing, in each iteration m, we stored epsim12 into epsim_cell, and compute f using another for loop. Using parfor on a 12-cores machine gives a four time speed up for computing msmf_corr_coeff.m (apology for not having a profiling screenshot here). Note that to avoid data-racing, in each iteration m, we stored epsim12 into epsim_cell, and compute f the actual speedup depends how many cores you are using another for loop.