Rust f64::asinh

Percentage Accurate: 30.0% → 99.1%

Time: 14.1s

Alternatives: 11

Speedup: 4.0×

Specification

Math

\[\begin{array}{l} \\ \sinh^{-1} x \end{array} \]

FPCore

(FPCore (x) :precision binary64 (asinh x))

C

double code(double x) {
	return asinh(x);
}

Python

def code(x):
	return math.asinh(x)

Julia

function code(x)
	return asinh(x)
end

MATLAB

function tmp = code(x)
	tmp = asinh(x);
end

Wolfram

code[x_] := N[ArcSinh[x], $MachinePrecision]

Initial Program: 30.0% accurate, 1.0× speedup

Math

\[\begin{array}{l} \\ \mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{x \cdot x + 1}\right), x\right) \end{array} \]

FPCore

(FPCore (x)
 :precision binary64
 (copysign (log (+ (fabs x) (sqrt (+ (* x x) 1.0)))) x))

C

double code(double x) {
	return copysign(log((fabs(x) + sqrt(((x * x) + 1.0)))), x);
}

Java

public static double code(double x) {
	return Math.copySign(Math.log((Math.abs(x) + Math.sqrt(((x * x) + 1.0)))), x);
}

Python

def code(x):
	return math.copysign(math.log((math.fabs(x) + math.sqrt(((x * x) + 1.0)))), x)

Julia

function code(x)
	return copysign(log(Float64(abs(x) + sqrt(Float64(Float64(x * x) + 1.0)))), x)
end

MATLAB

function tmp = code(x)
	tmp = sign(x) * abs(log((abs(x) + sqrt(((x * x) + 1.0)))));
end

Wolfram

