asinh fails for input 1.3407807929942597e154 (and higher). The actual result should be a reasonable 355 (approx). The problem is the z * z in the formula overflows giving infinity.
A good fix is to test for any input with abs value g.t. 9.490626562425155e7 (for double), since after this zz + 1 == zz due to loss of precision. The asinh formula then is equivalent to log 2 + log z.
This should have been fixed by https://sourceforge.net/p/mingw-w64/mingw-w64/ci/021dffb8a482eb9d1b39569cd1ea42b87226fdf7/.
Last edit: David James 2021-10-14