Average Error: 4.1 → 4.3
Time: 15.0s
Precision: 64
\[\frac{x}{x + y \cdot e^{2 \cdot \left(\frac{z \cdot \sqrt{t + a}}{t} - \left(b - c\right) \cdot \left(\left(a + \frac{5}{6}\right) - \frac{2}{t \cdot 3}\right)\right)}}\]
\[\begin{array}{l} \mathbf{if}\;t \le -1.871833520333572462735300804758002311224 \cdot 10^{-43} \lor \neg \left(t \le 1.196097185809923326285406708328933897124 \cdot 10^{-103}\right):\\ \;\;\;\;\frac{x}{x + y \cdot e^{2 \cdot \left(\frac{z \cdot \sqrt{t + a}}{t} - \left(b - c\right) \cdot \left(\left(a + \frac{5}{6}\right) - \sqrt[3]{{\left(\frac{2}{t \cdot 3}\right)}^{3}}\right)\right)}}\\ \mathbf{else}:\\ \;\;\;\;\frac{x}{x + y \cdot e^{2 \cdot \frac{\left(z \cdot \sqrt{t + a}\right) \cdot \left(\left(a - \frac{5}{6}\right) \cdot \left(t \cdot 3\right)\right) - t \cdot \left(\left(b - c\right) \cdot \left(\left(a \cdot a - \frac{5}{6} \cdot \frac{5}{6}\right) \cdot \left(t \cdot 3\right) - \left(a - \frac{5}{6}\right) \cdot 2\right)\right)}{t \cdot \left(\left(a - \frac{5}{6}\right) \cdot \left(t \cdot 3\right)\right)}}}\\ \end{array}\]
\frac{x}{x + y \cdot e^{2 \cdot \left(\frac{z \cdot \sqrt{t + a}}{t} - \left(b - c\right) \cdot \left(\left(a + \frac{5}{6}\right) - \frac{2}{t \cdot 3}\right)\right)}}
\begin{array}{l}
\mathbf{if}\;t \le -1.871833520333572462735300804758002311224 \cdot 10^{-43} \lor \neg \left(t \le 1.196097185809923326285406708328933897124 \cdot 10^{-103}\right):\\
\;\;\;\;\frac{x}{x + y \cdot e^{2 \cdot \left(\frac{z \cdot \sqrt{t + a}}{t} - \left(b - c\right) \cdot \left(\left(a + \frac{5}{6}\right) - \sqrt[3]{{\left(\frac{2}{t \cdot 3}\right)}^{3}}\right)\right)}}\\

\mathbf{else}:\\
\;\;\;\;\frac{x}{x + y \cdot e^{2 \cdot \frac{\left(z \cdot \sqrt{t + a}\right) \cdot \left(\left(a - \frac{5}{6}\right) \cdot \left(t \cdot 3\right)\right) - t \cdot \left(\left(b - c\right) \cdot \left(\left(a \cdot a - \frac{5}{6} \cdot \frac{5}{6}\right) \cdot \left(t \cdot 3\right) - \left(a - \frac{5}{6}\right) \cdot 2\right)\right)}{t \cdot \left(\left(a - \frac{5}{6}\right) \cdot \left(t \cdot 3\right)\right)}}}\\

\end{array}
double f(double x, double y, double z, double t, double a, double b, double c) {
        double r98265 = x;
        double r98266 = y;
        double r98267 = 2.0;
        double r98268 = z;
        double r98269 = t;
        double r98270 = a;
        double r98271 = r98269 + r98270;
        double r98272 = sqrt(r98271);
        double r98273 = r98268 * r98272;
        double r98274 = r98273 / r98269;
        double r98275 = b;
        double r98276 = c;
        double r98277 = r98275 - r98276;
        double r98278 = 5.0;
        double r98279 = 6.0;
        double r98280 = r98278 / r98279;
        double r98281 = r98270 + r98280;
        double r98282 = 3.0;
        double r98283 = r98269 * r98282;
        double r98284 = r98267 / r98283;
        double r98285 = r98281 - r98284;
        double r98286 = r98277 * r98285;
        double r98287 = r98274 - r98286;
        double r98288 = r98267 * r98287;
        double r98289 = exp(r98288);
        double r98290 = r98266 * r98289;
        double r98291 = r98265 + r98290;
        double r98292 = r98265 / r98291;
        return r98292;
}

