| Karim Belabas on Mon, 26 May 2014 15:47:17 +0200 |
[Date Prev] [Date Next] [Thread Prev] [Thread Next] [Date Index] [Thread Index]
| Re: bug in PARI/GP (Segmentation Fault), please report. |
* John W. Nicholson [2014-05-26 15:00]:
> But then, I wonder what I would need to do if I need to have prime
> values like prime(10^10) or greater as to make it where it don't take
> over 5 minutes to calculate one prime?
Our implementation of prime(n) runs in time O~(n), it's a very expensive
way to iterate over primes. Use forprime:
(15:45) gp > n = 0; forprime(p=10^50, 10^50+10^7, n++); n
time = 30,094 ms.
%1 = 87112
Cheers,
K.B.
--
Karim Belabas, IMB (UMR 5251) Tel: (+33) (0)5 40 00 26 17
Universite Bordeaux 1 Fax: (+33) (0)5 40 00 69 50
351, cours de la Liberation http://www.math.u-bordeaux1.fr/~kbelabas/
F-33405 Talence (France) http://pari.math.u-bordeaux1.fr/ [PARI/GP]
`