Eisenstein primes: Difference between revisions

m
Line 76:
lim = isqrt(graphlimitsquared)
arr = [Eisenstein(a, b) for a = -lim:lim, b = -lim:lim]
eprimes = sort!([filter(is_eisenstein_prime, arr[p[1]]), forlt p= in(x, y) -> enumeratenorm(arrx) if< is_eisenstein_primenorm(p[2]y)))],
lt = (x, y) -> norm(x) < norm(y))
for (i, c) in enumerate(eprimes)
if i <= printlimit
4,107

edits