mersenneforum.org

mersenneforum.org (https://www.mersenneforum.org/index.php)
-   Operazione Doppi Mersennes (https://www.mersenneforum.org/forumdisplay.php?f=99)
-   -   Trial division with CUDA (mmff) -- used, but runs like new! (https://www.mersenneforum.org/showthread.php?t=17162)

ET_ 2012-09-24 18:59

[QUOTE=Prime95;312590]Here we go again -- v 0.25:

[...]

My recommendation is to not upgrade until ATH, flashjh, and others have had time to try this version for a little bit. They have been quite effective in verifying the quality of recent releases.[/QUOTE]

Are you going to release a Linux 64bits binary of v 0.25?

Luigi

ATH 2012-09-24 22:45

Haven't found any issues with version 0.25 all my test cases works now :) Very nice.

It can now find 31 of the known fermat factors (remember PrintMode=1 in mmff.ini to avoid spam):

[CODE]worktodo.txt
FermatFactor=36,25709e6,25710e6
FermatFactor=33,5460e9,5470e9
FermatFactor=39,69,70
FermatFactor=45,11131e10,11132e10
FermatFactor=45,212e9,213e9
FermatFactor=50,2139e9,2140e9
FermatFactor=54,66,67
FermatFactor=54,78,79
FermatFactor=54,81900e9,81911e9
FermatFactor=61,67,68
FermatFactor=74,100,101
FermatFactor=77,98,99
FermatFactor=79,5e9,6e9
FermatFactor=87,1595e9,1596e9
FermatFactor=88,20018e9,20019e9
FermatFactor=90,119e9,120e9
FermatFactor=92,198e9,199e9
FermatFactor=93,103,104
FermatFactor=97,482e9,483e9
FermatFactor=101,3334e9,3335e9
FermatFactor=111,141,142
FermatFactor=120,3e9,4e9
FermatFactor=124,146,147
FermatFactor=127,129,130
FermatFactor=135,880e8,881e8
FermatFactor=145,167,168
FermatFactor=148,173,174
FermatFactor=149,160,161
FermatFactor=149,175,176
FermatFactor=154,166,167
FermatFactor=157,167,168

results.txt

F28 has a factor: 1766730974551267606529 [TF:70:71:mmff 0.25 mfaktc_barrett89_F32_63gs]
found 1 factor for k*2^36+1 in k range: 25709M to 25710M (71-bit factors) [mmff 0.25 mfaktc_barrett89_F32_63gs]
F31 has a factor: 46931635677864055013377 [TF:75:76:mmff 0.25 mfaktc_barrett89_F32_63gs]
found 1 factor for k*2^33+1 in k range: 5460G to 5470G (76-bit factors) [mmff 0.25 mfaktc_barrett89_F32_63gs]
F37 has a factor: 701179711390136401921 [TF:69:70:mmff 0.25 mfaktc_barrett89_F32_63gs]
found 1 factor for k*2^39+1 in k range: 1073741824 to 2147483647 (70-bit factors) [mmff 0.25 mfaktc_barrett89_F32_63gs]
F42 has a factor: 3916660235220715932328394753 [TF:91:92:mmff 0.25 mfaktc_barrett96_F32_63gs]
found 1 factor for k*2^45+1 in k range: 111310G to 111320G (92-bit factors) [mmff 0.25 mfaktc_barrett96_F32_63gs]
F43 has a factor: 7482850493766970889994241 [TF:82:83:mmff 0.25 mfaktc_barrett89_F32_63gs]
found 1 factor for k*2^45+1 in k range: 212G to 213G (83-bit factors) [mmff 0.25 mfaktc_barrett89_F32_63gs]
F48 has a factor: 2408911986953445595315961857 [TF:90:91:mmff 0.25 mfaktc_barrett96_F32_63gs]
found 1 factor for k*2^50+1 in k range: 2139G to 2140G (91-bit factors) [mmff 0.25 mfaktc_barrett96_F32_63gs]
F52 has a factor: 74201307460556292097 [TF:66:67:mmff 0.25 mfaktc_barrett89_F32_63gs]
found 1 factor for k*2^54+1 in k range: 4096 to 8191 (67-bit factors) [mmff 0.25 mfaktc_barrett89_F32_63gs]
F52 has a factor: 389591181597081096683521 [TF:78:79:mmff 0.25 mfaktc_barrett89_F32_63gs]
found 1 factor for k*2^54+1 in k range: 16777216 to 33554431 (79-bit factors) [mmff 0.25 mfaktc_barrett89_F32_63gs]
F52 has a factor: 1475547810493913550438096961537 [TF:100:101:mmff 0.25 mfaktc_barrett108_F32_63gs]
found 1 factor for k*2^54+1 in k range: 81900G to 81911G (101-bit factors) [mmff 0.25 mfaktc_barrett108_F32_63gs]
F58 has a factor: 219055085875300925441 [TF:67:68:mmff 0.25 mfaktc_barrett89_F32_63gs]
found 1 factor for k*2^61+1 in k range: 64 to 127 (68-bit factors) [mmff 0.25 mfaktc_barrett89_F32_63gs]
F72 has a factor: 1443765874709062348345951911937 [TF:100:101:mmff 0.25 mfaktc_barrett108_F64_95gs]
found 1 factor for k*2^74+1 in k range: 67108864 to 134217727 (101-bit factors) [mmff 0.25 mfaktc_barrett108_F64_95gs]
F75 has a factor: 520961043404985083798310879233 [TF:98:99:mmff 0.25 mfaktc_barrett108_F64_95gs]
found 1 factor for k*2^77+1 in k range: 2097152 to 4194303 (99-bit factors) [mmff 0.25 mfaktc_barrett108_F64_95gs]
F77 has a factor: 3590715923977960355577974656860161 [TF:111:112:mmff 0.25 mfaktc_barrett120_F64_95gs]
found 1 factor for k*2^79+1 in k range: 5G to 6G (112-bit factors) [mmff 0.25 mfaktc_barrett120_F64_95gs]
F83 has a factor: 246947940268608417020015902258307792897 [TF:127:128:mmff 0.25 mfaktc_barrett128_F64_95gs]
found 1 factor for k*2^87+1 in k range: 1595G to 1596G (128-bit factors) [mmff 0.25 mfaktc_barrett128_F64_95gs]
F86 has a factor: 6195449970597928748332522715641578258433 [TF:132:133:mmff 0.25 mfaktc_barrett140_F64_95gs]
found 1 factor for k*2^88+1 in k range: 20018G to 20019G (133-bit factors) [mmff 0.25 mfaktc_barrett140_F64_95gs]
F88 has a factor: 148481934042154969241780501829489000449 [TF:126:127:mmff 0.25 mfaktc_barrett128_F64_95gs]
found 1 factor for k*2^90+1 in k range: 119G to 120G (127-bit factors) [mmff 0.25 mfaktc_barrett128_F64_95gs]
F90 has a factor: 985016348367230226078056532654006730753 [TF:129:130:mmff 0.25 mfaktc_barrett140_F64_95gs]
found 1 factor for k*2^92+1 in k range: 198G to 199G (130-bit factors) [mmff 0.25 mfaktc_barrett140_F64_95gs]
F91 has a factor: 14072902366596202965053244178433 [TF:103:104:mmff 0.25 mfaktc_barrett108_F64_95gs]
found 1 factor for k*2^93+1 in k range: 1024 to 2047 (104-bit factors) [mmff 0.25 mfaktc_barrett108_F64_95gs]
F94 has a factor: 76459067246115642538831634131564386844673 [TF:135:136:mmff 0.25 mfaktc_barrett140_F96_127gs]
found 1 factor for k*2^97+1 in k range: 482G to 483G (136-bit factors) [mmff 0.25 mfaktc_barrett140_F96_127gs]
F96 has a factor: 8453027931784477309850388309101819121893377 [TF:142:143:mmff 0.25 mfaktc_barrett152_F96_127gs]
found 1 factor for k*2^101+1 in k range: 3334G to 3335G (143-bit factors) [mmff 0.25 mfaktc_barrett152_F96_127gs]
F107 has a factor: 3346902437331832346018436558958369334886401 [TF:141:142:mmff 0.25 mfaktc_barrett152_F96_127gs]
found 1 factor for k*2^111+1 in k range: 1073741824 to 2147483647 (142-bit factors) [mmff 0.25 mfaktc_barrett152_F96_127gs]
F116 has a factor: 4563438810603420826872624280490561141381005313 [TF:151:152:mmff 0.25 mfaktc_barrett152_F96_127gs]
found 1 factor for k*2^120+1 in k range: 3G to 4G (152-bit factors) [mmff 0.25 mfaktc_barrett152_F96_127gs]
F122 has a factor: 111331351706159727817280425663664652445286401 [TF:146:147:mmff 0.25 mfaktc_barrett152_F96_127gs]
found 1 factor for k*2^124+1 in k range: 4194304 to 8388607 (147-bit factors) [mmff 0.25 mfaktc_barrett152_F96_127gs]
F125 has a factor: 850705917302346158658436518579420528641 [TF:129:130:mmff 0.25 mfaktc_barrett140_F96_127gs]
found 1 factor for k*2^127+1 in k range: 4 to 7 (130-bit factors) [mmff 0.25 mfaktc_barrett140_F96_127gs]
F133 has a factor: 3836232386548105510567872577199319351015739156856833 [TF:171:172:mmff 0.25 mfaktc_barrett172_F128_159gs]
found 1 factor for k*2^135+1 in k range: 88000M to 88100M (172-bit factors) [mmff 0.25 mfaktc_barrett172_F128_159gs]
F142 has a factor: 363618066009591119386121910507749518730588867002369 [TF:167:168:mmff 0.25 mfaktc_barrett172_F128_159gs]
found 1 factor for k*2^145+1 in k range: 4194304 to 8388607 (168-bit factors) [mmff 0.25 mfaktc_barrett172_F128_159gs]
F146 has a factor: 13235038053749721162769301995307025251972223086886913 [TF:173:174:mmff 0.25 mfaktc_barrett183_F128_159gs]
found 1 factor for k*2^148+1 in k range: 33554432 to 67108863 (174-bit factors) [mmff 0.25 mfaktc_barrett183_F128_159gs]
F147 has a factor: 2230074519853062314153571827264836150598041600001 [TF:160:161:mmff 0.25 mfaktc_barrett172_F128_159gs]
found 1 factor for k*2^149+1 in k range: 2048 to 4095 (161-bit factors) [mmff 0.25 mfaktc_barrett172_F128_159gs]
F147 has a factor: 88894220732640180500173831441107513117330143465963521 [TF:175:176:mmff 0.25 mfaktc_barrett183_F128_159gs]
found 1 factor for k*2^149+1 in k range: 67108864 to 134217727 (176-bit factors) [mmff 0.25 mfaktc_barrett183_F128_159gs]
F150 has a factor: 124204803210043452689216278205372864748572142206977 [TF:166:167:mmff 0.25 mfaktc_barrett172_F128_159gs]
found 1 factor for k*2^154+1 in k range: 4096 to 8191 (167-bit factors) [mmff 0.25 mfaktc_barrett172_F128_159gs]
F150 has a factor: 287733134849521512021350451441018219494761719398401 [TF:167:168:mmff 0.25 mfaktc_barrett172_F128_159gs]
found 1 factor for k*2^157+1 in k range: 1024 to 2047 (168-bit factors) [mmff 0.25 mfaktc_barrett172_F128_159gs]
[/CODE]

Prime95 2012-09-25 02:06

1 Attachment(s)
Linux 64-bit build

firejuggler 2012-09-25 03:50

[code]

no factor for F46 from 2^96 to 2^97 (k range: 300000000000000 to 325000000000000) [mmff 0.20mmff mfaktc_barrett108_F30_61gs]
no factor for k*2^48+1 in k range: 325T to 350T (97-bit factors) [mmff 0.24 mfaktc_barrett108_F32_63gs]
no factor for k*2^48+1 in k range: 350T to 375T (97-bit factors) [mmff 0.25 mfaktc_barrett108_F32_63gs]
[/code]

Batalov 2012-09-25 11:52

The mmff-GFN patch
 
1 Attachment(s)
Because I can only dream of reaching George's level of generosity, the least I can do is spread the wealth.

Here's the patch the will turn mmff-0.25 into one of the five binaries for the Generalized Fermat factor search (bases 3,5,6,10,12 instead of 2). Read the included README, I will not repeat it here. It goes without saying that (if mmff is beta) this patch is alpha. Run the included tests!

Pick the k range and the base and good luck!
The reservations may be quite tricky (by email) but they are roughly summarized [URL="http://www1.uni-hamburg.de/RRZ/W.Keller/GFNsrch.txt"]here[/URL]; you may want to go above 10e12 (I've covered or in the process of covering the gap below 10e12).

The next stop will be xGFN (a[SUP]2[SUP]m[/SUP][/SUP] + b[SUP]2[SUP]m[/SUP][/SUP]). This search will be naturally slower, but could be fun, too!

-S

[B]EDIT: Caveat Emptor![/B] The mmff-gfn-0.25 binaries are [B]not [/B]capable of normal FermatFactoring. Keep the normal binary separately, and keep five GFN binaries separately - all in different folders.

Batalov 2012-09-25 12:24

1 Attachment(s)
P.S. The patch is a cleaner solution and its use is preferred (e.g. it can be applied to mmff-0.24 or maybe to mmff-0.26 later, with minimal hassle), but for the benefit for those who would prefer the source (and for Windows builders), here's the patched source, as well.

[B]EDIT: Caveat Emptor![/B] The mmff-gfn-0.25 binaries are [B]not [/B]capable of normal FermatFactoring. Keep the normal binary separately, and keep five GFN binaries separately - all in different folders.

The base-specific initializations are embedded into kernels, among other reasons - for speed. There are no switches. A future (xGFN) binary may be capable of all and any functions; if I will not think of a more elegant solution, there will be a sort of "fat binary" packing of a dozen of specialized kernels per each current kernel.

frmky 2012-09-26 07:08

[QUOTE=Batalov;312724](bases 3,5,6,10,12 instead of 2).[/QUOTE]
Any reason these bases were chosen, other than to check the same bases as Proth.exe?

Batalov 2012-09-26 07:20

I thought that too. It is easy to extend technically, but how would we know new factors from old? And who would volunteer to keep the factors and limits? (FactorDB has a limit on size; no pasaran).
It is probably because of the [URL="http://www.ams.org/journals/mcom/1998-67-221/S0025-5718-98-00891-6/S0025-5718-98-00891-6.pdf"]Bjorn/Riesel[/URL] legacy (and earlier Riesel references [3,4] therein).

frmky 2012-09-26 07:49

BR98 also includes bases 7, 8, and 11. If the idea is to only eliminate potential primes then these should be excluded, but so should 3 and 5. Perhaps the idea was to include the smallest bases, 3 & 5, plus bases that can include primes up to 12. Just a guess. :piggie:

Batalov 2012-09-26 08:00

Anyhow, W.Keller is keen on dealing with reservations on the 42[COLOR=red]*[/COLOR] (a,b) pairs and the five a's altogether. So, even though I would like to reserve the 25<=m<=150, k<=1e13, but I can only do that with xGFNs included. (reimplement "pfgw -gxo" in mmff-xgfn)

So, I plan to make mmff-xgfn. It would make sense: With five exponentiations (2,3,5,7,11), we will learn all GFN (including 7 and 11!) and xGFN factors in one sweep (some modular linear combinations need to be done; that's all)

GF'(3) and GF'(5) can be prime (after the /2)!
[I]F'[/I][SUB]6[/SUB](3) = 1716841910146256242328924544641 is prime, hey, larger than F(4)!

I also silently hope that every ~2000th found Proth prime at PrimeGrid could be a ticket, too. 2 found, ~1998 to go. ;-)
_______
[COLOR=red]*[/COLOR][SIZE=1]gotta love the number![/SIZE]

lalera 2012-09-29 17:36

hi,
i have mmff v 0.25 for linux x64 and the k - range is missing in the results.txt
is this a bug?
an example:
no factor for k*2^47+1 in (96-bit factors) [mmff 0.25 mfaktc_barrett96_F32_63gs]


All times are UTC. The time now is 04:47.

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