Fixed – thanks,

]]>Yes. For example, see https://www.ams.org/journals/mcom/1975-29-130/S0025-5718-1975-0379537-3/

]]>By (3) it is $\max_{i,j} |a_{ij}|$.

]]>>>> import mpmath

>>> mpmath.harmonic(65)

mpf(‘4.7592755190903846’)

>>> mpmath.harmonic(513)

mpf(‘6.8184658522885142’)

>>> mpmath.harmonic(2097152)

mpf(‘15.133306695078945’)

>>> mpmath.harmonic(2.81e14)

mpf(‘33.846591450163828’)

The Malone paper you reference discusses the errors for float64 in great detail.

]]>