Published In

Experimental Mathematics

Document Type

Article

Publication Date

10-2026

Subjects

Irrational rotation, discrepancy, Birkhoff sum, diophantine approximation, uniform distribution modulo one

Abstract

We present an algorithm for studying the sum๐บ๐‘โก(๐‘ฅ,๐œŒ)=๐‘โˆ‘๐‘˜=1[{๐‘ฅ+๐‘˜โข๐œŒ}โˆ’12],

which is closely related to the discrepancy of the irrational rotation. Computation of๐บ๐‘and associated statistics is notoriously difficult due to numerical instability and the failure of sampling to capture the discontinuities. Our algorithm instead fully defines๐บ๐‘through its discontinuities, which allows computation of๐บ๐‘โก(๐‘ฅ,๐œŒ)with errors on the order of๐‘times machine epsilon in๐‘‚โก(๐‘)time for practical๐‘. The associated probability density function can then be computed in๐‘‚โก(๐‘โขlogโก๐‘)time. Presented as a tool for experimental mathematics, the algorithm enables reliable plots of the descriptive statistics of๐บ๐‘over๐‘, in particular the length of the support of the pdf (the discrepancy), revealing extremal behavior related to the arithmetic properties of๐œŒ. We also observe a new pattern: when๐œŒis well approximated by a continued-fraction convergent๐‘๐‘›๐‘ž๐‘›, the pdf exhibits a predictable spiked-trapezoidal pattern when๐‘=๐‘˜โข๐‘ž๐‘›, eventually degrading at a rate that appears to depend on approximation quality.

Rights

Copyright (c) 2026 The Authors

Creative Commons License

This work is licensed under a Creative Commons Attribution 4.0 International License.

DOI

10.1080/10586458.2026.2742231

Persistent Identifier

https://archives.pdx.edu/ds/psu/45226

Share

COinS