
(FPCore (a k m) :precision binary64 (/ (* a (pow k m)) (+ (+ 1.0 (* 10.0 k)) (* k k))))
double code(double a, double k, double m) {
return (a * pow(k, m)) / ((1.0 + (10.0 * k)) + (k * k));
}
real(8) function code(a, k, m)
real(8), intent (in) :: a
real(8), intent (in) :: k
real(8), intent (in) :: m
code = (a * (k ** m)) / ((1.0d0 + (10.0d0 * k)) + (k * k))
end function
public static double code(double a, double k, double m) {
return (a * Math.pow(k, m)) / ((1.0 + (10.0 * k)) + (k * k));
}
def code(a, k, m): return (a * math.pow(k, m)) / ((1.0 + (10.0 * k)) + (k * k))
function code(a, k, m) return Float64(Float64(a * (k ^ m)) / Float64(Float64(1.0 + Float64(10.0 * k)) + Float64(k * k))) end
function tmp = code(a, k, m) tmp = (a * (k ^ m)) / ((1.0 + (10.0 * k)) + (k * k)); end
code[a_, k_, m_] := N[(N[(a * N[Power[k, m], $MachinePrecision]), $MachinePrecision] / N[(N[(1.0 + N[(10.0 * k), $MachinePrecision]), $MachinePrecision] + N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{a \cdot {k}^{m}}{\left(1 + 10 \cdot k\right) + k \cdot k}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 14 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (a k m) :precision binary64 (/ (* a (pow k m)) (+ (+ 1.0 (* 10.0 k)) (* k k))))
double code(double a, double k, double m) {
return (a * pow(k, m)) / ((1.0 + (10.0 * k)) + (k * k));
}
real(8) function code(a, k, m)
real(8), intent (in) :: a
real(8), intent (in) :: k
real(8), intent (in) :: m
code = (a * (k ** m)) / ((1.0d0 + (10.0d0 * k)) + (k * k))
end function
public static double code(double a, double k, double m) {
return (a * Math.pow(k, m)) / ((1.0 + (10.0 * k)) + (k * k));
}
def code(a, k, m): return (a * math.pow(k, m)) / ((1.0 + (10.0 * k)) + (k * k))
function code(a, k, m) return Float64(Float64(a * (k ^ m)) / Float64(Float64(1.0 + Float64(10.0 * k)) + Float64(k * k))) end
function tmp = code(a, k, m) tmp = (a * (k ^ m)) / ((1.0 + (10.0 * k)) + (k * k)); end
code[a_, k_, m_] := N[(N[(a * N[Power[k, m], $MachinePrecision]), $MachinePrecision] / N[(N[(1.0 + N[(10.0 * k), $MachinePrecision]), $MachinePrecision] + N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{a \cdot {k}^{m}}{\left(1 + 10 \cdot k\right) + k \cdot k}
\end{array}
(FPCore (a k m)
:precision binary64
(let* ((t_0 (* a (pow k m))))
(if (<= (/ t_0 (+ (+ 1.0 (* k 10.0)) (* k k))) 2e+270)
(* a (/ (pow k m) (+ 1.0 (* k (+ k 10.0)))))
t_0)))
double code(double a, double k, double m) {
double t_0 = a * pow(k, m);
double tmp;
if ((t_0 / ((1.0 + (k * 10.0)) + (k * k))) <= 2e+270) {
tmp = a * (pow(k, m) / (1.0 + (k * (k + 10.0))));
} else {
tmp = t_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 ** m)
if ((t_0 / ((1.0d0 + (k * 10.0d0)) + (k * k))) <= 2d+270) then
tmp = a * ((k ** m) / (1.0d0 + (k * (k + 10.0d0))))
else
tmp = t_0
end if
code = tmp
end function
public static double code(double a, double k, double m) {
double t_0 = a * Math.pow(k, m);
double tmp;
if ((t_0 / ((1.0 + (k * 10.0)) + (k * k))) <= 2e+270) {
tmp = a * (Math.pow(k, m) / (1.0 + (k * (k + 10.0))));
} else {
tmp = t_0;
}
return tmp;
}
def code(a, k, m): t_0 = a * math.pow(k, m) tmp = 0 if (t_0 / ((1.0 + (k * 10.0)) + (k * k))) <= 2e+270: tmp = a * (math.pow(k, m) / (1.0 + (k * (k + 10.0)))) else: tmp = t_0 return tmp
function code(a, k, m) t_0 = Float64(a * (k ^ m)) tmp = 0.0 if (Float64(t_0 / Float64(Float64(1.0 + Float64(k * 10.0)) + Float64(k * k))) <= 2e+270) tmp = Float64(a * Float64((k ^ m) / Float64(1.0 + Float64(k * Float64(k + 10.0))))); else tmp = t_0; end return tmp end
function tmp_2 = code(a, k, m) t_0 = a * (k ^ m); tmp = 0.0; if ((t_0 / ((1.0 + (k * 10.0)) + (k * k))) <= 2e+270) tmp = a * ((k ^ m) / (1.0 + (k * (k + 10.0)))); else tmp = t_0; end tmp_2 = tmp; end
code[a_, k_, m_] := Block[{t$95$0 = N[(a * N[Power[k, m], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[N[(t$95$0 / N[(N[(1.0 + N[(k * 10.0), $MachinePrecision]), $MachinePrecision] + N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2e+270], N[(a * N[(N[Power[k, m], $MachinePrecision] / N[(1.0 + N[(k * N[(k + 10.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$0]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := a \cdot {k}^{m}\\
\mathbf{if}\;\frac{t\_0}{\left(1 + k \cdot 10\right) + k \cdot k} \leq 2 \cdot 10^{+270}:\\
\;\;\;\;a \cdot \frac{{k}^{m}}{1 + k \cdot \left(k + 10\right)}\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if (/.f64 (*.f64 a (pow.f64 k m)) (+.f64 (+.f64 #s(literal 1 binary64) (*.f64 #s(literal 10 binary64) k)) (*.f64 k k))) < 2.0000000000000001e270Initial program 98.0%
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
pow-lowering-pow.f64N/A
associate-+l+N/A
+-lowering-+.f64N/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
+-lowering-+.f6498.0%
Applied egg-rr98.0%
if 2.0000000000000001e270 < (/.f64 (*.f64 a (pow.f64 k m)) (+.f64 (+.f64 #s(literal 1 binary64) (*.f64 #s(literal 10 binary64) k)) (*.f64 k k))) Initial program 60.3%
Taylor expanded in k around 0
*-lowering-*.f64N/A
pow-lowering-pow.f64100.0%
Simplified100.0%
Final simplification98.5%
(FPCore (a k m) :precision binary64 (if (<= m -0.00039) (* a (/ (pow k m) (* k k))) (if (<= m 2.2e-15) (/ a (+ 1.0 (* k (+ k 10.0)))) (* a (pow k m)))))
double code(double a, double k, double m) {
double tmp;
if (m <= -0.00039) {
tmp = a * (pow(k, m) / (k * k));
} else if (m <= 2.2e-15) {
tmp = a / (1.0 + (k * (k + 10.0)));
} else {
tmp = a * pow(k, m);
}
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.00039d0)) then
tmp = a * ((k ** m) / (k * k))
else if (m <= 2.2d-15) then
tmp = a / (1.0d0 + (k * (k + 10.0d0)))
else
tmp = a * (k ** m)
end if
code = tmp
end function
public static double code(double a, double k, double m) {
double tmp;
if (m <= -0.00039) {
tmp = a * (Math.pow(k, m) / (k * k));
} else if (m <= 2.2e-15) {
tmp = a / (1.0 + (k * (k + 10.0)));
} else {
tmp = a * Math.pow(k, m);
}
return tmp;
}
def code(a, k, m): tmp = 0 if m <= -0.00039: tmp = a * (math.pow(k, m) / (k * k)) elif m <= 2.2e-15: tmp = a / (1.0 + (k * (k + 10.0))) else: tmp = a * math.pow(k, m) return tmp
function code(a, k, m) tmp = 0.0 if (m <= -0.00039) tmp = Float64(a * Float64((k ^ m) / Float64(k * k))); elseif (m <= 2.2e-15) tmp = Float64(a / Float64(1.0 + Float64(k * Float64(k + 10.0)))); else tmp = Float64(a * (k ^ m)); end return tmp end
function tmp_2 = code(a, k, m) tmp = 0.0; if (m <= -0.00039) tmp = a * ((k ^ m) / (k * k)); elseif (m <= 2.2e-15) tmp = a / (1.0 + (k * (k + 10.0))); else tmp = a * (k ^ m); end tmp_2 = tmp; end
code[a_, k_, m_] := If[LessEqual[m, -0.00039], N[(a * N[(N[Power[k, m], $MachinePrecision] / N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[m, 2.2e-15], N[(a / N[(1.0 + N[(k * N[(k + 10.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(a * N[Power[k, m], $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq -0.00039:\\
\;\;\;\;a \cdot \frac{{k}^{m}}{k \cdot k}\\
\mathbf{elif}\;m \leq 2.2 \cdot 10^{-15}:\\
\;\;\;\;\frac{a}{1 + k \cdot \left(k + 10\right)}\\
\mathbf{else}:\\
\;\;\;\;a \cdot {k}^{m}\\
\end{array}
\end{array}
if m < -3.89999999999999993e-4Initial program 100.0%
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
pow-lowering-pow.f64N/A
associate-+l+N/A
+-lowering-+.f64N/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
+-lowering-+.f64100.0%
Applied egg-rr100.0%
Taylor expanded in k around inf
/-lowering-/.f64N/A
*-commutativeN/A
associate-*r*N/A
exp-prodN/A
neg-mul-1N/A
log-recN/A
remove-double-negN/A
exp-prodN/A
exp-to-powN/A
pow-lowering-pow.f64N/A
unpow2N/A
*-lowering-*.f64100.0%
Simplified100.0%
if -3.89999999999999993e-4 < m < 2.19999999999999986e-15Initial program 94.5%
Taylor expanded in m around 0
/-lowering-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
associate-*l*N/A
unpow2N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
Simplified92.5%
if 2.19999999999999986e-15 < m Initial program 75.3%
Taylor expanded in k around 0
*-lowering-*.f64N/A
pow-lowering-pow.f6499.0%
Simplified99.0%
Final simplification97.5%
(FPCore (a k m)
:precision binary64
(let* ((t_0 (* a (pow k m))))
(if (<= m -3.2e-10)
t_0
(if (<= m 2.2e-15) (/ a (+ 1.0 (* k (+ k 10.0)))) t_0))))
double code(double a, double k, double m) {
double t_0 = a * pow(k, m);
double tmp;
if (m <= -3.2e-10) {
tmp = t_0;
} else if (m <= 2.2e-15) {
tmp = a / (1.0 + (k * (k + 10.0)));
} else {
tmp = t_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 ** m)
if (m <= (-3.2d-10)) then
tmp = t_0
else if (m <= 2.2d-15) then
tmp = a / (1.0d0 + (k * (k + 10.0d0)))
else
tmp = t_0
end if
code = tmp
end function
public static double code(double a, double k, double m) {
double t_0 = a * Math.pow(k, m);
double tmp;
if (m <= -3.2e-10) {
tmp = t_0;
} else if (m <= 2.2e-15) {
tmp = a / (1.0 + (k * (k + 10.0)));
} else {
tmp = t_0;
}
return tmp;
}
def code(a, k, m): t_0 = a * math.pow(k, m) tmp = 0 if m <= -3.2e-10: tmp = t_0 elif m <= 2.2e-15: tmp = a / (1.0 + (k * (k + 10.0))) else: tmp = t_0 return tmp
function code(a, k, m) t_0 = Float64(a * (k ^ m)) tmp = 0.0 if (m <= -3.2e-10) tmp = t_0; elseif (m <= 2.2e-15) tmp = Float64(a / Float64(1.0 + Float64(k * Float64(k + 10.0)))); else tmp = t_0; end return tmp end
function tmp_2 = code(a, k, m) t_0 = a * (k ^ m); tmp = 0.0; if (m <= -3.2e-10) tmp = t_0; elseif (m <= 2.2e-15) tmp = a / (1.0 + (k * (k + 10.0))); else tmp = t_0; end tmp_2 = tmp; end
code[a_, k_, m_] := Block[{t$95$0 = N[(a * N[Power[k, m], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[m, -3.2e-10], t$95$0, If[LessEqual[m, 2.2e-15], N[(a / N[(1.0 + N[(k * N[(k + 10.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$0]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := a \cdot {k}^{m}\\
\mathbf{if}\;m \leq -3.2 \cdot 10^{-10}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;m \leq 2.2 \cdot 10^{-15}:\\
\;\;\;\;\frac{a}{1 + k \cdot \left(k + 10\right)}\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if m < -3.19999999999999981e-10 or 2.19999999999999986e-15 < m Initial program 87.8%
Taylor expanded in k around 0
*-lowering-*.f64N/A
pow-lowering-pow.f6498.5%
Simplified98.5%
if -3.19999999999999981e-10 < m < 2.19999999999999986e-15Initial program 94.3%
Taylor expanded in m around 0
/-lowering-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
associate-*l*N/A
unpow2N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
Simplified93.7%
Final simplification97.2%
(FPCore (a k m)
:precision binary64
(if (<= m -112000000.0)
(/ (- a (* (/ a k) (+ 10.0 (/ -99.0 k)))) (* k k))
(if (<= m 8e+18)
(/ a (+ 1.0 (* k (+ k 10.0))))
(/ a (- -99.0 (/ (+ (+ 1000.0 (/ 10000.0 k)) (/ 100000.0 (* k k))) k))))))
double code(double a, double k, double m) {
double tmp;
if (m <= -112000000.0) {
tmp = (a - ((a / k) * (10.0 + (-99.0 / k)))) / (k * k);
} else if (m <= 8e+18) {
tmp = a / (1.0 + (k * (k + 10.0)));
} else {
tmp = a / (-99.0 - (((1000.0 + (10000.0 / k)) + (100000.0 / (k * 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 (m <= (-112000000.0d0)) then
tmp = (a - ((a / k) * (10.0d0 + ((-99.0d0) / k)))) / (k * k)
else if (m <= 8d+18) then
tmp = a / (1.0d0 + (k * (k + 10.0d0)))
else
tmp = a / ((-99.0d0) - (((1000.0d0 + (10000.0d0 / k)) + (100000.0d0 / (k * k))) / k))
end if
code = tmp
end function
public static double code(double a, double k, double m) {
double tmp;
if (m <= -112000000.0) {
tmp = (a - ((a / k) * (10.0 + (-99.0 / k)))) / (k * k);
} else if (m <= 8e+18) {
tmp = a / (1.0 + (k * (k + 10.0)));
} else {
tmp = a / (-99.0 - (((1000.0 + (10000.0 / k)) + (100000.0 / (k * k))) / k));
}
return tmp;
}
def code(a, k, m): tmp = 0 if m <= -112000000.0: tmp = (a - ((a / k) * (10.0 + (-99.0 / k)))) / (k * k) elif m <= 8e+18: tmp = a / (1.0 + (k * (k + 10.0))) else: tmp = a / (-99.0 - (((1000.0 + (10000.0 / k)) + (100000.0 / (k * k))) / k)) return tmp
function code(a, k, m) tmp = 0.0 if (m <= -112000000.0) tmp = Float64(Float64(a - Float64(Float64(a / k) * Float64(10.0 + Float64(-99.0 / k)))) / Float64(k * k)); elseif (m <= 8e+18) tmp = Float64(a / Float64(1.0 + Float64(k * Float64(k + 10.0)))); else tmp = Float64(a / Float64(-99.0 - Float64(Float64(Float64(1000.0 + Float64(10000.0 / k)) + Float64(100000.0 / Float64(k * k))) / k))); end return tmp end
function tmp_2 = code(a, k, m) tmp = 0.0; if (m <= -112000000.0) tmp = (a - ((a / k) * (10.0 + (-99.0 / k)))) / (k * k); elseif (m <= 8e+18) tmp = a / (1.0 + (k * (k + 10.0))); else tmp = a / (-99.0 - (((1000.0 + (10000.0 / k)) + (100000.0 / (k * k))) / k)); end tmp_2 = tmp; end
code[a_, k_, m_] := If[LessEqual[m, -112000000.0], N[(N[(a - N[(N[(a / k), $MachinePrecision] * N[(10.0 + N[(-99.0 / k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(k * k), $MachinePrecision]), $MachinePrecision], If[LessEqual[m, 8e+18], N[(a / N[(1.0 + N[(k * N[(k + 10.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(a / N[(-99.0 - N[(N[(N[(1000.0 + N[(10000.0 / k), $MachinePrecision]), $MachinePrecision] + N[(100000.0 / N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq -112000000:\\
\;\;\;\;\frac{a - \frac{a}{k} \cdot \left(10 + \frac{-99}{k}\right)}{k \cdot k}\\
\mathbf{elif}\;m \leq 8 \cdot 10^{+18}:\\
\;\;\;\;\frac{a}{1 + k \cdot \left(k + 10\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{a}{-99 - \frac{\left(1000 + \frac{10000}{k}\right) + \frac{100000}{k \cdot k}}{k}}\\
\end{array}
\end{array}
if m < -1.12e8Initial program 100.0%
Taylor expanded in m around 0
/-lowering-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
associate-*l*N/A
unpow2N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
Simplified34.1%
flip3-+N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
clear-numN/A
flip3-+N/A
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f6434.1%
Applied egg-rr34.1%
Taylor expanded in k around inf
Simplified65.7%
if -1.12e8 < m < 8e18Initial program 93.8%
Taylor expanded in m around 0
/-lowering-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
associate-*l*N/A
unpow2N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
Simplified85.9%
if 8e18 < m Initial program 75.0%
Taylor expanded in m around 0
/-lowering-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
associate-*l*N/A
unpow2N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
Simplified2.9%
flip3-+N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
clear-numN/A
flip3-+N/A
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f642.9%
Applied egg-rr2.9%
Taylor expanded in k around 0
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f642.3%
Simplified2.3%
Taylor expanded in k around -inf
sub-negN/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
metadata-evalN/A
/-lowering-/.f64N/A
associate-+r+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6436.5%
Simplified36.5%
Final simplification61.9%
(FPCore (a k m)
:precision binary64
(if (<= m -112000000.0)
(/ (- a (* (/ a k) (+ 10.0 (/ -99.0 k)))) (* k k))
(if (<= m 8e+18)
(/ a (+ 1.0 (* k (+ k 10.0))))
(/ a (- -99.0 (/ (+ 1000.0 (/ 10000.0 k)) k))))))
double code(double a, double k, double m) {
double tmp;
if (m <= -112000000.0) {
tmp = (a - ((a / k) * (10.0 + (-99.0 / k)))) / (k * k);
} else if (m <= 8e+18) {
tmp = a / (1.0 + (k * (k + 10.0)));
} else {
tmp = a / (-99.0 - ((1000.0 + (10000.0 / 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 (m <= (-112000000.0d0)) then
tmp = (a - ((a / k) * (10.0d0 + ((-99.0d0) / k)))) / (k * k)
else if (m <= 8d+18) then
tmp = a / (1.0d0 + (k * (k + 10.0d0)))
else
tmp = a / ((-99.0d0) - ((1000.0d0 + (10000.0d0 / k)) / k))
end if
code = tmp
end function
public static double code(double a, double k, double m) {
double tmp;
if (m <= -112000000.0) {
tmp = (a - ((a / k) * (10.0 + (-99.0 / k)))) / (k * k);
} else if (m <= 8e+18) {
tmp = a / (1.0 + (k * (k + 10.0)));
} else {
tmp = a / (-99.0 - ((1000.0 + (10000.0 / k)) / k));
}
return tmp;
}
def code(a, k, m): tmp = 0 if m <= -112000000.0: tmp = (a - ((a / k) * (10.0 + (-99.0 / k)))) / (k * k) elif m <= 8e+18: tmp = a / (1.0 + (k * (k + 10.0))) else: tmp = a / (-99.0 - ((1000.0 + (10000.0 / k)) / k)) return tmp
function code(a, k, m) tmp = 0.0 if (m <= -112000000.0) tmp = Float64(Float64(a - Float64(Float64(a / k) * Float64(10.0 + Float64(-99.0 / k)))) / Float64(k * k)); elseif (m <= 8e+18) tmp = Float64(a / Float64(1.0 + Float64(k * Float64(k + 10.0)))); else tmp = Float64(a / Float64(-99.0 - Float64(Float64(1000.0 + Float64(10000.0 / k)) / k))); end return tmp end
function tmp_2 = code(a, k, m) tmp = 0.0; if (m <= -112000000.0) tmp = (a - ((a / k) * (10.0 + (-99.0 / k)))) / (k * k); elseif (m <= 8e+18) tmp = a / (1.0 + (k * (k + 10.0))); else tmp = a / (-99.0 - ((1000.0 + (10000.0 / k)) / k)); end tmp_2 = tmp; end
code[a_, k_, m_] := If[LessEqual[m, -112000000.0], N[(N[(a - N[(N[(a / k), $MachinePrecision] * N[(10.0 + N[(-99.0 / k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(k * k), $MachinePrecision]), $MachinePrecision], If[LessEqual[m, 8e+18], N[(a / N[(1.0 + N[(k * N[(k + 10.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(a / N[(-99.0 - N[(N[(1000.0 + N[(10000.0 / k), $MachinePrecision]), $MachinePrecision] / k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq -112000000:\\
\;\;\;\;\frac{a - \frac{a}{k} \cdot \left(10 + \frac{-99}{k}\right)}{k \cdot k}\\
\mathbf{elif}\;m \leq 8 \cdot 10^{+18}:\\
\;\;\;\;\frac{a}{1 + k \cdot \left(k + 10\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{a}{-99 - \frac{1000 + \frac{10000}{k}}{k}}\\
\end{array}
\end{array}
if m < -1.12e8Initial program 100.0%
Taylor expanded in m around 0
/-lowering-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
associate-*l*N/A
unpow2N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
Simplified34.1%
flip3-+N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
clear-numN/A
flip3-+N/A
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f6434.1%
Applied egg-rr34.1%
Taylor expanded in k around inf
Simplified65.7%
if -1.12e8 < m < 8e18Initial program 93.8%
Taylor expanded in m around 0
/-lowering-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
associate-*l*N/A
unpow2N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
Simplified85.9%
if 8e18 < m Initial program 75.0%
Taylor expanded in m around 0
/-lowering-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
associate-*l*N/A
unpow2N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
Simplified2.9%
flip3-+N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
clear-numN/A
flip3-+N/A
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f642.9%
Applied egg-rr2.9%
Taylor expanded in k around 0
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f642.3%
Simplified2.3%
Taylor expanded in k around inf
sub-negN/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
metadata-evalN/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f6429.3%
Simplified29.3%
Final simplification59.4%
(FPCore (a k m)
:precision binary64
(if (<= m -112000000.0)
(* a (/ 1.0 (* k k)))
(if (<= m 8e+18)
(/ a (+ 1.0 (* k (+ k 10.0))))
(/ a (- -99.0 (/ (+ 1000.0 (/ 10000.0 k)) k))))))
double code(double a, double k, double m) {
double tmp;
if (m <= -112000000.0) {
tmp = a * (1.0 / (k * k));
} else if (m <= 8e+18) {
tmp = a / (1.0 + (k * (k + 10.0)));
} else {
tmp = a / (-99.0 - ((1000.0 + (10000.0 / 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 (m <= (-112000000.0d0)) then
tmp = a * (1.0d0 / (k * k))
else if (m <= 8d+18) then
tmp = a / (1.0d0 + (k * (k + 10.0d0)))
else
tmp = a / ((-99.0d0) - ((1000.0d0 + (10000.0d0 / k)) / k))
end if
code = tmp
end function
public static double code(double a, double k, double m) {
double tmp;
if (m <= -112000000.0) {
tmp = a * (1.0 / (k * k));
} else if (m <= 8e+18) {
tmp = a / (1.0 + (k * (k + 10.0)));
} else {
tmp = a / (-99.0 - ((1000.0 + (10000.0 / k)) / k));
}
return tmp;
}
def code(a, k, m): tmp = 0 if m <= -112000000.0: tmp = a * (1.0 / (k * k)) elif m <= 8e+18: tmp = a / (1.0 + (k * (k + 10.0))) else: tmp = a / (-99.0 - ((1000.0 + (10000.0 / k)) / k)) return tmp
function code(a, k, m) tmp = 0.0 if (m <= -112000000.0) tmp = Float64(a * Float64(1.0 / Float64(k * k))); elseif (m <= 8e+18) tmp = Float64(a / Float64(1.0 + Float64(k * Float64(k + 10.0)))); else tmp = Float64(a / Float64(-99.0 - Float64(Float64(1000.0 + Float64(10000.0 / k)) / k))); end return tmp end
function tmp_2 = code(a, k, m) tmp = 0.0; if (m <= -112000000.0) tmp = a * (1.0 / (k * k)); elseif (m <= 8e+18) tmp = a / (1.0 + (k * (k + 10.0))); else tmp = a / (-99.0 - ((1000.0 + (10000.0 / k)) / k)); end tmp_2 = tmp; end
code[a_, k_, m_] := If[LessEqual[m, -112000000.0], N[(a * N[(1.0 / N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[m, 8e+18], N[(a / N[(1.0 + N[(k * N[(k + 10.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(a / N[(-99.0 - N[(N[(1000.0 + N[(10000.0 / k), $MachinePrecision]), $MachinePrecision] / k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq -112000000:\\
\;\;\;\;a \cdot \frac{1}{k \cdot k}\\
\mathbf{elif}\;m \leq 8 \cdot 10^{+18}:\\
\;\;\;\;\frac{a}{1 + k \cdot \left(k + 10\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{a}{-99 - \frac{1000 + \frac{10000}{k}}{k}}\\
\end{array}
\end{array}
if m < -1.12e8Initial program 100.0%
Taylor expanded in m around 0
/-lowering-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
associate-*l*N/A
unpow2N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
Simplified34.1%
Taylor expanded in k around inf
unpow2N/A
*-lowering-*.f6459.1%
Simplified59.1%
clear-numN/A
associate-/r/N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6461.2%
Applied egg-rr61.2%
if -1.12e8 < m < 8e18Initial program 93.8%
Taylor expanded in m around 0
/-lowering-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
associate-*l*N/A
unpow2N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
Simplified85.9%
if 8e18 < m Initial program 75.0%
Taylor expanded in m around 0
/-lowering-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
associate-*l*N/A
unpow2N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
Simplified2.9%
flip3-+N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
clear-numN/A
flip3-+N/A
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f642.9%
Applied egg-rr2.9%
Taylor expanded in k around 0
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f642.3%
Simplified2.3%
Taylor expanded in k around inf
sub-negN/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
metadata-evalN/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f6429.3%
Simplified29.3%
Final simplification57.9%
(FPCore (a k m)
:precision binary64
(if (<= m -112000000.0)
(* a (/ 1.0 (* k k)))
(if (<= m 1.2e+20)
(/ a (+ 1.0 (* k (+ k 10.0))))
(/ a (+ -99.0 (/ -1000.0 k))))))
double code(double a, double k, double m) {
double tmp;
if (m <= -112000000.0) {
tmp = a * (1.0 / (k * k));
} else if (m <= 1.2e+20) {
tmp = a / (1.0 + (k * (k + 10.0)));
} else {
tmp = a / (-99.0 + (-1000.0 / k));
}
return tmp;
}
real(8) function code(a, k, m)
real(8), intent (in) :: a
real(8), intent (in) :: k
real(8), intent (in) :: m
real(8) :: tmp
if (m <= (-112000000.0d0)) then
tmp = a * (1.0d0 / (k * k))
else if (m <= 1.2d+20) then
tmp = a / (1.0d0 + (k * (k + 10.0d0)))
else
tmp = a / ((-99.0d0) + ((-1000.0d0) / k))
end if
code = tmp
end function
public static double code(double a, double k, double m) {
double tmp;
if (m <= -112000000.0) {
tmp = a * (1.0 / (k * k));
} else if (m <= 1.2e+20) {
tmp = a / (1.0 + (k * (k + 10.0)));
} else {
tmp = a / (-99.0 + (-1000.0 / k));
}
return tmp;
}
def code(a, k, m): tmp = 0 if m <= -112000000.0: tmp = a * (1.0 / (k * k)) elif m <= 1.2e+20: tmp = a / (1.0 + (k * (k + 10.0))) else: tmp = a / (-99.0 + (-1000.0 / k)) return tmp
function code(a, k, m) tmp = 0.0 if (m <= -112000000.0) tmp = Float64(a * Float64(1.0 / Float64(k * k))); elseif (m <= 1.2e+20) tmp = Float64(a / Float64(1.0 + Float64(k * Float64(k + 10.0)))); else tmp = Float64(a / Float64(-99.0 + Float64(-1000.0 / k))); end return tmp end
function tmp_2 = code(a, k, m) tmp = 0.0; if (m <= -112000000.0) tmp = a * (1.0 / (k * k)); elseif (m <= 1.2e+20) tmp = a / (1.0 + (k * (k + 10.0))); else tmp = a / (-99.0 + (-1000.0 / k)); end tmp_2 = tmp; end
code[a_, k_, m_] := If[LessEqual[m, -112000000.0], N[(a * N[(1.0 / N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[m, 1.2e+20], N[(a / N[(1.0 + N[(k * N[(k + 10.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(a / N[(-99.0 + N[(-1000.0 / k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq -112000000:\\
\;\;\;\;a \cdot \frac{1}{k \cdot k}\\
\mathbf{elif}\;m \leq 1.2 \cdot 10^{+20}:\\
\;\;\;\;\frac{a}{1 + k \cdot \left(k + 10\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{a}{-99 + \frac{-1000}{k}}\\
\end{array}
\end{array}
if m < -1.12e8Initial program 100.0%
Taylor expanded in m around 0
/-lowering-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
associate-*l*N/A
unpow2N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
Simplified34.1%
Taylor expanded in k around inf
unpow2N/A
*-lowering-*.f6459.1%
Simplified59.1%
clear-numN/A
associate-/r/N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6461.2%
Applied egg-rr61.2%
if -1.12e8 < m < 1.2e20Initial program 93.8%
Taylor expanded in m around 0
/-lowering-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
associate-*l*N/A
unpow2N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
Simplified85.9%
if 1.2e20 < m Initial program 75.0%
Taylor expanded in m around 0
/-lowering-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
associate-*l*N/A
unpow2N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
Simplified2.9%
flip3-+N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
clear-numN/A
flip3-+N/A
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f642.9%
Applied egg-rr2.9%
Taylor expanded in k around 0
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f642.3%
Simplified2.3%
Taylor expanded in k around inf
mul-1-negN/A
distribute-neg-inN/A
+-lowering-+.f64N/A
metadata-evalN/A
associate-*r/N/A
metadata-evalN/A
distribute-neg-fracN/A
/-lowering-/.f64N/A
metadata-eval13.3%
Simplified13.3%
Final simplification52.4%
(FPCore (a k m) :precision binary64 (if (<= m -112000000.0) (* a (/ 1.0 (* k k))) (if (<= m 1.55e+19) (/ a (+ 1.0 (* k k))) (/ a (+ -99.0 (/ -1000.0 k))))))
double code(double a, double k, double m) {
double tmp;
if (m <= -112000000.0) {
tmp = a * (1.0 / (k * k));
} else if (m <= 1.55e+19) {
tmp = a / (1.0 + (k * k));
} else {
tmp = a / (-99.0 + (-1000.0 / k));
}
return tmp;
}
real(8) function code(a, k, m)
real(8), intent (in) :: a
real(8), intent (in) :: k
real(8), intent (in) :: m
real(8) :: tmp
if (m <= (-112000000.0d0)) then
tmp = a * (1.0d0 / (k * k))
else if (m <= 1.55d+19) then
tmp = a / (1.0d0 + (k * k))
else
tmp = a / ((-99.0d0) + ((-1000.0d0) / k))
end if
code = tmp
end function
public static double code(double a, double k, double m) {
double tmp;
if (m <= -112000000.0) {
tmp = a * (1.0 / (k * k));
} else if (m <= 1.55e+19) {
tmp = a / (1.0 + (k * k));
} else {
tmp = a / (-99.0 + (-1000.0 / k));
}
return tmp;
}
def code(a, k, m): tmp = 0 if m <= -112000000.0: tmp = a * (1.0 / (k * k)) elif m <= 1.55e+19: tmp = a / (1.0 + (k * k)) else: tmp = a / (-99.0 + (-1000.0 / k)) return tmp
function code(a, k, m) tmp = 0.0 if (m <= -112000000.0) tmp = Float64(a * Float64(1.0 / Float64(k * k))); elseif (m <= 1.55e+19) tmp = Float64(a / Float64(1.0 + Float64(k * k))); else tmp = Float64(a / Float64(-99.0 + Float64(-1000.0 / k))); end return tmp end
function tmp_2 = code(a, k, m) tmp = 0.0; if (m <= -112000000.0) tmp = a * (1.0 / (k * k)); elseif (m <= 1.55e+19) tmp = a / (1.0 + (k * k)); else tmp = a / (-99.0 + (-1000.0 / k)); end tmp_2 = tmp; end
code[a_, k_, m_] := If[LessEqual[m, -112000000.0], N[(a * N[(1.0 / N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[m, 1.55e+19], N[(a / N[(1.0 + N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(a / N[(-99.0 + N[(-1000.0 / k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq -112000000:\\
\;\;\;\;a \cdot \frac{1}{k \cdot k}\\
\mathbf{elif}\;m \leq 1.55 \cdot 10^{+19}:\\
\;\;\;\;\frac{a}{1 + k \cdot k}\\
\mathbf{else}:\\
\;\;\;\;\frac{a}{-99 + \frac{-1000}{k}}\\
\end{array}
\end{array}
if m < -1.12e8Initial program 100.0%
Taylor expanded in m around 0
/-lowering-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
associate-*l*N/A
unpow2N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
Simplified34.1%
Taylor expanded in k around inf
unpow2N/A
*-lowering-*.f6459.1%
Simplified59.1%
clear-numN/A
associate-/r/N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6461.2%
Applied egg-rr61.2%
if -1.12e8 < m < 1.55e19Initial program 93.8%
Taylor expanded in m around 0
/-lowering-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
associate-*l*N/A
unpow2N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
Simplified85.9%
Taylor expanded in k around inf
unpow2N/A
*-lowering-*.f6484.5%
Simplified84.5%
if 1.55e19 < m Initial program 75.0%
Taylor expanded in m around 0
/-lowering-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
associate-*l*N/A
unpow2N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
Simplified2.9%
flip3-+N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
clear-numN/A
flip3-+N/A
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f642.9%
Applied egg-rr2.9%
Taylor expanded in k around 0
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f642.3%
Simplified2.3%
Taylor expanded in k around inf
mul-1-negN/A
distribute-neg-inN/A
+-lowering-+.f64N/A
metadata-evalN/A
associate-*r/N/A
metadata-evalN/A
distribute-neg-fracN/A
/-lowering-/.f64N/A
metadata-eval13.3%
Simplified13.3%
Final simplification51.9%
(FPCore (a k m) :precision binary64 (if (<= m -4.2e-22) (* a (/ 1.0 (* k k))) (if (<= m 8e+18) (/ a (+ 1.0 (* k 10.0))) (/ a (+ -99.0 (/ -1000.0 k))))))
double code(double a, double k, double m) {
double tmp;
if (m <= -4.2e-22) {
tmp = a * (1.0 / (k * k));
} else if (m <= 8e+18) {
tmp = a / (1.0 + (k * 10.0));
} else {
tmp = a / (-99.0 + (-1000.0 / k));
}
return tmp;
}
real(8) function code(a, k, m)
real(8), intent (in) :: a
real(8), intent (in) :: k
real(8), intent (in) :: m
real(8) :: tmp
if (m <= (-4.2d-22)) then
tmp = a * (1.0d0 / (k * k))
else if (m <= 8d+18) then
tmp = a / (1.0d0 + (k * 10.0d0))
else
tmp = a / ((-99.0d0) + ((-1000.0d0) / k))
end if
code = tmp
end function
public static double code(double a, double k, double m) {
double tmp;
if (m <= -4.2e-22) {
tmp = a * (1.0 / (k * k));
} else if (m <= 8e+18) {
tmp = a / (1.0 + (k * 10.0));
} else {
tmp = a / (-99.0 + (-1000.0 / k));
}
return tmp;
}
def code(a, k, m): tmp = 0 if m <= -4.2e-22: tmp = a * (1.0 / (k * k)) elif m <= 8e+18: tmp = a / (1.0 + (k * 10.0)) else: tmp = a / (-99.0 + (-1000.0 / k)) return tmp
function code(a, k, m) tmp = 0.0 if (m <= -4.2e-22) tmp = Float64(a * Float64(1.0 / Float64(k * k))); elseif (m <= 8e+18) tmp = Float64(a / Float64(1.0 + Float64(k * 10.0))); else tmp = Float64(a / Float64(-99.0 + Float64(-1000.0 / k))); end return tmp end
function tmp_2 = code(a, k, m) tmp = 0.0; if (m <= -4.2e-22) tmp = a * (1.0 / (k * k)); elseif (m <= 8e+18) tmp = a / (1.0 + (k * 10.0)); else tmp = a / (-99.0 + (-1000.0 / k)); end tmp_2 = tmp; end
code[a_, k_, m_] := If[LessEqual[m, -4.2e-22], N[(a * N[(1.0 / N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[m, 8e+18], N[(a / N[(1.0 + N[(k * 10.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(a / N[(-99.0 + N[(-1000.0 / k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq -4.2 \cdot 10^{-22}:\\
\;\;\;\;a \cdot \frac{1}{k \cdot k}\\
\mathbf{elif}\;m \leq 8 \cdot 10^{+18}:\\
\;\;\;\;\frac{a}{1 + k \cdot 10}\\
\mathbf{else}:\\
\;\;\;\;\frac{a}{-99 + \frac{-1000}{k}}\\
\end{array}
\end{array}
if m < -4.20000000000000016e-22Initial program 100.0%
Taylor expanded in m around 0
/-lowering-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
associate-*l*N/A
unpow2N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
Simplified36.1%
Taylor expanded in k around inf
unpow2N/A
*-lowering-*.f6458.1%
Simplified58.1%
clear-numN/A
associate-/r/N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6460.1%
Applied egg-rr60.1%
if -4.20000000000000016e-22 < m < 8e18Initial program 93.2%
Taylor expanded in m around 0
/-lowering-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
associate-*l*N/A
unpow2N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
Simplified89.0%
Taylor expanded in k around 0
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6456.2%
Simplified56.2%
if 8e18 < m Initial program 75.0%
Taylor expanded in m around 0
/-lowering-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
associate-*l*N/A
unpow2N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
Simplified2.9%
flip3-+N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
clear-numN/A
flip3-+N/A
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f642.9%
Applied egg-rr2.9%
Taylor expanded in k around 0
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f642.3%
Simplified2.3%
Taylor expanded in k around inf
mul-1-negN/A
distribute-neg-inN/A
+-lowering-+.f64N/A
metadata-evalN/A
associate-*r/N/A
metadata-evalN/A
distribute-neg-fracN/A
/-lowering-/.f64N/A
metadata-eval13.3%
Simplified13.3%
Final simplification42.9%
(FPCore (a k m) :precision binary64 (if (<= k -4.7e-303) (* a (/ 1.0 (* k k))) (if (<= k 0.075) (* a (+ 1.0 (* k -10.0))) (/ a (* k (+ k 10.0))))))
double code(double a, double k, double m) {
double tmp;
if (k <= -4.7e-303) {
tmp = a * (1.0 / (k * k));
} else if (k <= 0.075) {
tmp = a * (1.0 + (k * -10.0));
} else {
tmp = a / (k * (k + 10.0));
}
return tmp;
}
real(8) function code(a, k, m)
real(8), intent (in) :: a
real(8), intent (in) :: k
real(8), intent (in) :: m
real(8) :: tmp
if (k <= (-4.7d-303)) then
tmp = a * (1.0d0 / (k * k))
else if (k <= 0.075d0) then
tmp = a * (1.0d0 + (k * (-10.0d0)))
else
tmp = a / (k * (k + 10.0d0))
end if
code = tmp
end function
public static double code(double a, double k, double m) {
double tmp;
if (k <= -4.7e-303) {
tmp = a * (1.0 / (k * k));
} else if (k <= 0.075) {
tmp = a * (1.0 + (k * -10.0));
} else {
tmp = a / (k * (k + 10.0));
}
return tmp;
}
def code(a, k, m): tmp = 0 if k <= -4.7e-303: tmp = a * (1.0 / (k * k)) elif k <= 0.075: tmp = a * (1.0 + (k * -10.0)) else: tmp = a / (k * (k + 10.0)) return tmp
function code(a, k, m) tmp = 0.0 if (k <= -4.7e-303) tmp = Float64(a * Float64(1.0 / Float64(k * k))); elseif (k <= 0.075) tmp = Float64(a * Float64(1.0 + Float64(k * -10.0))); else tmp = Float64(a / Float64(k * Float64(k + 10.0))); end return tmp end
function tmp_2 = code(a, k, m) tmp = 0.0; if (k <= -4.7e-303) tmp = a * (1.0 / (k * k)); elseif (k <= 0.075) tmp = a * (1.0 + (k * -10.0)); else tmp = a / (k * (k + 10.0)); end tmp_2 = tmp; end
code[a_, k_, m_] := If[LessEqual[k, -4.7e-303], N[(a * N[(1.0 / N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[k, 0.075], N[(a * N[(1.0 + N[(k * -10.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(a / N[(k * N[(k + 10.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;k \leq -4.7 \cdot 10^{-303}:\\
\;\;\;\;a \cdot \frac{1}{k \cdot k}\\
\mathbf{elif}\;k \leq 0.075:\\
\;\;\;\;a \cdot \left(1 + k \cdot -10\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{a}{k \cdot \left(k + 10\right)}\\
\end{array}
\end{array}
if k < -4.6999999999999997e-303Initial program 85.7%
Taylor expanded in m around 0
/-lowering-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
associate-*l*N/A
unpow2N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
Simplified20.6%
Taylor expanded in k around inf
unpow2N/A
*-lowering-*.f6428.2%
Simplified28.2%
clear-numN/A
associate-/r/N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6429.3%
Applied egg-rr29.3%
if -4.6999999999999997e-303 < k < 0.0749999999999999972Initial program 100.0%
Taylor expanded in m around 0
/-lowering-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
associate-*l*N/A
unpow2N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
Simplified42.3%
flip3-+N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
clear-numN/A
flip3-+N/A
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f6442.3%
Applied egg-rr42.3%
Taylor expanded in k around 0
*-commutativeN/A
associate-*r*N/A
distribute-rgt1-inN/A
+-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6442.3%
Simplified42.3%
if 0.0749999999999999972 < k Initial program 83.1%
Taylor expanded in m around 0
/-lowering-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
associate-*l*N/A
unpow2N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
Simplified54.5%
Taylor expanded in k around inf
unpow2N/A
associate-*l*N/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f6454.5%
Simplified54.5%
Final simplification42.2%
(FPCore (a k m) :precision binary64 (let* ((t_0 (* a (/ 1.0 (* k k))))) (if (<= k -4.7e-303) t_0 (if (<= k 0.1) (* a (+ 1.0 (* k -10.0))) t_0))))
double code(double a, double k, double m) {
double t_0 = a * (1.0 / (k * k));
double tmp;
if (k <= -4.7e-303) {
tmp = t_0;
} else if (k <= 0.1) {
tmp = a * (1.0 + (k * -10.0));
} else {
tmp = t_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 * (1.0d0 / (k * k))
if (k <= (-4.7d-303)) then
tmp = t_0
else if (k <= 0.1d0) then
tmp = a * (1.0d0 + (k * (-10.0d0)))
else
tmp = t_0
end if
code = tmp
end function
public static double code(double a, double k, double m) {
double t_0 = a * (1.0 / (k * k));
double tmp;
if (k <= -4.7e-303) {
tmp = t_0;
} else if (k <= 0.1) {
tmp = a * (1.0 + (k * -10.0));
} else {
tmp = t_0;
}
return tmp;
}
def code(a, k, m): t_0 = a * (1.0 / (k * k)) tmp = 0 if k <= -4.7e-303: tmp = t_0 elif k <= 0.1: tmp = a * (1.0 + (k * -10.0)) else: tmp = t_0 return tmp
function code(a, k, m) t_0 = Float64(a * Float64(1.0 / Float64(k * k))) tmp = 0.0 if (k <= -4.7e-303) tmp = t_0; elseif (k <= 0.1) tmp = Float64(a * Float64(1.0 + Float64(k * -10.0))); else tmp = t_0; end return tmp end
function tmp_2 = code(a, k, m) t_0 = a * (1.0 / (k * k)); tmp = 0.0; if (k <= -4.7e-303) tmp = t_0; elseif (k <= 0.1) tmp = a * (1.0 + (k * -10.0)); else tmp = t_0; end tmp_2 = tmp; end
code[a_, k_, m_] := Block[{t$95$0 = N[(a * N[(1.0 / N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[k, -4.7e-303], t$95$0, If[LessEqual[k, 0.1], N[(a * N[(1.0 + N[(k * -10.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$0]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := a \cdot \frac{1}{k \cdot k}\\
\mathbf{if}\;k \leq -4.7 \cdot 10^{-303}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;k \leq 0.1:\\
\;\;\;\;a \cdot \left(1 + k \cdot -10\right)\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if k < -4.6999999999999997e-303 or 0.10000000000000001 < k Initial program 84.4%
Taylor expanded in m around 0
/-lowering-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
associate-*l*N/A
unpow2N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
Simplified37.9%
Taylor expanded in k around inf
unpow2N/A
*-lowering-*.f6441.0%
Simplified41.0%
clear-numN/A
associate-/r/N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6441.6%
Applied egg-rr41.6%
if -4.6999999999999997e-303 < k < 0.10000000000000001Initial program 100.0%
Taylor expanded in m around 0
/-lowering-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
associate-*l*N/A
unpow2N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
Simplified42.3%
flip3-+N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
clear-numN/A
flip3-+N/A
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f6442.3%
Applied egg-rr42.3%
Taylor expanded in k around 0
*-commutativeN/A
associate-*r*N/A
distribute-rgt1-inN/A
+-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6442.3%
Simplified42.3%
Final simplification41.8%
(FPCore (a k m) :precision binary64 (let* ((t_0 (/ a (* k k)))) (if (<= k -4.7e-303) t_0 (if (<= k 0.1) (* a (+ 1.0 (* k -10.0))) t_0))))
double code(double a, double k, double m) {
double t_0 = a / (k * k);
double tmp;
if (k <= -4.7e-303) {
tmp = t_0;
} else if (k <= 0.1) {
tmp = a * (1.0 + (k * -10.0));
} else {
tmp = t_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 (k <= (-4.7d-303)) then
tmp = t_0
else if (k <= 0.1d0) then
tmp = a * (1.0d0 + (k * (-10.0d0)))
else
tmp = t_0
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 (k <= -4.7e-303) {
tmp = t_0;
} else if (k <= 0.1) {
tmp = a * (1.0 + (k * -10.0));
} else {
tmp = t_0;
}
return tmp;
}
def code(a, k, m): t_0 = a / (k * k) tmp = 0 if k <= -4.7e-303: tmp = t_0 elif k <= 0.1: tmp = a * (1.0 + (k * -10.0)) else: tmp = t_0 return tmp
function code(a, k, m) t_0 = Float64(a / Float64(k * k)) tmp = 0.0 if (k <= -4.7e-303) tmp = t_0; elseif (k <= 0.1) tmp = Float64(a * Float64(1.0 + Float64(k * -10.0))); else tmp = t_0; end return tmp end
function tmp_2 = code(a, k, m) t_0 = a / (k * k); tmp = 0.0; if (k <= -4.7e-303) tmp = t_0; elseif (k <= 0.1) tmp = a * (1.0 + (k * -10.0)); else tmp = t_0; end tmp_2 = tmp; end
code[a_, k_, m_] := Block[{t$95$0 = N[(a / N[(k * k), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[k, -4.7e-303], t$95$0, If[LessEqual[k, 0.1], N[(a * N[(1.0 + N[(k * -10.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$0]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{a}{k \cdot k}\\
\mathbf{if}\;k \leq -4.7 \cdot 10^{-303}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;k \leq 0.1:\\
\;\;\;\;a \cdot \left(1 + k \cdot -10\right)\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if k < -4.6999999999999997e-303 or 0.10000000000000001 < k Initial program 84.4%
Taylor expanded in m around 0
/-lowering-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
associate-*l*N/A
unpow2N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
Simplified37.9%
Taylor expanded in k around inf
unpow2N/A
*-lowering-*.f6441.0%
Simplified41.0%
if -4.6999999999999997e-303 < k < 0.10000000000000001Initial program 100.0%
Taylor expanded in m around 0
/-lowering-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
associate-*l*N/A
unpow2N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
Simplified42.3%
flip3-+N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
clear-numN/A
flip3-+N/A
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f6442.3%
Applied egg-rr42.3%
Taylor expanded in k around 0
*-commutativeN/A
associate-*r*N/A
distribute-rgt1-inN/A
+-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6442.3%
Simplified42.3%
Final simplification41.4%
(FPCore (a k m) :precision binary64 (let* ((t_0 (/ a (* k k)))) (if (<= k -4.7e-303) t_0 (if (<= k 235000.0) a t_0))))
double code(double a, double k, double m) {
double t_0 = a / (k * k);
double tmp;
if (k <= -4.7e-303) {
tmp = t_0;
} else if (k <= 235000.0) {
tmp = a;
} else {
tmp = t_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 (k <= (-4.7d-303)) then
tmp = t_0
else if (k <= 235000.0d0) then
tmp = a
else
tmp = t_0
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 (k <= -4.7e-303) {
tmp = t_0;
} else if (k <= 235000.0) {
tmp = a;
} else {
tmp = t_0;
}
return tmp;
}
def code(a, k, m): t_0 = a / (k * k) tmp = 0 if k <= -4.7e-303: tmp = t_0 elif k <= 235000.0: tmp = a else: tmp = t_0 return tmp
function code(a, k, m) t_0 = Float64(a / Float64(k * k)) tmp = 0.0 if (k <= -4.7e-303) tmp = t_0; elseif (k <= 235000.0) tmp = a; else tmp = t_0; end return tmp end
function tmp_2 = code(a, k, m) t_0 = a / (k * k); tmp = 0.0; if (k <= -4.7e-303) tmp = t_0; elseif (k <= 235000.0) tmp = a; else tmp = t_0; end tmp_2 = tmp; end
code[a_, k_, m_] := Block[{t$95$0 = N[(a / N[(k * k), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[k, -4.7e-303], t$95$0, If[LessEqual[k, 235000.0], a, t$95$0]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{a}{k \cdot k}\\
\mathbf{if}\;k \leq -4.7 \cdot 10^{-303}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;k \leq 235000:\\
\;\;\;\;a\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if k < -4.6999999999999997e-303 or 235000 < k Initial program 84.2%
Taylor expanded in m around 0
/-lowering-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
associate-*l*N/A
unpow2N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
Simplified38.4%
Taylor expanded in k around inf
unpow2N/A
*-lowering-*.f6441.5%
Simplified41.5%
if -4.6999999999999997e-303 < k < 235000Initial program 100.0%
Taylor expanded in m around 0
/-lowering-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
associate-*l*N/A
unpow2N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
Simplified41.4%
Taylor expanded in k around 0
Simplified41.3%
(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 89.5%
Taylor expanded in m around 0
/-lowering-/.f64N/A
unpow2N/A
distribute-rgt-inN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
*-lft-identityN/A
distribute-rgt-inN/A
+-commutativeN/A
associate-*l*N/A
unpow2N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
*-lft-identityN/A
Simplified39.4%
Taylor expanded in k around 0
Simplified16.7%
herbie shell --seed 2024191
(FPCore (a k m)
:name "Falkner and Boettcher, Appendix A"
:precision binary64
(/ (* a (pow k m)) (+ (+ 1.0 (* 10.0 k)) (* k k))))