
(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 13 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 (<= m -1.45e-8)
(/ (* a (pow k m)) (* k k))
(if (<= m 2.65e-14)
(* (pow (fma k k (fma 10.0 k 1.0)) -1.0) a)
(* (pow k m) a))))
double code(double a, double k, double m) {
double tmp;
if (m <= -1.45e-8) {
tmp = (a * pow(k, m)) / (k * k);
} else if (m <= 2.65e-14) {
tmp = pow(fma(k, k, fma(10.0, k, 1.0)), -1.0) * a;
} else {
tmp = pow(k, m) * a;
}
return tmp;
}
function code(a, k, m) tmp = 0.0 if (m <= -1.45e-8) tmp = Float64(Float64(a * (k ^ m)) / Float64(k * k)); elseif (m <= 2.65e-14) tmp = Float64((fma(k, k, fma(10.0, k, 1.0)) ^ -1.0) * a); else tmp = Float64((k ^ m) * a); end return tmp end
code[a_, k_, m_] := If[LessEqual[m, -1.45e-8], N[(N[(a * N[Power[k, m], $MachinePrecision]), $MachinePrecision] / N[(k * k), $MachinePrecision]), $MachinePrecision], If[LessEqual[m, 2.65e-14], N[(N[Power[N[(k * k + N[(10.0 * k + 1.0), $MachinePrecision]), $MachinePrecision], -1.0], $MachinePrecision] * a), $MachinePrecision], N[(N[Power[k, m], $MachinePrecision] * a), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq -1.45 \cdot 10^{-8}:\\
\;\;\;\;\frac{a \cdot {k}^{m}}{k \cdot k}\\
\mathbf{elif}\;m \leq 2.65 \cdot 10^{-14}:\\
\;\;\;\;{\left(\mathsf{fma}\left(k, k, \mathsf{fma}\left(10, k, 1\right)\right)\right)}^{-1} \cdot a\\
\mathbf{else}:\\
\;\;\;\;{k}^{m} \cdot a\\
\end{array}
\end{array}
if m < -1.4500000000000001e-8Initial program 100.0%
Taylor expanded in k around inf
unpow2N/A
lower-*.f64100.0
Applied rewrites100.0%
if -1.4500000000000001e-8 < m < 2.6500000000000001e-14Initial program 95.6%
lift-/.f64N/A
lift-*.f64N/A
associate-/l*N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f6495.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-+.f6495.6
Applied rewrites95.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
+-commutativeN/A
distribute-lft1-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
lower-+.f6495.6
Applied rewrites95.6%
Applied rewrites95.7%
if 2.6500000000000001e-14 < m Initial program 81.3%
lift-/.f64N/A
lift-*.f64N/A
associate-/l*N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f6481.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-+.f6481.3
Applied rewrites81.3%
Taylor expanded in k around 0
lower-pow.f64100.0
Applied rewrites100.0%
Final simplification98.5%
(FPCore (a k m) :precision binary64 (if (<= (/ (* a (pow k m)) (+ (+ 1.0 (* 10.0 k)) (* k k))) INFINITY) (* (/ (pow k m) (fma (+ k 10.0) 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 (((a * pow(k, m)) / ((1.0 + (10.0 * k)) + (k * k))) <= ((double) INFINITY)) {
tmp = (pow(k, m) / fma((k + 10.0), 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(a * (k ^ m)) / Float64(Float64(1.0 + Float64(10.0 * k)) + Float64(k * k))) <= Inf) tmp = Float64(Float64((k ^ m) / fma(Float64(k + 10.0), 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[(a * N[Power[k, m], $MachinePrecision]), $MachinePrecision] / N[(N[(1.0 + N[(10.0 * k), $MachinePrecision]), $MachinePrecision] + N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], Infinity], N[(N[(N[Power[k, m], $MachinePrecision] / N[(N[(k + 10.0), $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{a \cdot {k}^{m}}{\left(1 + 10 \cdot k\right) + k \cdot k} \leq \infty:\\
\;\;\;\;\frac{{k}^{m}}{\mathsf{fma}\left(k + 10, 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.4%
lift-/.f64N/A
lift-*.f64N/A
associate-/l*N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f6498.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-+.f6498.4
Applied rewrites98.4%
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
associate-+r+N/A
+-commutativeN/A
+-commutativeN/A
associate-+l+N/A
+-commutativeN/A
associate-+l+N/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
associate-*r*N/A
unpow2N/A
associate-+l+N/A
distribute-lft1-inN/A
+-commutativeN/A
unpow2N/A
associate-*r*N/A
lower-fma.f64N/A
Applied rewrites1.6%
Taylor expanded in k around 0
Applied rewrites73.1%
Applied rewrites100.0%
(FPCore (a k m) :precision binary64 (if (or (<= m -6e-7) (not (<= m 2.65e-14))) (* (pow k m) a) (* (pow (fma k k (fma 10.0 k 1.0)) -1.0) a)))
double code(double a, double k, double m) {
double tmp;
if ((m <= -6e-7) || !(m <= 2.65e-14)) {
tmp = pow(k, m) * a;
} else {
tmp = pow(fma(k, k, fma(10.0, k, 1.0)), -1.0) * a;
}
return tmp;
}
function code(a, k, m) tmp = 0.0 if ((m <= -6e-7) || !(m <= 2.65e-14)) tmp = Float64((k ^ m) * a); else tmp = Float64((fma(k, k, fma(10.0, k, 1.0)) ^ -1.0) * a); end return tmp end
code[a_, k_, m_] := If[Or[LessEqual[m, -6e-7], N[Not[LessEqual[m, 2.65e-14]], $MachinePrecision]], N[(N[Power[k, m], $MachinePrecision] * a), $MachinePrecision], N[(N[Power[N[(k * k + N[(10.0 * k + 1.0), $MachinePrecision]), $MachinePrecision], -1.0], $MachinePrecision] * a), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq -6 \cdot 10^{-7} \lor \neg \left(m \leq 2.65 \cdot 10^{-14}\right):\\
\;\;\;\;{k}^{m} \cdot a\\
\mathbf{else}:\\
\;\;\;\;{\left(\mathsf{fma}\left(k, k, \mathsf{fma}\left(10, k, 1\right)\right)\right)}^{-1} \cdot a\\
\end{array}
\end{array}
if m < -5.9999999999999997e-7 or 2.6500000000000001e-14 < m Initial program 91.5%
lift-/.f64N/A
lift-*.f64N/A
associate-/l*N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f6491.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-+.f6491.5
Applied rewrites91.5%
Taylor expanded in k around 0
lower-pow.f64100.0
Applied rewrites100.0%
if -5.9999999999999997e-7 < m < 2.6500000000000001e-14Initial program 95.7%
lift-/.f64N/A
lift-*.f64N/A
associate-/l*N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f6495.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-+.f6495.7
Applied rewrites95.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
+-commutativeN/A
distribute-lft1-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
lower-+.f6495.1
Applied rewrites95.1%
Applied rewrites95.1%
Final simplification98.2%
(FPCore (a k m)
:precision binary64
(if (<= m -0.245)
(/ (* (/ a (* k k)) 99.0) (* k k))
(if (<= m 1.3)
(* (pow (fma k k (fma 10.0 k 1.0)) -1.0) a)
(* (* (* a k) k) 99.0))))
double code(double a, double k, double m) {
double tmp;
if (m <= -0.245) {
tmp = ((a / (k * k)) * 99.0) / (k * k);
} else if (m <= 1.3) {
tmp = pow(fma(k, k, fma(10.0, k, 1.0)), -1.0) * a;
} else {
tmp = ((a * k) * k) * 99.0;
}
return tmp;
}
function code(a, k, m) tmp = 0.0 if (m <= -0.245) tmp = Float64(Float64(Float64(a / Float64(k * k)) * 99.0) / Float64(k * k)); elseif (m <= 1.3) tmp = Float64((fma(k, k, fma(10.0, k, 1.0)) ^ -1.0) * a); else tmp = Float64(Float64(Float64(a * k) * k) * 99.0); end return tmp end
code[a_, k_, m_] := If[LessEqual[m, -0.245], N[(N[(N[(a / N[(k * k), $MachinePrecision]), $MachinePrecision] * 99.0), $MachinePrecision] / N[(k * k), $MachinePrecision]), $MachinePrecision], If[LessEqual[m, 1.3], N[(N[Power[N[(k * k + N[(10.0 * k + 1.0), $MachinePrecision]), $MachinePrecision], -1.0], $MachinePrecision] * a), $MachinePrecision], N[(N[(N[(a * k), $MachinePrecision] * k), $MachinePrecision] * 99.0), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq -0.245:\\
\;\;\;\;\frac{\frac{a}{k \cdot k} \cdot 99}{k \cdot k}\\
\mathbf{elif}\;m \leq 1.3:\\
\;\;\;\;{\left(\mathsf{fma}\left(k, k, \mathsf{fma}\left(10, k, 1\right)\right)\right)}^{-1} \cdot a\\
\mathbf{else}:\\
\;\;\;\;\left(\left(a \cdot k\right) \cdot k\right) \cdot 99\\
\end{array}
\end{array}
if m < -0.245Initial program 100.0%
Taylor expanded in m around 0
lower-/.f64N/A
associate-+r+N/A
+-commutativeN/A
+-commutativeN/A
associate-+l+N/A
+-commutativeN/A
associate-+l+N/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
associate-*r*N/A
unpow2N/A
associate-+l+N/A
distribute-lft1-inN/A
+-commutativeN/A
unpow2N/A
associate-*r*N/A
lower-fma.f64N/A
Applied rewrites32.5%
Taylor expanded in k around inf
Applied rewrites64.9%
Taylor expanded in k around 0
Applied rewrites74.1%
if -0.245 < m < 1.30000000000000004Initial program 95.7%
lift-/.f64N/A
lift-*.f64N/A
associate-/l*N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f6495.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-+.f6495.7
Applied rewrites95.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
+-commutativeN/A
distribute-lft1-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
lower-+.f6494.8
Applied rewrites94.8%
Applied rewrites94.8%
if 1.30000000000000004 < m Initial program 81.1%
Taylor expanded in m around 0
lower-/.f64N/A
associate-+r+N/A
+-commutativeN/A
+-commutativeN/A
associate-+l+N/A
+-commutativeN/A
associate-+l+N/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
associate-*r*N/A
unpow2N/A
associate-+l+N/A
distribute-lft1-inN/A
+-commutativeN/A
unpow2N/A
associate-*r*N/A
lower-fma.f64N/A
Applied rewrites2.8%
Taylor expanded in k around 0
Applied rewrites19.5%
Taylor expanded in k around inf
Applied rewrites38.4%
Final simplification71.2%
(FPCore (a k m)
:precision binary64
(if (<= m -0.245)
(* (pow (* k k) -1.0) a)
(if (<= m 1.3)
(* (pow (fma k k (fma 10.0 k 1.0)) -1.0) a)
(* (* (* a k) k) 99.0))))
double code(double a, double k, double m) {
double tmp;
if (m <= -0.245) {
tmp = pow((k * k), -1.0) * a;
} else if (m <= 1.3) {
tmp = pow(fma(k, k, fma(10.0, k, 1.0)), -1.0) * a;
} else {
tmp = ((a * k) * k) * 99.0;
}
return tmp;
}
function code(a, k, m) tmp = 0.0 if (m <= -0.245) tmp = Float64((Float64(k * k) ^ -1.0) * a); elseif (m <= 1.3) tmp = Float64((fma(k, k, fma(10.0, k, 1.0)) ^ -1.0) * a); else tmp = Float64(Float64(Float64(a * k) * k) * 99.0); end return tmp end
code[a_, k_, m_] := If[LessEqual[m, -0.245], N[(N[Power[N[(k * k), $MachinePrecision], -1.0], $MachinePrecision] * a), $MachinePrecision], If[LessEqual[m, 1.3], N[(N[Power[N[(k * k + N[(10.0 * k + 1.0), $MachinePrecision]), $MachinePrecision], -1.0], $MachinePrecision] * a), $MachinePrecision], N[(N[(N[(a * k), $MachinePrecision] * k), $MachinePrecision] * 99.0), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq -0.245:\\
\;\;\;\;{\left(k \cdot k\right)}^{-1} \cdot a\\
\mathbf{elif}\;m \leq 1.3:\\
\;\;\;\;{\left(\mathsf{fma}\left(k, k, \mathsf{fma}\left(10, k, 1\right)\right)\right)}^{-1} \cdot a\\
\mathbf{else}:\\
\;\;\;\;\left(\left(a \cdot k\right) \cdot k\right) \cdot 99\\
\end{array}
\end{array}
if m < -0.245Initial program 100.0%
lift-/.f64N/A
lift-*.f64N/A
associate-/l*N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f64100.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-+.f64100.0
Applied rewrites100.0%
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
+-commutativeN/A
distribute-lft1-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
lower-+.f6432.5
Applied rewrites32.5%
Taylor expanded in k around inf
Applied rewrites56.3%
if -0.245 < m < 1.30000000000000004Initial program 95.7%
lift-/.f64N/A
lift-*.f64N/A
associate-/l*N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f6495.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-+.f6495.7
Applied rewrites95.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
+-commutativeN/A
distribute-lft1-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
lower-+.f6494.8
Applied rewrites94.8%
Applied rewrites94.8%
if 1.30000000000000004 < m Initial program 81.1%
Taylor expanded in m around 0
lower-/.f64N/A
associate-+r+N/A
+-commutativeN/A
+-commutativeN/A
associate-+l+N/A
+-commutativeN/A
associate-+l+N/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
associate-*r*N/A
unpow2N/A
associate-+l+N/A
distribute-lft1-inN/A
+-commutativeN/A
unpow2N/A
associate-*r*N/A
lower-fma.f64N/A
Applied rewrites2.8%
Taylor expanded in k around 0
Applied rewrites19.5%
Taylor expanded in k around inf
Applied rewrites38.4%
Final simplification65.0%
(FPCore (a k m)
:precision binary64
(if (<= m -0.245)
(* (pow (* k k) -1.0) a)
(if (<= m 1.3)
(* (pow (fma (+ 10.0 k) k 1.0) -1.0) a)
(* (* (* a k) k) 99.0))))
double code(double a, double k, double m) {
double tmp;
if (m <= -0.245) {
tmp = pow((k * k), -1.0) * a;
} else if (m <= 1.3) {
tmp = pow(fma((10.0 + k), k, 1.0), -1.0) * a;
} else {
tmp = ((a * k) * k) * 99.0;
}
return tmp;
}
function code(a, k, m) tmp = 0.0 if (m <= -0.245) tmp = Float64((Float64(k * k) ^ -1.0) * a); elseif (m <= 1.3) tmp = Float64((fma(Float64(10.0 + k), k, 1.0) ^ -1.0) * a); else tmp = Float64(Float64(Float64(a * k) * k) * 99.0); end return tmp end
code[a_, k_, m_] := If[LessEqual[m, -0.245], N[(N[Power[N[(k * k), $MachinePrecision], -1.0], $MachinePrecision] * a), $MachinePrecision], If[LessEqual[m, 1.3], N[(N[Power[N[(N[(10.0 + k), $MachinePrecision] * k + 1.0), $MachinePrecision], -1.0], $MachinePrecision] * a), $MachinePrecision], N[(N[(N[(a * k), $MachinePrecision] * k), $MachinePrecision] * 99.0), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq -0.245:\\
\;\;\;\;{\left(k \cdot k\right)}^{-1} \cdot a\\
\mathbf{elif}\;m \leq 1.3:\\
\;\;\;\;{\left(\mathsf{fma}\left(10 + k, k, 1\right)\right)}^{-1} \cdot a\\
\mathbf{else}:\\
\;\;\;\;\left(\left(a \cdot k\right) \cdot k\right) \cdot 99\\
\end{array}
\end{array}
if m < -0.245Initial program 100.0%
lift-/.f64N/A
lift-*.f64N/A
associate-/l*N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f64100.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-+.f64100.0
Applied rewrites100.0%
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
+-commutativeN/A
distribute-lft1-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
lower-+.f6432.5
Applied rewrites32.5%
Taylor expanded in k around inf
Applied rewrites56.3%
if -0.245 < m < 1.30000000000000004Initial program 95.7%
lift-/.f64N/A
lift-*.f64N/A
associate-/l*N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f6495.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-+.f6495.7
Applied rewrites95.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
+-commutativeN/A
distribute-lft1-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
lower-+.f6494.8
Applied rewrites94.8%
if 1.30000000000000004 < m Initial program 81.1%
Taylor expanded in m around 0
lower-/.f64N/A
associate-+r+N/A
+-commutativeN/A
+-commutativeN/A
associate-+l+N/A
+-commutativeN/A
associate-+l+N/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
associate-*r*N/A
unpow2N/A
associate-+l+N/A
distribute-lft1-inN/A
+-commutativeN/A
unpow2N/A
associate-*r*N/A
lower-fma.f64N/A
Applied rewrites2.8%
Taylor expanded in k around 0
Applied rewrites19.5%
Taylor expanded in k around inf
Applied rewrites38.4%
Final simplification65.0%
(FPCore (a k m) :precision binary64 (if (<= m -0.245) (* (pow (* k k) -1.0) a) (if (<= m 1.3) (/ a (fma (+ 10.0 k) k 1.0)) (* (* (* a k) k) 99.0))))
double code(double a, double k, double m) {
double tmp;
if (m <= -0.245) {
tmp = pow((k * k), -1.0) * a;
} else if (m <= 1.3) {
tmp = a / fma((10.0 + k), k, 1.0);
} else {
tmp = ((a * k) * k) * 99.0;
}
return tmp;
}
function code(a, k, m) tmp = 0.0 if (m <= -0.245) tmp = Float64((Float64(k * k) ^ -1.0) * a); elseif (m <= 1.3) tmp = Float64(a / fma(Float64(10.0 + k), k, 1.0)); else tmp = Float64(Float64(Float64(a * k) * k) * 99.0); end return tmp end
code[a_, k_, m_] := If[LessEqual[m, -0.245], N[(N[Power[N[(k * k), $MachinePrecision], -1.0], $MachinePrecision] * a), $MachinePrecision], If[LessEqual[m, 1.3], N[(a / N[(N[(10.0 + k), $MachinePrecision] * k + 1.0), $MachinePrecision]), $MachinePrecision], N[(N[(N[(a * k), $MachinePrecision] * k), $MachinePrecision] * 99.0), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq -0.245:\\
\;\;\;\;{\left(k \cdot k\right)}^{-1} \cdot a\\
\mathbf{elif}\;m \leq 1.3:\\
\;\;\;\;\frac{a}{\mathsf{fma}\left(10 + k, k, 1\right)}\\
\mathbf{else}:\\
\;\;\;\;\left(\left(a \cdot k\right) \cdot k\right) \cdot 99\\
\end{array}
\end{array}
if m < -0.245Initial program 100.0%
lift-/.f64N/A
lift-*.f64N/A
associate-/l*N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f64100.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-+.f64100.0
Applied rewrites100.0%
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
+-commutativeN/A
distribute-lft1-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
lower-+.f6432.5
Applied rewrites32.5%
Taylor expanded in k around inf
Applied rewrites56.3%
if -0.245 < m < 1.30000000000000004Initial program 95.7%
Taylor expanded in m around 0
lower-/.f64N/A
associate-+r+N/A
+-commutativeN/A
+-commutativeN/A
associate-+l+N/A
+-commutativeN/A
associate-+l+N/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
associate-*r*N/A
unpow2N/A
associate-+l+N/A
distribute-lft1-inN/A
+-commutativeN/A
unpow2N/A
associate-*r*N/A
lower-fma.f64N/A
Applied rewrites94.8%
if 1.30000000000000004 < m Initial program 81.1%
Taylor expanded in m around 0
lower-/.f64N/A
associate-+r+N/A
+-commutativeN/A
+-commutativeN/A
associate-+l+N/A
+-commutativeN/A
associate-+l+N/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
associate-*r*N/A
unpow2N/A
associate-+l+N/A
distribute-lft1-inN/A
+-commutativeN/A
unpow2N/A
associate-*r*N/A
lower-fma.f64N/A
Applied rewrites2.8%
Taylor expanded in k around 0
Applied rewrites19.5%
Taylor expanded in k around inf
Applied rewrites38.4%
Final simplification65.0%
(FPCore (a k m) :precision binary64 (if (<= m -0.245) (/ a (* k k)) (if (<= m 1.3) (/ a (fma (+ 10.0 k) k 1.0)) (* (* (* a k) k) 99.0))))
double code(double a, double k, double m) {
double tmp;
if (m <= -0.245) {
tmp = a / (k * k);
} else if (m <= 1.3) {
tmp = a / fma((10.0 + k), k, 1.0);
} else {
tmp = ((a * k) * k) * 99.0;
}
return tmp;
}
function code(a, k, m) tmp = 0.0 if (m <= -0.245) tmp = Float64(a / Float64(k * k)); elseif (m <= 1.3) tmp = Float64(a / fma(Float64(10.0 + k), k, 1.0)); else tmp = Float64(Float64(Float64(a * k) * k) * 99.0); end return tmp end
code[a_, k_, m_] := If[LessEqual[m, -0.245], N[(a / N[(k * k), $MachinePrecision]), $MachinePrecision], If[LessEqual[m, 1.3], N[(a / N[(N[(10.0 + k), $MachinePrecision] * k + 1.0), $MachinePrecision]), $MachinePrecision], N[(N[(N[(a * k), $MachinePrecision] * k), $MachinePrecision] * 99.0), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq -0.245:\\
\;\;\;\;\frac{a}{k \cdot k}\\
\mathbf{elif}\;m \leq 1.3:\\
\;\;\;\;\frac{a}{\mathsf{fma}\left(10 + k, k, 1\right)}\\
\mathbf{else}:\\
\;\;\;\;\left(\left(a \cdot k\right) \cdot k\right) \cdot 99\\
\end{array}
\end{array}
if m < -0.245Initial program 100.0%
Taylor expanded in m around 0
lower-/.f64N/A
associate-+r+N/A
+-commutativeN/A
+-commutativeN/A
associate-+l+N/A
+-commutativeN/A
associate-+l+N/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
associate-*r*N/A
unpow2N/A
associate-+l+N/A
distribute-lft1-inN/A
+-commutativeN/A
unpow2N/A
associate-*r*N/A
lower-fma.f64N/A
Applied rewrites32.5%
Taylor expanded in k around inf
Applied rewrites55.3%
if -0.245 < m < 1.30000000000000004Initial program 95.7%
Taylor expanded in m around 0
lower-/.f64N/A
associate-+r+N/A
+-commutativeN/A
+-commutativeN/A
associate-+l+N/A
+-commutativeN/A
associate-+l+N/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
associate-*r*N/A
unpow2N/A
associate-+l+N/A
distribute-lft1-inN/A
+-commutativeN/A
unpow2N/A
associate-*r*N/A
lower-fma.f64N/A
Applied rewrites94.8%
if 1.30000000000000004 < m Initial program 81.1%
Taylor expanded in m around 0
lower-/.f64N/A
associate-+r+N/A
+-commutativeN/A
+-commutativeN/A
associate-+l+N/A
+-commutativeN/A
associate-+l+N/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
associate-*r*N/A
unpow2N/A
associate-+l+N/A
distribute-lft1-inN/A
+-commutativeN/A
unpow2N/A
associate-*r*N/A
lower-fma.f64N/A
Applied rewrites2.8%
Taylor expanded in k around 0
Applied rewrites19.5%
Taylor expanded in k around inf
Applied rewrites38.4%
(FPCore (a k m) :precision binary64 (if (<= m -4.5e-120) (/ a (* k k)) (if (<= m 1.65) (/ a (fma 10.0 k 1.0)) (* (* (* a k) k) 99.0))))
double code(double a, double k, double m) {
double tmp;
if (m <= -4.5e-120) {
tmp = a / (k * k);
} else if (m <= 1.65) {
tmp = a / fma(10.0, k, 1.0);
} else {
tmp = ((a * k) * k) * 99.0;
}
return tmp;
}
function code(a, k, m) tmp = 0.0 if (m <= -4.5e-120) tmp = Float64(a / Float64(k * k)); elseif (m <= 1.65) tmp = Float64(a / fma(10.0, k, 1.0)); else tmp = Float64(Float64(Float64(a * k) * k) * 99.0); end return tmp end
code[a_, k_, m_] := If[LessEqual[m, -4.5e-120], N[(a / N[(k * k), $MachinePrecision]), $MachinePrecision], If[LessEqual[m, 1.65], N[(a / N[(10.0 * k + 1.0), $MachinePrecision]), $MachinePrecision], N[(N[(N[(a * k), $MachinePrecision] * k), $MachinePrecision] * 99.0), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq -4.5 \cdot 10^{-120}:\\
\;\;\;\;\frac{a}{k \cdot k}\\
\mathbf{elif}\;m \leq 1.65:\\
\;\;\;\;\frac{a}{\mathsf{fma}\left(10, k, 1\right)}\\
\mathbf{else}:\\
\;\;\;\;\left(\left(a \cdot k\right) \cdot k\right) \cdot 99\\
\end{array}
\end{array}
if m < -4.5e-120Initial program 100.0%
Taylor expanded in m around 0
lower-/.f64N/A
associate-+r+N/A
+-commutativeN/A
+-commutativeN/A
associate-+l+N/A
+-commutativeN/A
associate-+l+N/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
associate-*r*N/A
unpow2N/A
associate-+l+N/A
distribute-lft1-inN/A
+-commutativeN/A
unpow2N/A
associate-*r*N/A
lower-fma.f64N/A
Applied rewrites39.9%
Taylor expanded in k around inf
Applied rewrites57.2%
if -4.5e-120 < m < 1.6499999999999999Initial program 95.1%
Taylor expanded in m around 0
lower-/.f64N/A
associate-+r+N/A
+-commutativeN/A
+-commutativeN/A
associate-+l+N/A
+-commutativeN/A
associate-+l+N/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
associate-*r*N/A
unpow2N/A
associate-+l+N/A
distribute-lft1-inN/A
+-commutativeN/A
unpow2N/A
associate-*r*N/A
lower-fma.f64N/A
Applied rewrites94.8%
Taylor expanded in k around 0
Applied rewrites71.4%
if 1.6499999999999999 < m Initial program 81.1%
Taylor expanded in m around 0
lower-/.f64N/A
associate-+r+N/A
+-commutativeN/A
+-commutativeN/A
associate-+l+N/A
+-commutativeN/A
associate-+l+N/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
associate-*r*N/A
unpow2N/A
associate-+l+N/A
distribute-lft1-inN/A
+-commutativeN/A
unpow2N/A
associate-*r*N/A
lower-fma.f64N/A
Applied rewrites2.8%
Taylor expanded in k around 0
Applied rewrites19.5%
Taylor expanded in k around inf
Applied rewrites38.4%
(FPCore (a k m) :precision binary64 (if (<= m -2.35e-123) (/ a (* k k)) (if (<= m 0.45) (* 1.0 a) (* (* (* a k) k) 99.0))))
double code(double a, double k, double m) {
double tmp;
if (m <= -2.35e-123) {
tmp = a / (k * k);
} else if (m <= 0.45) {
tmp = 1.0 * a;
} else {
tmp = ((a * k) * k) * 99.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 <= (-2.35d-123)) then
tmp = a / (k * k)
else if (m <= 0.45d0) then
tmp = 1.0d0 * a
else
tmp = ((a * k) * k) * 99.0d0
end if
code = tmp
end function
public static double code(double a, double k, double m) {
double tmp;
if (m <= -2.35e-123) {
tmp = a / (k * k);
} else if (m <= 0.45) {
tmp = 1.0 * a;
} else {
tmp = ((a * k) * k) * 99.0;
}
return tmp;
}
def code(a, k, m): tmp = 0 if m <= -2.35e-123: tmp = a / (k * k) elif m <= 0.45: tmp = 1.0 * a else: tmp = ((a * k) * k) * 99.0 return tmp
function code(a, k, m) tmp = 0.0 if (m <= -2.35e-123) tmp = Float64(a / Float64(k * k)); elseif (m <= 0.45) tmp = Float64(1.0 * a); else tmp = Float64(Float64(Float64(a * k) * k) * 99.0); end return tmp end
function tmp_2 = code(a, k, m) tmp = 0.0; if (m <= -2.35e-123) tmp = a / (k * k); elseif (m <= 0.45) tmp = 1.0 * a; else tmp = ((a * k) * k) * 99.0; end tmp_2 = tmp; end
code[a_, k_, m_] := If[LessEqual[m, -2.35e-123], N[(a / N[(k * k), $MachinePrecision]), $MachinePrecision], If[LessEqual[m, 0.45], N[(1.0 * a), $MachinePrecision], N[(N[(N[(a * k), $MachinePrecision] * k), $MachinePrecision] * 99.0), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq -2.35 \cdot 10^{-123}:\\
\;\;\;\;\frac{a}{k \cdot k}\\
\mathbf{elif}\;m \leq 0.45:\\
\;\;\;\;1 \cdot a\\
\mathbf{else}:\\
\;\;\;\;\left(\left(a \cdot k\right) \cdot k\right) \cdot 99\\
\end{array}
\end{array}
if m < -2.3500000000000001e-123Initial program 100.0%
Taylor expanded in m around 0
lower-/.f64N/A
associate-+r+N/A
+-commutativeN/A
+-commutativeN/A
associate-+l+N/A
+-commutativeN/A
associate-+l+N/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
associate-*r*N/A
unpow2N/A
associate-+l+N/A
distribute-lft1-inN/A
+-commutativeN/A
unpow2N/A
associate-*r*N/A
lower-fma.f64N/A
Applied rewrites39.9%
Taylor expanded in k around inf
Applied rewrites57.2%
if -2.3500000000000001e-123 < m < 0.450000000000000011Initial program 95.1%
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.f6457.5
Applied rewrites57.5%
Taylor expanded in m around 0
Applied rewrites57.1%
if 0.450000000000000011 < m Initial program 81.1%
Taylor expanded in m around 0
lower-/.f64N/A
associate-+r+N/A
+-commutativeN/A
+-commutativeN/A
associate-+l+N/A
+-commutativeN/A
associate-+l+N/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
associate-*r*N/A
unpow2N/A
associate-+l+N/A
distribute-lft1-inN/A
+-commutativeN/A
unpow2N/A
associate-*r*N/A
lower-fma.f64N/A
Applied rewrites2.8%
Taylor expanded in k around 0
Applied rewrites19.5%
Taylor expanded in k around inf
Applied rewrites38.4%
(FPCore (a k m) :precision binary64 (if (<= m 0.45) (* 1.0 a) (* (* (* a k) k) 99.0)))
double code(double a, double k, double m) {
double tmp;
if (m <= 0.45) {
tmp = 1.0 * a;
} else {
tmp = ((a * k) * k) * 99.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.45d0) then
tmp = 1.0d0 * a
else
tmp = ((a * k) * k) * 99.0d0
end if
code = tmp
end function
public static double code(double a, double k, double m) {
double tmp;
if (m <= 0.45) {
tmp = 1.0 * a;
} else {
tmp = ((a * k) * k) * 99.0;
}
return tmp;
}
def code(a, k, m): tmp = 0 if m <= 0.45: tmp = 1.0 * a else: tmp = ((a * k) * k) * 99.0 return tmp
function code(a, k, m) tmp = 0.0 if (m <= 0.45) tmp = Float64(1.0 * a); else tmp = Float64(Float64(Float64(a * k) * k) * 99.0); end return tmp end
function tmp_2 = code(a, k, m) tmp = 0.0; if (m <= 0.45) tmp = 1.0 * a; else tmp = ((a * k) * k) * 99.0; end tmp_2 = tmp; end
code[a_, k_, m_] := If[LessEqual[m, 0.45], N[(1.0 * a), $MachinePrecision], N[(N[(N[(a * k), $MachinePrecision] * k), $MachinePrecision] * 99.0), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq 0.45:\\
\;\;\;\;1 \cdot a\\
\mathbf{else}:\\
\;\;\;\;\left(\left(a \cdot k\right) \cdot k\right) \cdot 99\\
\end{array}
\end{array}
if m < 0.450000000000000011Initial program 97.8%
lift-/.f64N/A
lift-*.f64N/A
associate-/l*N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f6497.9
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.8
Applied rewrites97.8%
Taylor expanded in k around 0
lower-pow.f6476.6
Applied rewrites76.6%
Taylor expanded in m around 0
Applied rewrites28.8%
if 0.450000000000000011 < m Initial program 81.1%
Taylor expanded in m around 0
lower-/.f64N/A
associate-+r+N/A
+-commutativeN/A
+-commutativeN/A
associate-+l+N/A
+-commutativeN/A
associate-+l+N/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
associate-*r*N/A
unpow2N/A
associate-+l+N/A
distribute-lft1-inN/A
+-commutativeN/A
unpow2N/A
associate-*r*N/A
lower-fma.f64N/A
Applied rewrites2.8%
Taylor expanded in k around 0
Applied rewrites19.5%
Taylor expanded in k around inf
Applied rewrites38.4%
(FPCore (a k m) :precision binary64 (if (<= m 1.9e+64) (* 1.0 a) (* (* a k) -10.0)))
double code(double a, double k, double m) {
double tmp;
if (m <= 1.9e+64) {
tmp = 1.0 * a;
} else {
tmp = (a * k) * -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 <= 1.9d+64) then
tmp = 1.0d0 * a
else
tmp = (a * k) * (-10.0d0)
end if
code = tmp
end function
public static double code(double a, double k, double m) {
double tmp;
if (m <= 1.9e+64) {
tmp = 1.0 * a;
} else {
tmp = (a * k) * -10.0;
}
return tmp;
}
def code(a, k, m): tmp = 0 if m <= 1.9e+64: tmp = 1.0 * a else: tmp = (a * k) * -10.0 return tmp
function code(a, k, m) tmp = 0.0 if (m <= 1.9e+64) tmp = Float64(1.0 * a); else tmp = Float64(Float64(a * k) * -10.0); end return tmp end
function tmp_2 = code(a, k, m) tmp = 0.0; if (m <= 1.9e+64) tmp = 1.0 * a; else tmp = (a * k) * -10.0; end tmp_2 = tmp; end
code[a_, k_, m_] := If[LessEqual[m, 1.9e+64], N[(1.0 * a), $MachinePrecision], N[(N[(a * k), $MachinePrecision] * -10.0), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq 1.9 \cdot 10^{+64}:\\
\;\;\;\;1 \cdot a\\
\mathbf{else}:\\
\;\;\;\;\left(a \cdot k\right) \cdot -10\\
\end{array}
\end{array}
if m < 1.9000000000000001e64Initial program 95.6%
lift-/.f64N/A
lift-*.f64N/A
associate-/l*N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f6495.6
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.6
Applied rewrites95.6%
Taylor expanded in k around 0
lower-pow.f6478.8
Applied rewrites78.8%
Taylor expanded in m around 0
Applied rewrites26.4%
if 1.9000000000000001e64 < m Initial program 83.6%
Taylor expanded in m around 0
lower-/.f64N/A
associate-+r+N/A
+-commutativeN/A
+-commutativeN/A
associate-+l+N/A
+-commutativeN/A
associate-+l+N/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
associate-*r*N/A
unpow2N/A
associate-+l+N/A
distribute-lft1-inN/A
+-commutativeN/A
unpow2N/A
associate-*r*N/A
lower-fma.f64N/A
Applied rewrites2.7%
Taylor expanded in k around 0
Applied rewrites6.3%
Taylor expanded in k around inf
Applied rewrites17.1%
(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 93.0%
lift-/.f64N/A
lift-*.f64N/A
associate-/l*N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f6493.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-+.f6493.0
Applied rewrites93.0%
Taylor expanded in k around 0
lower-pow.f6483.4
Applied rewrites83.4%
Taylor expanded in m around 0
Applied rewrites21.4%
herbie shell --seed 2024307
(FPCore (a k m)
:name "Falkner and Boettcher, Appendix A"
:precision binary64
(/ (* a (pow k m)) (+ (+ 1.0 (* 10.0 k)) (* k k))))