Gaussian primes: Difference between revisions

Content added Content deleted

Inline

@@ Line 169: / Line 169: @@
 """
 function isGaussianprime(n::T) where T <: Integer
-    r, c = real(n), imag(n)
+    r, c = abs(real(n)), imag(n)
+    return isprime(r * r + c * c) || r == abs(c) && isprime(r) && (r - 3) % 4 == 0
-    ar = abs(r)
-    return isprime(r * r + c * c) || ar == abs(c) && isprime(ar) && (ar - 3) % 4 == 0
 end