mersenneforum.org  

Go Back   mersenneforum.org > Prime Search Projects > And now for something completely different

Reply
 
Thread Tools
Old 2022-11-28, 19:16   #188
paulunderwood
 
paulunderwood's Avatar
 
Sep 2002
Database er0rr

5×29×31 Posts
Default

Quote:
Originally Posted by sweety439 View Post
OK, also you can see https://oeis.org/A092559, current the larger prime factor of 2^A092559(n)+1 have been proven prime for all a(n) up to a(52) = 117239, why you skip a(53) = 120788 and prove a(54) = 127031?
https://primes.utm.edu/top20/page.php?id=67 was used as a source, not OEIS.
paulunderwood is offline   Reply With Quote
Old 2022-11-28, 20:49   #189
R. Gerbicz
 
R. Gerbicz's Avatar
 
"Robert Gerbicz"
Oct 2005
Hungary

5×17×19 Posts
Default

Quote:
Originally Posted by paulunderwood View Post
This number failed at Size [8] 215618 bits with: mu_div_qr GNU MP assertion failed cy == 0.

I am now attempting to certify the cofactor of M221509 (prp66673).
Quote:
Originally Posted by paulunderwood View Post
The same error has occurred at Size [2]. Where is this bug?

Taking W127031.
Have you used for this my corn code ?
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].
R. Gerbicz is offline   Reply With Quote
Old 2022-11-28, 20:55   #190
paulunderwood
 
paulunderwood's Avatar
 
Sep 2002
Database er0rr

118F16 Posts
Default

Quote:
Originally Posted by R. Gerbicz View Post
Have you used for this my corn code ?
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].
Of course I am using your corn.c. I'd be a fool not to utilize the huge speed up.

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 2022-11-28 at 21:01
paulunderwood is offline   Reply With Quote
Old 2022-11-28, 21:50   #191
chalsall
If I May
 
chalsall's Avatar
 
"Chris Halsall"
Sep 2002
Barbados

255178 Posts
Default

Quote:
Originally Posted by paulunderwood View Post
Of course I am using your corn.c. I'd be a fool not to utilize the huge speed up.
Guys... Yoga breaths...

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.
chalsall is offline   Reply With Quote
Old 2022-11-28, 22:27   #192
paulunderwood
 
paulunderwood's Avatar
 
Sep 2002
Database er0rr

5×29×31 Posts
Default

Quote:
Originally Posted by chalsall View Post
Guys... Yoga breaths...

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.
There was no animosity intended. I was merely conveying how much I admire Robert's work that is corn.c, and agreeing with him that the error is most likely with GMP's mpn/mu_div_qr.c somewhere. That is one hell of a complicated bit of code and should be left to the authorities to debug.

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 2022-11-28 at 22:29
paulunderwood is offline   Reply With Quote
Old 2022-11-28, 22:35   #193
chalsall
If I May
 
chalsall's Avatar
 
"Chris Halsall"
Sep 2002
Barbados

11,087 Posts
Default

Quote:
Originally Posted by paulunderwood View Post
There was no animosity intended. I am hoping this bug is fixed because I really want to have another go at certifying R86453 even if that takes 2 years!
Super cool.

Some take questioning as an attack. When the exact opposite is intended.
chalsall is offline   Reply With Quote
Old 2022-11-28, 23:19   #194
chalsall
If I May
 
chalsall's Avatar
 
"Chris Halsall"
Sep 2002
Barbados

11,087 Posts
Default

Quote:
Originally Posted by chalsall View Post
Super cool. Some take questioning as an attack. When the exact opposite is intended.
Could I please ask a sincere question?

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...
chalsall is offline   Reply With Quote
Old 2022-11-29, 05:06   #195
VBCurtis
 
VBCurtis's Avatar
 
"Curtis"
Feb 2005
Riverside, CA

33×11×19 Posts
Default

Quote:
Originally Posted by chalsall View Post
Like... Seriously... Get over yourself... You are not that important...
Please take this part of your own advice, and stop posting noise in thread after thread. We're tired of your non sequiturs.
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.
VBCurtis is online now   Reply With Quote
Old 2022-11-29, 20:32   #196
gorignak
 
Oct 2022

5 Posts
Default

