Falkner and Boettcher, Appendix A

Percentage Accurate: 89.9% → 98.3%

Time: 12.1s

Alternatives: 15

Speedup: 1.1×

• 22.9% of points produce a very large (infinite) output. You may want to add a precondition.

Specification

Math

\[\begin{array}{l} \\ \frac{a \cdot {k}^{m}}{\left(1 + 10 \cdot k\right) + k \cdot k} \end{array} \]

FPCore

(FPCore (a k m)
 :precision binary64
 (/ (* a (pow k m)) (+ (+ 1.0 (* 10.0 k)) (* k k))))

C

double code(double a, double k, double m) {
	return (a * pow(k, m)) / ((1.0 + (10.0 * k)) + (k * k));
}

Fortran

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

Java

public static double code(double a, double k, double m) {
	return (a * Math.pow(k, m)) / ((1.0 + (10.0 * k)) + (k * k));
}

Python

def code(a, k, m):
	return (a * math.pow(k, m)) / ((1.0 + (10.0 * k)) + (k * k))

Julia

function code(a, k, m)
	return Float64(Float64(a * (k ^ m)) / Float64(Float64(1.0 + Float64(10.0 * k)) + Float64(k * k)))
end

MATLAB

function tmp = code(a, k, m)
	tmp = (a * (k ^ m)) / ((1.0 + (10.0 * k)) + (k * k));
end

Wolfram

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]

Initial Program: 89.9% accurate, 1.0× speedup

Math

\[\begin{array}{l} \\ \frac{a \cdot {k}^{m}}{\left(1 + 10 \cdot k\right) + k \cdot k} \end{array} \]

FPCore

(FPCore (a k m)
 :precision binary64
 (/ (* a (pow k m)) (+ (+ 1.0 (* 10.0 k)) (* k k))))

C

double code(double a, double k, double m) {
	return (a * pow(k, m)) / ((1.0 + (10.0 * k)) + (k * k));
}

Fortran

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

Java

public static double code(double a, double k, double m) {
	return (a * Math.pow(k, m)) / ((1.0 + (10.0 * k)) + (k * k));
}

Python

def code(a, k, m):
	return (a * math.pow(k, m)) / ((1.0 + (10.0 * k)) + (k * k))

Julia

function code(a, k, m)
	return Float64(Float64(a * (k ^ m)) / Float64(Float64(1.0 + Float64(10.0 * k)) + Float64(k * k)))
end

MATLAB

function tmp = code(a, k, m)
	tmp = (a * (k ^ m)) / ((1.0 + (10.0 * k)) + (k * k));
end

Wolfram

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]

Alternative 1: 98.3% accurate, 0.3× speedup

Math

\[\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 0.001:\\ \;\;\;\;a \cdot {k}^{m}\\ \mathbf{else}:\\ \;\;\;\;e^{\log \left(a \cdot \frac{1}{{k}^{\left(-m\right)}}\right) - 2 \cdot \log k}\\ \end{array} \end{array} \]

FPCore

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 0.001)
    (* a (pow k m))
    (exp (- (log (* a (/ 1.0 (pow k (- m))))) (* 2.0 (log k)))))))

C

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 <= 0.001) {
		tmp = a * pow(k, m);
	} else {
		tmp = exp((log((a * (1.0 / pow(k, -m)))) - (2.0 * log(k))));
	}
	return a_s * tmp;
}

Fortran

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 <= 0.001d0) then
        tmp = a * (k ** m)
    else
        tmp = exp((log((a * (1.0d0 / (k ** -m)))) - (2.0d0 * log(k))))
    end if
    code = a_s * tmp
end function

Java

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 <= 0.001) {
		tmp = a * Math.pow(k, m);
	} else {
		tmp = Math.exp((Math.log((a * (1.0 / Math.pow(k, -m)))) - (2.0 * Math.log(k))));
	}
	return a_s * tmp;
}

Python

a\_m = math.fabs(a)
a\_s = math.copysign(1.0, a)
def code(a_s, a_m, k, m):
	tmp = 0
	if k <= 0.001:
		tmp = a * math.pow(k, m)
	else:
		tmp = math.exp((math.log((a * (1.0 / math.pow(k, -m)))) - (2.0 * math.log(k))))
	return a_s * tmp

Julia

a\_m = abs(a)
a\_s = copysign(1.0, a)
function code(a_s, a_m, k, m)
	tmp = 0.0
	if (k <= 0.001)
		tmp = Float64(a * (k ^ m));
	else
		tmp = exp(Float64(log(Float64(a * Float64(1.0 / (k ^ Float64(-m))))) - Float64(2.0 * log(k))));
	end
	return Float64(a_s * tmp)
end

MATLAB

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 <= 0.001)
		tmp = a * (k ^ m);
	else
		tmp = exp((log((a * (1.0 / (k ^ -m)))) - (2.0 * log(k))));
	end
	tmp_2 = a_s * tmp;
end

Wolfram

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, 0.001], N[(a * N[Power[k, m], $MachinePrecision]), $MachinePrecision], N[Exp[N[(N[Log[N[(a * N[(1.0 / N[Power[k, (-m)], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] - N[(2.0 * N[Log[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]), $MachinePrecision]

Derivation

1. Split input into 2 regimes

2. if k < 1e-3

   1. Initial program 94.8%

      \[\frac{a \cdot {k}^{m}}{\left(1 + 10 \cdot k\right) + k \cdot k} \]
   2. Add Preprocessing

   3. Taylor expanded in k around 0

      \[\leadsto \color{blue}{a \cdot {k}^{m}} \]
   4. Step-by-step derivation

      1. lower-*.f64 N/A

         \[\leadsto \color{blue}{a \cdot {k}^{m}} \]

      2. lower-pow.f64 99.5

         \[\leadsto a \cdot \color{blue}{{k}^{m}} \]

   5. Applied rewrites 99.5%

      \[\leadsto \color{blue}{a \cdot {k}^{m}} \]

   if 1e-3 < k

   1. Initial program 79.2%

      \[\frac{a \cdot {k}^{m}}{\left(1 + 10 \cdot k\right) + k \cdot k} \]
   2. Add Preprocessing

   3. Taylor expanded in k around 0

      \[\leadsto \frac{a \cdot {k}^{m}}{\color{blue}{1 + 10 \cdot k}} \]
   4. Step-by-step derivation

      1. +-commutative N/A

         \[\leadsto \frac{a \cdot {k}^{m}}{\color{blue}{10 \cdot k + 1}} \]

      2. lower-fma.f64 61.7

         \[\leadsto \frac{a \cdot {k}^{m}}{\color{blue}{\mathsf{fma}\left(10, k, 1\right)}} \]

   5. Applied rewrites 61.7%

      \[\leadsto \frac{a \cdot {k}^{m}}{\color{blue}{\mathsf{fma}\left(10, k, 1\right)}} \]
   6. Step-by-step derivation

      1. lift-/.f64 N/A

         \[\leadsto \color{blue}{\frac{a \cdot {k}^{m}}{\mathsf{fma}\left(10, k, 1\right)}} \]

      2. clear-num N/A

         \[\leadsto \color{blue}{\frac{1}{\frac{\mathsf{fma}\left(10, k, 1\right)}{a \cdot {k}^{m}}}} \]

      3. inv-pow N/A

         \[\leadsto \color{blue}{{\left(\frac{\mathsf{fma}\left(10, k, 1\right)}{a \cdot {k}^{m}}\right)}^{-1}} \]

      4. pow-to-exp N/A

         \[\leadsto \color{blue}{e^{\log \left(\frac{\mathsf{fma}\left(10, k, 1\right)}{a \cdot {k}^{m}}\right) \cdot -1}} \]

      5. lower-exp.f64 N/A

         \[\leadsto \color{blue}{e^{\log \left(\frac{\mathsf{fma}\left(10, k, 1\right)}{a \cdot {k}^{m}}\right) \cdot -1}} \]

      6. lower-*.f64 N/A

         \[\leadsto e^{\color{blue}{\log \left(\frac{\mathsf{fma}\left(10, k, 1\right)}{a \cdot {k}^{m}}\right) \cdot -1}} \]

      7. lower-log.f64 N/A

         \[\leadsto e^{\color{blue}{\log \left(\frac{\mathsf{fma}\left(10, k, 1\right)}{a \cdot {k}^{m}}\right)} \cdot -1} \]

      8. lower-/.f64 34.8

         \[\leadsto e^{\log \color{blue}{\left(\frac{\mathsf{fma}\left(10, k, 1\right)}{a \cdot {k}^{m}}\right)} \cdot -1} \]

      9. lift-*.f64 N/A

         \[\leadsto e^{\log \left(\frac{\mathsf{Rewrite=>}\left(lower-fma.f64, \left(\mathsf{fma}\left(k, 10, 1\right)\right)\right)}{\color{blue}{a \cdot {k}^{m}}}\right) \cdot -1} \]

      10. *-commutative N/A

          \[\leadsto e^{\log \left(\frac{\mathsf{Rewrite=>}\left(lower-fma.f64, \left(\mathsf{fma}\left(k, 10, 1\right)\right)\right)}{\color{blue}{{k}^{m} \cdot a}}\right) \cdot -1} \]

      11. lower-*.f64 N/A

          \[\leadsto e^{\log \left(\frac{\mathsf{Rewrite=>}\left(lower-fma.f64, \left(\mathsf{fma}\left(k, 10, 1\right)\right)\right)}{\color{blue}{{k}^{m} \cdot a}}\right) \cdot -1} \]

   7. Applied rewrites 34.8%

      \[\leadsto \color{blue}{e^{\log \left(\frac{\mathsf{fma}\left(k, 10, 1\right)}{{k}^{m} \cdot a}\right) \cdot -1}} \]

   8. Taylor expanded in k around inf

      \[\leadsto e^{\color{blue}{-1 \cdot \left(\log \left(\frac{1}{a \cdot e^{-1 \cdot \left(m \cdot \log \left(\frac{1}{k}\right)\right)}}\right) + -2 \cdot \log \left(\frac{1}{k}\right)\right)}} \]
   9. Step-by-step derivation

      1. mul-1-neg N/A

         \[\leadsto e^{\color{blue}{\mathsf{neg}\left(\left(\log \left(\frac{1}{a \cdot e^{-1 \cdot \left(m \cdot \log \left(\frac{1}{k}\right)\right)}}\right) + -2 \cdot \log \left(\frac{1}{k}\right)\right)\right)}} \]

      2. lower-neg.f64 N/A

         \[\leadsto e^{\color{blue}{\mathsf{neg}\left(\left(\log \left(\frac{1}{a \cdot e^{-1 \cdot \left(m \cdot \log \left(\frac{1}{k}\right)\right)}}\right) + -2 \cdot \log \left(\frac{1}{k}\right)\right)\right)}} \]

      3. +-commutative N/A

         \[\leadsto e^{\mathsf{neg}\left(\color{blue}{\left(-2 \cdot \log \left(\frac{1}{k}\right) + \log \left(\frac{1}{a \cdot e^{-1 \cdot \left(m \cdot \log \left(\frac{1}{k}\right)\right)}}\right)\right)}\right)} \]

      4. log-rec N/A

         \[\leadsto e^{\mathsf{neg}\left(\left(-2 \cdot \log \left(\frac{1}{k}\right) + \color{blue}{\left(\mathsf{neg}\left(\log \left(a \cdot e^{-1 \cdot \left(m \cdot \log \left(\frac{1}{k}\right)\right)}\right)\right)\right)}\right)\right)} \]

      5. unsub-neg N/A

         \[\leadsto e^{\mathsf{neg}\left(\color{blue}{\left(-2 \cdot \log \left(\frac{1}{k}\right) - \log \left(a \cdot e^{-1 \cdot \left(m \cdot \log \left(\frac{1}{k}\right)\right)}\right)\right)}\right)} \]

      6. lower--.f64 N/A

         \[\leadsto e^{\mathsf{neg}\left(\color{blue}{\left(-2 \cdot \log \left(\frac{1}{k}\right) - \log \left(a \cdot e^{-1 \cdot \left(m \cdot \log \left(\frac{1}{k}\right)\right)}\right)\right)}\right)} \]

      7. log-rec N/A

         \[\leadsto e^{\mathsf{neg}\left(\left(-2 \cdot \color{blue}{\left(\mathsf{neg}\left(\log k\right)\right)} - \log \left(a \cdot e^{-1 \cdot \left(m \cdot \log \left(\frac{1}{k}\right)\right)}\right)\right)\right)} \]

      8. mul-1-neg N/A

         \[\leadsto e^{\mathsf{neg}\left(\left(-2 \cdot \color{blue}{\left(-1 \cdot \log k\right)} - \log \left(a \cdot e^{-1 \cdot \left(m \cdot \log \left(\frac{1}{k}\right)\right)}\right)\right)\right)} \]

      9. associate-*r* N/A

         \[\leadsto e^{\mathsf{neg}\left(\left(\color{blue}{\left(-2 \cdot -1\right) \cdot \log k} - \log \left(a \cdot e^{-1 \cdot \left(m \cdot \log \left(\frac{1}{k}\right)\right)}\right)\right)\right)} \]

      10. metadata-eval N/A

          \[\leadsto e^{\mathsf{neg}\left(\left(\color{blue}{2} \cdot \log k - \log \left(a \cdot e^{-1 \cdot \left(m \cdot \log \left(\frac{1}{k}\right)\right)}\right)\right)\right)} \]

      11. lower-*.f64 N/A

          \[\leadsto e^{\mathsf{neg}\left(\left(\color{blue}{2 \cdot \log k} - \log \left(a \cdot e^{-1 \cdot \left(m \cdot \log \left(\frac{1}{k}\right)\right)}\right)\right)\right)} \]

      12. lower-log.f64 N/A

          \[\leadsto e^{\mathsf{neg}\left(\left(2 \cdot \color{blue}{\log k} - \log \left(a \cdot e^{-1 \cdot \left(m \cdot \log \left(\frac{1}{k}\right)\right)}\right)\right)\right)} \]

      13. lower-log.f64 N/A

          \[\leadsto e^{\mathsf{neg}\left(\left(2 \cdot \log k - \color{blue}{\log \left(a \cdot e^{-1 \cdot \left(m \cdot \log \left(\frac{1}{k}\right)\right)}\right)}\right)\right)} \]

      14. lower-*.f64 N/A

          \[\leadsto e^{\mathsf{neg}\left(\left(2 \cdot \log k - \log \color{blue}{\left(a \cdot e^{-1 \cdot \left(m \cdot \log \left(\frac{1}{k}\right)\right)}\right)}\right)\right)} \]

      15. mul-1-neg N/A

          \[\leadsto e^{\mathsf{neg}\left(\left(2 \cdot \log k - \log \left(a \cdot e^{\color{blue}{\mathsf{neg}\left(m \cdot \log \left(\frac{1}{k}\right)\right)}}\right)\right)\right)} \]

   10. Applied rewrites 53.4%

       \[\leadsto e^{\color{blue}{-\left(2 \cdot \log k - \log \left(a \cdot \frac{1}{{k}^{\left(-m\right)}}\right)\right)}} \]
3. Recombined 2 regimes into one program.

4. Final simplification 81.0%

   \[\leadsto \begin{array}{l} \mathbf{if}\;k \leq 0.001:\\ \;\;\;\;a \cdot {k}^{m}\\ \mathbf{else}:\\ \;\;\;\;e^{\log \left(a \cdot \frac{1}{{k}^{\left(-m\right)}}\right) - 2 \cdot \log k}\\ \end{array} \]
5. Add Preprocessing

Alternative 2: 73.0% accurate, 0.3× speedup

Math

\[\begin{array}{l} a\_m = \left|a\right| \\ a\_s = \mathsf{copysign}\left(1, a\right) \\ \begin{array}{l} t_0 := \frac{a\_m \cdot {k}^{m}}{\left(1 + k \cdot 10\right) + k \cdot k}\\ a\_s \cdot \begin{array}{l} \mathbf{if}\;t\_0 \leq 0:\\ \;\;\;\;\frac{a}{\mathsf{fma}\left(k, -k, 100\right) \cdot \left(-1 - \frac{\left(10 + \frac{100}{k}\right) + \frac{1000}{k \cdot k}}{k}\right)}\\ \mathbf{elif}\;t\_0 \leq 5 \cdot 10^{+304}:\\ \;\;\;\;\frac{a}{\mathsf{fma}\left(k, k + 10, 1\right)}\\ \mathbf{elif}\;t\_0 \leq \infty:\\ \;\;\;\;\frac{a + \frac{\mathsf{fma}\left(a, -10, \frac{a \cdot 99}{k}\right)}{k}}{k \cdot k}\\ \mathbf{else}:\\ \;\;\;\;\mathsf{fma}\left(k, \mathsf{fma}\left(a, -10, k \cdot \left(a \cdot 99\right)\right), a\right)\\ \end{array} \end{array} \end{array} \]

FPCore

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)) (+ (+ 1.0 (* k 10.0)) (* k k)))))
   (*
    a_s
    (if (<= t_0 0.0)
      (/
       a
       (*
        (fma k (- k) 100.0)
        (- -1.0 (/ (+ (+ 10.0 (/ 100.0 k)) (/ 1000.0 (* k k))) k))))
      (if (<= t_0 5e+304)
        (/ a (fma k (+ k 10.0) 1.0))
        (if (<= t_0 INFINITY)
          (/ (+ a (/ (fma a -10.0 (/ (* a 99.0) k)) k)) (* k k))
          (fma k (fma a -10.0 (* k (* a 99.0))) a)))))))

C

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)) / ((1.0 + (k * 10.0)) + (k * k));
	double tmp;
	if (t_0 <= 0.0) {
		tmp = a / (fma(k, -k, 100.0) * (-1.0 - (((10.0 + (100.0 / k)) + (1000.0 / (k * k))) / k)));
	} else if (t_0 <= 5e+304) {
		tmp = a / fma(k, (k + 10.0), 1.0);
	} else if (t_0 <= ((double) INFINITY)) {
		tmp = (a + (fma(a, -10.0, ((a * 99.0) / k)) / k)) / (k * k);
	} else {
		tmp = fma(k, fma(a, -10.0, (k * (a * 99.0))), a);
	}
	return a_s * tmp;
}

Julia

a\_m = abs(a)
a\_s = copysign(1.0, a)
function code(a_s, a_m, k, m)
	t_0 = Float64(Float64(a_m * (k ^ m)) / Float64(Float64(1.0 + Float64(k * 10.0)) + Float64(k * k)))
	tmp = 0.0
	if (t_0 <= 0.0)
		tmp = Float64(a / Float64(fma(k, Float64(-k), 100.0) * Float64(-1.0 - Float64(Float64(Float64(10.0 + Float64(100.0 / k)) + Float64(1000.0 / Float64(k * k))) / k))));
	elseif (t_0 <= 5e+304)
		tmp = Float64(a / fma(k, Float64(k + 10.0), 1.0));
	elseif (t_0 <= Inf)
		tmp = Float64(Float64(a + Float64(fma(a, -10.0, Float64(Float64(a * 99.0) / k)) / k)) / Float64(k * k));
	else
		tmp = fma(k, fma(a, -10.0, Float64(k * Float64(a * 99.0))), a);
	end
	return Float64(a_s * tmp)
end

Wolfram

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[(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]}, N[(a$95$s * If[LessEqual[t$95$0, 0.0], N[(a / N[(N[(k * (-k) + 100.0), $MachinePrecision] * N[(-1.0 - N[(N[(N[(10.0 + N[(100.0 / k), $MachinePrecision]), $MachinePrecision] + N[(1000.0 / N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, 5e+304], N[(a / N[(k * N[(k + 10.0), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, Infinity], N[(N[(a + N[(N[(a * -10.0 + N[(N[(a * 99.0), $MachinePrecision] / k), $MachinePrecision]), $MachinePrecision] / k), $MachinePrecision]), $MachinePrecision] / N[(k * k), $MachinePrecision]), $MachinePrecision], N[(k * N[(a * -10.0 + N[(k * N[(a * 99.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + a), $MachinePrecision]]]]), $MachinePrecision]]

Derivation

1. Split input into 4 regimes

2. if (/.f64 (*.f64 a (pow.f64 k m)) (+.f64 (+.f64 #s(literal 1 binary64) (*.f64 #s(literal 10 binary64) k)) (*.f64 k k))) < 0.0

   1. Initial program 96.8%

      \[\frac{a \cdot {k}^{m}}{\left(1 + 10 \cdot k\right) + k \cdot k} \]
   2. Add Preprocessing

   3. Taylor expanded in m around 0

      \[\leadsto \color{blue}{\frac{a}{1 + \left(10 \cdot k + {k}^{2}\right)}} \]
   4. Step-by-step derivation

      1. lower-/.f64 N/A

         \[\leadsto \color{blue}{\frac{a}{1 + \left(10 \cdot k + {k}^{2}\right)}} \]

      2. unpow2 N/A

         \[\leadsto \frac{a}{1 + \left(10 \cdot k + \color{blue}{k \cdot k}\right)} \]

      3. distribute-rgt-in N/A

         \[\leadsto \frac{a}{1 + \color{blue}{k \cdot \left(10 + k\right)}} \]

      4. +-commutative N/A

         \[\leadsto \frac{a}{\color{blue}{k \cdot \left(10 + k\right) + 1}} \]

      5. metadata-eval N/A

         \[\leadsto \frac{a}{k \cdot \left(\color{blue}{10 \cdot 1} + k\right) + 1} \]

      6. lft-mult-inverse N/A

         \[\leadsto \frac{a}{k \cdot \left(10 \cdot \color{blue}{\left(\frac{1}{k} \cdot k\right)} + k\right) + 1} \]

      7. associate-*l* N/A

         \[\leadsto \frac{a}{k \cdot \left(\color{blue}{\left(10 \cdot \frac{1}{k}\right) \cdot k} + k\right) + 1} \]

      8. *-lft-identity N/A

         \[\leadsto \frac{a}{k \cdot \left(\left(10 \cdot \frac{1}{k}\right) \cdot k + \color{blue}{1 \cdot k}\right) + 1} \]

      9. distribute-rgt-in N/A

         \[\leadsto \frac{a}{k \cdot \color{blue}{\left(k \cdot \left(10 \cdot \frac{1}{k} + 1\right)\right)} + 1} \]

      10. +-commutative N/A

          \[\leadsto \frac{a}{k \cdot \left(k \cdot \color{blue}{\left(1 + 10 \cdot \frac{1}{k}\right)}\right) + 1} \]

      11. *-commutative N/A

          \[\leadsto \frac{a}{k \cdot \color{blue}{\left(\left(1 + 10 \cdot \frac{1}{k}\right) \cdot k\right)} + 1} \]

      12. lower-fma.f64 N/A

          \[\leadsto \frac{a}{\color{blue}{\mathsf{fma}\left(k, \left(1 + 10 \cdot \frac{1}{k}\right) \cdot k, 1\right)}} \]

      13. *-commutative N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{k \cdot \left(1 + 10 \cdot \frac{1}{k}\right)}, 1\right)} \]

      14. +-commutative N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, k \cdot \color{blue}{\left(10 \cdot \frac{1}{k} + 1\right)}, 1\right)} \]

      15. distribute-rgt-in N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{\left(10 \cdot \frac{1}{k}\right) \cdot k + 1 \cdot k}, 1\right)} \]

      16. associate-*l* N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{10 \cdot \left(\frac{1}{k} \cdot k\right)} + 1 \cdot k, 1\right)} \]

      17. lft-mult-inverse N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, 10 \cdot \color{blue}{1} + 1 \cdot k, 1\right)} \]

      18. metadata-eval N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{10} + 1 \cdot k, 1\right)} \]

      19. *-lft-identity N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, 10 + \color{blue}{k}, 1\right)} \]

      20. lower-+.f64 51.3

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{10 + k}, 1\right)} \]

   5. Applied rewrites 51.3%

      \[\leadsto \color{blue}{\frac{a}{\mathsf{fma}\left(k, 10 + k, 1\right)}} \]

   6. Taylor expanded in k around inf

      \[\leadsto \frac{a}{{k}^{2} \cdot \color{blue}{\left(1 + 10 \cdot \frac{1}{k}\right)}} \]
   7. Step-by-step derivation

   8. Applied rewrites 43.3%

      \[\leadsto \frac{a}{k \cdot \color{blue}{\left(k + 10\right)}} \]

   10. Applied rewrites 15.0%

       \[\leadsto \frac{a}{\mathsf{fma}\left(k, -k, 100\right) \cdot \left(\mathsf{fma}\left(k, k + 10, 100\right) \cdot \color{blue}{\frac{k}{\mathsf{fma}\left(k, k \cdot \left(-k\right), 1000\right)}}\right)} \]

   11. Taylor expanded in k around -inf

       \[\leadsto \frac{a}{\mathsf{fma}\left(k, \mathsf{neg}\left(k\right), 100\right) \cdot \left(-1 \cdot \frac{10 + \left(100 \cdot \frac{1}{k} + \frac{1000}{{k}^{2}}\right)}{k} - 1\right)} \]
   12. Step-by-step derivation

   13. Applied rewrites 56.8%

       \[\leadsto \frac{a}{\mathsf{fma}\left(k, -k, 100\right) \cdot \left(-1 + \frac{\left(10 + \frac{100}{k}\right) + \frac{1000}{k \cdot k}}{\color{blue}{-k}}\right)} \]

   if 0.0 < (/.f64 (*.f64 a (pow.f64 k m)) (+.f64 (+.f64 #s(literal 1 binary64) (*.f64 #s(literal 10 binary64) k)) (*.f64 k k))) < 4.9999999999999997e304

   1. Initial program 100.0%

      \[\frac{a \cdot {k}^{m}}{\left(1 + 10 \cdot k\right) + k \cdot k} \]
   2. Add Preprocessing

   3. Taylor expanded in m around 0

      \[\leadsto \color{blue}{\frac{a}{1 + \left(10 \cdot k + {k}^{2}\right)}} \]
   4. Step-by-step derivation

      1. lower-/.f64 N/A

         \[\leadsto \color{blue}{\frac{a}{1 + \left(10 \cdot k + {k}^{2}\right)}} \]

      2. unpow2 N/A

         \[\leadsto \frac{a}{1 + \left(10 \cdot k + \color{blue}{k \cdot k}\right)} \]

      3. distribute-rgt-in N/A

         \[\leadsto \frac{a}{1 + \color{blue}{k \cdot \left(10 + k\right)}} \]

      4. +-commutative N/A

         \[\leadsto \frac{a}{\color{blue}{k \cdot \left(10 + k\right) + 1}} \]

      5. metadata-eval N/A

         \[\leadsto \frac{a}{k \cdot \left(\color{blue}{10 \cdot 1} + k\right) + 1} \]

      6. lft-mult-inverse N/A

         \[\leadsto \frac{a}{k \cdot \left(10 \cdot \color{blue}{\left(\frac{1}{k} \cdot k\right)} + k\right) + 1} \]

      7. associate-*l* N/A

         \[\leadsto \frac{a}{k \cdot \left(\color{blue}{\left(10 \cdot \frac{1}{k}\right) \cdot k} + k\right) + 1} \]

      8. *-lft-identity N/A

         \[\leadsto \frac{a}{k \cdot \left(\left(10 \cdot \frac{1}{k}\right) \cdot k + \color{blue}{1 \cdot k}\right) + 1} \]

      9. distribute-rgt-in N/A

         \[\leadsto \frac{a}{k \cdot \color{blue}{\left(k \cdot \left(10 \cdot \frac{1}{k} + 1\right)\right)} + 1} \]

      10. +-commutative N/A

          \[\leadsto \frac{a}{k \cdot \left(k \cdot \color{blue}{\left(1 + 10 \cdot \frac{1}{k}\right)}\right) + 1} \]

      11. *-commutative N/A

          \[\leadsto \frac{a}{k \cdot \color{blue}{\left(\left(1 + 10 \cdot \frac{1}{k}\right) \cdot k\right)} + 1} \]

      12. lower-fma.f64 N/A

          \[\leadsto \frac{a}{\color{blue}{\mathsf{fma}\left(k, \left(1 + 10 \cdot \frac{1}{k}\right) \cdot k, 1\right)}} \]

      13. *-commutative N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{k \cdot \left(1 + 10 \cdot \frac{1}{k}\right)}, 1\right)} \]

      14. +-commutative N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, k \cdot \color{blue}{\left(10 \cdot \frac{1}{k} + 1\right)}, 1\right)} \]

      15. distribute-rgt-in N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{\left(10 \cdot \frac{1}{k}\right) \cdot k + 1 \cdot k}, 1\right)} \]

      16. associate-*l* N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{10 \cdot \left(\frac{1}{k} \cdot k\right)} + 1 \cdot k, 1\right)} \]

      17. lft-mult-inverse N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, 10 \cdot \color{blue}{1} + 1 \cdot k, 1\right)} \]

      18. metadata-eval N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{10} + 1 \cdot k, 1\right)} \]

      19. *-lft-identity N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, 10 + \color{blue}{k}, 1\right)} \]

      20. lower-+.f64 97.5

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{10 + k}, 1\right)} \]

   5. Applied rewrites 97.5%

      \[\leadsto \color{blue}{\frac{a}{\mathsf{fma}\left(k, 10 + k, 1\right)}} \]

   if 4.9999999999999997e304 < (/.f64 (*.f64 a (pow.f64 k m)) (+.f64 (+.f64 #s(literal 1 binary64) (*.f64 #s(literal 10 binary64) k)) (*.f64 k k))) < +inf.0

   1. Initial program 100.0%

      \[\frac{a \cdot {k}^{m}}{\left(1 + 10 \cdot k\right) + k \cdot k} \]
   2. Add Preprocessing

   3. Taylor expanded in m around 0

      \[\leadsto \color{blue}{\frac{a}{1 + \left(10 \cdot k + {k}^{2}\right)}} \]
   4. Step-by-step derivation

      1. lower-/.f64 N/A

         \[\leadsto \color{blue}{\frac{a}{1 + \left(10 \cdot k + {k}^{2}\right)}} \]

      2. unpow2 N/A

         \[\leadsto \frac{a}{1 + \left(10 \cdot k + \color{blue}{k \cdot k}\right)} \]

      3. distribute-rgt-in N/A

         \[\leadsto \frac{a}{1 + \color{blue}{k \cdot \left(10 + k\right)}} \]

      4. +-commutative N/A

         \[\leadsto \frac{a}{\color{blue}{k \cdot \left(10 + k\right) + 1}} \]

      5. metadata-eval N/A

         \[\leadsto \frac{a}{k \cdot \left(\color{blue}{10 \cdot 1} + k\right) + 1} \]

      6. lft-mult-inverse N/A

         \[\leadsto \frac{a}{k \cdot \left(10 \cdot \color{blue}{\left(\frac{1}{k} \cdot k\right)} + k\right) + 1} \]

      7. associate-*l* N/A

         \[\leadsto \frac{a}{k \cdot \left(\color{blue}{\left(10 \cdot \frac{1}{k}\right) \cdot k} + k\right) + 1} \]

      8. *-lft-identity N/A

         \[\leadsto \frac{a}{k \cdot \left(\left(10 \cdot \frac{1}{k}\right) \cdot k + \color{blue}{1 \cdot k}\right) + 1} \]

      9. distribute-rgt-in N/A

         \[\leadsto \frac{a}{k \cdot \color{blue}{\left(k \cdot \left(10 \cdot \frac{1}{k} + 1\right)\right)} + 1} \]

      10. +-commutative N/A

          \[\leadsto \frac{a}{k \cdot \left(k \cdot \color{blue}{\left(1 + 10 \cdot \frac{1}{k}\right)}\right) + 1} \]

      11. *-commutative N/A

          \[\leadsto \frac{a}{k \cdot \color{blue}{\left(\left(1 + 10 \cdot \frac{1}{k}\right) \cdot k\right)} + 1} \]

      12. lower-fma.f64 N/A

          \[\leadsto \frac{a}{\color{blue}{\mathsf{fma}\left(k, \left(1 + 10 \cdot \frac{1}{k}\right) \cdot k, 1\right)}} \]

      13. *-commutative N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{k \cdot \left(1 + 10 \cdot \frac{1}{k}\right)}, 1\right)} \]

      14. +-commutative N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, k \cdot \color{blue}{\left(10 \cdot \frac{1}{k} + 1\right)}, 1\right)} \]

      15. distribute-rgt-in N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{\left(10 \cdot \frac{1}{k}\right) \cdot k + 1 \cdot k}, 1\right)} \]

      16. associate-*l* N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{10 \cdot \left(\frac{1}{k} \cdot k\right)} + 1 \cdot k, 1\right)} \]

      17. lft-mult-inverse N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, 10 \cdot \color{blue}{1} + 1 \cdot k, 1\right)} \]

      18. metadata-eval N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{10} + 1 \cdot k, 1\right)} \]

      19. *-lft-identity N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, 10 + \color{blue}{k}, 1\right)} \]

      20. lower-+.f64 3.2

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{10 + k}, 1\right)} \]

   5. Applied rewrites 3.2%

      \[\leadsto \color{blue}{\frac{a}{\mathsf{fma}\left(k, 10 + k, 1\right)}} \]

   6. Taylor expanded in k around inf

      \[\leadsto \frac{\left(a + -1 \cdot \frac{a + -100 \cdot a}{{k}^{2}}\right) - 10 \cdot \frac{a}{k}}{\color{blue}{{k}^{2}}} \]

   8. Applied rewrites 37.5%

      \[\leadsto \frac{a + \frac{\mathsf{fma}\left(a, -10, \frac{a \cdot 99}{k}\right)}{k}}{\color{blue}{k \cdot k}} \]

   if +inf.0 < (/.f64 (*.f64 a (pow.f64 k m)) (+.f64 (+.f64 #s(literal 1 binary64) (*.f64 #s(literal 10 binary64) k)) (*.f64 k k)))

   1. Initial program 0.0%

      \[\frac{a \cdot {k}^{m}}{\left(1 + 10 \cdot k\right) + k \cdot k} \]
   2. Add Preprocessing

   3. Taylor expanded in m around 0

      \[\leadsto \color{blue}{\frac{a}{1 + \left(10 \cdot k + {k}^{2}\right)}} \]
   4. Step-by-step derivation

      1. lower-/.f64 N/A

         \[\leadsto \color{blue}{\frac{a}{1 + \left(10 \cdot k + {k}^{2}\right)}} \]

      2. unpow2 N/A

         \[\leadsto \frac{a}{1 + \left(10 \cdot k + \color{blue}{k \cdot k}\right)} \]

      3. distribute-rgt-in N/A

         \[\leadsto \frac{a}{1 + \color{blue}{k \cdot \left(10 + k\right)}} \]

      4. +-commutative N/A

         \[\leadsto \frac{a}{\color{blue}{k \cdot \left(10 + k\right) + 1}} \]

      5. metadata-eval N/A

         \[\leadsto \frac{a}{k \cdot \left(\color{blue}{10 \cdot 1} + k\right) + 1} \]

      6. lft-mult-inverse N/A

         \[\leadsto \frac{a}{k \cdot \left(10 \cdot \color{blue}{\left(\frac{1}{k} \cdot k\right)} + k\right) + 1} \]

      7. associate-*l* N/A

         \[\leadsto \frac{a}{k \cdot \left(\color{blue}{\left(10 \cdot \frac{1}{k}\right) \cdot k} + k\right) + 1} \]

      8. *-lft-identity N/A

         \[\leadsto \frac{a}{k \cdot \left(\left(10 \cdot \frac{1}{k}\right) \cdot k + \color{blue}{1 \cdot k}\right) + 1} \]

      9. distribute-rgt-in N/A

         \[\leadsto \frac{a}{k \cdot \color{blue}{\left(k \cdot \left(10 \cdot \frac{1}{k} + 1\right)\right)} + 1} \]

      10. +-commutative N/A

          \[\leadsto \frac{a}{k \cdot \left(k \cdot \color{blue}{\left(1 + 10 \cdot \frac{1}{k}\right)}\right) + 1} \]

      11. *-commutative N/A

          \[\leadsto \frac{a}{k \cdot \color{blue}{\left(\left(1 + 10 \cdot \frac{1}{k}\right) \cdot k\right)} + 1} \]

      12. lower-fma.f64 N/A

          \[\leadsto \frac{a}{\color{blue}{\mathsf{fma}\left(k, \left(1 + 10 \cdot \frac{1}{k}\right) \cdot k, 1\right)}} \]

      13. *-commutative N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{k \cdot \left(1 + 10 \cdot \frac{1}{k}\right)}, 1\right)} \]

      14. +-commutative N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, k \cdot \color{blue}{\left(10 \cdot \frac{1}{k} + 1\right)}, 1\right)} \]

      15. distribute-rgt-in N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{\left(10 \cdot \frac{1}{k}\right) \cdot k + 1 \cdot k}, 1\right)} \]

      16. associate-*l* N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{10 \cdot \left(\frac{1}{k} \cdot k\right)} + 1 \cdot k, 1\right)} \]

      17. lft-mult-inverse N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, 10 \cdot \color{blue}{1} + 1 \cdot k, 1\right)} \]

      18. metadata-eval N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{10} + 1 \cdot k, 1\right)} \]

      19. *-lft-identity N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, 10 + \color{blue}{k}, 1\right)} \]

      20. lower-+.f64 5.7

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{10 + k}, 1\right)} \]

   5. Applied rewrites 5.7%

      \[\leadsto \color{blue}{\frac{a}{\mathsf{fma}\left(k, 10 + k, 1\right)}} \]

   6. Taylor expanded in k around 0

      \[\leadsto a + \color{blue}{k \cdot \left(-1 \cdot \left(k \cdot \left(a + -100 \cdot a\right)\right) - 10 \cdot a\right)} \]
   7. Step-by-step derivation

   8. Applied rewrites 68.7%

      \[\leadsto \mathsf{fma}\left(k, \color{blue}{\mathsf{fma}\left(a, -10, k \cdot \left(a \cdot 99\right)\right)}, a\right) \]
