exp-w (used to crash)

Percentage Accurate: 99.4% → 99.4%

Time: 10.7s

Alternatives: 5

Speedup: 1.0×

• could not determine a ground truth

  1. w = 2.00000000000000007e154
  2. l = -2

Specification

Math

\[\begin{array}{l} \\ e^{-w} \cdot {\ell}^{\left(e^{w}\right)} \end{array} \]

FPCore

(FPCore (w l) :precision binary64 (* (exp (- w)) (pow l (exp w))))

C

double code(double w, double l) {
	return exp(-w) * pow(l, exp(w));
}

Fortran

real(8) function code(w, l)
    real(8), intent (in) :: w
    real(8), intent (in) :: l
    code = exp(-w) * (l ** exp(w))
end function

Java

public static double code(double w, double l) {
	return Math.exp(-w) * Math.pow(l, Math.exp(w));
}

Python

def code(w, l):
	return math.exp(-w) * math.pow(l, math.exp(w))

Julia

function code(w, l)
	return Float64(exp(Float64(-w)) * (l ^ exp(w)))
end

MATLAB

function tmp = code(w, l)
	tmp = exp(-w) * (l ^ exp(w));
end

Wolfram

code[w_, l_] := N[(N[Exp[(-w)], $MachinePrecision] * N[Power[l, N[Exp[w], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]

Initial Program: 99.4% accurate, 1.0× speedup

Math

\[\begin{array}{l} \\ e^{-w} \cdot {\ell}^{\left(e^{w}\right)} \end{array} \]

FPCore

(FPCore (w l) :precision binary64 (* (exp (- w)) (pow l (exp w))))

C

double code(double w, double l) {
	return exp(-w) * pow(l, exp(w));
}

Fortran

real(8) function code(w, l)
    real(8), intent (in) :: w
    real(8), intent (in) :: l
    code = exp(-w) * (l ** exp(w))
end function

Java

public static double code(double w, double l) {
	return Math.exp(-w) * Math.pow(l, Math.exp(w));
}

Python

def code(w, l):
	return math.exp(-w) * math.pow(l, math.exp(w))

Julia

function code(w, l)
	return Float64(exp(Float64(-w)) * (l ^ exp(w)))
end

MATLAB

function tmp = code(w, l)
	tmp = exp(-w) * (l ^ exp(w));
end

Wolfram

code[w_, l_] := N[(N[Exp[(-w)], $MachinePrecision] * N[Power[l, N[Exp[w], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]

Alternative 1: 99.4% accurate, 1.0× speedup

Math

\[\begin{array}{l} \\ e^{-w} \cdot {\ell}^{\left(e^{w}\right)} \end{array} \]

FPCore

(FPCore (w l) :precision binary64 (* (exp (- w)) (pow l (exp w))))

C

double code(double w, double l) {
	return exp(-w) * pow(l, exp(w));
}

Fortran

real(8) function code(w, l)
    real(8), intent (in) :: w
    real(8), intent (in) :: l
    code = exp(-w) * (l ** exp(w))
end function

Java

public static double code(double w, double l) {
	return Math.exp(-w) * Math.pow(l, Math.exp(w));
}

Python

def code(w, l):
	return math.exp(-w) * math.pow(l, math.exp(w))

Julia

function code(w, l)
	return Float64(exp(Float64(-w)) * (l ^ exp(w)))
end

MATLAB

function tmp = code(w, l)
	tmp = exp(-w) * (l ^ exp(w));
end

Wolfram

code[w_, l_] := N[(N[Exp[(-w)], $MachinePrecision] * N[Power[l, N[Exp[w], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]

Derivation

1. Initial program 99.4%

   \[e^{-w} \cdot {\ell}^{\left(e^{w}\right)} \]
2. Add Preprocessing

3. Final simplification 99.4%

   \[\leadsto e^{-w} \cdot {\ell}^{\left(e^{w}\right)} \]
4. Add Preprocessing

Alternative 2: 99.4% accurate, 1.0× speedup

Math

\[\begin{array}{l} \\ \frac{{\ell}^{\left(e^{w}\right)}}{e^{w}} \end{array} \]

FPCore

(FPCore (w l) :precision binary64 (/ (pow l (exp w)) (exp w)))

C

double code(double w, double l) {
	return pow(l, exp(w)) / exp(w);
}

Fortran

real(8) function code(w, l)
    real(8), intent (in) :: w
    real(8), intent (in) :: l
    code = (l ** exp(w)) / exp(w)
end function

Java

public static double code(double w, double l) {
	return Math.pow(l, Math.exp(w)) / Math.exp(w);
}

Python

def code(w, l):
	return math.pow(l, math.exp(w)) / math.exp(w)

Julia

function code(w, l)
	return Float64((l ^ exp(w)) / exp(w))
end

MATLAB

function tmp = code(w, l)
	tmp = (l ^ exp(w)) / exp(w);
end

Wolfram

code[w_, l_] := N[(N[Power[l, N[Exp[w], $MachinePrecision]], $MachinePrecision] / N[Exp[w], $MachinePrecision]), $MachinePrecision]

Derivation

1. Initial program 99.4%

   \[e^{-w} \cdot {\ell}^{\left(e^{w}\right)} \]
2. Step-by-step derivation

   1. exp-neg 99.4%

      \[\leadsto \color{blue}{\frac{1}{e^{w}}} \cdot {\ell}^{\left(e^{w}\right)} \]

   2. associate-*l/ 99.4%

      \[\leadsto \color{blue}{\frac{1 \cdot {\ell}^{\left(e^{w}\right)}}{e^{w}}} \]

   3. *-lft-identity 99.4%

      \[\leadsto \frac{\color{blue}{{\ell}^{\left(e^{w}\right)}}}{e^{w}} \]

3. Simplified 99.4%

   \[\leadsto \color{blue}{\frac{{\ell}^{\left(e^{w}\right)}}{e^{w}}} \]
4. Add Preprocessing

5. Final simplification 99.4%

   \[\leadsto \frac{{\ell}^{\left(e^{w}\right)}}{e^{w}} \]
6. Add Preprocessing

Alternative 3: 97.4% accurate, 3.0× speedup

Math

\[\begin{array}{l} \\ \frac{\ell}{e^{w}} \end{array} \]

FPCore

(FPCore (w l) :precision binary64 (/ l (exp w)))

C

double code(double w, double l) {
	return l / exp(w);
}

Fortran

real(8) function code(w, l)
    real(8), intent (in) :: w
    real(8), intent (in) :: l
    code = l / exp(w)
end function

Java

public static double code(double w, double l) {
	return l / Math.exp(w);
}

Python

def code(w, l):
	return l / math.exp(w)

Julia

function code(w, l)
	return Float64(l / exp(w))
end

MATLAB

function tmp = code(w, l)
	tmp = l / exp(w);
end

Wolfram

code[w_, l_] := N[(l / N[Exp[w], $MachinePrecision]), $MachinePrecision]

Derivation

1. Initial program 99.4%

   \[e^{-w} \cdot {\ell}^{\left(e^{w}\right)} \]
2. Step-by-step derivation

   1. exp-neg 99.4%

      \[\leadsto \color{blue}{\frac{1}{e^{w}}} \cdot {\ell}^{\left(e^{w}\right)} \]

   2. associate-*l/ 99.4%

      \[\leadsto \color{blue}{\frac{1 \cdot {\ell}^{\left(e^{w}\right)}}{e^{w}}} \]

   3. *-lft-identity 99.4%

      \[\leadsto \frac{\color{blue}{{\ell}^{\left(e^{w}\right)}}}{e^{w}} \]

3. Simplified 99.4%

   \[\leadsto \color{blue}{\frac{{\ell}^{\left(e^{w}\right)}}{e^{w}}} \]
4. Add Preprocessing
5. Step-by-step derivation

   1. add-sqr-sqrt 41.7%

      \[\leadsto \frac{{\ell}^{\left(e^{\color{blue}{\sqrt{w} \cdot \sqrt{w}}}\right)}}{e^{w}} \]

   2. sqrt-unprod 83.4%

      \[\leadsto \frac{{\ell}^{\left(e^{\color{blue}{\sqrt{w \cdot w}}}\right)}}{e^{w}} \]

   3. sqr-neg 83.4%

      \[\leadsto \frac{{\ell}^{\left(e^{\sqrt{\color{blue}{\left(-w\right) \cdot \left(-w\right)}}}\right)}}{e^{w}} \]

   4. sqrt-unprod 41.7%

      \[\leadsto \frac{{\ell}^{\left(e^{\color{blue}{\sqrt{-w} \cdot \sqrt{-w}}}\right)}}{e^{w}} \]

   5. add-sqr-sqrt 83.3%

      \[\leadsto \frac{{\ell}^{\left(e^{\color{blue}{-w}}\right)}}{e^{w}} \]

   6. add-sqr-sqrt 83.3%

      \[\leadsto \frac{{\ell}^{\color{blue}{\left(\sqrt{e^{-w}} \cdot \sqrt{e^{-w}}\right)}}}{e^{w}} \]

   7. sqrt-unprod 83.3%

      \[\leadsto \frac{{\ell}^{\color{blue}{\left(\sqrt{e^{-w} \cdot e^{-w}}\right)}}}{e^{w}} \]

   8. add-sqr-sqrt 41.7%

      \[\leadsto \frac{{\ell}^{\left(\sqrt{e^{\color{blue}{\sqrt{-w} \cdot \sqrt{-w}}} \cdot e^{-w}}\right)}}{e^{w}} \]

   9. sqrt-unprod 70.8%

      \[\leadsto \frac{{\ell}^{\left(\sqrt{e^{\color{blue}{\sqrt{\left(-w\right) \cdot \left(-w\right)}}} \cdot e^{-w}}\right)}}{e^{w}} \]

   10. sqr-neg 70.8%

       \[\leadsto \frac{{\ell}^{\left(\sqrt{e^{\sqrt{\color{blue}{w \cdot w}}} \cdot e^{-w}}\right)}}{e^{w}} \]

   11. sqrt-unprod 29.1%

       \[\leadsto \frac{{\ell}^{\left(\sqrt{e^{\color{blue}{\sqrt{w} \cdot \sqrt{w}}} \cdot e^{-w}}\right)}}{e^{w}} \]

   12. add-sqr-sqrt 54.8%

       \[\leadsto \frac{{\ell}^{\left(\sqrt{e^{\color{blue}{w}} \cdot e^{-w}}\right)}}{e^{w}} \]

   13. pow1 54.8%

       \[\leadsto \frac{{\ell}^{\left(\sqrt{\color{blue}{{\left(e^{w}\right)}^{1}} \cdot e^{-w}}\right)}}{e^{w}} \]

   14. exp-neg 54.8%

       \[\leadsto \frac{{\ell}^{\left(\sqrt{{\left(e^{w}\right)}^{1} \cdot \color{blue}{\frac{1}{e^{w}}}}\right)}}{e^{w}} \]

   15. inv-pow 54.8%

       \[\leadsto \frac{{\ell}^{\left(\sqrt{{\left(e^{w}\right)}^{1} \cdot \color{blue}{{\left(e^{w}\right)}^{-1}}}\right)}}{e^{w}} \]

   16. pow-prod-up 98.6%

       \[\leadsto \frac{{\ell}^{\left(\sqrt{\color{blue}{{\left(e^{w}\right)}^{\left(1 + -1\right)}}}\right)}}{e^{w}} \]

   17. metadata-eval 98.6%

       \[\leadsto \frac{{\ell}^{\left(\sqrt{{\left(e^{w}\right)}^{\color{blue}{0}}}\right)}}{e^{w}} \]

   18. metadata-eval 98.6%

       \[\leadsto \frac{{\ell}^{\left(\sqrt{\color{blue}{1}}\right)}}{e^{w}} \]

   19. metadata-eval 98.6%

       \[\leadsto \frac{{\ell}^{\color{blue}{1}}}{e^{w}} \]

   20. pow1 98.6%

       \[\leadsto \frac{\color{blue}{\ell}}{e^{w}} \]

   21. *-un-lft-identity 98.6%

       \[\leadsto \frac{\color{blue}{1 \cdot \ell}}{e^{w}} \]

   22. *-commutative 98.6%

       \[\leadsto \frac{\color{blue}{\ell \cdot 1}}{e^{w}} \]

6. Applied egg-rr 98.6%

   \[\leadsto \frac{\color{blue}{\ell \cdot 1}}{e^{w}} \]

7. Taylor expanded in l around 0 98.6%

   \[\leadsto \color{blue}{\frac{\ell}{e^{w}}} \]

8. Final simplification 98.6%

   \[\leadsto \frac{\ell}{e^{w}} \]
9. Add Preprocessing

Alternative 4: 63.4% accurate, 61.0× speedup

Math

\[\begin{array}{l} \\ \ell - w \cdot \ell \end{array} \]

FPCore

(FPCore (w l) :precision binary64 (- l (* w l)))

C

double code(double w, double l) {
	return l - (w * l);
}

Fortran

real(8) function code(w, l)
    real(8), intent (in) :: w
    real(8), intent (in) :: l
    code = l - (w * l)
end function

Java

public static double code(double w, double l) {
	return l - (w * l);
}

Python

def code(w, l):
	return l - (w * l)

Julia

function code(w, l)
	return Float64(l - Float64(w * l))
end

MATLAB

function tmp = code(w, l)
	tmp = l - (w * l);
end

Wolfram

code[w_, l_] := N[(l - N[(w * l), $MachinePrecision]), $MachinePrecision]

Derivation

1. Initial program 99.4%

   \[e^{-w} \cdot {\ell}^{\left(e^{w}\right)} \]
2. Step-by-step derivation

   1. exp-neg 99.4%

      \[\leadsto \color{blue}{\frac{1}{e^{w}}} \cdot {\ell}^{\left(e^{w}\right)} \]

   2. associate-*l/ 99.4%

      \[\leadsto \color{blue}{\frac{1 \cdot {\ell}^{\left(e^{w}\right)}}{e^{w}}} \]

   3. *-lft-identity 99.4%

      \[\leadsto \frac{\color{blue}{{\ell}^{\left(e^{w}\right)}}}{e^{w}} \]

3. Simplified 99.4%

   \[\leadsto \color{blue}{\frac{{\ell}^{\left(e^{w}\right)}}{e^{w}}} \]
4. Add Preprocessing
5. Step-by-step derivation

   1. add-sqr-sqrt 41.7%

      \[\leadsto \frac{{\ell}^{\left(e^{\color{blue}{\sqrt{w} \cdot \sqrt{w}}}\right)}}{e^{w}} \]

   2. sqrt-unprod 83.4%

      \[\leadsto \frac{{\ell}^{\left(e^{\color{blue}{\sqrt{w \cdot w}}}\right)}}{e^{w}} \]

   3. sqr-neg 83.4%

      \[\leadsto \frac{{\ell}^{\left(e^{\sqrt{\color{blue}{\left(-w\right) \cdot \left(-w\right)}}}\right)}}{e^{w}} \]

   4. sqrt-unprod 41.7%

      \[\leadsto \frac{{\ell}^{\left(e^{\color{blue}{\sqrt{-w} \cdot \sqrt{-w}}}\right)}}{e^{w}} \]

   5. add-sqr-sqrt 83.3%

      \[\leadsto \frac{{\ell}^{\left(e^{\color{blue}{-w}}\right)}}{e^{w}} \]

   6. add-sqr-sqrt 83.3%

      \[\leadsto \frac{{\ell}^{\color{blue}{\left(\sqrt{e^{-w}} \cdot \sqrt{e^{-w}}\right)}}}{e^{w}} \]

   7. sqrt-unprod 83.3%

      \[\leadsto \frac{{\ell}^{\color{blue}{\left(\sqrt{e^{-w} \cdot e^{-w}}\right)}}}{e^{w}} \]

   8. add-sqr-sqrt 41.7%

      \[\leadsto \frac{{\ell}^{\left(\sqrt{e^{\color{blue}{\sqrt{-w} \cdot \sqrt{-w}}} \cdot e^{-w}}\right)}}{e^{w}} \]

   9. sqrt-unprod 70.8%

      \[\leadsto \frac{{\ell}^{\left(\sqrt{e^{\color{blue}{\sqrt{\left(-w\right) \cdot \left(-w\right)}}} \cdot e^{-w}}\right)}}{e^{w}} \]

   10. sqr-neg 70.8%

       \[\leadsto \frac{{\ell}^{\left(\sqrt{e^{\sqrt{\color{blue}{w \cdot w}}} \cdot e^{-w}}\right)}}{e^{w}} \]

   11. sqrt-unprod 29.1%

       \[\leadsto \frac{{\ell}^{\left(\sqrt{e^{\color{blue}{\sqrt{w} \cdot \sqrt{w}}} \cdot e^{-w}}\right)}}{e^{w}} \]

   12. add-sqr-sqrt 54.8%

       \[\leadsto \frac{{\ell}^{\left(\sqrt{e^{\color{blue}{w}} \cdot e^{-w}}\right)}}{e^{w}} \]

   13. pow1 54.8%

       \[\leadsto \frac{{\ell}^{\left(\sqrt{\color{blue}{{\left(e^{w}\right)}^{1}} \cdot e^{-w}}\right)}}{e^{w}} \]

   14. exp-neg 54.8%

       \[\leadsto \frac{{\ell}^{\left(\sqrt{{\left(e^{w}\right)}^{1} \cdot \color{blue}{\frac{1}{e^{w}}}}\right)}}{e^{w}} \]

   15. inv-pow 54.8%

       \[\leadsto \frac{{\ell}^{\left(\sqrt{{\left(e^{w}\right)}^{1} \cdot \color{blue}{{\left(e^{w}\right)}^{-1}}}\right)}}{e^{w}} \]

   16. pow-prod-up 98.6%

       \[\leadsto \frac{{\ell}^{\left(\sqrt{\color{blue}{{\left(e^{w}\right)}^{\left(1 + -1\right)}}}\right)}}{e^{w}} \]

   17. metadata-eval 98.6%

       \[\leadsto \frac{{\ell}^{\left(\sqrt{{\left(e^{w}\right)}^{\color{blue}{0}}}\right)}}{e^{w}} \]

   18. metadata-eval 98.6%

       \[\leadsto \frac{{\ell}^{\left(\sqrt{\color{blue}{1}}\right)}}{e^{w}} \]

   19. metadata-eval 98.6%

       \[\leadsto \frac{{\ell}^{\color{blue}{1}}}{e^{w}} \]

   20. pow1 98.6%

       \[\leadsto \frac{\color{blue}{\ell}}{e^{w}} \]

   21. *-un-lft-identity 98.6%

       \[\leadsto \frac{\color{blue}{1 \cdot \ell}}{e^{w}} \]

   22. *-commutative 98.6%

       \[\leadsto \frac{\color{blue}{\ell \cdot 1}}{e^{w}} \]

6. Applied egg-rr 98.6%

   \[\leadsto \frac{\color{blue}{\ell \cdot 1}}{e^{w}} \]

7. Taylor expanded in w around 0 64.5%

   \[\leadsto \color{blue}{\ell + -1 \cdot \left(\ell \cdot w\right)} \]
8. Step-by-step derivation

   1. mul-1-neg 64.5%

      \[\leadsto \ell + \color{blue}{\left(-\ell \cdot w\right)} \]

   2. unsub-neg 64.5%

      \[\leadsto \color{blue}{\ell - \ell \cdot w} \]

9. Simplified 64.5%

   \[\leadsto \color{blue}{\ell - \ell \cdot w} \]

10. Final simplification 64.5%

    \[\leadsto \ell - w \cdot \ell \]
11. Add Preprocessing

Alternative 5: 56.0% accurate, 305.0× speedup

Math

\[\begin{array}{l} \\ \ell \end{array} \]

FPCore

(FPCore (w l) :precision binary64 l)

C

double code(double w, double l) {
	return l;
}

Fortran

real(8) function code(w, l)
    real(8), intent (in) :: w
    real(8), intent (in) :: l
    code = l
end function

Java

public static double code(double w, double l) {
	return l;
}

Python

def code(w, l):
	return l

Julia

function code(w, l)
	return l
end

MATLAB

function tmp = code(w, l)
	tmp = l;
end

Wolfram

code[w_, l_] := l

Derivation

1. Initial program 99.4%

   \[e^{-w} \cdot {\ell}^{\left(e^{w}\right)} \]
2. Add Preprocessing

3. Taylor expanded in w around 0 56.9%

   \[\leadsto \color{blue}{\ell} \]

4. Final simplification 56.9%

   \[\leadsto \ell \]
5. Add Preprocessing

Reproduce

herbie shell --seed 2024031 
(FPCore (w l)
  :name "exp-w (used to crash)"
  :precision binary64
  (* (exp (- w)) (pow l (exp w))))