20220620, 19:22  #166 
Sep 2002
Database er0rr
2^{2}·1,063 Posts 
I have just put out (post 148) code that does autoincrement of FFT size on excessive roundoff error. Attached is the concept test program.

20220620, 20:19  #167 
Jul 2003
So Cal
2×23×53 Posts 
Code:
*** Maximum error of 0.500000 excessive at iteration 108154. Increasing FFT size. prp real 0m41.429s Code:
*** Maximum error of 0.375000 excessive at iteration 108660. Increasing FFT size. prp real 0m41.536s 
20220622, 19:19  #169 
Sep 2002
Database er0rr
10234_{8} Posts 
The ultimate Ver_6 has been uploaded to post 148. It uses functions and a neat piece of code to calculate odd d for n1=d*2^r thanks to frmky. It might be a little bit quicker. As I state there, comment out the printf(); function messaging about FFT increases if you find them annoying at runtime.
frmky is working on 3SPRP GMP code to be used above a certain threshold based this UTM page. This is useful for ARM based computers. Last fiddled with by paulunderwood on 20220622 at 19:27 
20220622, 19:42  #170  
Jul 2003
So Cal
986_{16} Posts 
Quote:
This code is slower than Paul's version using gwnum but is faster than mpz_probab_prime_p(). It's useful where gwnum is not available. Code:
int cm_nt_is_prime (mpz_t n) { unsigned long j, r; int prime = 0; mpz_t d, a, x, n_sub_1; if (mpz_sizeinbase(n, 2) < 4096) return (mpz_probab_prime_p(n, 0) > 0); if (mpz_even_p(n)) return 0; mpz_inits(d, a, x, n_sub_1, NULL); mpz_sub_ui(n_sub_1, n, 1); r = mpz_scan1(n_sub_1, 0); mpz_fdiv_q_2exp(d, n_sub_1, r); mpz_set_ui(a, 3); mpz_powm(x, a, d, n); if ((mpz_cmp_ui(x, 1) == 0)  (mpz_cmp(x, n_sub_1) == 0)) prime = 1; else { for (j = 1; j < r; j++) { mpz_powm_ui(x, x, 2, n); if (mpz_cmp(x, n_sub_1) == 0) { prime = 1; break; } } } mpz_clears(d, a, x, n_sub_1, NULL); return prime; } 

20220630, 07:58  #171  
Aug 2020
79*6581e4;3*2539e3
11·53 Posts 
edit: running sudo ldconfig solved the issue...
Quote:
Last fiddled with by bur on 20220630 at 08:12 

20220630, 08:06  #172  
Bamboozled!
"๐บ๐๐ท๐ท๐ญ"
May 2003
Down not across
11,423 Posts 
Quote:
Code:
export LD_LIBRARY_PATH=/usr/local/lib:$LD_LIBRARY_PATH 

20220630, 08:13  #173 
Aug 2020
79*6581e4;3*2539e3
1001000111_{2} Posts 
Thanks, I did ldconfig which fixed it.

20220630, 11:22  #174 
Aug 2020
79*6581e4;3*2539e3
11×53 Posts 
A question regarding the progress with the v switch, I noticed values like qroot, Cornacchia and trial div increase over time, are those the values to look out for? If so, up to which value will they increase?

20220630, 12:24  #175  
Sep 2002
Database er0rr
2^{2}·1,063 Posts 
Quote:
I suggest that you start with some small test cases like (10^10311)/9. 

20220630, 17:35  #176  
"99(4^34019)99 palind"
Nov 2016
(P^81993)SZ base 36
7^{2}×73 Posts 
Quote:


Thread Tools  
Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
For which types of primes is GPU primality test software available?  bur  GPU Computing  6  20200828 06:20 
Fastest software for Mersenne primality test?  JonathanM  Information & Answers  25  20200616 02:47 
APRCL as primality proof  f1pokerspeed  FactorDB  14  20140109 21:06 
Proof of Primality Test for Fermat Numbers  princeps  Math  15  20120402 21:49 
PRIMALITY PROOF for Wagstaff numbers!  AntonVrba  Math  96  20090225 10:37 