3. Recombined 4 regimes into one program.

4. Final simplification 61.5%

   \[\leadsto \begin{array}{l} \mathbf{if}\;\frac{a \cdot {k}^{m}}{\left(1 + k \cdot 10\right) + k \cdot k} \leq 0:\\ \;\;\;\;\frac{a}{\mathsf{fma}\left(k, -k, 100\right) \cdot \left(-1 - \frac{\left(10 + \frac{100}{k}\right) + \frac{1000}{k \cdot k}}{k}\right)}\\ \mathbf{elif}\;\frac{a \cdot {k}^{m}}{\left(1 + k \cdot 10\right) + k \cdot k} \leq 5 \cdot 10^{+304}:\\ \;\;\;\;\frac{a}{\mathsf{fma}\left(k, k + 10, 1\right)}\\ \mathbf{elif}\;\frac{a \cdot {k}^{m}}{\left(1 + k \cdot 10\right) + k \cdot k} \leq \infty:\\ \;\;\;\;\frac{a + \frac{\mathsf{fma}\left(a, -10, \frac{a \cdot 99}{k}\right)}{k}}{k \cdot k}\\ \mathbf{else}:\\ \;\;\;\;\mathsf{fma}\left(k, \mathsf{fma}\left(a, -10, k \cdot \left(a \cdot 99\right)\right), a\right)\\ \end{array} \]
5. Add Preprocessing

Alternative 3: 97.1% accurate, 0.5× speedup

Math

\[\begin{array}{l} a\_m = \left|a\right| \\ a\_s = \mathsf{copysign}\left(1, a\right) \\ \begin{array}{l} t_0 := \frac{a\_m \cdot {k}^{m}}{\left(1 + k \cdot 10\right) + k \cdot k}\\ a\_s \cdot \begin{array}{l} \mathbf{if}\;t\_0 \leq 4 \cdot 10^{+155}:\\ \;\;\;\;t\_0\\ \mathbf{else}:\\ \;\;\;\;a \cdot {k}^{m}\\ \end{array} \end{array} \end{array} \]

FPCore

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)) (+ (+ 1.0 (* k 10.0)) (* k k)))))
   (* a_s (if (<= t_0 4e+155) t_0 (* a (pow k m))))))

C

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)) / ((1.0 + (k * 10.0)) + (k * k));
	double tmp;
	if (t_0 <= 4e+155) {
		tmp = t_0;
	} else {
		tmp = a * pow(k, m);
	}
	return a_s * tmp;
}

Fortran

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)) / ((1.0d0 + (k * 10.0d0)) + (k * k))
    if (t_0 <= 4d+155) then
        tmp = t_0
    else
        tmp = a * (k ** m)
    end if
    code = a_s * tmp
end function

Java

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)) / ((1.0 + (k * 10.0)) + (k * k));
	double tmp;
	if (t_0 <= 4e+155) {
		tmp = t_0;
	} else {
		tmp = a * Math.pow(k, m);
	}
	return a_s * tmp;
}

Python

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)) / ((1.0 + (k * 10.0)) + (k * k))
	tmp = 0
	if t_0 <= 4e+155:
		tmp = t_0
	else:
		tmp = a * math.pow(k, m)
	return a_s * tmp

Julia

a\_m = abs(a)
a\_s = copysign(1.0, a)
function code(a_s, a_m, k, m)
	t_0 = Float64(Float64(a_m * (k ^ m)) / Float64(Float64(1.0 + Float64(k * 10.0)) + Float64(k * k)))
	tmp = 0.0
	if (t_0 <= 4e+155)
		tmp = t_0;
	else
		tmp = Float64(a * (k ^ m));
	end
	return Float64(a_s * tmp)
end

MATLAB

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)) / ((1.0 + (k * 10.0)) + (k * k));
	tmp = 0.0;
	if (t_0 <= 4e+155)
		tmp = t_0;
	else
		tmp = a * (k ^ m);
	end
	tmp_2 = a_s * tmp;
end

