Go Back > Great Internet Mersenne Prime Search > PrimeNet > GPU to 72

Thread Tools
Old 2021-05-13, 21:01   #12
If I May
chalsall's Avatar
"Chris Halsall"
Sep 2002

22·7·359 Posts

Originally Posted by LOBES View Post
Doh, nevermind.
chalsall is offline   Reply With Quote
Old 2021-05-14, 02:11   #13
1976 Toyota Corona years forever!
petrw1's Avatar
Nov 2006
Saskatchewan, Canada

113748 Posts

Originally Posted by chalsall View Post
There's some excellent GPU P-1 software available now. But Prime95/mprime (CPU) is by far the easiest to work with.
I agree with "easiest".

I am running GPUOwl under colab.
When I can get a P100 (which is virtually every day) I can do a 8GHz P-1 in 30 minutes.
Running a similar P-1 on all 8 cores of my i7-7820X takes about 75 minutes.
petrw1 is offline   Reply With Quote
Old 2021-05-14, 22:36   #14
ewmayer's Avatar
Sep 2002
Rep├║blica de California

22·3·7·139 Posts

(Somehow missed this when PaulU OPed it:)

Originally Posted by paulunderwood View Post
Congrats! However:

? factor(499400852887245323683941126088449355702834653807158087 )

[       290582744822559701357207 1]

[1718618403140893608084889221841 1]
Ha! Real men (a.k.a. execrable masochists) factor numbers like this using the world's slowest bignum code - I mean of course *nix 'bc', which IIRC uses base-10 emulation of your CPU's base-2 instructions or some such ludicrously inefficient bignum implementation - and crappy bc-based functions of their own writing:

n = 499400852887245323683941126088449355702834653807158087; p = 105032111;
Stage 1 prime-powers seed = 105032111
Stage 1 residue A = 275242671610725931867172664303887659718570581548948384, gcd(A-1,n) = 1
Stage 2 interval = [10000,5000000]:
Using base= 3; Initializing M*24 = 120 [base^(A^(b^2)) % n] buffers for Stage 2...
Stage 2 q0 = 10080, k0 = 48
At q = 209790
At q = 419790
At q = 629790
At q = 839790
At q = 1049790
At q = 1259790
At q = 1469790
At q = 1679790
At q = 1889790
At q = 2099790
At q = 2309790
At q = 2519790
At q = 2729790
At q = 2939790
At q = 3149790
At q = 3359790
At q = 3569790
At q = 3779790
At q = 3989790
At q = 4199790
At q = 4409790
At q = 4619790
At q = 4829790
Stage 2: did 23762 loop passes. Residue B = -409059575611368065569985294315721920104412510081096652, gcd(B,n) = 290582744822559701357207
This factor is a probable prime.
Processed 581936 stage 2 primes, including 234294 prime-pairs and 113348 prime-singles [80.52 % paired].

Now back to work on my cutting-edge bc-based NFS implementation, with which I hope to someday factor numbers as large as the quantum-computer folks do: "the quantum factorization of the largest number to date, 56,153, smashing the previous record of 143 that was set in 2012."

Last fiddled with by ewmayer on 2021-05-14 at 22:37
ewmayer is offline   Reply With Quote
Old 2021-05-23, 22:12   #15
ewmayer's Avatar
Sep 2002
Rep├║blica de California

22×3×7×139 Posts

[We open our next scene with a hand slapping the owner's forehead, accompanied by the utterance "doh!"]

Re above: In fact it seems silly to use powerful general-modulus factoring machinery like ECM or QS on such (p-1)-found factor-product composites. Here's why: say we have some product of prime factors F = f1*f2*...*fn discovered by running p-1 to stage bounds b1 and b2 on an input Mersenne M(p) (or other bigum modulus with factors of a known form, allowing p-1 to be 'seeded' with a component of same). BY DEFINITION, each prime factor f1-fn will be b1/b2-smooth, in the sense than fj = 2*p*C + 1, where C is a composite all of whose prime factors are <= b1, save possibly one outlier-prime factor > b1 and <= b2. Thus if we again run p-1 to bounds b1/b2, but now with arithmetic modulo the relatively tiny factor product F, we are guaranteed to resolve all the prime factors f1-fn - the only trick is that we will need to do multiple GCDs along the way in order to capture the individual prime factors f1,...,fn, rather than have this secondary p-1 run modulo F again produce the same composite GCD = F which the original p-1 run mod M(p) did. Again, though, since in the followup p-1 run we are working mod F, all the arithmetic is trivially cheap, including the needed GCDs. And since the cost of a p-1 run is effectively akin to a single super-cheap ECM curve, we've reduced the work of resolving the composite F to just that equivalent.

Last fiddled with by ewmayer on 2021-05-23 at 22:17
ewmayer is offline   Reply With Quote

Thread Tools

Similar Threads
Thread Thread Starter Forum Replies Last Post
Largest k*2^n-1 Primes Found In 2019 Thomas11 Riesel Prime Search 2 2020-01-01 05:32
Largest k*2^n-1 Primes Found In 2015 Kosmaj Riesel Prime Search 0 2016-01-01 13:45
Largest k*2^n-1 Primes Found in 2013 Kosmaj Riesel Prime Search 3 2014-12-12 07:14
formula for largest prime found debasish Miscellaneous Math 20 2007-09-28 03:48
Largest prime found yet by 15k search TTn 15k Search 1 2004-05-11 03:04

All times are UTC. The time now is 09:35.

Wed Dec 8 09:35:01 UTC 2021 up 138 days, 4:04, 1 user, load averages: 1.37, 1.74, 1.66

Powered by vBulletin® Version 3.8.11
Copyright ©2000 - 2021, 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.