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