Daniel@0: function r = regularizeTwoDiag(W, X, gradient)
Daniel@0: %
Daniel@0: % r = regularizeTwoDiag(W, X, gradient)
Daniel@0: %
Daniel@0: %
Daniel@0:     d = length(W);
Daniel@0: 
Daniel@0:     if gradient
Daniel@0:         r = W / d;
Daniel@0:     else
Daniel@0:         r = sum(sum(W.^2)) / (2 * d);
Daniel@0:     end
Daniel@0: end