code[x_] := N[With[{TMP1 = Abs[N[Log[N[(N[Abs[x], $MachinePrecision] + N[Sqrt[N[(N[(x * x), $MachinePrecision] + 1.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], TMP2 = Sign[x]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]

Alternative 1: 99.1% accurate, 0.2× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} t_0 := \mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{x \cdot x + 1}\right), x\right)\\ t_1 := \left|x\right| + 1\\ t_2 := {t_1}^{2}\\ \mathbf{if}\;t_0 \leq -10:\\ \;\;\;\;\mathsf{copysign}\left(\log \left(\frac{1}{\mathsf{hypot}\left(1, x\right) - x}\right), x\right)\\ \mathbf{elif}\;t_0 \leq 0:\\ \;\;\;\;\mathsf{copysign}\left(\mathsf{fma}\left(-0.041666666666666664, \left(\frac{3}{t_1} + \frac{3}{t_2}\right) \cdot {x}^{4}, \mathsf{fma}\left(0.001388888888888889, \left(\frac{45}{t_1} + \left(\frac{30}{{t_1}^{3}} + \frac{45}{t_2}\right)\right) \cdot {x}^{6}, \mathsf{fma}\left(0.5, \frac{x \cdot x}{t_1}, \mathsf{log1p}\left(\left|x\right|\right)\right)\right)\right), x\right)\\ \mathbf{else}:\\ \;\;\;\;\mathsf{copysign}\left(\log \left(x + \mathsf{hypot}\left(1, x\right)\right), x\right)\\ \end{array} \end{array} \]

FPCore

(FPCore (x)
 :precision binary64
 (let* ((t_0 (copysign (log (+ (fabs x) (sqrt (+ (* x x) 1.0)))) x))
        (t_1 (+ (fabs x) 1.0))
        (t_2 (pow t_1 2.0)))
   (if (<= t_0 -10.0)
     (copysign (log (/ 1.0 (- (hypot 1.0 x) x))) x)
     (if (<= t_0 0.0)
       (copysign
        (fma
         -0.041666666666666664
         (* (+ (/ 3.0 t_1) (/ 3.0 t_2)) (pow x 4.0))
         (fma
          0.001388888888888889
          (*
           (+ (/ 45.0 t_1) (+ (/ 30.0 (pow t_1 3.0)) (/ 45.0 t_2)))
           (pow x 6.0))
          (fma 0.5 (/ (* x x) t_1) (log1p (fabs x)))))
        x)
       (copysign (log (+ x (hypot 1.0 x))) x)))))

C

double code(double x) {
	double t_0 = copysign(log((fabs(x) + sqrt(((x * x) + 1.0)))), x);
	double t_1 = fabs(x) + 1.0;
	double t_2 = pow(t_1, 2.0);
	double tmp;
	if (t_0 <= -10.0) {
		tmp = copysign(log((1.0 / (hypot(1.0, x) - x))), x);
	} else if (t_0 <= 0.0) {
		tmp = copysign(fma(-0.041666666666666664, (((3.0 / t_1) + (3.0 / t_2)) * pow(x, 4.0)), fma(0.001388888888888889, (((45.0 / t_1) + ((30.0 / pow(t_1, 3.0)) + (45.0 / t_2))) * pow(x, 6.0)), fma(0.5, ((x * x) / t_1), log1p(fabs(x))))), x);
	} else {
		tmp = copysign(log((x + hypot(1.0, x))), x);
	}
	return tmp;
}

Julia

function code(x)
	t_0 = copysign(log(Float64(abs(x) + sqrt(Float64(Float64(x * x) + 1.0)))), x)
	t_1 = Float64(abs(x) + 1.0)
	t_2 = t_1 ^ 2.0
	tmp = 0.0
	if (t_0 <= -10.0)
		tmp = copysign(log(Float64(1.0 / Float64(hypot(1.0, x) - x))), x);
	elseif (t_0 <= 0.0)
		tmp = copysign(fma(-0.041666666666666664, Float64(Float64(Float64(3.0 / t_1) + Float64(3.0 / t_2)) * (x ^ 4.0)), fma(0.001388888888888889, Float64(Float64(Float64(45.0 / t_1) + Float64(Float64(30.0 / (t_1 ^ 3.0)) + Float64(45.0 / t_2))) * (x ^ 6.0)), fma(0.5, Float64(Float64(x * x) / t_1), log1p(abs(x))))), x);
	else
		tmp = copysign(log(Float64(x + hypot(1.0, x))), x);
	end
	return tmp
end

Wolfram

code[x_] := Block[{t$95$0 = N[With[{TMP1 = Abs[N[Log[N[(N[Abs[x], $MachinePrecision] + N[Sqrt[N[(N[(x * x), $MachinePrecision] + 1.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], TMP2 = Sign[x]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]}, Block[{t$95$1 = N[(N[Abs[x], $MachinePrecision] + 1.0), $MachinePrecision]}, Block[{t$95$2 = N[Power[t$95$1, 2.0], $MachinePrecision]}, If[LessEqual[t$95$0, -10.0], N[With[{TMP1 = Abs[N[Log[N[(1.0 / N[(N[Sqrt[1.0 ^ 2 + x ^ 2], $MachinePrecision] - x), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], TMP2 = Sign[x]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision], If[LessEqual[t$95$0, 0.0], N[With[{TMP1 = Abs[N[(-0.041666666666666664 * N[(N[(N[(3.0 / t$95$1), $MachinePrecision] + N[(3.0 / t$95$2), $MachinePrecision]), $MachinePrecision] * N[Power[x, 4.0], $MachinePrecision]), $MachinePrecision] + N[(0.001388888888888889 * N[(N[(N[(45.0 / t$95$1), $MachinePrecision] + N[(N[(30.0 / N[Power[t$95$1, 3.0], $MachinePrecision]), $MachinePrecision] + N[(45.0 / t$95$2), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[Power[x, 6.0], $MachinePrecision]), $MachinePrecision] + N[(0.5 * N[(N[(x * x), $MachinePrecision] / t$95$1), $MachinePrecision] + N[Log[1 + N[Abs[x], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], TMP2 = Sign[x]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision], N[With[{TMP1 = Abs[N[Log[N[(x + N[Sqrt[1.0 ^ 2 + x ^ 2], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], TMP2 = Sign[x]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]]]]]]

Derivation

1. Split input into 3 regimes

2. if (copysign.f64 (log.f64 (+.f64 (fabs.f64 x) (sqrt.f64 (+.f64 (*.f64 x x) 1)))) x) < -10

   1. Initial program 50.8%

      \[\mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{x \cdot x + 1}\right), x\right) \]
   2. Step-by-step derivation

      1. +-commutative 50.8%

         \[\leadsto \mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{\color{blue}{1 + x \cdot x}}\right), x\right) \]

      2. hypot-1-def 100.0%

         \[\leadsto \mathsf{copysign}\left(\log \left(\left|x\right| + \color{blue}{\mathsf{hypot}\left(1, x\right)}\right), x\right) \]

   3. Simplified 100.0%

      \[\leadsto \color{blue}{\mathsf{copysign}\left(\log \left(\left|x\right| + \mathsf{hypot}\left(1, x\right)\right), x\right)} \]
   4. Step-by-step derivation

      1. flip-+ 0.9%

         \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(\frac{\left|x\right| \cdot \left|x\right| - \mathsf{hypot}\left(1, x\right) \cdot \mathsf{hypot}\left(1, x\right)}{\left|x\right| - \mathsf{hypot}\left(1, x\right)}\right)}, x\right) \]

      2. div-sub 0.9%

         \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(\frac{\left|x\right| \cdot \left|x\right|}{\left|x\right| - \mathsf{hypot}\left(1, x\right)} - \frac{\mathsf{hypot}\left(1, x\right) \cdot \mathsf{hypot}\left(1, x\right)}{\left|x\right| - \mathsf{hypot}\left(1, x\right)}\right)}, x\right) \]

      3. pow2 0.9%

         \[\leadsto \mathsf{copysign}\left(\log \left(\frac{\color{blue}{{\left(\left|x\right|\right)}^{2}}}{\left|x\right| - \mathsf{hypot}\left(1, x\right)} - \frac{\mathsf{hypot}\left(1, x\right) \cdot \mathsf{hypot}\left(1, x\right)}{\left|x\right| - \mathsf{hypot}\left(1, x\right)}\right), x\right) \]

      4. add-sqr-sqrt 0.0%

         \[\leadsto \mathsf{copysign}\left(\log \left(\frac{{\left(\left|\color{blue}{\sqrt{x} \cdot \sqrt{x}}\right|\right)}^{2}}{\left|x\right| - \mathsf{hypot}\left(1, x\right)} - \frac{\mathsf{hypot}\left(1, x\right) \cdot \mathsf{hypot}\left(1, x\right)}{\left|x\right| - \mathsf{hypot}\left(1, x\right)}\right), x\right) \]

      5. fabs-sqr 0.0%

         \[\leadsto \mathsf{copysign}\left(\log \left(\frac{{\color{blue}{\left(\sqrt{x} \cdot \sqrt{x}\right)}}^{2}}{\left|x\right| - \mathsf{hypot}\left(1, x\right)} - \frac{\mathsf{hypot}\left(1, x\right) \cdot \mathsf{hypot}\left(1, x\right)}{\left|x\right| - \mathsf{hypot}\left(1, x\right)}\right), x\right) \]

      6. add-sqr-sqrt 0.9%

         \[\leadsto \mathsf{copysign}\left(\log \left(\frac{{\color{blue}{x}}^{2}}{\left|x\right| - \mathsf{hypot}\left(1, x\right)} - \frac{\mathsf{hypot}\left(1, x\right) \cdot \mathsf{hypot}\left(1, x\right)}{\left|x\right| - \mathsf{hypot}\left(1, x\right)}\right), x\right) \]

      7. pow2 0.9%

         \[\leadsto \mathsf{copysign}\left(\log \left(\frac{\color{blue}{x \cdot x}}{\left|x\right| - \mathsf{hypot}\left(1, x\right)} - \frac{\mathsf{hypot}\left(1, x\right) \cdot \mathsf{hypot}\left(1, x\right)}{\left|x\right| - \mathsf{hypot}\left(1, x\right)}\right), x\right) \]

      8. add-sqr-sqrt 0.0%

         \[\leadsto \mathsf{copysign}\left(\log \left(\frac{x \cdot x}{\left|\color{blue}{\sqrt{x} \cdot \sqrt{x}}\right| - \mathsf{hypot}\left(1, x\right)} - \frac{\mathsf{hypot}\left(1, x\right) \cdot \mathsf{hypot}\left(1, x\right)}{\left|x\right| - \mathsf{hypot}\left(1, x\right)}\right), x\right) \]

      9. fabs-sqr 0.0%

         \[\leadsto \mathsf{copysign}\left(\log \left(\frac{x \cdot x}{\color{blue}{\sqrt{x} \cdot \sqrt{x}} - \mathsf{hypot}\left(1, x\right)} - \frac{\mathsf{hypot}\left(1, x\right) \cdot \mathsf{hypot}\left(1, x\right)}{\left|x\right| - \mathsf{hypot}\left(1, x\right)}\right), x\right) \]

      10. add-sqr-sqrt 0.3%

          \[\leadsto \mathsf{copysign}\left(\log \left(\frac{x \cdot x}{\color{blue}{x} - \mathsf{hypot}\left(1, x\right)} - \frac{\mathsf{hypot}\left(1, x\right) \cdot \mathsf{hypot}\left(1, x\right)}{\left|x\right| - \mathsf{hypot}\left(1, x\right)}\right), x\right) \]

      11. hypot-udef 0.3%

          \[\leadsto \mathsf{copysign}\left(\log \left(\frac{x \cdot x}{x - \mathsf{hypot}\left(1, x\right)} - \frac{\color{blue}{\sqrt{1 \cdot 1 + x \cdot x}} \cdot \mathsf{hypot}\left(1, x\right)}{\left|x\right| - \mathsf{hypot}\left(1, x\right)}\right), x\right) \]

      12. hypot-udef 0.3%

          \[\leadsto \mathsf{copysign}\left(\log \left(\frac{x \cdot x}{x - \mathsf{hypot}\left(1, x\right)} - \frac{\sqrt{1 \cdot 1 + x \cdot x} \cdot \color{blue}{\sqrt{1 \cdot 1 + x \cdot x}}}{\left|x\right| - \mathsf{hypot}\left(1, x\right)}\right), x\right) \]

      13. add-sqr-sqrt 0.3%

          \[\leadsto \mathsf{copysign}\left(\log \left(\frac{x \cdot x}{x - \mathsf{hypot}\left(1, x\right)} - \frac{\color{blue}{1 \cdot 1 + x \cdot x}}{\left|x\right| - \mathsf{hypot}\left(1, x\right)}\right), x\right) \]

      14. metadata-eval 0.3%

          \[\leadsto \mathsf{copysign}\left(\log \left(\frac{x \cdot x}{x - \mathsf{hypot}\left(1, x\right)} - \frac{\color{blue}{1} + x \cdot x}{\left|x\right| - \mathsf{hypot}\left(1, x\right)}\right), x\right) \]

      15. add-sqr-sqrt 0.0%

          \[\leadsto \mathsf{copysign}\left(\log \left(\frac{x \cdot x}{x - \mathsf{hypot}\left(1, x\right)} - \frac{1 + x \cdot x}{\left|\color{blue}{\sqrt{x} \cdot \sqrt{x}}\right| - \mathsf{hypot}\left(1, x\right)}\right), x\right) \]

      16. fabs-sqr 0.0%

          \[\leadsto \mathsf{copysign}\left(\log \left(\frac{x \cdot x}{x - \mathsf{hypot}\left(1, x\right)} - \frac{1 + x \cdot x}{\color{blue}{\sqrt{x} \cdot \sqrt{x}} - \mathsf{hypot}\left(1, x\right)}\right), x\right) \]

      17. add-sqr-sqrt 2.5%

          \[\leadsto \mathsf{copysign}\left(\log \left(\frac{x \cdot x}{x - \mathsf{hypot}\left(1, x\right)} - \frac{1 + x \cdot x}{\color{blue}{x} - \mathsf{hypot}\left(1, x\right)}\right), x\right) \]

   5. Applied egg-rr 2.5%

      \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(\frac{x \cdot x}{x - \mathsf{hypot}\left(1, x\right)} - \frac{1 + x \cdot x}{x - \mathsf{hypot}\left(1, x\right)}\right)}, x\right) \]
   6. Step-by-step derivation

      1. unpow2 2.5%

         \[\leadsto \mathsf{copysign}\left(\log \left(\frac{\color{blue}{{x}^{2}}}{x - \mathsf{hypot}\left(1, x\right)} - \frac{1 + x \cdot x}{x - \mathsf{hypot}\left(1, x\right)}\right), x\right) \]

      2. div-sub 3.0%

         \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(\frac{{x}^{2} - \left(1 + x \cdot x\right)}{x - \mathsf{hypot}\left(1, x\right)}\right)}, x\right) \]

      3. unpow2 3.0%

         \[\leadsto \mathsf{copysign}\left(\log \left(\frac{\color{blue}{x \cdot x} - \left(1 + x \cdot x\right)}{x - \mathsf{hypot}\left(1, x\right)}\right), x\right) \]

      4. unpow2 3.0%

         \[\leadsto \mathsf{copysign}\left(\log \left(\frac{\color{blue}{{x}^{2}} - \left(1 + x \cdot x\right)}{x - \mathsf{hypot}\left(1, x\right)}\right), x\right) \]

      5. unpow2 3.0%

         \[\leadsto \mathsf{copysign}\left(\log \left(\frac{{x}^{2} - \left(1 + \color{blue}{{x}^{2}}\right)}{x - \mathsf{hypot}\left(1, x\right)}\right), x\right) \]

      6. +-commutative 3.0%

         \[\leadsto \mathsf{copysign}\left(\log \left(\frac{{x}^{2} - \color{blue}{\left({x}^{2} + 1\right)}}{x - \mathsf{hypot}\left(1, x\right)}\right), x\right) \]

      7. associate--r+ 49.2%

         \[\leadsto \mathsf{copysign}\left(\log \left(\frac{\color{blue}{\left({x}^{2} - {x}^{2}\right) - 1}}{x - \mathsf{hypot}\left(1, x\right)}\right), x\right) \]

      8. +-inverses 100.0%

         \[\leadsto \mathsf{copysign}\left(\log \left(\frac{\color{blue}{0} - 1}{x - \mathsf{hypot}\left(1, x\right)}\right), x\right) \]

      9. metadata-eval 100.0%

         \[\leadsto \mathsf{copysign}\left(\log \left(\frac{\color{blue}{-1}}{x - \mathsf{hypot}\left(1, x\right)}\right), x\right) \]

      10. metadata-eval 100.0%

          \[\leadsto \mathsf{copysign}\left(\log \left(\frac{\color{blue}{\frac{1}{-1}}}{x - \mathsf{hypot}\left(1, x\right)}\right), x\right) \]

      11. associate-/r* 100.0%

          \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(\frac{1}{-1 \cdot \left(x - \mathsf{hypot}\left(1, x\right)\right)}\right)}, x\right) \]

      12. neg-mul-1 100.0%

          \[\leadsto \mathsf{copysign}\left(\log \left(\frac{1}{\color{blue}{-\left(x - \mathsf{hypot}\left(1, x\right)\right)}}\right), x\right) \]

      13. sub-neg 100.0%

          \[\leadsto \mathsf{copysign}\left(\log \left(\frac{1}{-\color{blue}{\left(x + \left(-\mathsf{hypot}\left(1, x\right)\right)\right)}}\right), x\right) \]

      14. +-commutative 100.0%

          \[\leadsto \mathsf{copysign}\left(\log \left(\frac{1}{-\color{blue}{\left(\left(-\mathsf{hypot}\left(1, x\right)\right) + x\right)}}\right), x\right) \]

      15. distribute-neg-in 100.0%

          \[\leadsto \mathsf{copysign}\left(\log \left(\frac{1}{\color{blue}{\left(-\left(-\mathsf{hypot}\left(1, x\right)\right)\right) + \left(-x\right)}}\right), x\right) \]

      16. remove-double-neg 100.0%

          \[\leadsto \mathsf{copysign}\left(\log \left(\frac{1}{\color{blue}{\mathsf{hypot}\left(1, x\right)} + \left(-x\right)}\right), x\right) \]

      17. sub-neg 100.0%

          \[\leadsto \mathsf{copysign}\left(\log \left(\frac{1}{\color{blue}{\mathsf{hypot}\left(1, x\right) - x}}\right), x\right) \]

   7. Simplified 100.0%

      \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(\frac{1}{\mathsf{hypot}\left(1, x\right) - x}\right)}, x\right) \]

   if -10 < (copysign.f64 (log.f64 (+.f64 (fabs.f64 x) (sqrt.f64 (+.f64 (*.f64 x x) 1)))) x) < -0.0

   1. Initial program 9.4%

      \[\mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{x \cdot x + 1}\right), x\right) \]
   2. Step-by-step derivation

      1. +-commutative 9.4%

         \[\leadsto \mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{\color{blue}{1 + x \cdot x}}\right), x\right) \]

      2. hypot-1-def 9.5%

         \[\leadsto \mathsf{copysign}\left(\log \left(\left|x\right| + \color{blue}{\mathsf{hypot}\left(1, x\right)}\right), x\right) \]

   3. Simplified 9.5%

      \[\leadsto \color{blue}{\mathsf{copysign}\left(\log \left(\left|x\right| + \mathsf{hypot}\left(1, x\right)\right), x\right)} \]

   4. Taylor expanded in x around 0 10.4%

      \[\leadsto \mathsf{copysign}\left(\color{blue}{-0.041666666666666664 \cdot \left(\left(3 \cdot \frac{1}{1 + \left|x\right|} + 3 \cdot \frac{1}{{\left(1 + \left|x\right|\right)}^{2}}\right) \cdot {x}^{4}\right) + \left(0.001388888888888889 \cdot \left(\left(45 \cdot \frac{1}{1 + \left|x\right|} + \left(30 \cdot \frac{1}{{\left(1 + \left|x\right|\right)}^{3}} + 45 \cdot \frac{1}{{\left(1 + \left|x\right|\right)}^{2}}\right)\right) \cdot {x}^{6}\right) + \left(0.5 \cdot \frac{{x}^{2}}{1 + \left|x\right|} + \log \left(1 + \left|x\right|\right)\right)\right)}, x\right) \]
   5. Step-by-step derivation

      1. fma-def 10.4%

         \[\leadsto \mathsf{copysign}\left(\color{blue}{\mathsf{fma}\left(-0.041666666666666664, \left(3 \cdot \frac{1}{1 + \left|x\right|} + 3 \cdot \frac{1}{{\left(1 + \left|x\right|\right)}^{2}}\right) \cdot {x}^{4}, 0.001388888888888889 \cdot \left(\left(45 \cdot \frac{1}{1 + \left|x\right|} + \left(30 \cdot \frac{1}{{\left(1 + \left|x\right|\right)}^{3}} + 45 \cdot \frac{1}{{\left(1 + \left|x\right|\right)}^{2}}\right)\right) \cdot {x}^{6}\right) + \left(0.5 \cdot \frac{{x}^{2}}{1 + \left|x\right|} + \log \left(1 + \left|x\right|\right)\right)\right)}, x\right) \]

      2. associate-*r/ 10.4%

         \[\leadsto \mathsf{copysign}\left(\mathsf{fma}\left(-0.041666666666666664, \left(\color{blue}{\frac{3 \cdot 1}{1 + \left|x\right|}} + 3 \cdot \frac{1}{{\left(1 + \left|x\right|\right)}^{2}}\right) \cdot {x}^{4}, 0.001388888888888889 \cdot \left(\left(45 \cdot \frac{1}{1 + \left|x\right|} + \left(30 \cdot \frac{1}{{\left(1 + \left|x\right|\right)}^{3}} + 45 \cdot \frac{1}{{\left(1 + \left|x\right|\right)}^{2}}\right)\right) \cdot {x}^{6}\right) + \left(0.5 \cdot \frac{{x}^{2}}{1 + \left|x\right|} + \log \left(1 + \left|x\right|\right)\right)\right), x\right) \]

      3. metadata-eval 10.4%

         \[\leadsto \mathsf{copysign}\left(\mathsf{fma}\left(-0.041666666666666664, \left(\frac{\color{blue}{3}}{1 + \left|x\right|} + 3 \cdot \frac{1}{{\left(1 + \left|x\right|\right)}^{2}}\right) \cdot {x}^{4}, 0.001388888888888889 \cdot \left(\left(45 \cdot \frac{1}{1 + \left|x\right|} + \left(30 \cdot \frac{1}{{\left(1 + \left|x\right|\right)}^{3}} + 45 \cdot \frac{1}{{\left(1 + \left|x\right|\right)}^{2}}\right)\right) \cdot {x}^{6}\right) + \left(0.5 \cdot \frac{{x}^{2}}{1 + \left|x\right|} + \log \left(1 + \left|x\right|\right)\right)\right), x\right) \]

      4. associate-*r/ 10.4%

         \[\leadsto \mathsf{copysign}\left(\mathsf{fma}\left(-0.041666666666666664, \left(\frac{3}{1 + \left|x\right|} + \color{blue}{\frac{3 \cdot 1}{{\left(1 + \left|x\right|\right)}^{2}}}\right) \cdot {x}^{4}, 0.001388888888888889 \cdot \left(\left(45 \cdot \frac{1}{1 + \left|x\right|} + \left(30 \cdot \frac{1}{{\left(1 + \left|x\right|\right)}^{3}} + 45 \cdot \frac{1}{{\left(1 + \left|x\right|\right)}^{2}}\right)\right) \cdot {x}^{6}\right) + \left(0.5 \cdot \frac{{x}^{2}}{1 + \left|x\right|} + \log \left(1 + \left|x\right|\right)\right)\right), x\right) \]

      5. metadata-eval 10.4%

         \[\leadsto \mathsf{copysign}\left(\mathsf{fma}\left(-0.041666666666666664, \left(\frac{3}{1 + \left|x\right|} + \frac{\color{blue}{3}}{{\left(1 + \left|x\right|\right)}^{2}}\right) \cdot {x}^{4}, 0.001388888888888889 \cdot \left(\left(45 \cdot \frac{1}{1 + \left|x\right|} + \left(30 \cdot \frac{1}{{\left(1 + \left|x\right|\right)}^{3}} + 45 \cdot \frac{1}{{\left(1 + \left|x\right|\right)}^{2}}\right)\right) \cdot {x}^{6}\right) + \left(0.5 \cdot \frac{{x}^{2}}{1 + \left|x\right|} + \log \left(1 + \left|x\right|\right)\right)\right), x\right) \]

   6. Simplified 100.0%

      \[\leadsto \mathsf{copysign}\left(\color{blue}{\mathsf{fma}\left(-0.041666666666666664, \left(\frac{3}{1 + \left|x\right|} + \frac{3}{{\left(1 + \left|x\right|\right)}^{2}}\right) \cdot {x}^{4}, \mathsf{fma}\left(0.001388888888888889, \left(\frac{45}{1 + \left|x\right|} + \left(\frac{30}{{\left(1 + \left|x\right|\right)}^{3}} + \frac{45}{{\left(1 + \left|x\right|\right)}^{2}}\right)\right) \cdot {x}^{6}, \mathsf{fma}\left(0.5, \frac{x \cdot x}{1 + \left|x\right|}, \mathsf{log1p}\left(\left|x\right|\right)\right)\right)\right)}, x\right) \]

   if -0.0 < (copysign.f64 (log.f64 (+.f64 (fabs.f64 x) (sqrt.f64 (+.f64 (*.f64 x x) 1)))) x)

   1. Initial program 52.3%

      \[\mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{x \cdot x + 1}\right), x\right) \]
   2. Step-by-step derivation

      1. +-commutative 52.3%

         \[\leadsto \mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{\color{blue}{1 + x \cdot x}}\right), x\right) \]

      2. hypot-1-def 99.9%

         \[\leadsto \mathsf{copysign}\left(\log \left(\left|x\right| + \color{blue}{\mathsf{hypot}\left(1, x\right)}\right), x\right) \]

   3. Simplified 99.9%

      \[\leadsto \color{blue}{\mathsf{copysign}\left(\log \left(\left|x\right| + \mathsf{hypot}\left(1, x\right)\right), x\right)} \]
   4. Step-by-step derivation

      1. *-un-lft-identity 99.9%

         \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(1 \cdot \left(\left|x\right| + \mathsf{hypot}\left(1, x\right)\right)\right)}, x\right) \]

      2. log-prod 99.9%

         \[\leadsto \mathsf{copysign}\left(\color{blue}{\log 1 + \log \left(\left|x\right| + \mathsf{hypot}\left(1, x\right)\right)}, x\right) \]

      3. metadata-eval 99.9%

         \[\leadsto \mathsf{copysign}\left(\color{blue}{0} + \log \left(\left|x\right| + \mathsf{hypot}\left(1, x\right)\right), x\right) \]

      4. *-un-lft-identity 99.9%

         \[\leadsto \mathsf{copysign}\left(0 + \log \color{blue}{\left(1 \cdot \left(\left|x\right| + \mathsf{hypot}\left(1, x\right)\right)\right)}, x\right) \]

      5. *-un-lft-identity 99.9%

         \[\leadsto \mathsf{copysign}\left(0 + \log \color{blue}{\left(\left|x\right| + \mathsf{hypot}\left(1, x\right)\right)}, x\right) \]

      6. add-sqr-sqrt 99.9%

         \[\leadsto \mathsf{copysign}\left(0 + \log \left(\left|\color{blue}{\sqrt{x} \cdot \sqrt{x}}\right| + \mathsf{hypot}\left(1, x\right)\right), x\right) \]

      7. fabs-sqr 99.9%

         \[\leadsto \mathsf{copysign}\left(0 + \log \left(\color{blue}{\sqrt{x} \cdot \sqrt{x}} + \mathsf{hypot}\left(1, x\right)\right), x\right) \]

      8. add-sqr-sqrt 99.9%

         \[\leadsto \mathsf{copysign}\left(0 + \log \left(\color{blue}{x} + \mathsf{hypot}\left(1, x\right)\right), x\right) \]

   5. Applied egg-rr 99.9%

      \[\leadsto \mathsf{copysign}\left(\color{blue}{0 + \log \left(x + \mathsf{hypot}\left(1, x\right)\right)}, x\right) \]
   6. Step-by-step derivation

      1. +-lft-identity 99.9%

         \[\leadsto \mathsf{copysign}\left(\color{blue}{\log \left(x + \mathsf{hypot}\left(1, x\right)\right)}, x\right) \]

   7. Simplified 99.9%

      \[\leadsto \mathsf{copysign}\left(\color{blue}{\log \left(x + \mathsf{hypot}\left(1, x\right)\right)}, x\right) \]
3. Recombined 3 regimes into one program.

4. Final simplification 100.0%

   \[\leadsto \begin{array}{l} \mathbf{if}\;\mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{x \cdot x + 1}\right), x\right) \leq -10:\\ \;\;\;\;\mathsf{copysign}\left(\log \left(\frac{1}{\mathsf{hypot}\left(1, x\right) - x}\right), x\right)\\ \mathbf{elif}\;\mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{x \cdot x + 1}\right), x\right) \leq 0:\\ \;\;\;\;\mathsf{copysign}\left(\mathsf{fma}\left(-0.041666666666666664, \left(\frac{3}{\left|x\right| + 1} + \frac{3}{{\left(\left|x\right| + 1\right)}^{2}}\right) \cdot {x}^{4}, \mathsf{fma}\left(0.001388888888888889, \left(\frac{45}{\left|x\right| + 1} + \left(\frac{30}{{\left(\left|x\right| + 1\right)}^{3}} + \frac{45}{{\left(\left|x\right| + 1\right)}^{2}}\right)\right) \cdot {x}^{6}, \mathsf{fma}\left(0.5, \frac{x \cdot x}{\left|x\right| + 1}, \mathsf{log1p}\left(\left|x\right|\right)\right)\right)\right), x\right)\\ \mathbf{else}:\\ \;\;\;\;\mathsf{copysign}\left(\log \left(x + \mathsf{hypot}\left(1, x\right)\right), x\right)\\ \end{array} \]

Alternative 2: 99.0% accurate, 0.3× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} t_0 := \mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{x \cdot x + 1}\right), x\right)\\ \mathbf{if}\;t_0 \leq -10:\\ \;\;\;\;\mathsf{copysign}\left(\log \left(\frac{1}{\mathsf{hypot}\left(1, x\right) - x}\right), x\right)\\ \mathbf{elif}\;t_0 \leq 0:\\ \;\;\;\;\mathsf{copysign}\left(-0.16666666666666666 \cdot {x}^{3} + \left(0.075 \cdot {x}^{5} + \left(x + {x}^{6} \cdot -0.16666666666666666\right)\right), x\right)\\ \mathbf{else}:\\ \;\;\;\;\mathsf{copysign}\left(\log \left(x + \mathsf{hypot}\left(1, x\right)\right), x\right)\\ \end{array} \end{array} \]

FPCore

(FPCore (x)
 :precision binary64
 (let* ((t_0 (copysign (log (+ (fabs x) (sqrt (+ (* x x) 1.0)))) x)))
   (if (<= t_0 -10.0)
     (copysign (log (/ 1.0 (- (hypot 1.0 x) x))) x)
     (if (<= t_0 0.0)
       (copysign
        (+
         (* -0.16666666666666666 (pow x 3.0))
         (+ (* 0.075 (pow x 5.0)) (+ x (* (pow x 6.0) -0.16666666666666666))))
        x)
       (copysign (log (+ x (hypot 1.0 x))) x)))))

C

double code(double x) {
	double t_0 = copysign(log((fabs(x) + sqrt(((x * x) + 1.0)))), x);
	double tmp;
	if (t_0 <= -10.0) {
		tmp = copysign(log((1.0 / (hypot(1.0, x) - x))), x);
	} else if (t_0 <= 0.0) {
		tmp = copysign(((-0.16666666666666666 * pow(x, 3.0)) + ((0.075 * pow(x, 5.0)) + (x + (pow(x, 6.0) * -0.16666666666666666)))), x);
	} else {
		tmp = copysign(log((x + hypot(1.0, x))), x);
	}
	return tmp;
}

Java

public static double code(double x) {
	double t_0 = Math.copySign(Math.log((Math.abs(x) + Math.sqrt(((x * x) + 1.0)))), x);
	double tmp;
	if (t_0 <= -10.0) {
		tmp = Math.copySign(Math.log((1.0 / (Math.hypot(1.0, x) - x))), x);
	} else if (t_0 <= 0.0) {
		tmp = Math.copySign(((-0.16666666666666666 * Math.pow(x, 3.0)) + ((0.075 * Math.pow(x, 5.0)) + (x + (Math.pow(x, 6.0) * -0.16666666666666666)))), x);
	} else {
		tmp = Math.copySign(Math.log((x + Math.hypot(1.0, x))), x);
	}
	return tmp;
}

Python

def code(x):
	t_0 = math.copysign(math.log((math.fabs(x) + math.sqrt(((x * x) + 1.0)))), x)
	tmp = 0
	if t_0 <= -10.0:
		tmp = math.copysign(math.log((1.0 / (math.hypot(1.0, x) - x))), x)
	elif t_0 <= 0.0:
		tmp = math.copysign(((-0.16666666666666666 * math.pow(x, 3.0)) + ((0.075 * math.pow(x, 5.0)) + (x + (math.pow(x, 6.0) * -0.16666666666666666)))), x)
	else:
		tmp = math.copysign(math.log((x + math.hypot(1.0, x))), x)
	return tmp

Julia

function code(x)
	t_0 = copysign(log(Float64(abs(x) + sqrt(Float64(Float64(x * x) + 1.0)))), x)
	tmp = 0.0
	if (t_0 <= -10.0)
		tmp = copysign(log(Float64(1.0 / Float64(hypot(1.0, x) - x))), x);
	elseif (t_0 <= 0.0)
		tmp = copysign(Float64(Float64(-0.16666666666666666 * (x ^ 3.0)) + Float64(Float64(0.075 * (x ^ 5.0)) + Float64(x + Float64((x ^ 6.0) * -0.16666666666666666)))), x);
	else
		tmp = copysign(log(Float64(x + hypot(1.0, x))), x);
	end
	return tmp
end

MATLAB

function tmp_2 = code(x)
	t_0 = sign(x) * abs(log((abs(x) + sqrt(((x * x) + 1.0)))));
	tmp = 0.0;
	if (t_0 <= -10.0)
		tmp = sign(x) * abs(log((1.0 / (hypot(1.0, x) - x))));
	elseif (t_0 <= 0.0)
		tmp = sign(x) * abs(((-0.16666666666666666 * (x ^ 3.0)) + ((0.075 * (x ^ 5.0)) + (x + ((x ^ 6.0) * -0.16666666666666666)))));
	else
		tmp = sign(x) * abs(log((x + hypot(1.0, x))));
	end
	tmp_2 = tmp;
end

Wolfram

code[x_] := Block[{t$95$0 = N[With[{TMP1 = Abs[N[Log[N[(N[Abs[x], $MachinePrecision] + N[Sqrt[N[(N[(x * x), $MachinePrecision] + 1.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], TMP2 = Sign[x]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]}, If[LessEqual[t$95$0, -10.0], N[With[{TMP1 = Abs[N[Log[N[(1.0 / N[(N[Sqrt[1.0 ^ 2 + x ^ 2], $MachinePrecision] - x), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], TMP2 = Sign[x]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision], If[LessEqual[t$95$0, 0.0], N[With[{TMP1 = Abs[N[(N[(-0.16666666666666666 * N[Power[x, 3.0], $MachinePrecision]), $MachinePrecision] + N[(N[(0.075 * N[Power[x, 5.0], $MachinePrecision]), $MachinePrecision] + N[(x + N[(N[Power[x, 6.0], $MachinePrecision] * -0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], TMP2 = Sign[x]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision], N[With[{TMP1 = Abs[N[Log[N[(x + N[Sqrt[1.0 ^ 2 + x ^ 2], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], TMP2 = Sign[x]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]]]]

Derivation

1. Split input into 3 regimes

2. if (copysign.f64 (log.f64 (+.f64 (fabs.f64 x) (sqrt.f64 (+.f64 (*.f64 x x) 1)))) x) < -10

   1. Initial program 50.8%

      \[\mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{x \cdot x + 1}\right), x\right) \]
   2. Step-by-step derivation

      1. +-commutative 50.8%

         \[\leadsto \mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{\color{blue}{1 + x \cdot x}}\right), x\right) \]

      2. hypot-1-def 100.0%

         \[\leadsto \mathsf{copysign}\left(\log \left(\left|x\right| + \color{blue}{\mathsf{hypot}\left(1, x\right)}\right), x\right) \]

   3. Simplified 100.0%

      \[\leadsto \color{blue}{\mathsf{copysign}\left(\log \left(\left|x\right| + \mathsf{hypot}\left(1, x\right)\right), x\right)} \]
   4. Step-by-step derivation

      1. flip-+ 0.9%

         \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(\frac{\left|x\right| \cdot \left|x\right| - \mathsf{hypot}\left(1, x\right) \cdot \mathsf{hypot}\left(1, x\right)}{\left|x\right| - \mathsf{hypot}\left(1, x\right)}\right)}, x\right) \]

      2. div-sub 0.9%

         \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(\frac{\left|x\right| \cdot \left|x\right|}{\left|x\right| - \mathsf{hypot}\left(1, x\right)} - \frac{\mathsf{hypot}\left(1, x\right) \cdot \mathsf{hypot}\left(1, x\right)}{\left|x\right| - \mathsf{hypot}\left(1, x\right)}\right)}, x\right) \]

      3. pow2 0.9%

         \[\leadsto \mathsf{copysign}\left(\log \left(\frac{\color{blue}{{\left(\left|x\right|\right)}^{2}}}{\left|x\right| - \mathsf{hypot}\left(1, x\right)} - \frac{\mathsf{hypot}\left(1, x\right) \cdot \mathsf{hypot}\left(1, x\right)}{\left|x\right| - \mathsf{hypot}\left(1, x\right)}\right), x\right) \]

      4. add-sqr-sqrt 0.0%

         \[\leadsto \mathsf{copysign}\left(\log \left(\frac{{\left(\left|\color{blue}{\sqrt{x} \cdot \sqrt{x}}\right|\right)}^{2}}{\left|x\right| - \mathsf{hypot}\left(1, x\right)} - \frac{\mathsf{hypot}\left(1, x\right) \cdot \mathsf{hypot}\left(1, x\right)}{\left|x\right| - \mathsf{hypot}\left(1, x\right)}\right), x\right) \]

      5. fabs-sqr 0.0%

         \[\leadsto \mathsf{copysign}\left(\log \left(\frac{{\color{blue}{\left(\sqrt{x} \cdot \sqrt{x}\right)}}^{2}}{\left|x\right| - \mathsf{hypot}\left(1, x\right)} - \frac{\mathsf{hypot}\left(1, x\right) \cdot \mathsf{hypot}\left(1, x\right)}{\left|x\right| - \mathsf{hypot}\left(1, x\right)}\right), x\right) \]

      6. add-sqr-sqrt 0.9%

         \[\leadsto \mathsf{copysign}\left(\log \left(\frac{{\color{blue}{x}}^{2}}{\left|x\right| - \mathsf{hypot}\left(1, x\right)} - \frac{\mathsf{hypot}\left(1, x\right) \cdot \mathsf{hypot}\left(1, x\right)}{\left|x\right| - \mathsf{hypot}\left(1, x\right)}\right), x\right) \]

      7. pow2 0.9%

         \[\leadsto \mathsf{copysign}\left(\log \left(\frac{\color{blue}{x \cdot x}}{\left|x\right| - \mathsf{hypot}\left(1, x\right)} - \frac{\mathsf{hypot}\left(1, x\right) \cdot \mathsf{hypot}\left(1, x\right)}{\left|x\right| - \mathsf{hypot}\left(1, x\right)}\right), x\right) \]

      8. add-sqr-sqrt 0.0%

         \[\leadsto \mathsf{copysign}\left(\log \left(\frac{x \cdot x}{\left|\color{blue}{\sqrt{x} \cdot \sqrt{x}}\right| - \mathsf{hypot}\left(1, x\right)} - \frac{\mathsf{hypot}\left(1, x\right) \cdot \mathsf{hypot}\left(1, x\right)}{\left|x\right| - \mathsf{hypot}\left(1, x\right)}\right), x\right) \]

      9. fabs-sqr 0.0%

         \[\leadsto \mathsf{copysign}\left(\log \left(\frac{x \cdot x}{\color{blue}{\sqrt{x} \cdot \sqrt{x}} - \mathsf{hypot}\left(1, x\right)} - \frac{\mathsf{hypot}\left(1, x\right) \cdot \mathsf{hypot}\left(1, x\right)}{\left|x\right| - \mathsf{hypot}\left(1, x\right)}\right), x\right) \]

      10. add-sqr-sqrt 0.3%

          \[\leadsto \mathsf{copysign}\left(\log \left(\frac{x \cdot x}{\color{blue}{x} - \mathsf{hypot}\left(1, x\right)} - \frac{\mathsf{hypot}\left(1, x\right) \cdot \mathsf{hypot}\left(1, x\right)}{\left|x\right| - \mathsf{hypot}\left(1, x\right)}\right), x\right) \]

      11. hypot-udef 0.3%

          \[\leadsto \mathsf{copysign}\left(\log \left(\frac{x \cdot x}{x - \mathsf{hypot}\left(1, x\right)} - \frac{\color{blue}{\sqrt{1 \cdot 1 + x \cdot x}} \cdot \mathsf{hypot}\left(1, x\right)}{\left|x\right| - \mathsf{hypot}\left(1, x\right)}\right), x\right) \]

      12. hypot-udef 0.3%

          \[\leadsto \mathsf{copysign}\left(\log \left(\frac{x \cdot x}{x - \mathsf{hypot}\left(1, x\right)} - \frac{\sqrt{1 \cdot 1 + x \cdot x} \cdot \color{blue}{\sqrt{1 \cdot 1 + x \cdot x}}}{\left|x\right| - \mathsf{hypot}\left(1, x\right)}\right), x\right) \]

      13. add-sqr-sqrt 0.3%

          \[\leadsto \mathsf{copysign}\left(\log \left(\frac{x \cdot x}{x - \mathsf{hypot}\left(1, x\right)} - \frac{\color{blue}{1 \cdot 1 + x \cdot x}}{\left|x\right| - \mathsf{hypot}\left(1, x\right)}\right), x\right) \]

      14. metadata-eval 0.3%

          \[\leadsto \mathsf{copysign}\left(\log \left(\frac{x \cdot x}{x - \mathsf{hypot}\left(1, x\right)} - \frac{\color{blue}{1} + x \cdot x}{\left|x\right| - \mathsf{hypot}\left(1, x\right)}\right), x\right) \]

      15. add-sqr-sqrt 0.0%

          \[\leadsto \mathsf{copysign}\left(\log \left(\frac{x \cdot x}{x - \mathsf{hypot}\left(1, x\right)} - \frac{1 + x \cdot x}{\left|\color{blue}{\sqrt{x} \cdot \sqrt{x}}\right| - \mathsf{hypot}\left(1, x\right)}\right), x\right) \]

      16. fabs-sqr 0.0%

          \[\leadsto \mathsf{copysign}\left(\log \left(\frac{x \cdot x}{x - \mathsf{hypot}\left(1, x\right)} - \frac{1 + x \cdot x}{\color{blue}{\sqrt{x} \cdot \sqrt{x}} - \mathsf{hypot}\left(1, x\right)}\right), x\right) \]

      17. add-sqr-sqrt 2.5%

          \[\leadsto \mathsf{copysign}\left(\log \left(\frac{x \cdot x}{x - \mathsf{hypot}\left(1, x\right)} - \frac{1 + x \cdot x}{\color{blue}{x} - \mathsf{hypot}\left(1, x\right)}\right), x\right) \]

   5. Applied egg-rr 2.5%

      \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(\frac{x \cdot x}{x - \mathsf{hypot}\left(1, x\right)} - \frac{1 + x \cdot x}{x - \mathsf{hypot}\left(1, x\right)}\right)}, x\right) \]
   6. Step-by-step derivation

      1. unpow2 2.5%

         \[\leadsto \mathsf{copysign}\left(\log \left(\frac{\color{blue}{{x}^{2}}}{x - \mathsf{hypot}\left(1, x\right)} - \frac{1 + x \cdot x}{x - \mathsf{hypot}\left(1, x\right)}\right), x\right) \]

      2. div-sub 3.0%

         \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(\frac{{x}^{2} - \left(1 + x \cdot x\right)}{x - \mathsf{hypot}\left(1, x\right)}\right)}, x\right) \]

      3. unpow2 3.0%

         \[\leadsto \mathsf{copysign}\left(\log \left(\frac{\color{blue}{x \cdot x} - \left(1 + x \cdot x\right)}{x - \mathsf{hypot}\left(1, x\right)}\right), x\right) \]

      4. unpow2 3.0%

         \[\leadsto \mathsf{copysign}\left(\log \left(\frac{\color{blue}{{x}^{2}} - \left(1 + x \cdot x\right)}{x - \mathsf{hypot}\left(1, x\right)}\right), x\right) \]

      5. unpow2 3.0%

         \[\leadsto \mathsf{copysign}\left(\log \left(\frac{{x}^{2} - \left(1 + \color{blue}{{x}^{2}}\right)}{x - \mathsf{hypot}\left(1, x\right)}\right), x\right) \]

      6. +-commutative 3.0%

         \[\leadsto \mathsf{copysign}\left(\log \left(\frac{{x}^{2} - \color{blue}{\left({x}^{2} + 1\right)}}{x - \mathsf{hypot}\left(1, x\right)}\right), x\right) \]

      7. associate--r+ 49.2%

         \[\leadsto \mathsf{copysign}\left(\log \left(\frac{\color{blue}{\left({x}^{2} - {x}^{2}\right) - 1}}{x - \mathsf{hypot}\left(1, x\right)}\right), x\right) \]

      8. +-inverses 100.0%

         \[\leadsto \mathsf{copysign}\left(\log \left(\frac{\color{blue}{0} - 1}{x - \mathsf{hypot}\left(1, x\right)}\right), x\right) \]

      9. metadata-eval 100.0%

         \[\leadsto \mathsf{copysign}\left(\log \left(\frac{\color{blue}{-1}}{x - \mathsf{hypot}\left(1, x\right)}\right), x\right) \]

      10. metadata-eval 100.0%

          \[\leadsto \mathsf{copysign}\left(\log \left(\frac{\color{blue}{\frac{1}{-1}}}{x - \mathsf{hypot}\left(1, x\right)}\right), x\right) \]

      11. associate-/r* 100.0%

          \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(\frac{1}{-1 \cdot \left(x - \mathsf{hypot}\left(1, x\right)\right)}\right)}, x\right) \]

      12. neg-mul-1 100.0%

          \[\leadsto \mathsf{copysign}\left(\log \left(\frac{1}{\color{blue}{-\left(x - \mathsf{hypot}\left(1, x\right)\right)}}\right), x\right) \]

      13. sub-neg 100.0%

          \[\leadsto \mathsf{copysign}\left(\log \left(\frac{1}{-\color{blue}{\left(x + \left(-\mathsf{hypot}\left(1, x\right)\right)\right)}}\right), x\right) \]

      14. +-commutative 100.0%

          \[\leadsto \mathsf{copysign}\left(\log \left(\frac{1}{-\color{blue}{\left(\left(-\mathsf{hypot}\left(1, x\right)\right) + x\right)}}\right), x\right) \]

      15. distribute-neg-in 100.0%

          \[\leadsto \mathsf{copysign}\left(\log \left(\frac{1}{\color{blue}{\left(-\left(-\mathsf{hypot}\left(1, x\right)\right)\right) + \left(-x\right)}}\right), x\right) \]

      16. remove-double-neg 100.0%

          \[\leadsto \mathsf{copysign}\left(\log \left(\frac{1}{\color{blue}{\mathsf{hypot}\left(1, x\right)} + \left(-x\right)}\right), x\right) \]

      17. sub-neg 100.0%

          \[\leadsto \mathsf{copysign}\left(\log \left(\frac{1}{\color{blue}{\mathsf{hypot}\left(1, x\right) - x}}\right), x\right) \]

   7. Simplified 100.0%

      \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(\frac{1}{\mathsf{hypot}\left(1, x\right) - x}\right)}, x\right) \]

   if -10 < (copysign.f64 (log.f64 (+.f64 (fabs.f64 x) (sqrt.f64 (+.f64 (*.f64 x x) 1)))) x) < -0.0

   1. Initial program 9.4%

      \[\mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{x \cdot x + 1}\right), x\right) \]
   2. Step-by-step derivation

      1. +-commutative 9.4%

         \[\leadsto \mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{\color{blue}{1 + x \cdot x}}\right), x\right) \]

      2. hypot-1-def 9.5%

         \[\leadsto \mathsf{copysign}\left(\log \left(\left|x\right| + \color{blue}{\mathsf{hypot}\left(1, x\right)}\right), x\right) \]

   3. Simplified 9.5%

      \[\leadsto \color{blue}{\mathsf{copysign}\left(\log \left(\left|x\right| + \mathsf{hypot}\left(1, x\right)\right), x\right)} \]

   4. Taylor expanded in x around 0 10.4%

      \[\leadsto \mathsf{copysign}\left(\color{blue}{-0.041666666666666664 \cdot \left(\left(3 \cdot \frac{1}{1 + \left|x\right|} + 3 \cdot \frac{1}{{\left(1 + \left|x\right|\right)}^{2}}\right) \cdot {x}^{4}\right) + \left(0.5 \cdot \frac{{x}^{2}}{1 + \left|x\right|} + \log \left(1 + \left|x\right|\right)\right)}, x\right) \]
   5. Step-by-step derivation

      1. associate-*r* 10.4%

         \[\leadsto \mathsf{copysign}\left(\color{blue}{\left(-0.041666666666666664 \cdot \left(3 \cdot \frac{1}{1 + \left|x\right|} + 3 \cdot \frac{1}{{\left(1 + \left|x\right|\right)}^{2}}\right)\right) \cdot {x}^{4}} + \left(0.5 \cdot \frac{{x}^{2}}{1 + \left|x\right|} + \log \left(1 + \left|x\right|\right)\right), x\right) \]

      2. *-commutative 10.4%

         \[\leadsto \mathsf{copysign}\left(\color{blue}{{x}^{4} \cdot \left(-0.041666666666666664 \cdot \left(3 \cdot \frac{1}{1 + \left|x\right|} + 3 \cdot \frac{1}{{\left(1 + \left|x\right|\right)}^{2}}\right)\right)} + \left(0.5 \cdot \frac{{x}^{2}}{1 + \left|x\right|} + \log \left(1 + \left|x\right|\right)\right), x\right) \]

      3. fma-def 10.4%

         \[\leadsto \mathsf{copysign}\left(\color{blue}{\mathsf{fma}\left({x}^{4}, -0.041666666666666664 \cdot \left(3 \cdot \frac{1}{1 + \left|x\right|} + 3 \cdot \frac{1}{{\left(1 + \left|x\right|\right)}^{2}}\right), 0.5 \cdot \frac{{x}^{2}}{1 + \left|x\right|} + \log \left(1 + \left|x\right|\right)\right)}, x\right) \]

   6. Simplified 100.0%

      \[\leadsto \mathsf{copysign}\left(\color{blue}{\mathsf{fma}\left({x}^{4}, \frac{-0.125}{1 + x} + \frac{-0.125}{{\left(1 + x\right)}^{2}}, \mathsf{fma}\left(0.5, \frac{x \cdot x}{1 + x}, \mathsf{log1p}\left(x\right)\right)\right)}, x\right) \]

   7. Taylor expanded in x around 0 100.0%

      \[\leadsto \mathsf{copysign}\left(\color{blue}{-0.16666666666666666 \cdot {x}^{3} + \left(0.075 \cdot {x}^{5} + \left(-0.16666666666666666 \cdot {x}^{6} + x\right)\right)}, x\right) \]

   if -0.0 < (copysign.f64 (log.f64 (+.f64 (fabs.f64 x) (sqrt.f64 (+.f64 (*.f64 x x) 1)))) x)

   1. Initial program 52.3%

      \[\mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{x \cdot x + 1}\right), x\right) \]
   2. Step-by-step derivation

      1. +-commutative 52.3%

         \[\leadsto \mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{\color{blue}{1 + x \cdot x}}\right), x\right) \]

      2. hypot-1-def 99.9%

         \[\leadsto \mathsf{copysign}\left(\log \left(\left|x\right| + \color{blue}{\mathsf{hypot}\left(1, x\right)}\right), x\right) \]

   3. Simplified 99.9%

      \[\leadsto \color{blue}{\mathsf{copysign}\left(\log \left(\left|x\right| + \mathsf{hypot}\left(1, x\right)\right), x\right)} \]
   4. Step-by-step derivation

      1. *-un-lft-identity 99.9%

         \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(1 \cdot \left(\left|x\right| + \mathsf{hypot}\left(1, x\right)\right)\right)}, x\right) \]

      2. log-prod 99.9%

         \[\leadsto \mathsf{copysign}\left(\color{blue}{\log 1 + \log \left(\left|x\right| + \mathsf{hypot}\left(1, x\right)\right)}, x\right) \]

      3. metadata-eval 99.9%

         \[\leadsto \mathsf{copysign}\left(\color{blue}{0} + \log \left(\left|x\right| + \mathsf{hypot}\left(1, x\right)\right), x\right) \]

      4. *-un-lft-identity 99.9%

         \[\leadsto \mathsf{copysign}\left(0 + \log \color{blue}{\left(1 \cdot \left(\left|x\right| + \mathsf{hypot}\left(1, x\right)\right)\right)}, x\right) \]

      5. *-un-lft-identity 99.9%

         \[\leadsto \mathsf{copysign}\left(0 + \log \color{blue}{\left(\left|x\right| + \mathsf{hypot}\left(1, x\right)\right)}, x\right) \]

      6. add-sqr-sqrt 99.9%

         \[\leadsto \mathsf{copysign}\left(0 + \log \left(\left|\color{blue}{\sqrt{x} \cdot \sqrt{x}}\right| + \mathsf{hypot}\left(1, x\right)\right), x\right) \]

      7. fabs-sqr 99.9%

         \[\leadsto \mathsf{copysign}\left(0 + \log \left(\color{blue}{\sqrt{x} \cdot \sqrt{x}} + \mathsf{hypot}\left(1, x\right)\right), x\right) \]

      8. add-sqr-sqrt 99.9%

         \[\leadsto \mathsf{copysign}\left(0 + \log \left(\color{blue}{x} + \mathsf{hypot}\left(1, x\right)\right), x\right) \]

   5. Applied egg-rr 99.9%

      \[\leadsto \mathsf{copysign}\left(\color{blue}{0 + \log \left(x + \mathsf{hypot}\left(1, x\right)\right)}, x\right) \]
   6. Step-by-step derivation

      1. +-lft-identity 99.9%

         \[\leadsto \mathsf{copysign}\left(\color{blue}{\log \left(x + \mathsf{hypot}\left(1, x\right)\right)}, x\right) \]

   7. Simplified 99.9%

      \[\leadsto \mathsf{copysign}\left(\color{blue}{\log \left(x + \mathsf{hypot}\left(1, x\right)\right)}, x\right) \]
3. Recombined 3 regimes into one program.

4. Final simplification 100.0%

   \[\leadsto \begin{array}{l} \mathbf{if}\;\mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{x \cdot x + 1}\right), x\right) \leq -10:\\ \;\;\;\;\mathsf{copysign}\left(\log \left(\frac{1}{\mathsf{hypot}\left(1, x\right) - x}\right), x\right)\\ \mathbf{elif}\;\mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{x \cdot x + 1}\right), x\right) \leq 0:\\ \;\;\;\;\mathsf{copysign}\left(-0.16666666666666666 \cdot {x}^{3} + \left(0.075 \cdot {x}^{5} + \left(x + {x}^{6} \cdot -0.16666666666666666\right)\right), x\right)\\ \mathbf{else}:\\ \;\;\;\;\mathsf{copysign}\left(\log \left(x + \mathsf{hypot}\left(1, x\right)\right), x\right)\\ \end{array} \]

Alternative 3: 99.4% accurate, 0.4× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} t_0 := \mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{x \cdot x + 1}\right), x\right)\\ \mathbf{if}\;t_0 \leq -0.0005:\\ \;\;\;\;\mathsf{copysign}\left(\log \left(\frac{1}{\mathsf{hypot}\left(1, x\right) - x}\right), x\right)\\ \mathbf{elif}\;t_0 \leq 0:\\ \;\;\;\;\mathsf{copysign}\left(0.5 \cdot \frac{x}{\frac{x + 1}{x}} + \mathsf{log1p}\left(x\right), x\right)\\ \mathbf{else}:\\ \;\;\;\;\mathsf{copysign}\left(\log \left(x + \mathsf{hypot}\left(1, x\right)\right), x\right)\\ \end{array} \end{array} \]

