function x = l2(A)
% ell-2 operator norm of A
% given by the maximal eigenvalue^.5 of A'*A
x = sqrt(max(eig(A'*A)));
return