double f(double x, double y, double z, double t, double a, double b, double c) {
        double r98293 = t;
        double r98294 = -1.8718335203335725e-43;
        bool r98295 = r98293 <= r98294;
        double r98296 = 1.1960971858099233e-103;
        bool r98297 = r98293 <= r98296;
        double r98298 = !r98297;
        bool r98299 = r98295 || r98298;
        double r98300 = x;
        double r98301 = y;
        double r98302 = 2.0;
        double r98303 = z;
        double r98304 = a;
        double r98305 = r98293 + r98304;
        double r98306 = sqrt(r98305);
        double r98307 = r98303 * r98306;
        double r98308 = r98307 / r98293;
        double r98309 = b;
        double r98310 = c;
        double r98311 = r98309 - r98310;
        double r98312 = 5.0;
        double r98313 = 6.0;
        double r98314 = r98312 / r98313;
        double r98315 = r98304 + r98314;
        double r98316 = 3.0;
        double r98317 = r98293 * r98316;
        double r98318 = r98302 / r98317;
        double r98319 = 3.0;
        double r98320 = pow(r98318, r98319);
        double r98321 = cbrt(r98320);
        double r98322 = r98315 - r98321;
        double r98323 = r98311 * r98322;
        double r98324 = r98308 - r98323;
        double r98325 = r98302 * r98324;
        double r98326 = exp(r98325);
        double r98327 = r98301 * r98326;
        double r98328 = r98300 + r98327;
        double r98329 = r98300 / r98328;
        double r98330 = r98304 - r98314;
        double r98331 = r98330 * r98317;
        double r98332 = r98307 * r98331;
        double r98333 = r98304 * r98304;
        double r98334 = r98314 * r98314;
        double r98335 = r98333 - r98334;
        double r98336 = r98335 * r98317;
        double r98337 = r98330 * r98302;
        double r98338 = r98336 - r98337;
        double r98339 = r98311 * r98338;
        double r98340 = r98293 * r98339;
        double r98341 = r98332 - r98340;
        double r98342 = r98293 * r98331;
        double r98343 = r98341 / r98342;
        double r98344 = r98302 * r98343;
        double r98345 = exp(r98344);
        double r98346 = r98301 * r98345;
        double r98347 = r98300 + r98346;
        double r98348 = r98300 / r98347;
        double r98349 = r98299 ? r98329 : r98348;
        return r98349;
}

Error

Bits error versus x

Bits error versus y

Bits error versus z

Bits error versus t

Bits error versus a

Bits error versus b

Bits error versus c

Try it out

Your Program's Arguments

Results

Enter valid numbers for all inputs

