| Alternative 1 | |
|---|---|
| Accuracy | 96.1% |
| Cost | 10704 |

(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 (fma (- t x) (/ (- y z) (- a z)) x))
(t_2 (+ x (* (- z y) (/ (- x t) (- a z))))))
(if (<= t_2 -1e-293)
t_1
(if (<= t_2 0.0)
(+ t (* (/ x z) (- y a)))
(if (<= t_2 0.5)
(+
(/ (* (- y z) t) (- a z))
(* x (- (+ (/ z (- a z)) 1.0) (/ y (- a z)))))
(if (<= t_2 4e+286) (fma (- y z) (/ (- t x) (- a z)) x) t_1))))))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 = fma((t - x), ((y - z) / (a - z)), x);
double t_2 = x + ((z - y) * ((x - t) / (a - z)));
double tmp;
if (t_2 <= -1e-293) {
tmp = t_1;
} else if (t_2 <= 0.0) {
tmp = t + ((x / z) * (y - a));
} else if (t_2 <= 0.5) {
tmp = (((y - z) * t) / (a - z)) + (x * (((z / (a - z)) + 1.0) - (y / (a - z))));
} else if (t_2 <= 4e+286) {
tmp = fma((y - z), ((t - x) / (a - z)), x);
} else {
tmp = t_1;
}
return tmp;
}
function code(x, y, z, t, a) return Float64(x + Float64(Float64(y - z) * Float64(Float64(t - x) / Float64(a - z)))) end
function code(x, y, z, t, a) t_1 = fma(Float64(t - x), Float64(Float64(y - z) / Float64(a - z)), x) t_2 = Float64(x + Float64(Float64(z - y) * Float64(Float64(x - t) / Float64(a - z)))) tmp = 0.0 if (t_2 <= -1e-293) tmp = t_1; elseif (t_2 <= 0.0) tmp = Float64(t + Float64(Float64(x / z) * Float64(y - a))); elseif (t_2 <= 0.5) tmp = Float64(Float64(Float64(Float64(y - z) * t) / Float64(a - z)) + Float64(x * Float64(Float64(Float64(z / Float64(a - z)) + 1.0) - Float64(y / Float64(a - z))))); elseif (t_2 <= 4e+286) tmp = fma(Float64(y - z), Float64(Float64(t - x) / Float64(a - z)), x); else tmp = t_1; end return tmp end
code[x_, y_, z_, t_, a_] := N[(x + N[(N[(y - z), $MachinePrecision] * N[(N[(t - x), $MachinePrecision] / N[(a - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(N[(t - x), $MachinePrecision] * N[(N[(y - z), $MachinePrecision] / N[(a - z), $MachinePrecision]), $MachinePrecision] + x), $MachinePrecision]}, Block[{t$95$2 = N[(x + N[(N[(z - y), $MachinePrecision] * N[(N[(x - t), $MachinePrecision] / N[(a - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$2, -1e-293], t$95$1, If[LessEqual[t$95$2, 0.0], N[(t + N[(N[(x / z), $MachinePrecision] * N[(y - a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$2, 0.5], N[(N[(N[(N[(y - z), $MachinePrecision] * t), $MachinePrecision] / N[(a - z), $MachinePrecision]), $MachinePrecision] + N[(x * N[(N[(N[(z / N[(a - z), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision] - N[(y / N[(a - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$2, 4e+286], N[(N[(y - z), $MachinePrecision] * N[(N[(t - x), $MachinePrecision] / N[(a - z), $MachinePrecision]), $MachinePrecision] + x), $MachinePrecision], t$95$1]]]]]]
x + \left(y - z\right) \cdot \frac{t - x}{a - z}
\begin{array}{l}
t_1 := \mathsf{fma}\left(t - x, \frac{y - z}{a - z}, x\right)\\
t_2 := x + \left(z - y\right) \cdot \frac{x - t}{a - z}\\
\mathbf{if}\;t_2 \leq -1 \cdot 10^{-293}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;t_2 \leq 0:\\
\;\;\;\;t + \frac{x}{z} \cdot \left(y - a\right)\\
\mathbf{elif}\;t_2 \leq 0.5:\\
\;\;\;\;\frac{\left(y - z\right) \cdot t}{a - z} + x \cdot \left(\left(\frac{z}{a - z} + 1\right) - \frac{y}{a - z}\right)\\
\mathbf{elif}\;t_2 \leq 4 \cdot 10^{+286}:\\
\;\;\;\;\mathsf{fma}\left(y - z, \frac{t - x}{a - z}, x\right)\\
\mathbf{else}:\\
\;\;\;\;t_1\\
\end{array}
Herbie found 25 alternatives:
| Alternative | Accuracy | Speedup |
|---|
if (+.f64 x (*.f64 (-.f64 y z) (/.f64 (-.f64 t x) (-.f64 a z)))) < -1.0000000000000001e-293 or 4.00000000000000013e286 < (+.f64 x (*.f64 (-.f64 y z) (/.f64 (-.f64 t x) (-.f64 a z)))) Initial program 90.5%
Simplified96.3%
[Start]90.5% | \[ x + \left(y - z\right) \cdot \frac{t - x}{a - z}
\] |
|---|---|
+-commutative [=>]90.5% | \[ \color{blue}{\left(y - z\right) \cdot \frac{t - x}{a - z} + x}
\] |
associate-*r/ [=>]76.4% | \[ \color{blue}{\frac{\left(y - z\right) \cdot \left(t - x\right)}{a - z}} + x
\] |
*-commutative [=>]76.4% | \[ \frac{\color{blue}{\left(t - x\right) \cdot \left(y - z\right)}}{a - z} + x
\] |
associate-*r/ [<=]96.2% | \[ \color{blue}{\left(t - x\right) \cdot \frac{y - z}{a - z}} + x
\] |
fma-def [=>]96.3% | \[ \color{blue}{\mathsf{fma}\left(t - x, \frac{y - z}{a - z}, x\right)}
\] |
if -1.0000000000000001e-293 < (+.f64 x (*.f64 (-.f64 y z) (/.f64 (-.f64 t x) (-.f64 a z)))) < 0.0Initial program 3.4%
Taylor expanded in z around inf 81.4%
Simplified81.4%
[Start]81.4% | \[ \left(-1 \cdot \frac{y \cdot \left(t - x\right)}{z} + t\right) - -1 \cdot \frac{a \cdot \left(t - x\right)}{z}
\] |
|---|---|
+-commutative [=>]81.4% | \[ \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+ [=>]81.4% | \[ \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/ [=>]81.4% | \[ 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/ [=>]81.4% | \[ 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 [<=]81.4% | \[ 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-- [=>]81.4% | \[ t + \frac{\color{blue}{-1 \cdot \left(y \cdot \left(t - x\right) - a \cdot \left(t - x\right)\right)}}{z}
\] |
mul-1-neg [=>]81.4% | \[ t + \frac{\color{blue}{-\left(y \cdot \left(t - x\right) - a \cdot \left(t - x\right)\right)}}{z}
\] |
distribute-neg-frac [<=]81.4% | \[ t + \color{blue}{\left(-\frac{y \cdot \left(t - x\right) - a \cdot \left(t - x\right)}{z}\right)}
\] |
unsub-neg [=>]81.4% | \[ \color{blue}{t - \frac{y \cdot \left(t - x\right) - a \cdot \left(t - x\right)}{z}}
\] |
distribute-rgt-out-- [=>]81.4% | \[ t - \frac{\color{blue}{\left(t - x\right) \cdot \left(y - a\right)}}{z}
\] |
Taylor expanded in t around 0 81.4%
Simplified99.8%
[Start]81.4% | \[ t - -1 \cdot \frac{\left(y - a\right) \cdot x}{z}
\] |
|---|---|
associate-*r/ [=>]81.4% | \[ t - \color{blue}{\frac{-1 \cdot \left(\left(y - a\right) \cdot x\right)}{z}}
\] |
mul-1-neg [=>]81.4% | \[ t - \frac{\color{blue}{-\left(y - a\right) \cdot x}}{z}
\] |
*-commutative [=>]81.4% | \[ t - \frac{-\color{blue}{x \cdot \left(y - a\right)}}{z}
\] |
distribute-rgt-neg-in [=>]81.4% | \[ t - \frac{\color{blue}{x \cdot \left(-\left(y - a\right)\right)}}{z}
\] |
neg-sub0 [=>]81.4% | \[ t - \frac{x \cdot \color{blue}{\left(0 - \left(y - a\right)\right)}}{z}
\] |
associate--r- [=>]81.4% | \[ t - \frac{x \cdot \color{blue}{\left(\left(0 - y\right) + a\right)}}{z}
\] |
neg-sub0 [<=]81.4% | \[ t - \frac{x \cdot \left(\color{blue}{\left(-y\right)} + a\right)}{z}
\] |
distribute-rgt-out [<=]81.4% | \[ t - \frac{\color{blue}{\left(-y\right) \cdot x + a \cdot x}}{z}
\] |
+-commutative [<=]81.4% | \[ t - \frac{\color{blue}{a \cdot x + \left(-y\right) \cdot x}}{z}
\] |
distribute-rgt-out [=>]81.4% | \[ t - \frac{\color{blue}{x \cdot \left(a + \left(-y\right)\right)}}{z}
\] |
sub-neg [<=]81.4% | \[ t - \frac{x \cdot \color{blue}{\left(a - y\right)}}{z}
\] |
*-commutative [=>]81.4% | \[ t - \frac{\color{blue}{\left(a - y\right) \cdot x}}{z}
\] |
associate-*r/ [<=]99.8% | \[ t - \color{blue}{\left(a - y\right) \cdot \frac{x}{z}}
\] |
if 0.0 < (+.f64 x (*.f64 (-.f64 y z) (/.f64 (-.f64 t x) (-.f64 a z)))) < 0.5Initial program 77.9%
Taylor expanded in x around -inf 99.5%
if 0.5 < (+.f64 x (*.f64 (-.f64 y z) (/.f64 (-.f64 t x) (-.f64 a z)))) < 4.00000000000000013e286Initial program 97.9%
Simplified97.9%
[Start]97.9% | \[ x + \left(y - z\right) \cdot \frac{t - x}{a - z}
\] |
|---|---|
+-commutative [=>]97.9% | \[ \color{blue}{\left(y - z\right) \cdot \frac{t - x}{a - z} + x}
\] |
fma-def [=>]97.9% | \[ \color{blue}{\mathsf{fma}\left(y - z, \frac{t - x}{a - z}, x\right)}
\] |
Final simplification97.3%
| Alternative 1 | |
|---|---|
| Accuracy | 96.1% |
| Cost | 10704 |
| Alternative 2 | |
|---|---|
| Accuracy | 96.0% |
| Cost | 10704 |
| Alternative 3 | |
|---|---|
| Accuracy | 93.2% |
| Cost | 4300 |
| Alternative 4 | |
|---|---|
| Accuracy | 91.1% |
| Cost | 2633 |
| Alternative 5 | |
|---|---|
| Accuracy | 67.7% |
| Cost | 1632 |
| Alternative 6 | |
|---|---|
| Accuracy | 50.9% |
| Cost | 1568 |
| Alternative 7 | |
|---|---|
| Accuracy | 49.9% |
| Cost | 1504 |
| Alternative 8 | |
|---|---|
| Accuracy | 49.1% |
| Cost | 1504 |
| Alternative 9 | |
|---|---|
| Accuracy | 70.5% |
| Cost | 1496 |
| Alternative 10 | |
|---|---|
| Accuracy | 57.8% |
| Cost | 1368 |
| Alternative 11 | |
|---|---|
| Accuracy | 69.2% |
| Cost | 1368 |
| Alternative 12 | |
|---|---|
| Accuracy | 69.1% |
| Cost | 1368 |
| Alternative 13 | |
|---|---|
| Accuracy | 52.2% |
| Cost | 1304 |
| Alternative 14 | |
|---|---|
| Accuracy | 48.5% |
| Cost | 1240 |
| Alternative 15 | |
|---|---|
| Accuracy | 73.3% |
| Cost | 1232 |
| Alternative 16 | |
|---|---|
| Accuracy | 69.5% |
| Cost | 1232 |
| Alternative 17 | |
|---|---|
| Accuracy | 67.8% |
| Cost | 1104 |
| Alternative 18 | |
|---|---|
| Accuracy | 36.5% |
| Cost | 976 |
| Alternative 19 | |
|---|---|
| Accuracy | 37.0% |
| Cost | 976 |
| Alternative 20 | |
|---|---|
| Accuracy | 37.6% |
| Cost | 848 |
| Alternative 21 | |
|---|---|
| Accuracy | 61.3% |
| Cost | 841 |
| Alternative 22 | |
|---|---|
| Accuracy | 61.3% |
| Cost | 840 |
| Alternative 23 | |
|---|---|
| Accuracy | 36.5% |
| Cost | 716 |
| Alternative 24 | |
|---|---|
| Accuracy | 38.6% |
| Cost | 328 |
| Alternative 25 | |
|---|---|
| Accuracy | 25.4% |
| Cost | 64 |
herbie shell --seed 2023171
(FPCore (x y z t a)
:name "Numeric.Signal:interpolate from hsignal-0.2.7.1"
:precision binary64
(+ x (* (- y z) (/ (- t x) (- a z)))))