Changed the way of smearing the covariance matrix: instead of using the smearing functions, just scale the matrix elements with (1+ SF^2) (as suggested by @tstreble)
@kleney please try this one