
(FPCore (x) :precision binary64 (asinh x))
double code(double x) {
return asinh(x);
}
def code(x): return math.asinh(x)
function code(x) return asinh(x) end
function tmp = code(x) tmp = asinh(x); end
code[x_] := N[ArcSinh[x], $MachinePrecision]
\begin{array}{l}
\\
\sinh^{-1} x
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 14 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x) :precision binary64 (copysign (log (+ (fabs x) (sqrt (+ (* x x) 1.0)))) x))
double code(double x) {
return copysign(log((fabs(x) + sqrt(((x * x) + 1.0)))), x);
}
public static double code(double x) {
return Math.copySign(Math.log((Math.abs(x) + Math.sqrt(((x * x) + 1.0)))), x);
}
def code(x): return math.copysign(math.log((math.fabs(x) + math.sqrt(((x * x) + 1.0)))), x)
function code(x) return copysign(log(Float64(abs(x) + sqrt(Float64(Float64(x * x) + 1.0)))), x) end
function tmp = code(x) tmp = sign(x) * abs(log((abs(x) + sqrt(((x * x) + 1.0))))); end
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]
\begin{array}{l}
\\
\mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{x \cdot x + 1}\right), x\right)
\end{array}
(FPCore (x)
:precision binary64
(let* ((t_0 (copysign (log (+ (fabs x) (sqrt (+ (* x x) 1.0)))) x)))
(if (<= t_0 -0.1)
(copysign (log (/ 1.0 (- (hypot 1.0 x) x))) x)
(if (<= t_0 0.01)
(copysign
(+
(* -0.16666666666666666 (pow x 3.0))
(+ (* 0.075 (pow x 5.0)) (+ x (* -0.044642857142857144 (pow x 7.0)))))
x)
(copysign (log (+ x (+ x (/ 0.5 x)))) x)))))
double code(double x) {
double t_0 = copysign(log((fabs(x) + sqrt(((x * x) + 1.0)))), x);
double tmp;
if (t_0 <= -0.1) {
tmp = copysign(log((1.0 / (hypot(1.0, x) - x))), x);
} else if (t_0 <= 0.01) {
tmp = copysign(((-0.16666666666666666 * pow(x, 3.0)) + ((0.075 * pow(x, 5.0)) + (x + (-0.044642857142857144 * pow(x, 7.0))))), x);
} else {
tmp = copysign(log((x + (x + (0.5 / x)))), x);
}
return tmp;
}
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.1) {
tmp = Math.copySign(Math.log((1.0 / (Math.hypot(1.0, x) - x))), x);
} else if (t_0 <= 0.01) {
tmp = Math.copySign(((-0.16666666666666666 * Math.pow(x, 3.0)) + ((0.075 * Math.pow(x, 5.0)) + (x + (-0.044642857142857144 * Math.pow(x, 7.0))))), x);
} else {
tmp = Math.copySign(Math.log((x + (x + (0.5 / x)))), x);
}
return tmp;
}
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.1: tmp = math.copysign(math.log((1.0 / (math.hypot(1.0, x) - x))), x) elif t_0 <= 0.01: tmp = math.copysign(((-0.16666666666666666 * math.pow(x, 3.0)) + ((0.075 * math.pow(x, 5.0)) + (x + (-0.044642857142857144 * math.pow(x, 7.0))))), x) else: tmp = math.copysign(math.log((x + (x + (0.5 / x)))), x) return tmp
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.1) tmp = copysign(log(Float64(1.0 / Float64(hypot(1.0, x) - x))), x); elseif (t_0 <= 0.01) tmp = copysign(Float64(Float64(-0.16666666666666666 * (x ^ 3.0)) + Float64(Float64(0.075 * (x ^ 5.0)) + Float64(x + Float64(-0.044642857142857144 * (x ^ 7.0))))), x); else tmp = copysign(log(Float64(x + Float64(x + Float64(0.5 / x)))), x); end return tmp end
function tmp_2 = code(x) t_0 = sign(x) * abs(log((abs(x) + sqrt(((x * x) + 1.0))))); tmp = 0.0; if (t_0 <= -0.1) tmp = sign(x) * abs(log((1.0 / (hypot(1.0, x) - x)))); elseif (t_0 <= 0.01) tmp = sign(x) * abs(((-0.16666666666666666 * (x ^ 3.0)) + ((0.075 * (x ^ 5.0)) + (x + (-0.044642857142857144 * (x ^ 7.0)))))); else tmp = sign(x) * abs(log((x + (x + (0.5 / x))))); end tmp_2 = tmp; end
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.1], 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.01], 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[(-0.044642857142857144 * N[Power[x, 7.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $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]]]]
\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.1:\\
\;\;\;\;\mathsf{copysign}\left(\log \left(\frac{1}{\mathsf{hypot}\left(1, x\right) - x}\right), x\right)\\
\mathbf{elif}\;t_0 \leq 0.01:\\
\;\;\;\;\mathsf{copysign}\left(-0.16666666666666666 \cdot {x}^{3} + \left(0.075 \cdot {x}^{5} + \left(x + -0.044642857142857144 \cdot {x}^{7}\right)\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}
if (copysign.f64 (log.f64 (+.f64 (fabs.f64 x) (sqrt.f64 (+.f64 (*.f64 x x) 1)))) x) < -0.10000000000000001Initial program 60.6%
+-commutative60.6%
hypot-1-def100.0%
Simplified100.0%
flip-+7.4%
div-sub7.4%
pow27.4%
add-sqr-sqrt0.0%
fabs-sqr0.0%
add-sqr-sqrt7.4%
pow27.4%
add-sqr-sqrt0.0%
fabs-sqr0.0%
add-sqr-sqrt1.7%
hypot-udef1.7%
hypot-udef1.7%
add-sqr-sqrt1.7%
metadata-eval1.7%
add-sqr-sqrt0.0%
fabs-sqr0.0%
add-sqr-sqrt8.9%
Applied egg-rr8.9%
unpow28.8%
div-sub10.0%
unpow210.0%
unpow210.0%
unpow210.0%
+-commutative10.0%
associate--r+59.3%
+-inverses100.0%
metadata-eval100.0%
metadata-eval100.0%
associate-/r*100.0%
neg-mul-1100.0%
sub-neg100.0%
+-commutative100.0%
distribute-neg-in100.0%
remove-double-neg100.0%
sub-neg100.0%
Simplified100.0%
if -0.10000000000000001 < (copysign.f64 (log.f64 (+.f64 (fabs.f64 x) (sqrt.f64 (+.f64 (*.f64 x x) 1)))) x) < 0.0100000000000000002Initial program 8.0%
+-commutative8.0%
hypot-1-def8.0%
Simplified8.0%
flip-+8.0%
div-sub8.0%
pow28.0%
add-sqr-sqrt4.6%
fabs-sqr4.6%
add-sqr-sqrt8.0%
pow28.0%
add-sqr-sqrt4.6%
fabs-sqr4.6%
add-sqr-sqrt8.0%
hypot-udef8.0%
hypot-udef8.0%
add-sqr-sqrt8.0%
metadata-eval8.0%
add-sqr-sqrt4.6%
fabs-sqr4.6%
add-sqr-sqrt8.0%
Applied egg-rr8.0%
unpow28.0%
div-sub8.0%
unpow28.0%
unpow28.0%
unpow28.0%
+-commutative8.0%
associate--r+8.0%
+-inverses8.0%
metadata-eval8.0%
metadata-eval8.0%
associate-/r*8.0%
neg-mul-18.0%
sub-neg8.0%
+-commutative8.0%
distribute-neg-in8.0%
remove-double-neg8.0%
sub-neg8.0%
Simplified8.0%
Taylor expanded in x around 0 100.0%
if 0.0100000000000000002 < (copysign.f64 (log.f64 (+.f64 (fabs.f64 x) (sqrt.f64 (+.f64 (*.f64 x x) 1)))) x) Initial program 61.0%
+-commutative61.0%
hypot-1-def100.0%
Simplified100.0%
Taylor expanded in x around inf 100.0%
associate-+r+100.0%
+-commutative100.0%
+-commutative100.0%
unpow1100.0%
sqr-pow100.0%
fabs-sqr100.0%
sqr-pow100.0%
unpow1100.0%
associate-*r/100.0%
metadata-eval100.0%
Simplified100.0%
Final simplification100.0%
(FPCore (x)
:precision binary64
(let* ((t_0 (copysign (log (+ (fabs x) (sqrt (+ (* x x) 1.0)))) x)))
(if (<= t_0 -0.1)
(copysign (log (/ 1.0 (- (hypot 1.0 x) x))) x)
(if (<= t_0 1e-5)
(copysign (+ x (* -0.16666666666666666 (pow x 3.0))) x)
(copysign (* 2.0 (log (sqrt (+ x (hypot 1.0 x))))) x)))))
double code(double x) {
double t_0 = copysign(log((fabs(x) + sqrt(((x * x) + 1.0)))), x);
double tmp;
if (t_0 <= -0.1) {
tmp = copysign(log((1.0 / (hypot(1.0, x) - x))), x);
} else if (t_0 <= 1e-5) {
tmp = copysign((x + (-0.16666666666666666 * pow(x, 3.0))), x);
} else {
tmp = copysign((2.0 * log(sqrt((x + hypot(1.0, x))))), x);
}
return tmp;
}
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.1) {
tmp = Math.copySign(Math.log((1.0 / (Math.hypot(1.0, x) - x))), x);
} else if (t_0 <= 1e-5) {
tmp = Math.copySign((x + (-0.16666666666666666 * Math.pow(x, 3.0))), x);
} else {
tmp = Math.copySign((2.0 * Math.log(Math.sqrt((x + Math.hypot(1.0, x))))), x);
}
return tmp;
}
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.1: tmp = math.copysign(math.log((1.0 / (math.hypot(1.0, x) - x))), x) elif t_0 <= 1e-5: tmp = math.copysign((x + (-0.16666666666666666 * math.pow(x, 3.0))), x) else: tmp = math.copysign((2.0 * math.log(math.sqrt((x + math.hypot(1.0, x))))), x) return tmp
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.1) tmp = copysign(log(Float64(1.0 / Float64(hypot(1.0, x) - x))), x); elseif (t_0 <= 1e-5) tmp = copysign(Float64(x + Float64(-0.16666666666666666 * (x ^ 3.0))), x); else tmp = copysign(Float64(2.0 * log(sqrt(Float64(x + hypot(1.0, x))))), x); end return tmp end
function tmp_2 = code(x) t_0 = sign(x) * abs(log((abs(x) + sqrt(((x * x) + 1.0))))); tmp = 0.0; if (t_0 <= -0.1) tmp = sign(x) * abs(log((1.0 / (hypot(1.0, x) - x)))); elseif (t_0 <= 1e-5) tmp = sign(x) * abs((x + (-0.16666666666666666 * (x ^ 3.0)))); else tmp = sign(x) * abs((2.0 * log(sqrt((x + hypot(1.0, x)))))); end tmp_2 = tmp; end
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.1], 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, 1e-5], N[With[{TMP1 = Abs[N[(x + N[(-0.16666666666666666 * N[Power[x, 3.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], TMP2 = Sign[x]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision], N[With[{TMP1 = Abs[N[(2.0 * N[Log[N[Sqrt[N[(x + N[Sqrt[1.0 ^ 2 + x ^ 2], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]], TMP2 = Sign[x]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]]]]
\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.1:\\
\;\;\;\;\mathsf{copysign}\left(\log \left(\frac{1}{\mathsf{hypot}\left(1, x\right) - x}\right), x\right)\\
\mathbf{elif}\;t_0 \leq 10^{-5}:\\
\;\;\;\;\mathsf{copysign}\left(x + -0.16666666666666666 \cdot {x}^{3}, x\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{copysign}\left(2 \cdot \log \left(\sqrt{x + \mathsf{hypot}\left(1, x\right)}\right), x\right)\\
\end{array}
\end{array}
if (copysign.f64 (log.f64 (+.f64 (fabs.f64 x) (sqrt.f64 (+.f64 (*.f64 x x) 1)))) x) < -0.10000000000000001Initial program 60.6%
+-commutative60.6%
hypot-1-def100.0%
Simplified100.0%
flip-+7.4%
div-sub7.4%
pow27.4%
add-sqr-sqrt0.0%
fabs-sqr0.0%
add-sqr-sqrt7.4%
pow27.4%
add-sqr-sqrt0.0%
fabs-sqr0.0%
add-sqr-sqrt1.7%
hypot-udef1.7%
hypot-udef1.7%
add-sqr-sqrt1.7%
metadata-eval1.7%
add-sqr-sqrt0.0%
fabs-sqr0.0%
add-sqr-sqrt8.9%
Applied egg-rr8.9%
unpow28.8%
div-sub10.0%
unpow210.0%
unpow210.0%
unpow210.0%
+-commutative10.0%
associate--r+59.3%
+-inverses100.0%
metadata-eval100.0%
metadata-eval100.0%
associate-/r*100.0%
neg-mul-1100.0%
sub-neg100.0%
+-commutative100.0%
distribute-neg-in100.0%
remove-double-neg100.0%
sub-neg100.0%
Simplified100.0%
if -0.10000000000000001 < (copysign.f64 (log.f64 (+.f64 (fabs.f64 x) (sqrt.f64 (+.f64 (*.f64 x x) 1)))) x) < 1.00000000000000008e-5Initial program 7.3%
+-commutative7.3%
hypot-1-def7.3%
Simplified7.3%
flip-+7.3%
div-sub7.3%
pow27.3%
add-sqr-sqrt3.9%
fabs-sqr3.9%
add-sqr-sqrt7.3%
pow27.3%
add-sqr-sqrt3.9%
fabs-sqr3.9%
add-sqr-sqrt7.3%
hypot-udef7.3%
hypot-udef7.3%
add-sqr-sqrt7.3%
metadata-eval7.3%
add-sqr-sqrt3.9%
fabs-sqr3.9%
add-sqr-sqrt7.3%
Applied egg-rr7.3%
unpow27.3%
div-sub7.3%
unpow27.3%
unpow27.3%
unpow27.3%
+-commutative7.3%
associate--r+7.3%
+-inverses7.3%
metadata-eval7.3%
metadata-eval7.3%
associate-/r*7.3%
neg-mul-17.3%
sub-neg7.3%
+-commutative7.3%
distribute-neg-in7.3%
remove-double-neg7.3%
sub-neg7.3%
Simplified7.3%
Taylor expanded in x around 0 100.0%
if 1.00000000000000008e-5 < (copysign.f64 (log.f64 (+.f64 (fabs.f64 x) (sqrt.f64 (+.f64 (*.f64 x x) 1)))) x) Initial program 61.4%
+-commutative61.4%
hypot-1-def99.9%
Simplified99.9%
*-un-lft-identity99.9%
log-prod99.9%
metadata-eval99.9%
*-un-lft-identity99.9%
*-un-lft-identity99.9%
add-sqr-sqrt99.9%
fabs-sqr99.9%
add-sqr-sqrt99.9%
Applied egg-rr99.9%
+-lft-identity99.9%
Simplified99.9%
add-sqr-sqrt99.9%
pow299.9%
log-pow100.0%
Applied egg-rr100.0%
Final simplification100.0%
(FPCore (x)
:precision binary64
(if (<= x -0.00105)
(copysign (log (/ 1.0 (- (hypot 1.0 x) x))) x)
(if (<= x 0.00078)
(copysign (+ x (* -0.16666666666666666 (pow x 3.0))) x)
(copysign (log (+ (fabs x) (hypot 1.0 x))) x))))
double code(double x) {
double tmp;
if (x <= -0.00105) {
tmp = copysign(log((1.0 / (hypot(1.0, x) - x))), x);
} else if (x <= 0.00078) {
tmp = copysign((x + (-0.16666666666666666 * pow(x, 3.0))), x);
} else {
tmp = copysign(log((fabs(x) + hypot(1.0, x))), x);
}
return tmp;
}
public static double code(double x) {
double tmp;
if (x <= -0.00105) {
tmp = Math.copySign(Math.log((1.0 / (Math.hypot(1.0, x) - x))), x);
} else if (x <= 0.00078) {
tmp = Math.copySign((x + (-0.16666666666666666 * Math.pow(x, 3.0))), x);
} else {
tmp = Math.copySign(Math.log((Math.abs(x) + Math.hypot(1.0, x))), x);
}
return tmp;
}
def code(x): tmp = 0 if x <= -0.00105: tmp = math.copysign(math.log((1.0 / (math.hypot(1.0, x) - x))), x) elif x <= 0.00078: tmp = math.copysign((x + (-0.16666666666666666 * math.pow(x, 3.0))), x) else: tmp = math.copysign(math.log((math.fabs(x) + math.hypot(1.0, x))), x) return tmp
function code(x) tmp = 0.0 if (x <= -0.00105) tmp = copysign(log(Float64(1.0 / Float64(hypot(1.0, x) - x))), x); elseif (x <= 0.00078) tmp = copysign(Float64(x + Float64(-0.16666666666666666 * (x ^ 3.0))), x); else tmp = copysign(log(Float64(abs(x) + hypot(1.0, x))), x); end return tmp end
function tmp_2 = code(x) tmp = 0.0; if (x <= -0.00105) tmp = sign(x) * abs(log((1.0 / (hypot(1.0, x) - x)))); elseif (x <= 0.00078) tmp = sign(x) * abs((x + (-0.16666666666666666 * (x ^ 3.0)))); else tmp = sign(x) * abs(log((abs(x) + hypot(1.0, x)))); end tmp_2 = tmp; end
code[x_] := If[LessEqual[x, -0.00105], 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[x, 0.00078], N[With[{TMP1 = Abs[N[(x + N[(-0.16666666666666666 * N[Power[x, 3.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], TMP2 = Sign[x]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision], N[With[{TMP1 = Abs[N[Log[N[(N[Abs[x], $MachinePrecision] + N[Sqrt[1.0 ^ 2 + x ^ 2], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], TMP2 = Sign[x]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -0.00105:\\
\;\;\;\;\mathsf{copysign}\left(\log \left(\frac{1}{\mathsf{hypot}\left(1, x\right) - x}\right), x\right)\\
\mathbf{elif}\;x \leq 0.00078:\\
\;\;\;\;\mathsf{copysign}\left(x + -0.16666666666666666 \cdot {x}^{3}, x\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{copysign}\left(\log \left(\left|x\right| + \mathsf{hypot}\left(1, x\right)\right), x\right)\\
\end{array}
\end{array}
if x < -0.00104999999999999994Initial program 60.6%
+-commutative60.6%
hypot-1-def100.0%
Simplified100.0%
flip-+7.4%
div-sub7.4%
pow27.4%
add-sqr-sqrt0.0%
fabs-sqr0.0%
add-sqr-sqrt7.4%
pow27.4%
add-sqr-sqrt0.0%
fabs-sqr0.0%
add-sqr-sqrt1.7%
hypot-udef1.7%
hypot-udef1.7%
add-sqr-sqrt1.7%
metadata-eval1.7%
add-sqr-sqrt0.0%
fabs-sqr0.0%
add-sqr-sqrt8.9%
Applied egg-rr8.9%
unpow28.8%
div-sub10.0%
unpow210.0%
unpow210.0%
unpow210.0%
+-commutative10.0%
associate--r+59.3%
+-inverses100.0%
metadata-eval100.0%
metadata-eval100.0%
associate-/r*100.0%
neg-mul-1100.0%
sub-neg100.0%
+-commutative100.0%
distribute-neg-in100.0%
remove-double-neg100.0%
sub-neg100.0%
Simplified100.0%
if -0.00104999999999999994 < x < 7.79999999999999986e-4Initial program 7.3%
+-commutative7.3%
hypot-1-def7.3%
Simplified7.3%
flip-+7.3%
div-sub7.3%
pow27.3%
add-sqr-sqrt3.9%
fabs-sqr3.9%
add-sqr-sqrt7.3%
pow27.3%
add-sqr-sqrt3.9%
fabs-sqr3.9%
add-sqr-sqrt7.3%
hypot-udef7.3%
hypot-udef7.3%
add-sqr-sqrt7.3%
metadata-eval7.3%
add-sqr-sqrt3.9%
fabs-sqr3.9%
add-sqr-sqrt7.3%
Applied egg-rr7.3%
unpow27.3%
div-sub7.3%
unpow27.3%
unpow27.3%
unpow27.3%
+-commutative7.3%
associate--r+7.3%
+-inverses7.3%
metadata-eval7.3%
metadata-eval7.3%
associate-/r*7.3%
neg-mul-17.3%
sub-neg7.3%
+-commutative7.3%
distribute-neg-in7.3%
remove-double-neg7.3%
sub-neg7.3%
Simplified7.3%
Taylor expanded in x around 0 100.0%
if 7.79999999999999986e-4 < x Initial program 61.4%
+-commutative61.4%
hypot-1-def99.9%
Simplified99.9%
Final simplification100.0%
(FPCore (x)
:precision binary64
(if (<= x -0.95)
(copysign (log (- (* x -2.0) (/ 0.5 x))) x)
(if (<= x 0.00078)
(copysign (+ x (* -0.16666666666666666 (pow x 3.0))) x)
(copysign (log (+ x (hypot 1.0 x))) x))))
double code(double x) {
double tmp;
if (x <= -0.95) {
tmp = copysign(log(((x * -2.0) - (0.5 / x))), x);
} else if (x <= 0.00078) {
tmp = copysign((x + (-0.16666666666666666 * pow(x, 3.0))), x);
} else {
tmp = copysign(log((x + hypot(1.0, x))), x);
}
return tmp;
}
public static double code(double x) {
double tmp;
if (x <= -0.95) {
tmp = Math.copySign(Math.log(((x * -2.0) - (0.5 / x))), x);
} else if (x <= 0.00078) {
tmp = Math.copySign((x + (-0.16666666666666666 * Math.pow(x, 3.0))), x);
} else {
tmp = Math.copySign(Math.log((x + Math.hypot(1.0, x))), x);
}
return tmp;
}
def code(x): tmp = 0 if x <= -0.95: tmp = math.copysign(math.log(((x * -2.0) - (0.5 / x))), x) elif x <= 0.00078: tmp = math.copysign((x + (-0.16666666666666666 * math.pow(x, 3.0))), x) else: tmp = math.copysign(math.log((x + math.hypot(1.0, x))), x) return tmp
function code(x) tmp = 0.0 if (x <= -0.95) tmp = copysign(log(Float64(Float64(x * -2.0) - Float64(0.5 / x))), x); elseif (x <= 0.00078) tmp = copysign(Float64(x + Float64(-0.16666666666666666 * (x ^ 3.0))), x); else tmp = copysign(log(Float64(x + hypot(1.0, x))), x); end return tmp end
function tmp_2 = code(x) tmp = 0.0; if (x <= -0.95) tmp = sign(x) * abs(log(((x * -2.0) - (0.5 / x)))); elseif (x <= 0.00078) tmp = sign(x) * abs((x + (-0.16666666666666666 * (x ^ 3.0)))); else tmp = sign(x) * abs(log((x + hypot(1.0, x)))); end tmp_2 = tmp; end
code[x_] := If[LessEqual[x, -0.95], N[With[{TMP1 = Abs[N[Log[N[(N[(x * -2.0), $MachinePrecision] - N[(0.5 / x), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], TMP2 = Sign[x]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision], If[LessEqual[x, 0.00078], N[With[{TMP1 = Abs[N[(x + N[(-0.16666666666666666 * N[Power[x, 3.0], $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]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -0.95:\\
\;\;\;\;\mathsf{copysign}\left(\log \left(x \cdot -2 - \frac{0.5}{x}\right), x\right)\\
\mathbf{elif}\;x \leq 0.00078:\\
\;\;\;\;\mathsf{copysign}\left(x + -0.16666666666666666 \cdot {x}^{3}, x\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{copysign}\left(\log \left(x + \mathsf{hypot}\left(1, x\right)\right), x\right)\\
\end{array}
\end{array}
if x < -0.94999999999999996Initial program 59.9%
+-commutative59.9%
hypot-1-def100.0%
Simplified100.0%
flip-+5.8%
div-sub5.8%
pow25.8%
add-sqr-sqrt0.0%
fabs-sqr0.0%
add-sqr-sqrt5.8%
pow25.8%
add-sqr-sqrt0.0%
fabs-sqr0.0%
add-sqr-sqrt1.2%
hypot-udef1.2%
hypot-udef1.2%
add-sqr-sqrt1.2%
metadata-eval1.2%
add-sqr-sqrt0.0%
fabs-sqr0.0%
add-sqr-sqrt7.4%
Applied egg-rr7.4%
unpow27.3%
div-sub8.4%
unpow28.5%
unpow28.4%
unpow28.5%
+-commutative8.5%
associate--r+58.6%
+-inverses100.0%
metadata-eval100.0%
metadata-eval100.0%
associate-/r*100.0%
neg-mul-1100.0%
sub-neg100.0%
+-commutative100.0%
distribute-neg-in100.0%
remove-double-neg100.0%
sub-neg100.0%
Simplified100.0%
add-sqr-sqrt0.0%
sqrt-unprod100.0%
pow1/2100.0%
log-rec100.0%
log-rec100.0%
sqr-neg100.0%
pow1100.0%
pow1100.0%
pow-prod-up100.0%
metadata-eval100.0%
Applied egg-rr100.0%
unpow1/2100.0%
unpow2100.0%
rem-sqrt-square100.0%
unpow1100.0%
sqr-pow99.5%
fabs-sqr99.5%
sqr-pow100.0%
unpow1100.0%
Simplified100.0%
Taylor expanded in x around -inf 97.7%
*-commutative97.7%
associate-*r/97.7%
metadata-eval97.7%
Simplified97.7%
if -0.94999999999999996 < x < 7.79999999999999986e-4Initial program 8.0%
+-commutative8.0%
hypot-1-def8.0%
Simplified8.0%
flip-+8.0%
div-sub8.0%
pow28.0%
add-sqr-sqrt3.9%
fabs-sqr3.9%
add-sqr-sqrt8.0%
pow28.0%
add-sqr-sqrt3.9%
fabs-sqr3.9%
add-sqr-sqrt7.5%
hypot-udef7.5%
hypot-udef7.4%
add-sqr-sqrt7.5%
metadata-eval7.5%
add-sqr-sqrt3.8%
fabs-sqr3.8%
add-sqr-sqrt8.0%
Applied egg-rr8.0%
unpow28.0%
div-sub8.0%
unpow28.0%
unpow28.0%
unpow28.0%
+-commutative8.0%
associate--r+8.0%
+-inverses8.0%
metadata-eval8.0%
metadata-eval8.0%
associate-/r*8.0%
neg-mul-18.0%
sub-neg8.0%
+-commutative8.0%
distribute-neg-in8.0%
remove-double-neg8.0%
sub-neg8.0%
Simplified8.0%
Taylor expanded in x around 0 99.5%
if 7.79999999999999986e-4 < x Initial program 61.4%
+-commutative61.4%
hypot-1-def99.9%
Simplified99.9%
*-un-lft-identity99.9%
log-prod99.9%
metadata-eval99.9%
*-un-lft-identity99.9%
*-un-lft-identity99.9%
add-sqr-sqrt99.9%
fabs-sqr99.9%
add-sqr-sqrt99.9%
Applied egg-rr99.9%
+-lft-identity99.9%
Simplified99.9%
Final simplification99.2%
(FPCore (x)
:precision binary64
(if (<= x -0.00105)
(copysign (log (- (hypot 1.0 x) x)) x)
(if (<= x 0.00078)
(copysign (+ x (* -0.16666666666666666 (pow x 3.0))) x)
(copysign (log (+ x (hypot 1.0 x))) x))))
double code(double x) {
double tmp;
if (x <= -0.00105) {
tmp = copysign(log((hypot(1.0, x) - x)), x);
} else if (x <= 0.00078) {
tmp = copysign((x + (-0.16666666666666666 * pow(x, 3.0))), x);
} else {
tmp = copysign(log((x + hypot(1.0, x))), x);
}
return tmp;
}
public static double code(double x) {
double tmp;
if (x <= -0.00105) {
tmp = Math.copySign(Math.log((Math.hypot(1.0, x) - x)), x);
} else if (x <= 0.00078) {
tmp = Math.copySign((x + (-0.16666666666666666 * Math.pow(x, 3.0))), x);
} else {
tmp = Math.copySign(Math.log((x + Math.hypot(1.0, x))), x);
}
return tmp;
}
def code(x): tmp = 0 if x <= -0.00105: tmp = math.copysign(math.log((math.hypot(1.0, x) - x)), x) elif x <= 0.00078: tmp = math.copysign((x + (-0.16666666666666666 * math.pow(x, 3.0))), x) else: tmp = math.copysign(math.log((x + math.hypot(1.0, x))), x) return tmp
function code(x) tmp = 0.0 if (x <= -0.00105) tmp = copysign(log(Float64(hypot(1.0, x) - x)), x); elseif (x <= 0.00078) tmp = copysign(Float64(x + Float64(-0.16666666666666666 * (x ^ 3.0))), x); else tmp = copysign(log(Float64(x + hypot(1.0, x))), x); end return tmp end
function tmp_2 = code(x) tmp = 0.0; if (x <= -0.00105) tmp = sign(x) * abs(log((hypot(1.0, x) - x))); elseif (x <= 0.00078) tmp = sign(x) * abs((x + (-0.16666666666666666 * (x ^ 3.0)))); else tmp = sign(x) * abs(log((x + hypot(1.0, x)))); end tmp_2 = tmp; end
code[x_] := If[LessEqual[x, -0.00105], N[With[{TMP1 = Abs[N[Log[N[(N[Sqrt[1.0 ^ 2 + x ^ 2], $MachinePrecision] - x), $MachinePrecision]], $MachinePrecision]], TMP2 = Sign[x]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision], If[LessEqual[x, 0.00078], N[With[{TMP1 = Abs[N[(x + N[(-0.16666666666666666 * N[Power[x, 3.0], $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]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -0.00105:\\
\;\;\;\;\mathsf{copysign}\left(\log \left(\mathsf{hypot}\left(1, x\right) - x\right), x\right)\\
\mathbf{elif}\;x \leq 0.00078:\\
\;\;\;\;\mathsf{copysign}\left(x + -0.16666666666666666 \cdot {x}^{3}, x\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{copysign}\left(\log \left(x + \mathsf{hypot}\left(1, x\right)\right), x\right)\\
\end{array}
\end{array}
if x < -0.00104999999999999994Initial program 60.6%
+-commutative60.6%
hypot-1-def100.0%
Simplified100.0%
flip-+7.4%
div-sub7.4%
pow27.4%
add-sqr-sqrt0.0%
fabs-sqr0.0%
add-sqr-sqrt7.4%
pow27.4%
add-sqr-sqrt0.0%
fabs-sqr0.0%
add-sqr-sqrt1.7%
hypot-udef1.7%
hypot-udef1.7%
add-sqr-sqrt1.7%
metadata-eval1.7%
add-sqr-sqrt0.0%
fabs-sqr0.0%
add-sqr-sqrt8.9%
Applied egg-rr8.9%
unpow28.8%
div-sub10.0%
unpow210.0%
unpow210.0%
unpow210.0%
+-commutative10.0%
associate--r+59.3%
+-inverses100.0%
metadata-eval100.0%
metadata-eval100.0%
associate-/r*100.0%
neg-mul-1100.0%
sub-neg100.0%
+-commutative100.0%
distribute-neg-in100.0%
remove-double-neg100.0%
sub-neg100.0%
Simplified100.0%
add-sqr-sqrt0.0%
sqrt-unprod100.0%
pow1/2100.0%
log-rec100.0%
log-rec100.0%
sqr-neg100.0%
pow1100.0%
pow1100.0%
pow-prod-up100.0%
metadata-eval100.0%
Applied egg-rr100.0%
unpow1/2100.0%
unpow2100.0%
rem-sqrt-square100.0%
unpow1100.0%
sqr-pow99.4%
fabs-sqr99.4%
sqr-pow100.0%
unpow1100.0%
Simplified100.0%
if -0.00104999999999999994 < x < 7.79999999999999986e-4Initial program 7.3%
+-commutative7.3%
hypot-1-def7.3%
Simplified7.3%
flip-+7.3%
div-sub7.3%
pow27.3%
add-sqr-sqrt3.9%
fabs-sqr3.9%
add-sqr-sqrt7.3%
pow27.3%
add-sqr-sqrt3.9%
fabs-sqr3.9%
add-sqr-sqrt7.3%
hypot-udef7.3%
hypot-udef7.3%
add-sqr-sqrt7.3%
metadata-eval7.3%
add-sqr-sqrt3.9%
fabs-sqr3.9%
add-sqr-sqrt7.3%
Applied egg-rr7.3%
unpow27.3%
div-sub7.3%
unpow27.3%
unpow27.3%
unpow27.3%
+-commutative7.3%
associate--r+7.3%
+-inverses7.3%
metadata-eval7.3%
metadata-eval7.3%
associate-/r*7.3%
neg-mul-17.3%
sub-neg7.3%
+-commutative7.3%
distribute-neg-in7.3%
remove-double-neg7.3%
sub-neg7.3%
Simplified7.3%
Taylor expanded in x around 0 100.0%
if 7.79999999999999986e-4 < x Initial program 61.4%
+-commutative61.4%
hypot-1-def99.9%
Simplified99.9%
*-un-lft-identity99.9%
log-prod99.9%
metadata-eval99.9%
*-un-lft-identity99.9%
*-un-lft-identity99.9%
add-sqr-sqrt99.9%
fabs-sqr99.9%
add-sqr-sqrt99.9%
Applied egg-rr99.9%
+-lft-identity99.9%
Simplified99.9%
Final simplification100.0%
(FPCore (x)
:precision binary64
(if (<= x -0.00105)
(copysign (log (/ 1.0 (- (hypot 1.0 x) x))) x)
(if (<= x 0.00078)
(copysign (+ x (* -0.16666666666666666 (pow x 3.0))) x)
(copysign (log (+ x (hypot 1.0 x))) x))))
double code(double x) {
double tmp;
if (x <= -0.00105) {
tmp = copysign(log((1.0 / (hypot(1.0, x) - x))), x);
} else if (x <= 0.00078) {
tmp = copysign((x + (-0.16666666666666666 * pow(x, 3.0))), x);
} else {
tmp = copysign(log((x + hypot(1.0, x))), x);
}
return tmp;
}
public static double code(double x) {
double tmp;
if (x <= -0.00105) {
tmp = Math.copySign(Math.log((1.0 / (Math.hypot(1.0, x) - x))), x);
} else if (x <= 0.00078) {
tmp = Math.copySign((x + (-0.16666666666666666 * Math.pow(x, 3.0))), x);
} else {
tmp = Math.copySign(Math.log((x + Math.hypot(1.0, x))), x);
}
return tmp;
}
def code(x): tmp = 0 if x <= -0.00105: tmp = math.copysign(math.log((1.0 / (math.hypot(1.0, x) - x))), x) elif x <= 0.00078: tmp = math.copysign((x + (-0.16666666666666666 * math.pow(x, 3.0))), x) else: tmp = math.copysign(math.log((x + math.hypot(1.0, x))), x) return tmp
function code(x) tmp = 0.0 if (x <= -0.00105) tmp = copysign(log(Float64(1.0 / Float64(hypot(1.0, x) - x))), x); elseif (x <= 0.00078) tmp = copysign(Float64(x + Float64(-0.16666666666666666 * (x ^ 3.0))), x); else tmp = copysign(log(Float64(x + hypot(1.0, x))), x); end return tmp end
function tmp_2 = code(x) tmp = 0.0; if (x <= -0.00105) tmp = sign(x) * abs(log((1.0 / (hypot(1.0, x) - x)))); elseif (x <= 0.00078) tmp = sign(x) * abs((x + (-0.16666666666666666 * (x ^ 3.0)))); else tmp = sign(x) * abs(log((x + hypot(1.0, x)))); end tmp_2 = tmp; end
code[x_] := If[LessEqual[x, -0.00105], 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[x, 0.00078], N[With[{TMP1 = Abs[N[(x + N[(-0.16666666666666666 * N[Power[x, 3.0], $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]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -0.00105:\\
\;\;\;\;\mathsf{copysign}\left(\log \left(\frac{1}{\mathsf{hypot}\left(1, x\right) - x}\right), x\right)\\
\mathbf{elif}\;x \leq 0.00078:\\
\;\;\;\;\mathsf{copysign}\left(x + -0.16666666666666666 \cdot {x}^{3}, x\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{copysign}\left(\log \left(x + \mathsf{hypot}\left(1, x\right)\right), x\right)\\
\end{array}
\end{array}
if x < -0.00104999999999999994Initial program 60.6%
+-commutative60.6%
hypot-1-def100.0%
Simplified100.0%
flip-+7.4%
div-sub7.4%
pow27.4%
add-sqr-sqrt0.0%
fabs-sqr0.0%
add-sqr-sqrt7.4%
pow27.4%
add-sqr-sqrt0.0%
fabs-sqr0.0%
add-sqr-sqrt1.7%
hypot-udef1.7%
hypot-udef1.7%
add-sqr-sqrt1.7%
metadata-eval1.7%
add-sqr-sqrt0.0%
fabs-sqr0.0%
add-sqr-sqrt8.9%
Applied egg-rr8.9%
unpow28.8%
div-sub10.0%
unpow210.0%
unpow210.0%
unpow210.0%
+-commutative10.0%
associate--r+59.3%
+-inverses100.0%
metadata-eval100.0%
metadata-eval100.0%
associate-/r*100.0%
neg-mul-1100.0%
sub-neg100.0%
+-commutative100.0%
distribute-neg-in100.0%
remove-double-neg100.0%
sub-neg100.0%
Simplified100.0%
if -0.00104999999999999994 < x < 7.79999999999999986e-4Initial program 7.3%
+-commutative7.3%
hypot-1-def7.3%
Simplified7.3%
flip-+7.3%
div-sub7.3%
pow27.3%
add-sqr-sqrt3.9%
fabs-sqr3.9%
add-sqr-sqrt7.3%
pow27.3%
add-sqr-sqrt3.9%
fabs-sqr3.9%
add-sqr-sqrt7.3%
hypot-udef7.3%
hypot-udef7.3%
add-sqr-sqrt7.3%
metadata-eval7.3%
add-sqr-sqrt3.9%
fabs-sqr3.9%
add-sqr-sqrt7.3%
Applied egg-rr7.3%
unpow27.3%
div-sub7.3%
unpow27.3%
unpow27.3%
unpow27.3%
+-commutative7.3%
associate--r+7.3%
+-inverses7.3%
metadata-eval7.3%
metadata-eval7.3%
associate-/r*7.3%
neg-mul-17.3%
sub-neg7.3%
+-commutative7.3%
distribute-neg-in7.3%
remove-double-neg7.3%
sub-neg7.3%
Simplified7.3%
Taylor expanded in x around 0 100.0%
if 7.79999999999999986e-4 < x Initial program 61.4%
+-commutative61.4%
hypot-1-def99.9%
Simplified99.9%
*-un-lft-identity99.9%
log-prod99.9%
metadata-eval99.9%
*-un-lft-identity99.9%
*-un-lft-identity99.9%
add-sqr-sqrt99.9%
fabs-sqr99.9%
add-sqr-sqrt99.9%
Applied egg-rr99.9%
+-lft-identity99.9%
Simplified99.9%
Final simplification100.0%
(FPCore (x)
:precision binary64
(if (<= x -1.25)
(copysign (log (/ -0.5 x)) x)
(if (<= x 0.96)
(copysign (+ x (* -0.16666666666666666 (pow x 3.0))) x)
(copysign (log (+ x (+ x (/ 0.5 x)))) x))))
double code(double x) {
double tmp;
if (x <= -1.25) {
tmp = copysign(log((-0.5 / x)), x);
} else if (x <= 0.96) {
tmp = copysign((x + (-0.16666666666666666 * pow(x, 3.0))), x);
} else {
tmp = copysign(log((x + (x + (0.5 / x)))), x);
}
return tmp;
}
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.96) {
tmp = Math.copySign((x + (-0.16666666666666666 * Math.pow(x, 3.0))), x);
} else {
tmp = Math.copySign(Math.log((x + (x + (0.5 / x)))), x);
}
return tmp;
}
def code(x): tmp = 0 if x <= -1.25: tmp = math.copysign(math.log((-0.5 / x)), x) elif x <= 0.96: tmp = math.copysign((x + (-0.16666666666666666 * math.pow(x, 3.0))), x) else: tmp = math.copysign(math.log((x + (x + (0.5 / x)))), x) return tmp
function code(x) tmp = 0.0 if (x <= -1.25) tmp = copysign(log(Float64(-0.5 / x)), x); elseif (x <= 0.96) tmp = copysign(Float64(x + Float64(-0.16666666666666666 * (x ^ 3.0))), x); else tmp = copysign(log(Float64(x + Float64(x + Float64(0.5 / x)))), x); end return tmp end
function tmp_2 = code(x) tmp = 0.0; if (x <= -1.25) tmp = sign(x) * abs(log((-0.5 / x))); elseif (x <= 0.96) tmp = sign(x) * abs((x + (-0.16666666666666666 * (x ^ 3.0)))); else tmp = sign(x) * abs(log((x + (x + (0.5 / x))))); end tmp_2 = tmp; end
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.96], N[With[{TMP1 = Abs[N[(x + N[(-0.16666666666666666 * N[Power[x, 3.0], $MachinePrecision]), $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]]]
\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.96:\\
\;\;\;\;\mathsf{copysign}\left(x + -0.16666666666666666 \cdot {x}^{3}, x\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{copysign}\left(\log \left(x + \left(x + \frac{0.5}{x}\right)\right), x\right)\\
\end{array}
\end{array}
if x < -1.25Initial program 59.9%
+-commutative59.9%
hypot-1-def100.0%
Simplified100.0%
Taylor expanded in x around -inf 97.7%
associate--l+97.7%
unpow197.7%
sqr-pow0.0%
fabs-sqr0.0%
sqr-pow6.1%
unpow16.1%
associate-+r-96.3%
mul-1-neg96.3%
sub-neg96.3%
+-inverses96.3%
neg-sub096.3%
associate-*r/96.3%
metadata-eval96.3%
distribute-neg-frac96.3%
metadata-eval96.3%
Simplified96.3%
if -1.25 < x < 0.95999999999999996Initial program 8.7%
+-commutative8.7%
hypot-1-def8.7%
Simplified8.7%
flip-+8.7%
div-sub8.7%
pow28.7%
add-sqr-sqrt4.6%
fabs-sqr4.6%
add-sqr-sqrt8.7%
pow28.7%
add-sqr-sqrt4.6%
fabs-sqr4.6%
add-sqr-sqrt8.1%
hypot-udef8.1%
hypot-udef8.1%
add-sqr-sqrt8.1%
metadata-eval8.1%
add-sqr-sqrt4.5%
fabs-sqr4.5%
add-sqr-sqrt8.7%
Applied egg-rr8.7%
unpow28.7%
div-sub8.7%
unpow28.7%
unpow28.7%
unpow28.7%
+-commutative8.7%
associate--r+8.7%
+-inverses8.7%
metadata-eval8.7%
metadata-eval8.7%
associate-/r*8.7%
neg-mul-18.7%
sub-neg8.7%
+-commutative8.7%
distribute-neg-in8.7%
remove-double-neg8.7%
sub-neg8.7%
Simplified8.7%
Taylor expanded in x around 0 99.3%
if 0.95999999999999996 < x Initial program 61.0%
+-commutative61.0%
hypot-1-def100.0%
Simplified100.0%
Taylor expanded in x around inf 100.0%
associate-+r+100.0%
+-commutative100.0%
+-commutative100.0%
unpow1100.0%
sqr-pow100.0%
fabs-sqr100.0%
sqr-pow100.0%
unpow1100.0%
associate-*r/100.0%
metadata-eval100.0%
Simplified100.0%
Final simplification98.8%
(FPCore (x)
:precision binary64
(if (<= x -0.95)
(copysign (log (- (* x -2.0) (/ 0.5 x))) x)
(if (<= x 0.96)
(copysign (+ x (* -0.16666666666666666 (pow x 3.0))) x)
(copysign (log (+ x (+ x (/ 0.5 x)))) x))))
double code(double x) {
double tmp;
if (x <= -0.95) {
tmp = copysign(log(((x * -2.0) - (0.5 / x))), x);
} else if (x <= 0.96) {
tmp = copysign((x + (-0.16666666666666666 * pow(x, 3.0))), x);
} else {
tmp = copysign(log((x + (x + (0.5 / x)))), x);
}
return tmp;
}
public static double code(double x) {
double tmp;
if (x <= -0.95) {
tmp = Math.copySign(Math.log(((x * -2.0) - (0.5 / x))), x);
} else if (x <= 0.96) {
tmp = Math.copySign((x + (-0.16666666666666666 * Math.pow(x, 3.0))), x);
} else {
tmp = Math.copySign(Math.log((x + (x + (0.5 / x)))), x);
}
return tmp;
}
def code(x): tmp = 0 if x <= -0.95: tmp = math.copysign(math.log(((x * -2.0) - (0.5 / x))), x) elif x <= 0.96: tmp = math.copysign((x + (-0.16666666666666666 * math.pow(x, 3.0))), x) else: tmp = math.copysign(math.log((x + (x + (0.5 / x)))), x) return tmp
function code(x) tmp = 0.0 if (x <= -0.95) tmp = copysign(log(Float64(Float64(x * -2.0) - Float64(0.5 / x))), x); elseif (x <= 0.96) tmp = copysign(Float64(x + Float64(-0.16666666666666666 * (x ^ 3.0))), x); else tmp = copysign(log(Float64(x + Float64(x + Float64(0.5 / x)))), x); end return tmp end
function tmp_2 = code(x) tmp = 0.0; if (x <= -0.95) tmp = sign(x) * abs(log(((x * -2.0) - (0.5 / x)))); elseif (x <= 0.96) tmp = sign(x) * abs((x + (-0.16666666666666666 * (x ^ 3.0)))); else tmp = sign(x) * abs(log((x + (x + (0.5 / x))))); end tmp_2 = tmp; end
code[x_] := If[LessEqual[x, -0.95], N[With[{TMP1 = Abs[N[Log[N[(N[(x * -2.0), $MachinePrecision] - N[(0.5 / x), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], TMP2 = Sign[x]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision], If[LessEqual[x, 0.96], N[With[{TMP1 = Abs[N[(x + N[(-0.16666666666666666 * N[Power[x, 3.0], $MachinePrecision]), $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]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -0.95:\\
\;\;\;\;\mathsf{copysign}\left(\log \left(x \cdot -2 - \frac{0.5}{x}\right), x\right)\\
\mathbf{elif}\;x \leq 0.96:\\
\;\;\;\;\mathsf{copysign}\left(x + -0.16666666666666666 \cdot {x}^{3}, x\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{copysign}\left(\log \left(x + \left(x + \frac{0.5}{x}\right)\right), x\right)\\
\end{array}
\end{array}
if x < -0.94999999999999996Initial program 59.9%
+-commutative59.9%
hypot-1-def100.0%
Simplified100.0%
flip-+5.8%
div-sub5.8%
pow25.8%
add-sqr-sqrt0.0%
fabs-sqr0.0%
add-sqr-sqrt5.8%
pow25.8%
add-sqr-sqrt0.0%
fabs-sqr0.0%
add-sqr-sqrt1.2%
hypot-udef1.2%
hypot-udef1.2%
add-sqr-sqrt1.2%
metadata-eval1.2%
add-sqr-sqrt0.0%
fabs-sqr0.0%
add-sqr-sqrt7.4%
Applied egg-rr7.4%
unpow27.3%
div-sub8.4%
unpow28.5%
unpow28.4%
unpow28.5%
+-commutative8.5%
associate--r+58.6%
+-inverses100.0%
metadata-eval100.0%
metadata-eval100.0%
associate-/r*100.0%
neg-mul-1100.0%
sub-neg100.0%
+-commutative100.0%
distribute-neg-in100.0%
remove-double-neg100.0%
sub-neg100.0%
Simplified100.0%
add-sqr-sqrt0.0%
sqrt-unprod100.0%
pow1/2100.0%
log-rec100.0%
log-rec100.0%
sqr-neg100.0%
pow1100.0%
pow1100.0%
pow-prod-up100.0%
metadata-eval100.0%
Applied egg-rr100.0%
unpow1/2100.0%
unpow2100.0%
rem-sqrt-square100.0%
unpow1100.0%
sqr-pow99.5%
fabs-sqr99.5%
sqr-pow100.0%
unpow1100.0%
Simplified100.0%
Taylor expanded in x around -inf 97.7%
*-commutative97.7%
associate-*r/97.7%
metadata-eval97.7%
Simplified97.7%
if -0.94999999999999996 < x < 0.95999999999999996Initial program 8.7%
+-commutative8.7%
hypot-1-def8.7%
Simplified8.7%
flip-+8.7%
div-sub8.7%
pow28.7%
add-sqr-sqrt4.6%
fabs-sqr4.6%
add-sqr-sqrt8.7%
pow28.7%
add-sqr-sqrt4.6%
fabs-sqr4.6%
add-sqr-sqrt8.1%
hypot-udef8.1%
hypot-udef8.1%
add-sqr-sqrt8.1%
metadata-eval8.1%
add-sqr-sqrt4.5%
fabs-sqr4.5%
add-sqr-sqrt8.7%
Applied egg-rr8.7%
unpow28.7%
div-sub8.7%
unpow28.7%
unpow28.7%
unpow28.7%
+-commutative8.7%
associate--r+8.7%
+-inverses8.7%
metadata-eval8.7%
metadata-eval8.7%
associate-/r*8.7%
neg-mul-18.7%
sub-neg8.7%
+-commutative8.7%
distribute-neg-in8.7%
remove-double-neg8.7%
sub-neg8.7%
Simplified8.7%
Taylor expanded in x around 0 99.3%
if 0.95999999999999996 < x Initial program 61.0%
+-commutative61.0%
hypot-1-def100.0%
Simplified100.0%
Taylor expanded in x around inf 100.0%
associate-+r+100.0%
+-commutative100.0%
+-commutative100.0%
unpow1100.0%
sqr-pow100.0%
fabs-sqr100.0%
sqr-pow100.0%
unpow1100.0%
associate-*r/100.0%
metadata-eval100.0%
Simplified100.0%
Final simplification99.1%
(FPCore (x)
:precision binary64
(if (<= x -1.25)
(copysign (log (/ -0.5 x)) x)
(if (<= x 1.3)
(copysign (+ x (* -0.16666666666666666 (pow x 3.0))) x)
(copysign (log (+ x x)) x))))
double code(double x) {
double tmp;
if (x <= -1.25) {
tmp = copysign(log((-0.5 / x)), x);
} else if (x <= 1.3) {
tmp = copysign((x + (-0.16666666666666666 * pow(x, 3.0))), x);
} else {
tmp = copysign(log((x + x)), x);
}
return tmp;
}
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.3) {
tmp = Math.copySign((x + (-0.16666666666666666 * Math.pow(x, 3.0))), x);
} else {
tmp = Math.copySign(Math.log((x + x)), x);
}
return tmp;
}
def code(x): tmp = 0 if x <= -1.25: tmp = math.copysign(math.log((-0.5 / x)), x) elif x <= 1.3: tmp = math.copysign((x + (-0.16666666666666666 * math.pow(x, 3.0))), x) else: tmp = math.copysign(math.log((x + x)), x) return tmp
function code(x) tmp = 0.0 if (x <= -1.25) tmp = copysign(log(Float64(-0.5 / x)), x); elseif (x <= 1.3) tmp = copysign(Float64(x + Float64(-0.16666666666666666 * (x ^ 3.0))), x); else tmp = copysign(log(Float64(x + x)), x); end return tmp end
function tmp_2 = code(x) tmp = 0.0; if (x <= -1.25) tmp = sign(x) * abs(log((-0.5 / x))); elseif (x <= 1.3) tmp = sign(x) * abs((x + (-0.16666666666666666 * (x ^ 3.0)))); else tmp = sign(x) * abs(log((x + x))); end tmp_2 = tmp; end
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.3], N[With[{TMP1 = Abs[N[(x + N[(-0.16666666666666666 * N[Power[x, 3.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], 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]]]
\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.3:\\
\;\;\;\;\mathsf{copysign}\left(x + -0.16666666666666666 \cdot {x}^{3}, x\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{copysign}\left(\log \left(x + x\right), x\right)\\
\end{array}
\end{array}
if x < -1.25Initial program 59.9%
+-commutative59.9%
hypot-1-def100.0%
Simplified100.0%
Taylor expanded in x around -inf 97.7%
associate--l+97.7%
unpow197.7%
sqr-pow0.0%
fabs-sqr0.0%
sqr-pow6.1%
unpow16.1%
associate-+r-96.3%
mul-1-neg96.3%
sub-neg96.3%
+-inverses96.3%
neg-sub096.3%
associate-*r/96.3%
metadata-eval96.3%
distribute-neg-frac96.3%
metadata-eval96.3%
Simplified96.3%
if -1.25 < x < 1.30000000000000004Initial program 8.7%
+-commutative8.7%
hypot-1-def8.7%
Simplified8.7%
flip-+8.7%
div-sub8.7%
pow28.7%
add-sqr-sqrt4.6%
fabs-sqr4.6%
add-sqr-sqrt8.7%
pow28.7%
add-sqr-sqrt4.6%
fabs-sqr4.6%
add-sqr-sqrt8.1%
hypot-udef8.1%
hypot-udef8.1%
add-sqr-sqrt8.1%
metadata-eval8.1%
add-sqr-sqrt4.5%
fabs-sqr4.5%
add-sqr-sqrt8.7%
Applied egg-rr8.7%
unpow28.7%
div-sub8.7%
unpow28.7%
unpow28.7%
unpow28.7%
+-commutative8.7%
associate--r+8.7%
+-inverses8.7%
metadata-eval8.7%
metadata-eval8.7%
associate-/r*8.7%
neg-mul-18.7%
sub-neg8.7%
+-commutative8.7%
distribute-neg-in8.7%
remove-double-neg8.7%
sub-neg8.7%
Simplified8.7%
Taylor expanded in x around 0 99.3%
if 1.30000000000000004 < x Initial program 61.0%
+-commutative61.0%
hypot-1-def100.0%
Simplified100.0%
Taylor expanded in x around inf 99.3%
unpow199.3%
sqr-pow99.3%
fabs-sqr99.3%
sqr-pow99.3%
unpow199.3%
Simplified99.3%
Final simplification98.6%
(FPCore (x) :precision binary64 (if (<= x -3.15) (copysign (log (- x)) x) (if (<= x 1.25) (copysign x x) (copysign (log (+ x x)) x))))
double code(double x) {
double tmp;
if (x <= -3.15) {
tmp = copysign(log(-x), x);
} else if (x <= 1.25) {
tmp = copysign(x, x);
} else {
tmp = copysign(log((x + x)), x);
}
return tmp;
}
public static double code(double x) {
double tmp;
if (x <= -3.15) {
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;
}
def code(x): tmp = 0 if x <= -3.15: 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
function code(x) tmp = 0.0 if (x <= -3.15) 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
function tmp_2 = code(x) tmp = 0.0; if (x <= -3.15) 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
code[x_] := If[LessEqual[x, -3.15], 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]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -3.15:\\
\;\;\;\;\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}
if x < -3.14999999999999991Initial program 59.2%
+-commutative59.2%
hypot-1-def100.0%
Simplified100.0%
Taylor expanded in x around -inf 31.1%
mul-1-neg31.1%
Simplified31.1%
if -3.14999999999999991 < x < 1.25Initial program 9.4%
+-commutative9.4%
hypot-1-def9.4%
Simplified9.4%
Taylor expanded in x around 0 7.5%
unpow17.5%
sqr-pow3.8%
fabs-sqr3.8%
sqr-pow7.3%
unpow17.3%
Simplified7.3%
Taylor expanded in x around 0 98.2%
if 1.25 < x Initial program 61.0%
+-commutative61.0%
hypot-1-def100.0%
Simplified100.0%
Taylor expanded in x around inf 99.3%
unpow199.3%
sqr-pow99.3%
fabs-sqr99.3%
sqr-pow99.3%
unpow199.3%
Simplified99.3%
Final simplification83.6%
(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))))
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;
}
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;
}
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
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
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
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]]]
\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}
if x < -1.25Initial program 59.9%
+-commutative59.9%
hypot-1-def100.0%
Simplified100.0%
Taylor expanded in x around -inf 97.7%
associate--l+97.7%
unpow197.7%
sqr-pow0.0%
fabs-sqr0.0%
sqr-pow6.1%
unpow16.1%
associate-+r-96.3%
mul-1-neg96.3%
sub-neg96.3%
+-inverses96.3%
neg-sub096.3%
associate-*r/96.3%
metadata-eval96.3%
distribute-neg-frac96.3%
metadata-eval96.3%
Simplified96.3%
if -1.25 < x < 1.25Initial program 8.7%
+-commutative8.7%
hypot-1-def8.7%
Simplified8.7%
Taylor expanded in x around 0 7.4%
unpow17.4%
sqr-pow3.9%
fabs-sqr3.9%
sqr-pow7.4%
unpow17.4%
Simplified7.4%
Taylor expanded in x around 0 98.8%
if 1.25 < x Initial program 61.0%
+-commutative61.0%
hypot-1-def100.0%
Simplified100.0%
Taylor expanded in x around inf 99.3%
unpow199.3%
sqr-pow99.3%
fabs-sqr99.3%
sqr-pow99.3%
unpow199.3%
Simplified99.3%
Final simplification98.4%
(FPCore (x) :precision binary64 (if (<= x -0.5) (copysign (log (- x)) x) (copysign (log1p x) x)))
double code(double x) {
double tmp;
if (x <= -0.5) {
tmp = copysign(log(-x), x);
} else {
tmp = copysign(log1p(x), x);
}
return tmp;
}
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;
}
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
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
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]]
\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}
if x < -0.5Initial program 59.9%
+-commutative59.9%
hypot-1-def100.0%
Simplified100.0%
Taylor expanded in x around -inf 30.9%
mul-1-neg30.9%
Simplified30.9%
if -0.5 < x Initial program 27.7%
+-commutative27.7%
hypot-1-def41.9%
Simplified41.9%
Taylor expanded in x around 0 16.0%
log1p-def73.4%
unpow173.4%
sqr-pow43.1%
fabs-sqr43.1%
sqr-pow73.4%
unpow173.4%
Simplified73.4%
Final simplification63.8%
(FPCore (x) :precision binary64 (if (<= x 1.6) (copysign x x) (copysign (log1p x) x)))
double code(double x) {
double tmp;
if (x <= 1.6) {
tmp = copysign(x, x);
} else {
tmp = copysign(log1p(x), x);
}
return tmp;
}
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;
}
def code(x): tmp = 0 if x <= 1.6: tmp = math.copysign(x, x) else: tmp = math.copysign(math.log1p(x), x) return tmp
function code(x) tmp = 0.0 if (x <= 1.6) tmp = copysign(x, x); else tmp = copysign(log1p(x), x); end return tmp end
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]]
\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}
if x < 1.6000000000000001Initial program 24.8%
+-commutative24.8%
hypot-1-def37.5%
Simplified37.5%
Taylor expanded in x around 0 14.8%
unpow114.8%
sqr-pow2.6%
fabs-sqr2.6%
sqr-pow5.1%
unpow15.1%
Simplified5.1%
Taylor expanded in x around 0 69.5%
if 1.6000000000000001 < x Initial program 61.0%
+-commutative61.0%
hypot-1-def100.0%
Simplified100.0%
Taylor expanded in x around 0 31.2%
log1p-def31.2%
unpow131.2%
sqr-pow31.2%
fabs-sqr31.2%
sqr-pow31.2%
unpow131.2%
Simplified31.2%
Final simplification58.8%
(FPCore (x) :precision binary64 (copysign x x))
double code(double x) {
return copysign(x, x);
}
public static double code(double x) {
return Math.copySign(x, x);
}
def code(x): return math.copysign(x, x)
function code(x) return copysign(x, x) end
function tmp = code(x) tmp = sign(x) * abs(x); end
code[x_] := N[With[{TMP1 = Abs[x], TMP2 = Sign[x]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
\begin{array}{l}
\\
\mathsf{copysign}\left(x, x\right)
\end{array}
Initial program 35.0%
+-commutative35.0%
hypot-1-def55.1%
Simplified55.1%
Taylor expanded in x around 0 19.4%
unpow119.4%
sqr-pow10.7%
fabs-sqr10.7%
sqr-pow12.4%
unpow112.4%
Simplified12.4%
Taylor expanded in x around 0 51.6%
Final simplification51.6%
(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)))
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);
}
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);
}
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)
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
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]]
\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}
herbie shell --seed 2023196
(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))