
(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 18 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}
a\_m = (fabs.f64 a)
a\_s = (copysign.f64 #s(literal 1 binary64) a)
(FPCore (a_s a_m k m)
:precision binary64
(*
a_s
(if (<= (/ (* a_m (pow k m)) (+ (+ 1.0 (* k 10.0)) (* k k))) 2e+63)
(/ a_m (/ (+ 1.0 (* k (+ k 10.0))) (pow k m)))
(/ a_m (pow k (- 0.0 m))))))a\_m = fabs(a);
a\_s = copysign(1.0, a);
double code(double a_s, double a_m, double k, double m) {
double tmp;
if (((a_m * pow(k, m)) / ((1.0 + (k * 10.0)) + (k * k))) <= 2e+63) {
tmp = a_m / ((1.0 + (k * (k + 10.0))) / pow(k, m));
} else {
tmp = a_m / pow(k, (0.0 - m));
}
return a_s * tmp;
}
a\_m = abs(a)
a\_s = copysign(1.0d0, a)
real(8) function code(a_s, a_m, k, m)
real(8), intent (in) :: a_s
real(8), intent (in) :: a_m
real(8), intent (in) :: k
real(8), intent (in) :: m
real(8) :: tmp
if (((a_m * (k ** m)) / ((1.0d0 + (k * 10.0d0)) + (k * k))) <= 2d+63) then
tmp = a_m / ((1.0d0 + (k * (k + 10.0d0))) / (k ** m))
else
tmp = a_m / (k ** (0.0d0 - m))
end if
code = a_s * tmp
end function
a\_m = Math.abs(a);
a\_s = Math.copySign(1.0, a);
public static double code(double a_s, double a_m, double k, double m) {
double tmp;
if (((a_m * Math.pow(k, m)) / ((1.0 + (k * 10.0)) + (k * k))) <= 2e+63) {
tmp = a_m / ((1.0 + (k * (k + 10.0))) / Math.pow(k, m));
} else {
tmp = a_m / Math.pow(k, (0.0 - m));
}
return a_s * tmp;
}
a\_m = math.fabs(a) a\_s = math.copysign(1.0, a) def code(a_s, a_m, k, m): tmp = 0 if ((a_m * math.pow(k, m)) / ((1.0 + (k * 10.0)) + (k * k))) <= 2e+63: tmp = a_m / ((1.0 + (k * (k + 10.0))) / math.pow(k, m)) else: tmp = a_m / math.pow(k, (0.0 - m)) return a_s * tmp
a\_m = abs(a) a\_s = copysign(1.0, a) function code(a_s, a_m, k, m) tmp = 0.0 if (Float64(Float64(a_m * (k ^ m)) / Float64(Float64(1.0 + Float64(k * 10.0)) + Float64(k * k))) <= 2e+63) tmp = Float64(a_m / Float64(Float64(1.0 + Float64(k * Float64(k + 10.0))) / (k ^ m))); else tmp = Float64(a_m / (k ^ Float64(0.0 - m))); end return Float64(a_s * tmp) end
a\_m = abs(a); a\_s = sign(a) * abs(1.0); function tmp_2 = code(a_s, a_m, k, m) tmp = 0.0; if (((a_m * (k ^ m)) / ((1.0 + (k * 10.0)) + (k * k))) <= 2e+63) tmp = a_m / ((1.0 + (k * (k + 10.0))) / (k ^ m)); else tmp = a_m / (k ^ (0.0 - m)); end tmp_2 = a_s * tmp; end
a\_m = N[Abs[a], $MachinePrecision]
a\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[a]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[a$95$s_, a$95$m_, k_, m_] := N[(a$95$s * If[LessEqual[N[(N[(a$95$m * N[Power[k, m], $MachinePrecision]), $MachinePrecision] / N[(N[(1.0 + N[(k * 10.0), $MachinePrecision]), $MachinePrecision] + N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2e+63], N[(a$95$m / N[(N[(1.0 + N[(k * N[(k + 10.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[Power[k, m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(a$95$m / N[Power[k, N[(0.0 - m), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
a\_m = \left|a\right|
\\
a\_s = \mathsf{copysign}\left(1, a\right)
\\
a\_s \cdot \begin{array}{l}
\mathbf{if}\;\frac{a\_m \cdot {k}^{m}}{\left(1 + k \cdot 10\right) + k \cdot k} \leq 2 \cdot 10^{+63}:\\
\;\;\;\;\frac{a\_m}{\frac{1 + k \cdot \left(k + 10\right)}{{k}^{m}}}\\
\mathbf{else}:\\
\;\;\;\;\frac{a\_m}{{k}^{\left(0 - m\right)}}\\
\end{array}
\end{array}
if (/.f64 (*.f64 a (pow.f64 k m)) (+.f64 (+.f64 #s(literal 1 binary64) (*.f64 #s(literal 10 binary64) k)) (*.f64 k k))) < 2.00000000000000012e63Initial program 99.9%
/-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
+-commutativeN/A
+-lowering-+.f6499.9%
Simplified99.9%
associate-/l*N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
pow-lowering-pow.f6499.9%
Applied egg-rr99.9%
if 2.00000000000000012e63 < (/.f64 (*.f64 a (pow.f64 k m)) (+.f64 (+.f64 #s(literal 1 binary64) (*.f64 #s(literal 10 binary64) k)) (*.f64 k k))) Initial program 64.9%
/-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
+-commutativeN/A
+-lowering-+.f6464.9%
Simplified64.9%
Taylor expanded in k around inf
mul-1-negN/A
exp-negN/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
exp-to-powN/A
pow-lowering-pow.f64N/A
/-lowering-/.f6464.9%
Simplified64.9%
Taylor expanded in k around 0
Simplified100.0%
/-rgt-identityN/A
*-rgt-identityN/A
/-lowering-/.f64N/A
inv-powN/A
pow-powN/A
pow-lowering-pow.f64N/A
*-lowering-*.f64100.0%
Applied egg-rr100.0%
mul-1-negN/A
neg-lowering-neg.f64100.0%
Applied egg-rr100.0%
Final simplification100.0%
a\_m = (fabs.f64 a)
a\_s = (copysign.f64 #s(literal 1 binary64) a)
(FPCore (a_s a_m k m)
:precision binary64
(*
a_s
(if (<= m -22000000000.0)
(* a_m (pow k m))
(if (<= m 7.5e-17)
(/ a_m (+ 1.0 (* k (+ k 10.0))))
(/ a_m (pow k (- 0.0 m)))))))a\_m = fabs(a);
a\_s = copysign(1.0, a);
double code(double a_s, double a_m, double k, double m) {
double tmp;
if (m <= -22000000000.0) {
tmp = a_m * pow(k, m);
} else if (m <= 7.5e-17) {
tmp = a_m / (1.0 + (k * (k + 10.0)));
} else {
tmp = a_m / pow(k, (0.0 - m));
}
return a_s * tmp;
}
a\_m = abs(a)
a\_s = copysign(1.0d0, a)
real(8) function code(a_s, a_m, k, m)
real(8), intent (in) :: a_s
real(8), intent (in) :: a_m
real(8), intent (in) :: k
real(8), intent (in) :: m
real(8) :: tmp
if (m <= (-22000000000.0d0)) then
tmp = a_m * (k ** m)
else if (m <= 7.5d-17) then
tmp = a_m / (1.0d0 + (k * (k + 10.0d0)))
else
tmp = a_m / (k ** (0.0d0 - m))
end if
code = a_s * tmp
end function
a\_m = Math.abs(a);
a\_s = Math.copySign(1.0, a);
public static double code(double a_s, double a_m, double k, double m) {
double tmp;
if (m <= -22000000000.0) {
tmp = a_m * Math.pow(k, m);
} else if (m <= 7.5e-17) {
tmp = a_m / (1.0 + (k * (k + 10.0)));
} else {
tmp = a_m / Math.pow(k, (0.0 - m));
}
return a_s * tmp;
}
a\_m = math.fabs(a) a\_s = math.copysign(1.0, a) def code(a_s, a_m, k, m): tmp = 0 if m <= -22000000000.0: tmp = a_m * math.pow(k, m) elif m <= 7.5e-17: tmp = a_m / (1.0 + (k * (k + 10.0))) else: tmp = a_m / math.pow(k, (0.0 - m)) return a_s * tmp
a\_m = abs(a) a\_s = copysign(1.0, a) function code(a_s, a_m, k, m) tmp = 0.0 if (m <= -22000000000.0) tmp = Float64(a_m * (k ^ m)); elseif (m <= 7.5e-17) tmp = Float64(a_m / Float64(1.0 + Float64(k * Float64(k + 10.0)))); else tmp = Float64(a_m / (k ^ Float64(0.0 - m))); end return Float64(a_s * tmp) end
a\_m = abs(a); a\_s = sign(a) * abs(1.0); function tmp_2 = code(a_s, a_m, k, m) tmp = 0.0; if (m <= -22000000000.0) tmp = a_m * (k ^ m); elseif (m <= 7.5e-17) tmp = a_m / (1.0 + (k * (k + 10.0))); else tmp = a_m / (k ^ (0.0 - m)); end tmp_2 = a_s * tmp; end
a\_m = N[Abs[a], $MachinePrecision]
a\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[a]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[a$95$s_, a$95$m_, k_, m_] := N[(a$95$s * If[LessEqual[m, -22000000000.0], N[(a$95$m * N[Power[k, m], $MachinePrecision]), $MachinePrecision], If[LessEqual[m, 7.5e-17], N[(a$95$m / N[(1.0 + N[(k * N[(k + 10.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(a$95$m / N[Power[k, N[(0.0 - m), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]
\begin{array}{l}
a\_m = \left|a\right|
\\
a\_s = \mathsf{copysign}\left(1, a\right)
\\
a\_s \cdot \begin{array}{l}
\mathbf{if}\;m \leq -22000000000:\\
\;\;\;\;a\_m \cdot {k}^{m}\\
\mathbf{elif}\;m \leq 7.5 \cdot 10^{-17}:\\
\;\;\;\;\frac{a\_m}{1 + k \cdot \left(k + 10\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{a\_m}{{k}^{\left(0 - m\right)}}\\
\end{array}
\end{array}
if m < -2.2e10Initial program 100.0%
/-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
+-commutativeN/A
+-lowering-+.f64100.0%
Simplified100.0%
Taylor expanded in k around 0
*-lowering-*.f64N/A
pow-lowering-pow.f64100.0%
Simplified100.0%
if -2.2e10 < m < 7.49999999999999984e-17Initial program 99.9%
/-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
+-commutativeN/A
+-lowering-+.f6499.9%
Simplified99.9%
Taylor expanded in m around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f6499.5%
Simplified99.5%
if 7.49999999999999984e-17 < m Initial program 77.2%
/-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
+-commutativeN/A
+-lowering-+.f6477.2%
Simplified77.2%
Taylor expanded in k around inf
mul-1-negN/A
exp-negN/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
exp-to-powN/A
pow-lowering-pow.f64N/A
/-lowering-/.f6477.3%
Simplified77.3%
Taylor expanded in k around 0
Simplified100.0%
/-rgt-identityN/A
*-rgt-identityN/A
/-lowering-/.f64N/A
inv-powN/A
pow-powN/A
pow-lowering-pow.f64N/A
*-lowering-*.f64100.0%
Applied egg-rr100.0%
mul-1-negN/A
neg-lowering-neg.f64100.0%
Applied egg-rr100.0%
Final simplification99.8%
a\_m = (fabs.f64 a)
a\_s = (copysign.f64 #s(literal 1 binary64) a)
(FPCore (a_s a_m k m)
:precision binary64
(let* ((t_0 (* a_m (pow k m))))
(*
a_s
(if (<= m -22000000000.0)
t_0
(if (<= m 7.5e-17) (/ a_m (+ 1.0 (* k (+ k 10.0)))) t_0)))))a\_m = fabs(a);
a\_s = copysign(1.0, a);
double code(double a_s, double a_m, double k, double m) {
double t_0 = a_m * pow(k, m);
double tmp;
if (m <= -22000000000.0) {
tmp = t_0;
} else if (m <= 7.5e-17) {
tmp = a_m / (1.0 + (k * (k + 10.0)));
} else {
tmp = t_0;
}
return a_s * tmp;
}
a\_m = abs(a)
a\_s = copysign(1.0d0, a)
real(8) function code(a_s, a_m, k, m)
real(8), intent (in) :: a_s
real(8), intent (in) :: a_m
real(8), intent (in) :: k
real(8), intent (in) :: m
real(8) :: t_0
real(8) :: tmp
t_0 = a_m * (k ** m)
if (m <= (-22000000000.0d0)) then
tmp = t_0
else if (m <= 7.5d-17) then
tmp = a_m / (1.0d0 + (k * (k + 10.0d0)))
else
tmp = t_0
end if
code = a_s * tmp
end function
a\_m = Math.abs(a);
a\_s = Math.copySign(1.0, a);
public static double code(double a_s, double a_m, double k, double m) {
double t_0 = a_m * Math.pow(k, m);
double tmp;
if (m <= -22000000000.0) {
tmp = t_0;
} else if (m <= 7.5e-17) {
tmp = a_m / (1.0 + (k * (k + 10.0)));
} else {
tmp = t_0;
}
return a_s * tmp;
}
a\_m = math.fabs(a) a\_s = math.copysign(1.0, a) def code(a_s, a_m, k, m): t_0 = a_m * math.pow(k, m) tmp = 0 if m <= -22000000000.0: tmp = t_0 elif m <= 7.5e-17: tmp = a_m / (1.0 + (k * (k + 10.0))) else: tmp = t_0 return a_s * tmp
a\_m = abs(a) a\_s = copysign(1.0, a) function code(a_s, a_m, k, m) t_0 = Float64(a_m * (k ^ m)) tmp = 0.0 if (m <= -22000000000.0) tmp = t_0; elseif (m <= 7.5e-17) tmp = Float64(a_m / Float64(1.0 + Float64(k * Float64(k + 10.0)))); else tmp = t_0; end return Float64(a_s * tmp) end
a\_m = abs(a); a\_s = sign(a) * abs(1.0); function tmp_2 = code(a_s, a_m, k, m) t_0 = a_m * (k ^ m); tmp = 0.0; if (m <= -22000000000.0) tmp = t_0; elseif (m <= 7.5e-17) tmp = a_m / (1.0 + (k * (k + 10.0))); else tmp = t_0; end tmp_2 = a_s * tmp; end
a\_m = N[Abs[a], $MachinePrecision]
a\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[a]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[a$95$s_, a$95$m_, k_, m_] := Block[{t$95$0 = N[(a$95$m * N[Power[k, m], $MachinePrecision]), $MachinePrecision]}, N[(a$95$s * If[LessEqual[m, -22000000000.0], t$95$0, If[LessEqual[m, 7.5e-17], N[(a$95$m / N[(1.0 + N[(k * N[(k + 10.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$0]]), $MachinePrecision]]
\begin{array}{l}
a\_m = \left|a\right|
\\
a\_s = \mathsf{copysign}\left(1, a\right)
\\
\begin{array}{l}
t_0 := a\_m \cdot {k}^{m}\\
a\_s \cdot \begin{array}{l}
\mathbf{if}\;m \leq -22000000000:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;m \leq 7.5 \cdot 10^{-17}:\\
\;\;\;\;\frac{a\_m}{1 + k \cdot \left(k + 10\right)}\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
\end{array}
if m < -2.2e10 or 7.49999999999999984e-17 < m Initial program 88.1%
/-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
+-commutativeN/A
+-lowering-+.f6488.1%
Simplified88.1%
Taylor expanded in k around 0
*-lowering-*.f64N/A
pow-lowering-pow.f64100.0%
Simplified100.0%
if -2.2e10 < m < 7.49999999999999984e-17Initial program 99.9%
/-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
+-commutativeN/A
+-lowering-+.f6499.9%
Simplified99.9%
Taylor expanded in m around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f6499.5%
Simplified99.5%
a\_m = (fabs.f64 a) a\_s = (copysign.f64 #s(literal 1 binary64) a) (FPCore (a_s a_m k m) :precision binary64 (* a_s (if (<= k 1.0) (/ a_m (pow k (- 0.0 m))) (/ (/ (* a_m (pow k m)) k) k))))
a\_m = fabs(a);
a\_s = copysign(1.0, a);
double code(double a_s, double a_m, double k, double m) {
double tmp;
if (k <= 1.0) {
tmp = a_m / pow(k, (0.0 - m));
} else {
tmp = ((a_m * pow(k, m)) / k) / k;
}
return a_s * tmp;
}
a\_m = abs(a)
a\_s = copysign(1.0d0, a)
real(8) function code(a_s, a_m, k, m)
real(8), intent (in) :: a_s
real(8), intent (in) :: a_m
real(8), intent (in) :: k
real(8), intent (in) :: m
real(8) :: tmp
if (k <= 1.0d0) then
tmp = a_m / (k ** (0.0d0 - m))
else
tmp = ((a_m * (k ** m)) / k) / k
end if
code = a_s * tmp
end function
a\_m = Math.abs(a);
a\_s = Math.copySign(1.0, a);
public static double code(double a_s, double a_m, double k, double m) {
double tmp;
if (k <= 1.0) {
tmp = a_m / Math.pow(k, (0.0 - m));
} else {
tmp = ((a_m * Math.pow(k, m)) / k) / k;
}
return a_s * tmp;
}
a\_m = math.fabs(a) a\_s = math.copysign(1.0, a) def code(a_s, a_m, k, m): tmp = 0 if k <= 1.0: tmp = a_m / math.pow(k, (0.0 - m)) else: tmp = ((a_m * math.pow(k, m)) / k) / k return a_s * tmp
a\_m = abs(a) a\_s = copysign(1.0, a) function code(a_s, a_m, k, m) tmp = 0.0 if (k <= 1.0) tmp = Float64(a_m / (k ^ Float64(0.0 - m))); else tmp = Float64(Float64(Float64(a_m * (k ^ m)) / k) / k); end return Float64(a_s * tmp) end
a\_m = abs(a); a\_s = sign(a) * abs(1.0); function tmp_2 = code(a_s, a_m, k, m) tmp = 0.0; if (k <= 1.0) tmp = a_m / (k ^ (0.0 - m)); else tmp = ((a_m * (k ^ m)) / k) / k; end tmp_2 = a_s * tmp; end
a\_m = N[Abs[a], $MachinePrecision]
a\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[a]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[a$95$s_, a$95$m_, k_, m_] := N[(a$95$s * If[LessEqual[k, 1.0], N[(a$95$m / N[Power[k, N[(0.0 - m), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], N[(N[(N[(a$95$m * N[Power[k, m], $MachinePrecision]), $MachinePrecision] / k), $MachinePrecision] / k), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
a\_m = \left|a\right|
\\
a\_s = \mathsf{copysign}\left(1, a\right)
\\
a\_s \cdot \begin{array}{l}
\mathbf{if}\;k \leq 1:\\
\;\;\;\;\frac{a\_m}{{k}^{\left(0 - m\right)}}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{a\_m \cdot {k}^{m}}{k}}{k}\\
\end{array}
\end{array}
if k < 1Initial program 94.8%
/-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
+-commutativeN/A
+-lowering-+.f6494.8%
Simplified94.8%
Taylor expanded in k around inf
mul-1-negN/A
exp-negN/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
exp-to-powN/A
pow-lowering-pow.f64N/A
/-lowering-/.f6494.8%
Simplified94.8%
Taylor expanded in k around 0
Simplified98.9%
/-rgt-identityN/A
*-rgt-identityN/A
/-lowering-/.f64N/A
inv-powN/A
pow-powN/A
pow-lowering-pow.f64N/A
*-lowering-*.f6498.9%
Applied egg-rr98.9%
mul-1-negN/A
neg-lowering-neg.f6498.9%
Applied egg-rr98.9%
if 1 < k Initial program 86.3%
/-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
+-commutativeN/A
+-lowering-+.f6486.3%
Simplified86.3%
Taylor expanded in k around inf
unpow2N/A
*-lowering-*.f6485.9%
Simplified85.9%
Taylor expanded in a around 0
unpow2N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
pow-lowering-pow.f6499.5%
Simplified99.5%
Final simplification99.1%
a\_m = (fabs.f64 a) a\_s = (copysign.f64 #s(literal 1 binary64) a) (FPCore (a_s a_m k m) :precision binary64 (* a_s (if (<= k 1.0) (/ a_m (pow k (- 0.0 m))) (* a_m (pow k (- m 2.0))))))
a\_m = fabs(a);
a\_s = copysign(1.0, a);
double code(double a_s, double a_m, double k, double m) {
double tmp;
if (k <= 1.0) {
tmp = a_m / pow(k, (0.0 - m));
} else {
tmp = a_m * pow(k, (m - 2.0));
}
return a_s * tmp;
}
a\_m = abs(a)
a\_s = copysign(1.0d0, a)
real(8) function code(a_s, a_m, k, m)
real(8), intent (in) :: a_s
real(8), intent (in) :: a_m
real(8), intent (in) :: k
real(8), intent (in) :: m
real(8) :: tmp
if (k <= 1.0d0) then
tmp = a_m / (k ** (0.0d0 - m))
else
tmp = a_m * (k ** (m - 2.0d0))
end if
code = a_s * tmp
end function
a\_m = Math.abs(a);
a\_s = Math.copySign(1.0, a);
public static double code(double a_s, double a_m, double k, double m) {
double tmp;
if (k <= 1.0) {
tmp = a_m / Math.pow(k, (0.0 - m));
} else {
tmp = a_m * Math.pow(k, (m - 2.0));
}
return a_s * tmp;
}
a\_m = math.fabs(a) a\_s = math.copysign(1.0, a) def code(a_s, a_m, k, m): tmp = 0 if k <= 1.0: tmp = a_m / math.pow(k, (0.0 - m)) else: tmp = a_m * math.pow(k, (m - 2.0)) return a_s * tmp
a\_m = abs(a) a\_s = copysign(1.0, a) function code(a_s, a_m, k, m) tmp = 0.0 if (k <= 1.0) tmp = Float64(a_m / (k ^ Float64(0.0 - m))); else tmp = Float64(a_m * (k ^ Float64(m - 2.0))); end return Float64(a_s * tmp) end
a\_m = abs(a); a\_s = sign(a) * abs(1.0); function tmp_2 = code(a_s, a_m, k, m) tmp = 0.0; if (k <= 1.0) tmp = a_m / (k ^ (0.0 - m)); else tmp = a_m * (k ^ (m - 2.0)); end tmp_2 = a_s * tmp; end
a\_m = N[Abs[a], $MachinePrecision]
a\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[a]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[a$95$s_, a$95$m_, k_, m_] := N[(a$95$s * If[LessEqual[k, 1.0], N[(a$95$m / N[Power[k, N[(0.0 - m), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], N[(a$95$m * N[Power[k, N[(m - 2.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
a\_m = \left|a\right|
\\
a\_s = \mathsf{copysign}\left(1, a\right)
\\
a\_s \cdot \begin{array}{l}
\mathbf{if}\;k \leq 1:\\
\;\;\;\;\frac{a\_m}{{k}^{\left(0 - m\right)}}\\
\mathbf{else}:\\
\;\;\;\;a\_m \cdot {k}^{\left(m - 2\right)}\\
\end{array}
\end{array}
if k < 1Initial program 94.8%
/-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
+-commutativeN/A
+-lowering-+.f6494.8%
Simplified94.8%
Taylor expanded in k around inf
mul-1-negN/A
exp-negN/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
exp-to-powN/A
pow-lowering-pow.f64N/A
/-lowering-/.f6494.8%
Simplified94.8%
Taylor expanded in k around 0
Simplified98.9%
/-rgt-identityN/A
*-rgt-identityN/A
/-lowering-/.f64N/A
inv-powN/A
pow-powN/A
pow-lowering-pow.f64N/A
*-lowering-*.f6498.9%
Applied egg-rr98.9%
mul-1-negN/A
neg-lowering-neg.f6498.9%
Applied egg-rr98.9%
if 1 < k Initial program 86.3%
/-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
+-commutativeN/A
+-lowering-+.f6486.3%
Simplified86.3%
Taylor expanded in k around inf
unpow2N/A
*-lowering-*.f6485.9%
Simplified85.9%
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
pow2N/A
pow-divN/A
pow-lowering-pow.f64N/A
--lowering--.f6499.3%
Applied egg-rr99.3%
Final simplification99.0%
a\_m = (fabs.f64 a)
a\_s = (copysign.f64 #s(literal 1 binary64) a)
(FPCore (a_s a_m k m)
:precision binary64
(let* ((t_0 (* k (* k k))))
(*
a_s
(if (<= m -22000000000.0)
(/ (* a_m (/ 99.0 (* k k))) (* k k))
(if (<= m 0.68)
(/ a_m (+ 1.0 (* k (+ k 10.0))))
(/
a_m
(+
1.0
(*
(/ (- (/ (+ 1000.0 (/ 1000000.0 t_0)) t_0) -1.0) t_0)
(/
(* k (+ (* k k) -100.0))
(/ 1.0 (- (+ (* k k) 100.0) (* k -10.0))))))))))))a\_m = fabs(a);
a\_s = copysign(1.0, a);
double code(double a_s, double a_m, double k, double m) {
double t_0 = k * (k * k);
double tmp;
if (m <= -22000000000.0) {
tmp = (a_m * (99.0 / (k * k))) / (k * k);
} else if (m <= 0.68) {
tmp = a_m / (1.0 + (k * (k + 10.0)));
} else {
tmp = a_m / (1.0 + (((((1000.0 + (1000000.0 / t_0)) / t_0) - -1.0) / t_0) * ((k * ((k * k) + -100.0)) / (1.0 / (((k * k) + 100.0) - (k * -10.0))))));
}
return a_s * tmp;
}
a\_m = abs(a)
a\_s = copysign(1.0d0, a)
real(8) function code(a_s, a_m, k, m)
real(8), intent (in) :: a_s
real(8), intent (in) :: a_m
real(8), intent (in) :: k
real(8), intent (in) :: m
real(8) :: t_0
real(8) :: tmp
t_0 = k * (k * k)
if (m <= (-22000000000.0d0)) then
tmp = (a_m * (99.0d0 / (k * k))) / (k * k)
else if (m <= 0.68d0) then
tmp = a_m / (1.0d0 + (k * (k + 10.0d0)))
else
tmp = a_m / (1.0d0 + (((((1000.0d0 + (1000000.0d0 / t_0)) / t_0) - (-1.0d0)) / t_0) * ((k * ((k * k) + (-100.0d0))) / (1.0d0 / (((k * k) + 100.0d0) - (k * (-10.0d0)))))))
end if
code = a_s * tmp
end function
a\_m = Math.abs(a);
a\_s = Math.copySign(1.0, a);
public static double code(double a_s, double a_m, double k, double m) {
double t_0 = k * (k * k);
double tmp;
if (m <= -22000000000.0) {
tmp = (a_m * (99.0 / (k * k))) / (k * k);
} else if (m <= 0.68) {
tmp = a_m / (1.0 + (k * (k + 10.0)));
} else {
tmp = a_m / (1.0 + (((((1000.0 + (1000000.0 / t_0)) / t_0) - -1.0) / t_0) * ((k * ((k * k) + -100.0)) / (1.0 / (((k * k) + 100.0) - (k * -10.0))))));
}
return a_s * tmp;
}
a\_m = math.fabs(a) a\_s = math.copysign(1.0, a) def code(a_s, a_m, k, m): t_0 = k * (k * k) tmp = 0 if m <= -22000000000.0: tmp = (a_m * (99.0 / (k * k))) / (k * k) elif m <= 0.68: tmp = a_m / (1.0 + (k * (k + 10.0))) else: tmp = a_m / (1.0 + (((((1000.0 + (1000000.0 / t_0)) / t_0) - -1.0) / t_0) * ((k * ((k * k) + -100.0)) / (1.0 / (((k * k) + 100.0) - (k * -10.0)))))) return a_s * tmp
a\_m = abs(a) a\_s = copysign(1.0, a) function code(a_s, a_m, k, m) t_0 = Float64(k * Float64(k * k)) tmp = 0.0 if (m <= -22000000000.0) tmp = Float64(Float64(a_m * Float64(99.0 / Float64(k * k))) / Float64(k * k)); elseif (m <= 0.68) tmp = Float64(a_m / Float64(1.0 + Float64(k * Float64(k + 10.0)))); else tmp = Float64(a_m / Float64(1.0 + Float64(Float64(Float64(Float64(Float64(1000.0 + Float64(1000000.0 / t_0)) / t_0) - -1.0) / t_0) * Float64(Float64(k * Float64(Float64(k * k) + -100.0)) / Float64(1.0 / Float64(Float64(Float64(k * k) + 100.0) - Float64(k * -10.0))))))); end return Float64(a_s * tmp) end
a\_m = abs(a); a\_s = sign(a) * abs(1.0); function tmp_2 = code(a_s, a_m, k, m) t_0 = k * (k * k); tmp = 0.0; if (m <= -22000000000.0) tmp = (a_m * (99.0 / (k * k))) / (k * k); elseif (m <= 0.68) tmp = a_m / (1.0 + (k * (k + 10.0))); else tmp = a_m / (1.0 + (((((1000.0 + (1000000.0 / t_0)) / t_0) - -1.0) / t_0) * ((k * ((k * k) + -100.0)) / (1.0 / (((k * k) + 100.0) - (k * -10.0)))))); end tmp_2 = a_s * tmp; end
a\_m = N[Abs[a], $MachinePrecision]
a\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[a]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[a$95$s_, a$95$m_, k_, m_] := Block[{t$95$0 = N[(k * N[(k * k), $MachinePrecision]), $MachinePrecision]}, N[(a$95$s * If[LessEqual[m, -22000000000.0], N[(N[(a$95$m * N[(99.0 / N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(k * k), $MachinePrecision]), $MachinePrecision], If[LessEqual[m, 0.68], N[(a$95$m / N[(1.0 + N[(k * N[(k + 10.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(a$95$m / N[(1.0 + N[(N[(N[(N[(N[(1000.0 + N[(1000000.0 / t$95$0), $MachinePrecision]), $MachinePrecision] / t$95$0), $MachinePrecision] - -1.0), $MachinePrecision] / t$95$0), $MachinePrecision] * N[(N[(k * N[(N[(k * k), $MachinePrecision] + -100.0), $MachinePrecision]), $MachinePrecision] / N[(1.0 / N[(N[(N[(k * k), $MachinePrecision] + 100.0), $MachinePrecision] - N[(k * -10.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]]
\begin{array}{l}
a\_m = \left|a\right|
\\
a\_s = \mathsf{copysign}\left(1, a\right)
\\
\begin{array}{l}
t_0 := k \cdot \left(k \cdot k\right)\\
a\_s \cdot \begin{array}{l}
\mathbf{if}\;m \leq -22000000000:\\
\;\;\;\;\frac{a\_m \cdot \frac{99}{k \cdot k}}{k \cdot k}\\
\mathbf{elif}\;m \leq 0.68:\\
\;\;\;\;\frac{a\_m}{1 + k \cdot \left(k + 10\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{a\_m}{1 + \frac{\frac{1000 + \frac{1000000}{t\_0}}{t\_0} - -1}{t\_0} \cdot \frac{k \cdot \left(k \cdot k + -100\right)}{\frac{1}{\left(k \cdot k + 100\right) - k \cdot -10}}}\\
\end{array}
\end{array}
\end{array}
if m < -2.2e10Initial program 100.0%
/-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
+-commutativeN/A
+-lowering-+.f64100.0%
Simplified100.0%
Taylor expanded in m around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f6433.6%
Simplified33.6%
Taylor expanded in k around -inf
/-lowering-/.f64N/A
Simplified61.3%
Taylor expanded in k around 0
*-commutativeN/A
associate-*l/N/A
associate-/l*N/A
metadata-evalN/A
associate-*r/N/A
*-lowering-*.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6468.7%
Simplified68.7%
if -2.2e10 < m < 0.680000000000000049Initial program 99.8%
/-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
+-commutativeN/A
+-lowering-+.f6499.9%
Simplified99.9%
Taylor expanded in m around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f6498.2%
Simplified98.2%
if 0.680000000000000049 < m Initial program 76.5%
/-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
+-commutativeN/A
+-lowering-+.f6476.5%
Simplified76.5%
Taylor expanded in m around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f642.8%
Simplified2.8%
metadata-evalN/A
sub-negN/A
flip--N/A
sub-negN/A
metadata-evalN/A
metadata-evalN/A
associate-*r/N/A
div-invN/A
*-commutativeN/A
associate-*l/N/A
flip3-+N/A
div-invN/A
times-fracN/A
*-lowering-*.f64N/A
Applied egg-rr2.3%
Taylor expanded in k around -inf
mul-1-negN/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
Simplified50.8%
Final simplification73.1%
a\_m = (fabs.f64 a)
a\_s = (copysign.f64 #s(literal 1 binary64) a)
(FPCore (a_s a_m k m)
:precision binary64
(let* ((t_0 (* k (* k k))))
(*
a_s
(if (<= m -22000000000.0)
(/ (* a_m (/ 99.0 (* k k))) (* k k))
(if (<= m 1.2)
(/ a_m (+ 1.0 (* k (+ k 10.0))))
(/
a_m
(+
1.0
(*
(/
(* k (+ (* k k) -100.0))
(/ 1.0 (- (+ (* k k) 100.0) (* k -10.0))))
(/ (+ 1.0 (/ 1000.0 t_0)) t_0)))))))))a\_m = fabs(a);
a\_s = copysign(1.0, a);
double code(double a_s, double a_m, double k, double m) {
double t_0 = k * (k * k);
double tmp;
if (m <= -22000000000.0) {
tmp = (a_m * (99.0 / (k * k))) / (k * k);
} else if (m <= 1.2) {
tmp = a_m / (1.0 + (k * (k + 10.0)));
} else {
tmp = a_m / (1.0 + (((k * ((k * k) + -100.0)) / (1.0 / (((k * k) + 100.0) - (k * -10.0)))) * ((1.0 + (1000.0 / t_0)) / t_0)));
}
return a_s * tmp;
}
a\_m = abs(a)
a\_s = copysign(1.0d0, a)
real(8) function code(a_s, a_m, k, m)
real(8), intent (in) :: a_s
real(8), intent (in) :: a_m
real(8), intent (in) :: k
real(8), intent (in) :: m
real(8) :: t_0
real(8) :: tmp
t_0 = k * (k * k)
if (m <= (-22000000000.0d0)) then
tmp = (a_m * (99.0d0 / (k * k))) / (k * k)
else if (m <= 1.2d0) then
tmp = a_m / (1.0d0 + (k * (k + 10.0d0)))
else
tmp = a_m / (1.0d0 + (((k * ((k * k) + (-100.0d0))) / (1.0d0 / (((k * k) + 100.0d0) - (k * (-10.0d0))))) * ((1.0d0 + (1000.0d0 / t_0)) / t_0)))
end if
code = a_s * tmp
end function
a\_m = Math.abs(a);
a\_s = Math.copySign(1.0, a);
public static double code(double a_s, double a_m, double k, double m) {
double t_0 = k * (k * k);
double tmp;
if (m <= -22000000000.0) {
tmp = (a_m * (99.0 / (k * k))) / (k * k);
} else if (m <= 1.2) {
tmp = a_m / (1.0 + (k * (k + 10.0)));
} else {
tmp = a_m / (1.0 + (((k * ((k * k) + -100.0)) / (1.0 / (((k * k) + 100.0) - (k * -10.0)))) * ((1.0 + (1000.0 / t_0)) / t_0)));
}
return a_s * tmp;
}
a\_m = math.fabs(a) a\_s = math.copysign(1.0, a) def code(a_s, a_m, k, m): t_0 = k * (k * k) tmp = 0 if m <= -22000000000.0: tmp = (a_m * (99.0 / (k * k))) / (k * k) elif m <= 1.2: tmp = a_m / (1.0 + (k * (k + 10.0))) else: tmp = a_m / (1.0 + (((k * ((k * k) + -100.0)) / (1.0 / (((k * k) + 100.0) - (k * -10.0)))) * ((1.0 + (1000.0 / t_0)) / t_0))) return a_s * tmp
a\_m = abs(a) a\_s = copysign(1.0, a) function code(a_s, a_m, k, m) t_0 = Float64(k * Float64(k * k)) tmp = 0.0 if (m <= -22000000000.0) tmp = Float64(Float64(a_m * Float64(99.0 / Float64(k * k))) / Float64(k * k)); elseif (m <= 1.2) tmp = Float64(a_m / Float64(1.0 + Float64(k * Float64(k + 10.0)))); else tmp = Float64(a_m / Float64(1.0 + Float64(Float64(Float64(k * Float64(Float64(k * k) + -100.0)) / Float64(1.0 / Float64(Float64(Float64(k * k) + 100.0) - Float64(k * -10.0)))) * Float64(Float64(1.0 + Float64(1000.0 / t_0)) / t_0)))); end return Float64(a_s * tmp) end
a\_m = abs(a); a\_s = sign(a) * abs(1.0); function tmp_2 = code(a_s, a_m, k, m) t_0 = k * (k * k); tmp = 0.0; if (m <= -22000000000.0) tmp = (a_m * (99.0 / (k * k))) / (k * k); elseif (m <= 1.2) tmp = a_m / (1.0 + (k * (k + 10.0))); else tmp = a_m / (1.0 + (((k * ((k * k) + -100.0)) / (1.0 / (((k * k) + 100.0) - (k * -10.0)))) * ((1.0 + (1000.0 / t_0)) / t_0))); end tmp_2 = a_s * tmp; end
a\_m = N[Abs[a], $MachinePrecision]
a\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[a]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[a$95$s_, a$95$m_, k_, m_] := Block[{t$95$0 = N[(k * N[(k * k), $MachinePrecision]), $MachinePrecision]}, N[(a$95$s * If[LessEqual[m, -22000000000.0], N[(N[(a$95$m * N[(99.0 / N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(k * k), $MachinePrecision]), $MachinePrecision], If[LessEqual[m, 1.2], N[(a$95$m / N[(1.0 + N[(k * N[(k + 10.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(a$95$m / N[(1.0 + N[(N[(N[(k * N[(N[(k * k), $MachinePrecision] + -100.0), $MachinePrecision]), $MachinePrecision] / N[(1.0 / N[(N[(N[(k * k), $MachinePrecision] + 100.0), $MachinePrecision] - N[(k * -10.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[(1.0 + N[(1000.0 / t$95$0), $MachinePrecision]), $MachinePrecision] / t$95$0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]]
\begin{array}{l}
a\_m = \left|a\right|
\\
a\_s = \mathsf{copysign}\left(1, a\right)
\\
\begin{array}{l}
t_0 := k \cdot \left(k \cdot k\right)\\
a\_s \cdot \begin{array}{l}
\mathbf{if}\;m \leq -22000000000:\\
\;\;\;\;\frac{a\_m \cdot \frac{99}{k \cdot k}}{k \cdot k}\\
\mathbf{elif}\;m \leq 1.2:\\
\;\;\;\;\frac{a\_m}{1 + k \cdot \left(k + 10\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{a\_m}{1 + \frac{k \cdot \left(k \cdot k + -100\right)}{\frac{1}{\left(k \cdot k + 100\right) - k \cdot -10}} \cdot \frac{1 + \frac{1000}{t\_0}}{t\_0}}\\
\end{array}
\end{array}
\end{array}
if m < -2.2e10Initial program 100.0%
/-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
+-commutativeN/A
+-lowering-+.f64100.0%
Simplified100.0%
Taylor expanded in m around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f6433.6%
Simplified33.6%
Taylor expanded in k around -inf
/-lowering-/.f64N/A
Simplified61.3%
Taylor expanded in k around 0
*-commutativeN/A
associate-*l/N/A
associate-/l*N/A
metadata-evalN/A
associate-*r/N/A
*-lowering-*.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6468.7%
Simplified68.7%
if -2.2e10 < m < 1.19999999999999996Initial program 99.8%
/-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
+-commutativeN/A
+-lowering-+.f6499.9%
Simplified99.9%
Taylor expanded in m around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f6498.2%
Simplified98.2%
if 1.19999999999999996 < m Initial program 76.5%
/-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
+-commutativeN/A
+-lowering-+.f6476.5%
Simplified76.5%
Taylor expanded in m around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f642.8%
Simplified2.8%
metadata-evalN/A
sub-negN/A
flip--N/A
sub-negN/A
metadata-evalN/A
metadata-evalN/A
associate-*r/N/A
div-invN/A
*-commutativeN/A
associate-*l/N/A
flip3-+N/A
div-invN/A
times-fracN/A
*-lowering-*.f64N/A
Applied egg-rr2.3%
Taylor expanded in k around inf
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6447.7%
Simplified47.7%
Final simplification72.1%
a\_m = (fabs.f64 a)
a\_s = (copysign.f64 #s(literal 1 binary64) a)
(FPCore (a_s a_m k m)
:precision binary64
(*
a_s
(if (<= m -22000000000.0)
(/ (* a_m (/ 99.0 (* k k))) (* k k))
(if (<= m 170000000000.0)
(/ a_m (+ 1.0 (* k (+ k 10.0))))
(/
a_m
(+
1.0
(*
(/ (* k (+ (* k k) -100.0)) (/ 1.0 (- (+ (* k k) 100.0) (* k -10.0))))
(/ 1.0 (* k (* k k))))))))))a\_m = fabs(a);
a\_s = copysign(1.0, a);
double code(double a_s, double a_m, double k, double m) {
double tmp;
if (m <= -22000000000.0) {
tmp = (a_m * (99.0 / (k * k))) / (k * k);
} else if (m <= 170000000000.0) {
tmp = a_m / (1.0 + (k * (k + 10.0)));
} else {
tmp = a_m / (1.0 + (((k * ((k * k) + -100.0)) / (1.0 / (((k * k) + 100.0) - (k * -10.0)))) * (1.0 / (k * (k * k)))));
}
return a_s * tmp;
}
a\_m = abs(a)
a\_s = copysign(1.0d0, a)
real(8) function code(a_s, a_m, k, m)
real(8), intent (in) :: a_s
real(8), intent (in) :: a_m
real(8), intent (in) :: k
real(8), intent (in) :: m
real(8) :: tmp
if (m <= (-22000000000.0d0)) then
tmp = (a_m * (99.0d0 / (k * k))) / (k * k)
else if (m <= 170000000000.0d0) then
tmp = a_m / (1.0d0 + (k * (k + 10.0d0)))
else
tmp = a_m / (1.0d0 + (((k * ((k * k) + (-100.0d0))) / (1.0d0 / (((k * k) + 100.0d0) - (k * (-10.0d0))))) * (1.0d0 / (k * (k * k)))))
end if
code = a_s * tmp
end function
a\_m = Math.abs(a);
a\_s = Math.copySign(1.0, a);
public static double code(double a_s, double a_m, double k, double m) {
double tmp;
if (m <= -22000000000.0) {
tmp = (a_m * (99.0 / (k * k))) / (k * k);
} else if (m <= 170000000000.0) {
tmp = a_m / (1.0 + (k * (k + 10.0)));
} else {
tmp = a_m / (1.0 + (((k * ((k * k) + -100.0)) / (1.0 / (((k * k) + 100.0) - (k * -10.0)))) * (1.0 / (k * (k * k)))));
}
return a_s * tmp;
}
a\_m = math.fabs(a) a\_s = math.copysign(1.0, a) def code(a_s, a_m, k, m): tmp = 0 if m <= -22000000000.0: tmp = (a_m * (99.0 / (k * k))) / (k * k) elif m <= 170000000000.0: tmp = a_m / (1.0 + (k * (k + 10.0))) else: tmp = a_m / (1.0 + (((k * ((k * k) + -100.0)) / (1.0 / (((k * k) + 100.0) - (k * -10.0)))) * (1.0 / (k * (k * k))))) return a_s * tmp
a\_m = abs(a) a\_s = copysign(1.0, a) function code(a_s, a_m, k, m) tmp = 0.0 if (m <= -22000000000.0) tmp = Float64(Float64(a_m * Float64(99.0 / Float64(k * k))) / Float64(k * k)); elseif (m <= 170000000000.0) tmp = Float64(a_m / Float64(1.0 + Float64(k * Float64(k + 10.0)))); else tmp = Float64(a_m / Float64(1.0 + Float64(Float64(Float64(k * Float64(Float64(k * k) + -100.0)) / Float64(1.0 / Float64(Float64(Float64(k * k) + 100.0) - Float64(k * -10.0)))) * Float64(1.0 / Float64(k * Float64(k * k)))))); end return Float64(a_s * tmp) end
a\_m = abs(a); a\_s = sign(a) * abs(1.0); function tmp_2 = code(a_s, a_m, k, m) tmp = 0.0; if (m <= -22000000000.0) tmp = (a_m * (99.0 / (k * k))) / (k * k); elseif (m <= 170000000000.0) tmp = a_m / (1.0 + (k * (k + 10.0))); else tmp = a_m / (1.0 + (((k * ((k * k) + -100.0)) / (1.0 / (((k * k) + 100.0) - (k * -10.0)))) * (1.0 / (k * (k * k))))); end tmp_2 = a_s * tmp; end
a\_m = N[Abs[a], $MachinePrecision]
a\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[a]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[a$95$s_, a$95$m_, k_, m_] := N[(a$95$s * If[LessEqual[m, -22000000000.0], N[(N[(a$95$m * N[(99.0 / N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(k * k), $MachinePrecision]), $MachinePrecision], If[LessEqual[m, 170000000000.0], N[(a$95$m / N[(1.0 + N[(k * N[(k + 10.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(a$95$m / N[(1.0 + N[(N[(N[(k * N[(N[(k * k), $MachinePrecision] + -100.0), $MachinePrecision]), $MachinePrecision] / N[(1.0 / N[(N[(N[(k * k), $MachinePrecision] + 100.0), $MachinePrecision] - N[(k * -10.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(1.0 / N[(k * N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]
\begin{array}{l}
a\_m = \left|a\right|
\\
a\_s = \mathsf{copysign}\left(1, a\right)
\\
a\_s \cdot \begin{array}{l}
\mathbf{if}\;m \leq -22000000000:\\
\;\;\;\;\frac{a\_m \cdot \frac{99}{k \cdot k}}{k \cdot k}\\
\mathbf{elif}\;m \leq 170000000000:\\
\;\;\;\;\frac{a\_m}{1 + k \cdot \left(k + 10\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{a\_m}{1 + \frac{k \cdot \left(k \cdot k + -100\right)}{\frac{1}{\left(k \cdot k + 100\right) - k \cdot -10}} \cdot \frac{1}{k \cdot \left(k \cdot k\right)}}\\
\end{array}
\end{array}
if m < -2.2e10Initial program 100.0%
/-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
+-commutativeN/A
+-lowering-+.f64100.0%
Simplified100.0%
Taylor expanded in m around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f6433.6%
Simplified33.6%
Taylor expanded in k around -inf
/-lowering-/.f64N/A
Simplified61.3%
Taylor expanded in k around 0
*-commutativeN/A
associate-*l/N/A
associate-/l*N/A
metadata-evalN/A
associate-*r/N/A
*-lowering-*.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6468.7%
Simplified68.7%
if -2.2e10 < m < 1.7e11Initial program 98.8%
/-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
+-commutativeN/A
+-lowering-+.f6498.8%
Simplified98.8%
Taylor expanded in m around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f6496.1%
Simplified96.1%
if 1.7e11 < m Initial program 77.1%
/-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
+-commutativeN/A
+-lowering-+.f6477.1%
Simplified77.1%
Taylor expanded in m around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f642.8%
Simplified2.8%
metadata-evalN/A
sub-negN/A
flip--N/A
sub-negN/A
metadata-evalN/A
metadata-evalN/A
associate-*r/N/A
div-invN/A
*-commutativeN/A
associate-*l/N/A
flip3-+N/A
div-invN/A
times-fracN/A
*-lowering-*.f64N/A
Applied egg-rr2.3%
Taylor expanded in k around inf
/-lowering-/.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6433.7%
Simplified33.7%
Final simplification67.2%
a\_m = (fabs.f64 a)
a\_s = (copysign.f64 #s(literal 1 binary64) a)
(FPCore (a_s a_m k m)
:precision binary64
(*
a_s
(if (<= m -22000000000.0)
(/ (* a_m (/ 99.0 (* k k))) (* k k))
(if (<= m 1.18e-5)
(/ a_m (+ 1.0 (* k (+ k 10.0))))
(+ a_m (* k (+ (* a_m -10.0) (* a_m (* k 99.0)))))))))a\_m = fabs(a);
a\_s = copysign(1.0, a);
double code(double a_s, double a_m, double k, double m) {
double tmp;
if (m <= -22000000000.0) {
tmp = (a_m * (99.0 / (k * k))) / (k * k);
} else if (m <= 1.18e-5) {
tmp = a_m / (1.0 + (k * (k + 10.0)));
} else {
tmp = a_m + (k * ((a_m * -10.0) + (a_m * (k * 99.0))));
}
return a_s * tmp;
}
a\_m = abs(a)
a\_s = copysign(1.0d0, a)
real(8) function code(a_s, a_m, k, m)
real(8), intent (in) :: a_s
real(8), intent (in) :: a_m
real(8), intent (in) :: k
real(8), intent (in) :: m
real(8) :: tmp
if (m <= (-22000000000.0d0)) then
tmp = (a_m * (99.0d0 / (k * k))) / (k * k)
else if (m <= 1.18d-5) then
tmp = a_m / (1.0d0 + (k * (k + 10.0d0)))
else
tmp = a_m + (k * ((a_m * (-10.0d0)) + (a_m * (k * 99.0d0))))
end if
code = a_s * tmp
end function
a\_m = Math.abs(a);
a\_s = Math.copySign(1.0, a);
public static double code(double a_s, double a_m, double k, double m) {
double tmp;
if (m <= -22000000000.0) {
tmp = (a_m * (99.0 / (k * k))) / (k * k);
} else if (m <= 1.18e-5) {
tmp = a_m / (1.0 + (k * (k + 10.0)));
} else {
tmp = a_m + (k * ((a_m * -10.0) + (a_m * (k * 99.0))));
}
return a_s * tmp;
}
a\_m = math.fabs(a) a\_s = math.copysign(1.0, a) def code(a_s, a_m, k, m): tmp = 0 if m <= -22000000000.0: tmp = (a_m * (99.0 / (k * k))) / (k * k) elif m <= 1.18e-5: tmp = a_m / (1.0 + (k * (k + 10.0))) else: tmp = a_m + (k * ((a_m * -10.0) + (a_m * (k * 99.0)))) return a_s * tmp
a\_m = abs(a) a\_s = copysign(1.0, a) function code(a_s, a_m, k, m) tmp = 0.0 if (m <= -22000000000.0) tmp = Float64(Float64(a_m * Float64(99.0 / Float64(k * k))) / Float64(k * k)); elseif (m <= 1.18e-5) tmp = Float64(a_m / Float64(1.0 + Float64(k * Float64(k + 10.0)))); else tmp = Float64(a_m + Float64(k * Float64(Float64(a_m * -10.0) + Float64(a_m * Float64(k * 99.0))))); end return Float64(a_s * tmp) end
a\_m = abs(a); a\_s = sign(a) * abs(1.0); function tmp_2 = code(a_s, a_m, k, m) tmp = 0.0; if (m <= -22000000000.0) tmp = (a_m * (99.0 / (k * k))) / (k * k); elseif (m <= 1.18e-5) tmp = a_m / (1.0 + (k * (k + 10.0))); else tmp = a_m + (k * ((a_m * -10.0) + (a_m * (k * 99.0)))); end tmp_2 = a_s * tmp; end
a\_m = N[Abs[a], $MachinePrecision]
a\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[a]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[a$95$s_, a$95$m_, k_, m_] := N[(a$95$s * If[LessEqual[m, -22000000000.0], N[(N[(a$95$m * N[(99.0 / N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(k * k), $MachinePrecision]), $MachinePrecision], If[LessEqual[m, 1.18e-5], N[(a$95$m / N[(1.0 + N[(k * N[(k + 10.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(a$95$m + N[(k * N[(N[(a$95$m * -10.0), $MachinePrecision] + N[(a$95$m * N[(k * 99.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]
\begin{array}{l}
a\_m = \left|a\right|
\\
a\_s = \mathsf{copysign}\left(1, a\right)
\\
a\_s \cdot \begin{array}{l}
\mathbf{if}\;m \leq -22000000000:\\
\;\;\;\;\frac{a\_m \cdot \frac{99}{k \cdot k}}{k \cdot k}\\
\mathbf{elif}\;m \leq 1.18 \cdot 10^{-5}:\\
\;\;\;\;\frac{a\_m}{1 + k \cdot \left(k + 10\right)}\\
\mathbf{else}:\\
\;\;\;\;a\_m + k \cdot \left(a\_m \cdot -10 + a\_m \cdot \left(k \cdot 99\right)\right)\\
\end{array}
\end{array}
if m < -2.2e10Initial program 100.0%
/-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
+-commutativeN/A
+-lowering-+.f64100.0%
Simplified100.0%
Taylor expanded in m around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f6433.6%
Simplified33.6%
Taylor expanded in k around -inf
/-lowering-/.f64N/A
Simplified61.3%
Taylor expanded in k around 0
*-commutativeN/A
associate-*l/N/A
associate-/l*N/A
metadata-evalN/A
associate-*r/N/A
*-lowering-*.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6468.7%
Simplified68.7%
if -2.2e10 < m < 1.18000000000000005e-5Initial program 99.8%
/-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
+-commutativeN/A
+-lowering-+.f6499.8%
Simplified99.8%
Taylor expanded in m around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f6498.9%
Simplified98.9%
if 1.18000000000000005e-5 < m Initial program 76.7%
/-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
+-commutativeN/A
+-lowering-+.f6476.7%
Simplified76.7%
Taylor expanded in m around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f643.2%
Simplified3.2%
+-commutativeN/A
fma-defineN/A
metadata-evalN/A
sub-negN/A
fma-defineN/A
flip--N/A
sub-negN/A
metadata-evalN/A
metadata-evalN/A
associate-*r/N/A
clear-numN/A
flip-+N/A
/-lowering-/.f64N/A
Applied egg-rr2.7%
Taylor expanded in k around 0
distribute-rgt-inN/A
+-commutativeN/A
sub-negN/A
mul-1-negN/A
distribute-neg-inN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
distribute-rgt-inN/A
metadata-evalN/A
cancel-sign-sub-invN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
cancel-sign-sub-invN/A
Simplified26.6%
Final simplification65.1%
a\_m = (fabs.f64 a)
a\_s = (copysign.f64 #s(literal 1 binary64) a)
(FPCore (a_s a_m k m)
:precision binary64
(*
a_s
(if (<= m -22000000000.0)
(/ (* a_m (/ 99.0 (* k k))) (* k k))
(if (<= m 1.55e+44)
(/ a_m (+ 1.0 (* k (+ k 10.0))))
(+ a_m (* a_m (* k -10.0)))))))a\_m = fabs(a);
a\_s = copysign(1.0, a);
double code(double a_s, double a_m, double k, double m) {
double tmp;
if (m <= -22000000000.0) {
tmp = (a_m * (99.0 / (k * k))) / (k * k);
} else if (m <= 1.55e+44) {
tmp = a_m / (1.0 + (k * (k + 10.0)));
} else {
tmp = a_m + (a_m * (k * -10.0));
}
return a_s * tmp;
}
a\_m = abs(a)
a\_s = copysign(1.0d0, a)
real(8) function code(a_s, a_m, k, m)
real(8), intent (in) :: a_s
real(8), intent (in) :: a_m
real(8), intent (in) :: k
real(8), intent (in) :: m
real(8) :: tmp
if (m <= (-22000000000.0d0)) then
tmp = (a_m * (99.0d0 / (k * k))) / (k * k)
else if (m <= 1.55d+44) then
tmp = a_m / (1.0d0 + (k * (k + 10.0d0)))
else
tmp = a_m + (a_m * (k * (-10.0d0)))
end if
code = a_s * tmp
end function
a\_m = Math.abs(a);
a\_s = Math.copySign(1.0, a);
public static double code(double a_s, double a_m, double k, double m) {
double tmp;
if (m <= -22000000000.0) {
tmp = (a_m * (99.0 / (k * k))) / (k * k);
} else if (m <= 1.55e+44) {
tmp = a_m / (1.0 + (k * (k + 10.0)));
} else {
tmp = a_m + (a_m * (k * -10.0));
}
return a_s * tmp;
}
a\_m = math.fabs(a) a\_s = math.copysign(1.0, a) def code(a_s, a_m, k, m): tmp = 0 if m <= -22000000000.0: tmp = (a_m * (99.0 / (k * k))) / (k * k) elif m <= 1.55e+44: tmp = a_m / (1.0 + (k * (k + 10.0))) else: tmp = a_m + (a_m * (k * -10.0)) return a_s * tmp
a\_m = abs(a) a\_s = copysign(1.0, a) function code(a_s, a_m, k, m) tmp = 0.0 if (m <= -22000000000.0) tmp = Float64(Float64(a_m * Float64(99.0 / Float64(k * k))) / Float64(k * k)); elseif (m <= 1.55e+44) tmp = Float64(a_m / Float64(1.0 + Float64(k * Float64(k + 10.0)))); else tmp = Float64(a_m + Float64(a_m * Float64(k * -10.0))); end return Float64(a_s * tmp) end
a\_m = abs(a); a\_s = sign(a) * abs(1.0); function tmp_2 = code(a_s, a_m, k, m) tmp = 0.0; if (m <= -22000000000.0) tmp = (a_m * (99.0 / (k * k))) / (k * k); elseif (m <= 1.55e+44) tmp = a_m / (1.0 + (k * (k + 10.0))); else tmp = a_m + (a_m * (k * -10.0)); end tmp_2 = a_s * tmp; end
a\_m = N[Abs[a], $MachinePrecision]
a\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[a]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[a$95$s_, a$95$m_, k_, m_] := N[(a$95$s * If[LessEqual[m, -22000000000.0], N[(N[(a$95$m * N[(99.0 / N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(k * k), $MachinePrecision]), $MachinePrecision], If[LessEqual[m, 1.55e+44], N[(a$95$m / N[(1.0 + N[(k * N[(k + 10.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(a$95$m + N[(a$95$m * N[(k * -10.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]
\begin{array}{l}
a\_m = \left|a\right|
\\
a\_s = \mathsf{copysign}\left(1, a\right)
\\
a\_s \cdot \begin{array}{l}
\mathbf{if}\;m \leq -22000000000:\\
\;\;\;\;\frac{a\_m \cdot \frac{99}{k \cdot k}}{k \cdot k}\\
\mathbf{elif}\;m \leq 1.55 \cdot 10^{+44}:\\
\;\;\;\;\frac{a\_m}{1 + k \cdot \left(k + 10\right)}\\
\mathbf{else}:\\
\;\;\;\;a\_m + a\_m \cdot \left(k \cdot -10\right)\\
\end{array}
\end{array}
if m < -2.2e10Initial program 100.0%
/-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
+-commutativeN/A
+-lowering-+.f64100.0%
Simplified100.0%
Taylor expanded in m around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f6433.6%
Simplified33.6%
Taylor expanded in k around -inf
/-lowering-/.f64N/A
Simplified61.3%
Taylor expanded in k around 0
*-commutativeN/A
associate-*l/N/A
associate-/l*N/A
metadata-evalN/A
associate-*r/N/A
*-lowering-*.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6468.7%
Simplified68.7%
if -2.2e10 < m < 1.54999999999999998e44Initial program 97.8%
/-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
+-commutativeN/A
+-lowering-+.f6497.8%
Simplified97.8%
Taylor expanded in m around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f6490.4%
Simplified90.4%
if 1.54999999999999998e44 < m Initial program 76.6%
/-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
+-commutativeN/A
+-lowering-+.f6476.6%
Simplified76.6%
Taylor expanded in m around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f642.8%
Simplified2.8%
+-commutativeN/A
fma-defineN/A
metadata-evalN/A
sub-negN/A
fma-defineN/A
flip--N/A
sub-negN/A
metadata-evalN/A
metadata-evalN/A
associate-*r/N/A
clear-numN/A
flip-+N/A
/-lowering-/.f64N/A
Applied egg-rr2.4%
Taylor expanded in k around 0
associate-*r*N/A
+-lowering-+.f64N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6411.7%
Simplified11.7%
a\_m = (fabs.f64 a)
a\_s = (copysign.f64 #s(literal 1 binary64) a)
(FPCore (a_s a_m k m)
:precision binary64
(*
a_s
(if (<= m -22000000000.0)
(/ a_m (* k k))
(if (<= m 1.38e+42)
(/ a_m (+ 1.0 (* k (+ k 10.0))))
(+ a_m (* a_m (* k -10.0)))))))a\_m = fabs(a);
a\_s = copysign(1.0, a);
double code(double a_s, double a_m, double k, double m) {
double tmp;
if (m <= -22000000000.0) {
tmp = a_m / (k * k);
} else if (m <= 1.38e+42) {
tmp = a_m / (1.0 + (k * (k + 10.0)));
} else {
tmp = a_m + (a_m * (k * -10.0));
}
return a_s * tmp;
}
a\_m = abs(a)
a\_s = copysign(1.0d0, a)
real(8) function code(a_s, a_m, k, m)
real(8), intent (in) :: a_s
real(8), intent (in) :: a_m
real(8), intent (in) :: k
real(8), intent (in) :: m
real(8) :: tmp
if (m <= (-22000000000.0d0)) then
tmp = a_m / (k * k)
else if (m <= 1.38d+42) then
tmp = a_m / (1.0d0 + (k * (k + 10.0d0)))
else
tmp = a_m + (a_m * (k * (-10.0d0)))
end if
code = a_s * tmp
end function
a\_m = Math.abs(a);
a\_s = Math.copySign(1.0, a);
public static double code(double a_s, double a_m, double k, double m) {
double tmp;
if (m <= -22000000000.0) {
tmp = a_m / (k * k);
} else if (m <= 1.38e+42) {
tmp = a_m / (1.0 + (k * (k + 10.0)));
} else {
tmp = a_m + (a_m * (k * -10.0));
}
return a_s * tmp;
}
a\_m = math.fabs(a) a\_s = math.copysign(1.0, a) def code(a_s, a_m, k, m): tmp = 0 if m <= -22000000000.0: tmp = a_m / (k * k) elif m <= 1.38e+42: tmp = a_m / (1.0 + (k * (k + 10.0))) else: tmp = a_m + (a_m * (k * -10.0)) return a_s * tmp
a\_m = abs(a) a\_s = copysign(1.0, a) function code(a_s, a_m, k, m) tmp = 0.0 if (m <= -22000000000.0) tmp = Float64(a_m / Float64(k * k)); elseif (m <= 1.38e+42) tmp = Float64(a_m / Float64(1.0 + Float64(k * Float64(k + 10.0)))); else tmp = Float64(a_m + Float64(a_m * Float64(k * -10.0))); end return Float64(a_s * tmp) end
a\_m = abs(a); a\_s = sign(a) * abs(1.0); function tmp_2 = code(a_s, a_m, k, m) tmp = 0.0; if (m <= -22000000000.0) tmp = a_m / (k * k); elseif (m <= 1.38e+42) tmp = a_m / (1.0 + (k * (k + 10.0))); else tmp = a_m + (a_m * (k * -10.0)); end tmp_2 = a_s * tmp; end
a\_m = N[Abs[a], $MachinePrecision]
a\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[a]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[a$95$s_, a$95$m_, k_, m_] := N[(a$95$s * If[LessEqual[m, -22000000000.0], N[(a$95$m / N[(k * k), $MachinePrecision]), $MachinePrecision], If[LessEqual[m, 1.38e+42], N[(a$95$m / N[(1.0 + N[(k * N[(k + 10.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(a$95$m + N[(a$95$m * N[(k * -10.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]
\begin{array}{l}
a\_m = \left|a\right|
\\
a\_s = \mathsf{copysign}\left(1, a\right)
\\
a\_s \cdot \begin{array}{l}
\mathbf{if}\;m \leq -22000000000:\\
\;\;\;\;\frac{a\_m}{k \cdot k}\\
\mathbf{elif}\;m \leq 1.38 \cdot 10^{+42}:\\
\;\;\;\;\frac{a\_m}{1 + k \cdot \left(k + 10\right)}\\
\mathbf{else}:\\
\;\;\;\;a\_m + a\_m \cdot \left(k \cdot -10\right)\\
\end{array}
\end{array}
if m < -2.2e10Initial program 100.0%
/-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
+-commutativeN/A
+-lowering-+.f64100.0%
Simplified100.0%
Taylor expanded in m around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f6433.6%
Simplified33.6%
Taylor expanded in k around inf
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6455.2%
Simplified55.2%
if -2.2e10 < m < 1.3800000000000001e42Initial program 97.8%
/-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
+-commutativeN/A
+-lowering-+.f6497.8%
Simplified97.8%
Taylor expanded in m around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f6491.3%
Simplified91.3%
if 1.3800000000000001e42 < m Initial program 76.9%
/-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
+-commutativeN/A
+-lowering-+.f6476.9%
Simplified76.9%
Taylor expanded in m around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f642.8%
Simplified2.8%
+-commutativeN/A
fma-defineN/A
metadata-evalN/A
sub-negN/A
fma-defineN/A
flip--N/A
sub-negN/A
metadata-evalN/A
metadata-evalN/A
associate-*r/N/A
clear-numN/A
flip-+N/A
/-lowering-/.f64N/A
Applied egg-rr2.4%
Taylor expanded in k around 0
associate-*r*N/A
+-lowering-+.f64N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6411.6%
Simplified11.6%
a\_m = (fabs.f64 a)
a\_s = (copysign.f64 #s(literal 1 binary64) a)
(FPCore (a_s a_m k m)
:precision binary64
(*
a_s
(if (<= m -22000000000.0)
(/ a_m (* k k))
(if (<= m 2.2e+44) (/ a_m (+ 1.0 (* k k))) (+ a_m (* a_m (* k -10.0)))))))a\_m = fabs(a);
a\_s = copysign(1.0, a);
double code(double a_s, double a_m, double k, double m) {
double tmp;
if (m <= -22000000000.0) {
tmp = a_m / (k * k);
} else if (m <= 2.2e+44) {
tmp = a_m / (1.0 + (k * k));
} else {
tmp = a_m + (a_m * (k * -10.0));
}
return a_s * tmp;
}
a\_m = abs(a)
a\_s = copysign(1.0d0, a)
real(8) function code(a_s, a_m, k, m)
real(8), intent (in) :: a_s
real(8), intent (in) :: a_m
real(8), intent (in) :: k
real(8), intent (in) :: m
real(8) :: tmp
if (m <= (-22000000000.0d0)) then
tmp = a_m / (k * k)
else if (m <= 2.2d+44) then
tmp = a_m / (1.0d0 + (k * k))
else
tmp = a_m + (a_m * (k * (-10.0d0)))
end if
code = a_s * tmp
end function
a\_m = Math.abs(a);
a\_s = Math.copySign(1.0, a);
public static double code(double a_s, double a_m, double k, double m) {
double tmp;
if (m <= -22000000000.0) {
tmp = a_m / (k * k);
} else if (m <= 2.2e+44) {
tmp = a_m / (1.0 + (k * k));
} else {
tmp = a_m + (a_m * (k * -10.0));
}
return a_s * tmp;
}
a\_m = math.fabs(a) a\_s = math.copysign(1.0, a) def code(a_s, a_m, k, m): tmp = 0 if m <= -22000000000.0: tmp = a_m / (k * k) elif m <= 2.2e+44: tmp = a_m / (1.0 + (k * k)) else: tmp = a_m + (a_m * (k * -10.0)) return a_s * tmp
a\_m = abs(a) a\_s = copysign(1.0, a) function code(a_s, a_m, k, m) tmp = 0.0 if (m <= -22000000000.0) tmp = Float64(a_m / Float64(k * k)); elseif (m <= 2.2e+44) tmp = Float64(a_m / Float64(1.0 + Float64(k * k))); else tmp = Float64(a_m + Float64(a_m * Float64(k * -10.0))); end return Float64(a_s * tmp) end
a\_m = abs(a); a\_s = sign(a) * abs(1.0); function tmp_2 = code(a_s, a_m, k, m) tmp = 0.0; if (m <= -22000000000.0) tmp = a_m / (k * k); elseif (m <= 2.2e+44) tmp = a_m / (1.0 + (k * k)); else tmp = a_m + (a_m * (k * -10.0)); end tmp_2 = a_s * tmp; end
a\_m = N[Abs[a], $MachinePrecision]
a\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[a]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[a$95$s_, a$95$m_, k_, m_] := N[(a$95$s * If[LessEqual[m, -22000000000.0], N[(a$95$m / N[(k * k), $MachinePrecision]), $MachinePrecision], If[LessEqual[m, 2.2e+44], N[(a$95$m / N[(1.0 + N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(a$95$m + N[(a$95$m * N[(k * -10.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]
\begin{array}{l}
a\_m = \left|a\right|
\\
a\_s = \mathsf{copysign}\left(1, a\right)
\\
a\_s \cdot \begin{array}{l}
\mathbf{if}\;m \leq -22000000000:\\
\;\;\;\;\frac{a\_m}{k \cdot k}\\
\mathbf{elif}\;m \leq 2.2 \cdot 10^{+44}:\\
\;\;\;\;\frac{a\_m}{1 + k \cdot k}\\
\mathbf{else}:\\
\;\;\;\;a\_m + a\_m \cdot \left(k \cdot -10\right)\\
\end{array}
\end{array}
if m < -2.2e10Initial program 100.0%
/-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
+-commutativeN/A
+-lowering-+.f64100.0%
Simplified100.0%
Taylor expanded in m around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f6433.6%
Simplified33.6%
Taylor expanded in k around inf
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6455.2%
Simplified55.2%
if -2.2e10 < m < 2.19999999999999996e44Initial program 97.8%
/-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
+-commutativeN/A
+-lowering-+.f6497.8%
Simplified97.8%
Taylor expanded in m around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f6490.4%
Simplified90.4%
Taylor expanded in k around inf
unpow2N/A
*-lowering-*.f6488.1%
Simplified88.1%
if 2.19999999999999996e44 < m Initial program 76.6%
/-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
+-commutativeN/A
+-lowering-+.f6476.6%
Simplified76.6%
Taylor expanded in m around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f642.8%
Simplified2.8%
+-commutativeN/A
fma-defineN/A
metadata-evalN/A
sub-negN/A
fma-defineN/A
flip--N/A
sub-negN/A
metadata-evalN/A
metadata-evalN/A
associate-*r/N/A
clear-numN/A
flip-+N/A
/-lowering-/.f64N/A
Applied egg-rr2.4%
Taylor expanded in k around 0
associate-*r*N/A
+-lowering-+.f64N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6411.7%
Simplified11.7%
a\_m = (fabs.f64 a)
a\_s = (copysign.f64 #s(literal 1 binary64) a)
(FPCore (a_s a_m k m)
:precision binary64
(*
a_s
(if (<= m -1.02e-84)
(/ a_m (* k k))
(if (<= m 1.8e+44)
(/ a_m (+ 1.0 (* k 10.0)))
(+ a_m (* a_m (* k -10.0)))))))a\_m = fabs(a);
a\_s = copysign(1.0, a);
double code(double a_s, double a_m, double k, double m) {
double tmp;
if (m <= -1.02e-84) {
tmp = a_m / (k * k);
} else if (m <= 1.8e+44) {
tmp = a_m / (1.0 + (k * 10.0));
} else {
tmp = a_m + (a_m * (k * -10.0));
}
return a_s * tmp;
}
a\_m = abs(a)
a\_s = copysign(1.0d0, a)
real(8) function code(a_s, a_m, k, m)
real(8), intent (in) :: a_s
real(8), intent (in) :: a_m
real(8), intent (in) :: k
real(8), intent (in) :: m
real(8) :: tmp
if (m <= (-1.02d-84)) then
tmp = a_m / (k * k)
else if (m <= 1.8d+44) then
tmp = a_m / (1.0d0 + (k * 10.0d0))
else
tmp = a_m + (a_m * (k * (-10.0d0)))
end if
code = a_s * tmp
end function
a\_m = Math.abs(a);
a\_s = Math.copySign(1.0, a);
public static double code(double a_s, double a_m, double k, double m) {
double tmp;
if (m <= -1.02e-84) {
tmp = a_m / (k * k);
} else if (m <= 1.8e+44) {
tmp = a_m / (1.0 + (k * 10.0));
} else {
tmp = a_m + (a_m * (k * -10.0));
}
return a_s * tmp;
}
a\_m = math.fabs(a) a\_s = math.copysign(1.0, a) def code(a_s, a_m, k, m): tmp = 0 if m <= -1.02e-84: tmp = a_m / (k * k) elif m <= 1.8e+44: tmp = a_m / (1.0 + (k * 10.0)) else: tmp = a_m + (a_m * (k * -10.0)) return a_s * tmp
a\_m = abs(a) a\_s = copysign(1.0, a) function code(a_s, a_m, k, m) tmp = 0.0 if (m <= -1.02e-84) tmp = Float64(a_m / Float64(k * k)); elseif (m <= 1.8e+44) tmp = Float64(a_m / Float64(1.0 + Float64(k * 10.0))); else tmp = Float64(a_m + Float64(a_m * Float64(k * -10.0))); end return Float64(a_s * tmp) end
a\_m = abs(a); a\_s = sign(a) * abs(1.0); function tmp_2 = code(a_s, a_m, k, m) tmp = 0.0; if (m <= -1.02e-84) tmp = a_m / (k * k); elseif (m <= 1.8e+44) tmp = a_m / (1.0 + (k * 10.0)); else tmp = a_m + (a_m * (k * -10.0)); end tmp_2 = a_s * tmp; end
a\_m = N[Abs[a], $MachinePrecision]
a\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[a]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[a$95$s_, a$95$m_, k_, m_] := N[(a$95$s * If[LessEqual[m, -1.02e-84], N[(a$95$m / N[(k * k), $MachinePrecision]), $MachinePrecision], If[LessEqual[m, 1.8e+44], N[(a$95$m / N[(1.0 + N[(k * 10.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(a$95$m + N[(a$95$m * N[(k * -10.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]
\begin{array}{l}
a\_m = \left|a\right|
\\
a\_s = \mathsf{copysign}\left(1, a\right)
\\
a\_s \cdot \begin{array}{l}
\mathbf{if}\;m \leq -1.02 \cdot 10^{-84}:\\
\;\;\;\;\frac{a\_m}{k \cdot k}\\
\mathbf{elif}\;m \leq 1.8 \cdot 10^{+44}:\\
\;\;\;\;\frac{a\_m}{1 + k \cdot 10}\\
\mathbf{else}:\\
\;\;\;\;a\_m + a\_m \cdot \left(k \cdot -10\right)\\
\end{array}
\end{array}
if m < -1.02000000000000004e-84Initial program 100.0%
/-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
+-commutativeN/A
+-lowering-+.f64100.0%
Simplified100.0%
Taylor expanded in m around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f6440.6%
Simplified40.6%
Taylor expanded in k around inf
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6457.7%
Simplified57.7%
if -1.02000000000000004e-84 < m < 1.8e44Initial program 97.6%
/-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
+-commutativeN/A
+-lowering-+.f6497.6%
Simplified97.6%
Taylor expanded in m around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f6489.7%
Simplified89.7%
Taylor expanded in k around 0
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6470.7%
Simplified70.7%
if 1.8e44 < m Initial program 76.6%
/-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
+-commutativeN/A
+-lowering-+.f6476.6%
Simplified76.6%
Taylor expanded in m around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f642.8%
Simplified2.8%
+-commutativeN/A
fma-defineN/A
metadata-evalN/A
sub-negN/A
fma-defineN/A
flip--N/A
sub-negN/A
metadata-evalN/A
metadata-evalN/A
associate-*r/N/A
clear-numN/A
flip-+N/A
/-lowering-/.f64N/A
Applied egg-rr2.4%
Taylor expanded in k around 0
associate-*r*N/A
+-lowering-+.f64N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6411.7%
Simplified11.7%
a\_m = (fabs.f64 a)
a\_s = (copysign.f64 #s(literal 1 binary64) a)
(FPCore (a_s a_m k m)
:precision binary64
(*
a_s
(if (<= k -8.5e-286)
(/ a_m (* k k))
(if (<= k 0.075) (+ a_m (* a_m (* k -10.0))) (/ a_m (* k (+ k 10.0)))))))a\_m = fabs(a);
a\_s = copysign(1.0, a);
double code(double a_s, double a_m, double k, double m) {
double tmp;
if (k <= -8.5e-286) {
tmp = a_m / (k * k);
} else if (k <= 0.075) {
tmp = a_m + (a_m * (k * -10.0));
} else {
tmp = a_m / (k * (k + 10.0));
}
return a_s * tmp;
}
a\_m = abs(a)
a\_s = copysign(1.0d0, a)
real(8) function code(a_s, a_m, k, m)
real(8), intent (in) :: a_s
real(8), intent (in) :: a_m
real(8), intent (in) :: k
real(8), intent (in) :: m
real(8) :: tmp
if (k <= (-8.5d-286)) then
tmp = a_m / (k * k)
else if (k <= 0.075d0) then
tmp = a_m + (a_m * (k * (-10.0d0)))
else
tmp = a_m / (k * (k + 10.0d0))
end if
code = a_s * tmp
end function
a\_m = Math.abs(a);
a\_s = Math.copySign(1.0, a);
public static double code(double a_s, double a_m, double k, double m) {
double tmp;
if (k <= -8.5e-286) {
tmp = a_m / (k * k);
} else if (k <= 0.075) {
tmp = a_m + (a_m * (k * -10.0));
} else {
tmp = a_m / (k * (k + 10.0));
}
return a_s * tmp;
}
a\_m = math.fabs(a) a\_s = math.copysign(1.0, a) def code(a_s, a_m, k, m): tmp = 0 if k <= -8.5e-286: tmp = a_m / (k * k) elif k <= 0.075: tmp = a_m + (a_m * (k * -10.0)) else: tmp = a_m / (k * (k + 10.0)) return a_s * tmp
a\_m = abs(a) a\_s = copysign(1.0, a) function code(a_s, a_m, k, m) tmp = 0.0 if (k <= -8.5e-286) tmp = Float64(a_m / Float64(k * k)); elseif (k <= 0.075) tmp = Float64(a_m + Float64(a_m * Float64(k * -10.0))); else tmp = Float64(a_m / Float64(k * Float64(k + 10.0))); end return Float64(a_s * tmp) end
a\_m = abs(a); a\_s = sign(a) * abs(1.0); function tmp_2 = code(a_s, a_m, k, m) tmp = 0.0; if (k <= -8.5e-286) tmp = a_m / (k * k); elseif (k <= 0.075) tmp = a_m + (a_m * (k * -10.0)); else tmp = a_m / (k * (k + 10.0)); end tmp_2 = a_s * tmp; end
a\_m = N[Abs[a], $MachinePrecision]
a\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[a]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[a$95$s_, a$95$m_, k_, m_] := N[(a$95$s * If[LessEqual[k, -8.5e-286], N[(a$95$m / N[(k * k), $MachinePrecision]), $MachinePrecision], If[LessEqual[k, 0.075], N[(a$95$m + N[(a$95$m * N[(k * -10.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(a$95$m / N[(k * N[(k + 10.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]
\begin{array}{l}
a\_m = \left|a\right|
\\
a\_s = \mathsf{copysign}\left(1, a\right)
\\
a\_s \cdot \begin{array}{l}
\mathbf{if}\;k \leq -8.5 \cdot 10^{-286}:\\
\;\;\;\;\frac{a\_m}{k \cdot k}\\
\mathbf{elif}\;k \leq 0.075:\\
\;\;\;\;a\_m + a\_m \cdot \left(k \cdot -10\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{a\_m}{k \cdot \left(k + 10\right)}\\
\end{array}
\end{array}
if k < -8.4999999999999998e-286Initial program 88.0%
/-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
+-commutativeN/A
+-lowering-+.f6488.0%
Simplified88.0%
Taylor expanded in m around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f6422.6%
Simplified22.6%
Taylor expanded in k around inf
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6429.9%
Simplified29.9%
if -8.4999999999999998e-286 < k < 0.0749999999999999972Initial program 99.9%
/-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
+-commutativeN/A
+-lowering-+.f64100.0%
Simplified100.0%
Taylor expanded in m around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f6451.0%
Simplified51.0%
+-commutativeN/A
fma-defineN/A
metadata-evalN/A
sub-negN/A
fma-defineN/A
flip--N/A
sub-negN/A
metadata-evalN/A
metadata-evalN/A
associate-*r/N/A
clear-numN/A
flip-+N/A
/-lowering-/.f64N/A
Applied egg-rr51.0%
Taylor expanded in k around 0
associate-*r*N/A
+-lowering-+.f64N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6450.3%
Simplified50.3%
if 0.0749999999999999972 < k Initial program 86.5%
/-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
+-commutativeN/A
+-lowering-+.f6486.5%
Simplified86.5%
Taylor expanded in m around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f6461.7%
Simplified61.7%
Taylor expanded in k around inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f6460.7%
Simplified60.7%
Taylor expanded in k around 0
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f6460.7%
Simplified60.7%
a\_m = (fabs.f64 a)
a\_s = (copysign.f64 #s(literal 1 binary64) a)
(FPCore (a_s a_m k m)
:precision binary64
(let* ((t_0 (/ a_m (* k k))))
(*
a_s
(if (<= k -8.5e-286)
t_0
(if (<= k 0.0076) (+ a_m (* a_m (* k -10.0))) t_0)))))a\_m = fabs(a);
a\_s = copysign(1.0, a);
double code(double a_s, double a_m, double k, double m) {
double t_0 = a_m / (k * k);
double tmp;
if (k <= -8.5e-286) {
tmp = t_0;
} else if (k <= 0.0076) {
tmp = a_m + (a_m * (k * -10.0));
} else {
tmp = t_0;
}
return a_s * tmp;
}
a\_m = abs(a)
a\_s = copysign(1.0d0, a)
real(8) function code(a_s, a_m, k, m)
real(8), intent (in) :: a_s
real(8), intent (in) :: a_m
real(8), intent (in) :: k
real(8), intent (in) :: m
real(8) :: t_0
real(8) :: tmp
t_0 = a_m / (k * k)
if (k <= (-8.5d-286)) then
tmp = t_0
else if (k <= 0.0076d0) then
tmp = a_m + (a_m * (k * (-10.0d0)))
else
tmp = t_0
end if
code = a_s * tmp
end function
a\_m = Math.abs(a);
a\_s = Math.copySign(1.0, a);
public static double code(double a_s, double a_m, double k, double m) {
double t_0 = a_m / (k * k);
double tmp;
if (k <= -8.5e-286) {
tmp = t_0;
} else if (k <= 0.0076) {
tmp = a_m + (a_m * (k * -10.0));
} else {
tmp = t_0;
}
return a_s * tmp;
}
a\_m = math.fabs(a) a\_s = math.copysign(1.0, a) def code(a_s, a_m, k, m): t_0 = a_m / (k * k) tmp = 0 if k <= -8.5e-286: tmp = t_0 elif k <= 0.0076: tmp = a_m + (a_m * (k * -10.0)) else: tmp = t_0 return a_s * tmp
a\_m = abs(a) a\_s = copysign(1.0, a) function code(a_s, a_m, k, m) t_0 = Float64(a_m / Float64(k * k)) tmp = 0.0 if (k <= -8.5e-286) tmp = t_0; elseif (k <= 0.0076) tmp = Float64(a_m + Float64(a_m * Float64(k * -10.0))); else tmp = t_0; end return Float64(a_s * tmp) end
a\_m = abs(a); a\_s = sign(a) * abs(1.0); function tmp_2 = code(a_s, a_m, k, m) t_0 = a_m / (k * k); tmp = 0.0; if (k <= -8.5e-286) tmp = t_0; elseif (k <= 0.0076) tmp = a_m + (a_m * (k * -10.0)); else tmp = t_0; end tmp_2 = a_s * tmp; end
a\_m = N[Abs[a], $MachinePrecision]
a\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[a]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[a$95$s_, a$95$m_, k_, m_] := Block[{t$95$0 = N[(a$95$m / N[(k * k), $MachinePrecision]), $MachinePrecision]}, N[(a$95$s * If[LessEqual[k, -8.5e-286], t$95$0, If[LessEqual[k, 0.0076], N[(a$95$m + N[(a$95$m * N[(k * -10.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$0]]), $MachinePrecision]]
\begin{array}{l}
a\_m = \left|a\right|
\\
a\_s = \mathsf{copysign}\left(1, a\right)
\\
\begin{array}{l}
t_0 := \frac{a\_m}{k \cdot k}\\
a\_s \cdot \begin{array}{l}
\mathbf{if}\;k \leq -8.5 \cdot 10^{-286}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;k \leq 0.0076:\\
\;\;\;\;a\_m + a\_m \cdot \left(k \cdot -10\right)\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
\end{array}
if k < -8.4999999999999998e-286 or 0.00759999999999999998 < k Initial program 87.2%
/-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
+-commutativeN/A
+-lowering-+.f6487.2%
Simplified87.2%
Taylor expanded in m around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f6443.0%
Simplified43.0%
Taylor expanded in k around inf
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6445.7%
Simplified45.7%
if -8.4999999999999998e-286 < k < 0.00759999999999999998Initial program 99.9%
/-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
+-commutativeN/A
+-lowering-+.f64100.0%
Simplified100.0%
Taylor expanded in m around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f6451.0%
Simplified51.0%
+-commutativeN/A
fma-defineN/A
metadata-evalN/A
sub-negN/A
fma-defineN/A
flip--N/A
sub-negN/A
metadata-evalN/A
metadata-evalN/A
associate-*r/N/A
clear-numN/A
flip-+N/A
/-lowering-/.f64N/A
Applied egg-rr51.0%
Taylor expanded in k around 0
associate-*r*N/A
+-lowering-+.f64N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6450.3%
Simplified50.3%
a\_m = (fabs.f64 a) a\_s = (copysign.f64 #s(literal 1 binary64) a) (FPCore (a_s a_m k m) :precision binary64 (let* ((t_0 (/ a_m (* k k)))) (* a_s (if (<= k -8.5e-286) t_0 (if (<= k 1.0) a_m t_0)))))
a\_m = fabs(a);
a\_s = copysign(1.0, a);
double code(double a_s, double a_m, double k, double m) {
double t_0 = a_m / (k * k);
double tmp;
if (k <= -8.5e-286) {
tmp = t_0;
} else if (k <= 1.0) {
tmp = a_m;
} else {
tmp = t_0;
}
return a_s * tmp;
}
a\_m = abs(a)
a\_s = copysign(1.0d0, a)
real(8) function code(a_s, a_m, k, m)
real(8), intent (in) :: a_s
real(8), intent (in) :: a_m
real(8), intent (in) :: k
real(8), intent (in) :: m
real(8) :: t_0
real(8) :: tmp
t_0 = a_m / (k * k)
if (k <= (-8.5d-286)) then
tmp = t_0
else if (k <= 1.0d0) then
tmp = a_m
else
tmp = t_0
end if
code = a_s * tmp
end function
a\_m = Math.abs(a);
a\_s = Math.copySign(1.0, a);
public static double code(double a_s, double a_m, double k, double m) {
double t_0 = a_m / (k * k);
double tmp;
if (k <= -8.5e-286) {
tmp = t_0;
} else if (k <= 1.0) {
tmp = a_m;
} else {
tmp = t_0;
}
return a_s * tmp;
}
a\_m = math.fabs(a) a\_s = math.copysign(1.0, a) def code(a_s, a_m, k, m): t_0 = a_m / (k * k) tmp = 0 if k <= -8.5e-286: tmp = t_0 elif k <= 1.0: tmp = a_m else: tmp = t_0 return a_s * tmp
a\_m = abs(a) a\_s = copysign(1.0, a) function code(a_s, a_m, k, m) t_0 = Float64(a_m / Float64(k * k)) tmp = 0.0 if (k <= -8.5e-286) tmp = t_0; elseif (k <= 1.0) tmp = a_m; else tmp = t_0; end return Float64(a_s * tmp) end
a\_m = abs(a); a\_s = sign(a) * abs(1.0); function tmp_2 = code(a_s, a_m, k, m) t_0 = a_m / (k * k); tmp = 0.0; if (k <= -8.5e-286) tmp = t_0; elseif (k <= 1.0) tmp = a_m; else tmp = t_0; end tmp_2 = a_s * tmp; end
a\_m = N[Abs[a], $MachinePrecision]
a\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[a]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[a$95$s_, a$95$m_, k_, m_] := Block[{t$95$0 = N[(a$95$m / N[(k * k), $MachinePrecision]), $MachinePrecision]}, N[(a$95$s * If[LessEqual[k, -8.5e-286], t$95$0, If[LessEqual[k, 1.0], a$95$m, t$95$0]]), $MachinePrecision]]
\begin{array}{l}
a\_m = \left|a\right|
\\
a\_s = \mathsf{copysign}\left(1, a\right)
\\
\begin{array}{l}
t_0 := \frac{a\_m}{k \cdot k}\\
a\_s \cdot \begin{array}{l}
\mathbf{if}\;k \leq -8.5 \cdot 10^{-286}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;k \leq 1:\\
\;\;\;\;a\_m\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
\end{array}
if k < -8.4999999999999998e-286 or 1 < k Initial program 87.1%
/-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
+-commutativeN/A
+-lowering-+.f6487.1%
Simplified87.1%
Taylor expanded in m around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f6442.6%
Simplified42.6%
Taylor expanded in k around inf
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6445.9%
Simplified45.9%
if -8.4999999999999998e-286 < k < 1Initial program 99.9%
/-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
+-commutativeN/A
+-lowering-+.f64100.0%
Simplified100.0%
Taylor expanded in m around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f6451.5%
Simplified51.5%
Taylor expanded in k around 0
Simplified49.5%
a\_m = (fabs.f64 a) a\_s = (copysign.f64 #s(literal 1 binary64) a) (FPCore (a_s a_m k m) :precision binary64 (let* ((t_0 (/ a_m (* k 10.0)))) (* a_s (if (<= k -2.52e+98) t_0 (if (<= k 0.1) a_m t_0)))))
a\_m = fabs(a);
a\_s = copysign(1.0, a);
double code(double a_s, double a_m, double k, double m) {
double t_0 = a_m / (k * 10.0);
double tmp;
if (k <= -2.52e+98) {
tmp = t_0;
} else if (k <= 0.1) {
tmp = a_m;
} else {
tmp = t_0;
}
return a_s * tmp;
}
a\_m = abs(a)
a\_s = copysign(1.0d0, a)
real(8) function code(a_s, a_m, k, m)
real(8), intent (in) :: a_s
real(8), intent (in) :: a_m
real(8), intent (in) :: k
real(8), intent (in) :: m
real(8) :: t_0
real(8) :: tmp
t_0 = a_m / (k * 10.0d0)
if (k <= (-2.52d+98)) then
tmp = t_0
else if (k <= 0.1d0) then
tmp = a_m
else
tmp = t_0
end if
code = a_s * tmp
end function
a\_m = Math.abs(a);
a\_s = Math.copySign(1.0, a);
public static double code(double a_s, double a_m, double k, double m) {
double t_0 = a_m / (k * 10.0);
double tmp;
if (k <= -2.52e+98) {
tmp = t_0;
} else if (k <= 0.1) {
tmp = a_m;
} else {
tmp = t_0;
}
return a_s * tmp;
}
a\_m = math.fabs(a) a\_s = math.copysign(1.0, a) def code(a_s, a_m, k, m): t_0 = a_m / (k * 10.0) tmp = 0 if k <= -2.52e+98: tmp = t_0 elif k <= 0.1: tmp = a_m else: tmp = t_0 return a_s * tmp
a\_m = abs(a) a\_s = copysign(1.0, a) function code(a_s, a_m, k, m) t_0 = Float64(a_m / Float64(k * 10.0)) tmp = 0.0 if (k <= -2.52e+98) tmp = t_0; elseif (k <= 0.1) tmp = a_m; else tmp = t_0; end return Float64(a_s * tmp) end
a\_m = abs(a); a\_s = sign(a) * abs(1.0); function tmp_2 = code(a_s, a_m, k, m) t_0 = a_m / (k * 10.0); tmp = 0.0; if (k <= -2.52e+98) tmp = t_0; elseif (k <= 0.1) tmp = a_m; else tmp = t_0; end tmp_2 = a_s * tmp; end
a\_m = N[Abs[a], $MachinePrecision]
a\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[a]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[a$95$s_, a$95$m_, k_, m_] := Block[{t$95$0 = N[(a$95$m / N[(k * 10.0), $MachinePrecision]), $MachinePrecision]}, N[(a$95$s * If[LessEqual[k, -2.52e+98], t$95$0, If[LessEqual[k, 0.1], a$95$m, t$95$0]]), $MachinePrecision]]
\begin{array}{l}
a\_m = \left|a\right|
\\
a\_s = \mathsf{copysign}\left(1, a\right)
\\
\begin{array}{l}
t_0 := \frac{a\_m}{k \cdot 10}\\
a\_s \cdot \begin{array}{l}
\mathbf{if}\;k \leq -2.52 \cdot 10^{+98}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;k \leq 0.1:\\
\;\;\;\;a\_m\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
\end{array}
if k < -2.51999999999999999e98 or 0.10000000000000001 < k Initial program 82.2%
/-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
+-commutativeN/A
+-lowering-+.f6482.2%
Simplified82.2%
Taylor expanded in m around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f6458.3%
Simplified58.3%
Taylor expanded in k around inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f6457.6%
Simplified57.6%
Taylor expanded in k around 0
*-commutativeN/A
*-lowering-*.f6428.6%
Simplified28.6%
if -2.51999999999999999e98 < k < 0.10000000000000001Initial program 100.0%
/-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
+-commutativeN/A
+-lowering-+.f64100.0%
Simplified100.0%
Taylor expanded in m around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f6436.4%
Simplified36.4%
Taylor expanded in k around 0
Simplified35.7%
a\_m = (fabs.f64 a) a\_s = (copysign.f64 #s(literal 1 binary64) a) (FPCore (a_s a_m k m) :precision binary64 (* a_s a_m))
a\_m = fabs(a);
a\_s = copysign(1.0, a);
double code(double a_s, double a_m, double k, double m) {
return a_s * a_m;
}
a\_m = abs(a)
a\_s = copysign(1.0d0, a)
real(8) function code(a_s, a_m, k, m)
real(8), intent (in) :: a_s
real(8), intent (in) :: a_m
real(8), intent (in) :: k
real(8), intent (in) :: m
code = a_s * a_m
end function
a\_m = Math.abs(a);
a\_s = Math.copySign(1.0, a);
public static double code(double a_s, double a_m, double k, double m) {
return a_s * a_m;
}
a\_m = math.fabs(a) a\_s = math.copysign(1.0, a) def code(a_s, a_m, k, m): return a_s * a_m
a\_m = abs(a) a\_s = copysign(1.0, a) function code(a_s, a_m, k, m) return Float64(a_s * a_m) end
a\_m = abs(a); a\_s = sign(a) * abs(1.0); function tmp = code(a_s, a_m, k, m) tmp = a_s * a_m; end
a\_m = N[Abs[a], $MachinePrecision]
a\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[a]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[a$95$s_, a$95$m_, k_, m_] := N[(a$95$s * a$95$m), $MachinePrecision]
\begin{array}{l}
a\_m = \left|a\right|
\\
a\_s = \mathsf{copysign}\left(1, a\right)
\\
a\_s \cdot a\_m
\end{array}
Initial program 92.1%
/-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
+-commutativeN/A
+-lowering-+.f6492.1%
Simplified92.1%
Taylor expanded in m around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f6446.1%
Simplified46.1%
Taylor expanded in k around 0
Simplified21.9%
herbie shell --seed 2024139
(FPCore (a k m)
:name "Falkner and Boettcher, Appendix A"
:precision binary64
(/ (* a (pow k m)) (+ (+ 1.0 (* 10.0 k)) (* k k))))