Wolfram

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[(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]}, N[(a$95$s * If[LessEqual[t$95$0, 4e+155], t$95$0, N[(a * N[Power[k, m], $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]]

Derivation

1. Split input into 2 regimes

2. if (/.f64 (*.f64 a (pow.f64 k m)) (+.f64 (+.f64 #s(literal 1 binary64) (*.f64 #s(literal 10 binary64) k)) (*.f64 k k))) < 4.00000000000000003e155

   1. Initial program 97.2%

      \[\frac{a \cdot {k}^{m}}{\left(1 + 10 \cdot k\right) + k \cdot k} \]
   2. Add Preprocessing

   if 4.00000000000000003e155 < (/.f64 (*.f64 a (pow.f64 k m)) (+.f64 (+.f64 #s(literal 1 binary64) (*.f64 #s(literal 10 binary64) k)) (*.f64 k k)))

   1. Initial program 61.9%

      \[\frac{a \cdot {k}^{m}}{\left(1 + 10 \cdot k\right) + k \cdot k} \]
   2. Add Preprocessing

   3. Taylor expanded in k around 0

      \[\leadsto \color{blue}{a \cdot {k}^{m}} \]
   4. Step-by-step derivation

      1. lower-*.f64 N/A

         \[\leadsto \color{blue}{a \cdot {k}^{m}} \]

      2. lower-pow.f64 100.0

         \[\leadsto a \cdot \color{blue}{{k}^{m}} \]

   5. Applied rewrites 100.0%

      \[\leadsto \color{blue}{a \cdot {k}^{m}} \]
3. Recombined 2 regimes into one program.

4. Final simplification 97.9%

   \[\leadsto \begin{array}{l} \mathbf{if}\;\frac{a \cdot {k}^{m}}{\left(1 + k \cdot 10\right) + k \cdot k} \leq 4 \cdot 10^{+155}:\\ \;\;\;\;\frac{a \cdot {k}^{m}}{\left(1 + k \cdot 10\right) + k \cdot k}\\ \mathbf{else}:\\ \;\;\;\;a \cdot {k}^{m}\\ \end{array} \]
5. Add Preprocessing

Alternative 4: 28.6% accurate, 0.9× speedup

Math

\[\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 5 \cdot 10^{-301}:\\ \;\;\;\;\frac{a}{k \cdot 10}\\ \mathbf{else}:\\ \;\;\;\;a \cdot \mathsf{fma}\left(k, -10, 1\right)\\ \end{array} \end{array} \]

FPCore

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))) 5e-301)
    (/ a (* k 10.0))
    (* a (fma k -10.0 1.0)))))

C

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))) <= 5e-301) {
		tmp = a / (k * 10.0);
	} else {
		tmp = a * fma(k, -10.0, 1.0);
	}
	return a_s * tmp;
}

Julia

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))) <= 5e-301)
		tmp = Float64(a / Float64(k * 10.0));
	else
		tmp = Float64(a * fma(k, -10.0, 1.0));
	end
	return Float64(a_s * tmp)
end

Wolfram

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], 5e-301], N[(a / N[(k * 10.0), $MachinePrecision]), $MachinePrecision], N[(a * N[(k * -10.0 + 1.0), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]

Derivation

1. Split input into 2 regimes

2. if (/.f64 (*.f64 a (pow.f64 k m)) (+.f64 (+.f64 #s(literal 1 binary64) (*.f64 #s(literal 10 binary64) k)) (*.f64 k k))) < 5.00000000000000013e-301

   1. Initial program 96.8%

      \[\frac{a \cdot {k}^{m}}{\left(1 + 10 \cdot k\right) + k \cdot k} \]
   2. Add Preprocessing

   3. Taylor expanded in m around 0

      \[\leadsto \color{blue}{\frac{a}{1 + \left(10 \cdot k + {k}^{2}\right)}} \]
   4. Step-by-step derivation

      1. lower-/.f64 N/A

         \[\leadsto \color{blue}{\frac{a}{1 + \left(10 \cdot k + {k}^{2}\right)}} \]

      2. unpow2 N/A

         \[\leadsto \frac{a}{1 + \left(10 \cdot k + \color{blue}{k \cdot k}\right)} \]

      3. distribute-rgt-in N/A

         \[\leadsto \frac{a}{1 + \color{blue}{k \cdot \left(10 + k\right)}} \]

      4. +-commutative N/A

         \[\leadsto \frac{a}{\color{blue}{k \cdot \left(10 + k\right) + 1}} \]

      5. metadata-eval N/A

         \[\leadsto \frac{a}{k \cdot \left(\color{blue}{10 \cdot 1} + k\right) + 1} \]

      6. lft-mult-inverse N/A

         \[\leadsto \frac{a}{k \cdot \left(10 \cdot \color{blue}{\left(\frac{1}{k} \cdot k\right)} + k\right) + 1} \]

      7. associate-*l* N/A

         \[\leadsto \frac{a}{k \cdot \left(\color{blue}{\left(10 \cdot \frac{1}{k}\right) \cdot k} + k\right) + 1} \]

      8. *-lft-identity N/A

         \[\leadsto \frac{a}{k \cdot \left(\left(10 \cdot \frac{1}{k}\right) \cdot k + \color{blue}{1 \cdot k}\right) + 1} \]

      9. distribute-rgt-in N/A

         \[\leadsto \frac{a}{k \cdot \color{blue}{\left(k \cdot \left(10 \cdot \frac{1}{k} + 1\right)\right)} + 1} \]

      10. +-commutative N/A

          \[\leadsto \frac{a}{k \cdot \left(k \cdot \color{blue}{\left(1 + 10 \cdot \frac{1}{k}\right)}\right) + 1} \]

      11. *-commutative N/A

          \[\leadsto \frac{a}{k \cdot \color{blue}{\left(\left(1 + 10 \cdot \frac{1}{k}\right) \cdot k\right)} + 1} \]

      12. lower-fma.f64 N/A

          \[\leadsto \frac{a}{\color{blue}{\mathsf{fma}\left(k, \left(1 + 10 \cdot \frac{1}{k}\right) \cdot k, 1\right)}} \]

      13. *-commutative N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{k \cdot \left(1 + 10 \cdot \frac{1}{k}\right)}, 1\right)} \]

      14. +-commutative N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, k \cdot \color{blue}{\left(10 \cdot \frac{1}{k} + 1\right)}, 1\right)} \]

      15. distribute-rgt-in N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{\left(10 \cdot \frac{1}{k}\right) \cdot k + 1 \cdot k}, 1\right)} \]

      16. associate-*l* N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{10 \cdot \left(\frac{1}{k} \cdot k\right)} + 1 \cdot k, 1\right)} \]

      17. lft-mult-inverse N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, 10 \cdot \color{blue}{1} + 1 \cdot k, 1\right)} \]

      18. metadata-eval N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{10} + 1 \cdot k, 1\right)} \]

      19. *-lft-identity N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, 10 + \color{blue}{k}, 1\right)} \]

      20. lower-+.f64 51.9

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{10 + k}, 1\right)} \]

   5. Applied rewrites 51.9%

      \[\leadsto \color{blue}{\frac{a}{\mathsf{fma}\left(k, 10 + k, 1\right)}} \]

   6. Taylor expanded in k around inf

      \[\leadsto \frac{a}{{k}^{2} \cdot \color{blue}{\left(1 + 10 \cdot \frac{1}{k}\right)}} \]
   7. Step-by-step derivation

   8. Applied rewrites 44.0%

      \[\leadsto \frac{a}{k \cdot \color{blue}{\left(k + 10\right)}} \]

   9. Taylor expanded in k around 0

      \[\leadsto \frac{a}{10 \cdot k} \]
   10. Step-by-step derivation

   11. Applied rewrites 19.6%

       \[\leadsto \frac{a}{k \cdot 10} \]

   if 5.00000000000000013e-301 < (/.f64 (*.f64 a (pow.f64 k m)) (+.f64 (+.f64 #s(literal 1 binary64) (*.f64 #s(literal 10 binary64) k)) (*.f64 k k)))

   1. Initial program 72.1%

      \[\frac{a \cdot {k}^{m}}{\left(1 + 10 \cdot k\right) + k \cdot k} \]
   2. Add Preprocessing

   3. Taylor expanded in m around 0

      \[\leadsto \color{blue}{\frac{a}{1 + \left(10 \cdot k + {k}^{2}\right)}} \]
   4. Step-by-step derivation

      1. lower-/.f64 N/A

         \[\leadsto \color{blue}{\frac{a}{1 + \left(10 \cdot k + {k}^{2}\right)}} \]

      2. unpow2 N/A

         \[\leadsto \frac{a}{1 + \left(10 \cdot k + \color{blue}{k \cdot k}\right)} \]

      3. distribute-rgt-in N/A

         \[\leadsto \frac{a}{1 + \color{blue}{k \cdot \left(10 + k\right)}} \]

      4. +-commutative N/A

         \[\leadsto \frac{a}{\color{blue}{k \cdot \left(10 + k\right) + 1}} \]

      5. metadata-eval N/A

         \[\leadsto \frac{a}{k \cdot \left(\color{blue}{10 \cdot 1} + k\right) + 1} \]

      6. lft-mult-inverse N/A

         \[\leadsto \frac{a}{k \cdot \left(10 \cdot \color{blue}{\left(\frac{1}{k} \cdot k\right)} + k\right) + 1} \]

      7. associate-*l* N/A

         \[\leadsto \frac{a}{k \cdot \left(\color{blue}{\left(10 \cdot \frac{1}{k}\right) \cdot k} + k\right) + 1} \]

      8. *-lft-identity N/A

         \[\leadsto \frac{a}{k \cdot \left(\left(10 \cdot \frac{1}{k}\right) \cdot k + \color{blue}{1 \cdot k}\right) + 1} \]

      9. distribute-rgt-in N/A

         \[\leadsto \frac{a}{k \cdot \color{blue}{\left(k \cdot \left(10 \cdot \frac{1}{k} + 1\right)\right)} + 1} \]

      10. +-commutative N/A

          \[\leadsto \frac{a}{k \cdot \left(k \cdot \color{blue}{\left(1 + 10 \cdot \frac{1}{k}\right)}\right) + 1} \]

      11. *-commutative N/A

          \[\leadsto \frac{a}{k \cdot \color{blue}{\left(\left(1 + 10 \cdot \frac{1}{k}\right) \cdot k\right)} + 1} \]

      12. lower-fma.f64 N/A

          \[\leadsto \frac{a}{\color{blue}{\mathsf{fma}\left(k, \left(1 + 10 \cdot \frac{1}{k}\right) \cdot k, 1\right)}} \]

      13. *-commutative N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{k \cdot \left(1 + 10 \cdot \frac{1}{k}\right)}, 1\right)} \]

      14. +-commutative N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, k \cdot \color{blue}{\left(10 \cdot \frac{1}{k} + 1\right)}, 1\right)} \]

      15. distribute-rgt-in N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{\left(10 \cdot \frac{1}{k}\right) \cdot k + 1 \cdot k}, 1\right)} \]

      16. associate-*l* N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{10 \cdot \left(\frac{1}{k} \cdot k\right)} + 1 \cdot k, 1\right)} \]

      17. lft-mult-inverse N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, 10 \cdot \color{blue}{1} + 1 \cdot k, 1\right)} \]

      18. metadata-eval N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{10} + 1 \cdot k, 1\right)} \]

      19. *-lft-identity N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, 10 + \color{blue}{k}, 1\right)} \]

      20. lower-+.f64 41.1

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{10 + k}, 1\right)} \]

   5. Applied rewrites 41.1%

      \[\leadsto \color{blue}{\frac{a}{\mathsf{fma}\left(k, 10 + k, 1\right)}} \]

   6. Taylor expanded in k around 0

      \[\leadsto a + \color{blue}{-10 \cdot \left(a \cdot k\right)} \]
   7. Step-by-step derivation

   8. Applied rewrites 35.3%

      \[\leadsto \mathsf{fma}\left(a, \color{blue}{k \cdot -10}, a\right) \]
   9. Step-by-step derivation

   10. Applied rewrites 35.3%

       \[\leadsto \mathsf{fma}\left(k, -10, 1\right) \cdot a \]
3. Recombined 2 regimes into one program.

4. Final simplification 24.8%

   \[\leadsto \begin{array}{l} \mathbf{if}\;\frac{a \cdot {k}^{m}}{\left(1 + k \cdot 10\right) + k \cdot k} \leq 5 \cdot 10^{-301}:\\ \;\;\;\;\frac{a}{k \cdot 10}\\ \mathbf{else}:\\ \;\;\;\;a \cdot \mathsf{fma}\left(k, -10, 1\right)\\ \end{array} \]
5. Add Preprocessing

Alternative 5: 97.5% accurate, 1.0× speedup

Math

\[\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 5 \cdot 10^{-6}:\\ \;\;\;\;a\_m \cdot \frac{{k}^{m}}{\mathsf{fma}\left(k, k + 10, 1\right)}\\ \mathbf{else}:\\ \;\;\;\;a \cdot {k}^{m}\\ \end{array} \end{array} \]

FPCore

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 5e-6)
    (* a_m (/ (pow k m) (fma k (+ k 10.0) 1.0)))
    (* a (pow k m)))))

C

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 <= 5e-6) {
		tmp = a_m * (pow(k, m) / fma(k, (k + 10.0), 1.0));
	} else {
		tmp = a * pow(k, m);
	}
	return a_s * tmp;
}

Julia

a\_m = abs(a)
a\_s = copysign(1.0, a)
function code(a_s, a_m, k, m)
	tmp = 0.0
	if (m <= 5e-6)
		tmp = Float64(a_m * Float64((k ^ m) / fma(k, Float64(k + 10.0), 1.0)));
	else
		tmp = Float64(a * (k ^ m));
	end
	return Float64(a_s * tmp)
end

Wolfram

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, 5e-6], N[(a$95$m * N[(N[Power[k, m], $MachinePrecision] / N[(k * N[(k + 10.0), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(a * N[Power[k, m], $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]

Derivation

1. Split input into 2 regimes

2. if m < 5.00000000000000041e-6

   1. Initial program 96.3%

      \[\frac{a \cdot {k}^{m}}{\left(1 + 10 \cdot k\right) + k \cdot k} \]
   2. Add Preprocessing
   3. Step-by-step derivation

      1. lift-/.f64 N/A

         \[\leadsto \color{blue}{\frac{a \cdot {k}^{m}}{\left(1 + 10 \cdot k\right) + k \cdot k}} \]

      2. lift-*.f64 N/A

         \[\leadsto \frac{\color{blue}{a \cdot {k}^{m}}}{\left(1 + 10 \cdot k\right) + k \cdot k} \]

      3. associate-/l* N/A

         \[\leadsto \color{blue}{a \cdot \frac{{k}^{m}}{\left(1 + 10 \cdot k\right) + k \cdot k}} \]

      4. *-commutative N/A

         \[\leadsto \color{blue}{\frac{{k}^{m}}{\left(1 + 10 \cdot k\right) + k \cdot k} \cdot a} \]

      5. lower-*.f64 N/A

         \[\leadsto \color{blue}{\frac{{k}^{m}}{\left(1 + 10 \cdot k\right) + k \cdot k} \cdot a} \]

      6. lower-/.f64 96.3

         \[\leadsto \color{blue}{\frac{{k}^{m}}{\left(1 + 10 \cdot k\right) + k \cdot k}} \cdot a \]

      7. lift-+.f64 N/A

         \[\leadsto \frac{{k}^{m}}{\color{blue}{\left(1 + 10 \cdot k\right) + k \cdot k}} \cdot a \]

      8. lift-+.f64 N/A

         \[\leadsto \frac{{k}^{m}}{\color{blue}{\left(1 + 10 \cdot k\right)} + k \cdot k} \cdot a \]

      9. associate-+l+ N/A

         \[\leadsto \frac{{k}^{m}}{\color{blue}{1 + \left(10 \cdot k + k \cdot k\right)}} \cdot a \]

      10. +-commutative N/A

          \[\leadsto \frac{{k}^{m}}{\color{blue}{\left(10 \cdot k + k \cdot k\right) + 1}} \cdot a \]

      11. lift-*.f64 N/A

          \[\leadsto \frac{{k}^{m}}{\left(\color{blue}{10 \cdot k} + k \cdot k\right) + 1} \cdot a \]

      12. lift-*.f64 N/A

          \[\leadsto \frac{{k}^{m}}{\left(10 \cdot k + \color{blue}{k \cdot k}\right) + 1} \cdot a \]

      13. distribute-rgt-out N/A

          \[\leadsto \frac{{k}^{m}}{\color{blue}{k \cdot \left(10 + k\right)} + 1} \cdot a \]

      14. lower-fma.f64 N/A

          \[\leadsto \frac{{k}^{m}}{\color{blue}{\mathsf{fma}\left(k, 10 + k, 1\right)}} \cdot a \]

      15. lower-+.f64 96.8

          \[\leadsto \frac{{k}^{m}}{\mathsf{fma}\left(k, \color{blue}{10 + k}, 1\right)} \cdot a \]

   4. Applied rewrites 96.8%

      \[\leadsto \color{blue}{\frac{{k}^{m}}{\mathsf{fma}\left(k, 10 + k, 1\right)} \cdot a} \]

   if 5.00000000000000041e-6 < m

   1. Initial program 72.0%

      \[\frac{a \cdot {k}^{m}}{\left(1 + 10 \cdot k\right) + k \cdot k} \]
   2. Add Preprocessing

   3. Taylor expanded in k around 0

      \[\leadsto \color{blue}{a \cdot {k}^{m}} \]
   4. Step-by-step derivation

      1. lower-*.f64 N/A

         \[\leadsto \color{blue}{a \cdot {k}^{m}} \]

      2. lower-pow.f64 100.0

         \[\leadsto a \cdot \color{blue}{{k}^{m}} \]

   5. Applied rewrites 100.0%

      \[\leadsto \color{blue}{a \cdot {k}^{m}} \]
3. Recombined 2 regimes into one program.

4. Final simplification 97.9%

   \[\leadsto \begin{array}{l} \mathbf{if}\;m \leq 5 \cdot 10^{-6}:\\ \;\;\;\;a \cdot \frac{{k}^{m}}{\mathsf{fma}\left(k, k + 10, 1\right)}\\ \mathbf{else}:\\ \;\;\;\;a \cdot {k}^{m}\\ \end{array} \]
5. Add Preprocessing

Alternative 6: 97.0% accurate, 1.0× speedup

Math

\[\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 -2 \cdot 10^{-17}:\\ \;\;\;\;\frac{a\_m \cdot {k}^{m}}{k \cdot k}\\ \mathbf{elif}\;m \leq 5 \cdot 10^{-6}:\\ \;\;\;\;\frac{a}{\mathsf{fma}\left(k, k + 10, 1\right)}\\ \mathbf{else}:\\ \;\;\;\;a \cdot {k}^{m}\\ \end{array} \end{array} \]

FPCore

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 -2e-17)
    (/ (* a_m (pow k m)) (* k k))
    (if (<= m 5e-6) (/ a (fma k (+ k 10.0) 1.0)) (* a (pow k m))))))

C

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 <= -2e-17) {
		tmp = (a_m * pow(k, m)) / (k * k);
	} else if (m <= 5e-6) {
		tmp = a / fma(k, (k + 10.0), 1.0);
	} else {
		tmp = a * pow(k, m);
	}
	return a_s * tmp;
}

Julia

a\_m = abs(a)
a\_s = copysign(1.0, a)
function code(a_s, a_m, k, m)
	tmp = 0.0
	if (m <= -2e-17)
		tmp = Float64(Float64(a_m * (k ^ m)) / Float64(k * k));
	elseif (m <= 5e-6)
		tmp = Float64(a / fma(k, Float64(k + 10.0), 1.0));
	else
		tmp = Float64(a * (k ^ m));
	end
	return Float64(a_s * tmp)
end

Wolfram

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, -2e-17], N[(N[(a$95$m * N[Power[k, m], $MachinePrecision]), $MachinePrecision] / N[(k * k), $MachinePrecision]), $MachinePrecision], If[LessEqual[m, 5e-6], N[(a / N[(k * N[(k + 10.0), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision], N[(a * N[Power[k, m], $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]

Derivation

1. Split input into 3 regimes

2. if m < -2.00000000000000014e-17

   1. Initial program 98.7%

      \[\frac{a \cdot {k}^{m}}{\left(1 + 10 \cdot k\right) + k \cdot k} \]
   2. Add Preprocessing

   3. Taylor expanded in k around inf

      \[\leadsto \frac{a \cdot {k}^{m}}{\color{blue}{{k}^{2}}} \]
   4. Step-by-step derivation

      1. unpow2 N/A

         \[\leadsto \frac{a \cdot {k}^{m}}{\color{blue}{k \cdot k}} \]

      2. lower-*.f64 100.0

         \[\leadsto \frac{a \cdot {k}^{m}}{\color{blue}{k \cdot k}} \]

   5. Applied rewrites 100.0%

      \[\leadsto \frac{a \cdot {k}^{m}}{\color{blue}{k \cdot k}} \]

   if -2.00000000000000014e-17 < m < 5.00000000000000041e-6

   1. Initial program 94.3%

      \[\frac{a \cdot {k}^{m}}{\left(1 + 10 \cdot k\right) + k \cdot k} \]
   2. Add Preprocessing

   3. Taylor expanded in m around 0

      \[\leadsto \color{blue}{\frac{a}{1 + \left(10 \cdot k + {k}^{2}\right)}} \]
   4. Step-by-step derivation

      1. lower-/.f64 N/A

         \[\leadsto \color{blue}{\frac{a}{1 + \left(10 \cdot k + {k}^{2}\right)}} \]

      2. unpow2 N/A

         \[\leadsto \frac{a}{1 + \left(10 \cdot k + \color{blue}{k \cdot k}\right)} \]

      3. distribute-rgt-in N/A

         \[\leadsto \frac{a}{1 + \color{blue}{k \cdot \left(10 + k\right)}} \]

      4. +-commutative N/A

         \[\leadsto \frac{a}{\color{blue}{k \cdot \left(10 + k\right) + 1}} \]

      5. metadata-eval N/A

         \[\leadsto \frac{a}{k \cdot \left(\color{blue}{10 \cdot 1} + k\right) + 1} \]

      6. lft-mult-inverse N/A

         \[\leadsto \frac{a}{k \cdot \left(10 \cdot \color{blue}{\left(\frac{1}{k} \cdot k\right)} + k\right) + 1} \]

      7. associate-*l* N/A

         \[\leadsto \frac{a}{k \cdot \left(\color{blue}{\left(10 \cdot \frac{1}{k}\right) \cdot k} + k\right) + 1} \]

      8. *-lft-identity N/A

         \[\leadsto \frac{a}{k \cdot \left(\left(10 \cdot \frac{1}{k}\right) \cdot k + \color{blue}{1 \cdot k}\right) + 1} \]

      9. distribute-rgt-in N/A

         \[\leadsto \frac{a}{k \cdot \color{blue}{\left(k \cdot \left(10 \cdot \frac{1}{k} + 1\right)\right)} + 1} \]

      10. +-commutative N/A

          \[\leadsto \frac{a}{k \cdot \left(k \cdot \color{blue}{\left(1 + 10 \cdot \frac{1}{k}\right)}\right) + 1} \]

      11. *-commutative N/A

          \[\leadsto \frac{a}{k \cdot \color{blue}{\left(\left(1 + 10 \cdot \frac{1}{k}\right) \cdot k\right)} + 1} \]

      12. lower-fma.f64 N/A

          \[\leadsto \frac{a}{\color{blue}{\mathsf{fma}\left(k, \left(1 + 10 \cdot \frac{1}{k}\right) \cdot k, 1\right)}} \]

      13. *-commutative N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{k \cdot \left(1 + 10 \cdot \frac{1}{k}\right)}, 1\right)} \]

      14. +-commutative N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, k \cdot \color{blue}{\left(10 \cdot \frac{1}{k} + 1\right)}, 1\right)} \]

      15. distribute-rgt-in N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{\left(10 \cdot \frac{1}{k}\right) \cdot k + 1 \cdot k}, 1\right)} \]

      16. associate-*l* N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{10 \cdot \left(\frac{1}{k} \cdot k\right)} + 1 \cdot k, 1\right)} \]

      17. lft-mult-inverse N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, 10 \cdot \color{blue}{1} + 1 \cdot k, 1\right)} \]

      18. metadata-eval N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{10} + 1 \cdot k, 1\right)} \]

      19. *-lft-identity N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, 10 + \color{blue}{k}, 1\right)} \]

      20. lower-+.f64 94.3

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{10 + k}, 1\right)} \]

   5. Applied rewrites 94.3%

      \[\leadsto \color{blue}{\frac{a}{\mathsf{fma}\left(k, 10 + k, 1\right)}} \]

   if 5.00000000000000041e-6 < m

   1. Initial program 72.0%

      \[\frac{a \cdot {k}^{m}}{\left(1 + 10 \cdot k\right) + k \cdot k} \]
   2. Add Preprocessing

   3. Taylor expanded in k around 0

      \[\leadsto \color{blue}{a \cdot {k}^{m}} \]
   4. Step-by-step derivation

      1. lower-*.f64 N/A

         \[\leadsto \color{blue}{a \cdot {k}^{m}} \]

      2. lower-pow.f64 100.0

         \[\leadsto a \cdot \color{blue}{{k}^{m}} \]

   5. Applied rewrites 100.0%

      \[\leadsto \color{blue}{a \cdot {k}^{m}} \]
