Average Error: 0.0 → 0.6
Time: 22.2s
Precision: 64
\[\left(\left(\sinh c\right) \bmod \left(c - {\left( -2.9807307601812193 \cdot 10^{+165} \right)}^{2}\right)\right)\]
↓
\[\left(\left(\mathsf{fma}\left(c \cdot \left(c \cdot c\right), \frac{1}{6}, \mathsf{fma}\left(\frac{1}{120}, {c}^{5}, c\right)\right)\right) \bmod \left(c - -2.9807307601812193 \cdot 10^{+165} \cdot -2.9807307601812193 \cdot 10^{+165}\right)\right)\]
\left(\left(\sinh c\right) \bmod \left(c - {\left( -2.9807307601812193 \cdot 10^{+165} \right)}^{2}\right)\right)
↓
\left(\left(\mathsf{fma}\left(c \cdot \left(c \cdot c\right), \frac{1}{6}, \mathsf{fma}\left(\frac{1}{120}, {c}^{5}, c\right)\right)\right) \bmod \left(c - -2.9807307601812193 \cdot 10^{+165} \cdot -2.9807307601812193 \cdot 10^{+165}\right)\right)
double f(double c) {
        double r677690 = c;
        double r677691 = sinh(r677690);
        double r677692 = -2.9807307601812193e+165;
        double r677693 = 2.0;
        double r677694 = pow(r677692, r677693);
        double r677695 = r677690 - r677694;
        double r677696 = fmod(r677691, r677695);
        return r677696;
}
↓
double f(double c) {
        double r677697 = c;
        double r677698 = r677697 * r677697;
        double r677699 = r677697 * r677698;
        double r677700 = 0.16666666666666666;
        double r677701 = 0.008333333333333333;
        double r677702 = 5.0;
        double r677703 = pow(r677697, r677702);
        double r677704 = fma(r677701, r677703, r677697);
        double r677705 = fma(r677699, r677700, r677704);
        double r677706 = -2.9807307601812193e+165;
        double r677707 = r677706 * r677706;
        double r677708 = r677697 - r677707;
        double r677709 = fmod(r677705, r677708);
        return r677709;
}

Error

Bits error versus c

Derivation

  1. Initial program 0.0

    \[\left(\left(\sinh c\right) \bmod \left(c - {\left( -2.9807307601812193 \cdot 10^{+165} \right)}^{2}\right)\right)\]

  2. Simplified0.0

    \[\leadsto \color{blue}{\left(\left(\sinh c\right) \bmod \left(c - -2.9807307601812193 \cdot 10^{+165} \cdot -2.9807307601812193 \cdot 10^{+165}\right)\right)}\]

  3. Taylor expanded around 0 0.6

    \[\leadsto \left(\color{blue}{\left(\frac{1}{6} \cdot {c}^{3} + \left(\frac{1}{120} \cdot {c}^{5} + c\right)\right)} \bmod \left(c - -2.9807307601812193 \cdot 10^{+165} \cdot -2.9807307601812193 \cdot 10^{+165}\right)\right)\]

  4. Simplified0.6

    \[\leadsto \left(\color{blue}{\left(\mathsf{fma}\left(c \cdot \left(c \cdot c\right), \frac{1}{6}, \mathsf{fma}\left(\frac{1}{120}, {c}^{5}, c\right)\right)\right)} \bmod \left(c - -2.9807307601812193 \cdot 10^{+165} \cdot -2.9807307601812193 \cdot 10^{+165}\right)\right)\]

  5. Final simplification0.6

    \[\leadsto \left(\left(\mathsf{fma}\left(c \cdot \left(c \cdot c\right), \frac{1}{6}, \mathsf{fma}\left(\frac{1}{120}, {c}^{5}, c\right)\right)\right) \bmod \left(c - -2.9807307601812193 \cdot 10^{+165} \cdot -2.9807307601812193 \cdot 10^{+165}\right)\right)\]

Reproduce

herbie shell --seed 2019153 +o rules:numerics
(FPCore (c)
  :name "Random Jason Timeout Test 002"
  (fmod (sinh c) (- c (pow -2.9807307601812193e+165 2))))