FPCore

(FPCore (x)
 :precision binary64
 (let* ((t_0 (copysign (log (+ (fabs x) (sqrt (+ (* x x) 1.0)))) x)))
   (if (<= t_0 -0.0005)
     (copysign (log (/ 1.0 (- (hypot 1.0 x) x))) x)
     (if (<= t_0 0.0)
       (copysign (+ (* 0.5 (/ x (/ (+ x 1.0) x))) (log1p x)) x)
       (copysign (log (+ x (hypot 1.0 x))) x)))))

C

double code(double x) {
	double t_0 = copysign(log((fabs(x) + sqrt(((x * x) + 1.0)))), x);
	double tmp;
	if (t_0 <= -0.0005) {
		tmp = copysign(log((1.0 / (hypot(1.0, x) - x))), x);
	} else if (t_0 <= 0.0) {
		tmp = copysign(((0.5 * (x / ((x + 1.0) / x))) + log1p(x)), x);
	} else {
		tmp = copysign(log((x + hypot(1.0, x))), x);
	}
	return tmp;
}

Java

public static double code(double x) {
	double t_0 = Math.copySign(Math.log((Math.abs(x) + Math.sqrt(((x * x) + 1.0)))), x);
	double tmp;
	if (t_0 <= -0.0005) {
		tmp = Math.copySign(Math.log((1.0 / (Math.hypot(1.0, x) - x))), x);
	} else if (t_0 <= 0.0) {
		tmp = Math.copySign(((0.5 * (x / ((x + 1.0) / x))) + Math.log1p(x)), x);
	} else {
		tmp = Math.copySign(Math.log((x + Math.hypot(1.0, x))), x);
	}
	return tmp;
}

Python

def code(x):
	t_0 = math.copysign(math.log((math.fabs(x) + math.sqrt(((x * x) + 1.0)))), x)
	tmp = 0
	if t_0 <= -0.0005:
		tmp = math.copysign(math.log((1.0 / (math.hypot(1.0, x) - x))), x)
	elif t_0 <= 0.0:
		tmp = math.copysign(((0.5 * (x / ((x + 1.0) / x))) + math.log1p(x)), x)
	else:
		tmp = math.copysign(math.log((x + math.hypot(1.0, x))), x)
	return tmp

Julia

function code(x)
	t_0 = copysign(log(Float64(abs(x) + sqrt(Float64(Float64(x * x) + 1.0)))), x)
	tmp = 0.0
	if (t_0 <= -0.0005)
		tmp = copysign(log(Float64(1.0 / Float64(hypot(1.0, x) - x))), x);
	elseif (t_0 <= 0.0)
		tmp = copysign(Float64(Float64(0.5 * Float64(x / Float64(Float64(x + 1.0) / x))) + log1p(x)), x);
	else
		tmp = copysign(log(Float64(x + hypot(1.0, x))), x);
	end
	return tmp
end

Wolfram

