Daniel@0: function dPsi = cpGradientDiag(X, S)
Daniel@0: 
Daniel@0:     dPsi    = diag(X * S * X');
Daniel@0: 
Daniel@0: end