function r = rms(u,v)
% 
% r = rms(u,v)
% Root mean square error between vectors u and v = sqrt( \sum_i (u_i - v_i)^2 )

u = u(:); v = v(:);
r = sqrt(sum((u - v) .* (u - v) ));
