
(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 16 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 (<= (/ (* a (pow k m)) (+ (+ 1.0 (* 10.0 k)) (* k k))) 5e+274) (* a (* (pow k m) (pow (fma (+ k 10.0) k 1.0) -1.0))) (* (pow k m) a)))
double code(double a, double k, double m) {
double tmp;
if (((a * pow(k, m)) / ((1.0 + (10.0 * k)) + (k * k))) <= 5e+274) {
tmp = a * (pow(k, m) * pow(fma((k + 10.0), k, 1.0), -1.0));
} else {
tmp = pow(k, m) * 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))) <= 5e+274) tmp = Float64(a * Float64((k ^ m) * (fma(Float64(k + 10.0), k, 1.0) ^ -1.0))); else tmp = Float64((k ^ m) * 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], 5e+274], N[(a * N[(N[Power[k, m], $MachinePrecision] * N[Power[N[(N[(k + 10.0), $MachinePrecision] * k + 1.0), $MachinePrecision], -1.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[Power[k, m], $MachinePrecision] * 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 5 \cdot 10^{+274}:\\
\;\;\;\;a \cdot \left({k}^{m} \cdot {\left(\mathsf{fma}\left(k + 10, k, 1\right)\right)}^{-1}\right)\\
\mathbf{else}:\\
\;\;\;\;{k}^{m} \cdot a\\
\end{array}
\end{array}
if (/.f64 (*.f64 a (pow.f64 k m)) (+.f64 (+.f64 #s(literal 1 binary64) (*.f64 #s(literal 10 binary64) k)) (*.f64 k k))) < 4.9999999999999998e274Initial program 97.3%
lift-/.f64N/A
frac-2negN/A
div-invN/A
lift-*.f64N/A
distribute-lft-neg-inN/A
associate-*l*N/A
lower-*.f64N/A
lower-neg.f64N/A
lower-*.f64N/A
neg-mul-1N/A
associate-/r*N/A
metadata-evalN/A
lower-/.f6497.3
lift-+.f64N/A
lift-+.f64N/A
associate-+l+N/A
+-commutativeN/A
Applied rewrites97.3%
if 4.9999999999999998e274 < (/.f64 (*.f64 a (pow.f64 k m)) (+.f64 (+.f64 #s(literal 1 binary64) (*.f64 #s(literal 10 binary64) k)) (*.f64 k k))) Initial program 47.8%
lift-/.f64N/A
frac-2negN/A
div-invN/A
lift-*.f64N/A
distribute-lft-neg-inN/A
associate-*l*N/A
lower-*.f64N/A
lower-neg.f64N/A
lower-*.f64N/A
neg-mul-1N/A
associate-/r*N/A
metadata-evalN/A
lower-/.f6447.8
lift-+.f64N/A
lift-+.f64N/A
associate-+l+N/A
+-commutativeN/A
Applied rewrites50.0%
Taylor expanded in k around 0
mul-1-negN/A
lower-neg.f64N/A
lower-pow.f64100.0
Applied rewrites100.0%
Taylor expanded in k around 0
*-commutativeN/A
lower-*.f64N/A
lower-pow.f64100.0
Applied rewrites100.0%
Final simplification97.8%
(FPCore (a k m) :precision binary64 (if (<= (/ (* a (pow k m)) (+ (+ 1.0 (* 10.0 k)) (* k k))) 5e+274) (* (/ (pow k m) (fma (+ k 10.0) k 1.0)) a) (* (pow k m) a)))
double code(double a, double k, double m) {
double tmp;
if (((a * pow(k, m)) / ((1.0 + (10.0 * k)) + (k * k))) <= 5e+274) {
tmp = (pow(k, m) / fma((k + 10.0), k, 1.0)) * a;
} else {
tmp = pow(k, m) * 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))) <= 5e+274) tmp = Float64(Float64((k ^ m) / fma(Float64(k + 10.0), k, 1.0)) * a); else tmp = Float64((k ^ m) * 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], 5e+274], N[(N[(N[Power[k, m], $MachinePrecision] / N[(N[(k + 10.0), $MachinePrecision] * k + 1.0), $MachinePrecision]), $MachinePrecision] * a), $MachinePrecision], N[(N[Power[k, m], $MachinePrecision] * 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 5 \cdot 10^{+274}:\\
\;\;\;\;\frac{{k}^{m}}{\mathsf{fma}\left(k + 10, k, 1\right)} \cdot a\\
\mathbf{else}:\\
\;\;\;\;{k}^{m} \cdot a\\
\end{array}
\end{array}
if (/.f64 (*.f64 a (pow.f64 k m)) (+.f64 (+.f64 #s(literal 1 binary64) (*.f64 #s(literal 10 binary64) k)) (*.f64 k k))) < 4.9999999999999998e274Initial program 97.3%
lift-/.f64N/A
lift-*.f64N/A
associate-/l*N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f6497.3
lift-+.f64N/A
lift-+.f64N/A
associate-+l+N/A
+-commutativeN/A
lift-*.f64N/A
lift-*.f64N/A
distribute-rgt-outN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-+.f6497.3
Applied rewrites97.3%
if 4.9999999999999998e274 < (/.f64 (*.f64 a (pow.f64 k m)) (+.f64 (+.f64 #s(literal 1 binary64) (*.f64 #s(literal 10 binary64) k)) (*.f64 k k))) Initial program 47.8%
lift-/.f64N/A
frac-2negN/A
div-invN/A
lift-*.f64N/A
distribute-lft-neg-inN/A
associate-*l*N/A
lower-*.f64N/A
lower-neg.f64N/A
lower-*.f64N/A
neg-mul-1N/A
associate-/r*N/A
metadata-evalN/A
lower-/.f6447.8
lift-+.f64N/A
lift-+.f64N/A
associate-+l+N/A
+-commutativeN/A
Applied rewrites50.0%
Taylor expanded in k around 0
mul-1-negN/A
lower-neg.f64N/A
lower-pow.f64100.0
Applied rewrites100.0%
Taylor expanded in k around 0
*-commutativeN/A
lower-*.f64N/A
lower-pow.f64100.0
Applied rewrites100.0%
(FPCore (a k m)
:precision binary64
(if (<= m -5.4e+182)
(/
a
(fma
(* (* k k) (- 100.0 (* k k)))
(/ (fma (fma (fma 0.0001 k 0.001) k 0.01) k 0.1) k)
1.0))
(if (<= m -6400000000.0)
(/ (- a (* (/ a k) (- (/ -99.0 k) -10.0))) (* k k))
(if (<= m 1.35)
(* (pow (fma (+ 10.0 k) k 1.0) -1.0) a)
(* (* (* k a) k) 99.0)))))
double code(double a, double k, double m) {
double tmp;
if (m <= -5.4e+182) {
tmp = a / fma(((k * k) * (100.0 - (k * k))), (fma(fma(fma(0.0001, k, 0.001), k, 0.01), k, 0.1) / k), 1.0);
} else if (m <= -6400000000.0) {
tmp = (a - ((a / k) * ((-99.0 / k) - -10.0))) / (k * k);
} else if (m <= 1.35) {
tmp = pow(fma((10.0 + k), k, 1.0), -1.0) * a;
} else {
tmp = ((k * a) * k) * 99.0;
}
return tmp;
}
function code(a, k, m) tmp = 0.0 if (m <= -5.4e+182) tmp = Float64(a / fma(Float64(Float64(k * k) * Float64(100.0 - Float64(k * k))), Float64(fma(fma(fma(0.0001, k, 0.001), k, 0.01), k, 0.1) / k), 1.0)); elseif (m <= -6400000000.0) tmp = Float64(Float64(a - Float64(Float64(a / k) * Float64(Float64(-99.0 / k) - -10.0))) / Float64(k * k)); elseif (m <= 1.35) tmp = Float64((fma(Float64(10.0 + k), k, 1.0) ^ -1.0) * a); else tmp = Float64(Float64(Float64(k * a) * k) * 99.0); end return tmp end
code[a_, k_, m_] := If[LessEqual[m, -5.4e+182], N[(a / N[(N[(N[(k * k), $MachinePrecision] * N[(100.0 - N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[(N[(N[(0.0001 * k + 0.001), $MachinePrecision] * k + 0.01), $MachinePrecision] * k + 0.1), $MachinePrecision] / k), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[m, -6400000000.0], N[(N[(a - N[(N[(a / k), $MachinePrecision] * N[(N[(-99.0 / k), $MachinePrecision] - -10.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(k * k), $MachinePrecision]), $MachinePrecision], If[LessEqual[m, 1.35], N[(N[Power[N[(N[(10.0 + k), $MachinePrecision] * k + 1.0), $MachinePrecision], -1.0], $MachinePrecision] * a), $MachinePrecision], N[(N[(N[(k * a), $MachinePrecision] * k), $MachinePrecision] * 99.0), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq -5.4 \cdot 10^{+182}:\\
\;\;\;\;\frac{a}{\mathsf{fma}\left(\left(k \cdot k\right) \cdot \left(100 - k \cdot k\right), \frac{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(0.0001, k, 0.001\right), k, 0.01\right), k, 0.1\right)}{k}, 1\right)}\\
\mathbf{elif}\;m \leq -6400000000:\\
\;\;\;\;\frac{a - \frac{a}{k} \cdot \left(\frac{-99}{k} - -10\right)}{k \cdot k}\\
\mathbf{elif}\;m \leq 1.35:\\
\;\;\;\;{\left(\mathsf{fma}\left(10 + k, k, 1\right)\right)}^{-1} \cdot a\\
\mathbf{else}:\\
\;\;\;\;\left(\left(k \cdot a\right) \cdot k\right) \cdot 99\\
\end{array}
\end{array}
if m < -5.4000000000000006e182Initial program 96.4%
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 rewrites47.7%
Applied rewrites26.5%
Taylor expanded in k around 0
Applied rewrites72.3%
if -5.4000000000000006e182 < m < -6.4e9Initial 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 rewrites38.3%
Taylor expanded in k around inf
Applied rewrites77.5%
if -6.4e9 < m < 1.3500000000000001Initial program 93.5%
lift-/.f64N/A
lift-*.f64N/A
associate-/l*N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f6493.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-+.f6493.6
Applied rewrites93.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
*-lft-identityN/A
distribute-rgt-inN/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-+.f6490.2
Applied rewrites90.2%
if 1.3500000000000001 < m Initial program 76.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 rewrites2.8%
Taylor expanded in k around 0
Applied rewrites25.7%
Taylor expanded in k around inf
Applied rewrites54.6%
Final simplification72.7%
(FPCore (a k m)
:precision binary64
(if (<= m -5.4e+182)
(/ a (fma (* (* k k) (- 100.0 (* k k))) (fma 0.001 k 0.01) 1.0))
(if (<= m -6400000000.0)
(/ (- a (* (/ a k) (- (/ -99.0 k) -10.0))) (* k k))
(if (<= m 1.35)
(* (pow (fma (+ 10.0 k) k 1.0) -1.0) a)
(* (* (* k a) k) 99.0)))))
double code(double a, double k, double m) {
double tmp;
if (m <= -5.4e+182) {
tmp = a / fma(((k * k) * (100.0 - (k * k))), fma(0.001, k, 0.01), 1.0);
} else if (m <= -6400000000.0) {
tmp = (a - ((a / k) * ((-99.0 / k) - -10.0))) / (k * k);
} else if (m <= 1.35) {
tmp = pow(fma((10.0 + k), k, 1.0), -1.0) * a;
} else {
tmp = ((k * a) * k) * 99.0;
}
return tmp;
}
function code(a, k, m) tmp = 0.0 if (m <= -5.4e+182) tmp = Float64(a / fma(Float64(Float64(k * k) * Float64(100.0 - Float64(k * k))), fma(0.001, k, 0.01), 1.0)); elseif (m <= -6400000000.0) tmp = Float64(Float64(a - Float64(Float64(a / k) * Float64(Float64(-99.0 / k) - -10.0))) / Float64(k * k)); elseif (m <= 1.35) tmp = Float64((fma(Float64(10.0 + k), k, 1.0) ^ -1.0) * a); else tmp = Float64(Float64(Float64(k * a) * k) * 99.0); end return tmp end
code[a_, k_, m_] := If[LessEqual[m, -5.4e+182], N[(a / N[(N[(N[(k * k), $MachinePrecision] * N[(100.0 - N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(0.001 * k + 0.01), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[m, -6400000000.0], N[(N[(a - N[(N[(a / k), $MachinePrecision] * N[(N[(-99.0 / k), $MachinePrecision] - -10.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(k * k), $MachinePrecision]), $MachinePrecision], If[LessEqual[m, 1.35], N[(N[Power[N[(N[(10.0 + k), $MachinePrecision] * k + 1.0), $MachinePrecision], -1.0], $MachinePrecision] * a), $MachinePrecision], N[(N[(N[(k * a), $MachinePrecision] * k), $MachinePrecision] * 99.0), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq -5.4 \cdot 10^{+182}:\\
\;\;\;\;\frac{a}{\mathsf{fma}\left(\left(k \cdot k\right) \cdot \left(100 - k \cdot k\right), \mathsf{fma}\left(0.001, k, 0.01\right), 1\right)}\\
\mathbf{elif}\;m \leq -6400000000:\\
\;\;\;\;\frac{a - \frac{a}{k} \cdot \left(\frac{-99}{k} - -10\right)}{k \cdot k}\\
\mathbf{elif}\;m \leq 1.35:\\
\;\;\;\;{\left(\mathsf{fma}\left(10 + k, k, 1\right)\right)}^{-1} \cdot a\\
\mathbf{else}:\\
\;\;\;\;\left(\left(k \cdot a\right) \cdot k\right) \cdot 99\\
\end{array}
\end{array}
if m < -5.4000000000000006e182Initial program 96.4%
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 rewrites47.7%
Applied rewrites26.5%
Taylor expanded in k around 0
Applied rewrites65.8%
Taylor expanded in k around inf
Applied rewrites65.8%
if -5.4000000000000006e182 < m < -6.4e9Initial 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 rewrites38.3%
Taylor expanded in k around inf
Applied rewrites77.5%
if -6.4e9 < m < 1.3500000000000001Initial program 93.5%
lift-/.f64N/A
lift-*.f64N/A
associate-/l*N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f6493.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-+.f6493.6
Applied rewrites93.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
*-lft-identityN/A
distribute-rgt-inN/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-+.f6490.2
Applied rewrites90.2%
if 1.3500000000000001 < m Initial program 76.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 rewrites2.8%
Taylor expanded in k around 0
Applied rewrites25.7%
Taylor expanded in k around inf
Applied rewrites54.6%
Final simplification72.0%
(FPCore (a k m)
:precision binary64
(if (<= m -5.4e+182)
(/ a (fma (* (* k k) (- 100.0 (* k k))) (fma 0.001 k 0.01) 1.0))
(if (<= m -6400000000.0)
(/ a (* k k))
(if (<= m 1.35)
(* (pow (fma (+ 10.0 k) k 1.0) -1.0) a)
(* (* (* k a) k) 99.0)))))
double code(double a, double k, double m) {
double tmp;
if (m <= -5.4e+182) {
tmp = a / fma(((k * k) * (100.0 - (k * k))), fma(0.001, k, 0.01), 1.0);
} else if (m <= -6400000000.0) {
tmp = a / (k * k);
} else if (m <= 1.35) {
tmp = pow(fma((10.0 + k), k, 1.0), -1.0) * a;
} else {
tmp = ((k * a) * k) * 99.0;
}
return tmp;
}
function code(a, k, m) tmp = 0.0 if (m <= -5.4e+182) tmp = Float64(a / fma(Float64(Float64(k * k) * Float64(100.0 - Float64(k * k))), fma(0.001, k, 0.01), 1.0)); elseif (m <= -6400000000.0) tmp = Float64(a / Float64(k * k)); elseif (m <= 1.35) tmp = Float64((fma(Float64(10.0 + k), k, 1.0) ^ -1.0) * a); else tmp = Float64(Float64(Float64(k * a) * k) * 99.0); end return tmp end
code[a_, k_, m_] := If[LessEqual[m, -5.4e+182], N[(a / N[(N[(N[(k * k), $MachinePrecision] * N[(100.0 - N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(0.001 * k + 0.01), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[m, -6400000000.0], N[(a / N[(k * k), $MachinePrecision]), $MachinePrecision], If[LessEqual[m, 1.35], N[(N[Power[N[(N[(10.0 + k), $MachinePrecision] * k + 1.0), $MachinePrecision], -1.0], $MachinePrecision] * a), $MachinePrecision], N[(N[(N[(k * a), $MachinePrecision] * k), $MachinePrecision] * 99.0), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq -5.4 \cdot 10^{+182}:\\
\;\;\;\;\frac{a}{\mathsf{fma}\left(\left(k \cdot k\right) \cdot \left(100 - k \cdot k\right), \mathsf{fma}\left(0.001, k, 0.01\right), 1\right)}\\
\mathbf{elif}\;m \leq -6400000000:\\
\;\;\;\;\frac{a}{k \cdot k}\\
\mathbf{elif}\;m \leq 1.35:\\
\;\;\;\;{\left(\mathsf{fma}\left(10 + k, k, 1\right)\right)}^{-1} \cdot a\\
\mathbf{else}:\\
\;\;\;\;\left(\left(k \cdot a\right) \cdot k\right) \cdot 99\\
\end{array}
\end{array}
if m < -5.4000000000000006e182Initial program 96.4%
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 rewrites47.7%
Applied rewrites26.5%
Taylor expanded in k around 0
Applied rewrites65.8%
Taylor expanded in k around inf
Applied rewrites65.8%
if -5.4000000000000006e182 < m < -6.4e9Initial 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 rewrites38.3%
Taylor expanded in k around inf
Applied rewrites71.2%
if -6.4e9 < m < 1.3500000000000001Initial program 93.5%
lift-/.f64N/A
lift-*.f64N/A
associate-/l*N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f6493.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-+.f6493.6
Applied rewrites93.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
*-lft-identityN/A
distribute-rgt-inN/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-+.f6490.2
Applied rewrites90.2%
if 1.3500000000000001 < m Initial program 76.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 rewrites2.8%
Taylor expanded in k around 0
Applied rewrites25.7%
Taylor expanded in k around inf
Applied rewrites54.6%
Final simplification70.9%
(FPCore (a k m)
:precision binary64
(if (<= m -5.4e+182)
(/ a (fma (* (* k k) (- 100.0 (* k k))) (* 0.001 k) 1.0))
(if (<= m -6400000000.0)
(/ a (* k k))
(if (<= m 1.35)
(* (pow (fma (+ 10.0 k) k 1.0) -1.0) a)
(* (* (* k a) k) 99.0)))))
double code(double a, double k, double m) {
double tmp;
if (m <= -5.4e+182) {
tmp = a / fma(((k * k) * (100.0 - (k * k))), (0.001 * k), 1.0);
} else if (m <= -6400000000.0) {
tmp = a / (k * k);
} else if (m <= 1.35) {
tmp = pow(fma((10.0 + k), k, 1.0), -1.0) * a;
} else {
tmp = ((k * a) * k) * 99.0;
}
return tmp;
}
function code(a, k, m) tmp = 0.0 if (m <= -5.4e+182) tmp = Float64(a / fma(Float64(Float64(k * k) * Float64(100.0 - Float64(k * k))), Float64(0.001 * k), 1.0)); elseif (m <= -6400000000.0) tmp = Float64(a / Float64(k * k)); elseif (m <= 1.35) tmp = Float64((fma(Float64(10.0 + k), k, 1.0) ^ -1.0) * a); else tmp = Float64(Float64(Float64(k * a) * k) * 99.0); end return tmp end
code[a_, k_, m_] := If[LessEqual[m, -5.4e+182], N[(a / N[(N[(N[(k * k), $MachinePrecision] * N[(100.0 - N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(0.001 * k), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[m, -6400000000.0], N[(a / N[(k * k), $MachinePrecision]), $MachinePrecision], If[LessEqual[m, 1.35], N[(N[Power[N[(N[(10.0 + k), $MachinePrecision] * k + 1.0), $MachinePrecision], -1.0], $MachinePrecision] * a), $MachinePrecision], N[(N[(N[(k * a), $MachinePrecision] * k), $MachinePrecision] * 99.0), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq -5.4 \cdot 10^{+182}:\\
\;\;\;\;\frac{a}{\mathsf{fma}\left(\left(k \cdot k\right) \cdot \left(100 - k \cdot k\right), 0.001 \cdot k, 1\right)}\\
\mathbf{elif}\;m \leq -6400000000:\\
\;\;\;\;\frac{a}{k \cdot k}\\
\mathbf{elif}\;m \leq 1.35:\\
\;\;\;\;{\left(\mathsf{fma}\left(10 + k, k, 1\right)\right)}^{-1} \cdot a\\
\mathbf{else}:\\
\;\;\;\;\left(\left(k \cdot a\right) \cdot k\right) \cdot 99\\
\end{array}
\end{array}
if m < -5.4000000000000006e182Initial program 96.4%
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 rewrites47.7%
Applied rewrites26.5%
Taylor expanded in k around 0
Applied rewrites65.8%
Taylor expanded in k around inf
Applied rewrites65.8%
if -5.4000000000000006e182 < m < -6.4e9Initial 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 rewrites38.3%
Taylor expanded in k around inf
Applied rewrites71.2%
if -6.4e9 < m < 1.3500000000000001Initial program 93.5%
lift-/.f64N/A
lift-*.f64N/A
associate-/l*N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f6493.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-+.f6493.6
Applied rewrites93.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
*-lft-identityN/A
distribute-rgt-inN/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-+.f6490.2
Applied rewrites90.2%
if 1.3500000000000001 < m Initial program 76.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 rewrites2.8%
Taylor expanded in k around 0
Applied rewrites25.7%
Taylor expanded in k around inf
Applied rewrites54.6%
Final simplification70.9%
(FPCore (a k m) :precision binary64 (if (or (<= m -9.2e-7) (not (<= m 7.5e-5))) (* (pow k m) a) (* (pow (fma (+ 10.0 k) k 1.0) -1.0) a)))
double code(double a, double k, double m) {
double tmp;
if ((m <= -9.2e-7) || !(m <= 7.5e-5)) {
tmp = pow(k, m) * a;
} else {
tmp = pow(fma((10.0 + k), k, 1.0), -1.0) * a;
}
return tmp;
}
function code(a, k, m) tmp = 0.0 if ((m <= -9.2e-7) || !(m <= 7.5e-5)) tmp = Float64((k ^ m) * a); else tmp = Float64((fma(Float64(10.0 + k), k, 1.0) ^ -1.0) * a); end return tmp end
code[a_, k_, m_] := If[Or[LessEqual[m, -9.2e-7], N[Not[LessEqual[m, 7.5e-5]], $MachinePrecision]], N[(N[Power[k, m], $MachinePrecision] * a), $MachinePrecision], N[(N[Power[N[(N[(10.0 + k), $MachinePrecision] * k + 1.0), $MachinePrecision], -1.0], $MachinePrecision] * a), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq -9.2 \cdot 10^{-7} \lor \neg \left(m \leq 7.5 \cdot 10^{-5}\right):\\
\;\;\;\;{k}^{m} \cdot a\\
\mathbf{else}:\\
\;\;\;\;{\left(\mathsf{fma}\left(10 + k, k, 1\right)\right)}^{-1} \cdot a\\
\end{array}
\end{array}
if m < -9.1999999999999998e-7 or 7.49999999999999934e-5 < m Initial program 86.0%
lift-/.f64N/A
frac-2negN/A
div-invN/A
lift-*.f64N/A
distribute-lft-neg-inN/A
associate-*l*N/A
lower-*.f64N/A
lower-neg.f64N/A
lower-*.f64N/A
neg-mul-1N/A
associate-/r*N/A
metadata-evalN/A
lower-/.f6486.0
lift-+.f64N/A
lift-+.f64N/A
associate-+l+N/A
+-commutativeN/A
Applied rewrites86.6%
Taylor expanded in k around 0
mul-1-negN/A
lower-neg.f64N/A
lower-pow.f64100.0
Applied rewrites100.0%
Taylor expanded in k around 0
*-commutativeN/A
lower-*.f64N/A
lower-pow.f64100.0
Applied rewrites100.0%
if -9.1999999999999998e-7 < m < 7.49999999999999934e-5Initial program 93.2%
lift-/.f64N/A
lift-*.f64N/A
associate-/l*N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f6493.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-+.f6493.3
Applied rewrites93.3%
Taylor expanded in m around 0
lower-/.f64N/A
+-commutativeN/A
*-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/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-+.f6492.7
Applied rewrites92.7%
Final simplification97.6%
(FPCore (a k m)
:precision binary64
(if (<= m -6400000000.0)
(/ a (* k k))
(if (<= m 1.35)
(* (pow (fma (+ 10.0 k) k 1.0) -1.0) a)
(* (* (* k a) k) 99.0))))
double code(double a, double k, double m) {
double tmp;
if (m <= -6400000000.0) {
tmp = a / (k * k);
} else if (m <= 1.35) {
tmp = pow(fma((10.0 + k), k, 1.0), -1.0) * a;
} else {
tmp = ((k * a) * k) * 99.0;
}
return tmp;
}
function code(a, k, m) tmp = 0.0 if (m <= -6400000000.0) tmp = Float64(a / Float64(k * k)); elseif (m <= 1.35) tmp = Float64((fma(Float64(10.0 + k), k, 1.0) ^ -1.0) * a); else tmp = Float64(Float64(Float64(k * a) * k) * 99.0); end return tmp end
code[a_, k_, m_] := If[LessEqual[m, -6400000000.0], N[(a / N[(k * k), $MachinePrecision]), $MachinePrecision], If[LessEqual[m, 1.35], N[(N[Power[N[(N[(10.0 + k), $MachinePrecision] * k + 1.0), $MachinePrecision], -1.0], $MachinePrecision] * a), $MachinePrecision], N[(N[(N[(k * a), $MachinePrecision] * k), $MachinePrecision] * 99.0), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq -6400000000:\\
\;\;\;\;\frac{a}{k \cdot k}\\
\mathbf{elif}\;m \leq 1.35:\\
\;\;\;\;{\left(\mathsf{fma}\left(10 + k, k, 1\right)\right)}^{-1} \cdot a\\
\mathbf{else}:\\
\;\;\;\;\left(\left(k \cdot a\right) \cdot k\right) \cdot 99\\
\end{array}
\end{array}
if m < -6.4e9Initial program 98.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 rewrites42.0%
Taylor expanded in k around inf
Applied rewrites63.4%
if -6.4e9 < m < 1.3500000000000001Initial program 93.5%
lift-/.f64N/A
lift-*.f64N/A
associate-/l*N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f6493.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-+.f6493.6
Applied rewrites93.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
*-lft-identityN/A
distribute-rgt-inN/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-+.f6490.2
Applied rewrites90.2%
if 1.3500000000000001 < m Initial program 76.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 rewrites2.8%
Taylor expanded in k around 0
Applied rewrites25.7%
Taylor expanded in k around inf
Applied rewrites54.6%
Final simplification69.3%
(FPCore (a k m) :precision binary64 (if (<= m -6400000000.0) (/ a (* k k)) (if (<= m 1.35) (/ a (fma (+ 10.0 k) k 1.0)) (* (* (* k a) k) 99.0))))
double code(double a, double k, double m) {
double tmp;
if (m <= -6400000000.0) {
tmp = a / (k * k);
} else if (m <= 1.35) {
tmp = a / fma((10.0 + k), k, 1.0);
} else {
tmp = ((k * a) * k) * 99.0;
}
return tmp;
}
function code(a, k, m) tmp = 0.0 if (m <= -6400000000.0) tmp = Float64(a / Float64(k * k)); elseif (m <= 1.35) tmp = Float64(a / fma(Float64(10.0 + k), k, 1.0)); else tmp = Float64(Float64(Float64(k * a) * k) * 99.0); end return tmp end
code[a_, k_, m_] := If[LessEqual[m, -6400000000.0], N[(a / N[(k * k), $MachinePrecision]), $MachinePrecision], If[LessEqual[m, 1.35], N[(a / N[(N[(10.0 + k), $MachinePrecision] * k + 1.0), $MachinePrecision]), $MachinePrecision], N[(N[(N[(k * a), $MachinePrecision] * k), $MachinePrecision] * 99.0), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq -6400000000:\\
\;\;\;\;\frac{a}{k \cdot k}\\
\mathbf{elif}\;m \leq 1.35:\\
\;\;\;\;\frac{a}{\mathsf{fma}\left(10 + k, k, 1\right)}\\
\mathbf{else}:\\
\;\;\;\;\left(\left(k \cdot a\right) \cdot k\right) \cdot 99\\
\end{array}
\end{array}
if m < -6.4e9Initial program 98.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 rewrites42.0%
Taylor expanded in k around inf
Applied rewrites63.4%
if -6.4e9 < m < 1.3500000000000001Initial program 93.5%
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 rewrites90.1%
if 1.3500000000000001 < m Initial program 76.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 rewrites2.8%
Taylor expanded in k around 0
Applied rewrites25.7%
Taylor expanded in k around inf
Applied rewrites54.6%
Final simplification69.3%
(FPCore (a k m) :precision binary64 (if (<= m -6.2e-16) (/ a (* k k)) (if (<= m 1.35) (/ a (fma 10.0 k 1.0)) (* (* (* k a) k) 99.0))))
double code(double a, double k, double m) {
double tmp;
if (m <= -6.2e-16) {
tmp = a / (k * k);
} else if (m <= 1.35) {
tmp = a / fma(10.0, k, 1.0);
} else {
tmp = ((k * a) * k) * 99.0;
}
return tmp;
}
function code(a, k, m) tmp = 0.0 if (m <= -6.2e-16) tmp = Float64(a / Float64(k * k)); elseif (m <= 1.35) tmp = Float64(a / fma(10.0, k, 1.0)); else tmp = Float64(Float64(Float64(k * a) * k) * 99.0); end return tmp end
code[a_, k_, m_] := If[LessEqual[m, -6.2e-16], N[(a / N[(k * k), $MachinePrecision]), $MachinePrecision], If[LessEqual[m, 1.35], N[(a / N[(10.0 * k + 1.0), $MachinePrecision]), $MachinePrecision], N[(N[(N[(k * a), $MachinePrecision] * k), $MachinePrecision] * 99.0), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq -6.2 \cdot 10^{-16}:\\
\;\;\;\;\frac{a}{k \cdot k}\\
\mathbf{elif}\;m \leq 1.35:\\
\;\;\;\;\frac{a}{\mathsf{fma}\left(10, k, 1\right)}\\
\mathbf{else}:\\
\;\;\;\;\left(\left(k \cdot a\right) \cdot k\right) \cdot 99\\
\end{array}
\end{array}
if m < -6.2000000000000002e-16Initial program 98.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 rewrites43.0%
Taylor expanded in k around inf
Applied rewrites63.5%
if -6.2000000000000002e-16 < m < 1.3500000000000001Initial program 93.3%
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 rewrites90.9%
Taylor expanded in k around 0
Applied rewrites62.8%
if 1.3500000000000001 < m Initial program 76.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 rewrites2.8%
Taylor expanded in k around 0
Applied rewrites25.7%
Taylor expanded in k around inf
Applied rewrites54.6%
Final simplification59.9%
(FPCore (a k m) :precision binary64 (if (<= m -5.4e-74) (/ a (* k k)) (if (<= m 0.56) (fma (fma 99.0 k -10.0) (* k a) a) (* (* (* k a) k) 99.0))))
double code(double a, double k, double m) {
double tmp;
if (m <= -5.4e-74) {
tmp = a / (k * k);
} else if (m <= 0.56) {
tmp = fma(fma(99.0, k, -10.0), (k * a), a);
} else {
tmp = ((k * a) * k) * 99.0;
}
return tmp;
}
function code(a, k, m) tmp = 0.0 if (m <= -5.4e-74) tmp = Float64(a / Float64(k * k)); elseif (m <= 0.56) tmp = fma(fma(99.0, k, -10.0), Float64(k * a), a); else tmp = Float64(Float64(Float64(k * a) * k) * 99.0); end return tmp end
code[a_, k_, m_] := If[LessEqual[m, -5.4e-74], N[(a / N[(k * k), $MachinePrecision]), $MachinePrecision], If[LessEqual[m, 0.56], N[(N[(99.0 * k + -10.0), $MachinePrecision] * N[(k * a), $MachinePrecision] + a), $MachinePrecision], N[(N[(N[(k * a), $MachinePrecision] * k), $MachinePrecision] * 99.0), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq -5.4 \cdot 10^{-74}:\\
\;\;\;\;\frac{a}{k \cdot k}\\
\mathbf{elif}\;m \leq 0.56:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(99, k, -10\right), k \cdot a, a\right)\\
\mathbf{else}:\\
\;\;\;\;\left(\left(k \cdot a\right) \cdot k\right) \cdot 99\\
\end{array}
\end{array}
if m < -5.40000000000000036e-74Initial program 98.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 rewrites47.2%
Taylor expanded in k around inf
Applied rewrites63.7%
if -5.40000000000000036e-74 < m < 0.56000000000000005Initial program 92.8%
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 rewrites90.2%
Taylor expanded in k around 0
Applied rewrites51.5%
Applied rewrites51.5%
if 0.56000000000000005 < m Initial program 76.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 rewrites2.8%
Taylor expanded in k around 0
Applied rewrites25.7%
Taylor expanded in k around inf
Applied rewrites54.6%
Final simplification56.6%
(FPCore (a k m) :precision binary64 (if (<= m -5.4e-74) (/ a (* k k)) (if (<= m 0.56) (fma (* (fma 99.0 k -10.0) k) a a) (* (* (* k a) k) 99.0))))
double code(double a, double k, double m) {
double tmp;
if (m <= -5.4e-74) {
tmp = a / (k * k);
} else if (m <= 0.56) {
tmp = fma((fma(99.0, k, -10.0) * k), a, a);
} else {
tmp = ((k * a) * k) * 99.0;
}
return tmp;
}
function code(a, k, m) tmp = 0.0 if (m <= -5.4e-74) tmp = Float64(a / Float64(k * k)); elseif (m <= 0.56) tmp = fma(Float64(fma(99.0, k, -10.0) * k), a, a); else tmp = Float64(Float64(Float64(k * a) * k) * 99.0); end return tmp end
code[a_, k_, m_] := If[LessEqual[m, -5.4e-74], N[(a / N[(k * k), $MachinePrecision]), $MachinePrecision], If[LessEqual[m, 0.56], N[(N[(N[(99.0 * k + -10.0), $MachinePrecision] * k), $MachinePrecision] * a + a), $MachinePrecision], N[(N[(N[(k * a), $MachinePrecision] * k), $MachinePrecision] * 99.0), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq -5.4 \cdot 10^{-74}:\\
\;\;\;\;\frac{a}{k \cdot k}\\
\mathbf{elif}\;m \leq 0.56:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(99, k, -10\right) \cdot k, a, a\right)\\
\mathbf{else}:\\
\;\;\;\;\left(\left(k \cdot a\right) \cdot k\right) \cdot 99\\
\end{array}
\end{array}
if m < -5.40000000000000036e-74Initial program 98.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 rewrites47.2%
Taylor expanded in k around inf
Applied rewrites63.7%
if -5.40000000000000036e-74 < m < 0.56000000000000005Initial program 92.8%
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 rewrites90.2%
Taylor expanded in k around 0
Applied rewrites51.5%
Applied rewrites51.5%
if 0.56000000000000005 < m Initial program 76.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 rewrites2.8%
Taylor expanded in k around 0
Applied rewrites25.7%
Taylor expanded in k around inf
Applied rewrites54.6%
Final simplification56.5%
(FPCore (a k m) :precision binary64 (if (<= m -2e-67) (/ a (* k k)) (if (<= m 0.56) (* (- a) -1.0) (* (* (* k a) k) 99.0))))
double code(double a, double k, double m) {
double tmp;
if (m <= -2e-67) {
tmp = a / (k * k);
} else if (m <= 0.56) {
tmp = -a * -1.0;
} else {
tmp = ((k * a) * 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 <= (-2d-67)) then
tmp = a / (k * k)
else if (m <= 0.56d0) then
tmp = -a * (-1.0d0)
else
tmp = ((k * a) * k) * 99.0d0
end if
code = tmp
end function
public static double code(double a, double k, double m) {
double tmp;
if (m <= -2e-67) {
tmp = a / (k * k);
} else if (m <= 0.56) {
tmp = -a * -1.0;
} else {
tmp = ((k * a) * k) * 99.0;
}
return tmp;
}
def code(a, k, m): tmp = 0 if m <= -2e-67: tmp = a / (k * k) elif m <= 0.56: tmp = -a * -1.0 else: tmp = ((k * a) * k) * 99.0 return tmp
function code(a, k, m) tmp = 0.0 if (m <= -2e-67) tmp = Float64(a / Float64(k * k)); elseif (m <= 0.56) tmp = Float64(Float64(-a) * -1.0); else tmp = Float64(Float64(Float64(k * a) * k) * 99.0); end return tmp end
function tmp_2 = code(a, k, m) tmp = 0.0; if (m <= -2e-67) tmp = a / (k * k); elseif (m <= 0.56) tmp = -a * -1.0; else tmp = ((k * a) * k) * 99.0; end tmp_2 = tmp; end
code[a_, k_, m_] := If[LessEqual[m, -2e-67], N[(a / N[(k * k), $MachinePrecision]), $MachinePrecision], If[LessEqual[m, 0.56], N[((-a) * -1.0), $MachinePrecision], N[(N[(N[(k * a), $MachinePrecision] * k), $MachinePrecision] * 99.0), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq -2 \cdot 10^{-67}:\\
\;\;\;\;\frac{a}{k \cdot k}\\
\mathbf{elif}\;m \leq 0.56:\\
\;\;\;\;\left(-a\right) \cdot -1\\
\mathbf{else}:\\
\;\;\;\;\left(\left(k \cdot a\right) \cdot k\right) \cdot 99\\
\end{array}
\end{array}
if m < -1.99999999999999989e-67Initial program 98.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 rewrites47.2%
Taylor expanded in k around inf
Applied rewrites63.7%
if -1.99999999999999989e-67 < m < 0.56000000000000005Initial program 92.8%
lift-/.f64N/A
frac-2negN/A
div-invN/A
lift-*.f64N/A
distribute-lft-neg-inN/A
associate-*l*N/A
lower-*.f64N/A
lower-neg.f64N/A
lower-*.f64N/A
neg-mul-1N/A
associate-/r*N/A
metadata-evalN/A
lower-/.f6492.9
lift-+.f64N/A
lift-+.f64N/A
associate-+l+N/A
+-commutativeN/A
Applied rewrites92.9%
Taylor expanded in k around 0
mul-1-negN/A
lower-neg.f64N/A
lower-pow.f6453.5
Applied rewrites53.5%
Taylor expanded in m around 0
Applied rewrites51.0%
if 0.56000000000000005 < m Initial program 76.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 rewrites2.8%
Taylor expanded in k around 0
Applied rewrites25.7%
Taylor expanded in k around inf
Applied rewrites54.6%
Final simplification56.4%
(FPCore (a k m) :precision binary64 (if (<= m 0.56) (* (- a) -1.0) (* (* (* k a) k) 99.0)))
double code(double a, double k, double m) {
double tmp;
if (m <= 0.56) {
tmp = -a * -1.0;
} else {
tmp = ((k * a) * 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.56d0) then
tmp = -a * (-1.0d0)
else
tmp = ((k * a) * k) * 99.0d0
end if
code = tmp
end function
public static double code(double a, double k, double m) {
double tmp;
if (m <= 0.56) {
tmp = -a * -1.0;
} else {
tmp = ((k * a) * k) * 99.0;
}
return tmp;
}
def code(a, k, m): tmp = 0 if m <= 0.56: tmp = -a * -1.0 else: tmp = ((k * a) * k) * 99.0 return tmp
function code(a, k, m) tmp = 0.0 if (m <= 0.56) tmp = Float64(Float64(-a) * -1.0); else tmp = Float64(Float64(Float64(k * a) * k) * 99.0); end return tmp end
function tmp_2 = code(a, k, m) tmp = 0.0; if (m <= 0.56) tmp = -a * -1.0; else tmp = ((k * a) * k) * 99.0; end tmp_2 = tmp; end
code[a_, k_, m_] := If[LessEqual[m, 0.56], N[((-a) * -1.0), $MachinePrecision], N[(N[(N[(k * a), $MachinePrecision] * k), $MachinePrecision] * 99.0), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq 0.56:\\
\;\;\;\;\left(-a\right) \cdot -1\\
\mathbf{else}:\\
\;\;\;\;\left(\left(k \cdot a\right) \cdot k\right) \cdot 99\\
\end{array}
\end{array}
if m < 0.56000000000000005Initial program 95.8%
lift-/.f64N/A
frac-2negN/A
div-invN/A
lift-*.f64N/A
distribute-lft-neg-inN/A
associate-*l*N/A
lower-*.f64N/A
lower-neg.f64N/A
lower-*.f64N/A
neg-mul-1N/A
associate-/r*N/A
metadata-evalN/A
lower-/.f6495.9
lift-+.f64N/A
lift-+.f64N/A
associate-+l+N/A
+-commutativeN/A
Applied rewrites96.5%
Taylor expanded in k around 0
mul-1-negN/A
lower-neg.f64N/A
lower-pow.f6474.1
Applied rewrites74.1%
Taylor expanded in m around 0
Applied rewrites28.4%
if 0.56000000000000005 < m Initial program 76.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 rewrites2.8%
Taylor expanded in k around 0
Applied rewrites25.7%
Taylor expanded in k around inf
Applied rewrites54.6%
Final simplification38.2%
(FPCore (a k m) :precision binary64 (if (<= m 1.46e+16) (* (- a) -1.0) (* (* -10.0 a) k)))
double code(double a, double k, double m) {
double tmp;
if (m <= 1.46e+16) {
tmp = -a * -1.0;
} else {
tmp = (-10.0 * a) * k;
}
return tmp;
}
real(8) function code(a, k, m)
real(8), intent (in) :: a
real(8), intent (in) :: k
real(8), intent (in) :: m
real(8) :: tmp
if (m <= 1.46d+16) then
tmp = -a * (-1.0d0)
else
tmp = ((-10.0d0) * a) * k
end if
code = tmp
end function
public static double code(double a, double k, double m) {
double tmp;
if (m <= 1.46e+16) {
tmp = -a * -1.0;
} else {
tmp = (-10.0 * a) * k;
}
return tmp;
}
def code(a, k, m): tmp = 0 if m <= 1.46e+16: tmp = -a * -1.0 else: tmp = (-10.0 * a) * k return tmp
function code(a, k, m) tmp = 0.0 if (m <= 1.46e+16) tmp = Float64(Float64(-a) * -1.0); else tmp = Float64(Float64(-10.0 * a) * k); end return tmp end
function tmp_2 = code(a, k, m) tmp = 0.0; if (m <= 1.46e+16) tmp = -a * -1.0; else tmp = (-10.0 * a) * k; end tmp_2 = tmp; end
code[a_, k_, m_] := If[LessEqual[m, 1.46e+16], N[((-a) * -1.0), $MachinePrecision], N[(N[(-10.0 * a), $MachinePrecision] * k), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq 1.46 \cdot 10^{+16}:\\
\;\;\;\;\left(-a\right) \cdot -1\\
\mathbf{else}:\\
\;\;\;\;\left(-10 \cdot a\right) \cdot k\\
\end{array}
\end{array}
if m < 1.46e16Initial program 94.6%
lift-/.f64N/A
frac-2negN/A
div-invN/A
lift-*.f64N/A
distribute-lft-neg-inN/A
associate-*l*N/A
lower-*.f64N/A
lower-neg.f64N/A
lower-*.f64N/A
neg-mul-1N/A
associate-/r*N/A
metadata-evalN/A
lower-/.f6494.7
lift-+.f64N/A
lift-+.f64N/A
associate-+l+N/A
+-commutativeN/A
Applied rewrites95.3%
Taylor expanded in k around 0
mul-1-negN/A
lower-neg.f64N/A
lower-pow.f6474.5
Applied rewrites74.5%
Taylor expanded in m around 0
Applied rewrites28.0%
if 1.46e16 < m Initial program 77.4%
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.9%
Taylor expanded in k around 0
Applied rewrites5.9%
Taylor expanded in k around inf
Applied rewrites21.7%
Final simplification25.7%
(FPCore (a k m) :precision binary64 (* (- a) -1.0))
double code(double a, double k, double m) {
return -a * -1.0;
}
real(8) function code(a, k, m)
real(8), intent (in) :: a
real(8), intent (in) :: k
real(8), intent (in) :: m
code = -a * (-1.0d0)
end function
public static double code(double a, double k, double m) {
return -a * -1.0;
}
def code(a, k, m): return -a * -1.0
function code(a, k, m) return Float64(Float64(-a) * -1.0) end
function tmp = code(a, k, m) tmp = -a * -1.0; end
code[a_, k_, m_] := N[((-a) * -1.0), $MachinePrecision]
\begin{array}{l}
\\
\left(-a\right) \cdot -1
\end{array}
Initial program 88.4%
lift-/.f64N/A
frac-2negN/A
div-invN/A
lift-*.f64N/A
distribute-lft-neg-inN/A
associate-*l*N/A
lower-*.f64N/A
lower-neg.f64N/A
lower-*.f64N/A
neg-mul-1N/A
associate-/r*N/A
metadata-evalN/A
lower-/.f6488.4
lift-+.f64N/A
lift-+.f64N/A
associate-+l+N/A
+-commutativeN/A
Applied rewrites88.8%
Taylor expanded in k around 0
mul-1-negN/A
lower-neg.f64N/A
lower-pow.f6483.8
Applied rewrites83.8%
Taylor expanded in m around 0
Applied rewrites19.2%
herbie shell --seed 2024314
(FPCore (a k m)
:name "Falkner and Boettcher, Appendix A"
:precision binary64
(/ (* a (pow k m)) (+ (+ 1.0 (* 10.0 k)) (* k k))))