3. Recombined 3 regimes into one program.

4. Final simplification 97.9%

   \[\leadsto \begin{array}{l} \mathbf{if}\;m \leq -2 \cdot 10^{-17}:\\ \;\;\;\;\frac{a \cdot {k}^{m}}{k \cdot k}\\ \mathbf{elif}\;m \leq 5 \cdot 10^{-6}:\\ \;\;\;\;\frac{a}{\mathsf{fma}\left(k, k + 10, 1\right)}\\ \mathbf{else}:\\ \;\;\;\;a \cdot {k}^{m}\\ \end{array} \]
5. Add Preprocessing

Alternative 7: 97.1% accurate, 1.1× speedup

Math

\[\begin{array}{l} a\_m = \left|a\right| \\ a\_s = \mathsf{copysign}\left(1, a\right) \\ \begin{array}{l} t_0 := a \cdot {k}^{m}\\ a\_s \cdot \begin{array}{l} \mathbf{if}\;m \leq -50:\\ \;\;\;\;t\_0\\ \mathbf{elif}\;m \leq 5 \cdot 10^{-6}:\\ \;\;\;\;\frac{a}{\mathsf{fma}\left(k, k + 10, 1\right)}\\ \mathbf{else}:\\ \;\;\;\;t\_0\\ \end{array} \end{array} \end{array} \]

FPCore

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 (pow k m))))
   (*
    a_s
    (if (<= m -50.0) t_0 (if (<= m 5e-6) (/ a (fma k (+ k 10.0) 1.0)) t_0)))))

C

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 * pow(k, m);
	double tmp;
	if (m <= -50.0) {
		tmp = t_0;
	} else if (m <= 5e-6) {
		tmp = a / fma(k, (k + 10.0), 1.0);
	} else {
		tmp = t_0;
	}
	return a_s * tmp;
}

Julia

a\_m = abs(a)
a\_s = copysign(1.0, a)
function code(a_s, a_m, k, m)
	t_0 = Float64(a * (k ^ m))
	tmp = 0.0
	if (m <= -50.0)
		tmp = t_0;
	elseif (m <= 5e-6)
		tmp = Float64(a / fma(k, Float64(k + 10.0), 1.0));
	else
		tmp = t_0;
	end
	return Float64(a_s * tmp)
end

Wolfram

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 * N[Power[k, m], $MachinePrecision]), $MachinePrecision]}, N[(a$95$s * If[LessEqual[m, -50.0], t$95$0, If[LessEqual[m, 5e-6], N[(a / N[(k * N[(k + 10.0), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision], t$95$0]]), $MachinePrecision]]

Derivation

1. Split input into 2 regimes

2. if m < -50 or 5.00000000000000041e-6 < m

   1. Initial program 84.9%

      \[\frac{a \cdot {k}^{m}}{\left(1 + 10 \cdot k\right) + k \cdot k} \]
   2. Add Preprocessing

   3. Taylor expanded in k around 0

      \[\leadsto \color{blue}{a \cdot {k}^{m}} \]
   4. Step-by-step derivation

      1. lower-*.f64 N/A

         \[\leadsto \color{blue}{a \cdot {k}^{m}} \]

      2. lower-pow.f64 100.0

         \[\leadsto a \cdot \color{blue}{{k}^{m}} \]

   5. Applied rewrites 100.0%

      \[\leadsto \color{blue}{a \cdot {k}^{m}} \]

   if -50 < m < 5.00000000000000041e-6

   1. Initial program 94.4%

      \[\frac{a \cdot {k}^{m}}{\left(1 + 10 \cdot k\right) + k \cdot k} \]
   2. Add Preprocessing

   3. Taylor expanded in m around 0

      \[\leadsto \color{blue}{\frac{a}{1 + \left(10 \cdot k + {k}^{2}\right)}} \]
   4. Step-by-step derivation

      1. lower-/.f64 N/A

         \[\leadsto \color{blue}{\frac{a}{1 + \left(10 \cdot k + {k}^{2}\right)}} \]

      2. unpow2 N/A

         \[\leadsto \frac{a}{1 + \left(10 \cdot k + \color{blue}{k \cdot k}\right)} \]

      3. distribute-rgt-in N/A

         \[\leadsto \frac{a}{1 + \color{blue}{k \cdot \left(10 + k\right)}} \]

      4. +-commutative N/A

         \[\leadsto \frac{a}{\color{blue}{k \cdot \left(10 + k\right) + 1}} \]

      5. metadata-eval N/A

         \[\leadsto \frac{a}{k \cdot \left(\color{blue}{10 \cdot 1} + k\right) + 1} \]

      6. lft-mult-inverse N/A

         \[\leadsto \frac{a}{k \cdot \left(10 \cdot \color{blue}{\left(\frac{1}{k} \cdot k\right)} + k\right) + 1} \]

      7. associate-*l* N/A

         \[\leadsto \frac{a}{k \cdot \left(\color{blue}{\left(10 \cdot \frac{1}{k}\right) \cdot k} + k\right) + 1} \]

      8. *-lft-identity N/A

         \[\leadsto \frac{a}{k \cdot \left(\left(10 \cdot \frac{1}{k}\right) \cdot k + \color{blue}{1 \cdot k}\right) + 1} \]

      9. distribute-rgt-in N/A

         \[\leadsto \frac{a}{k \cdot \color{blue}{\left(k \cdot \left(10 \cdot \frac{1}{k} + 1\right)\right)} + 1} \]

      10. +-commutative N/A

          \[\leadsto \frac{a}{k \cdot \left(k \cdot \color{blue}{\left(1 + 10 \cdot \frac{1}{k}\right)}\right) + 1} \]

      11. *-commutative N/A

          \[\leadsto \frac{a}{k \cdot \color{blue}{\left(\left(1 + 10 \cdot \frac{1}{k}\right) \cdot k\right)} + 1} \]

      12. lower-fma.f64 N/A

          \[\leadsto \frac{a}{\color{blue}{\mathsf{fma}\left(k, \left(1 + 10 \cdot \frac{1}{k}\right) \cdot k, 1\right)}} \]

      13. *-commutative N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{k \cdot \left(1 + 10 \cdot \frac{1}{k}\right)}, 1\right)} \]

      14. +-commutative N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, k \cdot \color{blue}{\left(10 \cdot \frac{1}{k} + 1\right)}, 1\right)} \]

      15. distribute-rgt-in N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{\left(10 \cdot \frac{1}{k}\right) \cdot k + 1 \cdot k}, 1\right)} \]

      16. associate-*l* N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{10 \cdot \left(\frac{1}{k} \cdot k\right)} + 1 \cdot k, 1\right)} \]

      17. lft-mult-inverse N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, 10 \cdot \color{blue}{1} + 1 \cdot k, 1\right)} \]

      18. metadata-eval N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{10} + 1 \cdot k, 1\right)} \]

      19. *-lft-identity N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, 10 + \color{blue}{k}, 1\right)} \]

      20. lower-+.f64 94.1

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{10 + k}, 1\right)} \]

   5. Applied rewrites 94.1%

      \[\leadsto \color{blue}{\frac{a}{\mathsf{fma}\left(k, 10 + k, 1\right)}} \]
3. Recombined 2 regimes into one program.

4. Final simplification 97.8%

   \[\leadsto \begin{array}{l} \mathbf{if}\;m \leq -50:\\ \;\;\;\;a \cdot {k}^{m}\\ \mathbf{elif}\;m \leq 5 \cdot 10^{-6}:\\ \;\;\;\;\frac{a}{\mathsf{fma}\left(k, k + 10, 1\right)}\\ \mathbf{else}:\\ \;\;\;\;a \cdot {k}^{m}\\ \end{array} \]
5. Add Preprocessing

Alternative 8: 71.3% accurate, 2.3× speedup

Math

\[\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 -0.34:\\ \;\;\;\;\frac{a + \frac{\mathsf{fma}\left(a, -10, \frac{a \cdot 99}{k}\right)}{k}}{k \cdot k}\\ \mathbf{elif}\;m \leq 210000:\\ \;\;\;\;\frac{a}{\mathsf{fma}\left(k, k + 10, 1\right)}\\ \mathbf{else}:\\ \;\;\;\;a \cdot \left(\left(k \cdot \left(k \cdot k\right)\right) \cdot -980\right)\\ \end{array} \end{array} \]

FPCore

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 -0.34)
    (/ (+ a (/ (fma a -10.0 (/ (* a 99.0) k)) k)) (* k k))
    (if (<= m 210000.0)
      (/ a (fma k (+ k 10.0) 1.0))
      (* a (* (* k (* k k)) -980.0))))))

C

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 <= -0.34) {
		tmp = (a + (fma(a, -10.0, ((a * 99.0) / k)) / k)) / (k * k);
	} else if (m <= 210000.0) {
		tmp = a / fma(k, (k + 10.0), 1.0);
	} else {
		tmp = a * ((k * (k * k)) * -980.0);
	}
	return a_s * tmp;
}

Julia

a\_m = abs(a)
a\_s = copysign(1.0, a)
function code(a_s, a_m, k, m)
	tmp = 0.0
	if (m <= -0.34)
		tmp = Float64(Float64(a + Float64(fma(a, -10.0, Float64(Float64(a * 99.0) / k)) / k)) / Float64(k * k));
	elseif (m <= 210000.0)
		tmp = Float64(a / fma(k, Float64(k + 10.0), 1.0));
	else
		tmp = Float64(a * Float64(Float64(k * Float64(k * k)) * -980.0));
	end
	return Float64(a_s * tmp)
end

Wolfram

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, -0.34], N[(N[(a + N[(N[(a * -10.0 + N[(N[(a * 99.0), $MachinePrecision] / k), $MachinePrecision]), $MachinePrecision] / k), $MachinePrecision]), $MachinePrecision] / N[(k * k), $MachinePrecision]), $MachinePrecision], If[LessEqual[m, 210000.0], N[(a / N[(k * N[(k + 10.0), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision], N[(a * N[(N[(k * N[(k * k), $MachinePrecision]), $MachinePrecision] * -980.0), $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]

Derivation

1. Split input into 3 regimes

2. if m < -0.340000000000000024

   1. Initial program 98.7%

      \[\frac{a \cdot {k}^{m}}{\left(1 + 10 \cdot k\right) + k \cdot k} \]
   2. Add Preprocessing

   3. Taylor expanded in m around 0

      \[\leadsto \color{blue}{\frac{a}{1 + \left(10 \cdot k + {k}^{2}\right)}} \]
   4. Step-by-step derivation

      1. lower-/.f64 N/A

         \[\leadsto \color{blue}{\frac{a}{1 + \left(10 \cdot k + {k}^{2}\right)}} \]

      2. unpow2 N/A

         \[\leadsto \frac{a}{1 + \left(10 \cdot k + \color{blue}{k \cdot k}\right)} \]

      3. distribute-rgt-in N/A

         \[\leadsto \frac{a}{1 + \color{blue}{k \cdot \left(10 + k\right)}} \]

      4. +-commutative N/A

         \[\leadsto \frac{a}{\color{blue}{k \cdot \left(10 + k\right) + 1}} \]

      5. metadata-eval N/A

         \[\leadsto \frac{a}{k \cdot \left(\color{blue}{10 \cdot 1} + k\right) + 1} \]

      6. lft-mult-inverse N/A

         \[\leadsto \frac{a}{k \cdot \left(10 \cdot \color{blue}{\left(\frac{1}{k} \cdot k\right)} + k\right) + 1} \]

      7. associate-*l* N/A

         \[\leadsto \frac{a}{k \cdot \left(\color{blue}{\left(10 \cdot \frac{1}{k}\right) \cdot k} + k\right) + 1} \]

      8. *-lft-identity N/A

         \[\leadsto \frac{a}{k \cdot \left(\left(10 \cdot \frac{1}{k}\right) \cdot k + \color{blue}{1 \cdot k}\right) + 1} \]

      9. distribute-rgt-in N/A

         \[\leadsto \frac{a}{k \cdot \color{blue}{\left(k \cdot \left(10 \cdot \frac{1}{k} + 1\right)\right)} + 1} \]

      10. +-commutative N/A

          \[\leadsto \frac{a}{k \cdot \left(k \cdot \color{blue}{\left(1 + 10 \cdot \frac{1}{k}\right)}\right) + 1} \]

      11. *-commutative N/A

          \[\leadsto \frac{a}{k \cdot \color{blue}{\left(\left(1 + 10 \cdot \frac{1}{k}\right) \cdot k\right)} + 1} \]

      12. lower-fma.f64 N/A

          \[\leadsto \frac{a}{\color{blue}{\mathsf{fma}\left(k, \left(1 + 10 \cdot \frac{1}{k}\right) \cdot k, 1\right)}} \]

      13. *-commutative N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{k \cdot \left(1 + 10 \cdot \frac{1}{k}\right)}, 1\right)} \]

      14. +-commutative N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, k \cdot \color{blue}{\left(10 \cdot \frac{1}{k} + 1\right)}, 1\right)} \]

      15. distribute-rgt-in N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{\left(10 \cdot \frac{1}{k}\right) \cdot k + 1 \cdot k}, 1\right)} \]

      16. associate-*l* N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{10 \cdot \left(\frac{1}{k} \cdot k\right)} + 1 \cdot k, 1\right)} \]

      17. lft-mult-inverse N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, 10 \cdot \color{blue}{1} + 1 \cdot k, 1\right)} \]

      18. metadata-eval N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{10} + 1 \cdot k, 1\right)} \]

      19. *-lft-identity N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, 10 + \color{blue}{k}, 1\right)} \]

      20. lower-+.f64 38.8

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{10 + k}, 1\right)} \]

   5. Applied rewrites 38.8%

      \[\leadsto \color{blue}{\frac{a}{\mathsf{fma}\left(k, 10 + k, 1\right)}} \]

   6. Taylor expanded in k around inf

      \[\leadsto \frac{\left(a + -1 \cdot \frac{a + -100 \cdot a}{{k}^{2}}\right) - 10 \cdot \frac{a}{k}}{\color{blue}{{k}^{2}}} \]

   8. Applied rewrites 67.6%

      \[\leadsto \frac{a + \frac{\mathsf{fma}\left(a, -10, \frac{a \cdot 99}{k}\right)}{k}}{\color{blue}{k \cdot k}} \]

   if -0.340000000000000024 < m < 2.1e5

   1. Initial program 93.6%

      \[\frac{a \cdot {k}^{m}}{\left(1 + 10 \cdot k\right) + k \cdot k} \]
   2. Add Preprocessing

   3. Taylor expanded in m around 0

      \[\leadsto \color{blue}{\frac{a}{1 + \left(10 \cdot k + {k}^{2}\right)}} \]
   4. Step-by-step derivation

      1. lower-/.f64 N/A

         \[\leadsto \color{blue}{\frac{a}{1 + \left(10 \cdot k + {k}^{2}\right)}} \]

      2. unpow2 N/A

         \[\leadsto \frac{a}{1 + \left(10 \cdot k + \color{blue}{k \cdot k}\right)} \]

      3. distribute-rgt-in N/A

         \[\leadsto \frac{a}{1 + \color{blue}{k \cdot \left(10 + k\right)}} \]

      4. +-commutative N/A

         \[\leadsto \frac{a}{\color{blue}{k \cdot \left(10 + k\right) + 1}} \]

      5. metadata-eval N/A

         \[\leadsto \frac{a}{k \cdot \left(\color{blue}{10 \cdot 1} + k\right) + 1} \]

      6. lft-mult-inverse N/A

         \[\leadsto \frac{a}{k \cdot \left(10 \cdot \color{blue}{\left(\frac{1}{k} \cdot k\right)} + k\right) + 1} \]

      7. associate-*l* N/A

         \[\leadsto \frac{a}{k \cdot \left(\color{blue}{\left(10 \cdot \frac{1}{k}\right) \cdot k} + k\right) + 1} \]

      8. *-lft-identity N/A

         \[\leadsto \frac{a}{k \cdot \left(\left(10 \cdot \frac{1}{k}\right) \cdot k + \color{blue}{1 \cdot k}\right) + 1} \]

      9. distribute-rgt-in N/A

         \[\leadsto \frac{a}{k \cdot \color{blue}{\left(k \cdot \left(10 \cdot \frac{1}{k} + 1\right)\right)} + 1} \]

      10. +-commutative N/A

          \[\leadsto \frac{a}{k \cdot \left(k \cdot \color{blue}{\left(1 + 10 \cdot \frac{1}{k}\right)}\right) + 1} \]

      11. *-commutative N/A

          \[\leadsto \frac{a}{k \cdot \color{blue}{\left(\left(1 + 10 \cdot \frac{1}{k}\right) \cdot k\right)} + 1} \]

      12. lower-fma.f64 N/A

          \[\leadsto \frac{a}{\color{blue}{\mathsf{fma}\left(k, \left(1 + 10 \cdot \frac{1}{k}\right) \cdot k, 1\right)}} \]

      13. *-commutative N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{k \cdot \left(1 + 10 \cdot \frac{1}{k}\right)}, 1\right)} \]

      14. +-commutative N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, k \cdot \color{blue}{\left(10 \cdot \frac{1}{k} + 1\right)}, 1\right)} \]

      15. distribute-rgt-in N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{\left(10 \cdot \frac{1}{k}\right) \cdot k + 1 \cdot k}, 1\right)} \]

      16. associate-*l* N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{10 \cdot \left(\frac{1}{k} \cdot k\right)} + 1 \cdot k, 1\right)} \]

      17. lft-mult-inverse N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, 10 \cdot \color{blue}{1} + 1 \cdot k, 1\right)} \]

      18. metadata-eval N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{10} + 1 \cdot k, 1\right)} \]

      19. *-lft-identity N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, 10 + \color{blue}{k}, 1\right)} \]

      20. lower-+.f64 91.5

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{10 + k}, 1\right)} \]

   5. Applied rewrites 91.5%

      \[\leadsto \color{blue}{\frac{a}{\mathsf{fma}\left(k, 10 + k, 1\right)}} \]

   if 2.1e5 < m

   1. Initial program 72.2%

      \[\frac{a \cdot {k}^{m}}{\left(1 + 10 \cdot k\right) + k \cdot k} \]
   2. Add Preprocessing

   3. Taylor expanded in m around 0

      \[\leadsto \color{blue}{\frac{a}{1 + \left(10 \cdot k + {k}^{2}\right)}} \]
   4. Step-by-step derivation

      1. lower-/.f64 N/A

         \[\leadsto \color{blue}{\frac{a}{1 + \left(10 \cdot k + {k}^{2}\right)}} \]

      2. unpow2 N/A

         \[\leadsto \frac{a}{1 + \left(10 \cdot k + \color{blue}{k \cdot k}\right)} \]

      3. distribute-rgt-in N/A

         \[\leadsto \frac{a}{1 + \color{blue}{k \cdot \left(10 + k\right)}} \]

      4. +-commutative N/A

         \[\leadsto \frac{a}{\color{blue}{k \cdot \left(10 + k\right) + 1}} \]

      5. metadata-eval N/A

         \[\leadsto \frac{a}{k \cdot \left(\color{blue}{10 \cdot 1} + k\right) + 1} \]

      6. lft-mult-inverse N/A

         \[\leadsto \frac{a}{k \cdot \left(10 \cdot \color{blue}{\left(\frac{1}{k} \cdot k\right)} + k\right) + 1} \]

      7. associate-*l* N/A

         \[\leadsto \frac{a}{k \cdot \left(\color{blue}{\left(10 \cdot \frac{1}{k}\right) \cdot k} + k\right) + 1} \]

      8. *-lft-identity N/A

         \[\leadsto \frac{a}{k \cdot \left(\left(10 \cdot \frac{1}{k}\right) \cdot k + \color{blue}{1 \cdot k}\right) + 1} \]

      9. distribute-rgt-in N/A

         \[\leadsto \frac{a}{k \cdot \color{blue}{\left(k \cdot \left(10 \cdot \frac{1}{k} + 1\right)\right)} + 1} \]

      10. +-commutative N/A

          \[\leadsto \frac{a}{k \cdot \left(k \cdot \color{blue}{\left(1 + 10 \cdot \frac{1}{k}\right)}\right) + 1} \]

      11. *-commutative N/A

          \[\leadsto \frac{a}{k \cdot \color{blue}{\left(\left(1 + 10 \cdot \frac{1}{k}\right) \cdot k\right)} + 1} \]

      12. lower-fma.f64 N/A

          \[\leadsto \frac{a}{\color{blue}{\mathsf{fma}\left(k, \left(1 + 10 \cdot \frac{1}{k}\right) \cdot k, 1\right)}} \]

      13. *-commutative N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{k \cdot \left(1 + 10 \cdot \frac{1}{k}\right)}, 1\right)} \]

      14. +-commutative N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, k \cdot \color{blue}{\left(10 \cdot \frac{1}{k} + 1\right)}, 1\right)} \]

      15. distribute-rgt-in N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{\left(10 \cdot \frac{1}{k}\right) \cdot k + 1 \cdot k}, 1\right)} \]

      16. associate-*l* N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{10 \cdot \left(\frac{1}{k} \cdot k\right)} + 1 \cdot k, 1\right)} \]

      17. lft-mult-inverse N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, 10 \cdot \color{blue}{1} + 1 \cdot k, 1\right)} \]

      18. metadata-eval N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{10} + 1 \cdot k, 1\right)} \]

      19. *-lft-identity N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, 10 + \color{blue}{k}, 1\right)} \]

      20. lower-+.f64 2.7

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{10 + k}, 1\right)} \]

   5. Applied rewrites 2.7%

      \[\leadsto \color{blue}{\frac{a}{\mathsf{fma}\left(k, 10 + k, 1\right)}} \]

   6. Taylor expanded in k around 0

      \[\leadsto a + \color{blue}{k \cdot \left(k \cdot \left(-1 \cdot \left(k \cdot \left(-10 \cdot a + -10 \cdot \left(a + -100 \cdot a\right)\right)\right) - \left(a + -100 \cdot a\right)\right) - 10 \cdot a\right)} \]
   7. Step-by-step derivation

   8. Applied rewrites 14.8%

      \[\leadsto \mathsf{fma}\left(k, \color{blue}{\mathsf{fma}\left(k, \mathsf{fma}\left(k, 10 \cdot \left(-98 \cdot a\right), a \cdot 99\right), a \cdot -10\right)}, a\right) \]

   9. Taylor expanded in k around inf

      \[\leadsto \mathsf{fma}\left(k, -980 \cdot \left(a \cdot \color{blue}{{k}^{2}}\right), a\right) \]
   10. Step-by-step derivation

   11. Applied rewrites 14.8%

       \[\leadsto \mathsf{fma}\left(k, a \cdot \left(\left(k \cdot k\right) \cdot \color{blue}{-980}\right), a\right) \]

   12. Taylor expanded in k around inf

       \[\leadsto -980 \cdot \left(a \cdot \color{blue}{{k}^{3}}\right) \]
   13. Step-by-step derivation

   14. Applied rewrites 49.9%

       \[\leadsto a \cdot \left(\left(k \cdot \left(k \cdot k\right)\right) \cdot \color{blue}{-980}\right) \]
