Mercurial > hg > octave-nkf

diff scripts/specfun/isprime.m @ 19237:920a400929ca

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
* isprime.m: Return an empty array for empty input
author Jordi Gutiérrez Hermoso <jordigh@octave.org>
date Sat, 23 Aug 2014 17:38:56 -0400
parents 0976f9fccbbd
children 89e275a4f6f6
line wrap: on
line diff
--- a/scripts/specfun/isprime.m
+++ b/scripts/specfun/isprime.m
@@ -43,6 +43,11 @@
     error ("isprime: X must be a non-negative integer");
   endif
 
+  if (isempty (x))
+    t = x;
+    return;
+  endif
+
   maxn = max (x(:));
   ## generate prime table of suitable length.
   maxp = min (maxn, max (sqrt (maxn), 1e7)); # FIXME: threshold not optimized.