code[x_] := Block[{t$95$0 = N[With[{TMP1 = Abs[N[Log[N[(N[Abs[x], $MachinePrecision] + N[Sqrt[N[(N[(x * x), $MachinePrecision] + 1.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], TMP2 = Sign[x]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]}, If[LessEqual[t$95$0, -0.0005], N[With[{TMP1 = Abs[N[Log[N[(1.0 / N[(N[Sqrt[1.0 ^ 2 + x ^ 2], $MachinePrecision] - x), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], TMP2 = Sign[x]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision], If[LessEqual[t$95$0, 0.0], N[With[{TMP1 = Abs[N[(N[(0.5 * N[(x / N[(N[(x + 1.0), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[Log[1 + x], $MachinePrecision]), $MachinePrecision]], TMP2 = Sign[x]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision], N[With[{TMP1 = Abs[N[Log[N[(x + N[Sqrt[1.0 ^ 2 + x ^ 2], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], TMP2 = Sign[x]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]]]]

Derivation

1. Split input into 3 regimes

2. if (copysign.f64 (log.f64 (+.f64 (fabs.f64 x) (sqrt.f64 (+.f64 (*.f64 x x) 1)))) x) < -5.0000000000000001e-4

   1. Initial program 51.9%

      \[\mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{x \cdot x + 1}\right), x\right) \]
   2. Step-by-step derivation

      1. +-commutative 51.9%

         \[\leadsto \mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{\color{blue}{1 + x \cdot x}}\right), x\right) \]

      2. hypot-1-def 99.6%

         \[\leadsto \mathsf{copysign}\left(\log \left(\left|x\right| + \color{blue}{\mathsf{hypot}\left(1, x\right)}\right), x\right) \]

   3. Simplified 99.6%

      \[\leadsto \color{blue}{\mathsf{copysign}\left(\log \left(\left|x\right| + \mathsf{hypot}\left(1, x\right)\right), x\right)} \]
   4. Step-by-step derivation

      1. flip-+ 3.4%

         \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(\frac{\left|x\right| \cdot \left|x\right| - \mathsf{hypot}\left(1, x\right) \cdot \mathsf{hypot}\left(1, x\right)}{\left|x\right| - \mathsf{hypot}\left(1, x\right)}\right)}, x\right) \]

      2. div-sub 3.4%

         \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(\frac{\left|x\right| \cdot \left|x\right|}{\left|x\right| - \mathsf{hypot}\left(1, x\right)} - \frac{\mathsf{hypot}\left(1, x\right) \cdot \mathsf{hypot}\left(1, x\right)}{\left|x\right| - \mathsf{hypot}\left(1, x\right)}\right)}, x\right) \]

      3. pow2 3.4%

         \[\leadsto \mathsf{copysign}\left(\log \left(\frac{\color{blue}{{\left(\left|x\right|\right)}^{2}}}{\left|x\right| - \mathsf{hypot}\left(1, x\right)} - \frac{\mathsf{hypot}\left(1, x\right) \cdot \mathsf{hypot}\left(1, x\right)}{\left|x\right| - \mathsf{hypot}\left(1, x\right)}\right), x\right) \]

      4. add-sqr-sqrt 0.0%

         \[\leadsto \mathsf{copysign}\left(\log \left(\frac{{\left(\left|\color{blue}{\sqrt{x} \cdot \sqrt{x}}\right|\right)}^{2}}{\left|x\right| - \mathsf{hypot}\left(1, x\right)} - \frac{\mathsf{hypot}\left(1, x\right) \cdot \mathsf{hypot}\left(1, x\right)}{\left|x\right| - \mathsf{hypot}\left(1, x\right)}\right), x\right) \]

      5. fabs-sqr 0.0%

         \[\leadsto \mathsf{copysign}\left(\log \left(\frac{{\color{blue}{\left(\sqrt{x} \cdot \sqrt{x}\right)}}^{2}}{\left|x\right| - \mathsf{hypot}\left(1, x\right)} - \frac{\mathsf{hypot}\left(1, x\right) \cdot \mathsf{hypot}\left(1, x\right)}{\left|x\right| - \mathsf{hypot}\left(1, x\right)}\right), x\right) \]

      6. add-sqr-sqrt 3.4%

         \[\leadsto \mathsf{copysign}\left(\log \left(\frac{{\color{blue}{x}}^{2}}{\left|x\right| - \mathsf{hypot}\left(1, x\right)} - \frac{\mathsf{hypot}\left(1, x\right) \cdot \mathsf{hypot}\left(1, x\right)}{\left|x\right| - \mathsf{hypot}\left(1, x\right)}\right), x\right) \]

      7. pow2 3.4%

         \[\leadsto \mathsf{copysign}\left(\log \left(\frac{\color{blue}{x \cdot x}}{\left|x\right| - \mathsf{hypot}\left(1, x\right)} - \frac{\mathsf{hypot}\left(1, x\right) \cdot \mathsf{hypot}\left(1, x\right)}{\left|x\right| - \mathsf{hypot}\left(1, x\right)}\right), x\right) \]

      8. add-sqr-sqrt 0.0%

         \[\leadsto \mathsf{copysign}\left(\log \left(\frac{x \cdot x}{\left|\color{blue}{\sqrt{x} \cdot \sqrt{x}}\right| - \mathsf{hypot}\left(1, x\right)} - \frac{\mathsf{hypot}\left(1, x\right) \cdot \mathsf{hypot}\left(1, x\right)}{\left|x\right| - \mathsf{hypot}\left(1, x\right)}\right), x\right) \]

      9. fabs-sqr 0.0%

         \[\leadsto \mathsf{copysign}\left(\log \left(\frac{x \cdot x}{\color{blue}{\sqrt{x} \cdot \sqrt{x}} - \mathsf{hypot}\left(1, x\right)} - \frac{\mathsf{hypot}\left(1, x\right) \cdot \mathsf{hypot}\left(1, x\right)}{\left|x\right| - \mathsf{hypot}\left(1, x\right)}\right), x\right) \]

      10. add-sqr-sqrt 1.7%

          \[\leadsto \mathsf{copysign}\left(\log \left(\frac{x \cdot x}{\color{blue}{x} - \mathsf{hypot}\left(1, x\right)} - \frac{\mathsf{hypot}\left(1, x\right) \cdot \mathsf{hypot}\left(1, x\right)}{\left|x\right| - \mathsf{hypot}\left(1, x\right)}\right), x\right) \]

      11. hypot-udef 1.7%

          \[\leadsto \mathsf{copysign}\left(\log \left(\frac{x \cdot x}{x - \mathsf{hypot}\left(1, x\right)} - \frac{\color{blue}{\sqrt{1 \cdot 1 + x \cdot x}} \cdot \mathsf{hypot}\left(1, x\right)}{\left|x\right| - \mathsf{hypot}\left(1, x\right)}\right), x\right) \]

      12. hypot-udef 1.7%

          \[\leadsto \mathsf{copysign}\left(\log \left(\frac{x \cdot x}{x - \mathsf{hypot}\left(1, x\right)} - \frac{\sqrt{1 \cdot 1 + x \cdot x} \cdot \color{blue}{\sqrt{1 \cdot 1 + x \cdot x}}}{\left|x\right| - \mathsf{hypot}\left(1, x\right)}\right), x\right) \]

      13. add-sqr-sqrt 1.7%

          \[\leadsto \mathsf{copysign}\left(\log \left(\frac{x \cdot x}{x - \mathsf{hypot}\left(1, x\right)} - \frac{\color{blue}{1 \cdot 1 + x \cdot x}}{\left|x\right| - \mathsf{hypot}\left(1, x\right)}\right), x\right) \]

      14. metadata-eval 1.7%

          \[\leadsto \mathsf{copysign}\left(\log \left(\frac{x \cdot x}{x - \mathsf{hypot}\left(1, x\right)} - \frac{\color{blue}{1} + x \cdot x}{\left|x\right| - \mathsf{hypot}\left(1, x\right)}\right), x\right) \]

      15. add-sqr-sqrt 0.0%

          \[\leadsto \mathsf{copysign}\left(\log \left(\frac{x \cdot x}{x - \mathsf{hypot}\left(1, x\right)} - \frac{1 + x \cdot x}{\left|\color{blue}{\sqrt{x} \cdot \sqrt{x}}\right| - \mathsf{hypot}\left(1, x\right)}\right), x\right) \]

      16. fabs-sqr 0.0%

          \[\leadsto \mathsf{copysign}\left(\log \left(\frac{x \cdot x}{x - \mathsf{hypot}\left(1, x\right)} - \frac{1 + x \cdot x}{\color{blue}{\sqrt{x} \cdot \sqrt{x}} - \mathsf{hypot}\left(1, x\right)}\right), x\right) \]

      17. add-sqr-sqrt 4.9%

          \[\leadsto \mathsf{copysign}\left(\log \left(\frac{x \cdot x}{x - \mathsf{hypot}\left(1, x\right)} - \frac{1 + x \cdot x}{\color{blue}{x} - \mathsf{hypot}\left(1, x\right)}\right), x\right) \]

   5. Applied egg-rr 4.9%

      \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(\frac{x \cdot x}{x - \mathsf{hypot}\left(1, x\right)} - \frac{1 + x \cdot x}{x - \mathsf{hypot}\left(1, x\right)}\right)}, x\right) \]
   6. Step-by-step derivation

      1. unpow2 4.9%

         \[\leadsto \mathsf{copysign}\left(\log \left(\frac{\color{blue}{{x}^{2}}}{x - \mathsf{hypot}\left(1, x\right)} - \frac{1 + x \cdot x}{x - \mathsf{hypot}\left(1, x\right)}\right), x\right) \]

      2. div-sub 5.5%

         \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(\frac{{x}^{2} - \left(1 + x \cdot x\right)}{x - \mathsf{hypot}\left(1, x\right)}\right)}, x\right) \]

      3. unpow2 5.5%

         \[\leadsto \mathsf{copysign}\left(\log \left(\frac{\color{blue}{x \cdot x} - \left(1 + x \cdot x\right)}{x - \mathsf{hypot}\left(1, x\right)}\right), x\right) \]

      4. unpow2 5.5%

         \[\leadsto \mathsf{copysign}\left(\log \left(\frac{\color{blue}{{x}^{2}} - \left(1 + x \cdot x\right)}{x - \mathsf{hypot}\left(1, x\right)}\right), x\right) \]

      5. unpow2 5.5%

         \[\leadsto \mathsf{copysign}\left(\log \left(\frac{{x}^{2} - \left(1 + \color{blue}{{x}^{2}}\right)}{x - \mathsf{hypot}\left(1, x\right)}\right), x\right) \]

      6. +-commutative 5.5%

         \[\leadsto \mathsf{copysign}\left(\log \left(\frac{{x}^{2} - \color{blue}{\left({x}^{2} + 1\right)}}{x - \mathsf{hypot}\left(1, x\right)}\right), x\right) \]

      7. associate--r+ 50.3%

         \[\leadsto \mathsf{copysign}\left(\log \left(\frac{\color{blue}{\left({x}^{2} - {x}^{2}\right) - 1}}{x - \mathsf{hypot}\left(1, x\right)}\right), x\right) \]

      8. +-inverses 99.6%

         \[\leadsto \mathsf{copysign}\left(\log \left(\frac{\color{blue}{0} - 1}{x - \mathsf{hypot}\left(1, x\right)}\right), x\right) \]

      9. metadata-eval 99.6%

         \[\leadsto \mathsf{copysign}\left(\log \left(\frac{\color{blue}{-1}}{x - \mathsf{hypot}\left(1, x\right)}\right), x\right) \]

      10. metadata-eval 99.6%

          \[\leadsto \mathsf{copysign}\left(\log \left(\frac{\color{blue}{\frac{1}{-1}}}{x - \mathsf{hypot}\left(1, x\right)}\right), x\right) \]

      11. associate-/r* 99.6%

          \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(\frac{1}{-1 \cdot \left(x - \mathsf{hypot}\left(1, x\right)\right)}\right)}, x\right) \]

      12. neg-mul-1 99.6%

          \[\leadsto \mathsf{copysign}\left(\log \left(\frac{1}{\color{blue}{-\left(x - \mathsf{hypot}\left(1, x\right)\right)}}\right), x\right) \]

      13. sub-neg 99.6%

          \[\leadsto \mathsf{copysign}\left(\log \left(\frac{1}{-\color{blue}{\left(x + \left(-\mathsf{hypot}\left(1, x\right)\right)\right)}}\right), x\right) \]

      14. +-commutative 99.6%

          \[\leadsto \mathsf{copysign}\left(\log \left(\frac{1}{-\color{blue}{\left(\left(-\mathsf{hypot}\left(1, x\right)\right) + x\right)}}\right), x\right) \]

      15. distribute-neg-in 99.6%

          \[\leadsto \mathsf{copysign}\left(\log \left(\frac{1}{\color{blue}{\left(-\left(-\mathsf{hypot}\left(1, x\right)\right)\right) + \left(-x\right)}}\right), x\right) \]

      16. remove-double-neg 99.6%

          \[\leadsto \mathsf{copysign}\left(\log \left(\frac{1}{\color{blue}{\mathsf{hypot}\left(1, x\right)} + \left(-x\right)}\right), x\right) \]

      17. sub-neg 99.6%

          \[\leadsto \mathsf{copysign}\left(\log \left(\frac{1}{\color{blue}{\mathsf{hypot}\left(1, x\right) - x}}\right), x\right) \]

   7. Simplified 99.6%

      \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(\frac{1}{\mathsf{hypot}\left(1, x\right) - x}\right)}, x\right) \]

   if -5.0000000000000001e-4 < (copysign.f64 (log.f64 (+.f64 (fabs.f64 x) (sqrt.f64 (+.f64 (*.f64 x x) 1)))) x) < -0.0

   1. Initial program 8.2%

      \[\mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{x \cdot x + 1}\right), x\right) \]
   2. Step-by-step derivation

      1. +-commutative 8.2%

         \[\leadsto \mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{\color{blue}{1 + x \cdot x}}\right), x\right) \]

      2. hypot-1-def 8.3%

         \[\leadsto \mathsf{copysign}\left(\log \left(\left|x\right| + \color{blue}{\mathsf{hypot}\left(1, x\right)}\right), x\right) \]

   3. Simplified 8.3%

      \[\leadsto \color{blue}{\mathsf{copysign}\left(\log \left(\left|x\right| + \mathsf{hypot}\left(1, x\right)\right), x\right)} \]

   4. Taylor expanded in x around 0 9.3%

      \[\leadsto \mathsf{copysign}\left(\color{blue}{0.5 \cdot \frac{{x}^{2}}{1 + \left|x\right|} + \log \left(1 + \left|x\right|\right)}, x\right) \]
   5. Step-by-step derivation

      1. fma-def 9.3%

         \[\leadsto \mathsf{copysign}\left(\color{blue}{\mathsf{fma}\left(0.5, \frac{{x}^{2}}{1 + \left|x\right|}, \log \left(1 + \left|x\right|\right)\right)}, x\right) \]

      2. unpow2 9.3%

         \[\leadsto \mathsf{copysign}\left(\mathsf{fma}\left(0.5, \frac{\color{blue}{x \cdot x}}{1 + \left|x\right|}, \log \left(1 + \left|x\right|\right)\right), x\right) \]

      3. log1p-def 99.8%

         \[\leadsto \mathsf{copysign}\left(\mathsf{fma}\left(0.5, \frac{x \cdot x}{1 + \left|x\right|}, \color{blue}{\mathsf{log1p}\left(\left|x\right|\right)}\right), x\right) \]

   6. Simplified 99.8%

      \[\leadsto \mathsf{copysign}\left(\color{blue}{\mathsf{fma}\left(0.5, \frac{x \cdot x}{1 + \left|x\right|}, \mathsf{log1p}\left(\left|x\right|\right)\right)}, x\right) \]
   7. Step-by-step derivation

      1. add-log-exp 8.4%

         \[\leadsto \mathsf{copysign}\left(\color{blue}{\log \left(e^{\mathsf{fma}\left(0.5, \frac{x \cdot x}{1 + \left|x\right|}, \mathsf{log1p}\left(\left|x\right|\right)\right)}\right)}, x\right) \]

      2. fma-udef 8.4%

         \[\leadsto \mathsf{copysign}\left(\log \left(e^{\color{blue}{0.5 \cdot \frac{x \cdot x}{1 + \left|x\right|} + \mathsf{log1p}\left(\left|x\right|\right)}}\right), x\right) \]

      3. exp-sum 8.3%

         \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(e^{0.5 \cdot \frac{x \cdot x}{1 + \left|x\right|}} \cdot e^{\mathsf{log1p}\left(\left|x\right|\right)}\right)}, x\right) \]

      4. add-sqr-sqrt 2.6%

         \[\leadsto \mathsf{copysign}\left(\log \left(e^{0.5 \cdot \frac{x \cdot x}{1 + \left|\color{blue}{\sqrt{x} \cdot \sqrt{x}}\right|}} \cdot e^{\mathsf{log1p}\left(\left|x\right|\right)}\right), x\right) \]

      5. fabs-sqr 2.6%

         \[\leadsto \mathsf{copysign}\left(\log \left(e^{0.5 \cdot \frac{x \cdot x}{1 + \color{blue}{\sqrt{x} \cdot \sqrt{x}}}} \cdot e^{\mathsf{log1p}\left(\left|x\right|\right)}\right), x\right) \]

      6. add-sqr-sqrt 8.0%

         \[\leadsto \mathsf{copysign}\left(\log \left(e^{0.5 \cdot \frac{x \cdot x}{1 + \color{blue}{x}}} \cdot e^{\mathsf{log1p}\left(\left|x\right|\right)}\right), x\right) \]

      7. log1p-udef 8.0%

         \[\leadsto \mathsf{copysign}\left(\log \left(e^{0.5 \cdot \frac{x \cdot x}{1 + x}} \cdot e^{\color{blue}{\log \left(1 + \left|x\right|\right)}}\right), x\right) \]

      8. add-exp-log 8.0%

         \[\leadsto \mathsf{copysign}\left(\log \left(e^{0.5 \cdot \frac{x \cdot x}{1 + x}} \cdot \color{blue}{\left(1 + \left|x\right|\right)}\right), x\right) \]

      9. add-sqr-sqrt 2.6%

         \[\leadsto \mathsf{copysign}\left(\log \left(e^{0.5 \cdot \frac{x \cdot x}{1 + x}} \cdot \left(1 + \left|\color{blue}{\sqrt{x} \cdot \sqrt{x}}\right|\right)\right), x\right) \]

      10. fabs-sqr 2.6%

          \[\leadsto \mathsf{copysign}\left(\log \left(e^{0.5 \cdot \frac{x \cdot x}{1 + x}} \cdot \left(1 + \color{blue}{\sqrt{x} \cdot \sqrt{x}}\right)\right), x\right) \]

      11. add-sqr-sqrt 8.4%

          \[\leadsto \mathsf{copysign}\left(\log \left(e^{0.5 \cdot \frac{x \cdot x}{1 + x}} \cdot \left(1 + \color{blue}{x}\right)\right), x\right) \]

      12. add-exp-log 8.4%

          \[\leadsto \mathsf{copysign}\left(\log \left(e^{0.5 \cdot \frac{x \cdot x}{1 + x}} \cdot \color{blue}{e^{\log \left(1 + x\right)}}\right), x\right) \]

      13. log1p-udef 8.4%

          \[\leadsto \mathsf{copysign}\left(\log \left(e^{0.5 \cdot \frac{x \cdot x}{1 + x}} \cdot e^{\color{blue}{\mathsf{log1p}\left(x\right)}}\right), x\right) \]

      14. exp-sum 8.4%

          \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(e^{0.5 \cdot \frac{x \cdot x}{1 + x} + \mathsf{log1p}\left(x\right)}\right)}, x\right) \]

   8. Applied egg-rr 99.8%

      \[\leadsto \mathsf{copysign}\left(\color{blue}{0.5 \cdot \frac{x}{\frac{x + 1}{x}} + \mathsf{log1p}\left(x\right)}, x\right) \]

   if -0.0 < (copysign.f64 (log.f64 (+.f64 (fabs.f64 x) (sqrt.f64 (+.f64 (*.f64 x x) 1)))) x)

   1. Initial program 52.3%

      \[\mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{x \cdot x + 1}\right), x\right) \]
   2. Step-by-step derivation

      1. +-commutative 52.3%

         \[\leadsto \mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{\color{blue}{1 + x \cdot x}}\right), x\right) \]

      2. hypot-1-def 99.9%

         \[\leadsto \mathsf{copysign}\left(\log \left(\left|x\right| + \color{blue}{\mathsf{hypot}\left(1, x\right)}\right), x\right) \]

   3. Simplified 99.9%

      \[\leadsto \color{blue}{\mathsf{copysign}\left(\log \left(\left|x\right| + \mathsf{hypot}\left(1, x\right)\right), x\right)} \]
   4. Step-by-step derivation

      1. *-un-lft-identity 99.9%

         \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(1 \cdot \left(\left|x\right| + \mathsf{hypot}\left(1, x\right)\right)\right)}, x\right) \]

      2. log-prod 99.9%

         \[\leadsto \mathsf{copysign}\left(\color{blue}{\log 1 + \log \left(\left|x\right| + \mathsf{hypot}\left(1, x\right)\right)}, x\right) \]

      3. metadata-eval 99.9%

         \[\leadsto \mathsf{copysign}\left(\color{blue}{0} + \log \left(\left|x\right| + \mathsf{hypot}\left(1, x\right)\right), x\right) \]

      4. *-un-lft-identity 99.9%

         \[\leadsto \mathsf{copysign}\left(0 + \log \color{blue}{\left(1 \cdot \left(\left|x\right| + \mathsf{hypot}\left(1, x\right)\right)\right)}, x\right) \]

      5. *-un-lft-identity 99.9%

         \[\leadsto \mathsf{copysign}\left(0 + \log \color{blue}{\left(\left|x\right| + \mathsf{hypot}\left(1, x\right)\right)}, x\right) \]

      6. add-sqr-sqrt 99.9%

         \[\leadsto \mathsf{copysign}\left(0 + \log \left(\left|\color{blue}{\sqrt{x} \cdot \sqrt{x}}\right| + \mathsf{hypot}\left(1, x\right)\right), x\right) \]

      7. fabs-sqr 99.9%

         \[\leadsto \mathsf{copysign}\left(0 + \log \left(\color{blue}{\sqrt{x} \cdot \sqrt{x}} + \mathsf{hypot}\left(1, x\right)\right), x\right) \]

      8. add-sqr-sqrt 99.9%

         \[\leadsto \mathsf{copysign}\left(0 + \log \left(\color{blue}{x} + \mathsf{hypot}\left(1, x\right)\right), x\right) \]

   5. Applied egg-rr 99.9%

      \[\leadsto \mathsf{copysign}\left(\color{blue}{0 + \log \left(x + \mathsf{hypot}\left(1, x\right)\right)}, x\right) \]
   6. Step-by-step derivation

      1. +-lft-identity 99.9%

         \[\leadsto \mathsf{copysign}\left(\color{blue}{\log \left(x + \mathsf{hypot}\left(1, x\right)\right)}, x\right) \]

   7. Simplified 99.9%

      \[\leadsto \mathsf{copysign}\left(\color{blue}{\log \left(x + \mathsf{hypot}\left(1, x\right)\right)}, x\right) \]
3. Recombined 3 regimes into one program.

4. Final simplification 99.8%

   \[\leadsto \begin{array}{l} \mathbf{if}\;\mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{x \cdot x + 1}\right), x\right) \leq -0.0005:\\ \;\;\;\;\mathsf{copysign}\left(\log \left(\frac{1}{\mathsf{hypot}\left(1, x\right) - x}\right), x\right)\\ \mathbf{elif}\;\mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{x \cdot x + 1}\right), x\right) \leq 0:\\ \;\;\;\;\mathsf{copysign}\left(0.5 \cdot \frac{x}{\frac{x + 1}{x}} + \mathsf{log1p}\left(x\right), x\right)\\ \mathbf{else}:\\ \;\;\;\;\mathsf{copysign}\left(\log \left(x + \mathsf{hypot}\left(1, x\right)\right), x\right)\\ \end{array} \]

