In[]:=
CloseKernels[]
Out[]=
{KernelObject[1,DESKTOP-TLQ2UBE,<defunct>],KernelObject[2,DESKTOP-TLQ2UBE,<defunct>],KernelObject[3,DESKTOP-TLQ2UBE,<defunct>],KernelObject[4,DESKTOP-TLQ2UBE,<defunct>],KernelObject[5,DESKTOP-TLQ2UBE,<defunct>],KernelObject[6,DESKTOP-TLQ2UBE,<defunct>],KernelObject[7,DESKTOP-TLQ2UBE,<defunct>],KernelObject[8,DESKTOP-TLQ2UBE,<defunct>],KernelObject[9,local,<defunct>],KernelObject[10,local,<defunct>],KernelObject[11,local,<defunct>],KernelObject[12,local,<defunct>],KernelObject[13,local,<defunct>],KernelObject[14,local,<defunct>],KernelObject[15,local,<defunct>],KernelObject[16,local,<defunct>]}
In[]:=
Needs["SubKernels`LocalKernels`"]Block[{$mathkernel=$mathkernel<>" -threadpriority=2"},LaunchKernels[]]
Out[]=
$Failed
Out[]=
{}
In[]:=
Out[]=
In[]:=
EntityInstance
Out[]=
EntityInstance
In[]:=
$Version
Out[]=
13.2.0 for Linux x86 (64-bit) (November 2, 2022)
In[]:=
5000 digits of the MRB constant
Out[]=
0.1878596424620671202485179340542732300559030949001387861720046840894772315646602137032966544331074969038423458562580190612313700947592266304389293488961841208373366260816136027381263793734352832125527639621714893217020762820621715167154084126804483635416719985197680252759893899391445798350556135096485210712078444230958681294976885269495642042555864836704410425279524710606660926339748341031157816786416689154600342222588380025455396892947114212218910509832871227730802003644521539053639505533220347062755115981282803951021926491467317629351619065981601866424582495069720338199295842093551516251439935760076459329128145170908242491588320416906640933443591480670556469280678700702811500938060693813938595336065798740556206234870432936073781956460310476395066489306136064552806751519350828083737671929686639810309494963749627738304984632456347931157530028921252329181619562697369707486576547607117801719578736830096590226066875365630551656736128815020143875613668655221067430537059103973575619148909369077798320355119336240463725349410542836369971702441855165483727935882200813448096105880203064781961959695375628783481233497638586301014072725292301472333336250918584024803704048881967676760119858111679169352796852044160027086137228688945101510291998853690572865928708687542549253379439534758970356331344038263888798665619598073351473990256577813317226107612797585272274277730898577492230597096257256271883675575297887925361687673940354321451362772549229313126276435732144621618778637715420542312822344629539653290332217147982028075984221065564890048536858707083268874877377635047689160983185536281667159108412193420164386000258508426556435006954832830120546193205155935040023508351261335921740897007329784277128967365161960225077117388084262325697885465378690462227085674874747093069357326668590856162823753865512432975647464914619179575869342996208149878536663170197264534260468378010759055148678719039578315060452444190757044511382058533398469219482879476486575931785958165274929778220959774409113714342169296245931753245373401299593995004917912983680848547143925846704238528608320053664510586678151196459676079196431734307671534498300497128869401656600427062111079053164721504556329943884005211152390168773115456961028369205036896108806031603660382896533239383524154510137534165673472607464891120088099838152046695415026377035573283592996630642717305158972116351999161135954670315408725287243998197872502746797388638897056867435377857981058556192492185716949135673462704077491448799682065482817465880642236348160780950777057939313495829806602825272128491688809230325290270059917755059615835919993190869393039736611646514858219972925337106768738686235047915879737968269847878082223410618789674667450680064404065553875213281494980700209858132220620109011265903449717410801063247564712834609549284370065147450218226120415643930308859826426256828126092491136733967235933714534216902560140050169469983875907342920361729301531400405936246406781400779475613077369732409923529464794580778164607696240864595666084141126399988575739429315226283898798436350719371486573491962025428443510411472841973814933807066222573191021481585745042886728477250434386718443149128948635448929492143259660847149600072534066215387561341325254274130158182476636432111506809477451406309160928029719327606796946860926362081763442272977546326737161110302220019498455407233859679972956745318490433826332931888160330054013690316104309973777863934393135621496549969937314205819065334661573835222280871390934331325238360305287172314811151029705856281299558991830381071966308132701670498617683168329529053798760030666570203435884960342104811488681216083611944605571913973292970683232645094571537170202325175198220852151188427534658891812172660333192804801174759046131898449720722068248919153933258002624627216176424468747789609764906070479605351740205227992021112876539835381740117955267003375831396078847726709215700142824833574188212706558826075722384346836425460624376294972557208442910901014692032297634008831340381640372911311314959883048662944964532193113574956453912385456662881566099860325447839862878014756977228191151861520827628240145565179425383818619479339940149759987240427388505644114375607952021157910255398765366539323158492006532957070559185522267542845736812670699565824674762111599567798871089943129103401442549765937135216907789286769520555697938362195599913809621446502826784457775198333561874495799884353199090550456424061091609429751367477262251350611686791529453643717739913922126483778806203872768088381032267284933085939078772734728272526200188560744654953143832715829738195578665687162281676008390755269474337748054271709338897975317018906895514393940484553829830198169253659797476778546865045859640238175930139996122667668117389670502033292349080513690067470281709256521198479569202858477768560978920686214307401593342577726291016675234940557410020728689850147927742910664313694152819958560981976523216255318623764260591320667123942398310748008900676975706071323330124919595164681348037017813465191315868
In[]:=
10000 digits of the MRB constant
Out[]=
0.1878596424620671202485179340542732300559030949001387861720046840894772315646602137032966544331074969038423458562580190612313700947592266304389293488961841208373366260816136027381263793734352832125527639621714893217020762820621715167154084126804483635416719985197680252759893899391445798350556135096485210712078444230958681294976885269495642042555864836704410425279524710606660926339748341031157816786416689154600342222588380025455396892947114212218910509832871227730802003644521539053639505533220347062755115981282803951021926491467317629351619065981601866424582495069720338199295842093551516251439935760076459329128145170908242491588320416906640933443591480670556469280678700702811500938060693813938595336065798740556206234870432936073781956460310476395066489306136064552806751519350828083737671929686639810309494963749627738304984632456347931157530028921252329181619562697369707486576547607117801719578736830096590226066875365630551656736128815020143875613668655221067430537059103973575619148909369077798320355119336240463725349410542836369971702441855165483727935882200813448096105880203064781961959695375628783481233497638586301014072725292301472333336250918584024803704048881967676760119858111679169352796852044160027086137228688945101510291998853690572865928708687542549253379439534758970356331344038263888798665619598073351473990256577813317226107612797585272274277730898577492230597096257256271883675575297887925361687673940354321451362772549229313126276435732144621618778637715420542312822344629539653290332217147982028075984221065564890048536858707083268874877377635047689160983185536281667159108412193420164386000258508426556435006954832830120546193205155935040023508351261335921740897007329784277128967365161960225077117388084262325697885465378690462227085674874747093069357326668590856162823753865512432975647464914619179575869342996208149878536663170197264534260468378010759055148678719039578315060452444190757044511382058533398469219482879476486575931785958165274929778220959774409113714342169296245931753245373401299593995004917912983680848547143925846704238528608320053664510586678151196459676079196431734307671534498300497128869401656600427062111079053164721504556329943884005211152390168773115456961028369205036896108806031603660382896533239383524154510137534165673472607464891120088099838152046695415026377035573283592996630642717305158972116351999161135954670315408725287243998197872502746797388638897056867435377857981058556192492185716949135673462704077491448799682065482817465880642236348160780950777057939313495829806602825272128491688809230325290270059917755059615835919993190869393039736611646514858219972925337106768738686235047915879737968269847878082223410618789674667450680064404065553875213281494980700209858132220620109011265903449717410801063247564712834609549284370065147450218226120415643930308859826426256828126092491136733967235933714534216902560140050169469983875907342920361729301531400405936246406781400779475613077369732409923529464794580778164607696240864595666084141126399988575739429315226283898798436350719371486573491962025428443510411472841973814933807066222573191021481585745042886728477250434386718443149128948635448929492143259660847149600072534066215387561341325254274130158182476636432111506809477451406309160928029719327606796946860926362081763442272977546326737161110302220019498455407233859679972956745318490433826332931888160330054013690316104309973777863934393135621496549969937314205819065334661573835222280871390934331325238360305287172314811151029705856281299558991830381071966308132701670498617683168329529053798760030666570203435884960342104811488681216083611944605571913973292970683232645094571537170202325175198220852151188427534658891812172660333192804801174759046131898449720722068248919153933258002624627216176424468747789609764906070479605351740205227992021112876539835381740117955267003375831396078847726709215700142824833574188212706558826075722384346836425460624376294972557208442910901014692032297634008831340381640372911311314959883048662944964532193113574956453912385456662881566099860325447839862878014756977228191151861520827628240145565179425383818619479339940149759987240427388505644114375607952021157910255398765366539323158492006532957070559185522267542845736812670699565824674762111599567798871089943129103401442549765937135216907789286769520555697938362195599913809621446502826784457775198333561874495799884353199090550456424061091609429751367477262251350611686791529453643717739913922126483778806203872768088381032267284933085939078772734728272526200188560744654953143832715829738195578665687162281676008390755269474337748054271709338897975317018906895514393940484553829830198169253659797476778546865045859640238175930139996122667668117389670502033292349080513690067470281709256521198479569202858477768560978920686214307401593342577726291016675234940557410020728689850147927742910664313694152819958560981976523216255318623764260591320667123942398310748008900676975706071323330124919595164681348037017813465191315868
In[]:=
Print["Start time is ",ds=DateString[],"."];prec=10000;(**Numberofrequireddecimals.*.*)ClearSystemCache[];T0=SessionTime[];expM[pre_]:=Module[{a,d,s,k,bb,c,end,iprec,xvals,x,pc,cores=16(*=4*numberofphysicalcores*),tsize=2^7,chunksize,start=1,ll,ctab,pr=Floor[1.005pre]},chunksize=cores*tsize;n=Floor[1.32pr];end=Ceiling[n/chunksize];Print["Iterations required: ",n];Print["Will give ",end," time estimates, each more accurate than the previous."];Print["Will stop at ",end*chunksize," iterations to ensure precsion of around ",pr," decimal places."];d=ChebyshevT[n,3];{b,c,s}={SetPrecision[-1,1.1*n],-d,0};iprec=Ceiling[pr/6912];Do[xvals=Flatten[Parallelize[Table[Table[ll=start+j*tsize+l;x=N[E^(Log[ll]/(ll)),iprec];pc=iprec;While[pc<pr/1024,pc=Min[3pc,pr/1024];x=SetPrecision[x,pc];y=x^ll-ll;x=x(1-2y/((ll+1)y+2llll));];(**N[Exp[Log[ll]/ll],pr/1024]**)x=SetPrecision[x,pr/256];xll=x^ll;z=(ll-xll)/xll;t=2ll-1;t2=t^2;x=x*(1+SetPrecision[4.5,pr/256](ll-1)/t2+(ll+1)z/(2llt)-SetPrecision[13.5,pr/256]ll(ll-1)1/(3llt2+t^3z));(*N[Exp[Log[ll]/ll],pr/256]*)x=SetPrecision[x,pr/64];xll=x^ll;z=(ll-xll)/xll;t=2ll-1;t2=t^2;x=x*(1+SetPrecision[4.5,pr/64](ll-1)/t2+(ll+1)z/(2llt)-SetPrecision[13.5,pr/64]ll(ll-1)1/(3llt2+t^3z));(**N[Exp[Log[ll]/ll],pr/64]**)x=SetPrecision[x,pr/16];xll=x^ll;z=(ll-xll)/xll;t=2ll-1;t2=t^2;x=x*(1+SetPrecision[4.5,pr/16](ll-1)/t2+(ll+1)z/(2llt)-SetPrecision[13.5,pr/16]ll(ll-1)1/(3llt2+t^3z));(**N[Exp[Log[ll]/ll],pr/16]**)x=SetPrecision[x,pr/4];xll=x^ll;z=(ll-xll)/xll;t=2ll-1;t2=t^2;x=x*(1+SetPrecision[4.5,pr/4](ll-1)/t2+(ll+1)z/(2llt)-SetPrecision[13.5,pr/4]ll(ll-1)1/(3llt2+t^3z));(**N[Exp[Log[ll]/ll],pr/4]**)x=SetPrecision[x,pr];xll=x^ll;z=(ll-xll)/xll;t=2ll-1;t2=t^2;x=x*(1+SetPrecision[4.5,pr](ll-1)/t2+(ll+1)z/(2llt)-SetPrecision[13.5,pr]ll(ll-1)1/(3llt2+t^3z));(*N[Exp[Log[ll]/ll],pr]*)x,{l,0,tsize-1}],{j,0,cores-1}]]];ctab=ParallelTable[Table[c=b-c;ll=start+l-2;b*=2(ll+n)(ll-n)/((ll+1)(2ll+1));c,{l,chunksize}],Method"Automatic"];s+=ctab.(xvals-1);start+=chunksize;st=SessionTime[]-T0;kc=k*chunksize;ti=(st)/(kc+10^-4)*(n)/(3600)/(24);Print[kc," iterations done in ",N[st,4]," seconds."," Should take ",N[ti,4]," days or ",N[ti*24*3600,4],"s, finish ",DatePlus[ds,ti],"."],{k,0,end-1}];N[-s/d,pr]];t2=Timing[MRBtest2=expM[prec];];Print["Finished on ",DateString[],". Proccessor time was ",t2[[1]]," s."];Print["Actual time was ",st];(*Print[*)MRBtest2(*]*)(*Remove(**)orenterMRBtest2toprintoutput*);Print["Enter MRBtest2 to print ",Floor[Precision[MRBtest2]]," digits"];Print["If you saved m3M, the difference between this and 3,014,991 known digits is ",N[MRBtest2-m3M,10]]