fixed gcd.i:is_prime function

LLNL · Mar 5, 2016 · bc24b71 · bc24b71
1 parent 01e228e
commit bc24b71
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/i/gcd.i b/i/gcd.i
@@ -72,7 +72,7 @@ func is_prime(x)
  */
 {
   x= long(abs(x));
-  if (x<2) return x==1;
+  if (x<4) return x>0;
   /* make a list of factors which includes 2, 3, and all larger
      odd numbers not divisible by 3 less or equal to sqrt(x) */
   top= min(long((sqrt(x)+8.)/3.+0.5), 20000);