3. Recombined 3 regimes into one program.

4. Final simplification 71.5%

   \[\leadsto \begin{array}{l} \mathbf{if}\;m \leq -0.34:\\ \;\;\;\;\frac{a + \frac{\mathsf{fma}\left(a, -10, \frac{a \cdot 99}{k}\right)}{k}}{k \cdot k}\\ \mathbf{elif}\;m \leq 210000:\\ \;\;\;\;\frac{a}{\mathsf{fma}\left(k, k + 10, 1\right)}\\ \mathbf{else}:\\ \;\;\;\;a \cdot \left(\left(k \cdot \left(k \cdot k\right)\right) \cdot -980\right)\\ \end{array} \]
5. Add Preprocessing

Alternative 9: 69.2% accurate, 4.1× speedup

Math

\[\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 -0.34:\\ \;\;\;\;\frac{a}{k \cdot k}\\ \mathbf{elif}\;m \leq 210000:\\ \;\;\;\;\frac{a}{\mathsf{fma}\left(k, k + 10, 1\right)}\\ \mathbf{else}:\\ \;\;\;\;a \cdot \left(\left(k \cdot \left(k \cdot k\right)\right) \cdot -980\right)\\ \end{array} \end{array} \]

FPCore

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 -0.34)
    (/ a (* k k))
    (if (<= m 210000.0)
      (/ a (fma k (+ k 10.0) 1.0))
      (* a (* (* k (* k k)) -980.0))))))

C

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 <= -0.34) {
		tmp = a / (k * k);
	} else if (m <= 210000.0) {
		tmp = a / fma(k, (k + 10.0), 1.0);
	} else {
		tmp = a * ((k * (k * k)) * -980.0);
	}
	return a_s * tmp;
}

Julia

a\_m = abs(a)
a\_s = copysign(1.0, a)
function code(a_s, a_m, k, m)
	tmp = 0.0
	if (m <= -0.34)
		tmp = Float64(a / Float64(k * k));
	elseif (m <= 210000.0)
		tmp = Float64(a / fma(k, Float64(k + 10.0), 1.0));
	else
		tmp = Float64(a * Float64(Float64(k * Float64(k * k)) * -980.0));
	end
	return Float64(a_s * tmp)
end

Wolfram

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, -0.34], N[(a / N[(k * k), $MachinePrecision]), $MachinePrecision], If[LessEqual[m, 210000.0], N[(a / N[(k * N[(k + 10.0), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision], N[(a * N[(N[(k * N[(k * k), $MachinePrecision]), $MachinePrecision] * -980.0), $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]

Derivation

1. Split input into 3 regimes

2. if m < -0.340000000000000024

   1. Initial program 98.7%

      \[\frac{a \cdot {k}^{m}}{\left(1 + 10 \cdot k\right) + k \cdot k} \]
   2. Add Preprocessing

   3. Taylor expanded in m around 0

      \[\leadsto \color{blue}{\frac{a}{1 + \left(10 \cdot k + {k}^{2}\right)}} \]
   4. Step-by-step derivation

      1. lower-/.f64 N/A

         \[\leadsto \color{blue}{\frac{a}{1 + \left(10 \cdot k + {k}^{2}\right)}} \]

      2. unpow2 N/A

         \[\leadsto \frac{a}{1 + \left(10 \cdot k + \color{blue}{k \cdot k}\right)} \]

      3. distribute-rgt-in N/A

         \[\leadsto \frac{a}{1 + \color{blue}{k \cdot \left(10 + k\right)}} \]

      4. +-commutative N/A

         \[\leadsto \frac{a}{\color{blue}{k \cdot \left(10 + k\right) + 1}} \]

      5. metadata-eval N/A

         \[\leadsto \frac{a}{k \cdot \left(\color{blue}{10 \cdot 1} + k\right) + 1} \]

      6. lft-mult-inverse N/A

         \[\leadsto \frac{a}{k \cdot \left(10 \cdot \color{blue}{\left(\frac{1}{k} \cdot k\right)} + k\right) + 1} \]

      7. associate-*l* N/A

         \[\leadsto \frac{a}{k \cdot \left(\color{blue}{\left(10 \cdot \frac{1}{k}\right) \cdot k} + k\right) + 1} \]

      8. *-lft-identity N/A

         \[\leadsto \frac{a}{k \cdot \left(\left(10 \cdot \frac{1}{k}\right) \cdot k + \color{blue}{1 \cdot k}\right) + 1} \]

      9. distribute-rgt-in N/A

         \[\leadsto \frac{a}{k \cdot \color{blue}{\left(k \cdot \left(10 \cdot \frac{1}{k} + 1\right)\right)} + 1} \]

      10. +-commutative N/A

          \[\leadsto \frac{a}{k \cdot \left(k \cdot \color{blue}{\left(1 + 10 \cdot \frac{1}{k}\right)}\right) + 1} \]

      11. *-commutative N/A

          \[\leadsto \frac{a}{k \cdot \color{blue}{\left(\left(1 + 10 \cdot \frac{1}{k}\right) \cdot k\right)} + 1} \]

      12. lower-fma.f64 N/A

          \[\leadsto \frac{a}{\color{blue}{\mathsf{fma}\left(k, \left(1 + 10 \cdot \frac{1}{k}\right) \cdot k, 1\right)}} \]

      13. *-commutative N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{k \cdot \left(1 + 10 \cdot \frac{1}{k}\right)}, 1\right)} \]

      14. +-commutative N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, k \cdot \color{blue}{\left(10 \cdot \frac{1}{k} + 1\right)}, 1\right)} \]

      15. distribute-rgt-in N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{\left(10 \cdot \frac{1}{k}\right) \cdot k + 1 \cdot k}, 1\right)} \]

      16. associate-*l* N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{10 \cdot \left(\frac{1}{k} \cdot k\right)} + 1 \cdot k, 1\right)} \]

      17. lft-mult-inverse N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, 10 \cdot \color{blue}{1} + 1 \cdot k, 1\right)} \]

      18. metadata-eval N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{10} + 1 \cdot k, 1\right)} \]

      19. *-lft-identity N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, 10 + \color{blue}{k}, 1\right)} \]

      20. lower-+.f64 38.8

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{10 + k}, 1\right)} \]

   5. Applied rewrites 38.8%

      \[\leadsto \color{blue}{\frac{a}{\mathsf{fma}\left(k, 10 + k, 1\right)}} \]

   6. Taylor expanded in k around inf

      \[\leadsto \frac{a}{{k}^{\color{blue}{2}}} \]
   7. Step-by-step derivation

   8. Applied rewrites 63.8%

      \[\leadsto \frac{a}{k \cdot \color{blue}{k}} \]

   if -0.340000000000000024 < m < 2.1e5

   1. Initial program 93.6%

      \[\frac{a \cdot {k}^{m}}{\left(1 + 10 \cdot k\right) + k \cdot k} \]
   2. Add Preprocessing

   3. Taylor expanded in m around 0

      \[\leadsto \color{blue}{\frac{a}{1 + \left(10 \cdot k + {k}^{2}\right)}} \]
   4. Step-by-step derivation

      1. lower-/.f64 N/A

         \[\leadsto \color{blue}{\frac{a}{1 + \left(10 \cdot k + {k}^{2}\right)}} \]

      2. unpow2 N/A

         \[\leadsto \frac{a}{1 + \left(10 \cdot k + \color{blue}{k \cdot k}\right)} \]

      3. distribute-rgt-in N/A

         \[\leadsto \frac{a}{1 + \color{blue}{k \cdot \left(10 + k\right)}} \]

      4. +-commutative N/A

         \[\leadsto \frac{a}{\color{blue}{k \cdot \left(10 + k\right) + 1}} \]

      5. metadata-eval N/A

         \[\leadsto \frac{a}{k \cdot \left(\color{blue}{10 \cdot 1} + k\right) + 1} \]

      6. lft-mult-inverse N/A

         \[\leadsto \frac{a}{k \cdot \left(10 \cdot \color{blue}{\left(\frac{1}{k} \cdot k\right)} + k\right) + 1} \]

      7. associate-*l* N/A

         \[\leadsto \frac{a}{k \cdot \left(\color{blue}{\left(10 \cdot \frac{1}{k}\right) \cdot k} + k\right) + 1} \]

      8. *-lft-identity N/A

         \[\leadsto \frac{a}{k \cdot \left(\left(10 \cdot \frac{1}{k}\right) \cdot k + \color{blue}{1 \cdot k}\right) + 1} \]

      9. distribute-rgt-in N/A

         \[\leadsto \frac{a}{k \cdot \color{blue}{\left(k \cdot \left(10 \cdot \frac{1}{k} + 1\right)\right)} + 1} \]

      10. +-commutative N/A

          \[\leadsto \frac{a}{k \cdot \left(k \cdot \color{blue}{\left(1 + 10 \cdot \frac{1}{k}\right)}\right) + 1} \]

      11. *-commutative N/A

          \[\leadsto \frac{a}{k \cdot \color{blue}{\left(\left(1 + 10 \cdot \frac{1}{k}\right) \cdot k\right)} + 1} \]

      12. lower-fma.f64 N/A

          \[\leadsto \frac{a}{\color{blue}{\mathsf{fma}\left(k, \left(1 + 10 \cdot \frac{1}{k}\right) \cdot k, 1\right)}} \]

      13. *-commutative N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{k \cdot \left(1 + 10 \cdot \frac{1}{k}\right)}, 1\right)} \]

      14. +-commutative N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, k \cdot \color{blue}{\left(10 \cdot \frac{1}{k} + 1\right)}, 1\right)} \]

      15. distribute-rgt-in N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{\left(10 \cdot \frac{1}{k}\right) \cdot k + 1 \cdot k}, 1\right)} \]

      16. associate-*l* N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{10 \cdot \left(\frac{1}{k} \cdot k\right)} + 1 \cdot k, 1\right)} \]

      17. lft-mult-inverse N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, 10 \cdot \color{blue}{1} + 1 \cdot k, 1\right)} \]

      18. metadata-eval N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{10} + 1 \cdot k, 1\right)} \]

      19. *-lft-identity N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, 10 + \color{blue}{k}, 1\right)} \]

      20. lower-+.f64 91.5

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{10 + k}, 1\right)} \]

   5. Applied rewrites 91.5%

      \[\leadsto \color{blue}{\frac{a}{\mathsf{fma}\left(k, 10 + k, 1\right)}} \]

   if 2.1e5 < m

   1. Initial program 72.2%

      \[\frac{a \cdot {k}^{m}}{\left(1 + 10 \cdot k\right) + k \cdot k} \]
   2. Add Preprocessing

   3. Taylor expanded in m around 0

      \[\leadsto \color{blue}{\frac{a}{1 + \left(10 \cdot k + {k}^{2}\right)}} \]
   4. Step-by-step derivation

      1. lower-/.f64 N/A

         \[\leadsto \color{blue}{\frac{a}{1 + \left(10 \cdot k + {k}^{2}\right)}} \]

      2. unpow2 N/A

         \[\leadsto \frac{a}{1 + \left(10 \cdot k + \color{blue}{k \cdot k}\right)} \]

      3. distribute-rgt-in N/A

         \[\leadsto \frac{a}{1 + \color{blue}{k \cdot \left(10 + k\right)}} \]

      4. +-commutative N/A

         \[\leadsto \frac{a}{\color{blue}{k \cdot \left(10 + k\right) + 1}} \]

      5. metadata-eval N/A

         \[\leadsto \frac{a}{k \cdot \left(\color{blue}{10 \cdot 1} + k\right) + 1} \]

      6. lft-mult-inverse N/A

         \[\leadsto \frac{a}{k \cdot \left(10 \cdot \color{blue}{\left(\frac{1}{k} \cdot k\right)} + k\right) + 1} \]

      7. associate-*l* N/A

         \[\leadsto \frac{a}{k \cdot \left(\color{blue}{\left(10 \cdot \frac{1}{k}\right) \cdot k} + k\right) + 1} \]

      8. *-lft-identity N/A

         \[\leadsto \frac{a}{k \cdot \left(\left(10 \cdot \frac{1}{k}\right) \cdot k + \color{blue}{1 \cdot k}\right) + 1} \]

      9. distribute-rgt-in N/A

         \[\leadsto \frac{a}{k \cdot \color{blue}{\left(k \cdot \left(10 \cdot \frac{1}{k} + 1\right)\right)} + 1} \]

      10. +-commutative N/A

          \[\leadsto \frac{a}{k \cdot \left(k \cdot \color{blue}{\left(1 + 10 \cdot \frac{1}{k}\right)}\right) + 1} \]

      11. *-commutative N/A

          \[\leadsto \frac{a}{k \cdot \color{blue}{\left(\left(1 + 10 \cdot \frac{1}{k}\right) \cdot k\right)} + 1} \]

      12. lower-fma.f64 N/A

          \[\leadsto \frac{a}{\color{blue}{\mathsf{fma}\left(k, \left(1 + 10 \cdot \frac{1}{k}\right) \cdot k, 1\right)}} \]

      13. *-commutative N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{k \cdot \left(1 + 10 \cdot \frac{1}{k}\right)}, 1\right)} \]

      14. +-commutative N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, k \cdot \color{blue}{\left(10 \cdot \frac{1}{k} + 1\right)}, 1\right)} \]

      15. distribute-rgt-in N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{\left(10 \cdot \frac{1}{k}\right) \cdot k + 1 \cdot k}, 1\right)} \]

      16. associate-*l* N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{10 \cdot \left(\frac{1}{k} \cdot k\right)} + 1 \cdot k, 1\right)} \]

      17. lft-mult-inverse N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, 10 \cdot \color{blue}{1} + 1 \cdot k, 1\right)} \]

      18. metadata-eval N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{10} + 1 \cdot k, 1\right)} \]

      19. *-lft-identity N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, 10 + \color{blue}{k}, 1\right)} \]

      20. lower-+.f64 2.7

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{10 + k}, 1\right)} \]

   5. Applied rewrites 2.7%

      \[\leadsto \color{blue}{\frac{a}{\mathsf{fma}\left(k, 10 + k, 1\right)}} \]

   6. Taylor expanded in k around 0

      \[\leadsto a + \color{blue}{k \cdot \left(k \cdot \left(-1 \cdot \left(k \cdot \left(-10 \cdot a + -10 \cdot \left(a + -100 \cdot a\right)\right)\right) - \left(a + -100 \cdot a\right)\right) - 10 \cdot a\right)} \]
   7. Step-by-step derivation

   8. Applied rewrites 14.8%

      \[\leadsto \mathsf{fma}\left(k, \color{blue}{\mathsf{fma}\left(k, \mathsf{fma}\left(k, 10 \cdot \left(-98 \cdot a\right), a \cdot 99\right), a \cdot -10\right)}, a\right) \]

   9. Taylor expanded in k around inf

      \[\leadsto \mathsf{fma}\left(k, -980 \cdot \left(a \cdot \color{blue}{{k}^{2}}\right), a\right) \]
   10. Step-by-step derivation

   11. Applied rewrites 14.8%

       \[\leadsto \mathsf{fma}\left(k, a \cdot \left(\left(k \cdot k\right) \cdot \color{blue}{-980}\right), a\right) \]

   12. Taylor expanded in k around inf

       \[\leadsto -980 \cdot \left(a \cdot \color{blue}{{k}^{3}}\right) \]
   13. Step-by-step derivation

   14. Applied rewrites 49.9%

       \[\leadsto a \cdot \left(\left(k \cdot \left(k \cdot k\right)\right) \cdot \color{blue}{-980}\right) \]
3. Recombined 3 regimes into one program.

4. Final simplification 70.3%

   \[\leadsto \begin{array}{l} \mathbf{if}\;m \leq -0.34:\\ \;\;\;\;\frac{a}{k \cdot k}\\ \mathbf{elif}\;m \leq 210000:\\ \;\;\;\;\frac{a}{\mathsf{fma}\left(k, k + 10, 1\right)}\\ \mathbf{else}:\\ \;\;\;\;a \cdot \left(\left(k \cdot \left(k \cdot k\right)\right) \cdot -980\right)\\ \end{array} \]
5. Add Preprocessing

Alternative 10: 58.3% accurate, 4.1× speedup

Math

\[\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 -3.8 \cdot 10^{-105}:\\ \;\;\;\;\frac{a}{k \cdot k}\\ \mathbf{elif}\;m \leq 210000:\\ \;\;\;\;\frac{a}{\mathsf{fma}\left(k, 10, 1\right)}\\ \mathbf{else}:\\ \;\;\;\;a \cdot \left(\left(k \cdot \left(k \cdot k\right)\right) \cdot -980\right)\\ \end{array} \end{array} \]

FPCore

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 -3.8e-105)
    (/ a (* k k))
    (if (<= m 210000.0)
      (/ a (fma k 10.0 1.0))
      (* a (* (* k (* k k)) -980.0))))))

