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