Alternative 4: 99.5% accurate, 1.3× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;x \leq -0.72:\\ \;\;\;\;\mathsf{copysign}\left(\log \left(\frac{-0.5}{x}\right), x\right)\\ \mathbf{elif}\;x \leq 0.000135:\\ \;\;\;\;\mathsf{copysign}\left(0.5 \cdot \frac{x}{\frac{x + 1}{x}} + \mathsf{log1p}\left(x\right), x\right)\\ \mathbf{else}:\\ \;\;\;\;\mathsf{copysign}\left(\log \left(x + \mathsf{hypot}\left(1, x\right)\right), x\right)\\ \end{array} \end{array} \]

FPCore

(FPCore (x)
 :precision binary64
 (if (<= x -0.72)
   (copysign (log (/ -0.5 x)) x)
   (if (<= x 0.000135)
     (copysign (+ (* 0.5 (/ x (/ (+ x 1.0) x))) (log1p x)) x)
     (copysign (log (+ x (hypot 1.0 x))) x))))

C

double code(double x) {
	double tmp;
	if (x <= -0.72) {
		tmp = copysign(log((-0.5 / x)), x);
	} else if (x <= 0.000135) {
		tmp = copysign(((0.5 * (x / ((x + 1.0) / x))) + log1p(x)), x);
	} else {
		tmp = copysign(log((x + hypot(1.0, x))), x);
	}
	return tmp;
}

Java

public static double code(double x) {
	double tmp;
	if (x <= -0.72) {
		tmp = Math.copySign(Math.log((-0.5 / x)), x);
	} else if (x <= 0.000135) {
		tmp = Math.copySign(((0.5 * (x / ((x + 1.0) / x))) + Math.log1p(x)), x);
	} else {
		tmp = Math.copySign(Math.log((x + Math.hypot(1.0, x))), x);
	}
	return tmp;
}

Python

def code(x):
	tmp = 0
	if x <= -0.72:
		tmp = math.copysign(math.log((-0.5 / x)), x)
	elif x <= 0.000135:
		tmp = math.copysign(((0.5 * (x / ((x + 1.0) / x))) + math.log1p(x)), x)
	else:
		tmp = math.copysign(math.log((x + math.hypot(1.0, x))), x)
	return tmp

Julia

function code(x)
	tmp = 0.0
	if (x <= -0.72)
		tmp = copysign(log(Float64(-0.5 / x)), x);
	elseif (x <= 0.000135)
		tmp = copysign(Float64(Float64(0.5 * Float64(x / Float64(Float64(x + 1.0) / x))) + log1p(x)), x);
	else
		tmp = copysign(log(Float64(x + hypot(1.0, x))), x);
	end
	return tmp
end

Wolfram

code[x_] := If[LessEqual[x, -0.72], N[With[{TMP1 = Abs[N[Log[N[(-0.5 / x), $MachinePrecision]], $MachinePrecision]], TMP2 = Sign[x]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision], If[LessEqual[x, 0.000135], N[With[{TMP1 = Abs[N[(N[(0.5 * N[(x / N[(N[(x + 1.0), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[Log[1 + x], $MachinePrecision]), $MachinePrecision]], TMP2 = Sign[x]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision], N[With[{TMP1 = Abs[N[Log[N[(x + N[Sqrt[1.0 ^ 2 + x ^ 2], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], TMP2 = Sign[x]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]]]

Derivation

1. Split input into 3 regimes

2. if x < -0.71999999999999997

   1. Initial program 50.8%

      \[\mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{x \cdot x + 1}\right), x\right) \]
   2. Step-by-step derivation

      1. +-commutative 50.8%

         \[\leadsto \mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{\color{blue}{1 + x \cdot x}}\right), x\right) \]

      2. hypot-1-def 100.0%

         \[\leadsto \mathsf{copysign}\left(\log \left(\left|x\right| + \color{blue}{\mathsf{hypot}\left(1, x\right)}\right), x\right) \]

   3. Simplified 100.0%

      \[\leadsto \color{blue}{\mathsf{copysign}\left(\log \left(\left|x\right| + \mathsf{hypot}\left(1, x\right)\right), x\right)} \]

   4. Taylor expanded in x around -inf 100.0%

      \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(\left(\left|x\right| + -1 \cdot x\right) - 0.5 \cdot \frac{1}{x}\right)}, x\right) \]
   5. Step-by-step derivation

      1. associate--l+ 100.0%

         \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(\left|x\right| + \left(-1 \cdot x - 0.5 \cdot \frac{1}{x}\right)\right)}, x\right) \]

      2. unpow1 100.0%

         \[\leadsto \mathsf{copysign}\left(\log \left(\left|\color{blue}{{x}^{1}}\right| + \left(-1 \cdot x - 0.5 \cdot \frac{1}{x}\right)\right), x\right) \]

      3. sqr-pow 0.0%

         \[\leadsto \mathsf{copysign}\left(\log \left(\left|\color{blue}{{x}^{\left(\frac{1}{2}\right)} \cdot {x}^{\left(\frac{1}{2}\right)}}\right| + \left(-1 \cdot x - 0.5 \cdot \frac{1}{x}\right)\right), x\right) \]

      4. fabs-sqr 0.0%

         \[\leadsto \mathsf{copysign}\left(\log \left(\color{blue}{{x}^{\left(\frac{1}{2}\right)} \cdot {x}^{\left(\frac{1}{2}\right)}} + \left(-1 \cdot x - 0.5 \cdot \frac{1}{x}\right)\right), x\right) \]

      5. sqr-pow 4.1%

         \[\leadsto \mathsf{copysign}\left(\log \left(\color{blue}{{x}^{1}} + \left(-1 \cdot x - 0.5 \cdot \frac{1}{x}\right)\right), x\right) \]

      6. unpow1 4.1%

         \[\leadsto \mathsf{copysign}\left(\log \left(\color{blue}{x} + \left(-1 \cdot x - 0.5 \cdot \frac{1}{x}\right)\right), x\right) \]

      7. associate-+r- 99.6%

         \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(\left(x + -1 \cdot x\right) - 0.5 \cdot \frac{1}{x}\right)}, x\right) \]

      8. mul-1-neg 99.6%

         \[\leadsto \mathsf{copysign}\left(\log \left(\left(x + \color{blue}{\left(-x\right)}\right) - 0.5 \cdot \frac{1}{x}\right), x\right) \]

      9. sub-neg 99.6%

         \[\leadsto \mathsf{copysign}\left(\log \left(\color{blue}{\left(x - x\right)} - 0.5 \cdot \frac{1}{x}\right), x\right) \]

      10. +-inverses 99.6%

          \[\leadsto \mathsf{copysign}\left(\log \left(\color{blue}{0} - 0.5 \cdot \frac{1}{x}\right), x\right) \]

      11. neg-sub0 99.6%

          \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(-0.5 \cdot \frac{1}{x}\right)}, x\right) \]

      12. associate-*r/ 99.6%

          \[\leadsto \mathsf{copysign}\left(\log \left(-\color{blue}{\frac{0.5 \cdot 1}{x}}\right), x\right) \]

      13. metadata-eval 99.6%

          \[\leadsto \mathsf{copysign}\left(\log \left(-\frac{\color{blue}{0.5}}{x}\right), x\right) \]

      14. distribute-neg-frac 99.6%

          \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(\frac{-0.5}{x}\right)}, x\right) \]

      15. metadata-eval 99.6%

          \[\leadsto \mathsf{copysign}\left(\log \left(\frac{\color{blue}{-0.5}}{x}\right), x\right) \]

   6. Simplified 99.6%

      \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(\frac{-0.5}{x}\right)}, x\right) \]

   if -0.71999999999999997 < x < 1.35000000000000002e-4

   1. Initial program 9.4%

      \[\mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{x \cdot x + 1}\right), x\right) \]
   2. Step-by-step derivation

      1. +-commutative 9.4%

         \[\leadsto \mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{\color{blue}{1 + x \cdot x}}\right), x\right) \]

      2. hypot-1-def 9.5%

         \[\leadsto \mathsf{copysign}\left(\log \left(\left|x\right| + \color{blue}{\mathsf{hypot}\left(1, x\right)}\right), x\right) \]

   3. Simplified 9.5%

      \[\leadsto \color{blue}{\mathsf{copysign}\left(\log \left(\left|x\right| + \mathsf{hypot}\left(1, x\right)\right), x\right)} \]

   4. Taylor expanded in x around 0 10.2%

      \[\leadsto \mathsf{copysign}\left(\color{blue}{0.5 \cdot \frac{{x}^{2}}{1 + \left|x\right|} + \log \left(1 + \left|x\right|\right)}, x\right) \]
   5. Step-by-step derivation

      1. fma-def 10.2%

         \[\leadsto \mathsf{copysign}\left(\color{blue}{\mathsf{fma}\left(0.5, \frac{{x}^{2}}{1 + \left|x\right|}, \log \left(1 + \left|x\right|\right)\right)}, x\right) \]

      2. unpow2 10.2%

         \[\leadsto \mathsf{copysign}\left(\mathsf{fma}\left(0.5, \frac{\color{blue}{x \cdot x}}{1 + \left|x\right|}, \log \left(1 + \left|x\right|\right)\right), x\right) \]

      3. log1p-def 99.3%

         \[\leadsto \mathsf{copysign}\left(\mathsf{fma}\left(0.5, \frac{x \cdot x}{1 + \left|x\right|}, \color{blue}{\mathsf{log1p}\left(\left|x\right|\right)}\right), x\right) \]

   6. Simplified 99.3%

      \[\leadsto \mathsf{copysign}\left(\color{blue}{\mathsf{fma}\left(0.5, \frac{x \cdot x}{1 + \left|x\right|}, \mathsf{log1p}\left(\left|x\right|\right)\right)}, x\right) \]
   7. Step-by-step derivation

      1. add-log-exp 9.3%

         \[\leadsto \mathsf{copysign}\left(\color{blue}{\log \left(e^{\mathsf{fma}\left(0.5, \frac{x \cdot x}{1 + \left|x\right|}, \mathsf{log1p}\left(\left|x\right|\right)\right)}\right)}, x\right) \]

      2. fma-udef 9.3%

         \[\leadsto \mathsf{copysign}\left(\log \left(e^{\color{blue}{0.5 \cdot \frac{x \cdot x}{1 + \left|x\right|} + \mathsf{log1p}\left(\left|x\right|\right)}}\right), x\right) \]

      3. exp-sum 9.2%

         \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(e^{0.5 \cdot \frac{x \cdot x}{1 + \left|x\right|}} \cdot e^{\mathsf{log1p}\left(\left|x\right|\right)}\right)}, x\right) \]

      4. add-sqr-sqrt 2.5%

         \[\leadsto \mathsf{copysign}\left(\log \left(e^{0.5 \cdot \frac{x \cdot x}{1 + \left|\color{blue}{\sqrt{x} \cdot \sqrt{x}}\right|}} \cdot e^{\mathsf{log1p}\left(\left|x\right|\right)}\right), x\right) \]

      5. fabs-sqr 2.5%

         \[\leadsto \mathsf{copysign}\left(\log \left(e^{0.5 \cdot \frac{x \cdot x}{1 + \color{blue}{\sqrt{x} \cdot \sqrt{x}}}} \cdot e^{\mathsf{log1p}\left(\left|x\right|\right)}\right), x\right) \]

      6. add-sqr-sqrt 8.6%

         \[\leadsto \mathsf{copysign}\left(\log \left(e^{0.5 \cdot \frac{x \cdot x}{1 + \color{blue}{x}}} \cdot e^{\mathsf{log1p}\left(\left|x\right|\right)}\right), x\right) \]

      7. log1p-udef 8.6%

         \[\leadsto \mathsf{copysign}\left(\log \left(e^{0.5 \cdot \frac{x \cdot x}{1 + x}} \cdot e^{\color{blue}{\log \left(1 + \left|x\right|\right)}}\right), x\right) \]

      8. add-exp-log 8.6%

         \[\leadsto \mathsf{copysign}\left(\log \left(e^{0.5 \cdot \frac{x \cdot x}{1 + x}} \cdot \color{blue}{\left(1 + \left|x\right|\right)}\right), x\right) \]

      9. add-sqr-sqrt 2.5%

         \[\leadsto \mathsf{copysign}\left(\log \left(e^{0.5 \cdot \frac{x \cdot x}{1 + x}} \cdot \left(1 + \left|\color{blue}{\sqrt{x} \cdot \sqrt{x}}\right|\right)\right), x\right) \]

      10. fabs-sqr 2.5%

          \[\leadsto \mathsf{copysign}\left(\log \left(e^{0.5 \cdot \frac{x \cdot x}{1 + x}} \cdot \left(1 + \color{blue}{\sqrt{x} \cdot \sqrt{x}}\right)\right), x\right) \]

      11. add-sqr-sqrt 9.4%

          \[\leadsto \mathsf{copysign}\left(\log \left(e^{0.5 \cdot \frac{x \cdot x}{1 + x}} \cdot \left(1 + \color{blue}{x}\right)\right), x\right) \]

      12. add-exp-log 9.4%

          \[\leadsto \mathsf{copysign}\left(\log \left(e^{0.5 \cdot \frac{x \cdot x}{1 + x}} \cdot \color{blue}{e^{\log \left(1 + x\right)}}\right), x\right) \]

      13. log1p-udef 9.4%

          \[\leadsto \mathsf{copysign}\left(\log \left(e^{0.5 \cdot \frac{x \cdot x}{1 + x}} \cdot e^{\color{blue}{\mathsf{log1p}\left(x\right)}}\right), x\right) \]

      14. exp-sum 9.4%

          \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(e^{0.5 \cdot \frac{x \cdot x}{1 + x} + \mathsf{log1p}\left(x\right)}\right)}, x\right) \]

   8. Applied egg-rr 99.3%

      \[\leadsto \mathsf{copysign}\left(\color{blue}{0.5 \cdot \frac{x}{\frac{x + 1}{x}} + \mathsf{log1p}\left(x\right)}, x\right) \]

   if 1.35000000000000002e-4 < x

   1. Initial program 52.3%

      \[\mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{x \cdot x + 1}\right), x\right) \]
   2. Step-by-step derivation

      1. +-commutative 52.3%

         \[\leadsto \mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{\color{blue}{1 + x \cdot x}}\right), x\right) \]

      2. hypot-1-def 99.9%

         \[\leadsto \mathsf{copysign}\left(\log \left(\left|x\right| + \color{blue}{\mathsf{hypot}\left(1, x\right)}\right), x\right) \]

   3. Simplified 99.9%

      \[\leadsto \color{blue}{\mathsf{copysign}\left(\log \left(\left|x\right| + \mathsf{hypot}\left(1, x\right)\right), x\right)} \]
   4. Step-by-step derivation

      1. *-un-lft-identity 99.9%

         \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(1 \cdot \left(\left|x\right| + \mathsf{hypot}\left(1, x\right)\right)\right)}, x\right) \]

      2. log-prod 99.9%

         \[\leadsto \mathsf{copysign}\left(\color{blue}{\log 1 + \log \left(\left|x\right| + \mathsf{hypot}\left(1, x\right)\right)}, x\right) \]

      3. metadata-eval 99.9%

         \[\leadsto \mathsf{copysign}\left(\color{blue}{0} + \log \left(\left|x\right| + \mathsf{hypot}\left(1, x\right)\right), x\right) \]

      4. *-un-lft-identity 99.9%

         \[\leadsto \mathsf{copysign}\left(0 + \log \color{blue}{\left(1 \cdot \left(\left|x\right| + \mathsf{hypot}\left(1, x\right)\right)\right)}, x\right) \]

      5. *-un-lft-identity 99.9%

         \[\leadsto \mathsf{copysign}\left(0 + \log \color{blue}{\left(\left|x\right| + \mathsf{hypot}\left(1, x\right)\right)}, x\right) \]

      6. add-sqr-sqrt 99.9%

         \[\leadsto \mathsf{copysign}\left(0 + \log \left(\left|\color{blue}{\sqrt{x} \cdot \sqrt{x}}\right| + \mathsf{hypot}\left(1, x\right)\right), x\right) \]

      7. fabs-sqr 99.9%

         \[\leadsto \mathsf{copysign}\left(0 + \log \left(\color{blue}{\sqrt{x} \cdot \sqrt{x}} + \mathsf{hypot}\left(1, x\right)\right), x\right) \]

      8. add-sqr-sqrt 99.9%

         \[\leadsto \mathsf{copysign}\left(0 + \log \left(\color{blue}{x} + \mathsf{hypot}\left(1, x\right)\right), x\right) \]

   5. Applied egg-rr 99.9%

      \[\leadsto \mathsf{copysign}\left(\color{blue}{0 + \log \left(x + \mathsf{hypot}\left(1, x\right)\right)}, x\right) \]
   6. Step-by-step derivation

      1. +-lft-identity 99.9%

         \[\leadsto \mathsf{copysign}\left(\color{blue}{\log \left(x + \mathsf{hypot}\left(1, x\right)\right)}, x\right) \]

   7. Simplified 99.9%

      \[\leadsto \mathsf{copysign}\left(\color{blue}{\log \left(x + \mathsf{hypot}\left(1, x\right)\right)}, x\right) \]
3. Recombined 3 regimes into one program.

4. Final simplification 99.5%

   \[\leadsto \begin{array}{l} \mathbf{if}\;x \leq -0.72:\\ \;\;\;\;\mathsf{copysign}\left(\log \left(\frac{-0.5}{x}\right), x\right)\\ \mathbf{elif}\;x \leq 0.000135:\\ \;\;\;\;\mathsf{copysign}\left(0.5 \cdot \frac{x}{\frac{x + 1}{x}} + \mathsf{log1p}\left(x\right), x\right)\\ \mathbf{else}:\\ \;\;\;\;\mathsf{copysign}\left(\log \left(x + \mathsf{hypot}\left(1, x\right)\right), x\right)\\ \end{array} \]

Alternative 5: 99.3% accurate, 1.9× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;x \leq -0.72:\\ \;\;\;\;\mathsf{copysign}\left(\log \left(\frac{-0.5}{x}\right), x\right)\\ \mathbf{elif}\;x \leq 0.91:\\ \;\;\;\;\mathsf{copysign}\left(0.5 \cdot \frac{x}{\frac{x + 1}{x}} + \mathsf{log1p}\left(x\right), x\right)\\ \mathbf{else}:\\ \;\;\;\;\mathsf{copysign}\left(\log \left(x + \left(x + \frac{0.5}{x}\right)\right), x\right)\\ \end{array} \end{array} \]

FPCore

(FPCore (x)
 :precision binary64
 (if (<= x -0.72)
   (copysign (log (/ -0.5 x)) x)
   (if (<= x 0.91)
     (copysign (+ (* 0.5 (/ x (/ (+ x 1.0) x))) (log1p x)) x)
     (copysign (log (+ x (+ x (/ 0.5 x)))) x))))

C

double code(double x) {
	double tmp;
	if (x <= -0.72) {
		tmp = copysign(log((-0.5 / x)), x);
	} else if (x <= 0.91) {
		tmp = copysign(((0.5 * (x / ((x + 1.0) / x))) + log1p(x)), x);
	} else {
		tmp = copysign(log((x + (x + (0.5 / x)))), x);
	}
	return tmp;
}

Java

public static double code(double x) {
	double tmp;
	if (x <= -0.72) {
		tmp = Math.copySign(Math.log((-0.5 / x)), x);
	} else if (x <= 0.91) {
		tmp = Math.copySign(((0.5 * (x / ((x + 1.0) / x))) + Math.log1p(x)), x);
	} else {
		tmp = Math.copySign(Math.log((x + (x + (0.5 / x)))), x);
	}
	return tmp;
}

Python

def code(x):
	tmp = 0
	if x <= -0.72:
		tmp = math.copysign(math.log((-0.5 / x)), x)
	elif x <= 0.91:
		tmp = math.copysign(((0.5 * (x / ((x + 1.0) / x))) + math.log1p(x)), x)
	else:
		tmp = math.copysign(math.log((x + (x + (0.5 / x)))), x)
	return tmp

Julia

function code(x)
	tmp = 0.0
	if (x <= -0.72)
		tmp = copysign(log(Float64(-0.5 / x)), x);
	elseif (x <= 0.91)
		tmp = copysign(Float64(Float64(0.5 * Float64(x / Float64(Float64(x + 1.0) / x))) + log1p(x)), x);
	else
		tmp = copysign(log(Float64(x + Float64(x + Float64(0.5 / x)))), x);
	end
	return tmp
end

Wolfram

