| Alternative 1 | |
|---|---|
| Accuracy | 89.3% |
| Cost | 2633 |
(FPCore (x y z t a) :precision binary64 (+ x (/ (* (- y z) (- t x)) (- a z))))
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (- x (/ (* (- y z) (- x t)) (- a z)))))
(if (<= t_1 -5e-294)
(- x (* (- t x) (/ (- z y) (- a z))))
(if (<= t_1 0.0)
(fma x (/ (- y a) z) (+ t (* (/ t z) (- a y))))
(+ x (/ (- t x) (/ (- a z) (- y z))))))))double code(double x, double y, double z, double t, double a) {
return x + (((y - z) * (t - x)) / (a - z));
}
double code(double x, double y, double z, double t, double a) {
double t_1 = x - (((y - z) * (x - t)) / (a - z));
double tmp;
if (t_1 <= -5e-294) {
tmp = x - ((t - x) * ((z - y) / (a - z)));
} else if (t_1 <= 0.0) {
tmp = fma(x, ((y - a) / z), (t + ((t / z) * (a - y))));
} else {
tmp = x + ((t - x) / ((a - z) / (y - z)));
}
return tmp;
}
function code(x, y, z, t, a) return Float64(x + Float64(Float64(Float64(y - z) * Float64(t - x)) / Float64(a - z))) end
function code(x, y, z, t, a) t_1 = Float64(x - Float64(Float64(Float64(y - z) * Float64(x - t)) / Float64(a - z))) tmp = 0.0 if (t_1 <= -5e-294) tmp = Float64(x - Float64(Float64(t - x) * Float64(Float64(z - y) / Float64(a - z)))); elseif (t_1 <= 0.0) tmp = fma(x, Float64(Float64(y - a) / z), Float64(t + Float64(Float64(t / z) * Float64(a - y)))); else tmp = Float64(x + Float64(Float64(t - x) / Float64(Float64(a - z) / Float64(y - z)))); end return tmp end
code[x_, y_, z_, t_, a_] := N[(x + N[(N[(N[(y - z), $MachinePrecision] * N[(t - x), $MachinePrecision]), $MachinePrecision] / N[(a - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(x - N[(N[(N[(y - z), $MachinePrecision] * N[(x - t), $MachinePrecision]), $MachinePrecision] / N[(a - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, -5e-294], N[(x - N[(N[(t - x), $MachinePrecision] * N[(N[(z - y), $MachinePrecision] / N[(a - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$1, 0.0], N[(x * N[(N[(y - a), $MachinePrecision] / z), $MachinePrecision] + N[(t + N[(N[(t / z), $MachinePrecision] * N[(a - y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x + N[(N[(t - x), $MachinePrecision] / N[(N[(a - z), $MachinePrecision] / N[(y - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
x + \frac{\left(y - z\right) \cdot \left(t - x\right)}{a - z}
\begin{array}{l}
t_1 := x - \frac{\left(y - z\right) \cdot \left(x - t\right)}{a - z}\\
\mathbf{if}\;t_1 \leq -5 \cdot 10^{-294}:\\
\;\;\;\;x - \left(t - x\right) \cdot \frac{z - y}{a - z}\\
\mathbf{elif}\;t_1 \leq 0:\\
\;\;\;\;\mathsf{fma}\left(x, \frac{y - a}{z}, t + \frac{t}{z} \cdot \left(a - y\right)\right)\\
\mathbf{else}:\\
\;\;\;\;x + \frac{t - x}{\frac{a - z}{y - z}}\\
\end{array}
| Original | 60.8% |
|---|---|
| Target | 81.5% |
| Herbie | 89.3% |
if (+.f64 x (/.f64 (*.f64 (-.f64 y z) (-.f64 t x)) (-.f64 a z))) < -5.0000000000000003e-294Initial program 65.3%
Simplified88.2%
[Start]65.3 | \[ x + \frac{\left(y - z\right) \cdot \left(t - x\right)}{a - z}
\] |
|---|---|
associate-*l/ [<=]88.2 | \[ x + \color{blue}{\frac{y - z}{a - z} \cdot \left(t - x\right)}
\] |
if -5.0000000000000003e-294 < (+.f64 x (/.f64 (*.f64 (-.f64 y z) (-.f64 t x)) (-.f64 a z))) < 0.0Initial program 5.9%
Simplified5.9%
[Start]5.9 | \[ x + \frac{\left(y - z\right) \cdot \left(t - x\right)}{a - z}
\] |
|---|---|
associate-*l/ [<=]5.9 | \[ x + \color{blue}{\frac{y - z}{a - z} \cdot \left(t - x\right)}
\] |
Applied egg-rr5.9%
[Start]5.9 | \[ x + \frac{y - z}{a - z} \cdot \left(t - x\right)
\] |
|---|---|
*-commutative [=>]5.9 | \[ x + \color{blue}{\left(t - x\right) \cdot \frac{y - z}{a - z}}
\] |
clear-num [=>]5.9 | \[ x + \left(t - x\right) \cdot \color{blue}{\frac{1}{\frac{a - z}{y - z}}}
\] |
un-div-inv [=>]5.9 | \[ x + \color{blue}{\frac{t - x}{\frac{a - z}{y - z}}}
\] |
Taylor expanded in z around inf 98.0%
Simplified98.1%
[Start]98.0 | \[ \left(-1 \cdot \frac{y \cdot \left(t - x\right)}{z} + t\right) - -1 \cdot \frac{a \cdot \left(t - x\right)}{z}
\] |
|---|---|
+-commutative [=>]98.0 | \[ \color{blue}{\left(t + -1 \cdot \frac{y \cdot \left(t - x\right)}{z}\right)} - -1 \cdot \frac{a \cdot \left(t - x\right)}{z}
\] |
associate--l+ [=>]98.0 | \[ \color{blue}{t + \left(-1 \cdot \frac{y \cdot \left(t - x\right)}{z} - -1 \cdot \frac{a \cdot \left(t - x\right)}{z}\right)}
\] |
associate-*r/ [=>]98.0 | \[ t + \left(\color{blue}{\frac{-1 \cdot \left(y \cdot \left(t - x\right)\right)}{z}} - -1 \cdot \frac{a \cdot \left(t - x\right)}{z}\right)
\] |
associate-*r/ [=>]98.0 | \[ t + \left(\frac{-1 \cdot \left(y \cdot \left(t - x\right)\right)}{z} - \color{blue}{\frac{-1 \cdot \left(a \cdot \left(t - x\right)\right)}{z}}\right)
\] |
div-sub [<=]98.0 | \[ t + \color{blue}{\frac{-1 \cdot \left(y \cdot \left(t - x\right)\right) - -1 \cdot \left(a \cdot \left(t - x\right)\right)}{z}}
\] |
distribute-lft-out-- [=>]98.0 | \[ t + \frac{\color{blue}{-1 \cdot \left(y \cdot \left(t - x\right) - a \cdot \left(t - x\right)\right)}}{z}
\] |
associate-*r/ [<=]98.0 | \[ t + \color{blue}{-1 \cdot \frac{y \cdot \left(t - x\right) - a \cdot \left(t - x\right)}{z}}
\] |
mul-1-neg [=>]98.0 | \[ t + \color{blue}{\left(-\frac{y \cdot \left(t - x\right) - a \cdot \left(t - x\right)}{z}\right)}
\] |
unsub-neg [=>]98.0 | \[ \color{blue}{t - \frac{y \cdot \left(t - x\right) - a \cdot \left(t - x\right)}{z}}
\] |
distribute-rgt-out-- [=>]98.0 | \[ t - \frac{\color{blue}{\left(t - x\right) \cdot \left(y - a\right)}}{z}
\] |
associate-/l* [=>]98.1 | \[ t - \color{blue}{\frac{t - x}{\frac{z}{y - a}}}
\] |
Taylor expanded in x around 0 98.0%
Simplified98.1%
[Start]98.0 | \[ \left(\frac{\left(y - a\right) \cdot x}{z} + t\right) - \frac{t \cdot \left(y - a\right)}{z}
\] |
|---|---|
associate--l+ [=>]98.0 | \[ \color{blue}{\frac{\left(y - a\right) \cdot x}{z} + \left(t - \frac{t \cdot \left(y - a\right)}{z}\right)}
\] |
*-commutative [=>]98.0 | \[ \frac{\color{blue}{x \cdot \left(y - a\right)}}{z} + \left(t - \frac{t \cdot \left(y - a\right)}{z}\right)
\] |
associate-*r/ [<=]98.2 | \[ \color{blue}{x \cdot \frac{y - a}{z}} + \left(t - \frac{t \cdot \left(y - a\right)}{z}\right)
\] |
associate-*l/ [<=]98.1 | \[ x \cdot \frac{y - a}{z} + \left(t - \color{blue}{\frac{t}{z} \cdot \left(y - a\right)}\right)
\] |
fma-def [=>]98.1 | \[ \color{blue}{\mathsf{fma}\left(x, \frac{y - a}{z}, t - \frac{t}{z} \cdot \left(y - a\right)\right)}
\] |
*-commutative [=>]98.1 | \[ \mathsf{fma}\left(x, \frac{y - a}{z}, t - \color{blue}{\left(y - a\right) \cdot \frac{t}{z}}\right)
\] |
if 0.0 < (+.f64 x (/.f64 (*.f64 (-.f64 y z) (-.f64 t x)) (-.f64 a z))) Initial program 66.6%
Simplified88.6%
[Start]66.6 | \[ x + \frac{\left(y - z\right) \cdot \left(t - x\right)}{a - z}
\] |
|---|---|
associate-*l/ [<=]88.6 | \[ x + \color{blue}{\frac{y - z}{a - z} \cdot \left(t - x\right)}
\] |
Applied egg-rr88.7%
[Start]88.6 | \[ x + \frac{y - z}{a - z} \cdot \left(t - x\right)
\] |
|---|---|
*-commutative [=>]88.6 | \[ x + \color{blue}{\left(t - x\right) \cdot \frac{y - z}{a - z}}
\] |
clear-num [=>]88.5 | \[ x + \left(t - x\right) \cdot \color{blue}{\frac{1}{\frac{a - z}{y - z}}}
\] |
un-div-inv [=>]88.7 | \[ x + \color{blue}{\frac{t - x}{\frac{a - z}{y - z}}}
\] |
Final simplification89.3%
| Alternative 1 | |
|---|---|
| Accuracy | 89.3% |
| Cost | 2633 |
| Alternative 2 | |
|---|---|
| Accuracy | 89.3% |
| Cost | 2632 |
| Alternative 3 | |
|---|---|
| Accuracy | 50.1% |
| Cost | 1900 |
| Alternative 4 | |
|---|---|
| Accuracy | 49.7% |
| Cost | 1900 |
| Alternative 5 | |
|---|---|
| Accuracy | 49.6% |
| Cost | 1900 |
| Alternative 6 | |
|---|---|
| Accuracy | 51.0% |
| Cost | 1764 |
| Alternative 7 | |
|---|---|
| Accuracy | 50.8% |
| Cost | 1636 |
| Alternative 8 | |
|---|---|
| Accuracy | 53.6% |
| Cost | 1636 |
| Alternative 9 | |
|---|---|
| Accuracy | 52.9% |
| Cost | 1636 |
| Alternative 10 | |
|---|---|
| Accuracy | 62.1% |
| Cost | 1633 |
| Alternative 11 | |
|---|---|
| Accuracy | 62.7% |
| Cost | 1633 |
| Alternative 12 | |
|---|---|
| Accuracy | 59.1% |
| Cost | 1632 |
| Alternative 13 | |
|---|---|
| Accuracy | 68.5% |
| Cost | 1500 |
| Alternative 14 | |
|---|---|
| Accuracy | 68.5% |
| Cost | 1500 |
| Alternative 15 | |
|---|---|
| Accuracy | 38.8% |
| Cost | 1372 |
| Alternative 16 | |
|---|---|
| Accuracy | 49.9% |
| Cost | 1372 |
| Alternative 17 | |
|---|---|
| Accuracy | 79.0% |
| Cost | 1232 |
| Alternative 18 | |
|---|---|
| Accuracy | 77.7% |
| Cost | 969 |
| Alternative 19 | |
|---|---|
| Accuracy | 48.7% |
| Cost | 712 |
| Alternative 20 | |
|---|---|
| Accuracy | 44.4% |
| Cost | 328 |
| Alternative 21 | |
|---|---|
| Accuracy | 29.5% |
| Cost | 64 |
herbie shell --seed 2023135
(FPCore (x y z t a)
:name "Graphics.Rendering.Chart.Axis.Types:invLinMap from Chart-1.5.3"
:precision binary64
:herbie-target
(if (< z -1.2536131056095036e+188) (- t (* (/ y z) (- t x))) (if (< z 4.446702369113811e+64) (+ x (/ (- y z) (/ (- a z) (- t x)))) (- t (* (/ y z) (- t x)))))
(+ x (/ (* (- y z) (- t x)) (- a z))))