| Alternative 1 | |
|---|---|
| Error | 9.3 |
| Cost | 4880 |
(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) (- t x)) (- a z))))
(t_2
(+ (/ (* t (- y z)) (- a z)) (* x (+ 1.0 (- (/ (- y z) (- a z))))))))
(if (<= t_1 (- INFINITY))
(* t (- (/ y (- a z)) (/ z (- a z))))
(if (<= t_1 -1e-299)
t_2
(if (<= t_1 0.0)
(+ t (- (/ (* (- t x) (- y a)) z)))
(if (<= t_1 4e+290)
t_2
(- (+ t (* y (- (/ x z) (/ t z)))) (- (/ (* a (- t x)) 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) * (t - x)) / (a - z));
double t_2 = ((t * (y - z)) / (a - z)) + (x * (1.0 + -((y - z) / (a - z))));
double tmp;
if (t_1 <= -((double) INFINITY)) {
tmp = t * ((y / (a - z)) - (z / (a - z)));
} else if (t_1 <= -1e-299) {
tmp = t_2;
} else if (t_1 <= 0.0) {
tmp = t + -(((t - x) * (y - a)) / z);
} else if (t_1 <= 4e+290) {
tmp = t_2;
} else {
tmp = (t + (y * ((x / z) - (t / z)))) - -((a * (t - x)) / z);
}
return tmp;
}
public static double code(double x, double y, double z, double t, double a) {
return x + (((y - z) * (t - x)) / (a - z));
}
public static double code(double x, double y, double z, double t, double a) {
double t_1 = x + (((y - z) * (t - x)) / (a - z));
double t_2 = ((t * (y - z)) / (a - z)) + (x * (1.0 + -((y - z) / (a - z))));
double tmp;
if (t_1 <= -Double.POSITIVE_INFINITY) {
tmp = t * ((y / (a - z)) - (z / (a - z)));
} else if (t_1 <= -1e-299) {
tmp = t_2;
} else if (t_1 <= 0.0) {
tmp = t + -(((t - x) * (y - a)) / z);
} else if (t_1 <= 4e+290) {
tmp = t_2;
} else {
tmp = (t + (y * ((x / z) - (t / z)))) - -((a * (t - x)) / z);
}
return tmp;
}
def code(x, y, z, t, a): return x + (((y - z) * (t - x)) / (a - z))
def code(x, y, z, t, a): t_1 = x + (((y - z) * (t - x)) / (a - z)) t_2 = ((t * (y - z)) / (a - z)) + (x * (1.0 + -((y - z) / (a - z)))) tmp = 0 if t_1 <= -math.inf: tmp = t * ((y / (a - z)) - (z / (a - z))) elif t_1 <= -1e-299: tmp = t_2 elif t_1 <= 0.0: tmp = t + -(((t - x) * (y - a)) / z) elif t_1 <= 4e+290: tmp = t_2 else: tmp = (t + (y * ((x / z) - (t / z)))) - -((a * (t - x)) / 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(t - x)) / Float64(a - z))) t_2 = Float64(Float64(Float64(t * Float64(y - z)) / Float64(a - z)) + Float64(x * Float64(1.0 + Float64(-Float64(Float64(y - z) / Float64(a - z)))))) tmp = 0.0 if (t_1 <= Float64(-Inf)) tmp = Float64(t * Float64(Float64(y / Float64(a - z)) - Float64(z / Float64(a - z)))); elseif (t_1 <= -1e-299) tmp = t_2; elseif (t_1 <= 0.0) tmp = Float64(t + Float64(-Float64(Float64(Float64(t - x) * Float64(y - a)) / z))); elseif (t_1 <= 4e+290) tmp = t_2; else tmp = Float64(Float64(t + Float64(y * Float64(Float64(x / z) - Float64(t / z)))) - Float64(-Float64(Float64(a * Float64(t - x)) / z))); end return tmp end
function tmp = code(x, y, z, t, a) tmp = x + (((y - z) * (t - x)) / (a - z)); end
function tmp_2 = code(x, y, z, t, a) t_1 = x + (((y - z) * (t - x)) / (a - z)); t_2 = ((t * (y - z)) / (a - z)) + (x * (1.0 + -((y - z) / (a - z)))); tmp = 0.0; if (t_1 <= -Inf) tmp = t * ((y / (a - z)) - (z / (a - z))); elseif (t_1 <= -1e-299) tmp = t_2; elseif (t_1 <= 0.0) tmp = t + -(((t - x) * (y - a)) / z); elseif (t_1 <= 4e+290) tmp = t_2; else tmp = (t + (y * ((x / z) - (t / z)))) - -((a * (t - x)) / z); end tmp_2 = 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[(t - x), $MachinePrecision]), $MachinePrecision] / N[(a - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(N[(N[(t * N[(y - z), $MachinePrecision]), $MachinePrecision] / N[(a - z), $MachinePrecision]), $MachinePrecision] + N[(x * N[(1.0 + (-N[(N[(y - z), $MachinePrecision] / N[(a - z), $MachinePrecision]), $MachinePrecision])), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, (-Infinity)], N[(t * N[(N[(y / N[(a - z), $MachinePrecision]), $MachinePrecision] - N[(z / N[(a - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$1, -1e-299], t$95$2, If[LessEqual[t$95$1, 0.0], N[(t + (-N[(N[(N[(t - x), $MachinePrecision] * N[(y - a), $MachinePrecision]), $MachinePrecision] / z), $MachinePrecision])), $MachinePrecision], If[LessEqual[t$95$1, 4e+290], t$95$2, N[(N[(t + N[(y * N[(N[(x / z), $MachinePrecision] - N[(t / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - (-N[(N[(a * N[(t - x), $MachinePrecision]), $MachinePrecision] / z), $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(t - x\right)}{a - z}\\
t_2 := \frac{t \cdot \left(y - z\right)}{a - z} + x \cdot \left(1 + \left(-\frac{y - z}{a - z}\right)\right)\\
\mathbf{if}\;t_1 \leq -\infty:\\
\;\;\;\;t \cdot \left(\frac{y}{a - z} - \frac{z}{a - z}\right)\\
\mathbf{elif}\;t_1 \leq -1 \cdot 10^{-299}:\\
\;\;\;\;t_2\\
\mathbf{elif}\;t_1 \leq 0:\\
\;\;\;\;t + \left(-\frac{\left(t - x\right) \cdot \left(y - a\right)}{z}\right)\\
\mathbf{elif}\;t_1 \leq 4 \cdot 10^{+290}:\\
\;\;\;\;t_2\\
\mathbf{else}:\\
\;\;\;\;\left(t + y \cdot \left(\frac{x}{z} - \frac{t}{z}\right)\right) - \left(-\frac{a \cdot \left(t - x\right)}{z}\right)\\
\end{array}
Results
| Original | 24.1 |
|---|---|
| Target | 11.8 |
| Herbie | 9.4 |
if (+.f64 x (/.f64 (*.f64 (-.f64 y z) (-.f64 t x)) (-.f64 a z))) < -inf.0Initial program 64.0
Taylor expanded in t around inf 25.9
if -inf.0 < (+.f64 x (/.f64 (*.f64 (-.f64 y z) (-.f64 t x)) (-.f64 a z))) < -9.99999999999999992e-300 or 0.0 < (+.f64 x (/.f64 (*.f64 (-.f64 y z) (-.f64 t x)) (-.f64 a z))) < 4.00000000000000025e290Initial program 2.0
Taylor expanded in x around 0 2.1
Simplified2.1
[Start]2.1 | \[ \left(1 + -1 \cdot \frac{y - z}{a - z}\right) \cdot x + \frac{t \cdot \left(y - z\right)}{a - z}
\] |
|---|---|
rational.json-simplify-1 [=>]2.1 | \[ \color{blue}{\frac{t \cdot \left(y - z\right)}{a - z} + \left(1 + -1 \cdot \frac{y - z}{a - z}\right) \cdot x}
\] |
rational.json-simplify-2 [=>]2.1 | \[ \frac{t \cdot \left(y - z\right)}{a - z} + \color{blue}{x \cdot \left(1 + -1 \cdot \frac{y - z}{a - z}\right)}
\] |
rational.json-simplify-2 [=>]2.1 | \[ \frac{t \cdot \left(y - z\right)}{a - z} + x \cdot \left(1 + \color{blue}{\frac{y - z}{a - z} \cdot -1}\right)
\] |
rational.json-simplify-9 [=>]2.1 | \[ \frac{t \cdot \left(y - z\right)}{a - z} + x \cdot \left(1 + \color{blue}{\left(-\frac{y - z}{a - z}\right)}\right)
\] |
if -9.99999999999999992e-300 < (+.f64 x (/.f64 (*.f64 (-.f64 y z) (-.f64 t x)) (-.f64 a z))) < 0.0Initial program 60.0
Taylor expanded in z around -inf 1.2
Simplified1.2
[Start]1.2 | \[ -1 \cdot \frac{y \cdot \left(t - x\right) - a \cdot \left(t - x\right)}{z} + t
\] |
|---|---|
rational.json-simplify-1 [=>]1.2 | \[ \color{blue}{t + -1 \cdot \frac{y \cdot \left(t - x\right) - a \cdot \left(t - x\right)}{z}}
\] |
rational.json-simplify-2 [=>]1.2 | \[ t + \color{blue}{\frac{y \cdot \left(t - x\right) - a \cdot \left(t - x\right)}{z} \cdot -1}
\] |
rational.json-simplify-9 [=>]1.2 | \[ t + \color{blue}{\left(-\frac{y \cdot \left(t - x\right) - a \cdot \left(t - x\right)}{z}\right)}
\] |
rational.json-simplify-2 [=>]1.2 | \[ t + \left(-\frac{\color{blue}{\left(t - x\right) \cdot y} - a \cdot \left(t - x\right)}{z}\right)
\] |
rational.json-simplify-48 [=>]1.2 | \[ t + \left(-\frac{\color{blue}{\left(t - x\right) \cdot \left(y - a\right)}}{z}\right)
\] |
if 4.00000000000000025e290 < (+.f64 x (/.f64 (*.f64 (-.f64 y z) (-.f64 t x)) (-.f64 a z))) Initial program 61.8
Taylor expanded in z around -inf 39.7
Simplified39.7
[Start]39.7 | \[ -1 \cdot \frac{y \cdot \left(t - x\right) - a \cdot \left(t - x\right)}{z} + t
\] |
|---|---|
rational.json-simplify-1 [=>]39.7 | \[ \color{blue}{t + -1 \cdot \frac{y \cdot \left(t - x\right) - a \cdot \left(t - x\right)}{z}}
\] |
rational.json-simplify-2 [=>]39.7 | \[ t + \color{blue}{\frac{y \cdot \left(t - x\right) - a \cdot \left(t - x\right)}{z} \cdot -1}
\] |
rational.json-simplify-9 [=>]39.7 | \[ t + \color{blue}{\left(-\frac{y \cdot \left(t - x\right) - a \cdot \left(t - x\right)}{z}\right)}
\] |
rational.json-simplify-2 [=>]39.7 | \[ t + \left(-\frac{\color{blue}{\left(t - x\right) \cdot y} - a \cdot \left(t - x\right)}{z}\right)
\] |
rational.json-simplify-48 [=>]39.7 | \[ t + \left(-\frac{\color{blue}{\left(t - x\right) \cdot \left(y - a\right)}}{z}\right)
\] |
Taylor expanded in y around 0 29.8
Simplified29.8
[Start]29.8 | \[ \left(y \cdot \left(\frac{x}{z} - \frac{t}{z}\right) + t\right) - -1 \cdot \frac{a \cdot \left(t - x\right)}{z}
\] |
|---|---|
rational.json-simplify-1 [=>]29.8 | \[ \color{blue}{\left(t + y \cdot \left(\frac{x}{z} - \frac{t}{z}\right)\right)} - -1 \cdot \frac{a \cdot \left(t - x\right)}{z}
\] |
rational.json-simplify-2 [=>]29.8 | \[ \left(t + y \cdot \left(\frac{x}{z} - \frac{t}{z}\right)\right) - \color{blue}{\frac{a \cdot \left(t - x\right)}{z} \cdot -1}
\] |
rational.json-simplify-9 [=>]29.8 | \[ \left(t + y \cdot \left(\frac{x}{z} - \frac{t}{z}\right)\right) - \color{blue}{\left(-\frac{a \cdot \left(t - x\right)}{z}\right)}
\] |
Final simplification9.4
| Alternative 1 | |
|---|---|
| Error | 9.3 |
| Cost | 4880 |
| Alternative 2 | |
|---|---|
| Error | 8.8 |
| Cost | 4432 |
| Alternative 3 | |
|---|---|
| Error | 23.2 |
| Cost | 2020 |
| Alternative 4 | |
|---|---|
| Error | 34.3 |
| Cost | 1900 |
| Alternative 5 | |
|---|---|
| Error | 33.0 |
| Cost | 1632 |
| Alternative 6 | |
|---|---|
| Error | 33.7 |
| Cost | 1372 |
| Alternative 7 | |
|---|---|
| Error | 33.9 |
| Cost | 1372 |
| Alternative 8 | |
|---|---|
| Error | 26.0 |
| Cost | 1300 |
| Alternative 9 | |
|---|---|
| Error | 25.9 |
| Cost | 1300 |
| Alternative 10 | |
|---|---|
| Error | 22.2 |
| Cost | 1296 |
| Alternative 11 | |
|---|---|
| Error | 31.2 |
| Cost | 1240 |
| Alternative 12 | |
|---|---|
| Error | 26.4 |
| Cost | 1172 |
| Alternative 13 | |
|---|---|
| Error | 24.4 |
| Cost | 1168 |
| Alternative 14 | |
|---|---|
| Error | 30.8 |
| Cost | 1108 |
| Alternative 15 | |
|---|---|
| Error | 37.7 |
| Cost | 976 |
| Alternative 16 | |
|---|---|
| Error | 37.6 |
| Cost | 976 |
| Alternative 17 | |
|---|---|
| Error | 35.9 |
| Cost | 844 |
| Alternative 18 | |
|---|---|
| Error | 32.9 |
| Cost | 844 |
| Alternative 19 | |
|---|---|
| Error | 27.2 |
| Cost | 840 |
| Alternative 20 | |
|---|---|
| Error | 35.9 |
| Cost | 712 |
| Alternative 21 | |
|---|---|
| Error | 35.8 |
| Cost | 328 |
| Alternative 22 | |
|---|---|
| Error | 45.9 |
| Cost | 64 |
herbie shell --seed 2023077
(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))))