root/gliderproc/trunk/MATLAB/matutil/identical.m

function theResult = identical(varargin)

% identical -- Nan-safe "isequal".
%  identical(x, y, ...) returns TRUE if all the
%   arguments are identical in shape and value,
%   even with respect to NaNs.  No distinction
%   is made between the number-like classes.
 
% Copyright (C) 2001 Dr. Charles R. Denham, ZYDECO.
%  All Rights Reserved.
%   Disclosure without explicit written consent from the
%    copyright owner does not constitute publication.
 
% Version of 05-Sep-2001 01:48:57.
% Updated    05-Sep-2001 02:32:14.

if nargout > 0, theResult = []; end
if nargin < 2, help(mfilename), return, end

% Check "isequal" first.

yes = isequal(varargin{:});
if yes
        if nargout > 0
                theResult = yes;
        else
                disp(yes)
        end
        return
end

% Not exactly equal.  Decompose the arguments,
%  checking for equality of the pieces.

x = varargin{1};

for k = 2:length(varargin)
        y = varargin{k};
        yes = isequal(x, y);
        if ~yes
                switch class(x)
                case 'cell'
                        if isequal(class(y), 'cell')
                                if isequal(size(x), size(y))
                                        for i = 1:prod(size(x))
                                                yes = feval(mfilename, x{i}, y{i});
                                                if ~yes, break, end
                                        end
                                end
                        end
                case 'struct'
                        if isequal(class(y), 'struct')
                                fx = fieldnames(x);
                                if isequal(fx, fieldnames(y))
                                        for i = 1:length(fx)
                                                gx = getfield(x, fx{i});
                                                gy = getfield(y, fx{i});
                                                yes = feval(mfilename, gx, gy);
                                                if yes, break, end
                                        end
                                end
                        end
                case {'char', 'double', 'sparse', 'uint8', 'uint32'}
                        nx = isnan(x);
                        if isequal(size(x), size(y)) & any(nx(:))
                                ny = isnan(y);
                                if isequal(nx, ny)
                                        yes = isequal(x(~nx), y(~ny));
                                end
                        end
                otherwise
                end
        end
        if ~yes, break, end
end

if nargout > 0
        theResult = yes;
else
        disp(yes)
end