| Alternative 1 | |
|---|---|
| Accuracy | 98.7% |
| Cost | 1609 |
(FPCore (x y z t) :precision binary64 (/ x (* (- y z) (- t z))))
(FPCore (x y z t)
:precision binary64
(let* ((t_1 (* (- y z) (- t z))))
(if (<= t_1 -5e+296)
(/ (/ 1.0 (- z y)) (/ (- z t) x))
(if (<= t_1 5e+154) (/ x t_1) (/ (/ x (- z t)) (- z y))))))double code(double x, double y, double z, double t) {
return x / ((y - z) * (t - z));
}
double code(double x, double y, double z, double t) {
double t_1 = (y - z) * (t - z);
double tmp;
if (t_1 <= -5e+296) {
tmp = (1.0 / (z - y)) / ((z - t) / x);
} else if (t_1 <= 5e+154) {
tmp = x / t_1;
} else {
tmp = (x / (z - t)) / (z - y);
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
code = x / ((y - z) * (t - z))
end function
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: t_1
real(8) :: tmp
t_1 = (y - z) * (t - z)
if (t_1 <= (-5d+296)) then
tmp = (1.0d0 / (z - y)) / ((z - t) / x)
else if (t_1 <= 5d+154) then
tmp = x / t_1
else
tmp = (x / (z - t)) / (z - y)
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
return x / ((y - z) * (t - z));
}
public static double code(double x, double y, double z, double t) {
double t_1 = (y - z) * (t - z);
double tmp;
if (t_1 <= -5e+296) {
tmp = (1.0 / (z - y)) / ((z - t) / x);
} else if (t_1 <= 5e+154) {
tmp = x / t_1;
} else {
tmp = (x / (z - t)) / (z - y);
}
return tmp;
}
def code(x, y, z, t): return x / ((y - z) * (t - z))
def code(x, y, z, t): t_1 = (y - z) * (t - z) tmp = 0 if t_1 <= -5e+296: tmp = (1.0 / (z - y)) / ((z - t) / x) elif t_1 <= 5e+154: tmp = x / t_1 else: tmp = (x / (z - t)) / (z - y) return tmp
function code(x, y, z, t) return Float64(x / Float64(Float64(y - z) * Float64(t - z))) end
function code(x, y, z, t) t_1 = Float64(Float64(y - z) * Float64(t - z)) tmp = 0.0 if (t_1 <= -5e+296) tmp = Float64(Float64(1.0 / Float64(z - y)) / Float64(Float64(z - t) / x)); elseif (t_1 <= 5e+154) tmp = Float64(x / t_1); else tmp = Float64(Float64(x / Float64(z - t)) / Float64(z - y)); end return tmp end
function tmp = code(x, y, z, t) tmp = x / ((y - z) * (t - z)); end
function tmp_2 = code(x, y, z, t) t_1 = (y - z) * (t - z); tmp = 0.0; if (t_1 <= -5e+296) tmp = (1.0 / (z - y)) / ((z - t) / x); elseif (t_1 <= 5e+154) tmp = x / t_1; else tmp = (x / (z - t)) / (z - y); end tmp_2 = tmp; end
code[x_, y_, z_, t_] := N[(x / N[(N[(y - z), $MachinePrecision] * N[(t - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
code[x_, y_, z_, t_] := Block[{t$95$1 = N[(N[(y - z), $MachinePrecision] * N[(t - z), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, -5e+296], N[(N[(1.0 / N[(z - y), $MachinePrecision]), $MachinePrecision] / N[(N[(z - t), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$1, 5e+154], N[(x / t$95$1), $MachinePrecision], N[(N[(x / N[(z - t), $MachinePrecision]), $MachinePrecision] / N[(z - y), $MachinePrecision]), $MachinePrecision]]]]
\frac{x}{\left(y - z\right) \cdot \left(t - z\right)}
\begin{array}{l}
t_1 := \left(y - z\right) \cdot \left(t - z\right)\\
\mathbf{if}\;t_1 \leq -5 \cdot 10^{+296}:\\
\;\;\;\;\frac{\frac{1}{z - y}}{\frac{z - t}{x}}\\
\mathbf{elif}\;t_1 \leq 5 \cdot 10^{+154}:\\
\;\;\;\;\frac{x}{t_1}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{x}{z - t}}{z - y}\\
\end{array}
Results
| Original | 88.7% |
|---|---|
| Target | 87.4% |
| Herbie | 98.7% |
if (*.f64 (-.f64 y z) (-.f64 t z)) < -5.0000000000000001e296Initial program 70.4%
Simplified99.9%
[Start]70.4 | \[ \frac{x}{\left(y - z\right) \cdot \left(t - z\right)}
\] |
|---|---|
sub-neg [=>]70.4 | \[ \frac{x}{\color{blue}{\left(y + \left(-z\right)\right)} \cdot \left(t - z\right)}
\] |
+-commutative [=>]70.4 | \[ \frac{x}{\color{blue}{\left(\left(-z\right) + y\right)} \cdot \left(t - z\right)}
\] |
neg-sub0 [=>]70.4 | \[ \frac{x}{\left(\color{blue}{\left(0 - z\right)} + y\right) \cdot \left(t - z\right)}
\] |
associate-+l- [=>]70.4 | \[ \frac{x}{\color{blue}{\left(0 - \left(z - y\right)\right)} \cdot \left(t - z\right)}
\] |
sub0-neg [=>]70.4 | \[ \frac{x}{\color{blue}{\left(-\left(z - y\right)\right)} \cdot \left(t - z\right)}
\] |
distribute-lft-neg-out [=>]70.4 | \[ \frac{x}{\color{blue}{-\left(z - y\right) \cdot \left(t - z\right)}}
\] |
distribute-rgt-neg-in [=>]70.4 | \[ \frac{x}{\color{blue}{\left(z - y\right) \cdot \left(-\left(t - z\right)\right)}}
\] |
neg-sub0 [=>]70.4 | \[ \frac{x}{\left(z - y\right) \cdot \color{blue}{\left(0 - \left(t - z\right)\right)}}
\] |
associate-+l- [<=]70.4 | \[ \frac{x}{\left(z - y\right) \cdot \color{blue}{\left(\left(0 - t\right) + z\right)}}
\] |
neg-sub0 [<=]70.4 | \[ \frac{x}{\left(z - y\right) \cdot \left(\color{blue}{\left(-t\right)} + z\right)}
\] |
+-commutative [<=]70.4 | \[ \frac{x}{\left(z - y\right) \cdot \color{blue}{\left(z + \left(-t\right)\right)}}
\] |
sub-neg [<=]70.4 | \[ \frac{x}{\left(z - y\right) \cdot \color{blue}{\left(z - t\right)}}
\] |
associate-/l/ [<=]99.9 | \[ \color{blue}{\frac{\frac{x}{z - t}}{z - y}}
\] |
Applied egg-rr72.6%
[Start]99.9 | \[ \frac{\frac{x}{z - t}}{z - y}
\] |
|---|---|
div-inv [=>]99.8 | \[ \color{blue}{\frac{x}{z - t} \cdot \frac{1}{z - y}}
\] |
clear-num [=>]99.8 | \[ \color{blue}{\frac{1}{\frac{z - t}{x}}} \cdot \frac{1}{z - y}
\] |
associate-*l/ [=>]99.8 | \[ \color{blue}{\frac{1 \cdot \frac{1}{z - y}}{\frac{z - t}{x}}}
\] |
*-un-lft-identity [<=]99.8 | \[ \frac{\color{blue}{\frac{1}{z - y}}}{\frac{z - t}{x}}
\] |
associate-/r/ [=>]72.6 | \[ \color{blue}{\frac{\frac{1}{z - y}}{z - t} \cdot x}
\] |
Applied egg-rr99.8%
[Start]72.6 | \[ \frac{\frac{1}{z - y}}{z - t} \cdot x
\] |
|---|---|
associate-/r/ [<=]99.8 | \[ \color{blue}{\frac{\frac{1}{z - y}}{\frac{z - t}{x}}}
\] |
if -5.0000000000000001e296 < (*.f64 (-.f64 y z) (-.f64 t z)) < 5.00000000000000004e154Initial program 97.7%
if 5.00000000000000004e154 < (*.f64 (-.f64 y z) (-.f64 t z)) Initial program 84.2%
Simplified99.5%
[Start]84.2 | \[ \frac{x}{\left(y - z\right) \cdot \left(t - z\right)}
\] |
|---|---|
sub-neg [=>]84.2 | \[ \frac{x}{\color{blue}{\left(y + \left(-z\right)\right)} \cdot \left(t - z\right)}
\] |
+-commutative [=>]84.2 | \[ \frac{x}{\color{blue}{\left(\left(-z\right) + y\right)} \cdot \left(t - z\right)}
\] |
neg-sub0 [=>]84.2 | \[ \frac{x}{\left(\color{blue}{\left(0 - z\right)} + y\right) \cdot \left(t - z\right)}
\] |
associate-+l- [=>]84.2 | \[ \frac{x}{\color{blue}{\left(0 - \left(z - y\right)\right)} \cdot \left(t - z\right)}
\] |
sub0-neg [=>]84.2 | \[ \frac{x}{\color{blue}{\left(-\left(z - y\right)\right)} \cdot \left(t - z\right)}
\] |
distribute-lft-neg-out [=>]84.2 | \[ \frac{x}{\color{blue}{-\left(z - y\right) \cdot \left(t - z\right)}}
\] |
distribute-rgt-neg-in [=>]84.2 | \[ \frac{x}{\color{blue}{\left(z - y\right) \cdot \left(-\left(t - z\right)\right)}}
\] |
neg-sub0 [=>]84.2 | \[ \frac{x}{\left(z - y\right) \cdot \color{blue}{\left(0 - \left(t - z\right)\right)}}
\] |
associate-+l- [<=]84.2 | \[ \frac{x}{\left(z - y\right) \cdot \color{blue}{\left(\left(0 - t\right) + z\right)}}
\] |
neg-sub0 [<=]84.2 | \[ \frac{x}{\left(z - y\right) \cdot \left(\color{blue}{\left(-t\right)} + z\right)}
\] |
+-commutative [<=]84.2 | \[ \frac{x}{\left(z - y\right) \cdot \color{blue}{\left(z + \left(-t\right)\right)}}
\] |
sub-neg [<=]84.2 | \[ \frac{x}{\left(z - y\right) \cdot \color{blue}{\left(z - t\right)}}
\] |
associate-/l/ [<=]99.5 | \[ \color{blue}{\frac{\frac{x}{z - t}}{z - y}}
\] |
Final simplification98.7%
| Alternative 1 | |
|---|---|
| Accuracy | 98.7% |
| Cost | 1609 |
| Alternative 2 | |
|---|---|
| Accuracy | 76.2% |
| Cost | 1108 |
| Alternative 3 | |
|---|---|
| Accuracy | 81.1% |
| Cost | 1104 |
| Alternative 4 | |
|---|---|
| Accuracy | 93.2% |
| Cost | 1104 |
| Alternative 5 | |
|---|---|
| Accuracy | 80.9% |
| Cost | 1040 |
| Alternative 6 | |
|---|---|
| Accuracy | 65.8% |
| Cost | 976 |
| Alternative 7 | |
|---|---|
| Accuracy | 65.9% |
| Cost | 976 |
| Alternative 8 | |
|---|---|
| Accuracy | 71.9% |
| Cost | 976 |
| Alternative 9 | |
|---|---|
| Accuracy | 78.4% |
| Cost | 976 |
| Alternative 10 | |
|---|---|
| Accuracy | 80.9% |
| Cost | 976 |
| Alternative 11 | |
|---|---|
| Accuracy | 65.9% |
| Cost | 912 |
| Alternative 12 | |
|---|---|
| Accuracy | 78.7% |
| Cost | 712 |
| Alternative 13 | |
|---|---|
| Accuracy | 78.8% |
| Cost | 712 |
| Alternative 14 | |
|---|---|
| Accuracy | 45.6% |
| Cost | 585 |
| Alternative 15 | |
|---|---|
| Accuracy | 62.1% |
| Cost | 585 |
| Alternative 16 | |
|---|---|
| Accuracy | 63.8% |
| Cost | 585 |
| Alternative 17 | |
|---|---|
| Accuracy | 67.5% |
| Cost | 585 |
| Alternative 18 | |
|---|---|
| Accuracy | 38.3% |
| Cost | 320 |
herbie shell --seed 2023140
(FPCore (x y z t)
:name "Data.Random.Distribution.Triangular:triangularCDF from random-fu-0.2.6.2, B"
:precision binary64
:herbie-target
(if (< (/ x (* (- y z) (- t z))) 0.0) (/ (/ x (- y z)) (- t z)) (* x (/ 1.0 (* (- y z) (- t z)))))
(/ x (* (- y z) (- t z))))