C

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 <= -3.8e-105) {
		tmp = a / (k * k);
	} else if (m <= 210000.0) {
		tmp = a / fma(k, 10.0, 1.0);
	} else {
		tmp = a * ((k * (k * k)) * -980.0);
	}
	return a_s * tmp;
}

Julia

a\_m = abs(a)
a\_s = copysign(1.0, a)
function code(a_s, a_m, k, m)
	tmp = 0.0
	if (m <= -3.8e-105)
		tmp = Float64(a / Float64(k * k));
	elseif (m <= 210000.0)
		tmp = Float64(a / fma(k, 10.0, 1.0));
	else
		tmp = Float64(a * Float64(Float64(k * Float64(k * k)) * -980.0));
	end
	return Float64(a_s * tmp)
end

Wolfram

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, -3.8e-105], N[(a / N[(k * k), $MachinePrecision]), $MachinePrecision], If[LessEqual[m, 210000.0], N[(a / N[(k * 10.0 + 1.0), $MachinePrecision]), $MachinePrecision], N[(a * N[(N[(k * N[(k * k), $MachinePrecision]), $MachinePrecision] * -980.0), $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]

Derivation

1. Split input into 3 regimes

2. if m < -3.7999999999999998e-105

   1. Initial program 98.0%

      \[\frac{a \cdot {k}^{m}}{\left(1 + 10 \cdot k\right) + k \cdot k} \]
   2. Add Preprocessing

   3. Taylor expanded in m around 0

      \[\leadsto \color{blue}{\frac{a}{1 + \left(10 \cdot k + {k}^{2}\right)}} \]
   4. Step-by-step derivation

      1. lower-/.f64 N/A

         \[\leadsto \color{blue}{\frac{a}{1 + \left(10 \cdot k + {k}^{2}\right)}} \]

      2. unpow2 N/A

         \[\leadsto \frac{a}{1 + \left(10 \cdot k + \color{blue}{k \cdot k}\right)} \]

      3. distribute-rgt-in N/A

         \[\leadsto \frac{a}{1 + \color{blue}{k \cdot \left(10 + k\right)}} \]

      4. +-commutative N/A

         \[\leadsto \frac{a}{\color{blue}{k \cdot \left(10 + k\right) + 1}} \]

      5. metadata-eval N/A

         \[\leadsto \frac{a}{k \cdot \left(\color{blue}{10 \cdot 1} + k\right) + 1} \]

      6. lft-mult-inverse N/A

         \[\leadsto \frac{a}{k \cdot \left(10 \cdot \color{blue}{\left(\frac{1}{k} \cdot k\right)} + k\right) + 1} \]

      7. associate-*l* N/A

         \[\leadsto \frac{a}{k \cdot \left(\color{blue}{\left(10 \cdot \frac{1}{k}\right) \cdot k} + k\right) + 1} \]

      8. *-lft-identity N/A

         \[\leadsto \frac{a}{k \cdot \left(\left(10 \cdot \frac{1}{k}\right) \cdot k + \color{blue}{1 \cdot k}\right) + 1} \]

      9. distribute-rgt-in N/A

         \[\leadsto \frac{a}{k \cdot \color{blue}{\left(k \cdot \left(10 \cdot \frac{1}{k} + 1\right)\right)} + 1} \]

      10. +-commutative N/A

          \[\leadsto \frac{a}{k \cdot \left(k \cdot \color{blue}{\left(1 + 10 \cdot \frac{1}{k}\right)}\right) + 1} \]

      11. *-commutative N/A

          \[\leadsto \frac{a}{k \cdot \color{blue}{\left(\left(1 + 10 \cdot \frac{1}{k}\right) \cdot k\right)} + 1} \]

      12. lower-fma.f64 N/A

          \[\leadsto \frac{a}{\color{blue}{\mathsf{fma}\left(k, \left(1 + 10 \cdot \frac{1}{k}\right) \cdot k, 1\right)}} \]

      13. *-commutative N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{k \cdot \left(1 + 10 \cdot \frac{1}{k}\right)}, 1\right)} \]

      14. +-commutative N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, k \cdot \color{blue}{\left(10 \cdot \frac{1}{k} + 1\right)}, 1\right)} \]

      15. distribute-rgt-in N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{\left(10 \cdot \frac{1}{k}\right) \cdot k + 1 \cdot k}, 1\right)} \]

      16. associate-*l* N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{10 \cdot \left(\frac{1}{k} \cdot k\right)} + 1 \cdot k, 1\right)} \]

      17. lft-mult-inverse N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, 10 \cdot \color{blue}{1} + 1 \cdot k, 1\right)} \]

      18. metadata-eval N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{10} + 1 \cdot k, 1\right)} \]

      19. *-lft-identity N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, 10 + \color{blue}{k}, 1\right)} \]

      20. lower-+.f64 49.7

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{10 + k}, 1\right)} \]

   5. Applied rewrites 49.7%

      \[\leadsto \color{blue}{\frac{a}{\mathsf{fma}\left(k, 10 + k, 1\right)}} \]

   6. Taylor expanded in k around inf

      \[\leadsto \frac{a}{{k}^{\color{blue}{2}}} \]
   7. Step-by-step derivation

   8. Applied rewrites 63.9%

      \[\leadsto \frac{a}{k \cdot \color{blue}{k}} \]

   if -3.7999999999999998e-105 < m < 2.1e5

   1. Initial program 93.2%

      \[\frac{a \cdot {k}^{m}}{\left(1 + 10 \cdot k\right) + k \cdot k} \]
   2. Add Preprocessing

   3. Taylor expanded in m around 0

      \[\leadsto \color{blue}{\frac{a}{1 + \left(10 \cdot k + {k}^{2}\right)}} \]
   4. Step-by-step derivation

      1. lower-/.f64 N/A

         \[\leadsto \color{blue}{\frac{a}{1 + \left(10 \cdot k + {k}^{2}\right)}} \]

      2. unpow2 N/A

         \[\leadsto \frac{a}{1 + \left(10 \cdot k + \color{blue}{k \cdot k}\right)} \]

      3. distribute-rgt-in N/A

         \[\leadsto \frac{a}{1 + \color{blue}{k \cdot \left(10 + k\right)}} \]

      4. +-commutative N/A

         \[\leadsto \frac{a}{\color{blue}{k \cdot \left(10 + k\right) + 1}} \]

      5. metadata-eval N/A

         \[\leadsto \frac{a}{k \cdot \left(\color{blue}{10 \cdot 1} + k\right) + 1} \]

      6. lft-mult-inverse N/A

         \[\leadsto \frac{a}{k \cdot \left(10 \cdot \color{blue}{\left(\frac{1}{k} \cdot k\right)} + k\right) + 1} \]

      7. associate-*l* N/A

         \[\leadsto \frac{a}{k \cdot \left(\color{blue}{\left(10 \cdot \frac{1}{k}\right) \cdot k} + k\right) + 1} \]

      8. *-lft-identity N/A

         \[\leadsto \frac{a}{k \cdot \left(\left(10 \cdot \frac{1}{k}\right) \cdot k + \color{blue}{1 \cdot k}\right) + 1} \]

      9. distribute-rgt-in N/A

         \[\leadsto \frac{a}{k \cdot \color{blue}{\left(k \cdot \left(10 \cdot \frac{1}{k} + 1\right)\right)} + 1} \]

      10. +-commutative N/A

          \[\leadsto \frac{a}{k \cdot \left(k \cdot \color{blue}{\left(1 + 10 \cdot \frac{1}{k}\right)}\right) + 1} \]

      11. *-commutative N/A

          \[\leadsto \frac{a}{k \cdot \color{blue}{\left(\left(1 + 10 \cdot \frac{1}{k}\right) \cdot k\right)} + 1} \]

      12. lower-fma.f64 N/A

          \[\leadsto \frac{a}{\color{blue}{\mathsf{fma}\left(k, \left(1 + 10 \cdot \frac{1}{k}\right) \cdot k, 1\right)}} \]

      13. *-commutative N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{k \cdot \left(1 + 10 \cdot \frac{1}{k}\right)}, 1\right)} \]

      14. +-commutative N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, k \cdot \color{blue}{\left(10 \cdot \frac{1}{k} + 1\right)}, 1\right)} \]

      15. distribute-rgt-in N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{\left(10 \cdot \frac{1}{k}\right) \cdot k + 1 \cdot k}, 1\right)} \]

      16. associate-*l* N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{10 \cdot \left(\frac{1}{k} \cdot k\right)} + 1 \cdot k, 1\right)} \]

      17. lft-mult-inverse N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, 10 \cdot \color{blue}{1} + 1 \cdot k, 1\right)} \]

      18. metadata-eval N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{10} + 1 \cdot k, 1\right)} \]

      19. *-lft-identity N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, 10 + \color{blue}{k}, 1\right)} \]

      20. lower-+.f64 91.0

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{10 + k}, 1\right)} \]

   5. Applied rewrites 91.0%

      \[\leadsto \color{blue}{\frac{a}{\mathsf{fma}\left(k, 10 + k, 1\right)}} \]

   6. Taylor expanded in k around 0

      \[\leadsto \frac{a}{\mathsf{fma}\left(k, 10, 1\right)} \]
   7. Step-by-step derivation

   8. Applied rewrites 62.2%

      \[\leadsto \frac{a}{\mathsf{fma}\left(k, 10, 1\right)} \]

   if 2.1e5 < m

   1. Initial program 72.2%

      \[\frac{a \cdot {k}^{m}}{\left(1 + 10 \cdot k\right) + k \cdot k} \]
   2. Add Preprocessing

   3. Taylor expanded in m around 0

      \[\leadsto \color{blue}{\frac{a}{1 + \left(10 \cdot k + {k}^{2}\right)}} \]
   4. Step-by-step derivation

      1. lower-/.f64 N/A

         \[\leadsto \color{blue}{\frac{a}{1 + \left(10 \cdot k + {k}^{2}\right)}} \]

      2. unpow2 N/A

         \[\leadsto \frac{a}{1 + \left(10 \cdot k + \color{blue}{k \cdot k}\right)} \]

      3. distribute-rgt-in N/A

         \[\leadsto \frac{a}{1 + \color{blue}{k \cdot \left(10 + k\right)}} \]

      4. +-commutative N/A

         \[\leadsto \frac{a}{\color{blue}{k \cdot \left(10 + k\right) + 1}} \]

      5. metadata-eval N/A

         \[\leadsto \frac{a}{k \cdot \left(\color{blue}{10 \cdot 1} + k\right) + 1} \]

      6. lft-mult-inverse N/A

         \[\leadsto \frac{a}{k \cdot \left(10 \cdot \color{blue}{\left(\frac{1}{k} \cdot k\right)} + k\right) + 1} \]

      7. associate-*l* N/A

         \[\leadsto \frac{a}{k \cdot \left(\color{blue}{\left(10 \cdot \frac{1}{k}\right) \cdot k} + k\right) + 1} \]

      8. *-lft-identity N/A

         \[\leadsto \frac{a}{k \cdot \left(\left(10 \cdot \frac{1}{k}\right) \cdot k + \color{blue}{1 \cdot k}\right) + 1} \]

      9. distribute-rgt-in N/A

         \[\leadsto \frac{a}{k \cdot \color{blue}{\left(k \cdot \left(10 \cdot \frac{1}{k} + 1\right)\right)} + 1} \]

      10. +-commutative N/A

          \[\leadsto \frac{a}{k \cdot \left(k \cdot \color{blue}{\left(1 + 10 \cdot \frac{1}{k}\right)}\right) + 1} \]

      11. *-commutative N/A

          \[\leadsto \frac{a}{k \cdot \color{blue}{\left(\left(1 + 10 \cdot \frac{1}{k}\right) \cdot k\right)} + 1} \]

      12. lower-fma.f64 N/A

          \[\leadsto \frac{a}{\color{blue}{\mathsf{fma}\left(k, \left(1 + 10 \cdot \frac{1}{k}\right) \cdot k, 1\right)}} \]

      13. *-commutative N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{k \cdot \left(1 + 10 \cdot \frac{1}{k}\right)}, 1\right)} \]

      14. +-commutative N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, k \cdot \color{blue}{\left(10 \cdot \frac{1}{k} + 1\right)}, 1\right)} \]

      15. distribute-rgt-in N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{\left(10 \cdot \frac{1}{k}\right) \cdot k + 1 \cdot k}, 1\right)} \]

      16. associate-*l* N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{10 \cdot \left(\frac{1}{k} \cdot k\right)} + 1 \cdot k, 1\right)} \]

      17. lft-mult-inverse N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, 10 \cdot \color{blue}{1} + 1 \cdot k, 1\right)} \]

      18. metadata-eval N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{10} + 1 \cdot k, 1\right)} \]

      19. *-lft-identity N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, 10 + \color{blue}{k}, 1\right)} \]

      20. lower-+.f64 2.7

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{10 + k}, 1\right)} \]

   5. Applied rewrites 2.7%

      \[\leadsto \color{blue}{\frac{a}{\mathsf{fma}\left(k, 10 + k, 1\right)}} \]

   6. Taylor expanded in k around 0

      \[\leadsto a + \color{blue}{k \cdot \left(k \cdot \left(-1 \cdot \left(k \cdot \left(-10 \cdot a + -10 \cdot \left(a + -100 \cdot a\right)\right)\right) - \left(a + -100 \cdot a\right)\right) - 10 \cdot a\right)} \]
   7. Step-by-step derivation

   8. Applied rewrites 14.8%

      \[\leadsto \mathsf{fma}\left(k, \color{blue}{\mathsf{fma}\left(k, \mathsf{fma}\left(k, 10 \cdot \left(-98 \cdot a\right), a \cdot 99\right), a \cdot -10\right)}, a\right) \]

   9. Taylor expanded in k around inf

      \[\leadsto \mathsf{fma}\left(k, -980 \cdot \left(a \cdot \color{blue}{{k}^{2}}\right), a\right) \]
   10. Step-by-step derivation

   11. Applied rewrites 14.8%

       \[\leadsto \mathsf{fma}\left(k, a \cdot \left(\left(k \cdot k\right) \cdot \color{blue}{-980}\right), a\right) \]

   12. Taylor expanded in k around inf

       \[\leadsto -980 \cdot \left(a \cdot \color{blue}{{k}^{3}}\right) \]
   13. Step-by-step derivation

   14. Applied rewrites 49.9%

       \[\leadsto a \cdot \left(\left(k \cdot \left(k \cdot k\right)\right) \cdot \color{blue}{-980}\right) \]
3. Recombined 3 regimes into one program.
4. Add Preprocessing

Alternative 11: 42.0% accurate, 4.2× speedup

Math

\[\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 10^{-170}:\\ \;\;\;\;\frac{a}{k \cdot k}\\ \mathbf{elif}\;k \leq 4 \cdot 10^{+19}:\\ \;\;\;\;\frac{a}{\mathsf{fma}\left(k, 10, 1\right)}\\ \mathbf{else}:\\ \;\;\;\;\frac{a}{k \cdot \left(k + 10\right)}\\ \end{array} \end{array} \]

FPCore

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 1e-170)
    (/ a (* k k))
    (if (<= k 4e+19) (/ a (fma k 10.0 1.0)) (/ a (* k (+ k 10.0)))))))

C

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 <= 1e-170) {
		tmp = a / (k * k);
	} else if (k <= 4e+19) {
		tmp = a / fma(k, 10.0, 1.0);
	} else {
		tmp = a / (k * (k + 10.0));
	}
	return a_s * tmp;
}

Julia

a\_m = abs(a)
a\_s = copysign(1.0, a)
function code(a_s, a_m, k, m)
	tmp = 0.0
	if (k <= 1e-170)
		tmp = Float64(a / Float64(k * k));
	elseif (k <= 4e+19)
		tmp = Float64(a / fma(k, 10.0, 1.0));
	else
		tmp = Float64(a / Float64(k * Float64(k + 10.0)));
	end
	return Float64(a_s * tmp)
end

