Last fix to computing C matrix, I promise. Ready for parallelI/O