
(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 (<= k 1.6e-11) (* (pow k m) a) (/ (* a (/ (pow k m) k)) k)))
double code(double a, double k, double m) {
double tmp;
if (k <= 1.6e-11) {
tmp = pow(k, m) * a;
} else {
tmp = (a * (pow(k, m) / k)) / 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 <= 1.6d-11) then
tmp = (k ** m) * a
else
tmp = (a * ((k ** m) / k)) / k
end if
code = tmp
end function
public static double code(double a, double k, double m) {
double tmp;
if (k <= 1.6e-11) {
tmp = Math.pow(k, m) * a;
} else {
tmp = (a * (Math.pow(k, m) / k)) / k;
}
return tmp;
}
def code(a, k, m): tmp = 0 if k <= 1.6e-11: tmp = math.pow(k, m) * a else: tmp = (a * (math.pow(k, m) / k)) / k return tmp
function code(a, k, m) tmp = 0.0 if (k <= 1.6e-11) tmp = Float64((k ^ m) * a); else tmp = Float64(Float64(a * Float64((k ^ m) / k)) / k); end return tmp end
function tmp_2 = code(a, k, m) tmp = 0.0; if (k <= 1.6e-11) tmp = (k ^ m) * a; else tmp = (a * ((k ^ m) / k)) / k; end tmp_2 = tmp; end
code[a_, k_, m_] := If[LessEqual[k, 1.6e-11], N[(N[Power[k, m], $MachinePrecision] * a), $MachinePrecision], N[(N[(a * N[(N[Power[k, m], $MachinePrecision] / k), $MachinePrecision]), $MachinePrecision] / k), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;k \leq 1.6 \cdot 10^{-11}:\\
\;\;\;\;{k}^{m} \cdot a\\
\mathbf{else}:\\
\;\;\;\;\frac{a \cdot \frac{{k}^{m}}{k}}{k}\\
\end{array}
\end{array}
if k < 1.59999999999999997e-11Initial program 92.0%
Taylor expanded in k around 0
rem-exp-logN/A
remove-double-negN/A
log-recN/A
exp-prodN/A
distribute-lft-neg-inN/A
*-commutativeN/A
mul-1-negN/A
+-rgt-identityN/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
Simplified99.1%
+-rgt-identityN/A
*-commutativeN/A
*-lowering-*.f64N/A
pow-lowering-pow.f6499.1
Applied egg-rr99.1%
if 1.59999999999999997e-11 < k Initial program 82.4%
Taylor expanded in k around inf
unpow2N/A
*-lowering-*.f6481.3
Simplified81.3%
*-commutativeN/A
times-fracN/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
pow-lowering-pow.f6498.8
Applied egg-rr98.8%
Final simplification99.0%
(FPCore (a k m) :precision binary64 (if (<= k 1.6e-11) (* (pow k m) a) (/ (* a (pow k (+ m -1.0))) k)))
double code(double a, double k, double m) {
double tmp;
if (k <= 1.6e-11) {
tmp = pow(k, m) * a;
} else {
tmp = (a * pow(k, (m + -1.0))) / k;
}
return tmp;
}
real(8) function code(a, k, m)
real(8), intent (in) :: a
real(8), intent (in) :: k
real(8), intent (in) :: m
real(8) :: tmp
if (k <= 1.6d-11) then
tmp = (k ** m) * a
else
tmp = (a * (k ** (m + (-1.0d0)))) / k
end if
code = tmp
end function
public static double code(double a, double k, double m) {
double tmp;
if (k <= 1.6e-11) {
tmp = Math.pow(k, m) * a;
} else {
tmp = (a * Math.pow(k, (m + -1.0))) / k;
}
return tmp;
}
def code(a, k, m): tmp = 0 if k <= 1.6e-11: tmp = math.pow(k, m) * a else: tmp = (a * math.pow(k, (m + -1.0))) / k return tmp
function code(a, k, m) tmp = 0.0 if (k <= 1.6e-11) tmp = Float64((k ^ m) * a); else tmp = Float64(Float64(a * (k ^ Float64(m + -1.0))) / k); end return tmp end
function tmp_2 = code(a, k, m) tmp = 0.0; if (k <= 1.6e-11) tmp = (k ^ m) * a; else tmp = (a * (k ^ (m + -1.0))) / k; end tmp_2 = tmp; end
code[a_, k_, m_] := If[LessEqual[k, 1.6e-11], N[(N[Power[k, m], $MachinePrecision] * a), $MachinePrecision], N[(N[(a * N[Power[k, N[(m + -1.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / k), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;k \leq 1.6 \cdot 10^{-11}:\\
\;\;\;\;{k}^{m} \cdot a\\
\mathbf{else}:\\
\;\;\;\;\frac{a \cdot {k}^{\left(m + -1\right)}}{k}\\
\end{array}
\end{array}
if k < 1.59999999999999997e-11Initial program 92.0%
Taylor expanded in k around 0
rem-exp-logN/A
remove-double-negN/A
log-recN/A
exp-prodN/A
distribute-lft-neg-inN/A
*-commutativeN/A
mul-1-negN/A
+-rgt-identityN/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
Simplified99.1%
+-rgt-identityN/A
*-commutativeN/A
*-lowering-*.f64N/A
pow-lowering-pow.f6499.1
Applied egg-rr99.1%
if 1.59999999999999997e-11 < k Initial program 82.4%
Taylor expanded in k around inf
unpow2N/A
*-lowering-*.f6481.3
Simplified81.3%
*-commutativeN/A
times-fracN/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
pow-lowering-pow.f6498.8
Applied egg-rr98.8%
/-lowering-/.f64N/A
*-lowering-*.f64N/A
div-invN/A
inv-powN/A
pow-prod-upN/A
pow-lowering-pow.f64N/A
+-lowering-+.f6498.6
Applied egg-rr98.6%
Final simplification98.9%
(FPCore (a k m)
:precision binary64
(let* ((t_0 (* (pow k m) a)))
(if (<= m -2.75e-17)
t_0
(if (<= m 0.98) (/ a (fma k (+ k 10.0) 1.0)) t_0))))
double code(double a, double k, double m) {
double t_0 = pow(k, m) * a;
double tmp;
if (m <= -2.75e-17) {
tmp = t_0;
} else if (m <= 0.98) {
tmp = a / fma(k, (k + 10.0), 1.0);
} else {
tmp = t_0;
}
return tmp;
}
function code(a, k, m) t_0 = Float64((k ^ m) * a) tmp = 0.0 if (m <= -2.75e-17) tmp = t_0; elseif (m <= 0.98) tmp = Float64(a / fma(k, Float64(k + 10.0), 1.0)); else tmp = t_0; end return tmp end
code[a_, k_, m_] := Block[{t$95$0 = N[(N[Power[k, m], $MachinePrecision] * a), $MachinePrecision]}, If[LessEqual[m, -2.75e-17], t$95$0, If[LessEqual[m, 0.98], N[(a / N[(k * N[(k + 10.0), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision], t$95$0]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := {k}^{m} \cdot a\\
\mathbf{if}\;m \leq -2.75 \cdot 10^{-17}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;m \leq 0.98:\\
\;\;\;\;\frac{a}{\mathsf{fma}\left(k, k + 10, 1\right)}\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if m < -2.75e-17 or 0.97999999999999998 < m Initial program 84.7%
Taylor expanded in k around 0
rem-exp-logN/A
remove-double-negN/A
log-recN/A
exp-prodN/A
distribute-lft-neg-inN/A
*-commutativeN/A
mul-1-negN/A
+-rgt-identityN/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
Simplified99.4%
+-rgt-identityN/A
*-commutativeN/A
*-lowering-*.f64N/A
pow-lowering-pow.f6499.4
Applied egg-rr99.4%
if -2.75e-17 < m < 0.97999999999999998Initial program 93.3%
Taylor expanded in m around 0
/-lowering-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
+-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
+-lowering-+.f6491.5
Simplified91.5%
Final simplification96.3%
(FPCore (a k m) :precision binary64 (if (<= k 1.6e-11) (* (pow k m) a) (* a (pow k (+ m -2.0)))))
double code(double a, double k, double m) {
double tmp;
if (k <= 1.6e-11) {
tmp = pow(k, m) * a;
} else {
tmp = a * pow(k, (m + -2.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 (k <= 1.6d-11) then
tmp = (k ** m) * a
else
tmp = a * (k ** (m + (-2.0d0)))
end if
code = tmp
end function
public static double code(double a, double k, double m) {
double tmp;
if (k <= 1.6e-11) {
tmp = Math.pow(k, m) * a;
} else {
tmp = a * Math.pow(k, (m + -2.0));
}
return tmp;
}
def code(a, k, m): tmp = 0 if k <= 1.6e-11: tmp = math.pow(k, m) * a else: tmp = a * math.pow(k, (m + -2.0)) return tmp
function code(a, k, m) tmp = 0.0 if (k <= 1.6e-11) tmp = Float64((k ^ m) * a); else tmp = Float64(a * (k ^ Float64(m + -2.0))); end return tmp end
function tmp_2 = code(a, k, m) tmp = 0.0; if (k <= 1.6e-11) tmp = (k ^ m) * a; else tmp = a * (k ^ (m + -2.0)); end tmp_2 = tmp; end
code[a_, k_, m_] := If[LessEqual[k, 1.6e-11], N[(N[Power[k, m], $MachinePrecision] * a), $MachinePrecision], N[(a * N[Power[k, N[(m + -2.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;k \leq 1.6 \cdot 10^{-11}:\\
\;\;\;\;{k}^{m} \cdot a\\
\mathbf{else}:\\
\;\;\;\;a \cdot {k}^{\left(m + -2\right)}\\
\end{array}
\end{array}
if k < 1.59999999999999997e-11Initial program 92.0%
Taylor expanded in k around 0
rem-exp-logN/A
remove-double-negN/A
log-recN/A
exp-prodN/A
distribute-lft-neg-inN/A
*-commutativeN/A
mul-1-negN/A
+-rgt-identityN/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
Simplified99.1%
+-rgt-identityN/A
*-commutativeN/A
*-lowering-*.f64N/A
pow-lowering-pow.f6499.1
Applied egg-rr99.1%
if 1.59999999999999997e-11 < k Initial program 82.4%
Taylor expanded in k around inf
unpow2N/A
*-lowering-*.f6481.3
Simplified81.3%
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
pow2N/A
pow-divN/A
pow-lowering-pow.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-eval92.8
Applied egg-rr92.8%
Final simplification96.5%
(FPCore (a k m)
:precision binary64
(if (<= m -1.52)
(/
(fma 9801.0 (/ a (* (* k k) (* k k))) (fma (/ a (* k k)) 99.0 a))
(* k k))
(if (<= m 1.4) (/ a (fma k (+ k 10.0) 1.0)) (* a (* (* k k) 99.0)))))
double code(double a, double k, double m) {
double tmp;
if (m <= -1.52) {
tmp = fma(9801.0, (a / ((k * k) * (k * k))), fma((a / (k * k)), 99.0, a)) / (k * k);
} else if (m <= 1.4) {
tmp = a / fma(k, (k + 10.0), 1.0);
} else {
tmp = a * ((k * k) * 99.0);
}
return tmp;
}
function code(a, k, m) tmp = 0.0 if (m <= -1.52) tmp = Float64(fma(9801.0, Float64(a / Float64(Float64(k * k) * Float64(k * k))), fma(Float64(a / Float64(k * k)), 99.0, a)) / Float64(k * k)); elseif (m <= 1.4) tmp = Float64(a / fma(k, Float64(k + 10.0), 1.0)); else tmp = Float64(a * Float64(Float64(k * k) * 99.0)); end return tmp end
code[a_, k_, m_] := If[LessEqual[m, -1.52], N[(N[(9801.0 * N[(a / N[(N[(k * k), $MachinePrecision] * N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(a / N[(k * k), $MachinePrecision]), $MachinePrecision] * 99.0 + a), $MachinePrecision]), $MachinePrecision] / N[(k * k), $MachinePrecision]), $MachinePrecision], If[LessEqual[m, 1.4], N[(a / N[(k * N[(k + 10.0), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision], N[(a * N[(N[(k * k), $MachinePrecision] * 99.0), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq -1.52:\\
\;\;\;\;\frac{\mathsf{fma}\left(9801, \frac{a}{\left(k \cdot k\right) \cdot \left(k \cdot k\right)}, \mathsf{fma}\left(\frac{a}{k \cdot k}, 99, a\right)\right)}{k \cdot k}\\
\mathbf{elif}\;m \leq 1.4:\\
\;\;\;\;\frac{a}{\mathsf{fma}\left(k, k + 10, 1\right)}\\
\mathbf{else}:\\
\;\;\;\;a \cdot \left(\left(k \cdot k\right) \cdot 99\right)\\
\end{array}
\end{array}
if m < -1.52Initial program 100.0%
Taylor expanded in m around 0
/-lowering-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
+-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
+-lowering-+.f6437.8
Simplified37.8%
*-commutativeN/A
flip-+N/A
associate-*l/N/A
div-invN/A
accelerator-lowering-fma.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
metadata-evalN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
--lowering--.f6441.6
Applied egg-rr41.6%
Taylor expanded in k around inf
/-lowering-/.f6440.3
Simplified40.3%
Taylor expanded in k around inf
Simplified69.3%
if -1.52 < m < 1.3999999999999999Initial program 93.6%
Taylor expanded in m around 0
/-lowering-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
+-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
+-lowering-+.f6490.7
Simplified90.7%
if 1.3999999999999999 < m Initial program 70.0%
Taylor expanded in m around 0
/-lowering-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
+-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
+-lowering-+.f642.6
Simplified2.6%
Taylor expanded in k around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
sub-negN/A
mul-1-negN/A
distribute-neg-outN/A
neg-sub0N/A
--lowering--.f64N/A
distribute-rgt1-inN/A
associate-*r*N/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
accelerator-lowering-fma.f64N/A
metadata-eval29.4
Simplified29.4%
Taylor expanded in k around inf
*-commutativeN/A
associate-*l*N/A
metadata-evalN/A
distribute-rgt-neg-inN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
distribute-rgt-neg-inN/A
metadata-evalN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6460.7
Simplified60.7%
Final simplification75.2%
(FPCore (a k m) :precision binary64 (if (<= m -1.52) (* a (/ (- 1.0 (/ (+ 10.0 (/ -99.0 k)) k)) (* k k))) (if (<= m 1.35) (/ a (fma k (+ k 10.0) 1.0)) (* a (* (* k k) 99.0)))))
double code(double a, double k, double m) {
double tmp;
if (m <= -1.52) {
tmp = a * ((1.0 - ((10.0 + (-99.0 / k)) / k)) / (k * k));
} else if (m <= 1.35) {
tmp = a / fma(k, (k + 10.0), 1.0);
} else {
tmp = a * ((k * k) * 99.0);
}
return tmp;
}
function code(a, k, m) tmp = 0.0 if (m <= -1.52) tmp = Float64(a * Float64(Float64(1.0 - Float64(Float64(10.0 + Float64(-99.0 / k)) / k)) / Float64(k * k))); elseif (m <= 1.35) tmp = Float64(a / fma(k, Float64(k + 10.0), 1.0)); else tmp = Float64(a * Float64(Float64(k * k) * 99.0)); end return tmp end
code[a_, k_, m_] := If[LessEqual[m, -1.52], N[(a * N[(N[(1.0 - N[(N[(10.0 + N[(-99.0 / k), $MachinePrecision]), $MachinePrecision] / k), $MachinePrecision]), $MachinePrecision] / N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[m, 1.35], N[(a / N[(k * N[(k + 10.0), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision], N[(a * N[(N[(k * k), $MachinePrecision] * 99.0), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq -1.52:\\
\;\;\;\;a \cdot \frac{1 - \frac{10 + \frac{-99}{k}}{k}}{k \cdot k}\\
\mathbf{elif}\;m \leq 1.35:\\
\;\;\;\;\frac{a}{\mathsf{fma}\left(k, k + 10, 1\right)}\\
\mathbf{else}:\\
\;\;\;\;a \cdot \left(\left(k \cdot k\right) \cdot 99\right)\\
\end{array}
\end{array}
if m < -1.52Initial program 100.0%
Taylor expanded in m around 0
/-lowering-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
+-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
+-lowering-+.f6437.8
Simplified37.8%
clear-numN/A
associate-/r/N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
+-lowering-+.f6437.8
Applied egg-rr37.8%
Taylor expanded in k around -inf
/-lowering-/.f64N/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
/-lowering-/.f64N/A
sub-negN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
distribute-neg-fracN/A
metadata-evalN/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6467.9
Simplified67.9%
if -1.52 < m < 1.3500000000000001Initial program 93.6%
Taylor expanded in m around 0
/-lowering-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
+-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
+-lowering-+.f6490.7
Simplified90.7%
if 1.3500000000000001 < m Initial program 70.0%
Taylor expanded in m around 0
/-lowering-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
+-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
+-lowering-+.f642.6
Simplified2.6%
Taylor expanded in k around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
sub-negN/A
mul-1-negN/A
distribute-neg-outN/A
neg-sub0N/A
--lowering--.f64N/A
distribute-rgt1-inN/A
associate-*r*N/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
accelerator-lowering-fma.f64N/A
metadata-eval29.4
Simplified29.4%
Taylor expanded in k around inf
*-commutativeN/A
associate-*l*N/A
metadata-evalN/A
distribute-rgt-neg-inN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
distribute-rgt-neg-inN/A
metadata-evalN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6460.7
Simplified60.7%
Final simplification74.8%
(FPCore (a k m) :precision binary64 (if (<= m -1.52) (/ (fma (/ a k) (- (/ 99.0 k) 10.0) a) (* k k)) (if (<= m 1.0) (/ a (fma k (+ k 10.0) 1.0)) (* a (* (* k k) 99.0)))))
double code(double a, double k, double m) {
double tmp;
if (m <= -1.52) {
tmp = fma((a / k), ((99.0 / k) - 10.0), a) / (k * k);
} else if (m <= 1.0) {
tmp = a / fma(k, (k + 10.0), 1.0);
} else {
tmp = a * ((k * k) * 99.0);
}
return tmp;
}
function code(a, k, m) tmp = 0.0 if (m <= -1.52) tmp = Float64(fma(Float64(a / k), Float64(Float64(99.0 / k) - 10.0), a) / Float64(k * k)); elseif (m <= 1.0) tmp = Float64(a / fma(k, Float64(k + 10.0), 1.0)); else tmp = Float64(a * Float64(Float64(k * k) * 99.0)); end return tmp end
code[a_, k_, m_] := If[LessEqual[m, -1.52], N[(N[(N[(a / k), $MachinePrecision] * N[(N[(99.0 / k), $MachinePrecision] - 10.0), $MachinePrecision] + a), $MachinePrecision] / N[(k * k), $MachinePrecision]), $MachinePrecision], If[LessEqual[m, 1.0], N[(a / N[(k * N[(k + 10.0), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision], N[(a * N[(N[(k * k), $MachinePrecision] * 99.0), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq -1.52:\\
\;\;\;\;\frac{\mathsf{fma}\left(\frac{a}{k}, \frac{99}{k} - 10, a\right)}{k \cdot k}\\
\mathbf{elif}\;m \leq 1:\\
\;\;\;\;\frac{a}{\mathsf{fma}\left(k, k + 10, 1\right)}\\
\mathbf{else}:\\
\;\;\;\;a \cdot \left(\left(k \cdot k\right) \cdot 99\right)\\
\end{array}
\end{array}
if m < -1.52Initial program 100.0%
Taylor expanded in m around 0
/-lowering-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
+-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
+-lowering-+.f6437.8
Simplified37.8%
clear-numN/A
associate-/r/N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
+-lowering-+.f6437.8
Applied egg-rr37.8%
Taylor expanded in k around inf
/-lowering-/.f64N/A
Simplified65.5%
if -1.52 < m < 1Initial program 93.6%
Taylor expanded in m around 0
/-lowering-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
+-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
+-lowering-+.f6490.7
Simplified90.7%
if 1 < m Initial program 70.0%
Taylor expanded in m around 0
/-lowering-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
+-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
+-lowering-+.f642.6
Simplified2.6%
Taylor expanded in k around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
sub-negN/A
mul-1-negN/A
distribute-neg-outN/A
neg-sub0N/A
--lowering--.f64N/A
distribute-rgt1-inN/A
associate-*r*N/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
accelerator-lowering-fma.f64N/A
metadata-eval29.4
Simplified29.4%
Taylor expanded in k around inf
*-commutativeN/A
associate-*l*N/A
metadata-evalN/A
distribute-rgt-neg-inN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
distribute-rgt-neg-inN/A
metadata-evalN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6460.7
Simplified60.7%
Final simplification74.2%
(FPCore (a k m) :precision binary64 (if (<= m -1.52) (/ (fma (/ a (* k k)) 99.0 a) (* k k)) (if (<= m 1.25) (/ a (fma k (+ k 10.0) 1.0)) (* a (* (* k k) 99.0)))))
double code(double a, double k, double m) {
double tmp;
if (m <= -1.52) {
tmp = fma((a / (k * k)), 99.0, a) / (k * k);
} else if (m <= 1.25) {
tmp = a / fma(k, (k + 10.0), 1.0);
} else {
tmp = a * ((k * k) * 99.0);
}
return tmp;
}
function code(a, k, m) tmp = 0.0 if (m <= -1.52) tmp = Float64(fma(Float64(a / Float64(k * k)), 99.0, a) / Float64(k * k)); elseif (m <= 1.25) tmp = Float64(a / fma(k, Float64(k + 10.0), 1.0)); else tmp = Float64(a * Float64(Float64(k * k) * 99.0)); end return tmp end
code[a_, k_, m_] := If[LessEqual[m, -1.52], N[(N[(N[(a / N[(k * k), $MachinePrecision]), $MachinePrecision] * 99.0 + a), $MachinePrecision] / N[(k * k), $MachinePrecision]), $MachinePrecision], If[LessEqual[m, 1.25], N[(a / N[(k * N[(k + 10.0), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision], N[(a * N[(N[(k * k), $MachinePrecision] * 99.0), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq -1.52:\\
\;\;\;\;\frac{\mathsf{fma}\left(\frac{a}{k \cdot k}, 99, a\right)}{k \cdot k}\\
\mathbf{elif}\;m \leq 1.25:\\
\;\;\;\;\frac{a}{\mathsf{fma}\left(k, k + 10, 1\right)}\\
\mathbf{else}:\\
\;\;\;\;a \cdot \left(\left(k \cdot k\right) \cdot 99\right)\\
\end{array}
\end{array}
if m < -1.52Initial program 100.0%
Taylor expanded in m around 0
/-lowering-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
+-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
+-lowering-+.f6437.8
Simplified37.8%
*-commutativeN/A
flip-+N/A
associate-*l/N/A
div-invN/A
accelerator-lowering-fma.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
metadata-evalN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
--lowering--.f6441.6
Applied egg-rr41.6%
Taylor expanded in k around inf
/-lowering-/.f6440.3
Simplified40.3%
Taylor expanded in k around inf
metadata-evalN/A
associate-*r*N/A
associate-*r/N/A
metadata-evalN/A
distribute-rgt1-inN/A
/-lowering-/.f64N/A
Simplified65.5%
if -1.52 < m < 1.25Initial program 93.6%
Taylor expanded in m around 0
/-lowering-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
+-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
+-lowering-+.f6490.7
Simplified90.7%
if 1.25 < m Initial program 70.0%
Taylor expanded in m around 0
/-lowering-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
+-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
+-lowering-+.f642.6
Simplified2.6%
Taylor expanded in k around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
sub-negN/A
mul-1-negN/A
distribute-neg-outN/A
neg-sub0N/A
--lowering--.f64N/A
distribute-rgt1-inN/A
associate-*r*N/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
accelerator-lowering-fma.f64N/A
metadata-eval29.4
Simplified29.4%
Taylor expanded in k around inf
*-commutativeN/A
associate-*l*N/A
metadata-evalN/A
distribute-rgt-neg-inN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
distribute-rgt-neg-inN/A
metadata-evalN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6460.7
Simplified60.7%
Final simplification74.2%
(FPCore (a k m)
:precision binary64
(let* ((t_0 (/ a (* k k))))
(if (<= m -1.66e-195)
t_0
(if (<= m 5e-82) a (if (<= m 1.05) t_0 (* a (* (* k k) 99.0)))))))
double code(double a, double k, double m) {
double t_0 = a / (k * k);
double tmp;
if (m <= -1.66e-195) {
tmp = t_0;
} else if (m <= 5e-82) {
tmp = a;
} else if (m <= 1.05) {
tmp = t_0;
} else {
tmp = a * ((k * k) * 99.0);
}
return tmp;
}
real(8) function code(a, k, m)
real(8), intent (in) :: a
real(8), intent (in) :: k
real(8), intent (in) :: m
real(8) :: t_0
real(8) :: tmp
t_0 = a / (k * k)
if (m <= (-1.66d-195)) then
tmp = t_0
else if (m <= 5d-82) then
tmp = a
else if (m <= 1.05d0) then
tmp = t_0
else
tmp = a * ((k * k) * 99.0d0)
end if
code = tmp
end function
public static double code(double a, double k, double m) {
double t_0 = a / (k * k);
double tmp;
if (m <= -1.66e-195) {
tmp = t_0;
} else if (m <= 5e-82) {
tmp = a;
} else if (m <= 1.05) {
tmp = t_0;
} else {
tmp = a * ((k * k) * 99.0);
}
return tmp;
}
def code(a, k, m): t_0 = a / (k * k) tmp = 0 if m <= -1.66e-195: tmp = t_0 elif m <= 5e-82: tmp = a elif m <= 1.05: tmp = t_0 else: tmp = a * ((k * k) * 99.0) return tmp
function code(a, k, m) t_0 = Float64(a / Float64(k * k)) tmp = 0.0 if (m <= -1.66e-195) tmp = t_0; elseif (m <= 5e-82) tmp = a; elseif (m <= 1.05) tmp = t_0; else tmp = Float64(a * Float64(Float64(k * k) * 99.0)); end return tmp end
function tmp_2 = code(a, k, m) t_0 = a / (k * k); tmp = 0.0; if (m <= -1.66e-195) tmp = t_0; elseif (m <= 5e-82) tmp = a; elseif (m <= 1.05) tmp = t_0; else tmp = a * ((k * k) * 99.0); end tmp_2 = tmp; end
code[a_, k_, m_] := Block[{t$95$0 = N[(a / N[(k * k), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[m, -1.66e-195], t$95$0, If[LessEqual[m, 5e-82], a, If[LessEqual[m, 1.05], t$95$0, N[(a * N[(N[(k * k), $MachinePrecision] * 99.0), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{a}{k \cdot k}\\
\mathbf{if}\;m \leq -1.66 \cdot 10^{-195}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;m \leq 5 \cdot 10^{-82}:\\
\;\;\;\;a\\
\mathbf{elif}\;m \leq 1.05:\\
\;\;\;\;t\_0\\
\mathbf{else}:\\
\;\;\;\;a \cdot \left(\left(k \cdot k\right) \cdot 99\right)\\
\end{array}
\end{array}
if m < -1.66e-195 or 4.9999999999999998e-82 < m < 1.05000000000000004Initial program 97.6%
Taylor expanded in m around 0
/-lowering-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
+-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
+-lowering-+.f6456.9
Simplified56.9%
Taylor expanded in k around inf
unpow2N/A
*-lowering-*.f6462.4
Simplified62.4%
if -1.66e-195 < m < 4.9999999999999998e-82Initial program 93.5%
Taylor expanded in m around 0
/-lowering-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
+-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
+-lowering-+.f6493.5
Simplified93.5%
Taylor expanded in k around 0
Simplified56.6%
if 1.05000000000000004 < m Initial program 70.0%
Taylor expanded in m around 0
/-lowering-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
+-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
+-lowering-+.f642.6
Simplified2.6%
Taylor expanded in k around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
sub-negN/A
mul-1-negN/A
distribute-neg-outN/A
neg-sub0N/A
--lowering--.f64N/A
distribute-rgt1-inN/A
associate-*r*N/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
accelerator-lowering-fma.f64N/A
metadata-eval29.4
Simplified29.4%
Taylor expanded in k around inf
*-commutativeN/A
associate-*l*N/A
metadata-evalN/A
distribute-rgt-neg-inN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
distribute-rgt-neg-inN/A
metadata-evalN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6460.7
Simplified60.7%
(FPCore (a k m) :precision binary64 (if (<= m -1.52) (/ a (* k k)) (if (<= m 1.1) (/ a (fma k (+ k 10.0) 1.0)) (* a (* (* k k) 99.0)))))
double code(double a, double k, double m) {
double tmp;
if (m <= -1.52) {
tmp = a / (k * k);
} else if (m <= 1.1) {
tmp = a / fma(k, (k + 10.0), 1.0);
} else {
tmp = a * ((k * k) * 99.0);
}
return tmp;
}
function code(a, k, m) tmp = 0.0 if (m <= -1.52) tmp = Float64(a / Float64(k * k)); elseif (m <= 1.1) tmp = Float64(a / fma(k, Float64(k + 10.0), 1.0)); else tmp = Float64(a * Float64(Float64(k * k) * 99.0)); end return tmp end
code[a_, k_, m_] := If[LessEqual[m, -1.52], N[(a / N[(k * k), $MachinePrecision]), $MachinePrecision], If[LessEqual[m, 1.1], N[(a / N[(k * N[(k + 10.0), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision], N[(a * N[(N[(k * k), $MachinePrecision] * 99.0), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq -1.52:\\
\;\;\;\;\frac{a}{k \cdot k}\\
\mathbf{elif}\;m \leq 1.1:\\
\;\;\;\;\frac{a}{\mathsf{fma}\left(k, k + 10, 1\right)}\\
\mathbf{else}:\\
\;\;\;\;a \cdot \left(\left(k \cdot k\right) \cdot 99\right)\\
\end{array}
\end{array}
if m < -1.52Initial program 100.0%
Taylor expanded in m around 0
/-lowering-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
+-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
+-lowering-+.f6437.8
Simplified37.8%
Taylor expanded in k around inf
unpow2N/A
*-lowering-*.f6462.6
Simplified62.6%
if -1.52 < m < 1.1000000000000001Initial program 93.6%
Taylor expanded in m around 0
/-lowering-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
+-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
+-lowering-+.f6490.7
Simplified90.7%
if 1.1000000000000001 < m Initial program 70.0%
Taylor expanded in m around 0
/-lowering-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
+-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
+-lowering-+.f642.6
Simplified2.6%
Taylor expanded in k around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
sub-negN/A
mul-1-negN/A
distribute-neg-outN/A
neg-sub0N/A
--lowering--.f64N/A
distribute-rgt1-inN/A
associate-*r*N/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
accelerator-lowering-fma.f64N/A
metadata-eval29.4
Simplified29.4%
Taylor expanded in k around inf
*-commutativeN/A
associate-*l*N/A
metadata-evalN/A
distribute-rgt-neg-inN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
distribute-rgt-neg-inN/A
metadata-evalN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6460.7
Simplified60.7%
Final simplification73.3%
(FPCore (a k m) :precision binary64 (if (<= m -1.52) (/ a (* k k)) (if (<= m 1.1) (/ a (fma k k 1.0)) (* a (* (* k k) 99.0)))))
double code(double a, double k, double m) {
double tmp;
if (m <= -1.52) {
tmp = a / (k * k);
} else if (m <= 1.1) {
tmp = a / fma(k, k, 1.0);
} else {
tmp = a * ((k * k) * 99.0);
}
return tmp;
}
function code(a, k, m) tmp = 0.0 if (m <= -1.52) tmp = Float64(a / Float64(k * k)); elseif (m <= 1.1) tmp = Float64(a / fma(k, k, 1.0)); else tmp = Float64(a * Float64(Float64(k * k) * 99.0)); end return tmp end
code[a_, k_, m_] := If[LessEqual[m, -1.52], N[(a / N[(k * k), $MachinePrecision]), $MachinePrecision], If[LessEqual[m, 1.1], N[(a / N[(k * k + 1.0), $MachinePrecision]), $MachinePrecision], N[(a * N[(N[(k * k), $MachinePrecision] * 99.0), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq -1.52:\\
\;\;\;\;\frac{a}{k \cdot k}\\
\mathbf{elif}\;m \leq 1.1:\\
\;\;\;\;\frac{a}{\mathsf{fma}\left(k, k, 1\right)}\\
\mathbf{else}:\\
\;\;\;\;a \cdot \left(\left(k \cdot k\right) \cdot 99\right)\\
\end{array}
\end{array}
if m < -1.52Initial program 100.0%
Taylor expanded in m around 0
/-lowering-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
+-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
+-lowering-+.f6437.8
Simplified37.8%
Taylor expanded in k around inf
unpow2N/A
*-lowering-*.f6462.6
Simplified62.6%
if -1.52 < m < 1.1000000000000001Initial program 93.6%
Taylor expanded in m around 0
/-lowering-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
+-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
+-lowering-+.f6490.7
Simplified90.7%
Taylor expanded in k around inf
Simplified89.3%
if 1.1000000000000001 < m Initial program 70.0%
Taylor expanded in m around 0
/-lowering-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
+-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
+-lowering-+.f642.6
Simplified2.6%
Taylor expanded in k around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
sub-negN/A
mul-1-negN/A
distribute-neg-outN/A
neg-sub0N/A
--lowering--.f64N/A
distribute-rgt1-inN/A
associate-*r*N/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
accelerator-lowering-fma.f64N/A
metadata-eval29.4
Simplified29.4%
Taylor expanded in k around inf
*-commutativeN/A
associate-*l*N/A
metadata-evalN/A
distribute-rgt-neg-inN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
distribute-rgt-neg-inN/A
metadata-evalN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6460.7
Simplified60.7%
(FPCore (a k m) :precision binary64 (if (<= m -5.8e-168) (/ a (* k k)) (if (<= m 1.05) (/ a (fma k 10.0 1.0)) (* a (* (* k k) 99.0)))))
double code(double a, double k, double m) {
double tmp;
if (m <= -5.8e-168) {
tmp = a / (k * k);
} else if (m <= 1.05) {
tmp = a / fma(k, 10.0, 1.0);
} else {
tmp = a * ((k * k) * 99.0);
}
return tmp;
}
function code(a, k, m) tmp = 0.0 if (m <= -5.8e-168) tmp = Float64(a / Float64(k * k)); elseif (m <= 1.05) tmp = Float64(a / fma(k, 10.0, 1.0)); else tmp = Float64(a * Float64(Float64(k * k) * 99.0)); end return tmp end
code[a_, k_, m_] := If[LessEqual[m, -5.8e-168], N[(a / N[(k * k), $MachinePrecision]), $MachinePrecision], If[LessEqual[m, 1.05], N[(a / N[(k * 10.0 + 1.0), $MachinePrecision]), $MachinePrecision], N[(a * N[(N[(k * k), $MachinePrecision] * 99.0), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq -5.8 \cdot 10^{-168}:\\
\;\;\;\;\frac{a}{k \cdot k}\\
\mathbf{elif}\;m \leq 1.05:\\
\;\;\;\;\frac{a}{\mathsf{fma}\left(k, 10, 1\right)}\\
\mathbf{else}:\\
\;\;\;\;a \cdot \left(\left(k \cdot k\right) \cdot 99\right)\\
\end{array}
\end{array}
if m < -5.7999999999999997e-168Initial program 99.0%
Taylor expanded in m around 0
/-lowering-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
+-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
+-lowering-+.f6453.4
Simplified53.4%
Taylor expanded in k around inf
unpow2N/A
*-lowering-*.f6462.2
Simplified62.2%
if -5.7999999999999997e-168 < m < 1.05000000000000004Initial program 92.3%
Taylor expanded in m around 0
/-lowering-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
+-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
+-lowering-+.f6489.9
Simplified89.9%
Taylor expanded in k around 0
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6464.9
Simplified64.9%
if 1.05000000000000004 < m Initial program 70.0%
Taylor expanded in m around 0
/-lowering-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
+-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
+-lowering-+.f642.6
Simplified2.6%
Taylor expanded in k around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
sub-negN/A
mul-1-negN/A
distribute-neg-outN/A
neg-sub0N/A
--lowering--.f64N/A
distribute-rgt1-inN/A
associate-*r*N/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
accelerator-lowering-fma.f64N/A
metadata-eval29.4
Simplified29.4%
Taylor expanded in k around inf
*-commutativeN/A
associate-*l*N/A
metadata-evalN/A
distribute-rgt-neg-inN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
distribute-rgt-neg-inN/A
metadata-evalN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6460.7
Simplified60.7%
(FPCore (a k m) :precision binary64 (if (<= m 1.7) a (* a (* (* k k) 99.0))))
double code(double a, double k, double m) {
double tmp;
if (m <= 1.7) {
tmp = a;
} else {
tmp = a * ((k * k) * 99.0);
}
return tmp;
}
real(8) function code(a, k, m)
real(8), intent (in) :: a
real(8), intent (in) :: k
real(8), intent (in) :: m
real(8) :: tmp
if (m <= 1.7d0) then
tmp = a
else
tmp = a * ((k * k) * 99.0d0)
end if
code = tmp
end function
public static double code(double a, double k, double m) {
double tmp;
if (m <= 1.7) {
tmp = a;
} else {
tmp = a * ((k * k) * 99.0);
}
return tmp;
}
def code(a, k, m): tmp = 0 if m <= 1.7: tmp = a else: tmp = a * ((k * k) * 99.0) return tmp
function code(a, k, m) tmp = 0.0 if (m <= 1.7) tmp = a; else tmp = Float64(a * Float64(Float64(k * k) * 99.0)); end return tmp end
function tmp_2 = code(a, k, m) tmp = 0.0; if (m <= 1.7) tmp = a; else tmp = a * ((k * k) * 99.0); end tmp_2 = tmp; end
code[a_, k_, m_] := If[LessEqual[m, 1.7], a, N[(a * N[(N[(k * k), $MachinePrecision] * 99.0), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq 1.7:\\
\;\;\;\;a\\
\mathbf{else}:\\
\;\;\;\;a \cdot \left(\left(k \cdot k\right) \cdot 99\right)\\
\end{array}
\end{array}
if m < 1.69999999999999996Initial program 96.2%
Taylor expanded in m around 0
/-lowering-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
+-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
+-lowering-+.f6468.8
Simplified68.8%
Taylor expanded in k around 0
Simplified27.7%
if 1.69999999999999996 < m Initial program 70.0%
Taylor expanded in m around 0
/-lowering-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
+-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
+-lowering-+.f642.6
Simplified2.6%
Taylor expanded in k around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
sub-negN/A
mul-1-negN/A
distribute-neg-outN/A
neg-sub0N/A
--lowering--.f64N/A
distribute-rgt1-inN/A
associate-*r*N/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
accelerator-lowering-fma.f64N/A
metadata-eval29.4
Simplified29.4%
Taylor expanded in k around inf
*-commutativeN/A
associate-*l*N/A
metadata-evalN/A
distribute-rgt-neg-inN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
distribute-rgt-neg-inN/A
metadata-evalN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6460.7
Simplified60.7%
(FPCore (a k m) :precision binary64 (if (<= m 7200.0) a (* a (* (* k k) -0.0001))))
double code(double a, double k, double m) {
double tmp;
if (m <= 7200.0) {
tmp = a;
} else {
tmp = a * ((k * k) * -0.0001);
}
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 <= 7200.0d0) then
tmp = a
else
tmp = a * ((k * k) * (-0.0001d0))
end if
code = tmp
end function
public static double code(double a, double k, double m) {
double tmp;
if (m <= 7200.0) {
tmp = a;
} else {
tmp = a * ((k * k) * -0.0001);
}
return tmp;
}
def code(a, k, m): tmp = 0 if m <= 7200.0: tmp = a else: tmp = a * ((k * k) * -0.0001) return tmp
function code(a, k, m) tmp = 0.0 if (m <= 7200.0) tmp = a; else tmp = Float64(a * Float64(Float64(k * k) * -0.0001)); end return tmp end
function tmp_2 = code(a, k, m) tmp = 0.0; if (m <= 7200.0) tmp = a; else tmp = a * ((k * k) * -0.0001); end tmp_2 = tmp; end
code[a_, k_, m_] := If[LessEqual[m, 7200.0], a, N[(a * N[(N[(k * k), $MachinePrecision] * -0.0001), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq 7200:\\
\;\;\;\;a\\
\mathbf{else}:\\
\;\;\;\;a \cdot \left(\left(k \cdot k\right) \cdot -0.0001\right)\\
\end{array}
\end{array}
if m < 7200Initial program 95.7%
Taylor expanded in m around 0
/-lowering-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
+-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
+-lowering-+.f6468.4
Simplified68.4%
Taylor expanded in k around 0
Simplified27.6%
if 7200 < m Initial program 70.9%
Taylor expanded in m around 0
/-lowering-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
+-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
+-lowering-+.f642.6
Simplified2.6%
*-commutativeN/A
flip-+N/A
associate-*l/N/A
div-invN/A
accelerator-lowering-fma.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
metadata-evalN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
--lowering--.f642.5
Applied egg-rr2.5%
Taylor expanded in k around inf
/-lowering-/.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f6430.8
Simplified30.8%
Taylor expanded in k around 0
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6427.4
Simplified27.4%
(FPCore (a k m) :precision binary64 (if (<= m 4.8e+42) a (* k (* a -10.0))))
double code(double a, double k, double m) {
double tmp;
if (m <= 4.8e+42) {
tmp = a;
} else {
tmp = k * (a * -10.0);
}
return tmp;
}
real(8) function code(a, k, m)
real(8), intent (in) :: a
real(8), intent (in) :: k
real(8), intent (in) :: m
real(8) :: tmp
if (m <= 4.8d+42) then
tmp = a
else
tmp = k * (a * (-10.0d0))
end if
code = tmp
end function
public static double code(double a, double k, double m) {
double tmp;
if (m <= 4.8e+42) {
tmp = a;
} else {
tmp = k * (a * -10.0);
}
return tmp;
}
def code(a, k, m): tmp = 0 if m <= 4.8e+42: tmp = a else: tmp = k * (a * -10.0) return tmp
function code(a, k, m) tmp = 0.0 if (m <= 4.8e+42) tmp = a; else tmp = Float64(k * Float64(a * -10.0)); end return tmp end
function tmp_2 = code(a, k, m) tmp = 0.0; if (m <= 4.8e+42) tmp = a; else tmp = k * (a * -10.0); end tmp_2 = tmp; end
code[a_, k_, m_] := If[LessEqual[m, 4.8e+42], a, N[(k * N[(a * -10.0), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq 4.8 \cdot 10^{+42}:\\
\;\;\;\;a\\
\mathbf{else}:\\
\;\;\;\;k \cdot \left(a \cdot -10\right)\\
\end{array}
\end{array}
if m < 4.7999999999999997e42Initial program 93.3%
Taylor expanded in m around 0
/-lowering-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
+-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
+-lowering-+.f6463.8
Simplified63.8%
Taylor expanded in k around 0
Simplified25.9%
if 4.7999999999999997e42 < m Initial program 72.7%
Taylor expanded in m around 0
/-lowering-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
+-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
+-lowering-+.f642.7
Simplified2.7%
Taylor expanded in k around 0
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
distribute-rgt1-inN/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
accelerator-lowering-fma.f6410.2
Simplified10.2%
Taylor expanded in k around inf
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6419.8
Simplified19.8%
(FPCore (a k m) :precision binary64 a)
double code(double a, double k, double m) {
return a;
}
real(8) function code(a, k, m)
real(8), intent (in) :: a
real(8), intent (in) :: k
real(8), intent (in) :: m
code = a
end function
public static double code(double a, double k, double m) {
return a;
}
def code(a, k, m): return a
function code(a, k, m) return a end
function tmp = code(a, k, m) tmp = a; end
code[a_, k_, m_] := a
\begin{array}{l}
\\
a
\end{array}
Initial program 88.0%
Taylor expanded in m around 0
/-lowering-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
+-commutativeN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
+-lowering-+.f6448.1
Simplified48.1%
Taylor expanded in k around 0
Simplified20.2%
herbie shell --seed 2024196
(FPCore (a k m)
:name "Falkner and Boettcher, Appendix A"
:precision binary64
(/ (* a (pow k m)) (+ (+ 1.0 (* 10.0 k)) (* k k))))