1st Run :
Start time is Sun 15 May 2022 16:25:21.
Iterations required: 13265998
Will give 6478 time estimates, each more accurate than the previous.
Will stop at 13266944 iterations to ensure precsion of around 10049999 decimal places.
As of Sun 15 May 2022 17:51:39 there were 2048 iterations done in 5178.5 seconds. That is 0.39548 iterations/s. 0.01543686% complete. It should take 388.237 days or 3.354×
7
10
s, and finish Wed 7 Jun 2023 22:07:18.
As of Sun 15 May 2022 18:35:36 there were 4096 iterations done in 7815.1 seconds. That is 0.52411 iterations/s. 0.03087373% complete. It should take 292.956 days or 2.531×
7
10
s, and finish Sat 4 Mar 2023 15:21:32.
As of Sun 15 May 2022 19:19:27 there were 6144 iterations done in 10447. seconds. That is 0.58814 iterations/s. 0.04631059% complete. It should take 261.065 days or 2.256×
7
10
s, and finish Tue 31 Jan 2023 17:58:33.
As of Sun 15 May 2022 20:03:13 there were 8192 iterations done in 13072. seconds. That is 0.62668 iterations/s. 0.06174745% complete. It should take 245.007 days or 2.117×
7
10
s, and finish Sun 15 Jan 2023 16:34:59.
As of Sun 15 May 2022 20:47:22 there were 10240 iterations done in 15721. seconds. That is 0.65136 iterations/s. 0.07718432% complete. It should take 235.726 days or 2.037×
7
10
s, and finish Fri 6 Jan 2023 09:50:38.
As of Sun 15 May 2022 21:31:16 there were 12288 iterations done in 18356. seconds. That is 0.66944 iterations/s. 0.09262118% complete. It should take 229.357 days or 1.982×
7
10
s, and finish Sat 31 Dec 2022 00:59:26.
As of Sun 15 May 2022 22:15:31 there were 14336 iterations done in 21010. seconds. That is 0.68233 iterations/s. 0.1080580% complete. It should take 225.024 days or 1.944×
7
10
s, and finish Mon 26 Dec 2022 16:59:57.
2ndrun
Start time is Sun 15 May 2022 23:20:17.
Iterations required: 13265998
Will give 6478 time estimates, each more accurate than the previous.
Will stop at 13266944 iterations to ensure precsion of around 10049999 decimal places.
As of Mon 16 May 2022 01:18:12 there were 2048 iterations done in 7074.4 seconds. That is 0.28949 iterations/s. 0.01543686% complete. It should take 530.378 days or 4.582×
7
10
s, and finish Sat 28 Oct 2023 08:25:04.
As of Mon 16 May 2022 02:01:49 there were 4096 iterations done in 9691.2 seconds. That is 0.42265 iterations/s. 0.03087373% complete. It should take 363.282 days or 3.139×
7
10
s, and finish Sun 14 May 2023 06:06:54.
$Aborted[]
Kernels
:Subkernel connected through KernelObject[10,DigitalStorm-PC] appears dead.
LinkObject
:Argument LinkObject[54002@192.168.1.212,54003@192.168.1.212,14674,20] in LinkClose[LinkObject[54002@192.168.1.212,54003@192.168.1.212,14674,20]] has an invalid LinkObject number; the link may be closed.
Kernels
:Subkernel connected through KernelObject[10,DigitalStorm-PC,<defunct>] appears dead.
LightweightGridClient`RemoteKernelClose
::nokernel
:Kernel could not be closed because no kernel was found for the given link "54002@192.168.1.212,54003@192.168.1.212.
​
LinkObject
:Unable to communicate with closed link LinkObject[54010@192.168.1.212,54011@192.168.1.212,14676,22].
Kernels
:Subkernel connected through KernelObject[12,DigitalStorm-PC] appears dead.
General
:Further output of Kernels::rdead will be suppressed during this calculation.
LinkObject
:Argument LinkObject[54010@192.168.1.212,54011@192.168.1.212,14676,22] in LinkClose[LinkObject[54010@192.168.1.212,54011@192.168.1.212,14676,22]] has an invalid LinkObject number; the link may be closed.
LightweightGridClient`RemoteKernelClose
::nokernel
:Kernel could not be closed because no kernel was found for the given link "54010@192.168.1.212,54011@192.168.1.212.
​
In[]:=
9
Out[]=
9
10milliondigitsoftheMRBconstantwith30kernelson30cores
​
In[]:=
CloseKernels[]
Out[]=
{}
In[]:=
Needs["SubKernels`LocalKernels`"]​​Block[{$mathkernel=$mathkernel<>" -threadpriority=2"},LaunchKernels[]]
Out[]=
{KernelObject[1,WIN-1AA39U1LQNT],KernelObject[2,WIN-1AA39U1LQNT],KernelObject[3,WIN-1AA39U1LQNT],KernelObject[4,WIN-1AA39U1LQNT],KernelObject[5,WIN-1AA39U1LQNT],KernelObject[6,WIN-1AA39U1LQNT],KernelObject[7,WIN-1AA39U1LQNT],KernelObject[8,WIN-1AA39U1LQNT],KernelObject[9,DigitalStorm-PC],KernelObject[10,DigitalStorm-PC],KernelObject[11,DigitalStorm-PC],KernelObject[12,DigitalStorm-PC],KernelObject[13,DigitalStorm-PC],KernelObject[14,DigitalStorm-PC],KernelObject[15,local],KernelObject[16,local],KernelObject[17,local],KernelObject[18,local],KernelObject[19,local],KernelObject[20,local],KernelObject[21,local],KernelObject[22,local],KernelObject[23,local],KernelObject[24,local],KernelObject[25,local],KernelObject[26,local],KernelObject[27,local],KernelObject[28,local],KernelObject[29,local],KernelObject[30,local]}
​
In[]:=
Print["Start time is ",ds=DateString[],"."];​​prec=10000000;​​(**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=pr/2^6;​​Do[xvals=Flatten[ParallelTable[Table[ll=start+j*tsize+l;​​x=N[E^(Log[ll]/(ll)),iprec];​​pc=iprec;​​While[pc<pr,pc=Min[4pc,pr];​​x=SetPrecision[x,pc];​​xll=x^ll;z=(ll-xll)/xll;​​t=2ll-1;t2=t^2;​​x=x*(1+SetPrecision[4.5,pc](ll-1)/t2+(ll+1)z/(2llt)-SetPrecision[13.5,2pc]ll(ll-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);​​If[kc>1,Print["As of ",DateString[]," there were ",kc," iterations done in ",N[st,5]," seconds. That is ",N[kc/st,5]," iterations/s. ",N[kc/(end*chunksize)*100,7],"% complete."," It should take ",N[ti,6]," days or ",N[ti*24*3600,4],"s, and finish ",DatePlus[ds,ti],"."]];​​Print[];,{k,0,end-1}];​​N[-s/d,pr]];​​t2=Timing[MRB1=expM[prec];];Print["Finished on ",DateString[],". Proccessor and actual time were ",t2[[1]]," and ",SessionTime[]-T0," s. respectively"];​​Print["Enter MRB1 to print ",Floor[Precision[MRB1]]," digits. The error from a 6,500,000 or more digit calculation that used a different method is "];N[MRB-MRB1,20]
Start time is Mon 16 May 2022 02:16:11.
Iterations required: 13265998
Will give 6478 time estimates, each more accurate than the previous.
Will stop at 13266944 iterations to ensure precsion of around 10049999 decimal places.
As of Mon 16 May 2022 03:42:17 there were 2048 iterations done in 5166.3 seconds. That is 0.39641 iterations/s. 0.01543686% complete. It should take 387.329 days or 3.347×
7
10
s, and finish Wed 7 Jun 2023 10:09:54.
As of Mon 16 May 2022 04:26:01 there were 4096 iterations done in 7790.1 seconds. That is 0.52580 iterations/s. 0.03087373% complete. It should take 292.017 days or 2.523×
7
10
s, and finish Sat 4 Mar 2023 02:40:14.
As of Mon 16 May 2022 05:10:14 there were 6144 iterations done in 10443. seconds. That is 0.58832 iterations/s. 0.04631059% complete. It should take 260.985 days or 2.255×
7
10
s, and finish Wed 1 Feb 2023 01:55:06.
As of Mon 16 May 2022 05:54:16 there were 8192 iterations done in 13085. seconds. That is 0.62605 iterations/s. 0.06174745% complete. It should take 245.254 days or 2.119×
7
10
s, and finish Mon 16 Jan 2023 08:22:02.
As of Mon 16 May 2022 06:38:31 there were 10240 iterations done in 15740. seconds. That is 0.65057 iterations/s. 0.07718432% complete. It should take 236.012 days or 2.039×
7
10
s, and finish Sat 7 Jan 2023 02:33:39.