| Alternative 1 | |
|---|---|
| Accuracy | 88.5% |
| Cost | 4432 |
(FPCore (x y z t a) :precision binary64 (+ x (/ (* (- y x) (- z t)) (- a t))))
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (* (/ (- y x) t) (- a z))))
(if (<= t -2.2e+164)
(+ (+ y t_1) (* (/ a t) t_1))
(if (<= t 3.6e+142)
(+ x (/ (- y x) (/ (- a t) (- z t))))
(+ y (/ (- x y) (/ t (- z a))))))))double code(double x, double y, double z, double t, double a) {
return x + (((y - x) * (z - t)) / (a - t));
}
double code(double x, double y, double z, double t, double a) {
double t_1 = ((y - x) / t) * (a - z);
double tmp;
if (t <= -2.2e+164) {
tmp = (y + t_1) + ((a / t) * t_1);
} else if (t <= 3.6e+142) {
tmp = x + ((y - x) / ((a - t) / (z - t)));
} else {
tmp = y + ((x - y) / (t / (z - a)));
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
code = x + (((y - x) * (z - t)) / (a - t))
end function
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: t_1
real(8) :: tmp
t_1 = ((y - x) / t) * (a - z)
if (t <= (-2.2d+164)) then
tmp = (y + t_1) + ((a / t) * t_1)
else if (t <= 3.6d+142) then
tmp = x + ((y - x) / ((a - t) / (z - t)))
else
tmp = y + ((x - y) / (t / (z - a)))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
return x + (((y - x) * (z - t)) / (a - t));
}
public static double code(double x, double y, double z, double t, double a) {
double t_1 = ((y - x) / t) * (a - z);
double tmp;
if (t <= -2.2e+164) {
tmp = (y + t_1) + ((a / t) * t_1);
} else if (t <= 3.6e+142) {
tmp = x + ((y - x) / ((a - t) / (z - t)));
} else {
tmp = y + ((x - y) / (t / (z - a)));
}
return tmp;
}
def code(x, y, z, t, a): return x + (((y - x) * (z - t)) / (a - t))
def code(x, y, z, t, a): t_1 = ((y - x) / t) * (a - z) tmp = 0 if t <= -2.2e+164: tmp = (y + t_1) + ((a / t) * t_1) elif t <= 3.6e+142: tmp = x + ((y - x) / ((a - t) / (z - t))) else: tmp = y + ((x - y) / (t / (z - a))) return tmp
function code(x, y, z, t, a) return Float64(x + Float64(Float64(Float64(y - x) * Float64(z - t)) / Float64(a - t))) end
function code(x, y, z, t, a) t_1 = Float64(Float64(Float64(y - x) / t) * Float64(a - z)) tmp = 0.0 if (t <= -2.2e+164) tmp = Float64(Float64(y + t_1) + Float64(Float64(a / t) * t_1)); elseif (t <= 3.6e+142) tmp = Float64(x + Float64(Float64(y - x) / Float64(Float64(a - t) / Float64(z - t)))); else tmp = Float64(y + Float64(Float64(x - y) / Float64(t / Float64(z - a)))); end return tmp end
function tmp = code(x, y, z, t, a) tmp = x + (((y - x) * (z - t)) / (a - t)); end
function tmp_2 = code(x, y, z, t, a) t_1 = ((y - x) / t) * (a - z); tmp = 0.0; if (t <= -2.2e+164) tmp = (y + t_1) + ((a / t) * t_1); elseif (t <= 3.6e+142) tmp = x + ((y - x) / ((a - t) / (z - t))); else tmp = y + ((x - y) / (t / (z - a))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := N[(x + N[(N[(N[(y - x), $MachinePrecision] * N[(z - t), $MachinePrecision]), $MachinePrecision] / N[(a - t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(N[(N[(y - x), $MachinePrecision] / t), $MachinePrecision] * N[(a - z), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t, -2.2e+164], N[(N[(y + t$95$1), $MachinePrecision] + N[(N[(a / t), $MachinePrecision] * t$95$1), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 3.6e+142], N[(x + N[(N[(y - x), $MachinePrecision] / N[(N[(a - t), $MachinePrecision] / N[(z - t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(y + N[(N[(x - y), $MachinePrecision] / N[(t / N[(z - a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
x + \frac{\left(y - x\right) \cdot \left(z - t\right)}{a - t}
\begin{array}{l}
t_1 := \frac{y - x}{t} \cdot \left(a - z\right)\\
\mathbf{if}\;t \leq -2.2 \cdot 10^{+164}:\\
\;\;\;\;\left(y + t_1\right) + \frac{a}{t} \cdot t_1\\
\mathbf{elif}\;t \leq 3.6 \cdot 10^{+142}:\\
\;\;\;\;x + \frac{y - x}{\frac{a - t}{z - t}}\\
\mathbf{else}:\\
\;\;\;\;y + \frac{x - y}{\frac{t}{z - a}}\\
\end{array}
Results
| Original | 62.9% |
|---|---|
| Target | 85.4% |
| Herbie | 87.8% |
if t < -2.20000000000000006e164Initial program 25.4%
Simplified62.0%
[Start]25.4 | \[ x + \frac{\left(y - x\right) \cdot \left(z - t\right)}{a - t}
\] |
|---|---|
+-commutative [=>]25.4 | \[ \color{blue}{\frac{\left(y - x\right) \cdot \left(z - t\right)}{a - t} + x}
\] |
associate-*r/ [<=]62.0 | \[ \color{blue}{\left(y - x\right) \cdot \frac{z - t}{a - t}} + x
\] |
*-commutative [<=]62.0 | \[ \color{blue}{\frac{z - t}{a - t} \cdot \left(y - x\right)} + x
\] |
fma-def [=>]62.0 | \[ \color{blue}{\mathsf{fma}\left(\frac{z - t}{a - t}, y - x, x\right)}
\] |
Taylor expanded in t around inf 54.7%
Simplified85.9%
[Start]54.7 | \[ \frac{a \cdot \left(\left(y - x\right) \cdot \left(-1 \cdot z - -1 \cdot a\right)\right)}{{t}^{2}} + \left(\frac{\left(y - x\right) \cdot \left(-1 \cdot z - -1 \cdot a\right)}{t} + y\right)
\] |
|---|---|
+-commutative [=>]54.7 | \[ \color{blue}{\left(\frac{\left(y - x\right) \cdot \left(-1 \cdot z - -1 \cdot a\right)}{t} + y\right) + \frac{a \cdot \left(\left(y - x\right) \cdot \left(-1 \cdot z - -1 \cdot a\right)\right)}{{t}^{2}}}
\] |
*-commutative [<=]54.7 | \[ \left(\frac{\left(y - x\right) \cdot \left(-1 \cdot z - -1 \cdot a\right)}{t} + y\right) + \frac{a \cdot \color{blue}{\left(\left(-1 \cdot z - -1 \cdot a\right) \cdot \left(y - x\right)\right)}}{{t}^{2}}
\] |
*-commutative [=>]54.7 | \[ \left(\frac{\left(y - x\right) \cdot \left(-1 \cdot z - -1 \cdot a\right)}{t} + y\right) + \frac{\color{blue}{\left(\left(-1 \cdot z - -1 \cdot a\right) \cdot \left(y - x\right)\right) \cdot a}}{{t}^{2}}
\] |
distribute-lft-out-- [=>]54.7 | \[ \left(\frac{\left(y - x\right) \cdot \left(-1 \cdot z - -1 \cdot a\right)}{t} + y\right) + \frac{\left(\color{blue}{\left(-1 \cdot \left(z - a\right)\right)} \cdot \left(y - x\right)\right) \cdot a}{{t}^{2}}
\] |
associate-*r* [<=]54.7 | \[ \left(\frac{\left(y - x\right) \cdot \left(-1 \cdot z - -1 \cdot a\right)}{t} + y\right) + \frac{\color{blue}{\left(-1 \cdot \left(\left(z - a\right) \cdot \left(y - x\right)\right)\right)} \cdot a}{{t}^{2}}
\] |
*-commutative [<=]54.7 | \[ \left(\frac{\left(y - x\right) \cdot \left(-1 \cdot z - -1 \cdot a\right)}{t} + y\right) + \frac{\left(-1 \cdot \color{blue}{\left(\left(y - x\right) \cdot \left(z - a\right)\right)}\right) \cdot a}{{t}^{2}}
\] |
associate-*r* [<=]54.7 | \[ \left(\frac{\left(y - x\right) \cdot \left(-1 \cdot z - -1 \cdot a\right)}{t} + y\right) + \frac{\color{blue}{-1 \cdot \left(\left(\left(y - x\right) \cdot \left(z - a\right)\right) \cdot a\right)}}{{t}^{2}}
\] |
*-commutative [<=]54.7 | \[ \left(\frac{\left(y - x\right) \cdot \left(-1 \cdot z - -1 \cdot a\right)}{t} + y\right) + \frac{-1 \cdot \color{blue}{\left(a \cdot \left(\left(y - x\right) \cdot \left(z - a\right)\right)\right)}}{{t}^{2}}
\] |
associate-*r/ [<=]54.7 | \[ \left(\frac{\left(y - x\right) \cdot \left(-1 \cdot z - -1 \cdot a\right)}{t} + y\right) + \color{blue}{-1 \cdot \frac{a \cdot \left(\left(y - x\right) \cdot \left(z - a\right)\right)}{{t}^{2}}}
\] |
associate-+r+ [<=]54.7 | \[ \color{blue}{\frac{\left(y - x\right) \cdot \left(-1 \cdot z - -1 \cdot a\right)}{t} + \left(y + -1 \cdot \frac{a \cdot \left(\left(y - x\right) \cdot \left(z - a\right)\right)}{{t}^{2}}\right)}
\] |
if -2.20000000000000006e164 < t < 3.6000000000000001e142Initial program 76.9%
Simplified88.8%
[Start]76.9 | \[ x + \frac{\left(y - x\right) \cdot \left(z - t\right)}{a - t}
\] |
|---|---|
associate-/l* [=>]88.8 | \[ x + \color{blue}{\frac{y - x}{\frac{a - t}{z - t}}}
\] |
if 3.6000000000000001e142 < t Initial program 28.6%
Simplified58.8%
[Start]28.6 | \[ x + \frac{\left(y - x\right) \cdot \left(z - t\right)}{a - t}
\] |
|---|---|
+-commutative [=>]28.6 | \[ \color{blue}{\frac{\left(y - x\right) \cdot \left(z - t\right)}{a - t} + x}
\] |
associate-*l/ [<=]58.8 | \[ \color{blue}{\frac{y - x}{a - t} \cdot \left(z - t\right)} + x
\] |
fma-def [=>]58.8 | \[ \color{blue}{\mathsf{fma}\left(\frac{y - x}{a - t}, z - t, x\right)}
\] |
Taylor expanded in t around -inf 62.0%
Simplified84.9%
[Start]62.0 | \[ y + -1 \cdot \frac{-1 \cdot \left(a \cdot \left(y - x\right)\right) + z \cdot \left(y - x\right)}{t}
\] |
|---|---|
mul-1-neg [=>]62.0 | \[ y + \color{blue}{\left(-\frac{-1 \cdot \left(a \cdot \left(y - x\right)\right) + z \cdot \left(y - x\right)}{t}\right)}
\] |
unsub-neg [=>]62.0 | \[ \color{blue}{y - \frac{-1 \cdot \left(a \cdot \left(y - x\right)\right) + z \cdot \left(y - x\right)}{t}}
\] |
+-commutative [=>]62.0 | \[ y - \frac{\color{blue}{z \cdot \left(y - x\right) + -1 \cdot \left(a \cdot \left(y - x\right)\right)}}{t}
\] |
associate-*r* [=>]62.0 | \[ y - \frac{z \cdot \left(y - x\right) + \color{blue}{\left(-1 \cdot a\right) \cdot \left(y - x\right)}}{t}
\] |
distribute-rgt-out [=>]62.0 | \[ y - \frac{\color{blue}{\left(y - x\right) \cdot \left(z + -1 \cdot a\right)}}{t}
\] |
associate-/l* [=>]84.9 | \[ y - \color{blue}{\frac{y - x}{\frac{t}{z + -1 \cdot a}}}
\] |
mul-1-neg [=>]84.9 | \[ y - \frac{y - x}{\frac{t}{z + \color{blue}{\left(-a\right)}}}
\] |
Final simplification87.8%
| Alternative 1 | |
|---|---|
| Accuracy | 88.5% |
| Cost | 4432 |
| Alternative 2 | |
|---|---|
| Accuracy | 87.8% |
| Cost | 1348 |
| Alternative 3 | |
|---|---|
| Accuracy | 50.3% |
| Cost | 1240 |
| Alternative 4 | |
|---|---|
| Accuracy | 48.1% |
| Cost | 1240 |
| Alternative 5 | |
|---|---|
| Accuracy | 70.3% |
| Cost | 1233 |
| Alternative 6 | |
|---|---|
| Accuracy | 70.5% |
| Cost | 1232 |
| Alternative 7 | |
|---|---|
| Accuracy | 73.8% |
| Cost | 1232 |
| Alternative 8 | |
|---|---|
| Accuracy | 74.5% |
| Cost | 1232 |
| Alternative 9 | |
|---|---|
| Accuracy | 68.3% |
| Cost | 1105 |
| Alternative 10 | |
|---|---|
| Accuracy | 53.9% |
| Cost | 1104 |
| Alternative 11 | |
|---|---|
| Accuracy | 63.1% |
| Cost | 1104 |
| Alternative 12 | |
|---|---|
| Accuracy | 65.2% |
| Cost | 1104 |
| Alternative 13 | |
|---|---|
| Accuracy | 85.1% |
| Cost | 1097 |
| Alternative 14 | |
|---|---|
| Accuracy | 88.1% |
| Cost | 1097 |
| Alternative 15 | |
|---|---|
| Accuracy | 48.0% |
| Cost | 976 |
| Alternative 16 | |
|---|---|
| Accuracy | 47.9% |
| Cost | 976 |
| Alternative 17 | |
|---|---|
| Accuracy | 48.2% |
| Cost | 712 |
| Alternative 18 | |
|---|---|
| Accuracy | 43.8% |
| Cost | 328 |
| Alternative 19 | |
|---|---|
| Accuracy | 3.0% |
| Cost | 64 |
| Alternative 20 | |
|---|---|
| Accuracy | 28.5% |
| Cost | 64 |
herbie shell --seed 2023133
(FPCore (x y z t a)
:name "Graphics.Rendering.Chart.Axis.Types:linMap from Chart-1.5.3"
:precision binary64
:herbie-target
(if (< a -1.6153062845442575e-142) (+ x (* (/ (- y x) 1.0) (/ (- z t) (- a t)))) (if (< a 3.774403170083174e-182) (- y (* (/ z t) (- y x))) (+ x (* (/ (- y x) 1.0) (/ (- z t) (- a t))))))
(+ x (/ (* (- y x) (- z t)) (- a t))))