
(FPCore (x y) :precision binary64 (exp (* (* x y) y)))
double code(double x, double y) {
return exp(((x * y) * y));
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = exp(((x * y) * y))
end function
public static double code(double x, double y) {
return Math.exp(((x * y) * y));
}
def code(x, y): return math.exp(((x * y) * y))
function code(x, y) return exp(Float64(Float64(x * y) * y)) end
function tmp = code(x, y) tmp = exp(((x * y) * y)); end
code[x_, y_] := N[Exp[N[(N[(x * y), $MachinePrecision] * y), $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
\\
e^{\left(x \cdot y\right) \cdot y}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 5 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x y) :precision binary64 (exp (* (* x y) y)))
double code(double x, double y) {
return exp(((x * y) * y));
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = exp(((x * y) * y))
end function
public static double code(double x, double y) {
return Math.exp(((x * y) * y));
}
def code(x, y): return math.exp(((x * y) * y))
function code(x, y) return exp(Float64(Float64(x * y) * y)) end
function tmp = code(x, y) tmp = exp(((x * y) * y)); end
code[x_, y_] := N[Exp[N[(N[(x * y), $MachinePrecision] * y), $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
\\
e^{\left(x \cdot y\right) \cdot y}
\end{array}
(FPCore (x y) :precision binary64 (* (exp (* x (* (* y y) 0.6666666666666666))) (pow E (* x (* (* y y) 0.3333333333333333)))))
double code(double x, double y) {
return exp((x * ((y * y) * 0.6666666666666666))) * pow(((double) M_E), (x * ((y * y) * 0.3333333333333333)));
}
public static double code(double x, double y) {
return Math.exp((x * ((y * y) * 0.6666666666666666))) * Math.pow(Math.E, (x * ((y * y) * 0.3333333333333333)));
}
def code(x, y): return math.exp((x * ((y * y) * 0.6666666666666666))) * math.pow(math.e, (x * ((y * y) * 0.3333333333333333)))
function code(x, y) return Float64(exp(Float64(x * Float64(Float64(y * y) * 0.6666666666666666))) * (exp(1) ^ Float64(x * Float64(Float64(y * y) * 0.3333333333333333)))) end
function tmp = code(x, y) tmp = exp((x * ((y * y) * 0.6666666666666666))) * (2.71828182845904523536 ^ (x * ((y * y) * 0.3333333333333333))); end
code[x_, y_] := N[(N[Exp[N[(x * N[(N[(y * y), $MachinePrecision] * 0.6666666666666666), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * N[Power[E, N[(x * N[(N[(y * y), $MachinePrecision] * 0.3333333333333333), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
e^{x \cdot \left(\left(y \cdot y\right) \cdot 0.6666666666666666\right)} \cdot {e}^{\left(x \cdot \left(\left(y \cdot y\right) \cdot 0.3333333333333333\right)\right)}
\end{array}
Initial program 100.0%
associate-*l*100.0%
Simplified100.0%
*-un-lft-identity100.0%
exp-prod100.0%
add-log-exp100.0%
add-cube-cbrt99.9%
log-prod99.9%
unpow-prod-up99.9%
Applied egg-rr99.9%
associate-*r*99.9%
metadata-eval99.9%
*-commutative99.9%
associate-*r*99.9%
*-commutative99.9%
*-commutative99.9%
*-commutative99.9%
associate-*l*99.9%
Simplified99.9%
e-exp-199.9%
pow-exp100.0%
*-un-lft-identity100.0%
Applied egg-rr100.0%
unpow2100.0%
Applied egg-rr100.0%
unpow2100.0%
Applied egg-rr100.0%
(FPCore (x y) :precision binary64 (pow (sqrt E) (* x (* (* y y) 2.0))))
double code(double x, double y) {
return pow(sqrt(((double) M_E)), (x * ((y * y) * 2.0)));
}
public static double code(double x, double y) {
return Math.pow(Math.sqrt(Math.E), (x * ((y * y) * 2.0)));
}
def code(x, y): return math.pow(math.sqrt(math.e), (x * ((y * y) * 2.0)))
function code(x, y) return sqrt(exp(1)) ^ Float64(x * Float64(Float64(y * y) * 2.0)) end
function tmp = code(x, y) tmp = sqrt(2.71828182845904523536) ^ (x * ((y * y) * 2.0)); end
code[x_, y_] := N[Power[N[Sqrt[E], $MachinePrecision], N[(x * N[(N[(y * y), $MachinePrecision] * 2.0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
\\
{\left(\sqrt{e}\right)}^{\left(x \cdot \left(\left(y \cdot y\right) \cdot 2\right)\right)}
\end{array}
Initial program 100.0%
associate-*l*100.0%
Simplified100.0%
*-un-lft-identity100.0%
exp-prod100.0%
add-log-exp100.0%
add-cube-cbrt99.9%
log-prod99.9%
unpow-prod-up99.9%
Applied egg-rr99.9%
associate-*r*99.9%
metadata-eval99.9%
*-commutative99.9%
associate-*r*99.9%
*-commutative99.9%
*-commutative99.9%
*-commutative99.9%
associate-*l*99.9%
Simplified99.9%
Applied egg-rr99.9%
fma-undefine99.9%
+-rgt-identity99.9%
pow-sqr100.0%
count-2100.0%
distribute-lft-out100.0%
count-2100.0%
Simplified100.0%
unpow2100.0%
Applied egg-rr100.0%
Final simplification100.0%
(FPCore (x y) :precision binary64 (+ 1.0 (expm1 (* x (* y y)))))
double code(double x, double y) {
return 1.0 + expm1((x * (y * y)));
}
public static double code(double x, double y) {
return 1.0 + Math.expm1((x * (y * y)));
}
def code(x, y): return 1.0 + math.expm1((x * (y * y)))
function code(x, y) return Float64(1.0 + expm1(Float64(x * Float64(y * y)))) end
code[x_, y_] := N[(1.0 + N[(Exp[N[(x * N[(y * y), $MachinePrecision]), $MachinePrecision]] - 1), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
1 + \mathsf{expm1}\left(x \cdot \left(y \cdot y\right)\right)
\end{array}
Initial program 100.0%
associate-*l*100.0%
Simplified100.0%
log1p-expm1-u100.0%
log1p-undefine100.0%
add-exp-log100.0%
pow2100.0%
Applied egg-rr100.0%
unpow2100.0%
Applied egg-rr100.0%
(FPCore (x y) :precision binary64 (exp (* x (* y y))))
double code(double x, double y) {
return exp((x * (y * y)));
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = exp((x * (y * y)))
end function
public static double code(double x, double y) {
return Math.exp((x * (y * y)));
}
def code(x, y): return math.exp((x * (y * y)))
function code(x, y) return exp(Float64(x * Float64(y * y))) end
function tmp = code(x, y) tmp = exp((x * (y * y))); end
code[x_, y_] := N[Exp[N[(x * N[(y * y), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
\\
e^{x \cdot \left(y \cdot y\right)}
\end{array}
Initial program 100.0%
associate-*l*100.0%
Simplified100.0%
(FPCore (x y) :precision binary64 1.0)
double code(double x, double y) {
return 1.0;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = 1.0d0
end function
public static double code(double x, double y) {
return 1.0;
}
def code(x, y): return 1.0
function code(x, y) return 1.0 end
function tmp = code(x, y) tmp = 1.0; end
code[x_, y_] := 1.0
\begin{array}{l}
\\
1
\end{array}
Initial program 100.0%
associate-*l*100.0%
Simplified100.0%
Taylor expanded in x around 0 55.3%
herbie shell --seed 2024170
(FPCore (x y)
:name "Data.Random.Distribution.Normal:normalF from random-fu-0.2.6.2"
:precision binary64
(exp (* (* x y) y)))