Daniel@0: function dPsi = cpGradientFullMKL(X, S)
Daniel@0: 
Daniel@0:     [d,n,m] = size(X);
Daniel@0: 
Daniel@0:     dPsi = zeros(d,d,m);
Daniel@0:     for i = 1:m
Daniel@0:         dPsi(:,:,i)    = X(:,:,i) * S * X(:,:,i)';
Daniel@0:     end
Daniel@0: 
Daniel@0: end