Damit lässt du die Primzahl "2" aber aus. Bzw. du müsstest sie spezialbehandeln.

Einen weiteren Performanceschub bekommt man, wenn man in is_prime nur bis zur Wurzel auf Teilbarkeit prüft.

Wenn man eh alle PZ von 1 bis N haben möchte, ist die Sieb-Variante auch sehr gut.