
(FPCore (a k m) :precision binary64 (/ (* a (pow k m)) (+ (+ 1.0 (* 10.0 k)) (* k k))))
double code(double a, double k, double m) {
return (a * pow(k, m)) / ((1.0 + (10.0 * k)) + (k * k));
}
real(8) function code(a, k, m)
real(8), intent (in) :: a
real(8), intent (in) :: k
real(8), intent (in) :: m
code = (a * (k ** m)) / ((1.0d0 + (10.0d0 * k)) + (k * k))
end function
public static double code(double a, double k, double m) {
return (a * Math.pow(k, m)) / ((1.0 + (10.0 * k)) + (k * k));
}
def code(a, k, m): return (a * math.pow(k, m)) / ((1.0 + (10.0 * k)) + (k * k))
function code(a, k, m) return Float64(Float64(a * (k ^ m)) / Float64(Float64(1.0 + Float64(10.0 * k)) + Float64(k * k))) end
function tmp = code(a, k, m) tmp = (a * (k ^ m)) / ((1.0 + (10.0 * k)) + (k * k)); end
code[a_, k_, m_] := N[(N[(a * N[Power[k, m], $MachinePrecision]), $MachinePrecision] / N[(N[(1.0 + N[(10.0 * k), $MachinePrecision]), $MachinePrecision] + N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{a \cdot {k}^{m}}{\left(1 + 10 \cdot k\right) + k \cdot k}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 14 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (a k m) :precision binary64 (/ (* a (pow k m)) (+ (+ 1.0 (* 10.0 k)) (* k k))))
double code(double a, double k, double m) {
return (a * pow(k, m)) / ((1.0 + (10.0 * k)) + (k * k));
}
real(8) function code(a, k, m)
real(8), intent (in) :: a
real(8), intent (in) :: k
real(8), intent (in) :: m
code = (a * (k ** m)) / ((1.0d0 + (10.0d0 * k)) + (k * k))
end function
public static double code(double a, double k, double m) {
return (a * Math.pow(k, m)) / ((1.0 + (10.0 * k)) + (k * k));
}
def code(a, k, m): return (a * math.pow(k, m)) / ((1.0 + (10.0 * k)) + (k * k))
function code(a, k, m) return Float64(Float64(a * (k ^ m)) / Float64(Float64(1.0 + Float64(10.0 * k)) + Float64(k * k))) end
function tmp = code(a, k, m) tmp = (a * (k ^ m)) / ((1.0 + (10.0 * k)) + (k * k)); end
code[a_, k_, m_] := N[(N[(a * N[Power[k, m], $MachinePrecision]), $MachinePrecision] / N[(N[(1.0 + N[(10.0 * k), $MachinePrecision]), $MachinePrecision] + N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{a \cdot {k}^{m}}{\left(1 + 10 \cdot k\right) + k \cdot k}
\end{array}
(FPCore (a k m)
:precision binary64
(let* ((t_0 (* a (pow k m))))
(if (<= (/ t_0 (+ (+ 1.0 (* k 10.0)) (* k k))) INFINITY)
(/ t_0 (fma (+ k 10.0) k 1.0))
(* (* k k) (* k (* a (fma k 1e-6 -1e-5)))))))
double code(double a, double k, double m) {
double t_0 = a * pow(k, m);
double tmp;
if ((t_0 / ((1.0 + (k * 10.0)) + (k * k))) <= ((double) INFINITY)) {
tmp = t_0 / fma((k + 10.0), k, 1.0);
} else {
tmp = (k * k) * (k * (a * fma(k, 1e-6, -1e-5)));
}
return tmp;
}
function code(a, k, m) t_0 = Float64(a * (k ^ m)) tmp = 0.0 if (Float64(t_0 / Float64(Float64(1.0 + Float64(k * 10.0)) + Float64(k * k))) <= Inf) tmp = Float64(t_0 / fma(Float64(k + 10.0), k, 1.0)); else tmp = Float64(Float64(k * k) * Float64(k * Float64(a * fma(k, 1e-6, -1e-5)))); end return tmp end
code[a_, k_, m_] := Block[{t$95$0 = N[(a * N[Power[k, m], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[N[(t$95$0 / N[(N[(1.0 + N[(k * 10.0), $MachinePrecision]), $MachinePrecision] + N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], Infinity], N[(t$95$0 / N[(N[(k + 10.0), $MachinePrecision] * k + 1.0), $MachinePrecision]), $MachinePrecision], N[(N[(k * k), $MachinePrecision] * N[(k * N[(a * N[(k * 1e-6 + -1e-5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := a \cdot {k}^{m}\\
\mathbf{if}\;\frac{t\_0}{\left(1 + k \cdot 10\right) + k \cdot k} \leq \infty:\\
\;\;\;\;\frac{t\_0}{\mathsf{fma}\left(k + 10, k, 1\right)}\\
\mathbf{else}:\\
\;\;\;\;\left(k \cdot k\right) \cdot \left(k \cdot \left(a \cdot \mathsf{fma}\left(k, 10^{-6}, -1 \cdot 10^{-5}\right)\right)\right)\\
\end{array}
\end{array}
if (/.f64 (*.f64 a (pow.f64 k m)) (+.f64 (+.f64 #s(literal 1 binary64) (*.f64 #s(literal 10 binary64) k)) (*.f64 k k))) < +inf.0Initial program 98.0%
lift-*.f64N/A
lift-*.f64N/A
associate-+l+N/A
+-commutativeN/A
lift-*.f64N/A
lift-*.f64N/A
distribute-rgt-outN/A
*-commutativeN/A
lower-fma.f64N/A
lower-+.f6498.0
Applied rewrites98.0%
if +inf.0 < (/.f64 (*.f64 a (pow.f64 k m)) (+.f64 (+.f64 #s(literal 1 binary64) (*.f64 #s(literal 10 binary64) k)) (*.f64 k k))) Initial program 0.0%
Taylor expanded in m around 0
lower-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
+-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
lower-+.f641.6
Applied rewrites1.6%
lift-+.f64N/A
*-commutativeN/A
lift-+.f64N/A
flip-+N/A
lift--.f64N/A
associate-*l/N/A
lift--.f64N/A
flip--N/A
flip-+N/A
lift--.f64N/A
associate-/r/N/A
times-fracN/A
lower-fma.f64N/A
Applied rewrites1.6%
Taylor expanded in k around -inf
sub-negN/A
metadata-evalN/A
+-commutativeN/A
lower-+.f64N/A
associate-*r/N/A
lower-/.f64N/A
Applied rewrites1.6%
Taylor expanded in k around 0
unpow3N/A
unpow2N/A
associate-*l*N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
distribute-lft-outN/A
lower-*.f64N/A
lower-fma.f64100.0
Applied rewrites100.0%
Final simplification98.2%
(FPCore (a k m) :precision binary64 (if (<= (/ (* a (pow k m)) (+ (+ 1.0 (* k 10.0)) (* k k))) INFINITY) (* a (/ (pow k m) (fma k (+ k 10.0) 1.0))) (* (* k k) (* k (* a (fma k 1e-6 -1e-5))))))
double code(double a, double k, double m) {
double tmp;
if (((a * pow(k, m)) / ((1.0 + (k * 10.0)) + (k * k))) <= ((double) INFINITY)) {
tmp = a * (pow(k, m) / fma(k, (k + 10.0), 1.0));
} else {
tmp = (k * k) * (k * (a * fma(k, 1e-6, -1e-5)));
}
return tmp;
}
function code(a, k, m) tmp = 0.0 if (Float64(Float64(a * (k ^ m)) / Float64(Float64(1.0 + Float64(k * 10.0)) + Float64(k * k))) <= Inf) tmp = Float64(a * Float64((k ^ m) / fma(k, Float64(k + 10.0), 1.0))); else tmp = Float64(Float64(k * k) * Float64(k * Float64(a * fma(k, 1e-6, -1e-5)))); end return tmp end
code[a_, k_, m_] := If[LessEqual[N[(N[(a * N[Power[k, m], $MachinePrecision]), $MachinePrecision] / N[(N[(1.0 + N[(k * 10.0), $MachinePrecision]), $MachinePrecision] + N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], Infinity], N[(a * N[(N[Power[k, m], $MachinePrecision] / N[(k * N[(k + 10.0), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(k * k), $MachinePrecision] * N[(k * N[(a * N[(k * 1e-6 + -1e-5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\frac{a \cdot {k}^{m}}{\left(1 + k \cdot 10\right) + k \cdot k} \leq \infty:\\
\;\;\;\;a \cdot \frac{{k}^{m}}{\mathsf{fma}\left(k, k + 10, 1\right)}\\
\mathbf{else}:\\
\;\;\;\;\left(k \cdot k\right) \cdot \left(k \cdot \left(a \cdot \mathsf{fma}\left(k, 10^{-6}, -1 \cdot 10^{-5}\right)\right)\right)\\
\end{array}
\end{array}
if (/.f64 (*.f64 a (pow.f64 k m)) (+.f64 (+.f64 #s(literal 1 binary64) (*.f64 #s(literal 10 binary64) k)) (*.f64 k k))) < +inf.0Initial program 98.0%
lift-pow.f64N/A
*-commutativeN/A
lift-*.f64N/A
lift-+.f64N/A
lift-*.f64N/A
lift-+.f64N/A
associate-*l/N/A
lower-*.f64N/A
lower-/.f6498.0
lift-+.f64N/A
lift-+.f64N/A
associate-+l+N/A
+-commutativeN/A
lift-*.f64N/A
lift-*.f64N/A
distribute-rgt-outN/A
lower-fma.f64N/A
lower-+.f6498.0
Applied rewrites98.0%
if +inf.0 < (/.f64 (*.f64 a (pow.f64 k m)) (+.f64 (+.f64 #s(literal 1 binary64) (*.f64 #s(literal 10 binary64) k)) (*.f64 k k))) Initial program 0.0%
Taylor expanded in m around 0
lower-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
+-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
lower-+.f641.6
Applied rewrites1.6%
lift-+.f64N/A
*-commutativeN/A
lift-+.f64N/A
flip-+N/A
lift--.f64N/A
associate-*l/N/A
lift--.f64N/A
flip--N/A
flip-+N/A
lift--.f64N/A
associate-/r/N/A
times-fracN/A
lower-fma.f64N/A
Applied rewrites1.6%
Taylor expanded in k around -inf
sub-negN/A
metadata-evalN/A
+-commutativeN/A
lower-+.f64N/A
associate-*r/N/A
lower-/.f64N/A
Applied rewrites1.6%
Taylor expanded in k around 0
unpow3N/A
unpow2N/A
associate-*l*N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
distribute-lft-outN/A
lower-*.f64N/A
lower-fma.f64100.0
Applied rewrites100.0%
Final simplification98.2%
(FPCore (a k m)
:precision binary64
(let* ((t_0 (* a (pow k m))))
(if (<= m -4.6e-5)
t_0
(if (<= m 1.02e-44) (/ a (fma k (+ k 10.0) 1.0)) t_0))))
double code(double a, double k, double m) {
double t_0 = a * pow(k, m);
double tmp;
if (m <= -4.6e-5) {
tmp = t_0;
} else if (m <= 1.02e-44) {
tmp = a / fma(k, (k + 10.0), 1.0);
} else {
tmp = t_0;
}
return tmp;
}
function code(a, k, m) t_0 = Float64(a * (k ^ m)) tmp = 0.0 if (m <= -4.6e-5) tmp = t_0; elseif (m <= 1.02e-44) tmp = Float64(a / fma(k, Float64(k + 10.0), 1.0)); else tmp = t_0; end return tmp end
code[a_, k_, m_] := Block[{t$95$0 = N[(a * N[Power[k, m], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[m, -4.6e-5], t$95$0, If[LessEqual[m, 1.02e-44], N[(a / N[(k * N[(k + 10.0), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision], t$95$0]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := a \cdot {k}^{m}\\
\mathbf{if}\;m \leq -4.6 \cdot 10^{-5}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;m \leq 1.02 \cdot 10^{-44}:\\
\;\;\;\;\frac{a}{\mathsf{fma}\left(k, k + 10, 1\right)}\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if m < -4.6e-5 or 1.0199999999999999e-44 < m Initial program 89.8%
Taylor expanded in k around 0
lower-*.f64N/A
lower-pow.f64100.0
Applied rewrites100.0%
if -4.6e-5 < m < 1.0199999999999999e-44Initial program 94.1%
Taylor expanded in m around 0
lower-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
+-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
lower-+.f6494.0
Applied rewrites94.0%
Final simplification98.1%
(FPCore (a k m)
:precision binary64
(if (<= m -6.8e+24)
(/ (fma (/ a k) (+ (/ 100.0 k) -10.0) a) (* k k))
(if (<= m 1.65)
(/ a (fma k (+ k 10.0) 1.0))
(*
k
(*
(* k k)
(fma
k
(- (* k (* a (fma -9.9e-9 k 1e-7))))
(* a (fma k 1e-6 -1e-5))))))))
double code(double a, double k, double m) {
double tmp;
if (m <= -6.8e+24) {
tmp = fma((a / k), ((100.0 / k) + -10.0), a) / (k * k);
} else if (m <= 1.65) {
tmp = a / fma(k, (k + 10.0), 1.0);
} else {
tmp = k * ((k * k) * fma(k, -(k * (a * fma(-9.9e-9, k, 1e-7))), (a * fma(k, 1e-6, -1e-5))));
}
return tmp;
}
function code(a, k, m) tmp = 0.0 if (m <= -6.8e+24) tmp = Float64(fma(Float64(a / k), Float64(Float64(100.0 / k) + -10.0), a) / Float64(k * k)); elseif (m <= 1.65) tmp = Float64(a / fma(k, Float64(k + 10.0), 1.0)); else tmp = Float64(k * Float64(Float64(k * k) * fma(k, Float64(-Float64(k * Float64(a * fma(-9.9e-9, k, 1e-7)))), Float64(a * fma(k, 1e-6, -1e-5))))); end return tmp end
code[a_, k_, m_] := If[LessEqual[m, -6.8e+24], N[(N[(N[(a / k), $MachinePrecision] * N[(N[(100.0 / k), $MachinePrecision] + -10.0), $MachinePrecision] + a), $MachinePrecision] / N[(k * k), $MachinePrecision]), $MachinePrecision], If[LessEqual[m, 1.65], N[(a / N[(k * N[(k + 10.0), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision], N[(k * N[(N[(k * k), $MachinePrecision] * N[(k * (-N[(k * N[(a * N[(-9.9e-9 * k + 1e-7), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]) + N[(a * N[(k * 1e-6 + -1e-5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq -6.8 \cdot 10^{+24}:\\
\;\;\;\;\frac{\mathsf{fma}\left(\frac{a}{k}, \frac{100}{k} + -10, a\right)}{k \cdot k}\\
\mathbf{elif}\;m \leq 1.65:\\
\;\;\;\;\frac{a}{\mathsf{fma}\left(k, k + 10, 1\right)}\\
\mathbf{else}:\\
\;\;\;\;k \cdot \left(\left(k \cdot k\right) \cdot \mathsf{fma}\left(k, -k \cdot \left(a \cdot \mathsf{fma}\left(-9.9 \cdot 10^{-9}, k, 10^{-7}\right)\right), a \cdot \mathsf{fma}\left(k, 10^{-6}, -1 \cdot 10^{-5}\right)\right)\right)\\
\end{array}
\end{array}
if m < -6.8000000000000001e24Initial program 100.0%
Taylor expanded in m around 0
lower-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
+-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
lower-+.f6433.1
Applied rewrites33.1%
Taylor expanded in k around inf
unpow2N/A
associate-*l*N/A
+-commutativeN/A
distribute-rgt-inN/A
*-lft-identityN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
lower-*.f64N/A
+-commutativeN/A
lower-+.f6443.2
Applied rewrites43.2%
Taylor expanded in k around inf
metadata-evalN/A
cancel-sign-sub-invN/A
associate--r+N/A
lower-/.f64N/A
Applied rewrites65.5%
if -6.8000000000000001e24 < m < 1.6499999999999999Initial program 94.4%
Taylor expanded in m around 0
lower-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
+-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
lower-+.f6492.2
Applied rewrites92.2%
if 1.6499999999999999 < m Initial program 80.4%
Taylor expanded in m around 0
lower-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
+-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
lower-+.f642.9
Applied rewrites2.9%
lift-+.f64N/A
*-commutativeN/A
lift-+.f64N/A
flip-+N/A
lift--.f64N/A
associate-*l/N/A
lift--.f64N/A
flip--N/A
flip-+N/A
lift--.f64N/A
associate-/r/N/A
times-fracN/A
lower-fma.f64N/A
Applied rewrites2.9%
Taylor expanded in k around -inf
sub-negN/A
metadata-evalN/A
+-commutativeN/A
lower-+.f64N/A
associate-*r/N/A
lower-/.f64N/A
Applied rewrites46.0%
Taylor expanded in k around 0
Applied rewrites76.0%
Final simplification78.1%
(FPCore (a k m)
:precision binary64
(if (<= m -6.8e+24)
(/ (fma (/ a k) (+ (/ 100.0 k) -10.0) a) (* k k))
(if (<= m 1.65)
(/ a (fma k (+ k 10.0) 1.0))
(* (* k k) (* k (* a (fma k 1e-6 -1e-5)))))))
double code(double a, double k, double m) {
double tmp;
if (m <= -6.8e+24) {
tmp = fma((a / k), ((100.0 / k) + -10.0), a) / (k * k);
} else if (m <= 1.65) {
tmp = a / fma(k, (k + 10.0), 1.0);
} else {
tmp = (k * k) * (k * (a * fma(k, 1e-6, -1e-5)));
}
return tmp;
}
function code(a, k, m) tmp = 0.0 if (m <= -6.8e+24) tmp = Float64(fma(Float64(a / k), Float64(Float64(100.0 / k) + -10.0), a) / Float64(k * k)); elseif (m <= 1.65) tmp = Float64(a / fma(k, Float64(k + 10.0), 1.0)); else tmp = Float64(Float64(k * k) * Float64(k * Float64(a * fma(k, 1e-6, -1e-5)))); end return tmp end
code[a_, k_, m_] := If[LessEqual[m, -6.8e+24], N[(N[(N[(a / k), $MachinePrecision] * N[(N[(100.0 / k), $MachinePrecision] + -10.0), $MachinePrecision] + a), $MachinePrecision] / N[(k * k), $MachinePrecision]), $MachinePrecision], If[LessEqual[m, 1.65], N[(a / N[(k * N[(k + 10.0), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision], N[(N[(k * k), $MachinePrecision] * N[(k * N[(a * N[(k * 1e-6 + -1e-5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq -6.8 \cdot 10^{+24}:\\
\;\;\;\;\frac{\mathsf{fma}\left(\frac{a}{k}, \frac{100}{k} + -10, a\right)}{k \cdot k}\\
\mathbf{elif}\;m \leq 1.65:\\
\;\;\;\;\frac{a}{\mathsf{fma}\left(k, k + 10, 1\right)}\\
\mathbf{else}:\\
\;\;\;\;\left(k \cdot k\right) \cdot \left(k \cdot \left(a \cdot \mathsf{fma}\left(k, 10^{-6}, -1 \cdot 10^{-5}\right)\right)\right)\\
\end{array}
\end{array}
if m < -6.8000000000000001e24Initial program 100.0%
Taylor expanded in m around 0
lower-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
+-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
lower-+.f6433.1
Applied rewrites33.1%
Taylor expanded in k around inf
unpow2N/A
associate-*l*N/A
+-commutativeN/A
distribute-rgt-inN/A
*-lft-identityN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
lower-*.f64N/A
+-commutativeN/A
lower-+.f6443.2
Applied rewrites43.2%
Taylor expanded in k around inf
metadata-evalN/A
cancel-sign-sub-invN/A
associate--r+N/A
lower-/.f64N/A
Applied rewrites65.5%
if -6.8000000000000001e24 < m < 1.6499999999999999Initial program 94.4%
Taylor expanded in m around 0
lower-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
+-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
lower-+.f6492.2
Applied rewrites92.2%
if 1.6499999999999999 < m Initial program 80.4%
Taylor expanded in m around 0
lower-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
+-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
lower-+.f642.9
Applied rewrites2.9%
lift-+.f64N/A
*-commutativeN/A
lift-+.f64N/A
flip-+N/A
lift--.f64N/A
associate-*l/N/A
lift--.f64N/A
flip--N/A
flip-+N/A
lift--.f64N/A
associate-/r/N/A
times-fracN/A
lower-fma.f64N/A
Applied rewrites2.9%
Taylor expanded in k around -inf
sub-negN/A
metadata-evalN/A
+-commutativeN/A
lower-+.f64N/A
associate-*r/N/A
lower-/.f64N/A
Applied rewrites46.0%
Taylor expanded in k around 0
unpow3N/A
unpow2N/A
associate-*l*N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
distribute-lft-outN/A
lower-*.f64N/A
lower-fma.f6474.1
Applied rewrites74.1%
Final simplification77.3%
(FPCore (a k m)
:precision binary64
(if (<= m -6.8e+24)
(/ a (* k k))
(if (<= m 1.65)
(/ a (fma k (+ k 10.0) 1.0))
(* (* k k) (* k (* a (fma k 1e-6 -1e-5)))))))
double code(double a, double k, double m) {
double tmp;
if (m <= -6.8e+24) {
tmp = a / (k * k);
} else if (m <= 1.65) {
tmp = a / fma(k, (k + 10.0), 1.0);
} else {
tmp = (k * k) * (k * (a * fma(k, 1e-6, -1e-5)));
}
return tmp;
}
function code(a, k, m) tmp = 0.0 if (m <= -6.8e+24) tmp = Float64(a / Float64(k * k)); elseif (m <= 1.65) tmp = Float64(a / fma(k, Float64(k + 10.0), 1.0)); else tmp = Float64(Float64(k * k) * Float64(k * Float64(a * fma(k, 1e-6, -1e-5)))); end return tmp end
code[a_, k_, m_] := If[LessEqual[m, -6.8e+24], N[(a / N[(k * k), $MachinePrecision]), $MachinePrecision], If[LessEqual[m, 1.65], N[(a / N[(k * N[(k + 10.0), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision], N[(N[(k * k), $MachinePrecision] * N[(k * N[(a * N[(k * 1e-6 + -1e-5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq -6.8 \cdot 10^{+24}:\\
\;\;\;\;\frac{a}{k \cdot k}\\
\mathbf{elif}\;m \leq 1.65:\\
\;\;\;\;\frac{a}{\mathsf{fma}\left(k, k + 10, 1\right)}\\
\mathbf{else}:\\
\;\;\;\;\left(k \cdot k\right) \cdot \left(k \cdot \left(a \cdot \mathsf{fma}\left(k, 10^{-6}, -1 \cdot 10^{-5}\right)\right)\right)\\
\end{array}
\end{array}
if m < -6.8000000000000001e24Initial program 100.0%
Taylor expanded in m around 0
lower-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
+-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
lower-+.f6433.1
Applied rewrites33.1%
Taylor expanded in k around inf
lower-/.f64N/A
unpow2N/A
lower-*.f6460.7
Applied rewrites60.7%
if -6.8000000000000001e24 < m < 1.6499999999999999Initial program 94.4%
Taylor expanded in m around 0
lower-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
+-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
lower-+.f6492.2
Applied rewrites92.2%
if 1.6499999999999999 < m Initial program 80.4%
Taylor expanded in m around 0
lower-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
+-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
lower-+.f642.9
Applied rewrites2.9%
lift-+.f64N/A
*-commutativeN/A
lift-+.f64N/A
flip-+N/A
lift--.f64N/A
associate-*l/N/A
lift--.f64N/A
flip--N/A
flip-+N/A
lift--.f64N/A
associate-/r/N/A
times-fracN/A
lower-fma.f64N/A
Applied rewrites2.9%
Taylor expanded in k around -inf
sub-negN/A
metadata-evalN/A
+-commutativeN/A
lower-+.f64N/A
associate-*r/N/A
lower-/.f64N/A
Applied rewrites46.0%
Taylor expanded in k around 0
unpow3N/A
unpow2N/A
associate-*l*N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
distribute-lft-outN/A
lower-*.f64N/A
lower-fma.f6474.1
Applied rewrites74.1%
Final simplification75.8%
(FPCore (a k m)
:precision binary64
(if (<= m -6.8e+24)
(/ a (* k k))
(if (<= m 1.65)
(/ a (fma k (+ k 10.0) 1.0))
(* -1e-5 (* a (* k (* k k)))))))
double code(double a, double k, double m) {
double tmp;
if (m <= -6.8e+24) {
tmp = a / (k * k);
} else if (m <= 1.65) {
tmp = a / fma(k, (k + 10.0), 1.0);
} else {
tmp = -1e-5 * (a * (k * (k * k)));
}
return tmp;
}
function code(a, k, m) tmp = 0.0 if (m <= -6.8e+24) tmp = Float64(a / Float64(k * k)); elseif (m <= 1.65) tmp = Float64(a / fma(k, Float64(k + 10.0), 1.0)); else tmp = Float64(-1e-5 * Float64(a * Float64(k * Float64(k * k)))); end return tmp end
code[a_, k_, m_] := If[LessEqual[m, -6.8e+24], N[(a / N[(k * k), $MachinePrecision]), $MachinePrecision], If[LessEqual[m, 1.65], N[(a / N[(k * N[(k + 10.0), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision], N[(-1e-5 * N[(a * N[(k * N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq -6.8 \cdot 10^{+24}:\\
\;\;\;\;\frac{a}{k \cdot k}\\
\mathbf{elif}\;m \leq 1.65:\\
\;\;\;\;\frac{a}{\mathsf{fma}\left(k, k + 10, 1\right)}\\
\mathbf{else}:\\
\;\;\;\;-1 \cdot 10^{-5} \cdot \left(a \cdot \left(k \cdot \left(k \cdot k\right)\right)\right)\\
\end{array}
\end{array}
if m < -6.8000000000000001e24Initial program 100.0%
Taylor expanded in m around 0
lower-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
+-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
lower-+.f6433.1
Applied rewrites33.1%
Taylor expanded in k around inf
lower-/.f64N/A
unpow2N/A
lower-*.f6460.7
Applied rewrites60.7%
if -6.8000000000000001e24 < m < 1.6499999999999999Initial program 94.4%
Taylor expanded in m around 0
lower-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
+-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
lower-+.f6492.2
Applied rewrites92.2%
if 1.6499999999999999 < m Initial program 80.4%
Taylor expanded in m around 0
lower-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
+-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
lower-+.f642.9
Applied rewrites2.9%
lift-+.f64N/A
*-commutativeN/A
lift-+.f64N/A
flip-+N/A
lift--.f64N/A
associate-*l/N/A
lift--.f64N/A
flip--N/A
flip-+N/A
lift--.f64N/A
associate-/r/N/A
times-fracN/A
lower-fma.f64N/A
Applied rewrites2.9%
Taylor expanded in k around -inf
sub-negN/A
metadata-evalN/A
+-commutativeN/A
lower-+.f64N/A
associate-*r/N/A
lower-/.f64N/A
Applied rewrites46.0%
Taylor expanded in k around 0
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
cube-multN/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6466.0
Applied rewrites66.0%
Final simplification72.9%
(FPCore (a k m) :precision binary64 (if (<= m -2.65e-30) (/ a (* k k)) (if (<= m 1.65) (/ a (fma k 10.0 1.0)) (* -1e-5 (* a (* k (* k k)))))))
double code(double a, double k, double m) {
double tmp;
if (m <= -2.65e-30) {
tmp = a / (k * k);
} else if (m <= 1.65) {
tmp = a / fma(k, 10.0, 1.0);
} else {
tmp = -1e-5 * (a * (k * (k * k)));
}
return tmp;
}
function code(a, k, m) tmp = 0.0 if (m <= -2.65e-30) tmp = Float64(a / Float64(k * k)); elseif (m <= 1.65) tmp = Float64(a / fma(k, 10.0, 1.0)); else tmp = Float64(-1e-5 * Float64(a * Float64(k * Float64(k * k)))); end return tmp end
code[a_, k_, m_] := If[LessEqual[m, -2.65e-30], N[(a / N[(k * k), $MachinePrecision]), $MachinePrecision], If[LessEqual[m, 1.65], N[(a / N[(k * 10.0 + 1.0), $MachinePrecision]), $MachinePrecision], N[(-1e-5 * N[(a * N[(k * N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq -2.65 \cdot 10^{-30}:\\
\;\;\;\;\frac{a}{k \cdot k}\\
\mathbf{elif}\;m \leq 1.65:\\
\;\;\;\;\frac{a}{\mathsf{fma}\left(k, 10, 1\right)}\\
\mathbf{else}:\\
\;\;\;\;-1 \cdot 10^{-5} \cdot \left(a \cdot \left(k \cdot \left(k \cdot k\right)\right)\right)\\
\end{array}
\end{array}
if m < -2.64999999999999987e-30Initial program 98.9%
Taylor expanded in m around 0
lower-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
+-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
lower-+.f6435.9
Applied rewrites35.9%
Taylor expanded in k around inf
lower-/.f64N/A
unpow2N/A
lower-*.f6459.1
Applied rewrites59.1%
if -2.64999999999999987e-30 < m < 1.6499999999999999Initial program 95.1%
Taylor expanded in m around 0
lower-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
+-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
lower-+.f6495.1
Applied rewrites95.1%
Taylor expanded in k around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f6460.4
Applied rewrites60.4%
if 1.6499999999999999 < m Initial program 80.4%
Taylor expanded in m around 0
lower-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
+-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
lower-+.f642.9
Applied rewrites2.9%
lift-+.f64N/A
*-commutativeN/A
lift-+.f64N/A
flip-+N/A
lift--.f64N/A
associate-*l/N/A
lift--.f64N/A
flip--N/A
flip-+N/A
lift--.f64N/A
associate-/r/N/A
times-fracN/A
lower-fma.f64N/A
Applied rewrites2.9%
Taylor expanded in k around -inf
sub-negN/A
metadata-evalN/A
+-commutativeN/A
lower-+.f64N/A
associate-*r/N/A
lower-/.f64N/A
Applied rewrites46.0%
Taylor expanded in k around 0
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
cube-multN/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6466.0
Applied rewrites66.0%
Final simplification62.0%
(FPCore (a k m) :precision binary64 (if (<= m -2.65e-30) (/ a (* k k)) (if (<= m 1.9) (/ a (fma k 10.0 1.0)) (* k (* a (fma 99.0 k -10.0))))))
double code(double a, double k, double m) {
double tmp;
if (m <= -2.65e-30) {
tmp = a / (k * k);
} else if (m <= 1.9) {
tmp = a / fma(k, 10.0, 1.0);
} else {
tmp = k * (a * fma(99.0, k, -10.0));
}
return tmp;
}
function code(a, k, m) tmp = 0.0 if (m <= -2.65e-30) tmp = Float64(a / Float64(k * k)); elseif (m <= 1.9) tmp = Float64(a / fma(k, 10.0, 1.0)); else tmp = Float64(k * Float64(a * fma(99.0, k, -10.0))); end return tmp end
code[a_, k_, m_] := If[LessEqual[m, -2.65e-30], N[(a / N[(k * k), $MachinePrecision]), $MachinePrecision], If[LessEqual[m, 1.9], N[(a / N[(k * 10.0 + 1.0), $MachinePrecision]), $MachinePrecision], N[(k * N[(a * N[(99.0 * k + -10.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq -2.65 \cdot 10^{-30}:\\
\;\;\;\;\frac{a}{k \cdot k}\\
\mathbf{elif}\;m \leq 1.9:\\
\;\;\;\;\frac{a}{\mathsf{fma}\left(k, 10, 1\right)}\\
\mathbf{else}:\\
\;\;\;\;k \cdot \left(a \cdot \mathsf{fma}\left(99, k, -10\right)\right)\\
\end{array}
\end{array}
if m < -2.64999999999999987e-30Initial program 98.9%
Taylor expanded in m around 0
lower-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
+-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
lower-+.f6435.9
Applied rewrites35.9%
Taylor expanded in k around inf
lower-/.f64N/A
unpow2N/A
lower-*.f6459.1
Applied rewrites59.1%
if -2.64999999999999987e-30 < m < 1.8999999999999999Initial program 95.1%
Taylor expanded in m around 0
lower-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
+-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
lower-+.f6495.1
Applied rewrites95.1%
Taylor expanded in k around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f6460.4
Applied rewrites60.4%
if 1.8999999999999999 < m Initial program 80.4%
Taylor expanded in m around 0
lower-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
+-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
lower-+.f642.9
Applied rewrites2.9%
Taylor expanded in k around 0
+-commutativeN/A
lower-fma.f64N/A
Applied rewrites23.9%
Taylor expanded in k around -inf
Applied rewrites36.4%
Taylor expanded in k around 0
lower-*.f64N/A
cancel-sign-sub-invN/A
metadata-evalN/A
*-commutativeN/A
associate-*r*N/A
distribute-rgt-outN/A
lower-*.f64N/A
lower-fma.f6436.3
Applied rewrites36.3%
(FPCore (a k m) :precision binary64 (if (<= m -0.0023) (/ a (* k 10.0)) (if (<= m 0.49) a (* k (* a (fma 99.0 k -10.0))))))
double code(double a, double k, double m) {
double tmp;
if (m <= -0.0023) {
tmp = a / (k * 10.0);
} else if (m <= 0.49) {
tmp = a;
} else {
tmp = k * (a * fma(99.0, k, -10.0));
}
return tmp;
}
function code(a, k, m) tmp = 0.0 if (m <= -0.0023) tmp = Float64(a / Float64(k * 10.0)); elseif (m <= 0.49) tmp = a; else tmp = Float64(k * Float64(a * fma(99.0, k, -10.0))); end return tmp end
code[a_, k_, m_] := If[LessEqual[m, -0.0023], N[(a / N[(k * 10.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[m, 0.49], a, N[(k * N[(a * N[(99.0 * k + -10.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq -0.0023:\\
\;\;\;\;\frac{a}{k \cdot 10}\\
\mathbf{elif}\;m \leq 0.49:\\
\;\;\;\;a\\
\mathbf{else}:\\
\;\;\;\;k \cdot \left(a \cdot \mathsf{fma}\left(99, k, -10\right)\right)\\
\end{array}
\end{array}
if m < -0.0023Initial program 100.0%
Taylor expanded in m around 0
lower-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
+-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
lower-+.f6434.3
Applied rewrites34.3%
Taylor expanded in k around inf
unpow2N/A
associate-*l*N/A
+-commutativeN/A
distribute-rgt-inN/A
*-lft-identityN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
lower-*.f64N/A
+-commutativeN/A
lower-+.f6444.1
Applied rewrites44.1%
Taylor expanded in k around 0
*-commutativeN/A
lower-*.f6426.0
Applied rewrites26.0%
if -0.0023 < m < 0.48999999999999999Initial program 94.2%
Taylor expanded in m around 0
lower-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
+-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
lower-+.f6493.1
Applied rewrites93.1%
Taylor expanded in k around 0
Applied rewrites45.6%
/-rgt-identity45.6
Applied rewrites45.6%
if 0.48999999999999999 < m Initial program 80.4%
Taylor expanded in m around 0
lower-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
+-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
lower-+.f642.9
Applied rewrites2.9%
Taylor expanded in k around 0
+-commutativeN/A
lower-fma.f64N/A
Applied rewrites23.9%
Taylor expanded in k around -inf
Applied rewrites36.4%
Taylor expanded in k around 0
lower-*.f64N/A
cancel-sign-sub-invN/A
metadata-evalN/A
*-commutativeN/A
associate-*r*N/A
distribute-rgt-outN/A
lower-*.f64N/A
lower-fma.f6436.3
Applied rewrites36.3%
(FPCore (a k m) :precision binary64 (if (<= m 1.9) (/ a (* k k)) (* k (* a (fma 99.0 k -10.0)))))
double code(double a, double k, double m) {
double tmp;
if (m <= 1.9) {
tmp = a / (k * k);
} else {
tmp = k * (a * fma(99.0, k, -10.0));
}
return tmp;
}
function code(a, k, m) tmp = 0.0 if (m <= 1.9) tmp = Float64(a / Float64(k * k)); else tmp = Float64(k * Float64(a * fma(99.0, k, -10.0))); end return tmp end
code[a_, k_, m_] := If[LessEqual[m, 1.9], N[(a / N[(k * k), $MachinePrecision]), $MachinePrecision], N[(k * N[(a * N[(99.0 * k + -10.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq 1.9:\\
\;\;\;\;\frac{a}{k \cdot k}\\
\mathbf{else}:\\
\;\;\;\;k \cdot \left(a \cdot \mathsf{fma}\left(99, k, -10\right)\right)\\
\end{array}
\end{array}
if m < 1.8999999999999999Initial program 97.1%
Taylor expanded in m around 0
lower-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
+-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
lower-+.f6463.4
Applied rewrites63.4%
Taylor expanded in k around inf
lower-/.f64N/A
unpow2N/A
lower-*.f6454.9
Applied rewrites54.9%
if 1.8999999999999999 < m Initial program 80.4%
Taylor expanded in m around 0
lower-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
+-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
lower-+.f642.9
Applied rewrites2.9%
Taylor expanded in k around 0
+-commutativeN/A
lower-fma.f64N/A
Applied rewrites23.9%
Taylor expanded in k around -inf
Applied rewrites36.4%
Taylor expanded in k around 0
lower-*.f64N/A
cancel-sign-sub-invN/A
metadata-evalN/A
*-commutativeN/A
associate-*r*N/A
distribute-rgt-outN/A
lower-*.f64N/A
lower-fma.f6436.3
Applied rewrites36.3%
(FPCore (a k m) :precision binary64 (if (<= m 0.49) a (* k (* a (fma 99.0 k -10.0)))))
double code(double a, double k, double m) {
double tmp;
if (m <= 0.49) {
tmp = a;
} else {
tmp = k * (a * fma(99.0, k, -10.0));
}
return tmp;
}
function code(a, k, m) tmp = 0.0 if (m <= 0.49) tmp = a; else tmp = Float64(k * Float64(a * fma(99.0, k, -10.0))); end return tmp end
code[a_, k_, m_] := If[LessEqual[m, 0.49], a, N[(k * N[(a * N[(99.0 * k + -10.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq 0.49:\\
\;\;\;\;a\\
\mathbf{else}:\\
\;\;\;\;k \cdot \left(a \cdot \mathsf{fma}\left(99, k, -10\right)\right)\\
\end{array}
\end{array}
if m < 0.48999999999999999Initial program 97.1%
Taylor expanded in m around 0
lower-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
+-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
lower-+.f6463.4
Applied rewrites63.4%
Taylor expanded in k around 0
Applied rewrites24.6%
/-rgt-identity24.6
Applied rewrites24.6%
if 0.48999999999999999 < m Initial program 80.4%
Taylor expanded in m around 0
lower-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
+-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
lower-+.f642.9
Applied rewrites2.9%
Taylor expanded in k around 0
+-commutativeN/A
lower-fma.f64N/A
Applied rewrites23.9%
Taylor expanded in k around -inf
Applied rewrites36.4%
Taylor expanded in k around 0
lower-*.f64N/A
cancel-sign-sub-invN/A
metadata-evalN/A
*-commutativeN/A
associate-*r*N/A
distribute-rgt-outN/A
lower-*.f64N/A
lower-fma.f6436.3
Applied rewrites36.3%
(FPCore (a k m) :precision binary64 (if (<= m 0.49) a (* k (* a -10.0))))
double code(double a, double k, double m) {
double tmp;
if (m <= 0.49) {
tmp = a;
} else {
tmp = k * (a * -10.0);
}
return tmp;
}
real(8) function code(a, k, m)
real(8), intent (in) :: a
real(8), intent (in) :: k
real(8), intent (in) :: m
real(8) :: tmp
if (m <= 0.49d0) then
tmp = a
else
tmp = k * (a * (-10.0d0))
end if
code = tmp
end function
public static double code(double a, double k, double m) {
double tmp;
if (m <= 0.49) {
tmp = a;
} else {
tmp = k * (a * -10.0);
}
return tmp;
}
def code(a, k, m): tmp = 0 if m <= 0.49: tmp = a else: tmp = k * (a * -10.0) return tmp
function code(a, k, m) tmp = 0.0 if (m <= 0.49) tmp = a; else tmp = Float64(k * Float64(a * -10.0)); end return tmp end
function tmp_2 = code(a, k, m) tmp = 0.0; if (m <= 0.49) tmp = a; else tmp = k * (a * -10.0); end tmp_2 = tmp; end
code[a_, k_, m_] := If[LessEqual[m, 0.49], a, N[(k * N[(a * -10.0), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq 0.49:\\
\;\;\;\;a\\
\mathbf{else}:\\
\;\;\;\;k \cdot \left(a \cdot -10\right)\\
\end{array}
\end{array}
if m < 0.48999999999999999Initial program 97.1%
Taylor expanded in m around 0
lower-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
+-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
lower-+.f6463.4
Applied rewrites63.4%
Taylor expanded in k around 0
Applied rewrites24.6%
/-rgt-identity24.6
Applied rewrites24.6%
if 0.48999999999999999 < m Initial program 80.4%
Taylor expanded in m around 0
lower-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
+-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
lower-+.f642.9
Applied rewrites2.9%
Taylor expanded in k around 0
+-commutativeN/A
lower-fma.f64N/A
Applied rewrites23.9%
Taylor expanded in k around -inf
Applied rewrites36.4%
Taylor expanded in k around 0
associate-*r*N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6421.9
Applied rewrites21.9%
(FPCore (a k m) :precision binary64 a)
double code(double a, double k, double m) {
return a;
}
real(8) function code(a, k, m)
real(8), intent (in) :: a
real(8), intent (in) :: k
real(8), intent (in) :: m
code = a
end function
public static double code(double a, double k, double m) {
return a;
}
def code(a, k, m): return a
function code(a, k, m) return a end
function tmp = code(a, k, m) tmp = a; end
code[a_, k_, m_] := a
\begin{array}{l}
\\
a
\end{array}
Initial program 91.1%
Taylor expanded in m around 0
lower-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
+-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
lower-+.f6441.6
Applied rewrites41.6%
Taylor expanded in k around 0
Applied rewrites17.1%
/-rgt-identity17.1
Applied rewrites17.1%
herbie shell --seed 2024219
(FPCore (a k m)
:name "Falkner and Boettcher, Appendix A"
:precision binary64
(/ (* a (pow k m)) (+ (+ 1.0 (* 10.0 k)) (* k k))))