code[x_] := If[LessEqual[x, -0.72], N[With[{TMP1 = Abs[N[Log[N[(-0.5 / x), $MachinePrecision]], $MachinePrecision]], TMP2 = Sign[x]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision], If[LessEqual[x, 0.91], N[With[{TMP1 = Abs[N[(N[(0.5 * N[(x / N[(N[(x + 1.0), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[Log[1 + x], $MachinePrecision]), $MachinePrecision]], TMP2 = Sign[x]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision], N[With[{TMP1 = Abs[N[Log[N[(x + N[(x + N[(0.5 / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], TMP2 = Sign[x]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]]]

Derivation

1. Split input into 3 regimes

2. if x < -0.71999999999999997

   1. Initial program 50.8%

      \[\mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{x \cdot x + 1}\right), x\right) \]
   2. Step-by-step derivation

      1. +-commutative 50.8%

         \[\leadsto \mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{\color{blue}{1 + x \cdot x}}\right), x\right) \]

      2. hypot-1-def 100.0%

         \[\leadsto \mathsf{copysign}\left(\log \left(\left|x\right| + \color{blue}{\mathsf{hypot}\left(1, x\right)}\right), x\right) \]

   3. Simplified 100.0%

      \[\leadsto \color{blue}{\mathsf{copysign}\left(\log \left(\left|x\right| + \mathsf{hypot}\left(1, x\right)\right), x\right)} \]

   4. Taylor expanded in x around -inf 100.0%

      \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(\left(\left|x\right| + -1 \cdot x\right) - 0.5 \cdot \frac{1}{x}\right)}, x\right) \]
   5. Step-by-step derivation

      1. associate--l+ 100.0%

         \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(\left|x\right| + \left(-1 \cdot x - 0.5 \cdot \frac{1}{x}\right)\right)}, x\right) \]

      2. unpow1 100.0%

         \[\leadsto \mathsf{copysign}\left(\log \left(\left|\color{blue}{{x}^{1}}\right| + \left(-1 \cdot x - 0.5 \cdot \frac{1}{x}\right)\right), x\right) \]

      3. sqr-pow 0.0%

         \[\leadsto \mathsf{copysign}\left(\log \left(\left|\color{blue}{{x}^{\left(\frac{1}{2}\right)} \cdot {x}^{\left(\frac{1}{2}\right)}}\right| + \left(-1 \cdot x - 0.5 \cdot \frac{1}{x}\right)\right), x\right) \]

      4. fabs-sqr 0.0%

         \[\leadsto \mathsf{copysign}\left(\log \left(\color{blue}{{x}^{\left(\frac{1}{2}\right)} \cdot {x}^{\left(\frac{1}{2}\right)}} + \left(-1 \cdot x - 0.5 \cdot \frac{1}{x}\right)\right), x\right) \]

      5. sqr-pow 4.1%

         \[\leadsto \mathsf{copysign}\left(\log \left(\color{blue}{{x}^{1}} + \left(-1 \cdot x - 0.5 \cdot \frac{1}{x}\right)\right), x\right) \]

      6. unpow1 4.1%

         \[\leadsto \mathsf{copysign}\left(\log \left(\color{blue}{x} + \left(-1 \cdot x - 0.5 \cdot \frac{1}{x}\right)\right), x\right) \]

      7. associate-+r- 99.6%

         \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(\left(x + -1 \cdot x\right) - 0.5 \cdot \frac{1}{x}\right)}, x\right) \]

      8. mul-1-neg 99.6%

         \[\leadsto \mathsf{copysign}\left(\log \left(\left(x + \color{blue}{\left(-x\right)}\right) - 0.5 \cdot \frac{1}{x}\right), x\right) \]

      9. sub-neg 99.6%

         \[\leadsto \mathsf{copysign}\left(\log \left(\color{blue}{\left(x - x\right)} - 0.5 \cdot \frac{1}{x}\right), x\right) \]

      10. +-inverses 99.6%

          \[\leadsto \mathsf{copysign}\left(\log \left(\color{blue}{0} - 0.5 \cdot \frac{1}{x}\right), x\right) \]

      11. neg-sub0 99.6%

          \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(-0.5 \cdot \frac{1}{x}\right)}, x\right) \]

      12. associate-*r/ 99.6%

          \[\leadsto \mathsf{copysign}\left(\log \left(-\color{blue}{\frac{0.5 \cdot 1}{x}}\right), x\right) \]

      13. metadata-eval 99.6%

          \[\leadsto \mathsf{copysign}\left(\log \left(-\frac{\color{blue}{0.5}}{x}\right), x\right) \]

      14. distribute-neg-frac 99.6%

          \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(\frac{-0.5}{x}\right)}, x\right) \]

      15. metadata-eval 99.6%

          \[\leadsto \mathsf{copysign}\left(\log \left(\frac{\color{blue}{-0.5}}{x}\right), x\right) \]

   6. Simplified 99.6%

      \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(\frac{-0.5}{x}\right)}, x\right) \]

   if -0.71999999999999997 < x < 0.910000000000000031

   1. Initial program 10.7%

      \[\mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{x \cdot x + 1}\right), x\right) \]
   2. Step-by-step derivation

      1. +-commutative 10.7%

         \[\leadsto \mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{\color{blue}{1 + x \cdot x}}\right), x\right) \]

      2. hypot-1-def 10.8%

         \[\leadsto \mathsf{copysign}\left(\log \left(\left|x\right| + \color{blue}{\mathsf{hypot}\left(1, x\right)}\right), x\right) \]

   3. Simplified 10.8%

      \[\leadsto \color{blue}{\mathsf{copysign}\left(\log \left(\left|x\right| + \mathsf{hypot}\left(1, x\right)\right), x\right)} \]

   4. Taylor expanded in x around 0 10.5%

      \[\leadsto \mathsf{copysign}\left(\color{blue}{0.5 \cdot \frac{{x}^{2}}{1 + \left|x\right|} + \log \left(1 + \left|x\right|\right)}, x\right) \]
   5. Step-by-step derivation

      1. fma-def 10.5%

         \[\leadsto \mathsf{copysign}\left(\color{blue}{\mathsf{fma}\left(0.5, \frac{{x}^{2}}{1 + \left|x\right|}, \log \left(1 + \left|x\right|\right)\right)}, x\right) \]

      2. unpow2 10.5%

         \[\leadsto \mathsf{copysign}\left(\mathsf{fma}\left(0.5, \frac{\color{blue}{x \cdot x}}{1 + \left|x\right|}, \log \left(1 + \left|x\right|\right)\right), x\right) \]

      3. log1p-def 98.3%

         \[\leadsto \mathsf{copysign}\left(\mathsf{fma}\left(0.5, \frac{x \cdot x}{1 + \left|x\right|}, \color{blue}{\mathsf{log1p}\left(\left|x\right|\right)}\right), x\right) \]

   6. Simplified 98.3%

      \[\leadsto \mathsf{copysign}\left(\color{blue}{\mathsf{fma}\left(0.5, \frac{x \cdot x}{1 + \left|x\right|}, \mathsf{log1p}\left(\left|x\right|\right)\right)}, x\right) \]
   7. Step-by-step derivation

      1. add-log-exp 9.6%

         \[\leadsto \mathsf{copysign}\left(\color{blue}{\log \left(e^{\mathsf{fma}\left(0.5, \frac{x \cdot x}{1 + \left|x\right|}, \mathsf{log1p}\left(\left|x\right|\right)\right)}\right)}, x\right) \]

      2. fma-udef 9.6%

         \[\leadsto \mathsf{copysign}\left(\log \left(e^{\color{blue}{0.5 \cdot \frac{x \cdot x}{1 + \left|x\right|} + \mathsf{log1p}\left(\left|x\right|\right)}}\right), x\right) \]

      3. exp-sum 9.5%

         \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(e^{0.5 \cdot \frac{x \cdot x}{1 + \left|x\right|}} \cdot e^{\mathsf{log1p}\left(\left|x\right|\right)}\right)}, x\right) \]

      4. add-sqr-sqrt 2.9%

         \[\leadsto \mathsf{copysign}\left(\log \left(e^{0.5 \cdot \frac{x \cdot x}{1 + \left|\color{blue}{\sqrt{x} \cdot \sqrt{x}}\right|}} \cdot e^{\mathsf{log1p}\left(\left|x\right|\right)}\right), x\right) \]

      5. fabs-sqr 2.9%

         \[\leadsto \mathsf{copysign}\left(\log \left(e^{0.5 \cdot \frac{x \cdot x}{1 + \color{blue}{\sqrt{x} \cdot \sqrt{x}}}} \cdot e^{\mathsf{log1p}\left(\left|x\right|\right)}\right), x\right) \]

      6. add-sqr-sqrt 8.9%

         \[\leadsto \mathsf{copysign}\left(\log \left(e^{0.5 \cdot \frac{x \cdot x}{1 + \color{blue}{x}}} \cdot e^{\mathsf{log1p}\left(\left|x\right|\right)}\right), x\right) \]

      7. log1p-udef 8.9%

         \[\leadsto \mathsf{copysign}\left(\log \left(e^{0.5 \cdot \frac{x \cdot x}{1 + x}} \cdot e^{\color{blue}{\log \left(1 + \left|x\right|\right)}}\right), x\right) \]

      8. add-exp-log 8.9%

         \[\leadsto \mathsf{copysign}\left(\log \left(e^{0.5 \cdot \frac{x \cdot x}{1 + x}} \cdot \color{blue}{\left(1 + \left|x\right|\right)}\right), x\right) \]

      9. add-sqr-sqrt 2.9%

         \[\leadsto \mathsf{copysign}\left(\log \left(e^{0.5 \cdot \frac{x \cdot x}{1 + x}} \cdot \left(1 + \left|\color{blue}{\sqrt{x} \cdot \sqrt{x}}\right|\right)\right), x\right) \]

      10. fabs-sqr 2.9%

          \[\leadsto \mathsf{copysign}\left(\log \left(e^{0.5 \cdot \frac{x \cdot x}{1 + x}} \cdot \left(1 + \color{blue}{\sqrt{x} \cdot \sqrt{x}}\right)\right), x\right) \]

      11. add-sqr-sqrt 9.6%

          \[\leadsto \mathsf{copysign}\left(\log \left(e^{0.5 \cdot \frac{x \cdot x}{1 + x}} \cdot \left(1 + \color{blue}{x}\right)\right), x\right) \]

      12. add-exp-log 9.6%

          \[\leadsto \mathsf{copysign}\left(\log \left(e^{0.5 \cdot \frac{x \cdot x}{1 + x}} \cdot \color{blue}{e^{\log \left(1 + x\right)}}\right), x\right) \]

      13. log1p-udef 9.6%

          \[\leadsto \mathsf{copysign}\left(\log \left(e^{0.5 \cdot \frac{x \cdot x}{1 + x}} \cdot e^{\color{blue}{\mathsf{log1p}\left(x\right)}}\right), x\right) \]

      14. exp-sum 9.7%

          \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(e^{0.5 \cdot \frac{x \cdot x}{1 + x} + \mathsf{log1p}\left(x\right)}\right)}, x\right) \]

   8. Applied egg-rr 98.3%

      \[\leadsto \mathsf{copysign}\left(\color{blue}{0.5 \cdot \frac{x}{\frac{x + 1}{x}} + \mathsf{log1p}\left(x\right)}, x\right) \]

   if 0.910000000000000031 < x

   1. Initial program 50.7%

      \[\mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{x \cdot x + 1}\right), x\right) \]
   2. Step-by-step derivation

      1. +-commutative 50.7%

         \[\leadsto \mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{\color{blue}{1 + x \cdot x}}\right), x\right) \]

      2. hypot-1-def 100.0%

         \[\leadsto \mathsf{copysign}\left(\log \left(\left|x\right| + \color{blue}{\mathsf{hypot}\left(1, x\right)}\right), x\right) \]

   3. Simplified 100.0%

      \[\leadsto \color{blue}{\mathsf{copysign}\left(\log \left(\left|x\right| + \mathsf{hypot}\left(1, x\right)\right), x\right)} \]

   4. Taylor expanded in x around inf 100.0%

      \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(0.5 \cdot \frac{1}{x} + \left(\left|x\right| + x\right)\right)}, x\right) \]
   5. Step-by-step derivation

      1. associate-+r+ 100.0%

         \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(\left(0.5 \cdot \frac{1}{x} + \left|x\right|\right) + x\right)}, x\right) \]

      2. +-commutative 100.0%

         \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(x + \left(0.5 \cdot \frac{1}{x} + \left|x\right|\right)\right)}, x\right) \]

      3. +-commutative 100.0%

         \[\leadsto \mathsf{copysign}\left(\log \left(x + \color{blue}{\left(\left|x\right| + 0.5 \cdot \frac{1}{x}\right)}\right), x\right) \]

      4. unpow1 100.0%

         \[\leadsto \mathsf{copysign}\left(\log \left(x + \left(\left|\color{blue}{{x}^{1}}\right| + 0.5 \cdot \frac{1}{x}\right)\right), x\right) \]

      5. sqr-pow 100.0%

         \[\leadsto \mathsf{copysign}\left(\log \left(x + \left(\left|\color{blue}{{x}^{\left(\frac{1}{2}\right)} \cdot {x}^{\left(\frac{1}{2}\right)}}\right| + 0.5 \cdot \frac{1}{x}\right)\right), x\right) \]

      6. fabs-sqr 100.0%

         \[\leadsto \mathsf{copysign}\left(\log \left(x + \left(\color{blue}{{x}^{\left(\frac{1}{2}\right)} \cdot {x}^{\left(\frac{1}{2}\right)}} + 0.5 \cdot \frac{1}{x}\right)\right), x\right) \]

      7. sqr-pow 100.0%

         \[\leadsto \mathsf{copysign}\left(\log \left(x + \left(\color{blue}{{x}^{1}} + 0.5 \cdot \frac{1}{x}\right)\right), x\right) \]

      8. unpow1 100.0%

         \[\leadsto \mathsf{copysign}\left(\log \left(x + \left(\color{blue}{x} + 0.5 \cdot \frac{1}{x}\right)\right), x\right) \]

      9. associate-*r/ 100.0%

         \[\leadsto \mathsf{copysign}\left(\log \left(x + \left(x + \color{blue}{\frac{0.5 \cdot 1}{x}}\right)\right), x\right) \]

      10. metadata-eval 100.0%

          \[\leadsto \mathsf{copysign}\left(\log \left(x + \left(x + \frac{\color{blue}{0.5}}{x}\right)\right), x\right) \]

   6. Simplified 100.0%

      \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(x + \left(x + \frac{0.5}{x}\right)\right)}, x\right) \]
3. Recombined 3 regimes into one program.

4. Final simplification 99.0%

   \[\leadsto \begin{array}{l} \mathbf{if}\;x \leq -0.72:\\ \;\;\;\;\mathsf{copysign}\left(\log \left(\frac{-0.5}{x}\right), x\right)\\ \mathbf{elif}\;x \leq 0.91:\\ \;\;\;\;\mathsf{copysign}\left(0.5 \cdot \frac{x}{\frac{x + 1}{x}} + \mathsf{log1p}\left(x\right), x\right)\\ \mathbf{else}:\\ \;\;\;\;\mathsf{copysign}\left(\log \left(x + \left(x + \frac{0.5}{x}\right)\right), x\right)\\ \end{array} \]

Alternative 6: 99.1% accurate, 1.9× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;x \leq -1.25:\\ \;\;\;\;\mathsf{copysign}\left(\log \left(\frac{-0.5}{x}\right), x\right)\\ \mathbf{elif}\;x \leq 0.8:\\ \;\;\;\;\mathsf{copysign}\left(x, x\right)\\ \mathbf{else}:\\ \;\;\;\;\mathsf{copysign}\left(\log \left(x + \left(x + \frac{0.5}{x}\right)\right), x\right)\\ \end{array} \end{array} \]

FPCore

(FPCore (x)
 :precision binary64
 (if (<= x -1.25)
   (copysign (log (/ -0.5 x)) x)
   (if (<= x 0.8) (copysign x x) (copysign (log (+ x (+ x (/ 0.5 x)))) x))))

C

double code(double x) {
	double tmp;
	if (x <= -1.25) {
		tmp = copysign(log((-0.5 / x)), x);
	} else if (x <= 0.8) {
		tmp = copysign(x, x);
	} else {
		tmp = copysign(log((x + (x + (0.5 / x)))), x);
	}
	return tmp;
}

Java

public static double code(double x) {
	double tmp;
	if (x <= -1.25) {
		tmp = Math.copySign(Math.log((-0.5 / x)), x);
	} else if (x <= 0.8) {
		tmp = Math.copySign(x, x);
	} else {
		tmp = Math.copySign(Math.log((x + (x + (0.5 / x)))), x);
	}
	return tmp;
}

Python

def code(x):
	tmp = 0
	if x <= -1.25:
		tmp = math.copysign(math.log((-0.5 / x)), x)
	elif x <= 0.8:
		tmp = math.copysign(x, x)
	else:
		tmp = math.copysign(math.log((x + (x + (0.5 / x)))), x)
	return tmp

Julia

function code(x)
	tmp = 0.0
	if (x <= -1.25)
		tmp = copysign(log(Float64(-0.5 / x)), x);
	elseif (x <= 0.8)
		tmp = copysign(x, x);
	else
		tmp = copysign(log(Float64(x + Float64(x + Float64(0.5 / x)))), x);
	end
	return tmp
end

MATLAB

function tmp_2 = code(x)
	tmp = 0.0;
	if (x <= -1.25)
		tmp = sign(x) * abs(log((-0.5 / x)));
	elseif (x <= 0.8)
		tmp = sign(x) * abs(x);
	else
		tmp = sign(x) * abs(log((x + (x + (0.5 / x)))));
	end
	tmp_2 = tmp;
end

Wolfram

