[Octave-cvsupdate] octave-forge/extra/NaN/inst nantest.m,1.2,1.3

[Alois S. <sch...@us...>]

Update of /cvsroot/octave/octave-forge/extra/NaN/inst
In directory sc8-pr-cvs3.sourceforge.net:/tmp/cvs-serv1938

Modified Files:
	nantest.m 
Log Message:
add test for quantile

Index: nantest.m
===================================================================
RCS file: /cvsroot/octave/octave-forge/extra/NaN/inst/nantest.m,v
retrieving revision 1.2
retrieving revision 1.3
diff -u -d -r1.2 -r1.3
--- nantest.m	23 Mar 2007 16:14:24 -0000	1.2
+++ nantest.m	28 Aug 2007 10:16:37 -0000	1.3
@@ -249,4 +249,28 @@
 
 tmp  = [tmp1;tmp2;tmp3;tmp4;tmp5;tmp6;tmp7;tmp8];
 
-warning(FLAG_WARNING);
+warning(FLAG_WARNING); 
+
+
+%%%%% QUANTILE TEST 
+d = [1 1 2 2 4 4 10 700]; 
+q = [-1,0,.05,.1,.25,.49,.5,.51,.75,.8, .999999,1,2];
+r = [ NaN, 1, 1, 1, 1.5, 2, 3, 4, 7, 10, 700, 700, NaN]; 
+if any( quantile(d, q) -  r>0)
+	fprintf(1,'Quantile(1): failed\n');
+else
+	fprintf(1,'Quantile(1): OK\n'); 
+end; 
+if exist('histo3','file')
+	H = histo3(d');
+else
+	H.X = [1;2;4;10;700];
+	H.H = [2;2;2;1;1];
+	H.datatype = 'histogram'; 
+end; 	 
+if any( quantile(H, q)' -  r>0)
+	fprintf(1,'Quantile(2): failed\n');
+else
+	fprintf(1,'Quantile(2): OK\n'); 
+end; 
+