Compute variance of the gradient estimator

[riberaborrell]

at each SGD iteration sample the gradient multiple times (10-50) and compute variance.