Derivation

  1. Split input into 2 regimes
  2. if t < -1.8718335203335725e-43 or 1.1960971858099233e-103 < t

    1. Initial program 2.3

      \[\frac{x}{x + y \cdot e^{2 \cdot \left(\frac{z \cdot \sqrt{t + a}}{t} - \left(b - c\right) \cdot \left(\left(a + \frac{5}{6}\right) - \frac{2}{t \cdot 3}\right)\right)}}\]
    2. Using strategy rm
    3. Applied add-cbrt-cube2.3

      \[\leadsto \frac{x}{x + y \cdot e^{2 \cdot \left(\frac{z \cdot \sqrt{t + a}}{t} - \left(b - c\right) \cdot \left(\left(a + \frac{5}{6}\right) - \frac{2}{t \cdot \color{blue}{\sqrt[3]{\left(3 \cdot 3\right) \cdot 3}}}\right)\right)}}\]
    4. Applied add-cbrt-cube2.3

      \[\leadsto \frac{x}{x + y \cdot e^{2 \cdot \left(\frac{z \cdot \sqrt{t + a}}{t} - \left(b - c\right) \cdot \left(\left(a + \frac{5}{6}\right) - \frac{2}{\color{blue}{\sqrt[3]{\left(t \cdot t\right) \cdot t}} \cdot \sqrt[3]{\left(3 \cdot 3\right) \cdot 3}}\right)\right)}}\]
    5. Applied cbrt-unprod2.3

      \[\leadsto \frac{x}{x + y \cdot e^{2 \cdot \left(\frac{z \cdot \sqrt{t + a}}{t} - \left(b - c\right) \cdot \left(\left(a + \frac{5}{6}\right) - \frac{2}{\color{blue}{\sqrt[3]{\left(\left(t \cdot t\right) \cdot t\right) \cdot \left(\left(3 \cdot 3\right) \cdot 3\right)}}}\right)\right)}}\]
    6. Applied add-cbrt-cube2.3

      \[\leadsto \frac{x}{x + y \cdot e^{2 \cdot \left(\frac{z \cdot \sqrt{t + a}}{t} - \left(b - c\right) \cdot \left(\left(a + \frac{5}{6}\right) - \frac{\color{blue}{\sqrt[3]{\left(2 \cdot 2\right) \cdot 2}}}{\sqrt[3]{\left(\left(t \cdot t\right) \cdot t\right) \cdot \left(\left(3 \cdot 3\right) \cdot 3\right)}}\right)\right)}}\]
    7. Applied cbrt-undiv2.3

      \[\leadsto \frac{x}{x + y \cdot e^{2 \cdot \left(\frac{z \cdot \sqrt{t + a}}{t} - \left(b - c\right) \cdot \left(\left(a + \frac{5}{6}\right) - \color{blue}{\sqrt[3]{\frac{\left(2 \cdot 2\right) \cdot 2}{\left(\left(t \cdot t\right) \cdot t\right) \cdot \left(\left(3 \cdot 3\right) \cdot 3\right)}}}\right)\right)}}\]
    8. Simplified2.3

      \[\leadsto \frac{x}{x + y \cdot e^{2 \cdot \left(\frac{z \cdot \sqrt{t + a}}{t} - \left(b - c\right) \cdot \left(\left(a + \frac{5}{6}\right) - \sqrt[3]{\color{blue}{{\left(\frac{2}{t \cdot 3}\right)}^{3}}}\right)\right)}}\]

    if -1.8718335203335725e-43 < t < 1.1960971858099233e-103

    1. Initial program 6.9

      \[\frac{x}{x + y \cdot e^{2 \cdot \left(\frac{z \cdot \sqrt{t + a}}{t} - \left(b - c\right) \cdot \left(\left(a + \frac{5}{6}\right) - \frac{2}{t \cdot 3}\right)\right)}}\]
    2. Using strategy rm
    3. Applied flip-+10.3

      \[\leadsto \frac{x}{x + y \cdot e^{2 \cdot \left(\frac{z \cdot \sqrt{t + a}}{t} - \left(b - c\right) \cdot \left(\color{blue}{\frac{a \cdot a - \frac{5}{6} \cdot \frac{5}{6}}{a - \frac{5}{6}}} - \frac{2}{t \cdot 3}\right)\right)}}\]
    4. Applied frac-sub10.3

      \[\leadsto \frac{x}{x + y \cdot e^{2 \cdot \left(\frac{z \cdot \sqrt{t + a}}{t} - \left(b - c\right) \cdot \color{blue}{\frac{\left(a \cdot a - \frac{5}{6} \cdot \frac{5}{6}\right) \cdot \left(t \cdot 3\right) - \left(a - \frac{5}{6}\right) \cdot 2}{\left(a - \frac{5}{6}\right) \cdot \left(t \cdot 3\right)}}\right)}}\]
    5. Applied associate-*r/10.3

      \[\leadsto \frac{x}{x + y \cdot e^{2 \cdot \left(\frac{z \cdot \sqrt{t + a}}{t} - \color{blue}{\frac{\left(b - c\right) \cdot \left(\left(a \cdot a - \frac{5}{6} \cdot \frac{5}{6}\right) \cdot \left(t \cdot 3\right) - \left(a - \frac{5}{6}\right) \cdot 2\right)}{\left(a - \frac{5}{6}\right) \cdot \left(t \cdot 3\right)}}\right)}}\]
    6. Applied frac-sub7.7

      \[\leadsto \frac{x}{x + y \cdot e^{2 \cdot \color{blue}{\frac{\left(z \cdot \sqrt{t + a}\right) \cdot \left(\left(a - \frac{5}{6}\right) \cdot \left(t \cdot 3\right)\right) - t \cdot \left(\left(b - c\right) \cdot \left(\left(a \cdot a - \frac{5}{6} \cdot \frac{5}{6}\right) \cdot \left(t \cdot 3\right) - \left(a - \frac{5}{6}\right) \cdot 2\right)\right)}{t \cdot \left(\left(a - \frac{5}{6}\right) \cdot \left(t \cdot 3\right)\right)}}}}\]
  3. Recombined 2 regimes into one program.
  4. Final simplification4.3

    \[\leadsto \begin{array}{l} \mathbf{if}\;t \le -1.871833520333572462735300804758002311224 \cdot 10^{-43} \lor \neg \left(t \le 1.196097185809923326285406708328933897124 \cdot 10^{-103}\right):\\ \;\;\;\;\frac{x}{x + y \cdot e^{2 \cdot \left(\frac{z \cdot \sqrt{t + a}}{t} - \left(b - c\right) \cdot \left(\left(a + \frac{5}{6}\right) - \sqrt[3]{{\left(\frac{2}{t \cdot 3}\right)}^{3}}\right)\right)}}\\ \mathbf{else}:\\ \;\;\;\;\frac{x}{x + y \cdot e^{2 \cdot \frac{\left(z \cdot \sqrt{t + a}\right) \cdot \left(\left(a - \frac{5}{6}\right) \cdot \left(t \cdot 3\right)\right) - t \cdot \left(\left(b - c\right) \cdot \left(\left(a \cdot a - \frac{5}{6} \cdot \frac{5}{6}\right) \cdot \left(t \cdot 3\right) - \left(a - \frac{5}{6}\right) \cdot 2\right)\right)}{t \cdot \left(\left(a - \frac{5}{6}\right) \cdot \left(t \cdot 3\right)\right)}}}\\ \end{array}\]

Reproduce

herbie shell --seed 2019362 
(FPCore (x y z t a b c)
  :name "Numeric.SpecFunctions:invIncompleteBetaWorker from math-functions-0.1.5.2"
  :precision binary64
  (/ x (+ x (* y (exp (* 2 (- (/ (* z (sqrt (+ t a))) t) (* (- b c) (- (+ a (/ 5 6)) (/ 2 (* t 3)))))))))))