
(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 (if (<= (/ (* (pow k m) a) (+ (* k k) (+ (* 10.0 k) 1.0))) INFINITY) (* (/ (pow k m) (fma (+ 10.0 k) k 1.0)) a) (fma (* (fma 99.0 k -10.0) k) a a)))
double code(double a, double k, double m) {
double tmp;
if (((pow(k, m) * a) / ((k * k) + ((10.0 * k) + 1.0))) <= ((double) INFINITY)) {
tmp = (pow(k, m) / fma((10.0 + k), k, 1.0)) * a;
} else {
tmp = fma((fma(99.0, k, -10.0) * k), a, a);
}
return tmp;
}
function code(a, k, m) tmp = 0.0 if (Float64(Float64((k ^ m) * a) / Float64(Float64(k * k) + Float64(Float64(10.0 * k) + 1.0))) <= Inf) tmp = Float64(Float64((k ^ m) / fma(Float64(10.0 + k), k, 1.0)) * a); else tmp = fma(Float64(fma(99.0, k, -10.0) * k), a, a); end return tmp end
code[a_, k_, m_] := If[LessEqual[N[(N[(N[Power[k, m], $MachinePrecision] * a), $MachinePrecision] / N[(N[(k * k), $MachinePrecision] + N[(N[(10.0 * k), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], Infinity], N[(N[(N[Power[k, m], $MachinePrecision] / N[(N[(10.0 + k), $MachinePrecision] * k + 1.0), $MachinePrecision]), $MachinePrecision] * a), $MachinePrecision], N[(N[(N[(99.0 * k + -10.0), $MachinePrecision] * k), $MachinePrecision] * a + a), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\frac{{k}^{m} \cdot a}{k \cdot k + \left(10 \cdot k + 1\right)} \leq \infty:\\
\;\;\;\;\frac{{k}^{m}}{\mathsf{fma}\left(10 + k, k, 1\right)} \cdot a\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(99, k, -10\right) \cdot k, a, a\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
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
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/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
associate-*l*N/A
unpow2N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
lower-fma.f64N/A
Applied rewrites1.6%
Taylor expanded in k around 0
Applied rewrites75.8%
Applied rewrites100.0%
Final simplification98.2%
(FPCore (a k m) :precision binary64 (if (<= (/ (* (pow k m) a) (+ (* k k) (+ (* 10.0 k) 1.0))) 0.0) (* (* -10.0 a) k) (* (fma -10.0 k 1.0) a)))
double code(double a, double k, double m) {
double tmp;
if (((pow(k, m) * a) / ((k * k) + ((10.0 * k) + 1.0))) <= 0.0) {
tmp = (-10.0 * a) * k;
} else {
tmp = fma(-10.0, k, 1.0) * a;
}
return tmp;
}
function code(a, k, m) tmp = 0.0 if (Float64(Float64((k ^ m) * a) / Float64(Float64(k * k) + Float64(Float64(10.0 * k) + 1.0))) <= 0.0) tmp = Float64(Float64(-10.0 * a) * k); else tmp = Float64(fma(-10.0, k, 1.0) * a); end return tmp end
code[a_, k_, m_] := If[LessEqual[N[(N[(N[Power[k, m], $MachinePrecision] * a), $MachinePrecision] / N[(N[(k * k), $MachinePrecision] + N[(N[(10.0 * k), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 0.0], N[(N[(-10.0 * a), $MachinePrecision] * k), $MachinePrecision], N[(N[(-10.0 * k + 1.0), $MachinePrecision] * a), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\frac{{k}^{m} \cdot a}{k \cdot k + \left(10 \cdot k + 1\right)} \leq 0:\\
\;\;\;\;\left(-10 \cdot a\right) \cdot k\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(-10, k, 1\right) \cdot a\\
\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))) < 0.0Initial program 97.5%
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
associate-*l*N/A
unpow2N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
lower-fma.f64N/A
Applied rewrites46.6%
Taylor expanded in k around 0
Applied rewrites14.0%
Taylor expanded in k around inf
Applied rewrites9.1%
if 0.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 73.5%
lift-/.f64N/A
lift-*.f64N/A
associate-/l*N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f6473.5
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
+-commutativeN/A
lower-+.f6473.5
Applied rewrites73.5%
Taylor expanded in m around 0
lower-/.f64N/A
+-commutativeN/A
*-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
distribute-lft1-inN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
distribute-rgt-inN/A
*-lft-identityN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
lower-+.f6444.4
Applied rewrites44.4%
Taylor expanded in k around 0
Applied rewrites35.7%
Final simplification16.6%
(FPCore (a k m) :precision binary64 (if (<= k 0.00165) (* (pow k m) a) (* (pow k (+ -1.0 m)) (/ a k))))
double code(double a, double k, double m) {
double tmp;
if (k <= 0.00165) {
tmp = pow(k, m) * a;
} else {
tmp = pow(k, (-1.0 + m)) * (a / k);
}
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 (k <= 0.00165d0) then
tmp = (k ** m) * a
else
tmp = (k ** ((-1.0d0) + m)) * (a / k)
end if
code = tmp
end function
public static double code(double a, double k, double m) {
double tmp;
if (k <= 0.00165) {
tmp = Math.pow(k, m) * a;
} else {
tmp = Math.pow(k, (-1.0 + m)) * (a / k);
}
return tmp;
}
def code(a, k, m): tmp = 0 if k <= 0.00165: tmp = math.pow(k, m) * a else: tmp = math.pow(k, (-1.0 + m)) * (a / k) return tmp
function code(a, k, m) tmp = 0.0 if (k <= 0.00165) tmp = Float64((k ^ m) * a); else tmp = Float64((k ^ Float64(-1.0 + m)) * Float64(a / k)); end return tmp end
function tmp_2 = code(a, k, m) tmp = 0.0; if (k <= 0.00165) tmp = (k ^ m) * a; else tmp = (k ^ (-1.0 + m)) * (a / k); end tmp_2 = tmp; end
code[a_, k_, m_] := If[LessEqual[k, 0.00165], N[(N[Power[k, m], $MachinePrecision] * a), $MachinePrecision], N[(N[Power[k, N[(-1.0 + m), $MachinePrecision]], $MachinePrecision] * N[(a / k), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;k \leq 0.00165:\\
\;\;\;\;{k}^{m} \cdot a\\
\mathbf{else}:\\
\;\;\;\;{k}^{\left(-1 + m\right)} \cdot \frac{a}{k}\\
\end{array}
\end{array}
if k < 0.00165Initial program 95.0%
lift-/.f64N/A
lift-*.f64N/A
associate-/l*N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f6495.1
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
+-commutativeN/A
lower-+.f6495.1
Applied rewrites95.1%
Taylor expanded in k around 0
lower-pow.f6498.9
Applied rewrites98.9%
if 0.00165 < k Initial program 83.2%
Taylor expanded in k around inf
*-commutativeN/A
associate-/l*N/A
*-commutativeN/A
*-commutativeN/A
associate-*r*N/A
exp-prodN/A
neg-mul-1N/A
log-recN/A
remove-double-negN/A
rem-exp-logN/A
lower-*.f64N/A
unpow2N/A
associate-/r*N/A
lower-/.f64N/A
lower-/.f64N/A
lower-pow.f6489.2
Applied rewrites89.2%
Applied rewrites96.5%
Final simplification98.0%
(FPCore (a k m)
:precision binary64
(let* ((t_0 (* (pow k m) a)))
(if (<= m -0.014)
t_0
(if (<= m 8.5e-35) (* (/ 1.0 (fma (+ 10.0 k) k 1.0)) a) t_0))))
double code(double a, double k, double m) {
double t_0 = pow(k, m) * a;
double tmp;
if (m <= -0.014) {
tmp = t_0;
} else if (m <= 8.5e-35) {
tmp = (1.0 / fma((10.0 + k), k, 1.0)) * a;
} else {
tmp = t_0;
}
return tmp;
}
function code(a, k, m) t_0 = Float64((k ^ m) * a) tmp = 0.0 if (m <= -0.014) tmp = t_0; elseif (m <= 8.5e-35) tmp = Float64(Float64(1.0 / fma(Float64(10.0 + k), k, 1.0)) * a); else tmp = t_0; end return tmp end
code[a_, k_, m_] := Block[{t$95$0 = N[(N[Power[k, m], $MachinePrecision] * a), $MachinePrecision]}, If[LessEqual[m, -0.014], t$95$0, If[LessEqual[m, 8.5e-35], N[(N[(1.0 / N[(N[(10.0 + k), $MachinePrecision] * k + 1.0), $MachinePrecision]), $MachinePrecision] * a), $MachinePrecision], t$95$0]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := {k}^{m} \cdot a\\
\mathbf{if}\;m \leq -0.014:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;m \leq 8.5 \cdot 10^{-35}:\\
\;\;\;\;\frac{1}{\mathsf{fma}\left(10 + k, k, 1\right)} \cdot a\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if m < -0.0140000000000000003 or 8.5000000000000001e-35 < m Initial program 88.5%
lift-/.f64N/A
lift-*.f64N/A
associate-/l*N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f6488.5
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
+-commutativeN/A
lower-+.f6488.5
Applied rewrites88.5%
Taylor expanded in k around 0
lower-pow.f6499.4
Applied rewrites99.4%
if -0.0140000000000000003 < m < 8.5000000000000001e-35Initial program 95.3%
lift-/.f64N/A
lift-*.f64N/A
associate-/l*N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f6495.4
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
+-commutativeN/A
lower-+.f6495.4
Applied rewrites95.4%
Taylor expanded in m around 0
lower-/.f64N/A
+-commutativeN/A
*-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
distribute-lft1-inN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
distribute-rgt-inN/A
*-lft-identityN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
lower-+.f6495.4
Applied rewrites95.4%
(FPCore (a k m)
:precision binary64
(if (<= m -1.0)
(/ (* (/ (/ a k) k) 99.0) (* k k))
(if (<= m 0.65)
(* (/ 1.0 (fma (+ 10.0 k) k 1.0)) a)
(* (* (* k a) 99.0) k))))
double code(double a, double k, double m) {
double tmp;
if (m <= -1.0) {
tmp = (((a / k) / k) * 99.0) / (k * k);
} else if (m <= 0.65) {
tmp = (1.0 / fma((10.0 + k), k, 1.0)) * a;
} else {
tmp = ((k * a) * 99.0) * k;
}
return tmp;
}
function code(a, k, m) tmp = 0.0 if (m <= -1.0) tmp = Float64(Float64(Float64(Float64(a / k) / k) * 99.0) / Float64(k * k)); elseif (m <= 0.65) tmp = Float64(Float64(1.0 / fma(Float64(10.0 + k), k, 1.0)) * a); else tmp = Float64(Float64(Float64(k * a) * 99.0) * k); end return tmp end
code[a_, k_, m_] := If[LessEqual[m, -1.0], N[(N[(N[(N[(a / k), $MachinePrecision] / k), $MachinePrecision] * 99.0), $MachinePrecision] / N[(k * k), $MachinePrecision]), $MachinePrecision], If[LessEqual[m, 0.65], N[(N[(1.0 / N[(N[(10.0 + k), $MachinePrecision] * k + 1.0), $MachinePrecision]), $MachinePrecision] * a), $MachinePrecision], N[(N[(N[(k * a), $MachinePrecision] * 99.0), $MachinePrecision] * k), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq -1:\\
\;\;\;\;\frac{\frac{\frac{a}{k}}{k} \cdot 99}{k \cdot k}\\
\mathbf{elif}\;m \leq 0.65:\\
\;\;\;\;\frac{1}{\mathsf{fma}\left(10 + k, k, 1\right)} \cdot a\\
\mathbf{else}:\\
\;\;\;\;\left(\left(k \cdot a\right) \cdot 99\right) \cdot k\\
\end{array}
\end{array}
if m < -1Initial 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
associate-*l*N/A
unpow2N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
lower-fma.f64N/A
Applied rewrites37.7%
Taylor expanded in k around inf
Applied rewrites68.7%
Taylor expanded in k around 0
Applied rewrites73.5%
if -1 < m < 0.650000000000000022Initial program 94.5%
lift-/.f64N/A
lift-*.f64N/A
associate-/l*N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f6494.5
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
+-commutativeN/A
lower-+.f6494.6
Applied rewrites94.6%
Taylor expanded in m around 0
lower-/.f64N/A
+-commutativeN/A
*-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
distribute-lft1-inN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
distribute-rgt-inN/A
*-lft-identityN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
lower-+.f6493.5
Applied rewrites93.5%
if 0.650000000000000022 < m Initial program 75.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
associate-*l*N/A
unpow2N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
lower-fma.f64N/A
Applied rewrites3.1%
Taylor expanded in k around 0
Applied rewrites24.9%
Taylor expanded in k around inf
Applied rewrites52.3%
(FPCore (a k m)
:precision binary64
(if (<= m -0.98)
(/ 1.0 (/ (* k k) a))
(if (<= m 0.65)
(* (/ 1.0 (fma (+ 10.0 k) k 1.0)) a)
(* (* (* k a) 99.0) k))))
double code(double a, double k, double m) {
double tmp;
if (m <= -0.98) {
tmp = 1.0 / ((k * k) / a);
} else if (m <= 0.65) {
tmp = (1.0 / fma((10.0 + k), k, 1.0)) * a;
} else {
tmp = ((k * a) * 99.0) * k;
}
return tmp;
}
function code(a, k, m) tmp = 0.0 if (m <= -0.98) tmp = Float64(1.0 / Float64(Float64(k * k) / a)); elseif (m <= 0.65) tmp = Float64(Float64(1.0 / fma(Float64(10.0 + k), k, 1.0)) * a); else tmp = Float64(Float64(Float64(k * a) * 99.0) * k); end return tmp end
code[a_, k_, m_] := If[LessEqual[m, -0.98], N[(1.0 / N[(N[(k * k), $MachinePrecision] / a), $MachinePrecision]), $MachinePrecision], If[LessEqual[m, 0.65], N[(N[(1.0 / N[(N[(10.0 + k), $MachinePrecision] * k + 1.0), $MachinePrecision]), $MachinePrecision] * a), $MachinePrecision], N[(N[(N[(k * a), $MachinePrecision] * 99.0), $MachinePrecision] * k), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq -0.98:\\
\;\;\;\;\frac{1}{\frac{k \cdot k}{a}}\\
\mathbf{elif}\;m \leq 0.65:\\
\;\;\;\;\frac{1}{\mathsf{fma}\left(10 + k, k, 1\right)} \cdot a\\
\mathbf{else}:\\
\;\;\;\;\left(\left(k \cdot a\right) \cdot 99\right) \cdot k\\
\end{array}
\end{array}
if m < -0.97999999999999998Initial 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
associate-*l*N/A
unpow2N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
lower-fma.f64N/A
Applied rewrites37.7%
Applied rewrites38.1%
Taylor expanded in k around inf
Applied rewrites64.8%
if -0.97999999999999998 < m < 0.650000000000000022Initial program 94.5%
lift-/.f64N/A
lift-*.f64N/A
associate-/l*N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f6494.5
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
+-commutativeN/A
lower-+.f6494.6
Applied rewrites94.6%
Taylor expanded in m around 0
lower-/.f64N/A
+-commutativeN/A
*-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
distribute-lft1-inN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
distribute-rgt-inN/A
*-lft-identityN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
lower-+.f6493.5
Applied rewrites93.5%
if 0.650000000000000022 < m Initial program 75.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
associate-*l*N/A
unpow2N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
lower-fma.f64N/A
Applied rewrites3.1%
Taylor expanded in k around 0
Applied rewrites24.9%
Taylor expanded in k around inf
Applied rewrites52.3%
(FPCore (a k m)
:precision binary64
(if (<= m -0.98)
(/ a (* k k))
(if (<= m 0.65)
(* (/ 1.0 (fma (+ 10.0 k) k 1.0)) a)
(* (* (* k a) 99.0) k))))
double code(double a, double k, double m) {
double tmp;
if (m <= -0.98) {
tmp = a / (k * k);
} else if (m <= 0.65) {
tmp = (1.0 / fma((10.0 + k), k, 1.0)) * a;
} else {
tmp = ((k * a) * 99.0) * k;
}
return tmp;
}
function code(a, k, m) tmp = 0.0 if (m <= -0.98) tmp = Float64(a / Float64(k * k)); elseif (m <= 0.65) tmp = Float64(Float64(1.0 / fma(Float64(10.0 + k), k, 1.0)) * a); else tmp = Float64(Float64(Float64(k * a) * 99.0) * k); end return tmp end
code[a_, k_, m_] := If[LessEqual[m, -0.98], N[(a / N[(k * k), $MachinePrecision]), $MachinePrecision], If[LessEqual[m, 0.65], N[(N[(1.0 / N[(N[(10.0 + k), $MachinePrecision] * k + 1.0), $MachinePrecision]), $MachinePrecision] * a), $MachinePrecision], N[(N[(N[(k * a), $MachinePrecision] * 99.0), $MachinePrecision] * k), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq -0.98:\\
\;\;\;\;\frac{a}{k \cdot k}\\
\mathbf{elif}\;m \leq 0.65:\\
\;\;\;\;\frac{1}{\mathsf{fma}\left(10 + k, k, 1\right)} \cdot a\\
\mathbf{else}:\\
\;\;\;\;\left(\left(k \cdot a\right) \cdot 99\right) \cdot k\\
\end{array}
\end{array}
if m < -0.97999999999999998Initial 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
associate-*l*N/A
unpow2N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
lower-fma.f64N/A
Applied rewrites37.7%
Taylor expanded in k around 0
Applied rewrites3.0%
Taylor expanded in k around inf
Applied rewrites64.4%
if -0.97999999999999998 < m < 0.650000000000000022Initial program 94.5%
lift-/.f64N/A
lift-*.f64N/A
associate-/l*N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f6494.5
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
+-commutativeN/A
lower-+.f6494.6
Applied rewrites94.6%
Taylor expanded in m around 0
lower-/.f64N/A
+-commutativeN/A
*-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
distribute-lft1-inN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
distribute-rgt-inN/A
*-lft-identityN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
lower-+.f6493.5
Applied rewrites93.5%
if 0.650000000000000022 < m Initial program 75.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
associate-*l*N/A
unpow2N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
lower-fma.f64N/A
Applied rewrites3.1%
Taylor expanded in k around 0
Applied rewrites24.9%
Taylor expanded in k around inf
Applied rewrites52.3%
(FPCore (a k m) :precision binary64 (if (<= m -0.98) (/ a (* k k)) (if (<= m 0.65) (/ a (fma (+ 10.0 k) k 1.0)) (* (* (* k a) 99.0) k))))
double code(double a, double k, double m) {
double tmp;
if (m <= -0.98) {
tmp = a / (k * k);
} else if (m <= 0.65) {
tmp = a / fma((10.0 + k), k, 1.0);
} else {
tmp = ((k * a) * 99.0) * k;
}
return tmp;
}
function code(a, k, m) tmp = 0.0 if (m <= -0.98) tmp = Float64(a / Float64(k * k)); elseif (m <= 0.65) tmp = Float64(a / fma(Float64(10.0 + k), k, 1.0)); else tmp = Float64(Float64(Float64(k * a) * 99.0) * k); end return tmp end
code[a_, k_, m_] := If[LessEqual[m, -0.98], N[(a / N[(k * k), $MachinePrecision]), $MachinePrecision], If[LessEqual[m, 0.65], N[(a / N[(N[(10.0 + k), $MachinePrecision] * k + 1.0), $MachinePrecision]), $MachinePrecision], N[(N[(N[(k * a), $MachinePrecision] * 99.0), $MachinePrecision] * k), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq -0.98:\\
\;\;\;\;\frac{a}{k \cdot k}\\
\mathbf{elif}\;m \leq 0.65:\\
\;\;\;\;\frac{a}{\mathsf{fma}\left(10 + k, k, 1\right)}\\
\mathbf{else}:\\
\;\;\;\;\left(\left(k \cdot a\right) \cdot 99\right) \cdot k\\
\end{array}
\end{array}
if m < -0.97999999999999998Initial 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
associate-*l*N/A
unpow2N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
lower-fma.f64N/A
Applied rewrites37.7%
Taylor expanded in k around 0
Applied rewrites3.0%
Taylor expanded in k around inf
Applied rewrites64.4%
if -0.97999999999999998 < m < 0.650000000000000022Initial program 94.5%
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
associate-*l*N/A
unpow2N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
lower-fma.f64N/A
Applied rewrites93.4%
if 0.650000000000000022 < m Initial program 75.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
associate-*l*N/A
unpow2N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
lower-fma.f64N/A
Applied rewrites3.1%
Taylor expanded in k around 0
Applied rewrites24.9%
Taylor expanded in k around inf
Applied rewrites52.3%
(FPCore (a k m) :precision binary64 (if (<= m -1.36e-16) (/ a (* k k)) (if (<= m 0.65) (/ a (fma 10.0 k 1.0)) (* (* (* k a) 99.0) k))))
double code(double a, double k, double m) {
double tmp;
if (m <= -1.36e-16) {
tmp = a / (k * k);
} else if (m <= 0.65) {
tmp = a / fma(10.0, k, 1.0);
} else {
tmp = ((k * a) * 99.0) * k;
}
return tmp;
}
function code(a, k, m) tmp = 0.0 if (m <= -1.36e-16) tmp = Float64(a / Float64(k * k)); elseif (m <= 0.65) tmp = Float64(a / fma(10.0, k, 1.0)); else tmp = Float64(Float64(Float64(k * a) * 99.0) * k); end return tmp end
code[a_, k_, m_] := If[LessEqual[m, -1.36e-16], N[(a / N[(k * k), $MachinePrecision]), $MachinePrecision], If[LessEqual[m, 0.65], N[(a / N[(10.0 * k + 1.0), $MachinePrecision]), $MachinePrecision], N[(N[(N[(k * a), $MachinePrecision] * 99.0), $MachinePrecision] * k), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq -1.36 \cdot 10^{-16}:\\
\;\;\;\;\frac{a}{k \cdot k}\\
\mathbf{elif}\;m \leq 0.65:\\
\;\;\;\;\frac{a}{\mathsf{fma}\left(10, k, 1\right)}\\
\mathbf{else}:\\
\;\;\;\;\left(\left(k \cdot a\right) \cdot 99\right) \cdot k\\
\end{array}
\end{array}
if m < -1.3599999999999999e-16Initial 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
associate-*l*N/A
unpow2N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
lower-fma.f64N/A
Applied rewrites39.0%
Taylor expanded in k around 0
Applied rewrites2.9%
Taylor expanded in k around inf
Applied rewrites65.2%
if -1.3599999999999999e-16 < m < 0.650000000000000022Initial 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
associate-*l*N/A
unpow2N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
lower-fma.f64N/A
Applied rewrites93.3%
Taylor expanded in k around 0
Applied rewrites65.4%
if 0.650000000000000022 < m Initial program 75.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
associate-*l*N/A
unpow2N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
lower-fma.f64N/A
Applied rewrites3.1%
Taylor expanded in k around 0
Applied rewrites24.9%
Taylor expanded in k around inf
Applied rewrites52.3%
(FPCore (a k m) :precision binary64 (if (<= m 3.2e-207) (/ a (* k k)) (if (<= m 0.52) (fma (* (fma 99.0 k -10.0) k) a a) (* (* (* k a) 99.0) k))))
double code(double a, double k, double m) {
double tmp;
if (m <= 3.2e-207) {
tmp = a / (k * k);
} else if (m <= 0.52) {
tmp = fma((fma(99.0, k, -10.0) * k), a, a);
} else {
tmp = ((k * a) * 99.0) * k;
}
return tmp;
}
function code(a, k, m) tmp = 0.0 if (m <= 3.2e-207) tmp = Float64(a / Float64(k * k)); elseif (m <= 0.52) tmp = fma(Float64(fma(99.0, k, -10.0) * k), a, a); else tmp = Float64(Float64(Float64(k * a) * 99.0) * k); end return tmp end
code[a_, k_, m_] := If[LessEqual[m, 3.2e-207], N[(a / N[(k * k), $MachinePrecision]), $MachinePrecision], If[LessEqual[m, 0.52], N[(N[(N[(99.0 * k + -10.0), $MachinePrecision] * k), $MachinePrecision] * a + a), $MachinePrecision], N[(N[(N[(k * a), $MachinePrecision] * 99.0), $MachinePrecision] * k), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq 3.2 \cdot 10^{-207}:\\
\;\;\;\;\frac{a}{k \cdot k}\\
\mathbf{elif}\;m \leq 0.52:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(99, k, -10\right) \cdot k, a, a\right)\\
\mathbf{else}:\\
\;\;\;\;\left(\left(k \cdot a\right) \cdot 99\right) \cdot k\\
\end{array}
\end{array}
if m < 3.2000000000000003e-207Initial program 98.6%
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
associate-*l*N/A
unpow2N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
lower-fma.f64N/A
Applied rewrites58.6%
Taylor expanded in k around 0
Applied rewrites17.9%
Taylor expanded in k around inf
Applied rewrites60.0%
if 3.2000000000000003e-207 < m < 0.52000000000000002Initial program 92.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
associate-*l*N/A
unpow2N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
lower-fma.f64N/A
Applied rewrites89.6%
Taylor expanded in k around 0
Applied rewrites59.9%
Applied rewrites59.9%
if 0.52000000000000002 < m Initial program 75.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
associate-*l*N/A
unpow2N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
lower-fma.f64N/A
Applied rewrites3.1%
Taylor expanded in k around 0
Applied rewrites24.9%
Taylor expanded in k around inf
Applied rewrites52.3%
(FPCore (a k m) :precision binary64 (if (<= m 3.2e-207) (/ a (* k k)) (if (<= m 0.52) (fma (* k a) -10.0 a) (* (* (* k a) 99.0) k))))
double code(double a, double k, double m) {
double tmp;
if (m <= 3.2e-207) {
tmp = a / (k * k);
} else if (m <= 0.52) {
tmp = fma((k * a), -10.0, a);
} else {
tmp = ((k * a) * 99.0) * k;
}
return tmp;
}
function code(a, k, m) tmp = 0.0 if (m <= 3.2e-207) tmp = Float64(a / Float64(k * k)); elseif (m <= 0.52) tmp = fma(Float64(k * a), -10.0, a); else tmp = Float64(Float64(Float64(k * a) * 99.0) * k); end return tmp end
code[a_, k_, m_] := If[LessEqual[m, 3.2e-207], N[(a / N[(k * k), $MachinePrecision]), $MachinePrecision], If[LessEqual[m, 0.52], N[(N[(k * a), $MachinePrecision] * -10.0 + a), $MachinePrecision], N[(N[(N[(k * a), $MachinePrecision] * 99.0), $MachinePrecision] * k), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq 3.2 \cdot 10^{-207}:\\
\;\;\;\;\frac{a}{k \cdot k}\\
\mathbf{elif}\;m \leq 0.52:\\
\;\;\;\;\mathsf{fma}\left(k \cdot a, -10, a\right)\\
\mathbf{else}:\\
\;\;\;\;\left(\left(k \cdot a\right) \cdot 99\right) \cdot k\\
\end{array}
\end{array}
if m < 3.2000000000000003e-207Initial program 98.6%
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
associate-*l*N/A
unpow2N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
lower-fma.f64N/A
Applied rewrites58.6%
Taylor expanded in k around 0
Applied rewrites17.9%
Taylor expanded in k around inf
Applied rewrites60.0%
if 3.2000000000000003e-207 < m < 0.52000000000000002Initial program 92.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
associate-*l*N/A
unpow2N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
lower-fma.f64N/A
Applied rewrites89.6%
Taylor expanded in k around 0
Applied rewrites59.0%
if 0.52000000000000002 < m Initial program 75.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
associate-*l*N/A
unpow2N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
lower-fma.f64N/A
Applied rewrites3.1%
Taylor expanded in k around 0
Applied rewrites24.9%
Taylor expanded in k around inf
Applied rewrites52.3%
Final simplification57.5%
(FPCore (a k m) :precision binary64 (if (<= m 0.52) (* 1.0 a) (* (* (* k a) 99.0) k)))
double code(double a, double k, double m) {
double tmp;
if (m <= 0.52) {
tmp = 1.0 * a;
} else {
tmp = ((k * a) * 99.0) * k;
}
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.52d0) then
tmp = 1.0d0 * a
else
tmp = ((k * a) * 99.0d0) * k
end if
code = tmp
end function
public static double code(double a, double k, double m) {
double tmp;
if (m <= 0.52) {
tmp = 1.0 * a;
} else {
tmp = ((k * a) * 99.0) * k;
}
return tmp;
}
def code(a, k, m): tmp = 0 if m <= 0.52: tmp = 1.0 * a else: tmp = ((k * a) * 99.0) * k return tmp
function code(a, k, m) tmp = 0.0 if (m <= 0.52) tmp = Float64(1.0 * a); else tmp = Float64(Float64(Float64(k * a) * 99.0) * k); end return tmp end
function tmp_2 = code(a, k, m) tmp = 0.0; if (m <= 0.52) tmp = 1.0 * a; else tmp = ((k * a) * 99.0) * k; end tmp_2 = tmp; end
code[a_, k_, m_] := If[LessEqual[m, 0.52], N[(1.0 * a), $MachinePrecision], N[(N[(N[(k * a), $MachinePrecision] * 99.0), $MachinePrecision] * k), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq 0.52:\\
\;\;\;\;1 \cdot a\\
\mathbf{else}:\\
\;\;\;\;\left(\left(k \cdot a\right) \cdot 99\right) \cdot k\\
\end{array}
\end{array}
if m < 0.52000000000000002Initial program 97.3%
lift-/.f64N/A
lift-*.f64N/A
associate-/l*N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f6497.3
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
+-commutativeN/A
lower-+.f6497.3
Applied rewrites97.3%
Taylor expanded in m around 0
lower-/.f64N/A
+-commutativeN/A
*-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
distribute-lft1-inN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
distribute-rgt-inN/A
*-lft-identityN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
lower-+.f6465.1
Applied rewrites65.1%
Taylor expanded in k around 0
Applied rewrites26.8%
if 0.52000000000000002 < m Initial program 75.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
associate-*l*N/A
unpow2N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
lower-fma.f64N/A
Applied rewrites3.1%
Taylor expanded in k around 0
Applied rewrites24.9%
Taylor expanded in k around inf
Applied rewrites52.3%
(FPCore (a k m) :precision binary64 (if (<= m 320000.0) (* 1.0 a) (* (* -10.0 a) k)))
double code(double a, double k, double m) {
double tmp;
if (m <= 320000.0) {
tmp = 1.0 * a;
} else {
tmp = (-10.0 * a) * k;
}
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 <= 320000.0d0) then
tmp = 1.0d0 * a
else
tmp = ((-10.0d0) * a) * k
end if
code = tmp
end function
public static double code(double a, double k, double m) {
double tmp;
if (m <= 320000.0) {
tmp = 1.0 * a;
} else {
tmp = (-10.0 * a) * k;
}
return tmp;
}
def code(a, k, m): tmp = 0 if m <= 320000.0: tmp = 1.0 * a else: tmp = (-10.0 * a) * k return tmp
function code(a, k, m) tmp = 0.0 if (m <= 320000.0) tmp = Float64(1.0 * a); else tmp = Float64(Float64(-10.0 * a) * k); end return tmp end
function tmp_2 = code(a, k, m) tmp = 0.0; if (m <= 320000.0) tmp = 1.0 * a; else tmp = (-10.0 * a) * k; end tmp_2 = tmp; end
code[a_, k_, m_] := If[LessEqual[m, 320000.0], N[(1.0 * a), $MachinePrecision], N[(N[(-10.0 * a), $MachinePrecision] * k), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq 320000:\\
\;\;\;\;1 \cdot a\\
\mathbf{else}:\\
\;\;\;\;\left(-10 \cdot a\right) \cdot k\\
\end{array}
\end{array}
if m < 3.2e5Initial program 96.8%
lift-/.f64N/A
lift-*.f64N/A
associate-/l*N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f6496.8
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
+-commutativeN/A
lower-+.f6496.8
Applied rewrites96.8%
Taylor expanded in m around 0
lower-/.f64N/A
+-commutativeN/A
*-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
distribute-lft1-inN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
distribute-rgt-inN/A
*-lft-identityN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
lower-+.f6463.7
Applied rewrites63.7%
Taylor expanded in k around 0
Applied rewrites26.3%
if 3.2e5 < m Initial program 76.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
associate-*l*N/A
unpow2N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
lower-fma.f64N/A
Applied rewrites3.1%
Taylor expanded in k around 0
Applied rewrites4.8%
Taylor expanded in k around inf
Applied rewrites20.3%
(FPCore (a k m) :precision binary64 (* 1.0 a))
double code(double a, double k, double m) {
return 1.0 * a;
}
real(8) function code(a, k, m)
real(8), intent (in) :: a
real(8), intent (in) :: k
real(8), intent (in) :: m
code = 1.0d0 * a
end function
public static double code(double a, double k, double m) {
return 1.0 * a;
}
def code(a, k, m): return 1.0 * a
function code(a, k, m) return Float64(1.0 * a) end
function tmp = code(a, k, m) tmp = 1.0 * a; end
code[a_, k_, m_] := N[(1.0 * a), $MachinePrecision]
\begin{array}{l}
\\
1 \cdot a
\end{array}
Initial program 90.7%
lift-/.f64N/A
lift-*.f64N/A
associate-/l*N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f6490.7
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
+-commutativeN/A
lower-+.f6490.7
Applied rewrites90.7%
Taylor expanded in m around 0
lower-/.f64N/A
+-commutativeN/A
*-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
distribute-lft1-inN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
distribute-rgt-inN/A
*-lft-identityN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
lower-+.f6446.0
Applied rewrites46.0%
Taylor expanded in k around 0
Applied rewrites19.7%
herbie shell --seed 2024297
(FPCore (a k m)
:name "Falkner and Boettcher, Appendix A"
:precision binary64
(/ (* a (pow k m)) (+ (+ 1.0 (* 10.0 k)) (* k k))))