4 minute read

Catalan’s Constant

Please cite:
Kim, S. Normality Analysis of Current World Record Computations for Catalan’s Constant and Arc Length of a Lemniscate with a=1. arXiv Preprint arXiv:1908.08925
if this article or the calculated digits were useful.

This world record computation of 600,000,000,100 digits was done from Sat May 25 22:37:01 2019 to Tue Jun 18 18:59:44 2019 using the Pilehrood (2010-short) algorithm. This time, I have also verified the calculation using the Guillera (2008) algorithm from Fri Jun 7 11:13:58 2019 to Tue Jul 16 10:29:12 2019. Validation file generated by y-cruncher v0.7.7 Build 9501 for computation, and y-cruncher v0.7.7 Build 9499 for the verification run:
Computation: https://web.archive.org/web/20190724102605/http://www.numberworld.org/y-cruncher/records/2019_6_18_catalan.txt
Verification: https://web.archive.org/web/20190724102625/http://www.numberworld.org/y-cruncher/records/2019_7_16_catalan.txt

\[G = \beta(2) = \sum_{n=0}^{\infty} \frac{(-1)^{n}}{(2n+1)^2} = \frac{1}{1^2} - \frac{1}{3^2} + \frac{1}{5^2} - \frac{1}{7^2} + \frac{1}{9^2} - \cdots\]

The definition of Catalan’s constant (Wikipedia)

Catalan’s constant is defined by the above equation, where β is the Dirichlet beta function, which is closely related to the Riemann zeta function (I think both of them are learned by undergrad students in mathematics). Even though it looks irrational in a numerical scope, it is unproven if it is transcendental, or even irrational. For more information on Catalan’s constant, see this Wolfram Mathworld entry.

Same as the Lemniscate Constant, I have used y-cruncher by Mr. Alexander J. Yee for this computation. This program is commonly used for stress testing and benchmarking overclocked PC builds (obviously this program performs a very rigorous computation), along with fellow mathematical computing program Prime95.
Maintaining server stability was challenging, as the computation pushed all components of the computer to their limits.

Computation:
System information:
Operating System: Linux 3.10.0-693.21.1.el7.x86_64 x86_64 (CentOS 7)
Processor(s):
Name: Intel(R) Xeon(R) Gold 6140 CPU @ 2.30GHz
Logical Cores: 72
Physical Cores: 36
Sockets: 2
NUMA Nodes: 2
Base Frequency: 2,294,522,799 Hz

I have used two CPUs that are compatible with AVX-512 operations that are used crucially in vector operations like y-cruncher. However, y-cruncher has severe I/O bottlenecks (although I managed to select a better algorithm than the Lemniscate Constant and optimized some important aspects of the computation), so the number of cores did not assist my operation well. I have used RAID scratch storage for my operations, but the speed of the disk was very slow compared to the required R/W operations.

Start Date: Sat May 25 22:37:01 2019
End Date: Tue Jun 18 18:59:44 2019
Total Computation Time: 2028121.582 seconds
Start-to-End Wall Time: 2060562.370 seconds
CPU Utilization: 564.31 % + 151.19 % kernel overhead
Multi-core Efficiency: 7.84 % + 2.10 % kernel overhead

Although the multi-core efficiency improved slightly compared to the Lemniscate Constant, it still did not reach optimal efficiency (Dr. Ian Cutress’s Lemniscate Constant calculation had an efficiency of 94.04 % and CPU utilization of 9027.61 %, which means that the CPU was not bottlenecked by other factors).

Memory:
Usable Memory: 201,159,380,992 ( 187 GiB)
Logical Peak Disk Usage: 3,962,541,437,992 (3.60 TiB)
Logical Disk Bytes Read: 543,482,162,425,752 ( 494 TiB)
Logical Disk Bytes Written: 474,688,611,298,000 ( 432 TiB)

The use of a more efficient algorithm reduced disk operations by half compared to the Lemniscate Constant, resulting in a faster computation.
It is worth noting that HDD I/O speeds can be a bottleneck for virtually any other component, and using Optane DIMMs or more RAM may greatly improve computation speed.

Verification:
System information:
Operating System: Linux 3.10.0-693.21.1.el7.x86_64 x86_64 (CentOS 7)
Processor(s):
Name: Intel(R) Xeon(R) Gold 6140 CPU @ 2.30GHz
Logical Cores: 72
Physical Cores: 36
Sockets: 2
NUMA Nodes: 2
Base Frequency: 2,294,527,776 Hz

Start Date: Fri Jun 7 11:13:58 2019
End Date: Tue Jul 16 10:29:12 2019
Total Computation Time: 3218837.554 seconds
Start-to-End Wall Time: 3366914.430 seconds
CPU Utilization: 639.34 % + 150.50 % kernel overhead
Multi-core Efficiency: 8.88 % + 2.09 % kernel overhead

Memory:
Usable Memory: 201,226,489,856 ( 187 GiB)
Logical Peak Disk Usage: 3,986,470,844,768 (3.63 TiB)
Logical Disk Bytes Read: 947,069,197,181,784 ( 861 TiB)
Logical Disk Bytes Written: 829,531,225,016,496 ( 754 TiB)

A less efficient algorithm led to substantially more disk I/O, resulting to be similar R/W to the Lemniscate Constant.

Overall, the simpler series of Catalan’s constant, compared to the AGM-Pi algorithm of the Lemniscate Constant, allowed for a faster computation with less I/O bottleneck.

To view the digits for the Catalan’s constant, you can download it from This Link (note that the total size is over 1 TB, but the link will simply redirect to a registry with a download link).

Note that the digits are released under an Attribution-NonCommercial-NoDerivatives 4.0 International License, which prohibits commercial use and distribution of remixed, transformed, or built upon versions without consent. Proper attribution and indication of changes are required even if it is not a prohibited use case.

Archive for computation results in the y-cruncher website:
https://web.archive.org/web/20190722034426/http://www.numberworld.org/y-cruncher/
Special thanks to Mr. Alexander J. Yee for developing and releasing y-cruncher and providing advice, and the Internet Archive for hosting the computed digits.