code[x_] := If[LessEqual[x, -1.25], N[With[{TMP1 = Abs[N[Log[N[(-0.5 / x), $MachinePrecision]], $MachinePrecision]], TMP2 = Sign[x]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision], If[LessEqual[x, 0.8], N[With[{TMP1 = Abs[x], TMP2 = Sign[x]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision], N[With[{TMP1 = Abs[N[Log[N[(x + N[(x + N[(0.5 / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], TMP2 = Sign[x]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]]]

Derivation

1. Split input into 3 regimes

2. if x < -1.25

   1. Initial program 50.8%

      \[\mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{x \cdot x + 1}\right), x\right) \]
   2. Step-by-step derivation

      1. +-commutative 50.8%

         \[\leadsto \mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{\color{blue}{1 + x \cdot x}}\right), x\right) \]

      2. hypot-1-def 100.0%

         \[\leadsto \mathsf{copysign}\left(\log \left(\left|x\right| + \color{blue}{\mathsf{hypot}\left(1, x\right)}\right), x\right) \]

   3. Simplified 100.0%

      \[\leadsto \color{blue}{\mathsf{copysign}\left(\log \left(\left|x\right| + \mathsf{hypot}\left(1, x\right)\right), x\right)} \]

   4. Taylor expanded in x around -inf 100.0%

      \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(\left(\left|x\right| + -1 \cdot x\right) - 0.5 \cdot \frac{1}{x}\right)}, x\right) \]
   5. Step-by-step derivation

      1. associate--l+ 100.0%

         \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(\left|x\right| + \left(-1 \cdot x - 0.5 \cdot \frac{1}{x}\right)\right)}, x\right) \]

      2. unpow1 100.0%

         \[\leadsto \mathsf{copysign}\left(\log \left(\left|\color{blue}{{x}^{1}}\right| + \left(-1 \cdot x - 0.5 \cdot \frac{1}{x}\right)\right), x\right) \]

      3. sqr-pow 0.0%

         \[\leadsto \mathsf{copysign}\left(\log \left(\left|\color{blue}{{x}^{\left(\frac{1}{2}\right)} \cdot {x}^{\left(\frac{1}{2}\right)}}\right| + \left(-1 \cdot x - 0.5 \cdot \frac{1}{x}\right)\right), x\right) \]

      4. fabs-sqr 0.0%

         \[\leadsto \mathsf{copysign}\left(\log \left(\color{blue}{{x}^{\left(\frac{1}{2}\right)} \cdot {x}^{\left(\frac{1}{2}\right)}} + \left(-1 \cdot x - 0.5 \cdot \frac{1}{x}\right)\right), x\right) \]

      5. sqr-pow 4.1%

         \[\leadsto \mathsf{copysign}\left(\log \left(\color{blue}{{x}^{1}} + \left(-1 \cdot x - 0.5 \cdot \frac{1}{x}\right)\right), x\right) \]

      6. unpow1 4.1%

         \[\leadsto \mathsf{copysign}\left(\log \left(\color{blue}{x} + \left(-1 \cdot x - 0.5 \cdot \frac{1}{x}\right)\right), x\right) \]

      7. associate-+r- 99.6%

         \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(\left(x + -1 \cdot x\right) - 0.5 \cdot \frac{1}{x}\right)}, x\right) \]

      8. mul-1-neg 99.6%

         \[\leadsto \mathsf{copysign}\left(\log \left(\left(x + \color{blue}{\left(-x\right)}\right) - 0.5 \cdot \frac{1}{x}\right), x\right) \]

      9. sub-neg 99.6%

         \[\leadsto \mathsf{copysign}\left(\log \left(\color{blue}{\left(x - x\right)} - 0.5 \cdot \frac{1}{x}\right), x\right) \]

      10. +-inverses 99.6%

          \[\leadsto \mathsf{copysign}\left(\log \left(\color{blue}{0} - 0.5 \cdot \frac{1}{x}\right), x\right) \]

      11. neg-sub0 99.6%

          \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(-0.5 \cdot \frac{1}{x}\right)}, x\right) \]

      12. associate-*r/ 99.6%

          \[\leadsto \mathsf{copysign}\left(\log \left(-\color{blue}{\frac{0.5 \cdot 1}{x}}\right), x\right) \]

      13. metadata-eval 99.6%

          \[\leadsto \mathsf{copysign}\left(\log \left(-\frac{\color{blue}{0.5}}{x}\right), x\right) \]

      14. distribute-neg-frac 99.6%

          \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(\frac{-0.5}{x}\right)}, x\right) \]

      15. metadata-eval 99.6%

          \[\leadsto \mathsf{copysign}\left(\log \left(\frac{\color{blue}{-0.5}}{x}\right), x\right) \]

   6. Simplified 99.6%

      \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(\frac{-0.5}{x}\right)}, x\right) \]

   if -1.25 < x < 0.80000000000000004

   1. Initial program 10.7%

      \[\mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{x \cdot x + 1}\right), x\right) \]
   2. Step-by-step derivation

      1. +-commutative 10.7%

         \[\leadsto \mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{\color{blue}{1 + x \cdot x}}\right), x\right) \]

      2. hypot-1-def 10.8%

         \[\leadsto \mathsf{copysign}\left(\log \left(\left|x\right| + \color{blue}{\mathsf{hypot}\left(1, x\right)}\right), x\right) \]

   3. Simplified 10.8%

      \[\leadsto \color{blue}{\mathsf{copysign}\left(\log \left(\left|x\right| + \mathsf{hypot}\left(1, x\right)\right), x\right)} \]

   4. Taylor expanded in x around 0 8.0%

      \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(1 + \left|x\right|\right)}, x\right) \]
   5. Step-by-step derivation

      1. unpow1 8.0%

         \[\leadsto \mathsf{copysign}\left(\log \left(1 + \left|\color{blue}{{x}^{1}}\right|\right), x\right) \]

      2. sqr-pow 2.8%

         \[\leadsto \mathsf{copysign}\left(\log \left(1 + \left|\color{blue}{{x}^{\left(\frac{1}{2}\right)} \cdot {x}^{\left(\frac{1}{2}\right)}}\right|\right), x\right) \]

      3. fabs-sqr 2.8%

         \[\leadsto \mathsf{copysign}\left(\log \left(1 + \color{blue}{{x}^{\left(\frac{1}{2}\right)} \cdot {x}^{\left(\frac{1}{2}\right)}}\right), x\right) \]

      4. sqr-pow 8.1%

         \[\leadsto \mathsf{copysign}\left(\log \left(1 + \color{blue}{{x}^{1}}\right), x\right) \]

      5. unpow1 8.1%

         \[\leadsto \mathsf{copysign}\left(\log \left(1 + \color{blue}{x}\right), x\right) \]

   6. Simplified 8.1%

      \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(1 + x\right)}, x\right) \]

   7. Taylor expanded in x around 0 97.4%

      \[\leadsto \mathsf{copysign}\left(\color{blue}{x}, x\right) \]

   if 0.80000000000000004 < x

   1. Initial program 50.7%

      \[\mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{x \cdot x + 1}\right), x\right) \]
   2. Step-by-step derivation

      1. +-commutative 50.7%

         \[\leadsto \mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{\color{blue}{1 + x \cdot x}}\right), x\right) \]

      2. hypot-1-def 100.0%

         \[\leadsto \mathsf{copysign}\left(\log \left(\left|x\right| + \color{blue}{\mathsf{hypot}\left(1, x\right)}\right), x\right) \]

   3. Simplified 100.0%

      \[\leadsto \color{blue}{\mathsf{copysign}\left(\log \left(\left|x\right| + \mathsf{hypot}\left(1, x\right)\right), x\right)} \]

   4. Taylor expanded in x around inf 100.0%

      \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(0.5 \cdot \frac{1}{x} + \left(\left|x\right| + x\right)\right)}, x\right) \]
   5. Step-by-step derivation

      1. associate-+r+ 100.0%

         \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(\left(0.5 \cdot \frac{1}{x} + \left|x\right|\right) + x\right)}, x\right) \]

      2. +-commutative 100.0%

         \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(x + \left(0.5 \cdot \frac{1}{x} + \left|x\right|\right)\right)}, x\right) \]

      3. +-commutative 100.0%

         \[\leadsto \mathsf{copysign}\left(\log \left(x + \color{blue}{\left(\left|x\right| + 0.5 \cdot \frac{1}{x}\right)}\right), x\right) \]

      4. unpow1 100.0%

         \[\leadsto \mathsf{copysign}\left(\log \left(x + \left(\left|\color{blue}{{x}^{1}}\right| + 0.5 \cdot \frac{1}{x}\right)\right), x\right) \]

      5. sqr-pow 100.0%

         \[\leadsto \mathsf{copysign}\left(\log \left(x + \left(\left|\color{blue}{{x}^{\left(\frac{1}{2}\right)} \cdot {x}^{\left(\frac{1}{2}\right)}}\right| + 0.5 \cdot \frac{1}{x}\right)\right), x\right) \]

      6. fabs-sqr 100.0%

         \[\leadsto \mathsf{copysign}\left(\log \left(x + \left(\color{blue}{{x}^{\left(\frac{1}{2}\right)} \cdot {x}^{\left(\frac{1}{2}\right)}} + 0.5 \cdot \frac{1}{x}\right)\right), x\right) \]

      7. sqr-pow 100.0%

         \[\leadsto \mathsf{copysign}\left(\log \left(x + \left(\color{blue}{{x}^{1}} + 0.5 \cdot \frac{1}{x}\right)\right), x\right) \]

      8. unpow1 100.0%

         \[\leadsto \mathsf{copysign}\left(\log \left(x + \left(\color{blue}{x} + 0.5 \cdot \frac{1}{x}\right)\right), x\right) \]

      9. associate-*r/ 100.0%

         \[\leadsto \mathsf{copysign}\left(\log \left(x + \left(x + \color{blue}{\frac{0.5 \cdot 1}{x}}\right)\right), x\right) \]

      10. metadata-eval 100.0%

          \[\leadsto \mathsf{copysign}\left(\log \left(x + \left(x + \frac{\color{blue}{0.5}}{x}\right)\right), x\right) \]

   6. Simplified 100.0%

      \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(x + \left(x + \frac{0.5}{x}\right)\right)}, x\right) \]
3. Recombined 3 regimes into one program.

4. Final simplification 98.6%

   \[\leadsto \begin{array}{l} \mathbf{if}\;x \leq -1.25:\\ \;\;\;\;\mathsf{copysign}\left(\log \left(\frac{-0.5}{x}\right), x\right)\\ \mathbf{elif}\;x \leq 0.8:\\ \;\;\;\;\mathsf{copysign}\left(x, x\right)\\ \mathbf{else}:\\ \;\;\;\;\mathsf{copysign}\left(\log \left(x + \left(x + \frac{0.5}{x}\right)\right), x\right)\\ \end{array} \]

Alternative 7: 82.6% accurate, 2.0× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;x \leq -3.2:\\ \;\;\;\;\mathsf{copysign}\left(\log \left(-x\right), x\right)\\ \mathbf{elif}\;x \leq 1.25:\\ \;\;\;\;\mathsf{copysign}\left(x, x\right)\\ \mathbf{else}:\\ \;\;\;\;\mathsf{copysign}\left(\log \left(x + x\right), x\right)\\ \end{array} \end{array} \]

FPCore

(FPCore (x)
 :precision binary64
 (if (<= x -3.2)
   (copysign (log (- x)) x)
   (if (<= x 1.25) (copysign x x) (copysign (log (+ x x)) x))))

C

double code(double x) {
	double tmp;
	if (x <= -3.2) {
		tmp = copysign(log(-x), x);
	} else if (x <= 1.25) {
		tmp = copysign(x, x);
	} else {
		tmp = copysign(log((x + x)), x);
	}
	return tmp;
}

Java

public static double code(double x) {
	double tmp;
	if (x <= -3.2) {
		tmp = Math.copySign(Math.log(-x), x);
	} else if (x <= 1.25) {
		tmp = Math.copySign(x, x);
	} else {
		tmp = Math.copySign(Math.log((x + x)), x);
	}
	return tmp;
}

Python

def code(x):
	tmp = 0
	if x <= -3.2:
		tmp = math.copysign(math.log(-x), x)
	elif x <= 1.25:
		tmp = math.copysign(x, x)
	else:
		tmp = math.copysign(math.log((x + x)), x)
	return tmp

Julia

function code(x)
	tmp = 0.0
	if (x <= -3.2)
		tmp = copysign(log(Float64(-x)), x);
	elseif (x <= 1.25)
		tmp = copysign(x, x);
	else
		tmp = copysign(log(Float64(x + x)), x);
	end
	return tmp
end

MATLAB

function tmp_2 = code(x)
	tmp = 0.0;
	if (x <= -3.2)
		tmp = sign(x) * abs(log(-x));
	elseif (x <= 1.25)
		tmp = sign(x) * abs(x);
	else
		tmp = sign(x) * abs(log((x + x)));
	end
	tmp_2 = tmp;
end

Wolfram

code[x_] := If[LessEqual[x, -3.2], N[With[{TMP1 = Abs[N[Log[(-x)], $MachinePrecision]], TMP2 = Sign[x]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision], If[LessEqual[x, 1.25], N[With[{TMP1 = Abs[x], TMP2 = Sign[x]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision], N[With[{TMP1 = Abs[N[Log[N[(x + x), $MachinePrecision]], $MachinePrecision]], TMP2 = Sign[x]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]]]

Derivation

1. Split input into 3 regimes

2. if x < -3.2000000000000002

   1. Initial program 50.8%

      \[\mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{x \cdot x + 1}\right), x\right) \]
   2. Step-by-step derivation

      1. +-commutative 50.8%

         \[\leadsto \mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{\color{blue}{1 + x \cdot x}}\right), x\right) \]

      2. hypot-1-def 100.0%

         \[\leadsto \mathsf{copysign}\left(\log \left(\left|x\right| + \color{blue}{\mathsf{hypot}\left(1, x\right)}\right), x\right) \]

   3. Simplified 100.0%

      \[\leadsto \color{blue}{\mathsf{copysign}\left(\log \left(\left|x\right| + \mathsf{hypot}\left(1, x\right)\right), x\right)} \]

   4. Taylor expanded in x around -inf 31.5%

      \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(-1 \cdot x\right)}, x\right) \]
   5. Step-by-step derivation

      1. mul-1-neg 31.5%

         \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(-x\right)}, x\right) \]

   6. Simplified 31.5%

      \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(-x\right)}, x\right) \]

   if -3.2000000000000002 < x < 1.25

   1. Initial program 10.7%

      \[\mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{x \cdot x + 1}\right), x\right) \]
   2. Step-by-step derivation

      1. +-commutative 10.7%

         \[\leadsto \mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{\color{blue}{1 + x \cdot x}}\right), x\right) \]

      2. hypot-1-def 10.8%

         \[\leadsto \mathsf{copysign}\left(\log \left(\left|x\right| + \color{blue}{\mathsf{hypot}\left(1, x\right)}\right), x\right) \]

   3. Simplified 10.8%

      \[\leadsto \color{blue}{\mathsf{copysign}\left(\log \left(\left|x\right| + \mathsf{hypot}\left(1, x\right)\right), x\right)} \]

   4. Taylor expanded in x around 0 8.0%

      \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(1 + \left|x\right|\right)}, x\right) \]
   5. Step-by-step derivation

      1. unpow1 8.0%

         \[\leadsto \mathsf{copysign}\left(\log \left(1 + \left|\color{blue}{{x}^{1}}\right|\right), x\right) \]

      2. sqr-pow 2.8%

         \[\leadsto \mathsf{copysign}\left(\log \left(1 + \left|\color{blue}{{x}^{\left(\frac{1}{2}\right)} \cdot {x}^{\left(\frac{1}{2}\right)}}\right|\right), x\right) \]

      3. fabs-sqr 2.8%

         \[\leadsto \mathsf{copysign}\left(\log \left(1 + \color{blue}{{x}^{\left(\frac{1}{2}\right)} \cdot {x}^{\left(\frac{1}{2}\right)}}\right), x\right) \]

      4. sqr-pow 8.1%

         \[\leadsto \mathsf{copysign}\left(\log \left(1 + \color{blue}{{x}^{1}}\right), x\right) \]

      5. unpow1 8.1%

         \[\leadsto \mathsf{copysign}\left(\log \left(1 + \color{blue}{x}\right), x\right) \]

   6. Simplified 8.1%

      \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(1 + x\right)}, x\right) \]

   7. Taylor expanded in x around 0 97.4%

      \[\leadsto \mathsf{copysign}\left(\color{blue}{x}, x\right) \]

   if 1.25 < x

   1. Initial program 50.7%

      \[\mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{x \cdot x + 1}\right), x\right) \]
   2. Step-by-step derivation

      1. +-commutative 50.7%

         \[\leadsto \mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{\color{blue}{1 + x \cdot x}}\right), x\right) \]

      2. hypot-1-def 100.0%

         \[\leadsto \mathsf{copysign}\left(\log \left(\left|x\right| + \color{blue}{\mathsf{hypot}\left(1, x\right)}\right), x\right) \]

   3. Simplified 100.0%

      \[\leadsto \color{blue}{\mathsf{copysign}\left(\log \left(\left|x\right| + \mathsf{hypot}\left(1, x\right)\right), x\right)} \]

   4. Taylor expanded in x around inf 99.7%

      \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(\left|x\right| + x\right)}, x\right) \]
   5. Step-by-step derivation

      1. unpow1 99.7%

         \[\leadsto \mathsf{copysign}\left(\log \left(\left|\color{blue}{{x}^{1}}\right| + x\right), x\right) \]

      2. sqr-pow 99.7%

         \[\leadsto \mathsf{copysign}\left(\log \left(\left|\color{blue}{{x}^{\left(\frac{1}{2}\right)} \cdot {x}^{\left(\frac{1}{2}\right)}}\right| + x\right), x\right) \]

      3. fabs-sqr 99.7%

         \[\leadsto \mathsf{copysign}\left(\log \left(\color{blue}{{x}^{\left(\frac{1}{2}\right)} \cdot {x}^{\left(\frac{1}{2}\right)}} + x\right), x\right) \]

      4. sqr-pow 99.7%

         \[\leadsto \mathsf{copysign}\left(\log \left(\color{blue}{{x}^{1}} + x\right), x\right) \]

      5. unpow1 99.7%

         \[\leadsto \mathsf{copysign}\left(\log \left(\color{blue}{x} + x\right), x\right) \]

   6. Simplified 99.7%

      \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(x + x\right)}, x\right) \]
3. Recombined 3 regimes into one program.

4. Final simplification 81.2%

   \[\leadsto \begin{array}{l} \mathbf{if}\;x \leq -3.2:\\ \;\;\;\;\mathsf{copysign}\left(\log \left(-x\right), x\right)\\ \mathbf{elif}\;x \leq 1.25:\\ \;\;\;\;\mathsf{copysign}\left(x, x\right)\\ \mathbf{else}:\\ \;\;\;\;\mathsf{copysign}\left(\log \left(x + x\right), x\right)\\ \end{array} \]

Alternative 8: 99.0% accurate, 2.0× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;x \leq -1.25:\\ \;\;\;\;\mathsf{copysign}\left(\log \left(\frac{-0.5}{x}\right), x\right)\\ \mathbf{elif}\;x \leq 1.25:\\ \;\;\;\;\mathsf{copysign}\left(x, x\right)\\ \mathbf{else}:\\ \;\;\;\;\mathsf{copysign}\left(\log \left(x + x\right), x\right)\\ \end{array} \end{array} \]

FPCore

(FPCore (x)
 :precision binary64
 (if (<= x -1.25)
   (copysign (log (/ -0.5 x)) x)
   (if (<= x 1.25) (copysign x x) (copysign (log (+ x x)) x))))

C

double code(double x) {
	double tmp;
	if (x <= -1.25) {
		tmp = copysign(log((-0.5 / x)), x);
	} else if (x <= 1.25) {
		tmp = copysign(x, x);
	} else {
		tmp = copysign(log((x + x)), x);
	}
	return tmp;
}

Java

public static double code(double x) {
	double tmp;
	if (x <= -1.25) {
		tmp = Math.copySign(Math.log((-0.5 / x)), x);
	} else if (x <= 1.25) {
		tmp = Math.copySign(x, x);
	} else {
		tmp = Math.copySign(Math.log((x + x)), x);
	}
	return tmp;
}

Python

def code(x):
	tmp = 0
	if x <= -1.25:
		tmp = math.copysign(math.log((-0.5 / x)), x)
	elif x <= 1.25:
		tmp = math.copysign(x, x)
	else:
		tmp = math.copysign(math.log((x + x)), x)
	return tmp

Julia

function code(x)
	tmp = 0.0
	if (x <= -1.25)
		tmp = copysign(log(Float64(-0.5 / x)), x);
	elseif (x <= 1.25)
		tmp = copysign(x, x);
	else
		tmp = copysign(log(Float64(x + x)), x);
	end
	return tmp
end

MATLAB

function tmp_2 = code(x)
	tmp = 0.0;
	if (x <= -1.25)
		tmp = sign(x) * abs(log((-0.5 / x)));
	elseif (x <= 1.25)
		tmp = sign(x) * abs(x);
	else
		tmp = sign(x) * abs(log((x + x)));
	end
	tmp_2 = tmp;
end

Wolfram

code[x_] := If[LessEqual[x, -1.25], N[With[{TMP1 = Abs[N[Log[N[(-0.5 / x), $MachinePrecision]], $MachinePrecision]], TMP2 = Sign[x]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision], If[LessEqual[x, 1.25], N[With[{TMP1 = Abs[x], TMP2 = Sign[x]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision], N[With[{TMP1 = Abs[N[Log[N[(x + x), $MachinePrecision]], $MachinePrecision]], TMP2 = Sign[x]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]]]

Derivation

1. Split input into 3 regimes

2. if x < -1.25

   1. Initial program 50.8%

      \[\mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{x \cdot x + 1}\right), x\right) \]
   2. Step-by-step derivation

      1. +-commutative 50.8%

         \[\leadsto \mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{\color{blue}{1 + x \cdot x}}\right), x\right) \]

      2. hypot-1-def 100.0%

         \[\leadsto \mathsf{copysign}\left(\log \left(\left|x\right| + \color{blue}{\mathsf{hypot}\left(1, x\right)}\right), x\right) \]

   3. Simplified 100.0%

      \[\leadsto \color{blue}{\mathsf{copysign}\left(\log \left(\left|x\right| + \mathsf{hypot}\left(1, x\right)\right), x\right)} \]

   4. Taylor expanded in x around -inf 100.0%

      \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(\left(\left|x\right| + -1 \cdot x\right) - 0.5 \cdot \frac{1}{x}\right)}, x\right) \]
   5. Step-by-step derivation

      1. associate--l+ 100.0%

         \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(\left|x\right| + \left(-1 \cdot x - 0.5 \cdot \frac{1}{x}\right)\right)}, x\right) \]

      2. unpow1 100.0%

         \[\leadsto \mathsf{copysign}\left(\log \left(\left|\color{blue}{{x}^{1}}\right| + \left(-1 \cdot x - 0.5 \cdot \frac{1}{x}\right)\right), x\right) \]

      3. sqr-pow 0.0%

         \[\leadsto \mathsf{copysign}\left(\log \left(\left|\color{blue}{{x}^{\left(\frac{1}{2}\right)} \cdot {x}^{\left(\frac{1}{2}\right)}}\right| + \left(-1 \cdot x - 0.5 \cdot \frac{1}{x}\right)\right), x\right) \]

      4. fabs-sqr 0.0%

         \[\leadsto \mathsf{copysign}\left(\log \left(\color{blue}{{x}^{\left(\frac{1}{2}\right)} \cdot {x}^{\left(\frac{1}{2}\right)}} + \left(-1 \cdot x - 0.5 \cdot \frac{1}{x}\right)\right), x\right) \]

      5. sqr-pow 4.1%

         \[\leadsto \mathsf{copysign}\left(\log \left(\color{blue}{{x}^{1}} + \left(-1 \cdot x - 0.5 \cdot \frac{1}{x}\right)\right), x\right) \]

      6. unpow1 4.1%

         \[\leadsto \mathsf{copysign}\left(\log \left(\color{blue}{x} + \left(-1 \cdot x - 0.5 \cdot \frac{1}{x}\right)\right), x\right) \]

      7. associate-+r- 99.6%

         \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(\left(x + -1 \cdot x\right) - 0.5 \cdot \frac{1}{x}\right)}, x\right) \]

      8. mul-1-neg 99.6%

         \[\leadsto \mathsf{copysign}\left(\log \left(\left(x + \color{blue}{\left(-x\right)}\right) - 0.5 \cdot \frac{1}{x}\right), x\right) \]

      9. sub-neg 99.6%

         \[\leadsto \mathsf{copysign}\left(\log \left(\color{blue}{\left(x - x\right)} - 0.5 \cdot \frac{1}{x}\right), x\right) \]

      10. +-inverses 99.6%

          \[\leadsto \mathsf{copysign}\left(\log \left(\color{blue}{0} - 0.5 \cdot \frac{1}{x}\right), x\right) \]

      11. neg-sub0 99.6%

          \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(-0.5 \cdot \frac{1}{x}\right)}, x\right) \]

      12. associate-*r/ 99.6%

          \[\leadsto \mathsf{copysign}\left(\log \left(-\color{blue}{\frac{0.5 \cdot 1}{x}}\right), x\right) \]

      13. metadata-eval 99.6%

          \[\leadsto \mathsf{copysign}\left(\log \left(-\frac{\color{blue}{0.5}}{x}\right), x\right) \]

      14. distribute-neg-frac 99.6%

          \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(\frac{-0.5}{x}\right)}, x\right) \]

      15. metadata-eval 99.6%

          \[\leadsto \mathsf{copysign}\left(\log \left(\frac{\color{blue}{-0.5}}{x}\right), x\right) \]

   6. Simplified 99.6%

      \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(\frac{-0.5}{x}\right)}, x\right) \]

   if -1.25 < x < 1.25

   1. Initial program 10.7%

      \[\mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{x \cdot x + 1}\right), x\right) \]
   2. Step-by-step derivation

      1. +-commutative 10.7%

         \[\leadsto \mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{\color{blue}{1 + x \cdot x}}\right), x\right) \]

      2. hypot-1-def 10.8%

         \[\leadsto \mathsf{copysign}\left(\log \left(\left|x\right| + \color{blue}{\mathsf{hypot}\left(1, x\right)}\right), x\right) \]

   3. Simplified 10.8%

      \[\leadsto \color{blue}{\mathsf{copysign}\left(\log \left(\left|x\right| + \mathsf{hypot}\left(1, x\right)\right), x\right)} \]

   4. Taylor expanded in x around 0 8.0%

      \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(1 + \left|x\right|\right)}, x\right) \]
   5. Step-by-step derivation

      1. unpow1 8.0%

         \[\leadsto \mathsf{copysign}\left(\log \left(1 + \left|\color{blue}{{x}^{1}}\right|\right), x\right) \]

      2. sqr-pow 2.8%

         \[\leadsto \mathsf{copysign}\left(\log \left(1 + \left|\color{blue}{{x}^{\left(\frac{1}{2}\right)} \cdot {x}^{\left(\frac{1}{2}\right)}}\right|\right), x\right) \]

      3. fabs-sqr 2.8%

         \[\leadsto \mathsf{copysign}\left(\log \left(1 + \color{blue}{{x}^{\left(\frac{1}{2}\right)} \cdot {x}^{\left(\frac{1}{2}\right)}}\right), x\right) \]

      4. sqr-pow 8.1%

         \[\leadsto \mathsf{copysign}\left(\log \left(1 + \color{blue}{{x}^{1}}\right), x\right) \]

      5. unpow1 8.1%

         \[\leadsto \mathsf{copysign}\left(\log \left(1 + \color{blue}{x}\right), x\right) \]

   6. Simplified 8.1%

      \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(1 + x\right)}, x\right) \]

   7. Taylor expanded in x around 0 97.4%

      \[\leadsto \mathsf{copysign}\left(\color{blue}{x}, x\right) \]

   if 1.25 < x

   1. Initial program 50.7%

      \[\mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{x \cdot x + 1}\right), x\right) \]
   2. Step-by-step derivation

      1. +-commutative 50.7%

         \[\leadsto \mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{\color{blue}{1 + x \cdot x}}\right), x\right) \]

      2. hypot-1-def 100.0%

         \[\leadsto \mathsf{copysign}\left(\log \left(\left|x\right| + \color{blue}{\mathsf{hypot}\left(1, x\right)}\right), x\right) \]

   3. Simplified 100.0%

      \[\leadsto \color{blue}{\mathsf{copysign}\left(\log \left(\left|x\right| + \mathsf{hypot}\left(1, x\right)\right), x\right)} \]

   4. Taylor expanded in x around inf 99.7%

      \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(\left|x\right| + x\right)}, x\right) \]
   5. Step-by-step derivation

      1. unpow1 99.7%

         \[\leadsto \mathsf{copysign}\left(\log \left(\left|\color{blue}{{x}^{1}}\right| + x\right), x\right) \]

      2. sqr-pow 99.7%

         \[\leadsto \mathsf{copysign}\left(\log \left(\left|\color{blue}{{x}^{\left(\frac{1}{2}\right)} \cdot {x}^{\left(\frac{1}{2}\right)}}\right| + x\right), x\right) \]

      3. fabs-sqr 99.7%

         \[\leadsto \mathsf{copysign}\left(\log \left(\color{blue}{{x}^{\left(\frac{1}{2}\right)} \cdot {x}^{\left(\frac{1}{2}\right)}} + x\right), x\right) \]

      4. sqr-pow 99.7%

         \[\leadsto \mathsf{copysign}\left(\log \left(\color{blue}{{x}^{1}} + x\right), x\right) \]

      5. unpow1 99.7%

         \[\leadsto \mathsf{copysign}\left(\log \left(\color{blue}{x} + x\right), x\right) \]

   6. Simplified 99.7%

      \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(x + x\right)}, x\right) \]