Wolfram

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, 1e-170], N[(a / N[(k * k), $MachinePrecision]), $MachinePrecision], If[LessEqual[k, 4e+19], N[(a / N[(k * 10.0 + 1.0), $MachinePrecision]), $MachinePrecision], N[(a / N[(k * N[(k + 10.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]

Derivation

1. Split input into 3 regimes

2. if k < 9.99999999999999983e-171

   1. Initial program 92.3%

      \[\frac{a \cdot {k}^{m}}{\left(1 + 10 \cdot k\right) + k \cdot k} \]
   2. Add Preprocessing

   3. Taylor expanded in m around 0

      \[\leadsto \color{blue}{\frac{a}{1 + \left(10 \cdot k + {k}^{2}\right)}} \]
   4. Step-by-step derivation

      1. lower-/.f64 N/A

         \[\leadsto \color{blue}{\frac{a}{1 + \left(10 \cdot k + {k}^{2}\right)}} \]

      2. unpow2 N/A

         \[\leadsto \frac{a}{1 + \left(10 \cdot k + \color{blue}{k \cdot k}\right)} \]

      3. distribute-rgt-in N/A

         \[\leadsto \frac{a}{1 + \color{blue}{k \cdot \left(10 + k\right)}} \]

      4. +-commutative N/A

         \[\leadsto \frac{a}{\color{blue}{k \cdot \left(10 + k\right) + 1}} \]

      5. metadata-eval N/A

         \[\leadsto \frac{a}{k \cdot \left(\color{blue}{10 \cdot 1} + k\right) + 1} \]

      6. lft-mult-inverse N/A

         \[\leadsto \frac{a}{k \cdot \left(10 \cdot \color{blue}{\left(\frac{1}{k} \cdot k\right)} + k\right) + 1} \]

      7. associate-*l* N/A

         \[\leadsto \frac{a}{k \cdot \left(\color{blue}{\left(10 \cdot \frac{1}{k}\right) \cdot k} + k\right) + 1} \]

      8. *-lft-identity N/A

         \[\leadsto \frac{a}{k \cdot \left(\left(10 \cdot \frac{1}{k}\right) \cdot k + \color{blue}{1 \cdot k}\right) + 1} \]

      9. distribute-rgt-in N/A

         \[\leadsto \frac{a}{k \cdot \color{blue}{\left(k \cdot \left(10 \cdot \frac{1}{k} + 1\right)\right)} + 1} \]

      10. +-commutative N/A

          \[\leadsto \frac{a}{k \cdot \left(k \cdot \color{blue}{\left(1 + 10 \cdot \frac{1}{k}\right)}\right) + 1} \]

      11. *-commutative N/A

          \[\leadsto \frac{a}{k \cdot \color{blue}{\left(\left(1 + 10 \cdot \frac{1}{k}\right) \cdot k\right)} + 1} \]

      12. lower-fma.f64 N/A

          \[\leadsto \frac{a}{\color{blue}{\mathsf{fma}\left(k, \left(1 + 10 \cdot \frac{1}{k}\right) \cdot k, 1\right)}} \]

      13. *-commutative N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{k \cdot \left(1 + 10 \cdot \frac{1}{k}\right)}, 1\right)} \]

      14. +-commutative N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, k \cdot \color{blue}{\left(10 \cdot \frac{1}{k} + 1\right)}, 1\right)} \]

      15. distribute-rgt-in N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{\left(10 \cdot \frac{1}{k}\right) \cdot k + 1 \cdot k}, 1\right)} \]

      16. associate-*l* N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{10 \cdot \left(\frac{1}{k} \cdot k\right)} + 1 \cdot k, 1\right)} \]

      17. lft-mult-inverse N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, 10 \cdot \color{blue}{1} + 1 \cdot k, 1\right)} \]

      18. metadata-eval N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{10} + 1 \cdot k, 1\right)} \]

      19. *-lft-identity N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, 10 + \color{blue}{k}, 1\right)} \]

      20. lower-+.f64 28.4

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{10 + k}, 1\right)} \]

   5. Applied rewrites 28.4%

      \[\leadsto \color{blue}{\frac{a}{\mathsf{fma}\left(k, 10 + k, 1\right)}} \]

   6. Taylor expanded in k around inf

      \[\leadsto \frac{a}{{k}^{\color{blue}{2}}} \]
   7. Step-by-step derivation

   8. Applied rewrites 34.5%

      \[\leadsto \frac{a}{k \cdot \color{blue}{k}} \]

   if 9.99999999999999983e-171 < k < 4e19

   1. Initial program 99.9%

      \[\frac{a \cdot {k}^{m}}{\left(1 + 10 \cdot k\right) + k \cdot k} \]
   2. Add Preprocessing

   3. Taylor expanded in m around 0

      \[\leadsto \color{blue}{\frac{a}{1 + \left(10 \cdot k + {k}^{2}\right)}} \]
   4. Step-by-step derivation

      1. lower-/.f64 N/A

         \[\leadsto \color{blue}{\frac{a}{1 + \left(10 \cdot k + {k}^{2}\right)}} \]

      2. unpow2 N/A

         \[\leadsto \frac{a}{1 + \left(10 \cdot k + \color{blue}{k \cdot k}\right)} \]

      3. distribute-rgt-in N/A

         \[\leadsto \frac{a}{1 + \color{blue}{k \cdot \left(10 + k\right)}} \]

      4. +-commutative N/A

         \[\leadsto \frac{a}{\color{blue}{k \cdot \left(10 + k\right) + 1}} \]

      5. metadata-eval N/A

         \[\leadsto \frac{a}{k \cdot \left(\color{blue}{10 \cdot 1} + k\right) + 1} \]

      6. lft-mult-inverse N/A

         \[\leadsto \frac{a}{k \cdot \left(10 \cdot \color{blue}{\left(\frac{1}{k} \cdot k\right)} + k\right) + 1} \]

      7. associate-*l* N/A

         \[\leadsto \frac{a}{k \cdot \left(\color{blue}{\left(10 \cdot \frac{1}{k}\right) \cdot k} + k\right) + 1} \]

      8. *-lft-identity N/A

         \[\leadsto \frac{a}{k \cdot \left(\left(10 \cdot \frac{1}{k}\right) \cdot k + \color{blue}{1 \cdot k}\right) + 1} \]

      9. distribute-rgt-in N/A

         \[\leadsto \frac{a}{k \cdot \color{blue}{\left(k \cdot \left(10 \cdot \frac{1}{k} + 1\right)\right)} + 1} \]

      10. +-commutative N/A

          \[\leadsto \frac{a}{k \cdot \left(k \cdot \color{blue}{\left(1 + 10 \cdot \frac{1}{k}\right)}\right) + 1} \]

      11. *-commutative N/A

          \[\leadsto \frac{a}{k \cdot \color{blue}{\left(\left(1 + 10 \cdot \frac{1}{k}\right) \cdot k\right)} + 1} \]

      12. lower-fma.f64 N/A

          \[\leadsto \frac{a}{\color{blue}{\mathsf{fma}\left(k, \left(1 + 10 \cdot \frac{1}{k}\right) \cdot k, 1\right)}} \]

      13. *-commutative N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{k \cdot \left(1 + 10 \cdot \frac{1}{k}\right)}, 1\right)} \]

      14. +-commutative N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, k \cdot \color{blue}{\left(10 \cdot \frac{1}{k} + 1\right)}, 1\right)} \]

      15. distribute-rgt-in N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{\left(10 \cdot \frac{1}{k}\right) \cdot k + 1 \cdot k}, 1\right)} \]

      16. associate-*l* N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{10 \cdot \left(\frac{1}{k} \cdot k\right)} + 1 \cdot k, 1\right)} \]

      17. lft-mult-inverse N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, 10 \cdot \color{blue}{1} + 1 \cdot k, 1\right)} \]

      18. metadata-eval N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{10} + 1 \cdot k, 1\right)} \]

      19. *-lft-identity N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, 10 + \color{blue}{k}, 1\right)} \]

      20. lower-+.f64 62.6

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{10 + k}, 1\right)} \]

   5. Applied rewrites 62.6%

      \[\leadsto \color{blue}{\frac{a}{\mathsf{fma}\left(k, 10 + k, 1\right)}} \]

   6. Taylor expanded in k around 0

      \[\leadsto \frac{a}{\mathsf{fma}\left(k, 10, 1\right)} \]
   7. Step-by-step derivation

   8. Applied rewrites 61.6%

      \[\leadsto \frac{a}{\mathsf{fma}\left(k, 10, 1\right)} \]

   if 4e19 < k

   1. Initial program 78.8%

      \[\frac{a \cdot {k}^{m}}{\left(1 + 10 \cdot k\right) + k \cdot k} \]
   2. Add Preprocessing

   3. Taylor expanded in m around 0

      \[\leadsto \color{blue}{\frac{a}{1 + \left(10 \cdot k + {k}^{2}\right)}} \]
   4. Step-by-step derivation

      1. lower-/.f64 N/A

         \[\leadsto \color{blue}{\frac{a}{1 + \left(10 \cdot k + {k}^{2}\right)}} \]

      2. unpow2 N/A

         \[\leadsto \frac{a}{1 + \left(10 \cdot k + \color{blue}{k \cdot k}\right)} \]

      3. distribute-rgt-in N/A

         \[\leadsto \frac{a}{1 + \color{blue}{k \cdot \left(10 + k\right)}} \]

      4. +-commutative N/A

         \[\leadsto \frac{a}{\color{blue}{k \cdot \left(10 + k\right) + 1}} \]

      5. metadata-eval N/A

         \[\leadsto \frac{a}{k \cdot \left(\color{blue}{10 \cdot 1} + k\right) + 1} \]

      6. lft-mult-inverse N/A

         \[\leadsto \frac{a}{k \cdot \left(10 \cdot \color{blue}{\left(\frac{1}{k} \cdot k\right)} + k\right) + 1} \]

      7. associate-*l* N/A

         \[\leadsto \frac{a}{k \cdot \left(\color{blue}{\left(10 \cdot \frac{1}{k}\right) \cdot k} + k\right) + 1} \]

      8. *-lft-identity N/A

         \[\leadsto \frac{a}{k \cdot \left(\left(10 \cdot \frac{1}{k}\right) \cdot k + \color{blue}{1 \cdot k}\right) + 1} \]

      9. distribute-rgt-in N/A

         \[\leadsto \frac{a}{k \cdot \color{blue}{\left(k \cdot \left(10 \cdot \frac{1}{k} + 1\right)\right)} + 1} \]

      10. +-commutative N/A

          \[\leadsto \frac{a}{k \cdot \left(k \cdot \color{blue}{\left(1 + 10 \cdot \frac{1}{k}\right)}\right) + 1} \]

      11. *-commutative N/A

          \[\leadsto \frac{a}{k \cdot \color{blue}{\left(\left(1 + 10 \cdot \frac{1}{k}\right) \cdot k\right)} + 1} \]

      12. lower-fma.f64 N/A

          \[\leadsto \frac{a}{\color{blue}{\mathsf{fma}\left(k, \left(1 + 10 \cdot \frac{1}{k}\right) \cdot k, 1\right)}} \]

      13. *-commutative N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{k \cdot \left(1 + 10 \cdot \frac{1}{k}\right)}, 1\right)} \]

      14. +-commutative N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, k \cdot \color{blue}{\left(10 \cdot \frac{1}{k} + 1\right)}, 1\right)} \]

      15. distribute-rgt-in N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{\left(10 \cdot \frac{1}{k}\right) \cdot k + 1 \cdot k}, 1\right)} \]

      16. associate-*l* N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{10 \cdot \left(\frac{1}{k} \cdot k\right)} + 1 \cdot k, 1\right)} \]

      17. lft-mult-inverse N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, 10 \cdot \color{blue}{1} + 1 \cdot k, 1\right)} \]

      18. metadata-eval N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{10} + 1 \cdot k, 1\right)} \]

      19. *-lft-identity N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, 10 + \color{blue}{k}, 1\right)} \]

      20. lower-+.f64 61.5

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{10 + k}, 1\right)} \]

   5. Applied rewrites 61.5%

      \[\leadsto \color{blue}{\frac{a}{\mathsf{fma}\left(k, 10 + k, 1\right)}} \]

   6. Taylor expanded in k around inf

      \[\leadsto \frac{a}{{k}^{2} \cdot \color{blue}{\left(1 + 10 \cdot \frac{1}{k}\right)}} \]
   7. Step-by-step derivation

   8. Applied rewrites 61.5%

      \[\leadsto \frac{a}{k \cdot \color{blue}{\left(k + 10\right)}} \]
3. Recombined 3 regimes into one program.
4. Add Preprocessing

Alternative 12: 42.0% accurate, 4.5× speedup

Math

\[\begin{array}{l} a\_m = \left|a\right| \\ a\_s = \mathsf{copysign}\left(1, a\right) \\ \begin{array}{l} t_0 := \frac{a}{k \cdot k}\\ a\_s \cdot \begin{array}{l} \mathbf{if}\;k \leq 10^{-170}:\\ \;\;\;\;t\_0\\ \mathbf{elif}\;k \leq 4 \cdot 10^{+19}:\\ \;\;\;\;\frac{a}{\mathsf{fma}\left(k, 10, 1\right)}\\ \mathbf{else}:\\ \;\;\;\;t\_0\\ \end{array} \end{array} \end{array} \]

FPCore

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 (* k k))))
   (*
    a_s
    (if (<= k 1e-170) t_0 (if (<= k 4e+19) (/ a (fma k 10.0 1.0)) t_0)))))

C

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 / (k * k);
	double tmp;
	if (k <= 1e-170) {
		tmp = t_0;
	} else if (k <= 4e+19) {
		tmp = a / fma(k, 10.0, 1.0);
	} else {
		tmp = t_0;
	}
	return a_s * tmp;
}

Julia

a\_m = abs(a)
a\_s = copysign(1.0, a)
function code(a_s, a_m, k, m)
	t_0 = Float64(a / Float64(k * k))
	tmp = 0.0
	if (k <= 1e-170)
		tmp = t_0;
	elseif (k <= 4e+19)
		tmp = Float64(a / fma(k, 10.0, 1.0));
	else
		tmp = t_0;
	end
	return Float64(a_s * tmp)
end

Wolfram

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 / N[(k * k), $MachinePrecision]), $MachinePrecision]}, N[(a$95$s * If[LessEqual[k, 1e-170], t$95$0, If[LessEqual[k, 4e+19], N[(a / N[(k * 10.0 + 1.0), $MachinePrecision]), $MachinePrecision], t$95$0]]), $MachinePrecision]]

Derivation

1. Split input into 2 regimes

2. if k < 9.99999999999999983e-171 or 4e19 < k

   1. Initial program 85.7%

      \[\frac{a \cdot {k}^{m}}{\left(1 + 10 \cdot k\right) + k \cdot k} \]
   2. Add Preprocessing

   3. Taylor expanded in m around 0

      \[\leadsto \color{blue}{\frac{a}{1 + \left(10 \cdot k + {k}^{2}\right)}} \]
   4. Step-by-step derivation

      1. lower-/.f64 N/A

         \[\leadsto \color{blue}{\frac{a}{1 + \left(10 \cdot k + {k}^{2}\right)}} \]

      2. unpow2 N/A

         \[\leadsto \frac{a}{1 + \left(10 \cdot k + \color{blue}{k \cdot k}\right)} \]

      3. distribute-rgt-in N/A

         \[\leadsto \frac{a}{1 + \color{blue}{k \cdot \left(10 + k\right)}} \]

      4. +-commutative N/A

         \[\leadsto \frac{a}{\color{blue}{k \cdot \left(10 + k\right) + 1}} \]

      5. metadata-eval N/A

         \[\leadsto \frac{a}{k \cdot \left(\color{blue}{10 \cdot 1} + k\right) + 1} \]

      6. lft-mult-inverse N/A

         \[\leadsto \frac{a}{k \cdot \left(10 \cdot \color{blue}{\left(\frac{1}{k} \cdot k\right)} + k\right) + 1} \]

      7. associate-*l* N/A

         \[\leadsto \frac{a}{k \cdot \left(\color{blue}{\left(10 \cdot \frac{1}{k}\right) \cdot k} + k\right) + 1} \]

      8. *-lft-identity N/A

         \[\leadsto \frac{a}{k \cdot \left(\left(10 \cdot \frac{1}{k}\right) \cdot k + \color{blue}{1 \cdot k}\right) + 1} \]

      9. distribute-rgt-in N/A

         \[\leadsto \frac{a}{k \cdot \color{blue}{\left(k \cdot \left(10 \cdot \frac{1}{k} + 1\right)\right)} + 1} \]

      10. +-commutative N/A

          \[\leadsto \frac{a}{k \cdot \left(k \cdot \color{blue}{\left(1 + 10 \cdot \frac{1}{k}\right)}\right) + 1} \]

      11. *-commutative N/A

          \[\leadsto \frac{a}{k \cdot \color{blue}{\left(\left(1 + 10 \cdot \frac{1}{k}\right) \cdot k\right)} + 1} \]

      12. lower-fma.f64 N/A

          \[\leadsto \frac{a}{\color{blue}{\mathsf{fma}\left(k, \left(1 + 10 \cdot \frac{1}{k}\right) \cdot k, 1\right)}} \]

      13. *-commutative N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{k \cdot \left(1 + 10 \cdot \frac{1}{k}\right)}, 1\right)} \]

      14. +-commutative N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, k \cdot \color{blue}{\left(10 \cdot \frac{1}{k} + 1\right)}, 1\right)} \]

      15. distribute-rgt-in N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{\left(10 \cdot \frac{1}{k}\right) \cdot k + 1 \cdot k}, 1\right)} \]

      16. associate-*l* N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{10 \cdot \left(\frac{1}{k} \cdot k\right)} + 1 \cdot k, 1\right)} \]

      17. lft-mult-inverse N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, 10 \cdot \color{blue}{1} + 1 \cdot k, 1\right)} \]

      18. metadata-eval N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{10} + 1 \cdot k, 1\right)} \]

      19. *-lft-identity N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, 10 + \color{blue}{k}, 1\right)} \]

      20. lower-+.f64 44.7

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{10 + k}, 1\right)} \]

   5. Applied rewrites 44.7%

      \[\leadsto \color{blue}{\frac{a}{\mathsf{fma}\left(k, 10 + k, 1\right)}} \]

   6. Taylor expanded in k around inf

      \[\leadsto \frac{a}{{k}^{\color{blue}{2}}} \]
   7. Step-by-step derivation

   8. Applied rewrites 47.8%

      \[\leadsto \frac{a}{k \cdot \color{blue}{k}} \]

   if 9.99999999999999983e-171 < k < 4e19

   1. Initial program 99.9%

      \[\frac{a \cdot {k}^{m}}{\left(1 + 10 \cdot k\right) + k \cdot k} \]
   2. Add Preprocessing

   3. Taylor expanded in m around 0

      \[\leadsto \color{blue}{\frac{a}{1 + \left(10 \cdot k + {k}^{2}\right)}} \]
   4. Step-by-step derivation

      1. lower-/.f64 N/A

         \[\leadsto \color{blue}{\frac{a}{1 + \left(10 \cdot k + {k}^{2}\right)}} \]

      2. unpow2 N/A

         \[\leadsto \frac{a}{1 + \left(10 \cdot k + \color{blue}{k \cdot k}\right)} \]

      3. distribute-rgt-in N/A

         \[\leadsto \frac{a}{1 + \color{blue}{k \cdot \left(10 + k\right)}} \]

      4. +-commutative N/A

         \[\leadsto \frac{a}{\color{blue}{k \cdot \left(10 + k\right) + 1}} \]

      5. metadata-eval N/A

         \[\leadsto \frac{a}{k \cdot \left(\color{blue}{10 \cdot 1} + k\right) + 1} \]

      6. lft-mult-inverse N/A

         \[\leadsto \frac{a}{k \cdot \left(10 \cdot \color{blue}{\left(\frac{1}{k} \cdot k\right)} + k\right) + 1} \]

      7. associate-*l* N/A

         \[\leadsto \frac{a}{k \cdot \left(\color{blue}{\left(10 \cdot \frac{1}{k}\right) \cdot k} + k\right) + 1} \]

      8. *-lft-identity N/A

         \[\leadsto \frac{a}{k \cdot \left(\left(10 \cdot \frac{1}{k}\right) \cdot k + \color{blue}{1 \cdot k}\right) + 1} \]

      9. distribute-rgt-in N/A

         \[\leadsto \frac{a}{k \cdot \color{blue}{\left(k \cdot \left(10 \cdot \frac{1}{k} + 1\right)\right)} + 1} \]

      10. +-commutative N/A

          \[\leadsto \frac{a}{k \cdot \left(k \cdot \color{blue}{\left(1 + 10 \cdot \frac{1}{k}\right)}\right) + 1} \]

      11. *-commutative N/A

          \[\leadsto \frac{a}{k \cdot \color{blue}{\left(\left(1 + 10 \cdot \frac{1}{k}\right) \cdot k\right)} + 1} \]

      12. lower-fma.f64 N/A

          \[\leadsto \frac{a}{\color{blue}{\mathsf{fma}\left(k, \left(1 + 10 \cdot \frac{1}{k}\right) \cdot k, 1\right)}} \]

      13. *-commutative N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{k \cdot \left(1 + 10 \cdot \frac{1}{k}\right)}, 1\right)} \]

      14. +-commutative N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, k \cdot \color{blue}{\left(10 \cdot \frac{1}{k} + 1\right)}, 1\right)} \]

      15. distribute-rgt-in N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{\left(10 \cdot \frac{1}{k}\right) \cdot k + 1 \cdot k}, 1\right)} \]

      16. associate-*l* N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{10 \cdot \left(\frac{1}{k} \cdot k\right)} + 1 \cdot k, 1\right)} \]

      17. lft-mult-inverse N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, 10 \cdot \color{blue}{1} + 1 \cdot k, 1\right)} \]

      18. metadata-eval N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{10} + 1 \cdot k, 1\right)} \]

      19. *-lft-identity N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, 10 + \color{blue}{k}, 1\right)} \]

      20. lower-+.f64 62.6

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{10 + k}, 1\right)} \]

   5. Applied rewrites 62.6%

      \[\leadsto \color{blue}{\frac{a}{\mathsf{fma}\left(k, 10 + k, 1\right)}} \]

   6. Taylor expanded in k around 0

      \[\leadsto \frac{a}{\mathsf{fma}\left(k, 10, 1\right)} \]
   7. Step-by-step derivation

   8. Applied rewrites 61.6%

      \[\leadsto \frac{a}{\mathsf{fma}\left(k, 10, 1\right)} \]
3. Recombined 2 regimes into one program.
4. Add Preprocessing

Alternative 13: 42.5% accurate, 4.6× speedup

Math

\[\begin{array}{l} a\_m = \left|a\right| \\ a\_s = \mathsf{copysign}\left(1, a\right) \\ \begin{array}{l} t_0 := \frac{a}{k \cdot k}\\ a\_s \cdot \begin{array}{l} \mathbf{if}\;k \leq 10^{-170}:\\ \;\;\;\;t\_0\\ \mathbf{elif}\;k \leq 0.001:\\ \;\;\;\;a \cdot \mathsf{fma}\left(k, -10, 1\right)\\ \mathbf{else}:\\ \;\;\;\;t\_0\\ \end{array} \end{array} \end{array} \]

FPCore

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 (* k k))))
   (*
    a_s
    (if (<= k 1e-170) t_0 (if (<= k 0.001) (* a (fma k -10.0 1.0)) t_0)))))

C

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 / (k * k);
	double tmp;
	if (k <= 1e-170) {
		tmp = t_0;
	} else if (k <= 0.001) {
		tmp = a * fma(k, -10.0, 1.0);
	} else {
		tmp = t_0;
	}
	return a_s * tmp;
}

Julia

a\_m = abs(a)
a\_s = copysign(1.0, a)
function code(a_s, a_m, k, m)
	t_0 = Float64(a / Float64(k * k))
	tmp = 0.0
	if (k <= 1e-170)
		tmp = t_0;
	elseif (k <= 0.001)
		tmp = Float64(a * fma(k, -10.0, 1.0));
	else
		tmp = t_0;
	end
	return Float64(a_s * tmp)
end

Wolfram

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 / N[(k * k), $MachinePrecision]), $MachinePrecision]}, N[(a$95$s * If[LessEqual[k, 1e-170], t$95$0, If[LessEqual[k, 0.001], N[(a * N[(k * -10.0 + 1.0), $MachinePrecision]), $MachinePrecision], t$95$0]]), $MachinePrecision]]

Derivation

1. Split input into 2 regimes

