wolffd@0: function dPsi = cpGradientFull(X, S, batchSize)
wolffd@0: 
wolffd@0:     dPsi    = X * S * X' / batchSize;
wolffd@0: 
wolffd@0: end