\[\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.5 s
Input Error: 24.1
Output Error: 18.2
Log: ⚲
Profile: 🕒
\(\begin{cases} {\left(\sqrt{\left(\left((e^{\log_* (1 + (e^{d} - 1)^* \cdot c)} - 1)^*\right) \bmod \left(\tan^{-1}_* \frac{8.61318337292339 \cdot 10^{-131}}{d}\right)\right)}\right)}^2 & \text{when } d \le 2.076811f-39 \\ \left(\left(\left(-c\right) \cdot (e^{-d} - 1)^*\right) \bmod \left(\tan^{-1}_* \frac{8.61318337292339 \cdot 10^{-131}}{-d}\right)\right) & \text{otherwise} \end{cases}\)

    if d < 2.076811f-39

    1. Started with
      \[\left(\left((e^{d} - 1)^* \cdot c\right) \bmod \left(\tan^{-1}_* \frac{8.61318337292339 \cdot 10^{-131}}{d}\right)\right)\]
      16.5
    2. Using strategy rm
      16.5
    3. Applied add-sqr-sqrt to get
      \[\color{red}{\left(\left((e^{d} - 1)^* \cdot c\right) \bmod \left(\tan^{-1}_* \frac{8.61318337292339 \cdot 10^{-131}}{d}\right)\right)} \leadsto \color{blue}{{\left(\sqrt{\left(\left((e^{d} - 1)^* \cdot c\right) \bmod \left(\tan^{-1}_* \frac{8.61318337292339 \cdot 10^{-131}}{d}\right)\right)}\right)}^2}\]
      16.5
    4. Using strategy rm
      16.5
    5. Applied expm1-log1p-u to get
      \[{\left(\sqrt{\left(\color{red}{\left((e^{d} - 1)^* \cdot c\right)} \bmod \left(\tan^{-1}_* \frac{8.61318337292339 \cdot 10^{-131}}{d}\right)\right)}\right)}^2 \leadsto {\left(\sqrt{\left(\color{blue}{\left((e^{\log_* (1 + (e^{d} - 1)^* \cdot c)} - 1)^*\right)} \bmod \left(\tan^{-1}_* \frac{8.61318337292339 \cdot 10^{-131}}{d}\right)\right)}\right)}^2\]
      16.6

    if 2.076811f-39 < d

    1. Started with
      \[\left(\left((e^{d} - 1)^* \cdot c\right) \bmod \left(\tan^{-1}_* \frac{8.61318337292339 \cdot 10^{-131}}{d}\right)\right)\]
      30.2
    2. Applied taylor to get
      \[\left(\left((e^{d} - 1)^* \cdot c\right) \bmod \left(\tan^{-1}_* \frac{8.61318337292339 \cdot 10^{-131}}{d}\right)\right) \leadsto \left(\left(-1 \cdot \frac{(e^{\frac{-1}{d}} - 1)^*}{c}\right) \bmod \left(\tan^{-1}_* \frac{8.61318337292339 \cdot 10^{-131}}{\frac{-1}{d}}\right)\right)\]
      29.9
    3. Taylor expanded around -inf to get
      \[\color{red}{\left(\left(-1 \cdot \frac{(e^{\frac{-1}{d}} - 1)^*}{c}\right) \bmod \left(\tan^{-1}_* \frac{8.61318337292339 \cdot 10^{-131}}{\frac{-1}{d}}\right)\right)} \leadsto \color{blue}{\left(\left(-1 \cdot \frac{(e^{\frac{-1}{d}} - 1)^*}{c}\right) \bmod \left(\tan^{-1}_* \frac{8.61318337292339 \cdot 10^{-131}}{\frac{-1}{d}}\right)\right)}\]
      29.9
    4. Applied taylor to get
      \[\left(\left(-1 \cdot \frac{(e^{\frac{-1}{d}} - 1)^*}{c}\right) \bmod \left(\tan^{-1}_* \frac{8.61318337292339 \cdot 10^{-131}}{\frac{-1}{d}}\right)\right) \leadsto \left(\left(-1 \cdot \left(c \cdot (e^{-1 \cdot d} - 1)^*\right)\right) \bmod \left(\tan^{-1}_* \frac{8.61318337292339 \cdot 10^{-131}}{-1 \cdot d}\right)\right)\]
      19.6
    5. Taylor expanded around inf to get
      \[\color{red}{\left(\left(-1 \cdot \left(c \cdot (e^{-1 \cdot d} - 1)^*\right)\right) \bmod \left(\tan^{-1}_* \frac{8.61318337292339 \cdot 10^{-131}}{-1 \cdot d}\right)\right)} \leadsto \color{blue}{\left(\left(-1 \cdot \left(c \cdot (e^{-1 \cdot d} - 1)^*\right)\right) \bmod \left(\tan^{-1}_* \frac{8.61318337292339 \cdot 10^{-131}}{-1 \cdot d}\right)\right)}\]
      19.6
    6. Applied simplify to get
      \[\left(\left(-1 \cdot \left(c \cdot (e^{-1 \cdot d} - 1)^*\right)\right) \bmod \left(\tan^{-1}_* \frac{8.61318337292339 \cdot 10^{-131}}{-1 \cdot d}\right)\right) \leadsto \left(\left(\left(-c\right) \cdot (e^{-d} - 1)^*\right) \bmod \left(\tan^{-1}_* \frac{8.61318337292339 \cdot 10^{-131}}{-d}\right)\right)\]
      19.6
    7. Applied final simplification
  1. Removed slow pow expressions

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)))