20221128, 19:16  #188  
Sep 2002
Database er0rr
5×29×31 Posts 
Quote:


20221128, 20:49  #189  
"Robert Gerbicz"
Oct 2005
Hungary
5×17×19 Posts 
Quote:
Quote:
Notice that even if in some stages the code fails then likely it is a persistent bug in gmp and will fail later. In later stages you are doing the same algorithms just with a barely random smaller/larger numbers [Where the first stage could be a little exceptional if N has a special form]. 

20221128, 20:55  #190  
Sep 2002
Database er0rr
4495_{10} Posts 
Quote:
I had R86453 fail with this error besides these two, before corn.c. I think the problem lies with GMP and sheer size of the numbers being certified, but I could be wrong. Last fiddled with by paulunderwood on 20221128 at 21:01 

20221128, 21:50  #191  
If I May
"Chris Halsall"
Sep 2002
Barbados
11,087 Posts 
Quote:
Could we please just try to press the reset button from time to time? We are all on the same side. Or, at least, most of us. Correct? On the same page. Noise is something to be expected. Correct? P.S. I sometimes talk to the Universe. I find it catharitic. 

20221128, 22:27  #192  
Sep 2002
Database er0rr
5×29×31 Posts 
Quote:
I am hoping this bug is fixed because I really want to have another go at certifying R86453 even if that takes 2 years! Last fiddled with by paulunderwood on 20221128 at 22:29 

20221128, 22:35  #193 
If I May
"Chris Halsall"
Sep 2002
Barbados
11,087 Posts 

20221128, 23:19  #194  
If I May
"Chris Halsall"
Sep 2002
Barbados
10101101001111_{2} Posts 
Quote:
Why do people assume (Ass U Me) a question means an attack? Like... Seriously... Get over yourself... You are not that important... I almost feel like I resonate with Ricky Gervais... 

20221129, 05:06  #195  
"Curtis"
Feb 2005
Riverside, CA
3^{3}×11×19 Posts 
Quote:
Give us a break for a month: If your post doesn't have any math content, delete it instead of posting it. Give it a try. Help us cut the noise in the forum. Muzzle yourself. 

20221129, 20:32  #196 
Oct 2022
5_{16} Posts 
I'd like to reserve M130439 cofactor PRP39261.

20221130, 17:58  #197  
"Robert Gerbicz"
Oct 2005
Hungary
3117_{8} Posts 
Quote:
And most likely only on mpz_mod. Could you run this tiny test, testing mpz_mod (directly not on random numbers, but it is fine since there is no special code for p^e mod q^f, where p,q is odd): Code:
#include <stdio.h> #include <time.h> #include <math.h> #include "gmp.h" #define mult 1.4 // for the next bitlevel: nbits=mult*nbits. #define div 2 // to divide n by roughly n^(1/d) #define L 3 // to see more tests increase it, for given nbits it will do L^2L tests, for L<=7 #define maxnbits 12e9 int main(void){ int ptab[]={3,5,7,11,13,17,19}; mpz_t n,d,r; FILE* fout; fout=fopen("tests.txt","a+"); mpz_init(n); mpz_init(d); mpz_init(r); for(long nbits=1024;nbits<=maxnbits;nbits=(long)nbits*mult){ printf("Testing %ld bits n numbers.\n",nbits); for(int i=0;i<L&&i<7;i++){ long E1=(double)nbits*(log(2)/log(ptab[i])); mpz_ui_pow_ui(n,ptab[i],E1); for(int j=0;j<L&&j<7;j++)if(i!=j){ long E2=((double)1.0/div)*(double)nbits*(log(2)/log(ptab[j])); mpz_ui_pow_ui(d,ptab[j],E2); time_t sec=time(NULL); mpz_mod(r,n,d); int parity=(mpz_even_p(r)!=0?0:1); printf("calc parity(%d^%ld mod %d^%ld)=%d, time=%ld sec.\n", ptab[i],E1,ptab[j],E2,parity,time(NULL)sec); fprintf(fout,"calc parity(%d^%ld mod %d^%ld)=%d, time=%ld sec.\n", ptab[i],E1,ptab[j],E2,parity,time(NULL)sec); fflush(fout); } } } fclose(fout); mpz_clear(n); mpz_clear(d); mpz_clear(r); return 0; } Last fiddled with by R. Gerbicz on 20221130 at 18:24 Reason: correction 

20221130, 19:26  #198 
Sep 2002
Database er0rr
5·29·31 Posts 
The results of the tests are attached.

Thread Tools  
Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
15*2^n1, n>1M Reservation Thread  Kosmaj  Riesel Prime Search  707  20211120 19:13 
Where on Earth can I download FastECPP? How fast is it compared with Primo?  Stargate38  Miscellaneous Math  5  20211116 17:33 
5*2^n1 Reservation Thread  VBCurtis  Riesel Prime Search  679  20211009 17:33 
Octoproth Reservation Thread  Greenbank  Octoproth Search  2  20071226 09:58 
Dodecaproth Reservation Thread  Greenbank  Octoproth Search  30  20060209 00:33 