2. if k < 9.99999999999999983e-171 or 1e-3 < k

   1. Initial program 85.8%

      \[\frac{a \cdot {k}^{m}}{\left(1 + 10 \cdot k\right) + k \cdot k} \]
   2. Add Preprocessing

   3. Taylor expanded in m around 0

      \[\leadsto \color{blue}{\frac{a}{1 + \left(10 \cdot k + {k}^{2}\right)}} \]
   4. Step-by-step derivation

      1. lower-/.f64 N/A

         \[\leadsto \color{blue}{\frac{a}{1 + \left(10 \cdot k + {k}^{2}\right)}} \]

      2. unpow2 N/A

         \[\leadsto \frac{a}{1 + \left(10 \cdot k + \color{blue}{k \cdot k}\right)} \]

      3. distribute-rgt-in N/A

         \[\leadsto \frac{a}{1 + \color{blue}{k \cdot \left(10 + k\right)}} \]

      4. +-commutative N/A

         \[\leadsto \frac{a}{\color{blue}{k \cdot \left(10 + k\right) + 1}} \]

      5. metadata-eval N/A

         \[\leadsto \frac{a}{k \cdot \left(\color{blue}{10 \cdot 1} + k\right) + 1} \]

      6. lft-mult-inverse N/A

         \[\leadsto \frac{a}{k \cdot \left(10 \cdot \color{blue}{\left(\frac{1}{k} \cdot k\right)} + k\right) + 1} \]

      7. associate-*l* N/A

         \[\leadsto \frac{a}{k \cdot \left(\color{blue}{\left(10 \cdot \frac{1}{k}\right) \cdot k} + k\right) + 1} \]

      8. *-lft-identity N/A

         \[\leadsto \frac{a}{k \cdot \left(\left(10 \cdot \frac{1}{k}\right) \cdot k + \color{blue}{1 \cdot k}\right) + 1} \]

      9. distribute-rgt-in N/A

         \[\leadsto \frac{a}{k \cdot \color{blue}{\left(k \cdot \left(10 \cdot \frac{1}{k} + 1\right)\right)} + 1} \]

      10. +-commutative N/A

          \[\leadsto \frac{a}{k \cdot \left(k \cdot \color{blue}{\left(1 + 10 \cdot \frac{1}{k}\right)}\right) + 1} \]

      11. *-commutative N/A

          \[\leadsto \frac{a}{k \cdot \color{blue}{\left(\left(1 + 10 \cdot \frac{1}{k}\right) \cdot k\right)} + 1} \]

      12. lower-fma.f64 N/A

          \[\leadsto \frac{a}{\color{blue}{\mathsf{fma}\left(k, \left(1 + 10 \cdot \frac{1}{k}\right) \cdot k, 1\right)}} \]

      13. *-commutative N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{k \cdot \left(1 + 10 \cdot \frac{1}{k}\right)}, 1\right)} \]

      14. +-commutative N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, k \cdot \color{blue}{\left(10 \cdot \frac{1}{k} + 1\right)}, 1\right)} \]

      15. distribute-rgt-in N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{\left(10 \cdot \frac{1}{k}\right) \cdot k + 1 \cdot k}, 1\right)} \]

      16. associate-*l* N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{10 \cdot \left(\frac{1}{k} \cdot k\right)} + 1 \cdot k, 1\right)} \]

      17. lft-mult-inverse N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, 10 \cdot \color{blue}{1} + 1 \cdot k, 1\right)} \]

      18. metadata-eval N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{10} + 1 \cdot k, 1\right)} \]

      19. *-lft-identity N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, 10 + \color{blue}{k}, 1\right)} \]

      20. lower-+.f64 44.3

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{10 + k}, 1\right)} \]

   5. Applied rewrites 44.3%

      \[\leadsto \color{blue}{\frac{a}{\mathsf{fma}\left(k, 10 + k, 1\right)}} \]

   6. Taylor expanded in k around inf

      \[\leadsto \frac{a}{{k}^{\color{blue}{2}}} \]
   7. Step-by-step derivation

   8. Applied rewrites 47.4%

      \[\leadsto \frac{a}{k \cdot \color{blue}{k}} \]

   if 9.99999999999999983e-171 < k < 1e-3

   1. Initial program 99.9%

      \[\frac{a \cdot {k}^{m}}{\left(1 + 10 \cdot k\right) + k \cdot k} \]
   2. Add Preprocessing

   3. Taylor expanded in m around 0

      \[\leadsto \color{blue}{\frac{a}{1 + \left(10 \cdot k + {k}^{2}\right)}} \]
   4. Step-by-step derivation

      1. lower-/.f64 N/A

         \[\leadsto \color{blue}{\frac{a}{1 + \left(10 \cdot k + {k}^{2}\right)}} \]

      2. unpow2 N/A

         \[\leadsto \frac{a}{1 + \left(10 \cdot k + \color{blue}{k \cdot k}\right)} \]

      3. distribute-rgt-in N/A

         \[\leadsto \frac{a}{1 + \color{blue}{k \cdot \left(10 + k\right)}} \]

      4. +-commutative N/A

         \[\leadsto \frac{a}{\color{blue}{k \cdot \left(10 + k\right) + 1}} \]

      5. metadata-eval N/A

         \[\leadsto \frac{a}{k \cdot \left(\color{blue}{10 \cdot 1} + k\right) + 1} \]

      6. lft-mult-inverse N/A

         \[\leadsto \frac{a}{k \cdot \left(10 \cdot \color{blue}{\left(\frac{1}{k} \cdot k\right)} + k\right) + 1} \]

      7. associate-*l* N/A

         \[\leadsto \frac{a}{k \cdot \left(\color{blue}{\left(10 \cdot \frac{1}{k}\right) \cdot k} + k\right) + 1} \]

      8. *-lft-identity N/A

         \[\leadsto \frac{a}{k \cdot \left(\left(10 \cdot \frac{1}{k}\right) \cdot k + \color{blue}{1 \cdot k}\right) + 1} \]

      9. distribute-rgt-in N/A

         \[\leadsto \frac{a}{k \cdot \color{blue}{\left(k \cdot \left(10 \cdot \frac{1}{k} + 1\right)\right)} + 1} \]

      10. +-commutative N/A

          \[\leadsto \frac{a}{k \cdot \left(k \cdot \color{blue}{\left(1 + 10 \cdot \frac{1}{k}\right)}\right) + 1} \]

      11. *-commutative N/A

          \[\leadsto \frac{a}{k \cdot \color{blue}{\left(\left(1 + 10 \cdot \frac{1}{k}\right) \cdot k\right)} + 1} \]

      12. lower-fma.f64 N/A

          \[\leadsto \frac{a}{\color{blue}{\mathsf{fma}\left(k, \left(1 + 10 \cdot \frac{1}{k}\right) \cdot k, 1\right)}} \]

      13. *-commutative N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{k \cdot \left(1 + 10 \cdot \frac{1}{k}\right)}, 1\right)} \]

      14. +-commutative N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, k \cdot \color{blue}{\left(10 \cdot \frac{1}{k} + 1\right)}, 1\right)} \]

      15. distribute-rgt-in N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{\left(10 \cdot \frac{1}{k}\right) \cdot k + 1 \cdot k}, 1\right)} \]

      16. associate-*l* N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{10 \cdot \left(\frac{1}{k} \cdot k\right)} + 1 \cdot k, 1\right)} \]

      17. lft-mult-inverse N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, 10 \cdot \color{blue}{1} + 1 \cdot k, 1\right)} \]

      18. metadata-eval N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{10} + 1 \cdot k, 1\right)} \]

      19. *-lft-identity N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, 10 + \color{blue}{k}, 1\right)} \]

      20. lower-+.f64 65.1

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{10 + k}, 1\right)} \]

   5. Applied rewrites 65.1%

      \[\leadsto \color{blue}{\frac{a}{\mathsf{fma}\left(k, 10 + k, 1\right)}} \]

   6. Taylor expanded in k around 0

      \[\leadsto a + \color{blue}{-10 \cdot \left(a \cdot k\right)} \]
   7. Step-by-step derivation

   8. Applied rewrites 63.9%

      \[\leadsto \mathsf{fma}\left(a, \color{blue}{k \cdot -10}, a\right) \]
   9. Step-by-step derivation

   10. Applied rewrites 63.9%

       \[\leadsto \mathsf{fma}\left(k, -10, 1\right) \cdot a \]
3. Recombined 2 regimes into one program.

4. Final simplification 50.6%

   \[\leadsto \begin{array}{l} \mathbf{if}\;k \leq 10^{-170}:\\ \;\;\;\;\frac{a}{k \cdot k}\\ \mathbf{elif}\;k \leq 0.001:\\ \;\;\;\;a \cdot \mathsf{fma}\left(k, -10, 1\right)\\ \mathbf{else}:\\ \;\;\;\;\frac{a}{k \cdot k}\\ \end{array} \]
5. Add Preprocessing

Alternative 14: 25.0% accurate, 7.9× speedup

Math

\[\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 220000:\\ \;\;\;\;a\_m \cdot 1\\ \mathbf{else}:\\ \;\;\;\;a \cdot \left(k \cdot -10\right)\\ \end{array} \end{array} \]

FPCore

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 220000.0) (* a_m 1.0) (* a (* k -10.0)))))

C

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 <= 220000.0) {
		tmp = a_m * 1.0;
	} else {
		tmp = a * (k * -10.0);
	}
	return a_s * tmp;
}

Fortran

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 <= 220000.0d0) then
        tmp = a_m * 1.0d0
    else
        tmp = a * (k * (-10.0d0))
    end if
    code = a_s * tmp
end function

Java

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 <= 220000.0) {
		tmp = a_m * 1.0;
	} else {
		tmp = a * (k * -10.0);
	}
	return a_s * tmp;
}

Python

a\_m = math.fabs(a)
a\_s = math.copysign(1.0, a)
def code(a_s, a_m, k, m):
	tmp = 0
	if m <= 220000.0:
		tmp = a_m * 1.0
	else:
		tmp = a * (k * -10.0)
	return a_s * tmp

Julia

a\_m = abs(a)
a\_s = copysign(1.0, a)
function code(a_s, a_m, k, m)
	tmp = 0.0
	if (m <= 220000.0)
		tmp = Float64(a_m * 1.0);
	else
		tmp = Float64(a * Float64(k * -10.0));
	end
	return Float64(a_s * tmp)
end

MATLAB

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 <= 220000.0)
		tmp = a_m * 1.0;
	else
		tmp = a * (k * -10.0);
	end
	tmp_2 = a_s * tmp;
end

Wolfram

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, 220000.0], N[(a$95$m * 1.0), $MachinePrecision], N[(a * N[(k * -10.0), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]

Derivation

1. Split input into 2 regimes

2. if m < 2.2e5

   1. Initial program 95.8%

      \[\frac{a \cdot {k}^{m}}{\left(1 + 10 \cdot k\right) + k \cdot k} \]
   2. Add Preprocessing
   3. Step-by-step derivation

      1. lift-/.f64 N/A

         \[\leadsto \color{blue}{\frac{a \cdot {k}^{m}}{\left(1 + 10 \cdot k\right) + k \cdot k}} \]

      2. lift-*.f64 N/A

         \[\leadsto \frac{\color{blue}{a \cdot {k}^{m}}}{\left(1 + 10 \cdot k\right) + k \cdot k} \]

      3. associate-/l* N/A

         \[\leadsto \color{blue}{a \cdot \frac{{k}^{m}}{\left(1 + 10 \cdot k\right) + k \cdot k}} \]

      4. *-commutative N/A

         \[\leadsto \color{blue}{\frac{{k}^{m}}{\left(1 + 10 \cdot k\right) + k \cdot k} \cdot a} \]

      5. lower-*.f64 N/A

         \[\leadsto \color{blue}{\frac{{k}^{m}}{\left(1 + 10 \cdot k\right) + k \cdot k} \cdot a} \]

      6. lower-/.f64 95.8

         \[\leadsto \color{blue}{\frac{{k}^{m}}{\left(1 + 10 \cdot k\right) + k \cdot k}} \cdot a \]

      7. lift-+.f64 N/A

         \[\leadsto \frac{{k}^{m}}{\color{blue}{\left(1 + 10 \cdot k\right) + k \cdot k}} \cdot a \]

      8. lift-+.f64 N/A

         \[\leadsto \frac{{k}^{m}}{\color{blue}{\left(1 + 10 \cdot k\right)} + k \cdot k} \cdot a \]

      9. associate-+l+ N/A

         \[\leadsto \frac{{k}^{m}}{\color{blue}{1 + \left(10 \cdot k + k \cdot k\right)}} \cdot a \]

      10. +-commutative N/A

          \[\leadsto \frac{{k}^{m}}{\color{blue}{\left(10 \cdot k + k \cdot k\right) + 1}} \cdot a \]

      11. lift-*.f64 N/A

          \[\leadsto \frac{{k}^{m}}{\left(\color{blue}{10 \cdot k} + k \cdot k\right) + 1} \cdot a \]

      12. lift-*.f64 N/A

          \[\leadsto \frac{{k}^{m}}{\left(10 \cdot k + \color{blue}{k \cdot k}\right) + 1} \cdot a \]

      13. distribute-rgt-out N/A

          \[\leadsto \frac{{k}^{m}}{\color{blue}{k \cdot \left(10 + k\right)} + 1} \cdot a \]

      14. lower-fma.f64 N/A

          \[\leadsto \frac{{k}^{m}}{\color{blue}{\mathsf{fma}\left(k, 10 + k, 1\right)}} \cdot a \]

      15. lower-+.f64 96.3

          \[\leadsto \frac{{k}^{m}}{\mathsf{fma}\left(k, \color{blue}{10 + k}, 1\right)} \cdot a \]

   4. Applied rewrites 96.3%

      \[\leadsto \color{blue}{\frac{{k}^{m}}{\mathsf{fma}\left(k, 10 + k, 1\right)} \cdot a} \]

   5. Taylor expanded in k around 0

      \[\leadsto \color{blue}{{k}^{m}} \cdot a \]
   6. Step-by-step derivation

      1. lower-pow.f64 70.5

         \[\leadsto \color{blue}{{k}^{m}} \cdot a \]

   7. Applied rewrites 70.5%

      \[\leadsto \color{blue}{{k}^{m}} \cdot a \]

   8. Taylor expanded in m around 0

      \[\leadsto 1 \cdot a \]
   9. Step-by-step derivation

   10. Applied rewrites 27.1%

       \[\leadsto 1 \cdot a \]

   if 2.2e5 < m

   1. Initial program 72.2%

      \[\frac{a \cdot {k}^{m}}{\left(1 + 10 \cdot k\right) + k \cdot k} \]
   2. Add Preprocessing

   3. Taylor expanded in m around 0

      \[\leadsto \color{blue}{\frac{a}{1 + \left(10 \cdot k + {k}^{2}\right)}} \]
   4. Step-by-step derivation

      1. lower-/.f64 N/A

         \[\leadsto \color{blue}{\frac{a}{1 + \left(10 \cdot k + {k}^{2}\right)}} \]

      2. unpow2 N/A

         \[\leadsto \frac{a}{1 + \left(10 \cdot k + \color{blue}{k \cdot k}\right)} \]

      3. distribute-rgt-in N/A

         \[\leadsto \frac{a}{1 + \color{blue}{k \cdot \left(10 + k\right)}} \]

      4. +-commutative N/A

         \[\leadsto \frac{a}{\color{blue}{k \cdot \left(10 + k\right) + 1}} \]

      5. metadata-eval N/A

         \[\leadsto \frac{a}{k \cdot \left(\color{blue}{10 \cdot 1} + k\right) + 1} \]

      6. lft-mult-inverse N/A

         \[\leadsto \frac{a}{k \cdot \left(10 \cdot \color{blue}{\left(\frac{1}{k} \cdot k\right)} + k\right) + 1} \]

      7. associate-*l* N/A

         \[\leadsto \frac{a}{k \cdot \left(\color{blue}{\left(10 \cdot \frac{1}{k}\right) \cdot k} + k\right) + 1} \]

      8. *-lft-identity N/A

         \[\leadsto \frac{a}{k \cdot \left(\left(10 \cdot \frac{1}{k}\right) \cdot k + \color{blue}{1 \cdot k}\right) + 1} \]

      9. distribute-rgt-in N/A

         \[\leadsto \frac{a}{k \cdot \color{blue}{\left(k \cdot \left(10 \cdot \frac{1}{k} + 1\right)\right)} + 1} \]

      10. +-commutative N/A

          \[\leadsto \frac{a}{k \cdot \left(k \cdot \color{blue}{\left(1 + 10 \cdot \frac{1}{k}\right)}\right) + 1} \]

      11. *-commutative N/A

          \[\leadsto \frac{a}{k \cdot \color{blue}{\left(\left(1 + 10 \cdot \frac{1}{k}\right) \cdot k\right)} + 1} \]

      12. lower-fma.f64 N/A

          \[\leadsto \frac{a}{\color{blue}{\mathsf{fma}\left(k, \left(1 + 10 \cdot \frac{1}{k}\right) \cdot k, 1\right)}} \]

      13. *-commutative N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{k \cdot \left(1 + 10 \cdot \frac{1}{k}\right)}, 1\right)} \]

      14. +-commutative N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, k \cdot \color{blue}{\left(10 \cdot \frac{1}{k} + 1\right)}, 1\right)} \]

      15. distribute-rgt-in N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{\left(10 \cdot \frac{1}{k}\right) \cdot k + 1 \cdot k}, 1\right)} \]

      16. associate-*l* N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{10 \cdot \left(\frac{1}{k} \cdot k\right)} + 1 \cdot k, 1\right)} \]

      17. lft-mult-inverse N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, 10 \cdot \color{blue}{1} + 1 \cdot k, 1\right)} \]

      18. metadata-eval N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{10} + 1 \cdot k, 1\right)} \]

      19. *-lft-identity N/A

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, 10 + \color{blue}{k}, 1\right)} \]

      20. lower-+.f64 2.7

          \[\leadsto \frac{a}{\mathsf{fma}\left(k, \color{blue}{10 + k}, 1\right)} \]

   5. Applied rewrites 2.7%

      \[\leadsto \color{blue}{\frac{a}{\mathsf{fma}\left(k, 10 + k, 1\right)}} \]

   6. Taylor expanded in k around 0

      \[\leadsto a + \color{blue}{-10 \cdot \left(a \cdot k\right)} \]
   7. Step-by-step derivation

   8. Applied rewrites 8.7%

      \[\leadsto \mathsf{fma}\left(a, \color{blue}{k \cdot -10}, a\right) \]

   9. Taylor expanded in k around inf

      \[\leadsto -10 \cdot \left(a \cdot \color{blue}{k}\right) \]
   10. Step-by-step derivation

   11. Applied rewrites 18.5%

       \[\leadsto a \cdot \left(k \cdot \color{blue}{-10}\right) \]
3. Recombined 2 regimes into one program.

4. Final simplification 24.5%

   \[\leadsto \begin{array}{l} \mathbf{if}\;m \leq 220000:\\ \;\;\;\;a \cdot 1\\ \mathbf{else}:\\ \;\;\;\;a \cdot \left(k \cdot -10\right)\\ \end{array} \]
5. Add Preprocessing

Alternative 15: 19.8% accurate, 22.3× speedup

Math

\[\begin{array}{l} a\_m = \left|a\right| \\ a\_s = \mathsf{copysign}\left(1, a\right) \\ a\_s \cdot \left(a\_m \cdot 1\right) \end{array} \]

FPCore

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 1.0)))

C

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 * 1.0);
}

Fortran

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 * 1.0d0)
end function

Java

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 * 1.0);
}

Python

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 * 1.0)

Julia

a\_m = abs(a)
a\_s = copysign(1.0, a)
function code(a_s, a_m, k, m)
	return Float64(a_s * Float64(a_m * 1.0))
end

MATLAB

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 * 1.0);
end

Wolfram

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 * N[(a$95$m * 1.0), $MachinePrecision]), $MachinePrecision]

Derivation

1. Initial program 88.5%

   \[\frac{a \cdot {k}^{m}}{\left(1 + 10 \cdot k\right) + k \cdot k} \]
2. Add Preprocessing
3. Step-by-step derivation

   1. lift-/.f64 N/A

      \[\leadsto \color{blue}{\frac{a \cdot {k}^{m}}{\left(1 + 10 \cdot k\right) + k \cdot k}} \]

   2. lift-*.f64 N/A

      \[\leadsto \frac{\color{blue}{a \cdot {k}^{m}}}{\left(1 + 10 \cdot k\right) + k \cdot k} \]

   3. associate-/l* N/A

      \[\leadsto \color{blue}{a \cdot \frac{{k}^{m}}{\left(1 + 10 \cdot k\right) + k \cdot k}} \]

   4. *-commutative N/A

      \[\leadsto \color{blue}{\frac{{k}^{m}}{\left(1 + 10 \cdot k\right) + k \cdot k} \cdot a} \]

   5. lower-*.f64 N/A

      \[\leadsto \color{blue}{\frac{{k}^{m}}{\left(1 + 10 \cdot k\right) + k \cdot k} \cdot a} \]

   6. lower-/.f64 88.5

      \[\leadsto \color{blue}{\frac{{k}^{m}}{\left(1 + 10 \cdot k\right) + k \cdot k}} \cdot a \]

   7. lift-+.f64 N/A

      \[\leadsto \frac{{k}^{m}}{\color{blue}{\left(1 + 10 \cdot k\right) + k \cdot k}} \cdot a \]

   8. lift-+.f64 N/A

      \[\leadsto \frac{{k}^{m}}{\color{blue}{\left(1 + 10 \cdot k\right)} + k \cdot k} \cdot a \]

   9. associate-+l+ N/A

      \[\leadsto \frac{{k}^{m}}{\color{blue}{1 + \left(10 \cdot k + k \cdot k\right)}} \cdot a \]

   10. +-commutative N/A

       \[\leadsto \frac{{k}^{m}}{\color{blue}{\left(10 \cdot k + k \cdot k\right) + 1}} \cdot a \]

   11. lift-*.f64 N/A

       \[\leadsto \frac{{k}^{m}}{\left(\color{blue}{10 \cdot k} + k \cdot k\right) + 1} \cdot a \]

   12. lift-*.f64 N/A

       \[\leadsto \frac{{k}^{m}}{\left(10 \cdot k + \color{blue}{k \cdot k}\right) + 1} \cdot a \]

   13. distribute-rgt-out N/A

       \[\leadsto \frac{{k}^{m}}{\color{blue}{k \cdot \left(10 + k\right)} + 1} \cdot a \]

   14. lower-fma.f64 N/A

       \[\leadsto \frac{{k}^{m}}{\color{blue}{\mathsf{fma}\left(k, 10 + k, 1\right)}} \cdot a \]

   15. lower-+.f64 88.9

       \[\leadsto \frac{{k}^{m}}{\mathsf{fma}\left(k, \color{blue}{10 + k}, 1\right)} \cdot a \]

4. Applied rewrites 88.9%

   \[\leadsto \color{blue}{\frac{{k}^{m}}{\mathsf{fma}\left(k, 10 + k, 1\right)} \cdot a} \]

5. Taylor expanded in k around 0

   \[\leadsto \color{blue}{{k}^{m}} \cdot a \]
6. Step-by-step derivation

   1. lower-pow.f64 79.6

      \[\leadsto \color{blue}{{k}^{m}} \cdot a \]

7. Applied rewrites 79.6%

   \[\leadsto \color{blue}{{k}^{m}} \cdot a \]

8. Taylor expanded in m around 0

   \[\leadsto 1 \cdot a \]
9. Step-by-step derivation

10. Applied rewrites 20.0%

    \[\leadsto 1 \cdot a \]

11. Final simplification 20.0%

    \[\leadsto a \cdot 1 \]
12. Add Preprocessing

Reproduce

herbie shell --seed 2024214 
(FPCore (a k m)
  :name "Falkner and Boettcher, Appendix A"
  :precision binary64
  (/ (* a (pow k m)) (+ (+ 1.0 (* 10.0 k)) (* k k))))