| Alternative 1 | |
|---|---|
| Accuracy | 97.6% |
| Cost | 20160 |
\[\mathsf{fma}\left(a, 1 - t, \mathsf{fma}\left(y + \left(t + -2\right), b, \mathsf{fma}\left(z, 1 - y, x\right)\right)\right)
\]

(FPCore (x y z t a b) :precision binary64 (+ (- (- x (* (- y 1.0) z)) (* (- t 1.0) a)) (* (- (+ y t) 2.0) b)))
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1
(+ (+ (+ x (* z (- 1.0 y))) (* a (- 1.0 t))) (* b (- (+ t y) 2.0)))))
(if (<= t_1 INFINITY) t_1 (* y (- b z)))))double code(double x, double y, double z, double t, double a, double b) {
return ((x - ((y - 1.0) * z)) - ((t - 1.0) * a)) + (((y + t) - 2.0) * b);
}
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = ((x + (z * (1.0 - y))) + (a * (1.0 - t))) + (b * ((t + y) - 2.0));
double tmp;
if (t_1 <= ((double) INFINITY)) {
tmp = t_1;
} else {
tmp = y * (b - z);
}
return tmp;
}
public static double code(double x, double y, double z, double t, double a, double b) {
return ((x - ((y - 1.0) * z)) - ((t - 1.0) * a)) + (((y + t) - 2.0) * b);
}
public static double code(double x, double y, double z, double t, double a, double b) {
double t_1 = ((x + (z * (1.0 - y))) + (a * (1.0 - t))) + (b * ((t + y) - 2.0));
double tmp;
if (t_1 <= Double.POSITIVE_INFINITY) {
tmp = t_1;
} else {
tmp = y * (b - z);
}
return tmp;
}
def code(x, y, z, t, a, b): return ((x - ((y - 1.0) * z)) - ((t - 1.0) * a)) + (((y + t) - 2.0) * b)
def code(x, y, z, t, a, b): t_1 = ((x + (z * (1.0 - y))) + (a * (1.0 - t))) + (b * ((t + y) - 2.0)) tmp = 0 if t_1 <= math.inf: tmp = t_1 else: tmp = y * (b - z) return tmp
function code(x, y, z, t, a, b) return Float64(Float64(Float64(x - Float64(Float64(y - 1.0) * z)) - Float64(Float64(t - 1.0) * a)) + Float64(Float64(Float64(y + t) - 2.0) * b)) end
function code(x, y, z, t, a, b) t_1 = Float64(Float64(Float64(x + Float64(z * Float64(1.0 - y))) + Float64(a * Float64(1.0 - t))) + Float64(b * Float64(Float64(t + y) - 2.0))) tmp = 0.0 if (t_1 <= Inf) tmp = t_1; else tmp = Float64(y * Float64(b - z)); end return tmp end
function tmp = code(x, y, z, t, a, b) tmp = ((x - ((y - 1.0) * z)) - ((t - 1.0) * a)) + (((y + t) - 2.0) * b); end
function tmp_2 = code(x, y, z, t, a, b) t_1 = ((x + (z * (1.0 - y))) + (a * (1.0 - t))) + (b * ((t + y) - 2.0)); tmp = 0.0; if (t_1 <= Inf) tmp = t_1; else tmp = y * (b - z); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := N[(N[(N[(x - N[(N[(y - 1.0), $MachinePrecision] * z), $MachinePrecision]), $MachinePrecision] - N[(N[(t - 1.0), $MachinePrecision] * a), $MachinePrecision]), $MachinePrecision] + N[(N[(N[(y + t), $MachinePrecision] - 2.0), $MachinePrecision] * b), $MachinePrecision]), $MachinePrecision]
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(N[(x + N[(z * N[(1.0 - y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(a * N[(1.0 - t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(b * N[(N[(t + y), $MachinePrecision] - 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, Infinity], t$95$1, N[(y * N[(b - z), $MachinePrecision]), $MachinePrecision]]]
\left(\left(x - \left(y - 1\right) \cdot z\right) - \left(t - 1\right) \cdot a\right) + \left(\left(y + t\right) - 2\right) \cdot b
\begin{array}{l}
t_1 := \left(\left(x + z \cdot \left(1 - y\right)\right) + a \cdot \left(1 - t\right)\right) + b \cdot \left(\left(t + y\right) - 2\right)\\
\mathbf{if}\;t_1 \leq \infty:\\
\;\;\;\;t_1\\
\mathbf{else}:\\
\;\;\;\;y \cdot \left(b - z\right)\\
\end{array}
Herbie found 27 alternatives:
| Alternative | Accuracy | Speedup |
|---|
Results
if (+.f64 (-.f64 (-.f64 x (*.f64 (-.f64 y 1) z)) (*.f64 (-.f64 t 1) a)) (*.f64 (-.f64 (+.f64 y t) 2) b)) < +inf.0Initial program 100.0%
if +inf.0 < (+.f64 (-.f64 (-.f64 x (*.f64 (-.f64 y 1) z)) (*.f64 (-.f64 t 1) a)) (*.f64 (-.f64 (+.f64 y t) 2) b)) Initial program 0.0%
Simplified0.0%
[Start]0.0% | \[ \left(\left(x - \left(y - 1\right) \cdot z\right) - \left(t - 1\right) \cdot a\right) + \left(\left(y + t\right) - 2\right) \cdot b
\] |
|---|---|
associate-+l- [=>]0.0% | \[ \color{blue}{\left(x - \left(y - 1\right) \cdot z\right) - \left(\left(t - 1\right) \cdot a - \left(\left(y + t\right) - 2\right) \cdot b\right)}
\] |
*-commutative [=>]0.0% | \[ \left(x - \color{blue}{z \cdot \left(y - 1\right)}\right) - \left(\left(t - 1\right) \cdot a - \left(\left(y + t\right) - 2\right) \cdot b\right)
\] |
*-commutative [<=]0.0% | \[ \left(x - \color{blue}{\left(y - 1\right) \cdot z}\right) - \left(\left(t - 1\right) \cdot a - \left(\left(y + t\right) - 2\right) \cdot b\right)
\] |
sub-neg [=>]0.0% | \[ \left(x - \color{blue}{\left(y + \left(-1\right)\right)} \cdot z\right) - \left(\left(t - 1\right) \cdot a - \left(\left(y + t\right) - 2\right) \cdot b\right)
\] |
metadata-eval [=>]0.0% | \[ \left(x - \left(y + \color{blue}{-1}\right) \cdot z\right) - \left(\left(t - 1\right) \cdot a - \left(\left(y + t\right) - 2\right) \cdot b\right)
\] |
remove-double-neg [<=]0.0% | \[ \left(x - \left(y + -1\right) \cdot z\right) - \left(\color{blue}{\left(-\left(-\left(t - 1\right)\right)\right)} \cdot a - \left(\left(y + t\right) - 2\right) \cdot b\right)
\] |
remove-double-neg [=>]0.0% | \[ \left(x - \left(y + -1\right) \cdot z\right) - \left(\color{blue}{\left(t - 1\right)} \cdot a - \left(\left(y + t\right) - 2\right) \cdot b\right)
\] |
sub-neg [=>]0.0% | \[ \left(x - \left(y + -1\right) \cdot z\right) - \left(\color{blue}{\left(t + \left(-1\right)\right)} \cdot a - \left(\left(y + t\right) - 2\right) \cdot b\right)
\] |
metadata-eval [=>]0.0% | \[ \left(x - \left(y + -1\right) \cdot z\right) - \left(\left(t + \color{blue}{-1}\right) \cdot a - \left(\left(y + t\right) - 2\right) \cdot b\right)
\] |
associate--l+ [=>]0.0% | \[ \left(x - \left(y + -1\right) \cdot z\right) - \left(\left(t + -1\right) \cdot a - \color{blue}{\left(y + \left(t - 2\right)\right)} \cdot b\right)
\] |
Taylor expanded in y around inf 86.1%
Final simplification99.6%
| Alternative 1 | |
|---|---|
| Accuracy | 97.6% |
| Cost | 20160 |
| Alternative 2 | |
|---|---|
| Accuracy | 97.4% |
| Cost | 13888 |
| Alternative 3 | |
|---|---|
| Accuracy | 79.5% |
| Cost | 1756 |
| Alternative 4 | |
|---|---|
| Accuracy | 59.7% |
| Cost | 1505 |
| Alternative 5 | |
|---|---|
| Accuracy | 71.5% |
| Cost | 1496 |
| Alternative 6 | |
|---|---|
| Accuracy | 71.5% |
| Cost | 1364 |
| Alternative 7 | |
|---|---|
| Accuracy | 72.3% |
| Cost | 1360 |
| Alternative 8 | |
|---|---|
| Accuracy | 59.2% |
| Cost | 1241 |
| Alternative 9 | |
|---|---|
| Accuracy | 57.1% |
| Cost | 1240 |
| Alternative 10 | |
|---|---|
| Accuracy | 70.3% |
| Cost | 1236 |
| Alternative 11 | |
|---|---|
| Accuracy | 85.3% |
| Cost | 1225 |
| Alternative 12 | |
|---|---|
| Accuracy | 68.5% |
| Cost | 1104 |
| Alternative 13 | |
|---|---|
| Accuracy | 35.7% |
| Cost | 1049 |
| Alternative 14 | |
|---|---|
| Accuracy | 60.6% |
| Cost | 977 |
| Alternative 15 | |
|---|---|
| Accuracy | 26.6% |
| Cost | 852 |
| Alternative 16 | |
|---|---|
| Accuracy | 47.0% |
| Cost | 848 |
| Alternative 17 | |
|---|---|
| Accuracy | 45.6% |
| Cost | 848 |
| Alternative 18 | |
|---|---|
| Accuracy | 50.3% |
| Cost | 848 |
| Alternative 19 | |
|---|---|
| Accuracy | 50.5% |
| Cost | 848 |
| Alternative 20 | |
|---|---|
| Accuracy | 69.1% |
| Cost | 841 |
| Alternative 21 | |
|---|---|
| Accuracy | 36.0% |
| Cost | 785 |
| Alternative 22 | |
|---|---|
| Accuracy | 34.4% |
| Cost | 720 |
| Alternative 23 | |
|---|---|
| Accuracy | 23.2% |
| Cost | 460 |
| Alternative 24 | |
|---|---|
| Accuracy | 21.3% |
| Cost | 328 |
| Alternative 25 | |
|---|---|
| Accuracy | 21.2% |
| Cost | 328 |
| Alternative 26 | |
|---|---|
| Accuracy | 11.3% |
| Cost | 64 |
herbie shell --seed 2023164
(FPCore (x y z t a b)
:name "Statistics.Distribution.Beta:$centropy from math-functions-0.1.5.2"
:precision binary64
(+ (- (- x (* (- y 1.0) z)) (* (- t 1.0) a)) (* (- (+ y t) 2.0) b)))