I'd like to reserve M130439 cofactor PRP39261.
gorignak is offline   Reply With Quote
Old 2022-11-30, 17:58   #197
R. Gerbicz
 
R. Gerbicz's Avatar
 
"Robert Gerbicz"
Oct 2005
Hungary

64F16 Posts
Default

Quote:
Originally Posted by paulunderwood View Post
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.
What I can see even running on an 8k digits input that mpn_preinv_mu_div_qr is called only in trial factorization step.
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^2-L 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;
}
It will save also the output. Maybe this is already an overshoot because on the code you'll use smaller numbers due to constants. (With 16 GB Ram it could run, for smaller machines it could die).

Last fiddled with by R. Gerbicz on 2022-11-30 at 18:24 Reason: correction
R. Gerbicz is offline   Reply With Quote
Old 2022-11-30, 19:26   #198
paulunderwood
 
paulunderwood's Avatar
 
Sep 2002
Database er0rr

5×29×31 Posts
Default

The results of the tests are attached.
Attached Files
File Type: txt tests.txt (14.7 KB, 8 views)
paulunderwood is offline   Reply With Quote
Reply

Thread Tools


Similar Threads
Thread Thread Starter Forum Replies Last Post
15*2^n-1, n>1M Reservation Thread Kosmaj Riesel Prime Search 707 2021-11-20 19:13
Where on Earth can I download FastECPP? How fast is it compared with Primo? Stargate38 Miscellaneous Math 5 2021-11-16 17:33
5*2^n-1 Reservation Thread VBCurtis Riesel Prime Search 679 2021-10-09 17:33
Octoproth Reservation Thread Greenbank Octoproth Search 2 2007-12-26 09:58
Dodecaproth Reservation Thread Greenbank Octoproth Search 30 2006-02-09 00:33

All times are UTC. The time now is 21:48.


Mon Feb 6 21:48:47 UTC 2023 up 172 days, 19:17, 1 user, load averages: 1.05, 1.13, 1.14

Powered by vBulletin® Version 3.8.11
Copyright ©2000 - 2023, Jelsoft Enterprises Ltd.

This forum has received and complied with 0 (zero) government requests for information.

Permission is granted to copy, distribute and/or modify this document under the terms of the GNU Free Documentation License, Version 1.2 or any later version published by the Free Software Foundation.
A copy of the license is included in the FAQ.

≠ ± ∓ ÷ × · − √ ‰ ⊗ ⊕ ⊖ ⊘ ⊙ ≤ ≥ ≦ ≧ ≨ ≩ ≺ ≻ ≼ ≽ ⊏ ⊐ ⊑ ⊒ ² ³ °
∠ ∟ ° ≅ ~ ‖ ⟂ ⫛
≡ ≜ ≈ ∝ ∞ ≪ ≫ ⌊⌋ ⌈⌉ ∘ ∏ ∐ ∑ ∧ ∨ ∩ ∪ ⨀ ⊕ ⊗ 𝖕 𝖖 𝖗 ⊲ ⊳
∅ ∖ ∁ ↦ ↣ ∩ ∪ ⊆ ⊂ ⊄ ⊊ ⊇ ⊃ ⊅ ⊋ ⊖ ∈ ∉ ∋ ∌ ℕ ℤ ℚ ℝ ℂ ℵ ℶ ℷ ℸ 𝓟
¬ ∨ ∧ ⊕ → ← ⇒ ⇐ ⇔ ∀ ∃ ∄ ∴ ∵ ⊤ ⊥ ⊢ ⊨ ⫤ ⊣ … ⋯ ⋮ ⋰ ⋱
∫ ∬ ∭ ∮ ∯ ∰ ∇ ∆ δ ∂ ℱ ℒ ℓ
𝛢𝛼 𝛣𝛽 𝛤𝛾 𝛥𝛿 𝛦𝜀𝜖 𝛧𝜁 𝛨𝜂 𝛩𝜃𝜗 𝛪𝜄 𝛫𝜅 𝛬𝜆 𝛭𝜇 𝛮𝜈 𝛯𝜉 𝛰𝜊 𝛱𝜋 𝛲𝜌 𝛴𝜎𝜍 𝛵𝜏 𝛶𝜐 𝛷𝜙𝜑 𝛸𝜒 𝛹𝜓 𝛺𝜔