utils/matlab/stats.m

function [varargout] = stats(A)
% stats(A) writes the basic statistics of the A to the terminal which
% are i) the minimum finite value
%    ii) the maximum finite value
%   iii) the mean of the finite values
%    iv) the S.D. of the finite values ( RMS of [A-mean] )
%     v) the fraction of non-finite elements excluded from calculations
%
% e.g.
%  >> stats(topo)
% Min -4555     Max 0         Mean -2331.07  SD 1207.3441   N-Z 1024/1024
%
% [Min Max Mean SD]=stats(topo); returns the statistics in Min, Max,
% Mean and SD and does not write to the terminal.
%
% $Header: /u/gcmpack/models/MITgcmUV/utils/matlab/Attic/stats.m,v 1.1.2.1 2001/04/13 18:28:23 adcroft Exp $

%ii=find(A~=0);
ii=find(isfinite(A));
if isempty(ii)
 ii=1;
end
sZ=prod(size(A));
minA=min(A(ii));
maxA=max(A(ii));
meanA=mean(A(ii));
sdA=sqrt( mean( (A(ii)-meanA).^2 ) );
nZ=prod(sum(sum(sum(~isfinite(A)))));
switch max(nargout)
 case {0}
  disp( ...
   sprintf('Min %0.5g  Max %0.5g  Mean %0.5g  SD %0.5g   NaN %i/%i',...
         minA,maxA,meanA,sdA,nZ,sZ) );
 case {1}
  varargout(1)={minA};
 case {2}
  varargout(1)={minA};
  varargout(2)={maxA};
 case {3}
  varargout(1)={minA};
  varargout(2)={maxA};
  varargout(3)={meanA};
 case {4}
  varargout(1)={minA};
  varargout(2)={maxA};
  varargout(3)={meanA};
  varargout(4)={sdA};
 otherwise
  error('Too many return arguments requested')
end
1	function [varargout] = stats(A)
2	% stats(A) writes the basic statistics of the A to the terminal which
3	% are i) the minimum finite value
4	% ii) the maximum finite value
5	% iii) the mean of the finite values
6	% iv) the S.D. of the finite values ( RMS of [A-mean] )
7	% v) the fraction of non-finite elements excluded from calculations
8	%
9	% e.g.
10	% >> stats(topo)
11	% Min -4555 Max 0 Mean -2331.07 SD 1207.3441 N-Z 1024/1024
12	%
13	% [Min Max Mean SD]=stats(topo); returns the statistics in Min, Max,
14	% Mean and SD and does not write to the terminal.
15	%
16	% $Header: /u/gcmpack/models/MITgcmUV/utils/matlab/Attic/stats.m,v 1.1.2.1 2001/04/13 18:28:23 adcroft Exp $
17
18	%ii=find(A~=0);
19	ii=find(isfinite(A));
20	if isempty(ii)
21	ii=1;
22	end
23	sZ=prod(size(A));
24	minA=min(A(ii));
25	maxA=max(A(ii));
26	meanA=mean(A(ii));
27	sdA=sqrt( mean( (A(ii)-meanA).^2 ) );
28	nZ=prod(sum(sum(sum(~isfinite(A)))));
29	switch max(nargout)
30	case {0}
31	disp( ...
32	sprintf('Min %0.5g Max %0.5g Mean %0.5g SD %0.5g NaN %i/%i',...
33	minA,maxA,meanA,sdA,nZ,sZ) );
34	case {1}
35	varargout(1)={minA};
36	case {2}
37	varargout(1)={minA};
38	varargout(2)={maxA};
39	case {3}
40	varargout(1)={minA};
41	varargout(2)={maxA};
42	varargout(3)={meanA};
43	case {4}
44	varargout(1)={minA};
45	varargout(2)={maxA};
46	varargout(3)={meanA};
47	varargout(4)={sdA};
48	otherwise
49	error('Too many return arguments requested')
50	end