
(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 17 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 (* (pow k m) a)))
(if (<= (/ t_0 (- (* k k) (- -1.0 (* 10.0 k)))) 2e+249)
(* (/ (pow k m) (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 ((t_0 / ((k * k) - (-1.0 - (10.0 * k)))) <= 2e+249) {
tmp = (pow(k, m) / 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 (Float64(t_0 / Float64(Float64(k * k) - Float64(-1.0 - Float64(10.0 * k)))) <= 2e+249) tmp = Float64(Float64((k ^ m) / 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[N[(t$95$0 / N[(N[(k * k), $MachinePrecision] - N[(-1.0 - N[(10.0 * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2e+249], N[(N[(N[Power[k, m], $MachinePrecision] / 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}\;\frac{t\_0}{k \cdot k - \left(-1 - 10 \cdot k\right)} \leq 2 \cdot 10^{+249}:\\
\;\;\;\;\frac{{k}^{m}}{\mathsf{fma}\left(10 + k, k, 1\right)} \cdot a\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\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))) < 1.9999999999999998e249Initial program 97.1%
lift-/.f64N/A
lift-*.f64N/A
associate-/l*N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f6497.2
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.2
Applied rewrites97.2%
if 1.9999999999999998e249 < (/.f64 (*.f64 a (pow.f64 k m)) (+.f64 (+.f64 #s(literal 1 binary64) (*.f64 #s(literal 10 binary64) k)) (*.f64 k k))) Initial program 64.9%
lift-/.f64N/A
lift-*.f64N/A
associate-/l*N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f6464.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-+.f6464.9
Applied rewrites64.9%
Taylor expanded in k around 0
lower-pow.f64100.0
Applied rewrites100.0%
Final simplification97.8%
(FPCore (a k m)
:precision binary64
(let* ((t_0 (/ (* (pow k m) a) (- (* k k) (- -1.0 (* 10.0 k))))))
(if (<= t_0 1e-299)
(/
1.0
(/
(fma
(* (+ -100.0 (* k k)) k)
(/ (- (/ (+ (/ 100.0 k) 10.0) k) -1.0) k)
1.0)
a))
(if (<= t_0 2e+277)
(* (/ 1.0 (fma (+ 10.0 k) k 1.0)) a)
(if (<= t_0 INFINITY)
(/ (- a (* (- (/ -99.0 k) -10.0) (/ a k))) (* k k))
(fma (* (fma (- k) -99.0 -10.0) a) k a))))))
double code(double a, double k, double m) {
double t_0 = (pow(k, m) * a) / ((k * k) - (-1.0 - (10.0 * k)));
double tmp;
if (t_0 <= 1e-299) {
tmp = 1.0 / (fma(((-100.0 + (k * k)) * k), (((((100.0 / k) + 10.0) / k) - -1.0) / k), 1.0) / a);
} else if (t_0 <= 2e+277) {
tmp = (1.0 / fma((10.0 + k), k, 1.0)) * a;
} else if (t_0 <= ((double) INFINITY)) {
tmp = (a - (((-99.0 / k) - -10.0) * (a / k))) / (k * k);
} else {
tmp = fma((fma(-k, -99.0, -10.0) * a), k, a);
}
return tmp;
}
function code(a, k, m) t_0 = Float64(Float64((k ^ m) * a) / Float64(Float64(k * k) - Float64(-1.0 - Float64(10.0 * k)))) tmp = 0.0 if (t_0 <= 1e-299) tmp = Float64(1.0 / Float64(fma(Float64(Float64(-100.0 + Float64(k * k)) * k), Float64(Float64(Float64(Float64(Float64(100.0 / k) + 10.0) / k) - -1.0) / k), 1.0) / a)); elseif (t_0 <= 2e+277) tmp = Float64(Float64(1.0 / fma(Float64(10.0 + k), k, 1.0)) * a); elseif (t_0 <= Inf) tmp = Float64(Float64(a - Float64(Float64(Float64(-99.0 / k) - -10.0) * Float64(a / k))) / Float64(k * k)); else tmp = fma(Float64(fma(Float64(-k), -99.0, -10.0) * a), k, a); end return tmp end
code[a_, k_, m_] := Block[{t$95$0 = N[(N[(N[Power[k, m], $MachinePrecision] * a), $MachinePrecision] / N[(N[(k * k), $MachinePrecision] - N[(-1.0 - N[(10.0 * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, 1e-299], N[(1.0 / N[(N[(N[(N[(-100.0 + N[(k * k), $MachinePrecision]), $MachinePrecision] * k), $MachinePrecision] * N[(N[(N[(N[(N[(100.0 / k), $MachinePrecision] + 10.0), $MachinePrecision] / k), $MachinePrecision] - -1.0), $MachinePrecision] / k), $MachinePrecision] + 1.0), $MachinePrecision] / a), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, 2e+277], N[(N[(1.0 / N[(N[(10.0 + k), $MachinePrecision] * k + 1.0), $MachinePrecision]), $MachinePrecision] * a), $MachinePrecision], If[LessEqual[t$95$0, Infinity], N[(N[(a - N[(N[(N[(-99.0 / k), $MachinePrecision] - -10.0), $MachinePrecision] * N[(a / k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(k * k), $MachinePrecision]), $MachinePrecision], N[(N[(N[((-k) * -99.0 + -10.0), $MachinePrecision] * a), $MachinePrecision] * k + a), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{{k}^{m} \cdot a}{k \cdot k - \left(-1 - 10 \cdot k\right)}\\
\mathbf{if}\;t\_0 \leq 10^{-299}:\\
\;\;\;\;\frac{1}{\frac{\mathsf{fma}\left(\left(-100 + k \cdot k\right) \cdot k, \frac{\frac{\frac{100}{k} + 10}{k} - -1}{k}, 1\right)}{a}}\\
\mathbf{elif}\;t\_0 \leq 2 \cdot 10^{+277}:\\
\;\;\;\;\frac{1}{\mathsf{fma}\left(10 + k, k, 1\right)} \cdot a\\
\mathbf{elif}\;t\_0 \leq \infty:\\
\;\;\;\;\frac{a - \left(\frac{-99}{k} - -10\right) \cdot \frac{a}{k}}{k \cdot k}\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(-k, -99, -10\right) \cdot a, k, 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))) < 9.99999999999999992e-300Initial program 96.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 rewrites41.1%
Applied rewrites41.1%
Taylor expanded in k around -inf
Applied rewrites53.6%
Applied rewrites56.8%
if 9.99999999999999992e-300 < (/.f64 (*.f64 a (pow.f64 k m)) (+.f64 (+.f64 #s(literal 1 binary64) (*.f64 #s(literal 10 binary64) k)) (*.f64 k k))) < 2.00000000000000001e277Initial program 99.9%
lift-/.f64N/A
lift-*.f64N/A
associate-/l*N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f6499.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-+.f6499.9
Applied rewrites99.9%
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-+.f6496.0
Applied rewrites96.0%
if 2.00000000000000001e277 < (/.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 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 rewrites3.8%
Taylor expanded in k around inf
Applied rewrites61.2%
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 rewrites71.6%
Final simplification63.8%
(FPCore (a k m)
:precision binary64
(let* ((t_0 (/ (* (pow k m) a) (- (* k k) (- -1.0 (* 10.0 k))))))
(if (<= t_0 1e-299)
(/
a
(fma
(* (+ -100.0 (* k k)) k)
(/ (- (/ (+ (/ 100.0 k) 10.0) k) -1.0) k)
1.0))
(if (<= t_0 2e+277)
(* (/ 1.0 (fma (+ 10.0 k) k 1.0)) a)
(if (<= t_0 INFINITY)
(/ (- a (* (- (/ -99.0 k) -10.0) (/ a k))) (* k k))
(fma (* (fma (- k) -99.0 -10.0) a) k a))))))
double code(double a, double k, double m) {
double t_0 = (pow(k, m) * a) / ((k * k) - (-1.0 - (10.0 * k)));
double tmp;
if (t_0 <= 1e-299) {
tmp = a / fma(((-100.0 + (k * k)) * k), (((((100.0 / k) + 10.0) / k) - -1.0) / k), 1.0);
} else if (t_0 <= 2e+277) {
tmp = (1.0 / fma((10.0 + k), k, 1.0)) * a;
} else if (t_0 <= ((double) INFINITY)) {
tmp = (a - (((-99.0 / k) - -10.0) * (a / k))) / (k * k);
} else {
tmp = fma((fma(-k, -99.0, -10.0) * a), k, a);
}
return tmp;
}
function code(a, k, m) t_0 = Float64(Float64((k ^ m) * a) / Float64(Float64(k * k) - Float64(-1.0 - Float64(10.0 * k)))) tmp = 0.0 if (t_0 <= 1e-299) tmp = Float64(a / fma(Float64(Float64(-100.0 + Float64(k * k)) * k), Float64(Float64(Float64(Float64(Float64(100.0 / k) + 10.0) / k) - -1.0) / k), 1.0)); elseif (t_0 <= 2e+277) tmp = Float64(Float64(1.0 / fma(Float64(10.0 + k), k, 1.0)) * a); elseif (t_0 <= Inf) tmp = Float64(Float64(a - Float64(Float64(Float64(-99.0 / k) - -10.0) * Float64(a / k))) / Float64(k * k)); else tmp = fma(Float64(fma(Float64(-k), -99.0, -10.0) * a), k, a); end return tmp end
code[a_, k_, m_] := Block[{t$95$0 = N[(N[(N[Power[k, m], $MachinePrecision] * a), $MachinePrecision] / N[(N[(k * k), $MachinePrecision] - N[(-1.0 - N[(10.0 * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, 1e-299], N[(a / N[(N[(N[(-100.0 + N[(k * k), $MachinePrecision]), $MachinePrecision] * k), $MachinePrecision] * N[(N[(N[(N[(N[(100.0 / k), $MachinePrecision] + 10.0), $MachinePrecision] / k), $MachinePrecision] - -1.0), $MachinePrecision] / k), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, 2e+277], N[(N[(1.0 / N[(N[(10.0 + k), $MachinePrecision] * k + 1.0), $MachinePrecision]), $MachinePrecision] * a), $MachinePrecision], If[LessEqual[t$95$0, Infinity], N[(N[(a - N[(N[(N[(-99.0 / k), $MachinePrecision] - -10.0), $MachinePrecision] * N[(a / k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(k * k), $MachinePrecision]), $MachinePrecision], N[(N[(N[((-k) * -99.0 + -10.0), $MachinePrecision] * a), $MachinePrecision] * k + a), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{{k}^{m} \cdot a}{k \cdot k - \left(-1 - 10 \cdot k\right)}\\
\mathbf{if}\;t\_0 \leq 10^{-299}:\\
\;\;\;\;\frac{a}{\mathsf{fma}\left(\left(-100 + k \cdot k\right) \cdot k, \frac{\frac{\frac{100}{k} + 10}{k} - -1}{k}, 1\right)}\\
\mathbf{elif}\;t\_0 \leq 2 \cdot 10^{+277}:\\
\;\;\;\;\frac{1}{\mathsf{fma}\left(10 + k, k, 1\right)} \cdot a\\
\mathbf{elif}\;t\_0 \leq \infty:\\
\;\;\;\;\frac{a - \left(\frac{-99}{k} - -10\right) \cdot \frac{a}{k}}{k \cdot k}\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(-k, -99, -10\right) \cdot a, k, 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))) < 9.99999999999999992e-300Initial program 96.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 rewrites41.1%
Applied rewrites41.1%
Taylor expanded in k around -inf
Applied rewrites53.6%
Applied rewrites56.8%
if 9.99999999999999992e-300 < (/.f64 (*.f64 a (pow.f64 k m)) (+.f64 (+.f64 #s(literal 1 binary64) (*.f64 #s(literal 10 binary64) k)) (*.f64 k k))) < 2.00000000000000001e277Initial program 99.9%
lift-/.f64N/A
lift-*.f64N/A
associate-/l*N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f6499.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-+.f6499.9
Applied rewrites99.9%
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-+.f6496.0
Applied rewrites96.0%
if 2.00000000000000001e277 < (/.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 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 rewrites3.8%
Taylor expanded in k around inf
Applied rewrites61.2%
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 rewrites71.6%
Final simplification63.8%
(FPCore (a k m)
:precision binary64
(let* ((t_0 (/ (* (pow k m) a) (- (* k k) (- -1.0 (* 10.0 k))))))
(if (<= t_0 1e-299)
(/ a (fma (* (/ (- (/ 10.0 k) -1.0) k) (- (* k k) 100.0)) k 1.0))
(if (<= t_0 2e+277)
(* (/ 1.0 (fma (+ 10.0 k) k 1.0)) a)
(if (<= t_0 INFINITY)
(/ (- a (* (- (/ -99.0 k) -10.0) (/ a k))) (* k k))
(fma (* (fma (- k) -99.0 -10.0) a) k a))))))
double code(double a, double k, double m) {
double t_0 = (pow(k, m) * a) / ((k * k) - (-1.0 - (10.0 * k)));
double tmp;
if (t_0 <= 1e-299) {
tmp = a / fma(((((10.0 / k) - -1.0) / k) * ((k * k) - 100.0)), k, 1.0);
} else if (t_0 <= 2e+277) {
tmp = (1.0 / fma((10.0 + k), k, 1.0)) * a;
} else if (t_0 <= ((double) INFINITY)) {
tmp = (a - (((-99.0 / k) - -10.0) * (a / k))) / (k * k);
} else {
tmp = fma((fma(-k, -99.0, -10.0) * a), k, a);
}
return tmp;
}
function code(a, k, m) t_0 = Float64(Float64((k ^ m) * a) / Float64(Float64(k * k) - Float64(-1.0 - Float64(10.0 * k)))) tmp = 0.0 if (t_0 <= 1e-299) tmp = Float64(a / fma(Float64(Float64(Float64(Float64(10.0 / k) - -1.0) / k) * Float64(Float64(k * k) - 100.0)), k, 1.0)); elseif (t_0 <= 2e+277) tmp = Float64(Float64(1.0 / fma(Float64(10.0 + k), k, 1.0)) * a); elseif (t_0 <= Inf) tmp = Float64(Float64(a - Float64(Float64(Float64(-99.0 / k) - -10.0) * Float64(a / k))) / Float64(k * k)); else tmp = fma(Float64(fma(Float64(-k), -99.0, -10.0) * a), k, a); end return tmp end
code[a_, k_, m_] := Block[{t$95$0 = N[(N[(N[Power[k, m], $MachinePrecision] * a), $MachinePrecision] / N[(N[(k * k), $MachinePrecision] - N[(-1.0 - N[(10.0 * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, 1e-299], N[(a / N[(N[(N[(N[(N[(10.0 / k), $MachinePrecision] - -1.0), $MachinePrecision] / k), $MachinePrecision] * N[(N[(k * k), $MachinePrecision] - 100.0), $MachinePrecision]), $MachinePrecision] * k + 1.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, 2e+277], N[(N[(1.0 / N[(N[(10.0 + k), $MachinePrecision] * k + 1.0), $MachinePrecision]), $MachinePrecision] * a), $MachinePrecision], If[LessEqual[t$95$0, Infinity], N[(N[(a - N[(N[(N[(-99.0 / k), $MachinePrecision] - -10.0), $MachinePrecision] * N[(a / k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(k * k), $MachinePrecision]), $MachinePrecision], N[(N[(N[((-k) * -99.0 + -10.0), $MachinePrecision] * a), $MachinePrecision] * k + a), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{{k}^{m} \cdot a}{k \cdot k - \left(-1 - 10 \cdot k\right)}\\
\mathbf{if}\;t\_0 \leq 10^{-299}:\\
\;\;\;\;\frac{a}{\mathsf{fma}\left(\frac{\frac{10}{k} - -1}{k} \cdot \left(k \cdot k - 100\right), k, 1\right)}\\
\mathbf{elif}\;t\_0 \leq 2 \cdot 10^{+277}:\\
\;\;\;\;\frac{1}{\mathsf{fma}\left(10 + k, k, 1\right)} \cdot a\\
\mathbf{elif}\;t\_0 \leq \infty:\\
\;\;\;\;\frac{a - \left(\frac{-99}{k} - -10\right) \cdot \frac{a}{k}}{k \cdot k}\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(-k, -99, -10\right) \cdot a, k, 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))) < 9.99999999999999992e-300Initial program 96.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 rewrites41.1%
Applied rewrites41.1%
Taylor expanded in k around inf
Applied rewrites47.9%
if 9.99999999999999992e-300 < (/.f64 (*.f64 a (pow.f64 k m)) (+.f64 (+.f64 #s(literal 1 binary64) (*.f64 #s(literal 10 binary64) k)) (*.f64 k k))) < 2.00000000000000001e277Initial program 99.9%
lift-/.f64N/A
lift-*.f64N/A
associate-/l*N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f6499.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-+.f6499.9
Applied rewrites99.9%
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-+.f6496.0
Applied rewrites96.0%
if 2.00000000000000001e277 < (/.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 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 rewrites3.8%
Taylor expanded in k around inf
Applied rewrites61.2%
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 rewrites71.6%
Final simplification58.0%
(FPCore (a k m)
:precision binary64
(let* ((t_0 (/ (* (pow k m) a) (- (* k k) (- -1.0 (* 10.0 k))))))
(if (<= t_0 0.0)
(/
a
(fma
(* (fma (fma (fma -0.0001 k -0.001) k -0.01) k -0.1) (- (* k k) 100.0))
k
1.0))
(if (<= t_0 2e+277)
(/ a (fma (+ 10.0 k) k 1.0))
(if (<= t_0 INFINITY)
(/ (- a (* (- (/ -99.0 k) -10.0) (/ a k))) (* k k))
(fma (* (fma (- k) -99.0 -10.0) a) k a))))))
double code(double a, double k, double m) {
double t_0 = (pow(k, m) * a) / ((k * k) - (-1.0 - (10.0 * k)));
double tmp;
if (t_0 <= 0.0) {
tmp = a / fma((fma(fma(fma(-0.0001, k, -0.001), k, -0.01), k, -0.1) * ((k * k) - 100.0)), k, 1.0);
} else if (t_0 <= 2e+277) {
tmp = a / fma((10.0 + k), k, 1.0);
} else if (t_0 <= ((double) INFINITY)) {
tmp = (a - (((-99.0 / k) - -10.0) * (a / k))) / (k * k);
} else {
tmp = fma((fma(-k, -99.0, -10.0) * a), k, a);
}
return tmp;
}
function code(a, k, m) t_0 = Float64(Float64((k ^ m) * a) / Float64(Float64(k * k) - Float64(-1.0 - Float64(10.0 * k)))) tmp = 0.0 if (t_0 <= 0.0) tmp = Float64(a / fma(Float64(fma(fma(fma(-0.0001, k, -0.001), k, -0.01), k, -0.1) * Float64(Float64(k * k) - 100.0)), k, 1.0)); elseif (t_0 <= 2e+277) tmp = Float64(a / fma(Float64(10.0 + k), k, 1.0)); elseif (t_0 <= Inf) tmp = Float64(Float64(a - Float64(Float64(Float64(-99.0 / k) - -10.0) * Float64(a / k))) / Float64(k * k)); else tmp = fma(Float64(fma(Float64(-k), -99.0, -10.0) * a), k, a); end return tmp end
code[a_, k_, m_] := Block[{t$95$0 = N[(N[(N[Power[k, m], $MachinePrecision] * a), $MachinePrecision] / N[(N[(k * k), $MachinePrecision] - N[(-1.0 - N[(10.0 * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, 0.0], N[(a / N[(N[(N[(N[(N[(-0.0001 * k + -0.001), $MachinePrecision] * k + -0.01), $MachinePrecision] * k + -0.1), $MachinePrecision] * N[(N[(k * k), $MachinePrecision] - 100.0), $MachinePrecision]), $MachinePrecision] * k + 1.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, 2e+277], N[(a / N[(N[(10.0 + k), $MachinePrecision] * k + 1.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, Infinity], N[(N[(a - N[(N[(N[(-99.0 / k), $MachinePrecision] - -10.0), $MachinePrecision] * N[(a / k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(k * k), $MachinePrecision]), $MachinePrecision], N[(N[(N[((-k) * -99.0 + -10.0), $MachinePrecision] * a), $MachinePrecision] * k + a), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{{k}^{m} \cdot a}{k \cdot k - \left(-1 - 10 \cdot k\right)}\\
\mathbf{if}\;t\_0 \leq 0:\\
\;\;\;\;\frac{a}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(-0.0001, k, -0.001\right), k, -0.01\right), k, -0.1\right) \cdot \left(k \cdot k - 100\right), k, 1\right)}\\
\mathbf{elif}\;t\_0 \leq 2 \cdot 10^{+277}:\\
\;\;\;\;\frac{a}{\mathsf{fma}\left(10 + k, k, 1\right)}\\
\mathbf{elif}\;t\_0 \leq \infty:\\
\;\;\;\;\frac{a - \left(\frac{-99}{k} - -10\right) \cdot \frac{a}{k}}{k \cdot k}\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(-k, -99, -10\right) \cdot a, k, 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))) < 0.0Initial program 96.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 rewrites41.1%
Applied rewrites41.1%
Taylor expanded in k around 0
Applied rewrites45.7%
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))) < 2.00000000000000001e277Initial program 99.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 rewrites96.0%
if 2.00000000000000001e277 < (/.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 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 rewrites3.8%
Taylor expanded in k around inf
Applied rewrites61.2%
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 rewrites71.6%
Final simplification56.5%
(FPCore (a k m)
:precision binary64
(let* ((t_0 (/ (* (pow k m) a) (- (* k k) (- -1.0 (* 10.0 k))))))
(if (<= t_0 0.0)
(/ a (fma (* (fma (fma -0.001 k -0.01) k -0.1) (fma k k -100.0)) k 1.0))
(if (<= t_0 2e+277)
(/ a (fma (+ 10.0 k) k 1.0))
(if (<= t_0 INFINITY)
(/ (- a (* (- (/ -99.0 k) -10.0) (/ a k))) (* k k))
(fma (* (fma (- k) -99.0 -10.0) a) k a))))))
double code(double a, double k, double m) {
double t_0 = (pow(k, m) * a) / ((k * k) - (-1.0 - (10.0 * k)));
double tmp;
if (t_0 <= 0.0) {
tmp = a / fma((fma(fma(-0.001, k, -0.01), k, -0.1) * fma(k, k, -100.0)), k, 1.0);
} else if (t_0 <= 2e+277) {
tmp = a / fma((10.0 + k), k, 1.0);
} else if (t_0 <= ((double) INFINITY)) {
tmp = (a - (((-99.0 / k) - -10.0) * (a / k))) / (k * k);
} else {
tmp = fma((fma(-k, -99.0, -10.0) * a), k, a);
}
return tmp;
}
function code(a, k, m) t_0 = Float64(Float64((k ^ m) * a) / Float64(Float64(k * k) - Float64(-1.0 - Float64(10.0 * k)))) tmp = 0.0 if (t_0 <= 0.0) tmp = Float64(a / fma(Float64(fma(fma(-0.001, k, -0.01), k, -0.1) * fma(k, k, -100.0)), k, 1.0)); elseif (t_0 <= 2e+277) tmp = Float64(a / fma(Float64(10.0 + k), k, 1.0)); elseif (t_0 <= Inf) tmp = Float64(Float64(a - Float64(Float64(Float64(-99.0 / k) - -10.0) * Float64(a / k))) / Float64(k * k)); else tmp = fma(Float64(fma(Float64(-k), -99.0, -10.0) * a), k, a); end return tmp end
code[a_, k_, m_] := Block[{t$95$0 = N[(N[(N[Power[k, m], $MachinePrecision] * a), $MachinePrecision] / N[(N[(k * k), $MachinePrecision] - N[(-1.0 - N[(10.0 * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, 0.0], N[(a / N[(N[(N[(N[(-0.001 * k + -0.01), $MachinePrecision] * k + -0.1), $MachinePrecision] * N[(k * k + -100.0), $MachinePrecision]), $MachinePrecision] * k + 1.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, 2e+277], N[(a / N[(N[(10.0 + k), $MachinePrecision] * k + 1.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, Infinity], N[(N[(a - N[(N[(N[(-99.0 / k), $MachinePrecision] - -10.0), $MachinePrecision] * N[(a / k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(k * k), $MachinePrecision]), $MachinePrecision], N[(N[(N[((-k) * -99.0 + -10.0), $MachinePrecision] * a), $MachinePrecision] * k + a), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{{k}^{m} \cdot a}{k \cdot k - \left(-1 - 10 \cdot k\right)}\\
\mathbf{if}\;t\_0 \leq 0:\\
\;\;\;\;\frac{a}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(-0.001, k, -0.01\right), k, -0.1\right) \cdot \mathsf{fma}\left(k, k, -100\right), k, 1\right)}\\
\mathbf{elif}\;t\_0 \leq 2 \cdot 10^{+277}:\\
\;\;\;\;\frac{a}{\mathsf{fma}\left(10 + k, k, 1\right)}\\
\mathbf{elif}\;t\_0 \leq \infty:\\
\;\;\;\;\frac{a - \left(\frac{-99}{k} - -10\right) \cdot \frac{a}{k}}{k \cdot k}\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(-k, -99, -10\right) \cdot a, k, 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))) < 0.0Initial program 96.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 rewrites41.1%
Applied rewrites41.1%
Taylor expanded in k around 0
Applied rewrites45.2%
Taylor expanded in k around 0
Applied rewrites45.2%
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))) < 2.00000000000000001e277Initial program 99.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 rewrites96.0%
if 2.00000000000000001e277 < (/.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 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 rewrites3.8%
Taylor expanded in k around inf
Applied rewrites61.2%
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 rewrites71.6%
Final simplification56.2%
(FPCore (a k m)
:precision binary64
(let* ((t_0 (/ (* (pow k m) a) (- (* k k) (- -1.0 (* 10.0 k))))))
(if (<= t_0 0.0)
(/ a (fma (* (fma (fma -0.001 k -0.01) k -0.1) (fma k k -100.0)) k 1.0))
(if (<= t_0 2e+277)
(/ a (fma (+ 10.0 k) k 1.0))
(if (<= t_0 INFINITY)
(/ 1.0 (/ (* k k) a))
(fma (* (fma (- k) -99.0 -10.0) a) k a))))))
double code(double a, double k, double m) {
double t_0 = (pow(k, m) * a) / ((k * k) - (-1.0 - (10.0 * k)));
double tmp;
if (t_0 <= 0.0) {
tmp = a / fma((fma(fma(-0.001, k, -0.01), k, -0.1) * fma(k, k, -100.0)), k, 1.0);
} else if (t_0 <= 2e+277) {
tmp = a / fma((10.0 + k), k, 1.0);
} else if (t_0 <= ((double) INFINITY)) {
tmp = 1.0 / ((k * k) / a);
} else {
tmp = fma((fma(-k, -99.0, -10.0) * a), k, a);
}
return tmp;
}
function code(a, k, m) t_0 = Float64(Float64((k ^ m) * a) / Float64(Float64(k * k) - Float64(-1.0 - Float64(10.0 * k)))) tmp = 0.0 if (t_0 <= 0.0) tmp = Float64(a / fma(Float64(fma(fma(-0.001, k, -0.01), k, -0.1) * fma(k, k, -100.0)), k, 1.0)); elseif (t_0 <= 2e+277) tmp = Float64(a / fma(Float64(10.0 + k), k, 1.0)); elseif (t_0 <= Inf) tmp = Float64(1.0 / Float64(Float64(k * k) / a)); else tmp = fma(Float64(fma(Float64(-k), -99.0, -10.0) * a), k, a); end return tmp end
code[a_, k_, m_] := Block[{t$95$0 = N[(N[(N[Power[k, m], $MachinePrecision] * a), $MachinePrecision] / N[(N[(k * k), $MachinePrecision] - N[(-1.0 - N[(10.0 * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, 0.0], N[(a / N[(N[(N[(N[(-0.001 * k + -0.01), $MachinePrecision] * k + -0.1), $MachinePrecision] * N[(k * k + -100.0), $MachinePrecision]), $MachinePrecision] * k + 1.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, 2e+277], N[(a / N[(N[(10.0 + k), $MachinePrecision] * k + 1.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, Infinity], N[(1.0 / N[(N[(k * k), $MachinePrecision] / a), $MachinePrecision]), $MachinePrecision], N[(N[(N[((-k) * -99.0 + -10.0), $MachinePrecision] * a), $MachinePrecision] * k + a), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{{k}^{m} \cdot a}{k \cdot k - \left(-1 - 10 \cdot k\right)}\\
\mathbf{if}\;t\_0 \leq 0:\\
\;\;\;\;\frac{a}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(-0.001, k, -0.01\right), k, -0.1\right) \cdot \mathsf{fma}\left(k, k, -100\right), k, 1\right)}\\
\mathbf{elif}\;t\_0 \leq 2 \cdot 10^{+277}:\\
\;\;\;\;\frac{a}{\mathsf{fma}\left(10 + k, k, 1\right)}\\
\mathbf{elif}\;t\_0 \leq \infty:\\
\;\;\;\;\frac{1}{\frac{k \cdot k}{a}}\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(-k, -99, -10\right) \cdot a, k, 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))) < 0.0Initial program 96.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 rewrites41.1%
Applied rewrites41.1%
Taylor expanded in k around 0
Applied rewrites45.2%
Taylor expanded in k around 0
Applied rewrites45.2%
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))) < 2.00000000000000001e277Initial program 99.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 rewrites96.0%
if 2.00000000000000001e277 < (/.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 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 rewrites3.8%
Applied rewrites3.8%
Taylor expanded in k around inf
Applied rewrites50.2%
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 rewrites71.6%
Final simplification54.7%
(FPCore (a k m)
:precision binary64
(let* ((t_0 (/ (* (pow k m) a) (- (* k k) (- -1.0 (* 10.0 k))))))
(if (<= t_0 0.0)
(/ a (fma (* (fma -0.01 k -0.1) (- (* k k) 100.0)) k 1.0))
(if (<= t_0 2e+277)
(/ a (fma (+ 10.0 k) k 1.0))
(if (<= t_0 INFINITY)
(/ 1.0 (/ (* k k) a))
(fma (* (fma (- k) -99.0 -10.0) a) k a))))))
double code(double a, double k, double m) {
double t_0 = (pow(k, m) * a) / ((k * k) - (-1.0 - (10.0 * k)));
double tmp;
if (t_0 <= 0.0) {
tmp = a / fma((fma(-0.01, k, -0.1) * ((k * k) - 100.0)), k, 1.0);
} else if (t_0 <= 2e+277) {
tmp = a / fma((10.0 + k), k, 1.0);
} else if (t_0 <= ((double) INFINITY)) {
tmp = 1.0 / ((k * k) / a);
} else {
tmp = fma((fma(-k, -99.0, -10.0) * a), k, a);
}
return tmp;
}
function code(a, k, m) t_0 = Float64(Float64((k ^ m) * a) / Float64(Float64(k * k) - Float64(-1.0 - Float64(10.0 * k)))) tmp = 0.0 if (t_0 <= 0.0) tmp = Float64(a / fma(Float64(fma(-0.01, k, -0.1) * Float64(Float64(k * k) - 100.0)), k, 1.0)); elseif (t_0 <= 2e+277) tmp = Float64(a / fma(Float64(10.0 + k), k, 1.0)); elseif (t_0 <= Inf) tmp = Float64(1.0 / Float64(Float64(k * k) / a)); else tmp = fma(Float64(fma(Float64(-k), -99.0, -10.0) * a), k, a); end return tmp end
code[a_, k_, m_] := Block[{t$95$0 = N[(N[(N[Power[k, m], $MachinePrecision] * a), $MachinePrecision] / N[(N[(k * k), $MachinePrecision] - N[(-1.0 - N[(10.0 * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, 0.0], N[(a / N[(N[(N[(-0.01 * k + -0.1), $MachinePrecision] * N[(N[(k * k), $MachinePrecision] - 100.0), $MachinePrecision]), $MachinePrecision] * k + 1.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, 2e+277], N[(a / N[(N[(10.0 + k), $MachinePrecision] * k + 1.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, Infinity], N[(1.0 / N[(N[(k * k), $MachinePrecision] / a), $MachinePrecision]), $MachinePrecision], N[(N[(N[((-k) * -99.0 + -10.0), $MachinePrecision] * a), $MachinePrecision] * k + a), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{{k}^{m} \cdot a}{k \cdot k - \left(-1 - 10 \cdot k\right)}\\
\mathbf{if}\;t\_0 \leq 0:\\
\;\;\;\;\frac{a}{\mathsf{fma}\left(\mathsf{fma}\left(-0.01, k, -0.1\right) \cdot \left(k \cdot k - 100\right), k, 1\right)}\\
\mathbf{elif}\;t\_0 \leq 2 \cdot 10^{+277}:\\
\;\;\;\;\frac{a}{\mathsf{fma}\left(10 + k, k, 1\right)}\\
\mathbf{elif}\;t\_0 \leq \infty:\\
\;\;\;\;\frac{1}{\frac{k \cdot k}{a}}\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(-k, -99, -10\right) \cdot a, k, 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))) < 0.0Initial program 96.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 rewrites41.1%
Applied rewrites41.1%
Taylor expanded in k around 0
Applied rewrites42.3%
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))) < 2.00000000000000001e277Initial program 99.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 rewrites96.0%
if 2.00000000000000001e277 < (/.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 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 rewrites3.8%
Applied rewrites3.8%
Taylor expanded in k around inf
Applied rewrites50.2%
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 rewrites71.6%
Final simplification52.8%
(FPCore (a k m) :precision binary64 (if (<= k 3.8e-15) (* (pow k m) a) (* (pow k (+ -1.0 m)) (/ a k))))
double code(double a, double k, double m) {
double tmp;
if (k <= 3.8e-15) {
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 <= 3.8d-15) 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 <= 3.8e-15) {
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 <= 3.8e-15: 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 <= 3.8e-15) 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 <= 3.8e-15) tmp = (k ^ m) * a; else tmp = (k ^ (-1.0 + m)) * (a / k); end tmp_2 = tmp; end
code[a_, k_, m_] := If[LessEqual[k, 3.8e-15], 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 3.8 \cdot 10^{-15}:\\
\;\;\;\;{k}^{m} \cdot a\\
\mathbf{else}:\\
\;\;\;\;{k}^{\left(-1 + m\right)} \cdot \frac{a}{k}\\
\end{array}
\end{array}
if k < 3.8000000000000002e-15Initial program 96.9%
lift-/.f64N/A
lift-*.f64N/A
associate-/l*N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f6496.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-+.f6496.9
Applied rewrites96.9%
Taylor expanded in k around 0
lower-pow.f64100.0
Applied rewrites100.0%
if 3.8000000000000002e-15 < k Initial program 78.5%
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.f6487.4
Applied rewrites87.4%
Applied rewrites91.2%
Final simplification96.7%
(FPCore (a k m)
:precision binary64
(let* ((t_0 (* (pow k m) a)))
(if (<= m -1.25e-8)
t_0
(if (<= m 3.5e-6) (/ a (fma (+ 10.0 k) k 1.0)) t_0))))
double code(double a, double k, double m) {
double t_0 = pow(k, m) * a;
double tmp;
if (m <= -1.25e-8) {
tmp = t_0;
} else if (m <= 3.5e-6) {
tmp = a / fma((10.0 + k), k, 1.0);
} else {
tmp = t_0;
}
return tmp;
}
function code(a, k, m) t_0 = Float64((k ^ m) * a) tmp = 0.0 if (m <= -1.25e-8) tmp = t_0; elseif (m <= 3.5e-6) tmp = Float64(a / fma(Float64(10.0 + k), k, 1.0)); 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, -1.25e-8], t$95$0, If[LessEqual[m, 3.5e-6], N[(a / N[(N[(10.0 + k), $MachinePrecision] * k + 1.0), $MachinePrecision]), $MachinePrecision], t$95$0]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := {k}^{m} \cdot a\\
\mathbf{if}\;m \leq -1.25 \cdot 10^{-8}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;m \leq 3.5 \cdot 10^{-6}:\\
\;\;\;\;\frac{a}{\mathsf{fma}\left(10 + k, k, 1\right)}\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if m < -1.2499999999999999e-8 or 3.49999999999999995e-6 < m Initial program 88.7%
lift-/.f64N/A
lift-*.f64N/A
associate-/l*N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f6488.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-+.f6488.7
Applied rewrites88.7%
Taylor expanded in k around 0
lower-pow.f64100.0
Applied rewrites100.0%
if -1.2499999999999999e-8 < m < 3.49999999999999995e-6Initial program 92.8%
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 rewrites90.5%
(FPCore (a k m)
:precision binary64
(if (<= m -2.35e-11)
(/ 1.0 (/ (* k k) a))
(if (<= m 2.4e-5)
(/ a (fma (+ 10.0 k) k 1.0))
(fma (* (fma (- k) -99.0 -10.0) a) k a))))
double code(double a, double k, double m) {
double tmp;
if (m <= -2.35e-11) {
tmp = 1.0 / ((k * k) / a);
} else if (m <= 2.4e-5) {
tmp = a / fma((10.0 + k), k, 1.0);
} else {
tmp = fma((fma(-k, -99.0, -10.0) * a), k, a);
}
return tmp;
}
function code(a, k, m) tmp = 0.0 if (m <= -2.35e-11) tmp = Float64(1.0 / Float64(Float64(k * k) / a)); elseif (m <= 2.4e-5) tmp = Float64(a / fma(Float64(10.0 + k), k, 1.0)); else tmp = fma(Float64(fma(Float64(-k), -99.0, -10.0) * a), k, a); end return tmp end
code[a_, k_, m_] := If[LessEqual[m, -2.35e-11], N[(1.0 / N[(N[(k * k), $MachinePrecision] / a), $MachinePrecision]), $MachinePrecision], If[LessEqual[m, 2.4e-5], N[(a / N[(N[(10.0 + k), $MachinePrecision] * k + 1.0), $MachinePrecision]), $MachinePrecision], N[(N[(N[((-k) * -99.0 + -10.0), $MachinePrecision] * a), $MachinePrecision] * k + a), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq -2.35 \cdot 10^{-11}:\\
\;\;\;\;\frac{1}{\frac{k \cdot k}{a}}\\
\mathbf{elif}\;m \leq 2.4 \cdot 10^{-5}:\\
\;\;\;\;\frac{a}{\mathsf{fma}\left(10 + k, k, 1\right)}\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(-k, -99, -10\right) \cdot a, k, a\right)\\
\end{array}
\end{array}
if m < -2.34999999999999996e-11Initial 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 rewrites32.4%
Applied rewrites32.4%
Taylor expanded in k around inf
Applied rewrites61.1%
if -2.34999999999999996e-11 < m < 2.4000000000000001e-5Initial program 92.7%
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 rewrites90.9%
if 2.4000000000000001e-5 < m Initial program 77.3%
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.2%
Taylor expanded in k around 0
Applied rewrites25.0%
Final simplification57.8%
(FPCore (a k m)
:precision binary64
(if (<= m -0.3)
(/ a (* k k))
(if (<= m 2.4e-5)
(/ a (fma (+ 10.0 k) k 1.0))
(fma (* (fma (- k) -99.0 -10.0) a) k a))))
double code(double a, double k, double m) {
double tmp;
if (m <= -0.3) {
tmp = a / (k * k);
} else if (m <= 2.4e-5) {
tmp = a / fma((10.0 + k), k, 1.0);
} else {
tmp = fma((fma(-k, -99.0, -10.0) * a), k, a);
}
return tmp;
}
function code(a, k, m) tmp = 0.0 if (m <= -0.3) tmp = Float64(a / Float64(k * k)); elseif (m <= 2.4e-5) tmp = Float64(a / fma(Float64(10.0 + k), k, 1.0)); else tmp = fma(Float64(fma(Float64(-k), -99.0, -10.0) * a), k, a); end return tmp end
code[a_, k_, m_] := If[LessEqual[m, -0.3], N[(a / N[(k * k), $MachinePrecision]), $MachinePrecision], If[LessEqual[m, 2.4e-5], N[(a / N[(N[(10.0 + k), $MachinePrecision] * k + 1.0), $MachinePrecision]), $MachinePrecision], N[(N[(N[((-k) * -99.0 + -10.0), $MachinePrecision] * a), $MachinePrecision] * k + a), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq -0.3:\\
\;\;\;\;\frac{a}{k \cdot k}\\
\mathbf{elif}\;m \leq 2.4 \cdot 10^{-5}:\\
\;\;\;\;\frac{a}{\mathsf{fma}\left(10 + k, k, 1\right)}\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(-k, -99, -10\right) \cdot a, k, a\right)\\
\end{array}
\end{array}
if m < -0.299999999999999989Initial 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 rewrites32.1%
Taylor expanded in k around 0
Applied rewrites3.2%
Taylor expanded in k around inf
Applied rewrites61.1%
if -0.299999999999999989 < m < 2.4000000000000001e-5Initial program 92.8%
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 rewrites90.5%
if 2.4000000000000001e-5 < m Initial program 77.3%
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.2%
Taylor expanded in k around 0
Applied rewrites25.0%
Final simplification57.8%
(FPCore (a k m)
:precision binary64
(if (<= m -2e-11)
(/ a (* k k))
(if (<= m 2.4e-5)
(/ a (fma 10.0 k 1.0))
(fma (* (fma (- k) -99.0 -10.0) a) k a))))
double code(double a, double k, double m) {
double tmp;
if (m <= -2e-11) {
tmp = a / (k * k);
} else if (m <= 2.4e-5) {
tmp = a / fma(10.0, k, 1.0);
} else {
tmp = fma((fma(-k, -99.0, -10.0) * a), k, a);
}
return tmp;
}
function code(a, k, m) tmp = 0.0 if (m <= -2e-11) tmp = Float64(a / Float64(k * k)); elseif (m <= 2.4e-5) tmp = Float64(a / fma(10.0, k, 1.0)); else tmp = fma(Float64(fma(Float64(-k), -99.0, -10.0) * a), k, a); end return tmp end
code[a_, k_, m_] := If[LessEqual[m, -2e-11], N[(a / N[(k * k), $MachinePrecision]), $MachinePrecision], If[LessEqual[m, 2.4e-5], N[(a / N[(10.0 * k + 1.0), $MachinePrecision]), $MachinePrecision], N[(N[(N[((-k) * -99.0 + -10.0), $MachinePrecision] * a), $MachinePrecision] * k + a), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq -2 \cdot 10^{-11}:\\
\;\;\;\;\frac{a}{k \cdot k}\\
\mathbf{elif}\;m \leq 2.4 \cdot 10^{-5}:\\
\;\;\;\;\frac{a}{\mathsf{fma}\left(10, k, 1\right)}\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(-k, -99, -10\right) \cdot a, k, a\right)\\
\end{array}
\end{array}
if m < -1.99999999999999988e-11Initial 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 rewrites32.4%
Taylor expanded in k around 0
Applied rewrites3.1%
Taylor expanded in k around inf
Applied rewrites61.1%
if -1.99999999999999988e-11 < m < 2.4000000000000001e-5Initial program 92.7%
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 rewrites90.9%
Taylor expanded in k around 0
Applied rewrites67.3%
if 2.4000000000000001e-5 < m Initial program 77.3%
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.2%
Taylor expanded in k around 0
Applied rewrites25.0%
Final simplification50.6%
(FPCore (a k m) :precision binary64 (if (<= m -2e-11) (/ a (* k k)) (if (<= m 52000000.0) (/ a (fma 10.0 k 1.0)) (* (* -10.0 a) k))))
double code(double a, double k, double m) {
double tmp;
if (m <= -2e-11) {
tmp = a / (k * k);
} else if (m <= 52000000.0) {
tmp = a / fma(10.0, k, 1.0);
} else {
tmp = (-10.0 * a) * k;
}
return tmp;
}
function code(a, k, m) tmp = 0.0 if (m <= -2e-11) tmp = Float64(a / Float64(k * k)); elseif (m <= 52000000.0) tmp = Float64(a / fma(10.0, k, 1.0)); else tmp = Float64(Float64(-10.0 * a) * k); end return tmp end
code[a_, k_, m_] := If[LessEqual[m, -2e-11], N[(a / N[(k * k), $MachinePrecision]), $MachinePrecision], If[LessEqual[m, 52000000.0], N[(a / N[(10.0 * k + 1.0), $MachinePrecision]), $MachinePrecision], N[(N[(-10.0 * a), $MachinePrecision] * k), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq -2 \cdot 10^{-11}:\\
\;\;\;\;\frac{a}{k \cdot k}\\
\mathbf{elif}\;m \leq 52000000:\\
\;\;\;\;\frac{a}{\mathsf{fma}\left(10, k, 1\right)}\\
\mathbf{else}:\\
\;\;\;\;\left(-10 \cdot a\right) \cdot k\\
\end{array}
\end{array}
if m < -1.99999999999999988e-11Initial 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 rewrites32.4%
Taylor expanded in k around 0
Applied rewrites3.1%
Taylor expanded in k around inf
Applied rewrites61.1%
if -1.99999999999999988e-11 < m < 5.2e7Initial program 91.7%
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.0%
Taylor expanded in k around 0
Applied rewrites66.0%
if 5.2e7 < m Initial program 77.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 rewrites2.9%
Taylor expanded in k around 0
Applied rewrites3.7%
Taylor expanded in k around inf
Applied rewrites16.6%
(FPCore (a k m) :precision binary64 (let* ((t_0 (/ a (* k k)))) (if (<= k 2.7e-281) t_0 (if (<= k 3.8e-15) (* (fma -10.0 k 1.0) a) t_0))))
double code(double a, double k, double m) {
double t_0 = a / (k * k);
double tmp;
if (k <= 2.7e-281) {
tmp = t_0;
} else if (k <= 3.8e-15) {
tmp = fma(-10.0, k, 1.0) * a;
} else {
tmp = t_0;
}
return tmp;
}
function code(a, k, m) t_0 = Float64(a / Float64(k * k)) tmp = 0.0 if (k <= 2.7e-281) tmp = t_0; elseif (k <= 3.8e-15) tmp = Float64(fma(-10.0, k, 1.0) * a); else tmp = t_0; end return tmp end
code[a_, k_, m_] := Block[{t$95$0 = N[(a / N[(k * k), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[k, 2.7e-281], t$95$0, If[LessEqual[k, 3.8e-15], N[(N[(-10.0 * k + 1.0), $MachinePrecision] * a), $MachinePrecision], t$95$0]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{a}{k \cdot k}\\
\mathbf{if}\;k \leq 2.7 \cdot 10^{-281}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;k \leq 3.8 \cdot 10^{-15}:\\
\;\;\;\;\mathsf{fma}\left(-10, k, 1\right) \cdot a\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if k < 2.6999999999999999e-281 or 3.8000000000000002e-15 < k Initial program 85.7%
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 rewrites34.4%
Taylor expanded in k around 0
Applied rewrites4.1%
Taylor expanded in k around inf
Applied rewrites41.0%
if 2.6999999999999999e-281 < k < 3.8000000000000002e-15Initial 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 rewrites53.6%
Taylor expanded in k around 0
Applied rewrites53.6%
Taylor expanded in k around 0
Applied rewrites53.6%
(FPCore (a k m) :precision binary64 (if (<= m 1550000.0) (* (fma -10.0 k 1.0) a) (* (* -10.0 a) k)))
double code(double a, double k, double m) {
double tmp;
if (m <= 1550000.0) {
tmp = fma(-10.0, k, 1.0) * a;
} else {
tmp = (-10.0 * a) * k;
}
return tmp;
}
function code(a, k, m) tmp = 0.0 if (m <= 1550000.0) tmp = Float64(fma(-10.0, k, 1.0) * a); else tmp = Float64(Float64(-10.0 * a) * k); end return tmp end
code[a_, k_, m_] := If[LessEqual[m, 1550000.0], N[(N[(-10.0 * k + 1.0), $MachinePrecision] * a), $MachinePrecision], N[(N[(-10.0 * a), $MachinePrecision] * k), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq 1550000:\\
\;\;\;\;\mathsf{fma}\left(-10, k, 1\right) \cdot a\\
\mathbf{else}:\\
\;\;\;\;\left(-10 \cdot a\right) \cdot k\\
\end{array}
\end{array}
if m < 1.55e6Initial program 96.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
associate-*l*N/A
unpow2N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
lower-fma.f64N/A
Applied rewrites59.0%
Taylor expanded in k around 0
Applied rewrites26.7%
Taylor expanded in k around 0
Applied rewrites26.7%
if 1.55e6 < m Initial program 77.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 rewrites2.9%
Taylor expanded in k around 0
Applied rewrites3.7%
Taylor expanded in k around inf
Applied rewrites16.6%
(FPCore (a k m) :precision binary64 (* (* -10.0 a) k))
double code(double a, double k, double m) {
return (-10.0 * a) * k;
}
real(8) function code(a, k, m)
real(8), intent (in) :: a
real(8), intent (in) :: k
real(8), intent (in) :: m
code = ((-10.0d0) * a) * k
end function
public static double code(double a, double k, double m) {
return (-10.0 * a) * k;
}
def code(a, k, m): return (-10.0 * a) * k
function code(a, k, m) return Float64(Float64(-10.0 * a) * k) end
function tmp = code(a, k, m) tmp = (-10.0 * a) * k; end
code[a_, k_, m_] := N[(N[(-10.0 * a), $MachinePrecision] * k), $MachinePrecision]
\begin{array}{l}
\\
\left(-10 \cdot a\right) \cdot k
\end{array}
Initial program 90.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 rewrites40.2%
Taylor expanded in k around 0
Applied rewrites19.0%
Taylor expanded in k around inf
Applied rewrites7.0%
herbie shell --seed 2024288
(FPCore (a k m)
:name "Falkner and Boettcher, Appendix A"
:precision binary64
(/ (* a (pow k m)) (+ (+ 1.0 (* 10.0 k)) (* k k))))