\[\left(\left((e^{d} - 1)^* \cdot c\right) \bmod \left(\tan^{-1}_* \frac{8.61318337292339 \cdot 10^{-131}}{d}\right)\right)\]
Test:
(fmod (* (expm1 d) c) (atan2 8.61318337292339e-131 d))
Bits:
128 bits
Bits error versus a
Bits error versus b
Bits error versus c
Bits error versus d
Time: 24.6 s
Input Error: 40.6
Output Error: 40.6
Log: ⚲
Profile: 🕒
\(\left(\left((e^{d} - 1)^* \cdot c\right) \bmod \left(\tan^{-1}_* \frac{8.61318337292339 \cdot 10^{-131}}{d}\right)\right)\)
  1. Started with
    \[\left(\left((e^{d} - 1)^* \cdot c\right) \bmod \left(\tan^{-1}_* \frac{8.61318337292339 \cdot 10^{-131}}{d}\right)\right)\]
    40.6

Original test:


(lambda ((a default) (b default) (c default) (d default))
  #:name "(fmod (* (expm1 d) c) (atan2 8.61318337292339e-131 d))"
  (fmod (* (expm1 d) c) (atan2 8.61318337292339e-131 d)))