Daniel@0: function r = regularizeTraceDiag(W, X, gradient)
Daniel@0: %
Daniel@0: % r = regularizeTraceDiag(W, X, gradient)
Daniel@0: %
Daniel@0: %
Daniel@0:     d = length(W);
Daniel@0: 
Daniel@0:     if gradient
Daniel@0:         r = ones(d,1) / d;
Daniel@0:     else
Daniel@0:         r = sum(W) / d;
Daniel@0:     end
Daniel@0: end