view general/numerical/matrix/issym.m @ 61:eff6bddf82e3 tip

Finally implemented perceptual brightness thing.

function z=issym(x,th)
% issym - Test if matrix is symmetric with optional tolerance
%
% issym :: [[N,N]] -> bool.
% issym :: [[N,N]], nonneg -> bool.

if nargin<2,
	z=all(all(x==x'));
else
	z=all(all(abs(x-x')<=th));
end