function y=rectify(x)
% rectify - 'half wave' rectification, rectify(x) x>0 ? x : 0
%
% rectify :: [Size->real] -> [Size->nonneg].

y = (x>0) .* x;
% alternative implementation?
% i=(x<0); y=x; y(i)=0;
