
(FPCore (x y z) :precision binary64 (* (/ 1.0 2.0) (+ x (* y (sqrt z)))))
double code(double x, double y, double z) {
return (1.0 / 2.0) * (x + (y * sqrt(z)));
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
code = (1.0d0 / 2.0d0) * (x + (y * sqrt(z)))
end function
public static double code(double x, double y, double z) {
return (1.0 / 2.0) * (x + (y * Math.sqrt(z)));
}
def code(x, y, z): return (1.0 / 2.0) * (x + (y * math.sqrt(z)))
function code(x, y, z) return Float64(Float64(1.0 / 2.0) * Float64(x + Float64(y * sqrt(z)))) end
function tmp = code(x, y, z) tmp = (1.0 / 2.0) * (x + (y * sqrt(z))); end
code[x_, y_, z_] := N[(N[(1.0 / 2.0), $MachinePrecision] * N[(x + N[(y * N[Sqrt[z], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{1}{2} \cdot \left(x + y \cdot \sqrt{z}\right)
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 4 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x y z) :precision binary64 (* (/ 1.0 2.0) (+ x (* y (sqrt z)))))
double code(double x, double y, double z) {
return (1.0 / 2.0) * (x + (y * sqrt(z)));
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
code = (1.0d0 / 2.0d0) * (x + (y * sqrt(z)))
end function
public static double code(double x, double y, double z) {
return (1.0 / 2.0) * (x + (y * Math.sqrt(z)));
}
def code(x, y, z): return (1.0 / 2.0) * (x + (y * math.sqrt(z)))
function code(x, y, z) return Float64(Float64(1.0 / 2.0) * Float64(x + Float64(y * sqrt(z)))) end
function tmp = code(x, y, z) tmp = (1.0 / 2.0) * (x + (y * sqrt(z))); end
code[x_, y_, z_] := N[(N[(1.0 / 2.0), $MachinePrecision] * N[(x + N[(y * N[Sqrt[z], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{1}{2} \cdot \left(x + y \cdot \sqrt{z}\right)
\end{array}
(FPCore (x y z) :precision binary64 (* 0.5 (+ x (/ y (pow z -0.5)))))
double code(double x, double y, double z) {
return 0.5 * (x + (y / pow(z, -0.5)));
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
code = 0.5d0 * (x + (y / (z ** (-0.5d0))))
end function
public static double code(double x, double y, double z) {
return 0.5 * (x + (y / Math.pow(z, -0.5)));
}
def code(x, y, z): return 0.5 * (x + (y / math.pow(z, -0.5)))
function code(x, y, z) return Float64(0.5 * Float64(x + Float64(y / (z ^ -0.5)))) end
function tmp = code(x, y, z) tmp = 0.5 * (x + (y / (z ^ -0.5))); end
code[x_, y_, z_] := N[(0.5 * N[(x + N[(y / N[Power[z, -0.5], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
0.5 \cdot \left(x + \frac{y}{{z}^{-0.5}}\right)
\end{array}
Initial program 99.8%
metadata-eval99.8%
Simplified99.8%
add099.8%
flip-+70.1%
metadata-eval70.1%
fma-neg70.1%
metadata-eval70.1%
fma-define70.1%
add070.1%
*-commutative70.1%
*-commutative70.1%
swap-sqr61.2%
add-sqr-sqrt61.2%
pow261.2%
fma-neg61.2%
metadata-eval61.2%
fma-define61.2%
add061.2%
Applied egg-rr61.2%
associate-/l*68.4%
unpow268.4%
times-frac92.2%
*-inverses92.2%
associate-*r/92.2%
*-lft-identity92.2%
Simplified92.2%
clear-num92.2%
associate-/r/91.5%
clear-num91.5%
Applied egg-rr91.5%
associate-*l/92.7%
associate-/l*99.7%
Applied egg-rr99.7%
*-un-lft-identity99.7%
add-sqr-sqrt99.4%
times-frac99.5%
pow1/299.5%
pow199.5%
pow-div99.5%
metadata-eval99.5%
pow1/299.5%
pow199.5%
pow-div99.5%
metadata-eval99.5%
Applied egg-rr99.5%
associate-*l/99.5%
*-lft-identity99.5%
associate-/r*99.5%
rem-square-sqrt99.8%
Simplified99.8%
Final simplification99.8%
(FPCore (x y z) :precision binary64 (if (or (<= x -6.8e+16) (not (<= x 9.5e-21))) (* 0.5 x) (* 0.5 (* y (sqrt z)))))
double code(double x, double y, double z) {
double tmp;
if ((x <= -6.8e+16) || !(x <= 9.5e-21)) {
tmp = 0.5 * x;
} else {
tmp = 0.5 * (y * sqrt(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) :: tmp
if ((x <= (-6.8d+16)) .or. (.not. (x <= 9.5d-21))) then
tmp = 0.5d0 * x
else
tmp = 0.5d0 * (y * sqrt(z))
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double tmp;
if ((x <= -6.8e+16) || !(x <= 9.5e-21)) {
tmp = 0.5 * x;
} else {
tmp = 0.5 * (y * Math.sqrt(z));
}
return tmp;
}
def code(x, y, z): tmp = 0 if (x <= -6.8e+16) or not (x <= 9.5e-21): tmp = 0.5 * x else: tmp = 0.5 * (y * math.sqrt(z)) return tmp
function code(x, y, z) tmp = 0.0 if ((x <= -6.8e+16) || !(x <= 9.5e-21)) tmp = Float64(0.5 * x); else tmp = Float64(0.5 * Float64(y * sqrt(z))); end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if ((x <= -6.8e+16) || ~((x <= 9.5e-21))) tmp = 0.5 * x; else tmp = 0.5 * (y * sqrt(z)); end tmp_2 = tmp; end
code[x_, y_, z_] := If[Or[LessEqual[x, -6.8e+16], N[Not[LessEqual[x, 9.5e-21]], $MachinePrecision]], N[(0.5 * x), $MachinePrecision], N[(0.5 * N[(y * N[Sqrt[z], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -6.8 \cdot 10^{+16} \lor \neg \left(x \leq 9.5 \cdot 10^{-21}\right):\\
\;\;\;\;0.5 \cdot x\\
\mathbf{else}:\\
\;\;\;\;0.5 \cdot \left(y \cdot \sqrt{z}\right)\\
\end{array}
\end{array}
if x < -6.8e16 or 9.4999999999999994e-21 < x Initial program 99.9%
metadata-eval99.9%
Simplified99.9%
Taylor expanded in x around inf 78.3%
if -6.8e16 < x < 9.4999999999999994e-21Initial program 99.7%
metadata-eval99.7%
Simplified99.7%
Taylor expanded in x around 0 77.0%
Final simplification77.6%
(FPCore (x y z) :precision binary64 (* 0.5 (+ x (* y (sqrt z)))))
double code(double x, double y, double z) {
return 0.5 * (x + (y * sqrt(z)));
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
code = 0.5d0 * (x + (y * sqrt(z)))
end function
public static double code(double x, double y, double z) {
return 0.5 * (x + (y * Math.sqrt(z)));
}
def code(x, y, z): return 0.5 * (x + (y * math.sqrt(z)))
function code(x, y, z) return Float64(0.5 * Float64(x + Float64(y * sqrt(z)))) end
function tmp = code(x, y, z) tmp = 0.5 * (x + (y * sqrt(z))); end
code[x_, y_, z_] := N[(0.5 * N[(x + N[(y * N[Sqrt[z], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
0.5 \cdot \left(x + y \cdot \sqrt{z}\right)
\end{array}
Initial program 99.8%
metadata-eval99.8%
Simplified99.8%
Final simplification99.8%
(FPCore (x y z) :precision binary64 (* 0.5 x))
double code(double x, double y, double z) {
return 0.5 * x;
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
code = 0.5d0 * x
end function
public static double code(double x, double y, double z) {
return 0.5 * x;
}
def code(x, y, z): return 0.5 * x
function code(x, y, z) return Float64(0.5 * x) end
function tmp = code(x, y, z) tmp = 0.5 * x; end
code[x_, y_, z_] := N[(0.5 * x), $MachinePrecision]
\begin{array}{l}
\\
0.5 \cdot x
\end{array}
Initial program 99.8%
metadata-eval99.8%
Simplified99.8%
Taylor expanded in x around inf 51.3%
Final simplification51.3%
herbie shell --seed 2024034
(FPCore (x y z)
:name "Diagrams.Solve.Polynomial:quadForm from diagrams-solve-0.1, B"
:precision binary64
(* (/ 1.0 2.0) (+ x (* y (sqrt z)))))