mersenneforum.org

mersenneforum.org (https://www.mersenneforum.org/index.php)
-   Data (https://www.mersenneforum.org/forumdisplay.php?f=21)
-   -   Newer milestone thread (https://www.mersenneforum.org/showthread.php?t=13871)

NBtarheel_33 2015-02-04 20:00

February 3, 2015 All exponents below 33,219,253 (10 million digits) double-checked.

Just because I'm feeling extra nitpicky this afternoon: 33,219,2[B]53[/B] in the above should really be 33,219,2[B]81[/B].

Uncwilly 2015-02-05 01:13

[QUOTE=ATH;394443]Checking the milestones list against some old status files:

Milestones graph: [URL="www.hoegge.dk/mersenne/milestones.png"]milestones.png[/URL][/QUOTE]

It seems that you have some data points on that that I don't have (2008 and 2009 especially.) Would you be willing to share what you have? I can PM you my e-mail address.

Madpoo 2015-02-05 05:49

[QUOTE=NBtarheel_33;394455]February 3, 2015 All exponents below 33,219,253 (10 million digits) double-checked.

Just because I'm feeling extra nitpicky this afternoon: 33,219,2[B]53[/B] in the above should really be 33,219,2[B]81[/B].[/QUOTE]

I thought the first 2^n-1 that resulted in 10M digits is n=33219280 (which isn't prime, but doesn't matter). 2^33,219,281 - 1 would have more than 10M digits, wouldn't it?

I don't remember where or how I came up with 33219280 as the magic number, but that's what I was using as a basis. I don't have the energy to double-check right now... past my bedtime. :smile:

Madpoo 2015-02-05 05:54

[QUOTE=Uncwilly;394498]It seems that you have some data points on that that I don't have (2008 and 2009 especially.) Would you be willing to share what you have? I can PM you my e-mail address.[/QUOTE]

Maybe I can look at the milestone list and double-check, add the missing n-millionth items and verify the ones listed. They're probably okay.

axn 2015-02-05 06:06

[QUOTE=Madpoo;394509]I thought the first 2^n-1 that resulted in 10M digits is n=33219280 (which isn't prime, but doesn't matter). 2^33,219,281 - 1 would have more than 10M digits, wouldn't it?

I don't remember where or how I came up with 33219280 as the magic number, but that's what I was using as a basis. I don't have the energy to double-check right now... past my bedtime. :smile:[/QUOTE]

All exponents [B][I][U]below[/U][/I][/B] (I bolded, underlined & italicized for added effect :smile:)

retina 2015-02-05 06:13

[QUOTE=axn;394511]All exponents [B][I][U]below[/U][/I][/B] (I bolded, underlined & italicized for added effect :smile:)[/QUOTE]All exponents [color=red][B][I][U][size=5]below[/size][/U][/I][/B][/color] (I bolded, underlined, italicised coloured & enbiggened for added effect :smile:)

LaurV 2015-02-05 06:18

[QUOTE=Madpoo;394509]I thought the first 2^n-1 that resulted in 10M digits is n=33219280 (which isn't prime, but doesn't matter). 2^33,219,281 - 1 would have more than 10M digits, wouldn't it?

I don't remember where or how I came up with 33219280 as the magic number, but that's what I was using as a basis. I don't have the energy to double-check right now... past my bedtime. :smile:[/QUOTE]
Actually if we are to be nitpicking...:razz:
[CODE]
gp > ceil(332192[COLOR=Red][B]77[/B][/COLOR]*log(2)/log(10))
[COLOR=Magenta]time = 1 ms.
[/COLOR]%1 = 9999999
gp > ceil(332192[B][COLOR=Red]78[/COLOR][/B]*log(2)/log(10))
[COLOR=Magenta]time = 1 ms.
[/COLOR]%2 = 10000000
gp > length(Str(1<<13-1))
%3 = 4
gp > length(Str(1<<14-1))
%4 = 5
gp > length(Str(1<<332192[COLOR=Red][B]77[/B][/COLOR]-1))
[COLOR=Magenta]time = 14,928 ms.[/COLOR]
%5 = [B]9999999[/B]
gp > length(Str(1<<332192[COLOR=Red][B]78[/B][/COLOR]-1))
[COLOR=Magenta]time = 15,241 ms.[/COLOR]
%6 = [B]10000000[/B]
gp > precprime(332192[COLOR=Red][B]78[/B][/COLOR])
[COLOR=Magenta]time = 1 ms.[/COLOR]
%7 = 332192[B][COLOR=Red]53[/COLOR][/B]
gp > nextprime(332192[COLOR=Red][B]78[/B][/COLOR])
%8 = 332192[B][COLOR=Red]81[/COLOR][/B]
gp >
[/CODE]

(the log/ceil stuff is sometime tricky, if not enough precision in floats, here the precision is high enough, but usually "printing it" in a decimal string will resolve all the arguments, even if this is very slow, especially in pari/gp)

retina 2015-02-05 06:33

[QUOTE=LaurV;394513](the log/ceil stuff is sometime tricky, if not enough precision in floats, here the precision is high enough, but usually "printing it" in a decimal string will resolve all the arguments, even if this is very slow, especially in pari/gp)[/QUOTE]I would have thought the standard log function would suffice:

(10^7-1)/log10(2) = 33,219,277.62.....

And no more than 10 digits of precision needed.

The next decade is: (10^8-1)/log10(2) = 332,192,806.16...

LaurV 2015-02-05 06:39

[QUOTE=retina;394512]All exponents [SIZE=4][COLOR=Red][B]ʍo|ǝq[/B][/COLOR][/SIZE] (I [STRIKE]bolded, underlined, italicised coloured & enbiggened[/STRIKE][B] uʍop ǝpısdn ʇı pǝddı|ɟ[/B] for added effect :smile:)[/QUOTE]

noʎ ɹoɟ ʇɐɥʇ pǝxıɟ...

LaurV 2015-02-05 06:45

[QUOTE=retina;394514]I would have thought the standard log function would suffice:

(10^7-1)/log10(2) = 33,219,277.62.....

And no more than 10 digits of precision needed.

The next decade is: (10^8-1)/log10(2) = 332,192,806.16...[/QUOTE]

Try: (10^7)/log10(2) = ?
(and then start thinking what happens with the exponents in between...) :razz:

Madpoo 2015-02-05 07:02

[QUOTE=axn;394511]All exponents [B][I][U]below[/U][/I][/B] (I bolded, underlined & italicized for added effect :smile:)[/QUOTE]

I was confused what y'all were talking about, but I see now you're just being super precise (which I should appreciate, but I missed it).

So I should say "All exponents below *or equal to* blah blah have been double-checked"

I guess that's one thing in favor of using the language "up to" which is inclusive.

I'm looking at the milestone page and with all the stuff in there for the n-millionth stuff it's a little cluttered as time goes on, so maybe I'll try to organize it better... I'll make the change then.

FYI, looking at the database to try and figure out milestone dates is a little more daunting since triple-checks or other things are sometimes done, so I can't just look at the last date of a result in a certain exponent range. I'll have to do a little SQL magic to work out the actual date a single or double-check came in. The smaller exponents won't necessarily have that info in the v5 database either...maybe in the v4 database. :)


All times are UTC. The time now is 23:16.

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