3. Recombined 3 regimes into one program.

4. Final simplification 98.5%

   \[\leadsto \begin{array}{l} \mathbf{if}\;x \leq -1.25:\\ \;\;\;\;\mathsf{copysign}\left(\log \left(\frac{-0.5}{x}\right), x\right)\\ \mathbf{elif}\;x \leq 1.25:\\ \;\;\;\;\mathsf{copysign}\left(x, x\right)\\ \mathbf{else}:\\ \;\;\;\;\mathsf{copysign}\left(\log \left(x + x\right), x\right)\\ \end{array} \]

Alternative 9: 64.6% accurate, 2.0× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;x \leq -0.5:\\ \;\;\;\;\mathsf{copysign}\left(\log \left(-x\right), x\right)\\ \mathbf{else}:\\ \;\;\;\;\mathsf{copysign}\left(\mathsf{log1p}\left(x\right), x\right)\\ \end{array} \end{array} \]

FPCore

(FPCore (x)
 :precision binary64
 (if (<= x -0.5) (copysign (log (- x)) x) (copysign (log1p x) x)))

C

double code(double x) {
	double tmp;
	if (x <= -0.5) {
		tmp = copysign(log(-x), x);
	} else {
		tmp = copysign(log1p(x), x);
	}
	return tmp;
}

Java

public static double code(double x) {
	double tmp;
	if (x <= -0.5) {
		tmp = Math.copySign(Math.log(-x), x);
	} else {
		tmp = Math.copySign(Math.log1p(x), x);
	}
	return tmp;
}

Python

def code(x):
	tmp = 0
	if x <= -0.5:
		tmp = math.copysign(math.log(-x), x)
	else:
		tmp = math.copysign(math.log1p(x), x)
	return tmp

Julia

function code(x)
	tmp = 0.0
	if (x <= -0.5)
		tmp = copysign(log(Float64(-x)), x);
	else
		tmp = copysign(log1p(x), x);
	end
	return tmp
end

Wolfram

code[x_] := If[LessEqual[x, -0.5], N[With[{TMP1 = Abs[N[Log[(-x)], $MachinePrecision]], TMP2 = Sign[x]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision], N[With[{TMP1 = Abs[N[Log[1 + x], $MachinePrecision]], TMP2 = Sign[x]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]]

Derivation

1. Split input into 2 regimes

2. if x < -0.5

   1. Initial program 50.8%

      \[\mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{x \cdot x + 1}\right), x\right) \]
   2. Step-by-step derivation

      1. +-commutative 50.8%

         \[\leadsto \mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{\color{blue}{1 + x \cdot x}}\right), x\right) \]

      2. hypot-1-def 100.0%

         \[\leadsto \mathsf{copysign}\left(\log \left(\left|x\right| + \color{blue}{\mathsf{hypot}\left(1, x\right)}\right), x\right) \]

   3. Simplified 100.0%

      \[\leadsto \color{blue}{\mathsf{copysign}\left(\log \left(\left|x\right| + \mathsf{hypot}\left(1, x\right)\right), x\right)} \]

   4. Taylor expanded in x around -inf 31.5%

      \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(-1 \cdot x\right)}, x\right) \]
   5. Step-by-step derivation

      1. mul-1-neg 31.5%

         \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(-x\right)}, x\right) \]

   6. Simplified 31.5%

      \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(-x\right)}, x\right) \]

   if -0.5 < x

   1. Initial program 22.7%

      \[\mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{x \cdot x + 1}\right), x\right) \]
   2. Step-by-step derivation

      1. +-commutative 22.7%

         \[\leadsto \mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{\color{blue}{1 + x \cdot x}}\right), x\right) \]

      2. hypot-1-def 37.4%

         \[\leadsto \mathsf{copysign}\left(\log \left(\left|x\right| + \color{blue}{\mathsf{hypot}\left(1, x\right)}\right), x\right) \]

   3. Simplified 37.4%

      \[\leadsto \color{blue}{\mathsf{copysign}\left(\log \left(\left|x\right| + \mathsf{hypot}\left(1, x\right)\right), x\right)} \]

   4. Taylor expanded in x around 0 15.0%

      \[\leadsto \mathsf{copysign}\left(\color{blue}{\log \left(1 + \left|x\right|\right)}, x\right) \]
   5. Step-by-step derivation

      1. log1p-def 76.7%

         \[\leadsto \mathsf{copysign}\left(\color{blue}{\mathsf{log1p}\left(\left|x\right|\right)}, x\right) \]

      2. unpow1 76.7%

         \[\leadsto \mathsf{copysign}\left(\mathsf{log1p}\left(\left|\color{blue}{{x}^{1}}\right|\right), x\right) \]

      3. sqr-pow 41.8%

         \[\leadsto \mathsf{copysign}\left(\mathsf{log1p}\left(\left|\color{blue}{{x}^{\left(\frac{1}{2}\right)} \cdot {x}^{\left(\frac{1}{2}\right)}}\right|\right), x\right) \]

      4. fabs-sqr 41.8%

         \[\leadsto \mathsf{copysign}\left(\mathsf{log1p}\left(\color{blue}{{x}^{\left(\frac{1}{2}\right)} \cdot {x}^{\left(\frac{1}{2}\right)}}\right), x\right) \]

      5. sqr-pow 76.7%

         \[\leadsto \mathsf{copysign}\left(\mathsf{log1p}\left(\color{blue}{{x}^{1}}\right), x\right) \]

      6. unpow1 76.7%

         \[\leadsto \mathsf{copysign}\left(\mathsf{log1p}\left(\color{blue}{x}\right), x\right) \]

   6. Simplified 76.7%

      \[\leadsto \mathsf{copysign}\left(\color{blue}{\mathsf{log1p}\left(x\right)}, x\right) \]
3. Recombined 2 regimes into one program.

4. Final simplification 65.2%

   \[\leadsto \begin{array}{l} \mathbf{if}\;x \leq -0.5:\\ \;\;\;\;\mathsf{copysign}\left(\log \left(-x\right), x\right)\\ \mathbf{else}:\\ \;\;\;\;\mathsf{copysign}\left(\mathsf{log1p}\left(x\right), x\right)\\ \end{array} \]

Alternative 10: 59.0% accurate, 2.0× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;x \leq 1.6:\\ \;\;\;\;\mathsf{copysign}\left(x, x\right)\\ \mathbf{else}:\\ \;\;\;\;\mathsf{copysign}\left(\mathsf{log1p}\left(x\right), x\right)\\ \end{array} \end{array} \]

FPCore

(FPCore (x)
 :precision binary64
 (if (<= x 1.6) (copysign x x) (copysign (log1p x) x)))

C

double code(double x) {
	double tmp;
	if (x <= 1.6) {
		tmp = copysign(x, x);
	} else {
		tmp = copysign(log1p(x), x);
	}
	return tmp;
}

Java

public static double code(double x) {
	double tmp;
	if (x <= 1.6) {
		tmp = Math.copySign(x, x);
	} else {
		tmp = Math.copySign(Math.log1p(x), x);
	}
	return tmp;
}

Python

def code(x):
	tmp = 0
	if x <= 1.6:
		tmp = math.copysign(x, x)
	else:
		tmp = math.copysign(math.log1p(x), x)
	return tmp

Julia

function code(x)
	tmp = 0.0
	if (x <= 1.6)
		tmp = copysign(x, x);
	else
		tmp = copysign(log1p(x), x);
	end
	return tmp
end

Wolfram

code[x_] := If[LessEqual[x, 1.6], N[With[{TMP1 = Abs[x], TMP2 = Sign[x]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision], N[With[{TMP1 = Abs[N[Log[1 + x], $MachinePrecision]], TMP2 = Sign[x]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]]

Derivation

1. Split input into 2 regimes

2. if x < 1.6000000000000001

   1. Initial program 23.8%

      \[\mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{x \cdot x + 1}\right), x\right) \]
   2. Step-by-step derivation

      1. +-commutative 23.8%

         \[\leadsto \mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{\color{blue}{1 + x \cdot x}}\right), x\right) \]

      2. hypot-1-def 40.0%

         \[\leadsto \mathsf{copysign}\left(\log \left(\left|x\right| + \color{blue}{\mathsf{hypot}\left(1, x\right)}\right), x\right) \]

   3. Simplified 40.0%

      \[\leadsto \color{blue}{\mathsf{copysign}\left(\log \left(\left|x\right| + \mathsf{hypot}\left(1, x\right)\right), x\right)} \]

   4. Taylor expanded in x around 0 15.7%

      \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(1 + \left|x\right|\right)}, x\right) \]
   5. Step-by-step derivation

      1. unpow1 15.7%

         \[\leadsto \mathsf{copysign}\left(\log \left(1 + \left|\color{blue}{{x}^{1}}\right|\right), x\right) \]

      2. sqr-pow 1.9%

         \[\leadsto \mathsf{copysign}\left(\log \left(1 + \left|\color{blue}{{x}^{\left(\frac{1}{2}\right)} \cdot {x}^{\left(\frac{1}{2}\right)}}\right|\right), x\right) \]

      3. fabs-sqr 1.9%

         \[\leadsto \mathsf{copysign}\left(\log \left(1 + \color{blue}{{x}^{\left(\frac{1}{2}\right)} \cdot {x}^{\left(\frac{1}{2}\right)}}\right), x\right) \]

      4. sqr-pow 5.5%

         \[\leadsto \mathsf{copysign}\left(\log \left(1 + \color{blue}{{x}^{1}}\right), x\right) \]

      5. unpow1 5.5%

         \[\leadsto \mathsf{copysign}\left(\log \left(1 + \color{blue}{x}\right), x\right) \]

   6. Simplified 5.5%

      \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(1 + x\right)}, x\right) \]

   7. Taylor expanded in x around 0 67.4%

      \[\leadsto \mathsf{copysign}\left(\color{blue}{x}, x\right) \]

   if 1.6000000000000001 < x

   1. Initial program 50.7%

      \[\mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{x \cdot x + 1}\right), x\right) \]
   2. Step-by-step derivation

      1. +-commutative 50.7%

         \[\leadsto \mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{\color{blue}{1 + x \cdot x}}\right), x\right) \]

      2. hypot-1-def 100.0%

         \[\leadsto \mathsf{copysign}\left(\log \left(\left|x\right| + \color{blue}{\mathsf{hypot}\left(1, x\right)}\right), x\right) \]

   3. Simplified 100.0%

      \[\leadsto \color{blue}{\mathsf{copysign}\left(\log \left(\left|x\right| + \mathsf{hypot}\left(1, x\right)\right), x\right)} \]

   4. Taylor expanded in x around 0 31.4%

      \[\leadsto \mathsf{copysign}\left(\color{blue}{\log \left(1 + \left|x\right|\right)}, x\right) \]
   5. Step-by-step derivation

      1. log1p-def 31.4%

         \[\leadsto \mathsf{copysign}\left(\color{blue}{\mathsf{log1p}\left(\left|x\right|\right)}, x\right) \]

      2. unpow1 31.4%

         \[\leadsto \mathsf{copysign}\left(\mathsf{log1p}\left(\left|\color{blue}{{x}^{1}}\right|\right), x\right) \]

      3. sqr-pow 31.4%

         \[\leadsto \mathsf{copysign}\left(\mathsf{log1p}\left(\left|\color{blue}{{x}^{\left(\frac{1}{2}\right)} \cdot {x}^{\left(\frac{1}{2}\right)}}\right|\right), x\right) \]

      4. fabs-sqr 31.4%

         \[\leadsto \mathsf{copysign}\left(\mathsf{log1p}\left(\color{blue}{{x}^{\left(\frac{1}{2}\right)} \cdot {x}^{\left(\frac{1}{2}\right)}}\right), x\right) \]

      5. sqr-pow 31.4%

         \[\leadsto \mathsf{copysign}\left(\mathsf{log1p}\left(\color{blue}{{x}^{1}}\right), x\right) \]

      6. unpow1 31.4%

         \[\leadsto \mathsf{copysign}\left(\mathsf{log1p}\left(\color{blue}{x}\right), x\right) \]

   6. Simplified 31.4%

      \[\leadsto \mathsf{copysign}\left(\color{blue}{\mathsf{log1p}\left(x\right)}, x\right) \]
3. Recombined 2 regimes into one program.

4. Final simplification 59.4%

   \[\leadsto \begin{array}{l} \mathbf{if}\;x \leq 1.6:\\ \;\;\;\;\mathsf{copysign}\left(x, x\right)\\ \mathbf{else}:\\ \;\;\;\;\mathsf{copysign}\left(\mathsf{log1p}\left(x\right), x\right)\\ \end{array} \]

Alternative 11: 52.2% accurate, 4.0× speedup

Math

\[\begin{array}{l} \\ \mathsf{copysign}\left(x, x\right) \end{array} \]

FPCore

(FPCore (x) :precision binary64 (copysign x x))

C

double code(double x) {
	return copysign(x, x);
}

Java

public static double code(double x) {
	return Math.copySign(x, x);
}

Python

def code(x):
	return math.copysign(x, x)

Julia

function code(x)
	return copysign(x, x)
end

MATLAB

function tmp = code(x)
	tmp = sign(x) * abs(x);
end

Wolfram

code[x_] := N[With[{TMP1 = Abs[x], TMP2 = Sign[x]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]

Derivation

1. Initial program 29.8%

   \[\mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{x \cdot x + 1}\right), x\right) \]
2. Step-by-step derivation

   1. +-commutative 29.8%

      \[\leadsto \mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{\color{blue}{1 + x \cdot x}}\right), x\right) \]

   2. hypot-1-def 53.3%

      \[\leadsto \mathsf{copysign}\left(\log \left(\left|x\right| + \color{blue}{\mathsf{hypot}\left(1, x\right)}\right), x\right) \]

3. Simplified 53.3%

   \[\leadsto \color{blue}{\mathsf{copysign}\left(\log \left(\left|x\right| + \mathsf{hypot}\left(1, x\right)\right), x\right)} \]

4. Taylor expanded in x around 0 19.2%

   \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(1 + \left|x\right|\right)}, x\right) \]
5. Step-by-step derivation

   1. unpow1 19.2%

      \[\leadsto \mathsf{copysign}\left(\log \left(1 + \left|\color{blue}{{x}^{1}}\right|\right), x\right) \]

   2. sqr-pow 8.5%

      \[\leadsto \mathsf{copysign}\left(\log \left(1 + \left|\color{blue}{{x}^{\left(\frac{1}{2}\right)} \cdot {x}^{\left(\frac{1}{2}\right)}}\right|\right), x\right) \]

   3. fabs-sqr 8.5%

      \[\leadsto \mathsf{copysign}\left(\log \left(1 + \color{blue}{{x}^{\left(\frac{1}{2}\right)} \cdot {x}^{\left(\frac{1}{2}\right)}}\right), x\right) \]

   4. sqr-pow 11.2%

      \[\leadsto \mathsf{copysign}\left(\log \left(1 + \color{blue}{{x}^{1}}\right), x\right) \]

   5. unpow1 11.2%

      \[\leadsto \mathsf{copysign}\left(\log \left(1 + \color{blue}{x}\right), x\right) \]

6. Simplified 11.2%

   \[\leadsto \mathsf{copysign}\left(\log \color{blue}{\left(1 + x\right)}, x\right) \]

7. Taylor expanded in x around 0 53.6%

   \[\leadsto \mathsf{copysign}\left(\color{blue}{x}, x\right) \]

8. Final simplification 53.6%

   \[\leadsto \mathsf{copysign}\left(x, x\right) \]

Developer target: 99.9% accurate, 0.6× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} t_0 := \frac{1}{\left|x\right|}\\ \mathsf{copysign}\left(\mathsf{log1p}\left(\left|x\right| + \frac{\left|x\right|}{\mathsf{hypot}\left(1, t_0\right) + t_0}\right), x\right) \end{array} \end{array} \]

FPCore

(FPCore (x)
 :precision binary64
 (let* ((t_0 (/ 1.0 (fabs x))))
   (copysign (log1p (+ (fabs x) (/ (fabs x) (+ (hypot 1.0 t_0) t_0)))) x)))

C

double code(double x) {
	double t_0 = 1.0 / fabs(x);
	return copysign(log1p((fabs(x) + (fabs(x) / (hypot(1.0, t_0) + t_0)))), x);
}

Java

public static double code(double x) {
	double t_0 = 1.0 / Math.abs(x);
	return Math.copySign(Math.log1p((Math.abs(x) + (Math.abs(x) / (Math.hypot(1.0, t_0) + t_0)))), x);
}

Python

def code(x):
	t_0 = 1.0 / math.fabs(x)
	return math.copysign(math.log1p((math.fabs(x) + (math.fabs(x) / (math.hypot(1.0, t_0) + t_0)))), x)

Julia

function code(x)
	t_0 = Float64(1.0 / abs(x))
	return copysign(log1p(Float64(abs(x) + Float64(abs(x) / Float64(hypot(1.0, t_0) + t_0)))), x)
end

Wolfram

code[x_] := Block[{t$95$0 = N[(1.0 / N[Abs[x], $MachinePrecision]), $MachinePrecision]}, N[With[{TMP1 = Abs[N[Log[1 + N[(N[Abs[x], $MachinePrecision] + N[(N[Abs[x], $MachinePrecision] / N[(N[Sqrt[1.0 ^ 2 + t$95$0 ^ 2], $MachinePrecision] + t$95$0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], TMP2 = Sign[x]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]]

Reproduce

herbie shell --seed 2023278 
(FPCore (x)
  :name "Rust f64::asinh"
  :precision binary64

  :herbie-target
  (copysign (log1p (+ (fabs x) (/ (fabs x) (+ (hypot 1.0 (/ 1.0 (fabs x))) (/ 1.0 (fabs x)))))) x)

  (copysign (log (+ (fabs x) (sqrt (+ (* x x) 1.0)))) x))