
(FPCore (x y z) :precision binary64 (/ (+ x y) (- 1.0 (/ y z))))
double code(double x, double y, double z) {
return (x + y) / (1.0 - (y / z));
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
code = (x + y) / (1.0d0 - (y / z))
end function
public static double code(double x, double y, double z) {
return (x + y) / (1.0 - (y / z));
}
def code(x, y, z): return (x + y) / (1.0 - (y / z))
function code(x, y, z) return Float64(Float64(x + y) / Float64(1.0 - Float64(y / z))) end
function tmp = code(x, y, z) tmp = (x + y) / (1.0 - (y / z)); end
code[x_, y_, z_] := N[(N[(x + y), $MachinePrecision] / N[(1.0 - N[(y / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{x + y}{1 - \frac{y}{z}}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 9 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x y z) :precision binary64 (/ (+ x y) (- 1.0 (/ y z))))
double code(double x, double y, double z) {
return (x + y) / (1.0 - (y / z));
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
code = (x + y) / (1.0d0 - (y / z))
end function
public static double code(double x, double y, double z) {
return (x + y) / (1.0 - (y / z));
}
def code(x, y, z): return (x + y) / (1.0 - (y / z))
function code(x, y, z) return Float64(Float64(x + y) / Float64(1.0 - Float64(y / z))) end
function tmp = code(x, y, z) tmp = (x + y) / (1.0 - (y / z)); end
code[x_, y_, z_] := N[(N[(x + y), $MachinePrecision] / N[(1.0 - N[(y / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{x + y}{1 - \frac{y}{z}}
\end{array}
(FPCore (x y z) :precision binary64 (let* ((t_0 (/ (+ y x) (- 1.0 (/ y z))))) (if (<= t_0 -5e-179) t_0 (if (<= t_0 2e-256) (* (- -1.0 (/ x y)) z) t_0))))
double code(double x, double y, double z) {
double t_0 = (y + x) / (1.0 - (y / z));
double tmp;
if (t_0 <= -5e-179) {
tmp = t_0;
} else if (t_0 <= 2e-256) {
tmp = (-1.0 - (x / y)) * z;
} else {
tmp = t_0;
}
return tmp;
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8) :: t_0
real(8) :: tmp
t_0 = (y + x) / (1.0d0 - (y / z))
if (t_0 <= (-5d-179)) then
tmp = t_0
else if (t_0 <= 2d-256) then
tmp = ((-1.0d0) - (x / y)) * z
else
tmp = t_0
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double t_0 = (y + x) / (1.0 - (y / z));
double tmp;
if (t_0 <= -5e-179) {
tmp = t_0;
} else if (t_0 <= 2e-256) {
tmp = (-1.0 - (x / y)) * z;
} else {
tmp = t_0;
}
return tmp;
}
def code(x, y, z): t_0 = (y + x) / (1.0 - (y / z)) tmp = 0 if t_0 <= -5e-179: tmp = t_0 elif t_0 <= 2e-256: tmp = (-1.0 - (x / y)) * z else: tmp = t_0 return tmp
function code(x, y, z) t_0 = Float64(Float64(y + x) / Float64(1.0 - Float64(y / z))) tmp = 0.0 if (t_0 <= -5e-179) tmp = t_0; elseif (t_0 <= 2e-256) tmp = Float64(Float64(-1.0 - Float64(x / y)) * z); else tmp = t_0; end return tmp end
function tmp_2 = code(x, y, z) t_0 = (y + x) / (1.0 - (y / z)); tmp = 0.0; if (t_0 <= -5e-179) tmp = t_0; elseif (t_0 <= 2e-256) tmp = (-1.0 - (x / y)) * z; else tmp = t_0; end tmp_2 = tmp; end
code[x_, y_, z_] := Block[{t$95$0 = N[(N[(y + x), $MachinePrecision] / N[(1.0 - N[(y / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, -5e-179], t$95$0, If[LessEqual[t$95$0, 2e-256], N[(N[(-1.0 - N[(x / y), $MachinePrecision]), $MachinePrecision] * z), $MachinePrecision], t$95$0]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{y + x}{1 - \frac{y}{z}}\\
\mathbf{if}\;t\_0 \leq -5 \cdot 10^{-179}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;t\_0 \leq 2 \cdot 10^{-256}:\\
\;\;\;\;\left(-1 - \frac{x}{y}\right) \cdot z\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if (/.f64 (+.f64 x y) (-.f64 #s(literal 1 binary64) (/.f64 y z))) < -4.9999999999999998e-179 or 1.99999999999999995e-256 < (/.f64 (+.f64 x y) (-.f64 #s(literal 1 binary64) (/.f64 y z))) Initial program 99.9%
if -4.9999999999999998e-179 < (/.f64 (+.f64 x y) (-.f64 #s(literal 1 binary64) (/.f64 y z))) < 1.99999999999999995e-256Initial program 23.0%
Taylor expanded in z around 0
associate-/l*N/A
*-commutativeN/A
associate-*l*N/A
lower-*.f64N/A
associate-*r/N/A
+-commutativeN/A
distribute-lft-inN/A
mul-1-negN/A
unsub-negN/A
div-subN/A
associate-*l/N/A
metadata-evalN/A
distribute-neg-fracN/A
distribute-lft-neg-outN/A
lft-mult-inverseN/A
metadata-evalN/A
lower--.f64N/A
lower-/.f6499.9
Applied rewrites99.9%
Final simplification99.9%
(FPCore (x y z)
:precision binary64
(let* ((t_0 (* (/ y (- z y)) z)))
(if (<= y -2.9e+97)
t_0
(if (<= y 1.35e-90)
(/ x (- 1.0 (/ y z)))
(if (<= y 4.8e+100) t_0 (* (- -1.0 (/ x y)) z))))))
double code(double x, double y, double z) {
double t_0 = (y / (z - y)) * z;
double tmp;
if (y <= -2.9e+97) {
tmp = t_0;
} else if (y <= 1.35e-90) {
tmp = x / (1.0 - (y / z));
} else if (y <= 4.8e+100) {
tmp = t_0;
} else {
tmp = (-1.0 - (x / y)) * z;
}
return tmp;
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8) :: t_0
real(8) :: tmp
t_0 = (y / (z - y)) * z
if (y <= (-2.9d+97)) then
tmp = t_0
else if (y <= 1.35d-90) then
tmp = x / (1.0d0 - (y / z))
else if (y <= 4.8d+100) then
tmp = t_0
else
tmp = ((-1.0d0) - (x / y)) * z
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double t_0 = (y / (z - y)) * z;
double tmp;
if (y <= -2.9e+97) {
tmp = t_0;
} else if (y <= 1.35e-90) {
tmp = x / (1.0 - (y / z));
} else if (y <= 4.8e+100) {
tmp = t_0;
} else {
tmp = (-1.0 - (x / y)) * z;
}
return tmp;
}
def code(x, y, z): t_0 = (y / (z - y)) * z tmp = 0 if y <= -2.9e+97: tmp = t_0 elif y <= 1.35e-90: tmp = x / (1.0 - (y / z)) elif y <= 4.8e+100: tmp = t_0 else: tmp = (-1.0 - (x / y)) * z return tmp
function code(x, y, z) t_0 = Float64(Float64(y / Float64(z - y)) * z) tmp = 0.0 if (y <= -2.9e+97) tmp = t_0; elseif (y <= 1.35e-90) tmp = Float64(x / Float64(1.0 - Float64(y / z))); elseif (y <= 4.8e+100) tmp = t_0; else tmp = Float64(Float64(-1.0 - Float64(x / y)) * z); end return tmp end
function tmp_2 = code(x, y, z) t_0 = (y / (z - y)) * z; tmp = 0.0; if (y <= -2.9e+97) tmp = t_0; elseif (y <= 1.35e-90) tmp = x / (1.0 - (y / z)); elseif (y <= 4.8e+100) tmp = t_0; else tmp = (-1.0 - (x / y)) * z; end tmp_2 = tmp; end
code[x_, y_, z_] := Block[{t$95$0 = N[(N[(y / N[(z - y), $MachinePrecision]), $MachinePrecision] * z), $MachinePrecision]}, If[LessEqual[y, -2.9e+97], t$95$0, If[LessEqual[y, 1.35e-90], N[(x / N[(1.0 - N[(y / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 4.8e+100], t$95$0, N[(N[(-1.0 - N[(x / y), $MachinePrecision]), $MachinePrecision] * z), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{y}{z - y} \cdot z\\
\mathbf{if}\;y \leq -2.9 \cdot 10^{+97}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;y \leq 1.35 \cdot 10^{-90}:\\
\;\;\;\;\frac{x}{1 - \frac{y}{z}}\\
\mathbf{elif}\;y \leq 4.8 \cdot 10^{+100}:\\
\;\;\;\;t\_0\\
\mathbf{else}:\\
\;\;\;\;\left(-1 - \frac{x}{y}\right) \cdot z\\
\end{array}
\end{array}
if y < -2.89999999999999987e97 or 1.34999999999999998e-90 < y < 4.80000000000000023e100Initial program 81.6%
Taylor expanded in x around 0
*-inversesN/A
div-subN/A
associate-/r/N/A
lower-*.f64N/A
lower-/.f64N/A
lower--.f6467.4
Applied rewrites67.4%
if -2.89999999999999987e97 < y < 1.34999999999999998e-90Initial program 98.7%
Taylor expanded in y around inf
associate--l+N/A
+-commutativeN/A
associate-/l*N/A
associate-*r*N/A
unpow2N/A
associate-/l*N/A
distribute-rgt-out--N/A
lower-fma.f64N/A
lower-/.f64N/A
lower--.f64N/A
mul-1-negN/A
lower-neg.f64N/A
mul-1-negN/A
lower-neg.f6429.2
Applied rewrites29.2%
Taylor expanded in x around inf
lower-/.f64N/A
lower--.f64N/A
lower-/.f6472.0
Applied rewrites72.0%
if 4.80000000000000023e100 < y Initial program 68.1%
Taylor expanded in z around 0
associate-/l*N/A
*-commutativeN/A
associate-*l*N/A
lower-*.f64N/A
associate-*r/N/A
+-commutativeN/A
distribute-lft-inN/A
mul-1-negN/A
unsub-negN/A
div-subN/A
associate-*l/N/A
metadata-evalN/A
distribute-neg-fracN/A
distribute-lft-neg-outN/A
lft-mult-inverseN/A
metadata-evalN/A
lower--.f64N/A
lower-/.f6482.8
Applied rewrites82.8%
(FPCore (x y z)
:precision binary64
(let* ((t_0 (* (/ (+ y x) (- y)) z)))
(if (< y -3.7429310762689856e+171)
t_0
(if (< y 3.5534662456086734e+168) (/ (+ x y) (- 1.0 (/ y z))) t_0))))
double code(double x, double y, double z) {
double t_0 = ((y + x) / -y) * z;
double tmp;
if (y < -3.7429310762689856e+171) {
tmp = t_0;
} else if (y < 3.5534662456086734e+168) {
tmp = (x + y) / (1.0 - (y / z));
} else {
tmp = t_0;
}
return tmp;
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8) :: t_0
real(8) :: tmp
t_0 = ((y + x) / -y) * z
if (y < (-3.7429310762689856d+171)) then
tmp = t_0
else if (y < 3.5534662456086734d+168) then
tmp = (x + y) / (1.0d0 - (y / z))
else
tmp = t_0
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double t_0 = ((y + x) / -y) * z;
double tmp;
if (y < -3.7429310762689856e+171) {
tmp = t_0;
} else if (y < 3.5534662456086734e+168) {
tmp = (x + y) / (1.0 - (y / z));
} else {
tmp = t_0;
}
return tmp;
}
def code(x, y, z): t_0 = ((y + x) / -y) * z tmp = 0 if y < -3.7429310762689856e+171: tmp = t_0 elif y < 3.5534662456086734e+168: tmp = (x + y) / (1.0 - (y / z)) else: tmp = t_0 return tmp
function code(x, y, z) t_0 = Float64(Float64(Float64(y + x) / Float64(-y)) * z) tmp = 0.0 if (y < -3.7429310762689856e+171) tmp = t_0; elseif (y < 3.5534662456086734e+168) tmp = Float64(Float64(x + y) / Float64(1.0 - Float64(y / z))); else tmp = t_0; end return tmp end
function tmp_2 = code(x, y, z) t_0 = ((y + x) / -y) * z; tmp = 0.0; if (y < -3.7429310762689856e+171) tmp = t_0; elseif (y < 3.5534662456086734e+168) tmp = (x + y) / (1.0 - (y / z)); else tmp = t_0; end tmp_2 = tmp; end
code[x_, y_, z_] := Block[{t$95$0 = N[(N[(N[(y + x), $MachinePrecision] / (-y)), $MachinePrecision] * z), $MachinePrecision]}, If[Less[y, -3.7429310762689856e+171], t$95$0, If[Less[y, 3.5534662456086734e+168], N[(N[(x + y), $MachinePrecision] / N[(1.0 - N[(y / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$0]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{y + x}{-y} \cdot z\\
\mathbf{if}\;y < -3.7429310762689856 \cdot 10^{+171}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;y < 3.5534662456086734 \cdot 10^{+168}:\\
\;\;\;\;\frac{x + y}{1 - \frac{y}{z}}\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
herbie shell --seed 2024230
(FPCore (x y z)
:name "Graphics.Rendering.Chart.Backend.Diagrams:calcFontMetrics from Chart-diagrams-1.5.1, A"
:precision binary64
:alt
(! :herbie-platform default (if (< y -3742931076268985600000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000) (* (/ (+ y x) (- y)) z) (if (< y 3553466245608673400000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000) (/ (+ x y) (- 1 (/ y z))) (* (/ (+ y x) (- y)) z))))
(/ (+ x y) (- 1.0 (/ y z))))