# HG changeset patch
# User schloegl
# Date 1442043235 0
# Node ID 7c1bc8d8c4063cd212444f81ba62dbc5f0151d75
# Parent  f26b1170ea90e36e24c15c42b35d58f81ac06896
fix range check of k-index

diff -r f26b1170ea90 -r 7c1bc8d8c406 extra/NaN/src/kth_element.cpp
--- a/extra/NaN/src/kth_element.cpp	Sat Sep 12 07:15:01 2015 +0000
+++ b/extra/NaN/src/kth_element.cpp	Sat Sep 12 07:33:55 2015 +0000
@@ -190,10 +190,10 @@
 	POutput[0] = mxCreateDoubleMatrix(mxGetM(PInputs[1]),mxGetN(PInputs[1]),mxREAL);
 	Y = (double*) mxGetData(POutput[0]);
 	for (k=0; k < szK; k++) {
-		n = K[k]-1;       // convert to zero-based indexing 
-		if (n >= szX || n < 0)
+		if (K[k] > szX || K[k] < 1)
 			Y[k] = 0.0/0.0;	// NaN: result undefined
 		else {
+        		n = (size_t)(K[k]-1);       // convert to zero-based indexing, round towards 0 
         		findFirstK(T, 0, szX-1, n);
         		Y[k] = T[n];
 		}