
(FPCore (wj x) :precision binary64 (let* ((t_0 (* wj (exp wj)))) (- wj (/ (- t_0 x) (+ (exp wj) t_0)))))
double code(double wj, double x) {
double t_0 = wj * exp(wj);
return wj - ((t_0 - x) / (exp(wj) + t_0));
}
real(8) function code(wj, x)
real(8), intent (in) :: wj
real(8), intent (in) :: x
real(8) :: t_0
t_0 = wj * exp(wj)
code = wj - ((t_0 - x) / (exp(wj) + t_0))
end function
public static double code(double wj, double x) {
double t_0 = wj * Math.exp(wj);
return wj - ((t_0 - x) / (Math.exp(wj) + t_0));
}
def code(wj, x): t_0 = wj * math.exp(wj) return wj - ((t_0 - x) / (math.exp(wj) + t_0))
function code(wj, x) t_0 = Float64(wj * exp(wj)) return Float64(wj - Float64(Float64(t_0 - x) / Float64(exp(wj) + t_0))) end
function tmp = code(wj, x) t_0 = wj * exp(wj); tmp = wj - ((t_0 - x) / (exp(wj) + t_0)); end
code[wj_, x_] := Block[{t$95$0 = N[(wj * N[Exp[wj], $MachinePrecision]), $MachinePrecision]}, N[(wj - N[(N[(t$95$0 - x), $MachinePrecision] / N[(N[Exp[wj], $MachinePrecision] + t$95$0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := wj \cdot e^{wj}\\
wj - \frac{t\_0 - x}{e^{wj} + t\_0}
\end{array}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 20 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (wj x) :precision binary64 (let* ((t_0 (* wj (exp wj)))) (- wj (/ (- t_0 x) (+ (exp wj) t_0)))))
double code(double wj, double x) {
double t_0 = wj * exp(wj);
return wj - ((t_0 - x) / (exp(wj) + t_0));
}
real(8) function code(wj, x)
real(8), intent (in) :: wj
real(8), intent (in) :: x
real(8) :: t_0
t_0 = wj * exp(wj)
code = wj - ((t_0 - x) / (exp(wj) + t_0))
end function
public static double code(double wj, double x) {
double t_0 = wj * Math.exp(wj);
return wj - ((t_0 - x) / (Math.exp(wj) + t_0));
}
def code(wj, x): t_0 = wj * math.exp(wj) return wj - ((t_0 - x) / (math.exp(wj) + t_0))
function code(wj, x) t_0 = Float64(wj * exp(wj)) return Float64(wj - Float64(Float64(t_0 - x) / Float64(exp(wj) + t_0))) end
function tmp = code(wj, x) t_0 = wj * exp(wj); tmp = wj - ((t_0 - x) / (exp(wj) + t_0)); end
code[wj_, x_] := Block[{t$95$0 = N[(wj * N[Exp[wj], $MachinePrecision]), $MachinePrecision]}, N[(wj - N[(N[(t$95$0 - x), $MachinePrecision] / N[(N[Exp[wj], $MachinePrecision] + t$95$0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := wj \cdot e^{wj}\\
wj - \frac{t\_0 - x}{e^{wj} + t\_0}
\end{array}
\end{array}
(FPCore (wj x)
:precision binary64
(if (<= wj -6.8e-6)
(*
x
(+
(- (/ wj x) (/ (/ 1.0 (exp wj)) (- -1.0 wj)))
(* (/ wj x) (/ (+ 1.0 (* wj (+ wj -1.0))) (- -1.0 (* wj (* wj wj)))))))
(if (<= wj 1.3e-8)
(+
x
(*
wj
(+
(* x -2.0)
(*
wj
(+
(+
1.0
(*
wj
(- (- -1.0 (* x 0.6666666666666666)) (+ (* x -3.0) (* x 5.0)))))
(* x 2.5))))))
(+ wj (* (/ 1.0 (- -1.0 wj)) (- wj (/ x (exp wj))))))))
double code(double wj, double x) {
double tmp;
if (wj <= -6.8e-6) {
tmp = x * (((wj / x) - ((1.0 / exp(wj)) / (-1.0 - wj))) + ((wj / x) * ((1.0 + (wj * (wj + -1.0))) / (-1.0 - (wj * (wj * wj))))));
} else if (wj <= 1.3e-8) {
tmp = x + (wj * ((x * -2.0) + (wj * ((1.0 + (wj * ((-1.0 - (x * 0.6666666666666666)) - ((x * -3.0) + (x * 5.0))))) + (x * 2.5)))));
} else {
tmp = wj + ((1.0 / (-1.0 - wj)) * (wj - (x / exp(wj))));
}
return tmp;
}
real(8) function code(wj, x)
real(8), intent (in) :: wj
real(8), intent (in) :: x
real(8) :: tmp
if (wj <= (-6.8d-6)) then
tmp = x * (((wj / x) - ((1.0d0 / exp(wj)) / ((-1.0d0) - wj))) + ((wj / x) * ((1.0d0 + (wj * (wj + (-1.0d0)))) / ((-1.0d0) - (wj * (wj * wj))))))
else if (wj <= 1.3d-8) then
tmp = x + (wj * ((x * (-2.0d0)) + (wj * ((1.0d0 + (wj * (((-1.0d0) - (x * 0.6666666666666666d0)) - ((x * (-3.0d0)) + (x * 5.0d0))))) + (x * 2.5d0)))))
else
tmp = wj + ((1.0d0 / ((-1.0d0) - wj)) * (wj - (x / exp(wj))))
end if
code = tmp
end function
public static double code(double wj, double x) {
double tmp;
if (wj <= -6.8e-6) {
tmp = x * (((wj / x) - ((1.0 / Math.exp(wj)) / (-1.0 - wj))) + ((wj / x) * ((1.0 + (wj * (wj + -1.0))) / (-1.0 - (wj * (wj * wj))))));
} else if (wj <= 1.3e-8) {
tmp = x + (wj * ((x * -2.0) + (wj * ((1.0 + (wj * ((-1.0 - (x * 0.6666666666666666)) - ((x * -3.0) + (x * 5.0))))) + (x * 2.5)))));
} else {
tmp = wj + ((1.0 / (-1.0 - wj)) * (wj - (x / Math.exp(wj))));
}
return tmp;
}
def code(wj, x): tmp = 0 if wj <= -6.8e-6: tmp = x * (((wj / x) - ((1.0 / math.exp(wj)) / (-1.0 - wj))) + ((wj / x) * ((1.0 + (wj * (wj + -1.0))) / (-1.0 - (wj * (wj * wj)))))) elif wj <= 1.3e-8: tmp = x + (wj * ((x * -2.0) + (wj * ((1.0 + (wj * ((-1.0 - (x * 0.6666666666666666)) - ((x * -3.0) + (x * 5.0))))) + (x * 2.5))))) else: tmp = wj + ((1.0 / (-1.0 - wj)) * (wj - (x / math.exp(wj)))) return tmp
function code(wj, x) tmp = 0.0 if (wj <= -6.8e-6) tmp = Float64(x * Float64(Float64(Float64(wj / x) - Float64(Float64(1.0 / exp(wj)) / Float64(-1.0 - wj))) + Float64(Float64(wj / x) * Float64(Float64(1.0 + Float64(wj * Float64(wj + -1.0))) / Float64(-1.0 - Float64(wj * Float64(wj * wj))))))); elseif (wj <= 1.3e-8) tmp = Float64(x + Float64(wj * Float64(Float64(x * -2.0) + Float64(wj * Float64(Float64(1.0 + Float64(wj * Float64(Float64(-1.0 - Float64(x * 0.6666666666666666)) - Float64(Float64(x * -3.0) + Float64(x * 5.0))))) + Float64(x * 2.5)))))); else tmp = Float64(wj + Float64(Float64(1.0 / Float64(-1.0 - wj)) * Float64(wj - Float64(x / exp(wj))))); end return tmp end
function tmp_2 = code(wj, x) tmp = 0.0; if (wj <= -6.8e-6) tmp = x * (((wj / x) - ((1.0 / exp(wj)) / (-1.0 - wj))) + ((wj / x) * ((1.0 + (wj * (wj + -1.0))) / (-1.0 - (wj * (wj * wj)))))); elseif (wj <= 1.3e-8) tmp = x + (wj * ((x * -2.0) + (wj * ((1.0 + (wj * ((-1.0 - (x * 0.6666666666666666)) - ((x * -3.0) + (x * 5.0))))) + (x * 2.5))))); else tmp = wj + ((1.0 / (-1.0 - wj)) * (wj - (x / exp(wj)))); end tmp_2 = tmp; end
code[wj_, x_] := If[LessEqual[wj, -6.8e-6], N[(x * N[(N[(N[(wj / x), $MachinePrecision] - N[(N[(1.0 / N[Exp[wj], $MachinePrecision]), $MachinePrecision] / N[(-1.0 - wj), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(wj / x), $MachinePrecision] * N[(N[(1.0 + N[(wj * N[(wj + -1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(-1.0 - N[(wj * N[(wj * wj), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[wj, 1.3e-8], N[(x + N[(wj * N[(N[(x * -2.0), $MachinePrecision] + N[(wj * N[(N[(1.0 + N[(wj * N[(N[(-1.0 - N[(x * 0.6666666666666666), $MachinePrecision]), $MachinePrecision] - N[(N[(x * -3.0), $MachinePrecision] + N[(x * 5.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(x * 2.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(wj + N[(N[(1.0 / N[(-1.0 - wj), $MachinePrecision]), $MachinePrecision] * N[(wj - N[(x / N[Exp[wj], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;wj \leq -6.8 \cdot 10^{-6}:\\
\;\;\;\;x \cdot \left(\left(\frac{wj}{x} - \frac{\frac{1}{e^{wj}}}{-1 - wj}\right) + \frac{wj}{x} \cdot \frac{1 + wj \cdot \left(wj + -1\right)}{-1 - wj \cdot \left(wj \cdot wj\right)}\right)\\
\mathbf{elif}\;wj \leq 1.3 \cdot 10^{-8}:\\
\;\;\;\;x + wj \cdot \left(x \cdot -2 + wj \cdot \left(\left(1 + wj \cdot \left(\left(-1 - x \cdot 0.6666666666666666\right) - \left(x \cdot -3 + x \cdot 5\right)\right)\right) + x \cdot 2.5\right)\right)\\
\mathbf{else}:\\
\;\;\;\;wj + \frac{1}{-1 - wj} \cdot \left(wj - \frac{x}{e^{wj}}\right)\\
\end{array}
\end{array}
if wj < -6.80000000000000012e-6Initial program 51.4%
sub-negN/A
+-lowering-+.f64N/A
distribute-rgt1-inN/A
associate-/l/N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
div-subN/A
associate-/l*N/A
*-inversesN/A
*-rgt-identityN/A
--lowering--.f64N/A
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
+-commutativeN/A
*-inversesN/A
distribute-neg-inN/A
Simplified96.1%
+-commutativeN/A
div-subN/A
associate-+l-N/A
flip3--N/A
associate-/r/N/A
fmm-defN/A
fma-lowering-fma.f64N/A
Applied egg-rr96.1%
Taylor expanded in x around inf
*-lowering-*.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
Simplified96.8%
if -6.80000000000000012e-6 < wj < 1.3000000000000001e-8Initial program 79.4%
Taylor expanded in wj around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
cancel-sign-sub-invN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
Simplified100.0%
if 1.3000000000000001e-8 < wj Initial program 49.7%
sub-negN/A
+-lowering-+.f64N/A
distribute-rgt1-inN/A
associate-/l/N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
div-subN/A
associate-/l*N/A
*-inversesN/A
*-rgt-identityN/A
--lowering--.f64N/A
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
+-commutativeN/A
*-inversesN/A
distribute-neg-inN/A
Simplified99.7%
div-invN/A
*-commutativeN/A
*-lowering-*.f64N/A
metadata-evalN/A
/-lowering-/.f64N/A
metadata-evalN/A
--lowering--.f64N/A
--lowering--.f64N/A
/-lowering-/.f64N/A
exp-lowering-exp.f64100.0%
Applied egg-rr100.0%
Final simplification99.8%
(FPCore (wj x)
:precision binary64
(let* ((t_0 (- wj (/ x (exp wj)))))
(if (<= wj -6.5e-6)
(+ wj (/ t_0 (- -1.0 wj)))
(if (<= wj 1.3e-8)
(+
x
(*
wj
(+
(* x -2.0)
(*
wj
(+
(+
1.0
(*
wj
(- (- -1.0 (* x 0.6666666666666666)) (+ (* x -3.0) (* x 5.0)))))
(* x 2.5))))))
(+ wj (* (/ 1.0 (- -1.0 wj)) t_0))))))
double code(double wj, double x) {
double t_0 = wj - (x / exp(wj));
double tmp;
if (wj <= -6.5e-6) {
tmp = wj + (t_0 / (-1.0 - wj));
} else if (wj <= 1.3e-8) {
tmp = x + (wj * ((x * -2.0) + (wj * ((1.0 + (wj * ((-1.0 - (x * 0.6666666666666666)) - ((x * -3.0) + (x * 5.0))))) + (x * 2.5)))));
} else {
tmp = wj + ((1.0 / (-1.0 - wj)) * t_0);
}
return tmp;
}
real(8) function code(wj, x)
real(8), intent (in) :: wj
real(8), intent (in) :: x
real(8) :: t_0
real(8) :: tmp
t_0 = wj - (x / exp(wj))
if (wj <= (-6.5d-6)) then
tmp = wj + (t_0 / ((-1.0d0) - wj))
else if (wj <= 1.3d-8) then
tmp = x + (wj * ((x * (-2.0d0)) + (wj * ((1.0d0 + (wj * (((-1.0d0) - (x * 0.6666666666666666d0)) - ((x * (-3.0d0)) + (x * 5.0d0))))) + (x * 2.5d0)))))
else
tmp = wj + ((1.0d0 / ((-1.0d0) - wj)) * t_0)
end if
code = tmp
end function
public static double code(double wj, double x) {
double t_0 = wj - (x / Math.exp(wj));
double tmp;
if (wj <= -6.5e-6) {
tmp = wj + (t_0 / (-1.0 - wj));
} else if (wj <= 1.3e-8) {
tmp = x + (wj * ((x * -2.0) + (wj * ((1.0 + (wj * ((-1.0 - (x * 0.6666666666666666)) - ((x * -3.0) + (x * 5.0))))) + (x * 2.5)))));
} else {
tmp = wj + ((1.0 / (-1.0 - wj)) * t_0);
}
return tmp;
}
def code(wj, x): t_0 = wj - (x / math.exp(wj)) tmp = 0 if wj <= -6.5e-6: tmp = wj + (t_0 / (-1.0 - wj)) elif wj <= 1.3e-8: tmp = x + (wj * ((x * -2.0) + (wj * ((1.0 + (wj * ((-1.0 - (x * 0.6666666666666666)) - ((x * -3.0) + (x * 5.0))))) + (x * 2.5))))) else: tmp = wj + ((1.0 / (-1.0 - wj)) * t_0) return tmp
function code(wj, x) t_0 = Float64(wj - Float64(x / exp(wj))) tmp = 0.0 if (wj <= -6.5e-6) tmp = Float64(wj + Float64(t_0 / Float64(-1.0 - wj))); elseif (wj <= 1.3e-8) tmp = Float64(x + Float64(wj * Float64(Float64(x * -2.0) + Float64(wj * Float64(Float64(1.0 + Float64(wj * Float64(Float64(-1.0 - Float64(x * 0.6666666666666666)) - Float64(Float64(x * -3.0) + Float64(x * 5.0))))) + Float64(x * 2.5)))))); else tmp = Float64(wj + Float64(Float64(1.0 / Float64(-1.0 - wj)) * t_0)); end return tmp end
function tmp_2 = code(wj, x) t_0 = wj - (x / exp(wj)); tmp = 0.0; if (wj <= -6.5e-6) tmp = wj + (t_0 / (-1.0 - wj)); elseif (wj <= 1.3e-8) tmp = x + (wj * ((x * -2.0) + (wj * ((1.0 + (wj * ((-1.0 - (x * 0.6666666666666666)) - ((x * -3.0) + (x * 5.0))))) + (x * 2.5))))); else tmp = wj + ((1.0 / (-1.0 - wj)) * t_0); end tmp_2 = tmp; end
code[wj_, x_] := Block[{t$95$0 = N[(wj - N[(x / N[Exp[wj], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[wj, -6.5e-6], N[(wj + N[(t$95$0 / N[(-1.0 - wj), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[wj, 1.3e-8], N[(x + N[(wj * N[(N[(x * -2.0), $MachinePrecision] + N[(wj * N[(N[(1.0 + N[(wj * N[(N[(-1.0 - N[(x * 0.6666666666666666), $MachinePrecision]), $MachinePrecision] - N[(N[(x * -3.0), $MachinePrecision] + N[(x * 5.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(x * 2.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(wj + N[(N[(1.0 / N[(-1.0 - wj), $MachinePrecision]), $MachinePrecision] * t$95$0), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := wj - \frac{x}{e^{wj}}\\
\mathbf{if}\;wj \leq -6.5 \cdot 10^{-6}:\\
\;\;\;\;wj + \frac{t\_0}{-1 - wj}\\
\mathbf{elif}\;wj \leq 1.3 \cdot 10^{-8}:\\
\;\;\;\;x + wj \cdot \left(x \cdot -2 + wj \cdot \left(\left(1 + wj \cdot \left(\left(-1 - x \cdot 0.6666666666666666\right) - \left(x \cdot -3 + x \cdot 5\right)\right)\right) + x \cdot 2.5\right)\right)\\
\mathbf{else}:\\
\;\;\;\;wj + \frac{1}{-1 - wj} \cdot t\_0\\
\end{array}
\end{array}
if wj < -6.4999999999999996e-6Initial program 51.4%
sub-negN/A
+-lowering-+.f64N/A
distribute-rgt1-inN/A
associate-/l/N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
div-subN/A
associate-/l*N/A
*-inversesN/A
*-rgt-identityN/A
--lowering--.f64N/A
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
+-commutativeN/A
*-inversesN/A
distribute-neg-inN/A
Simplified96.1%
if -6.4999999999999996e-6 < wj < 1.3000000000000001e-8Initial program 79.4%
Taylor expanded in wj around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
cancel-sign-sub-invN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
Simplified100.0%
if 1.3000000000000001e-8 < wj Initial program 49.7%
sub-negN/A
+-lowering-+.f64N/A
distribute-rgt1-inN/A
associate-/l/N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
div-subN/A
associate-/l*N/A
*-inversesN/A
*-rgt-identityN/A
--lowering--.f64N/A
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
+-commutativeN/A
*-inversesN/A
distribute-neg-inN/A
Simplified99.7%
div-invN/A
*-commutativeN/A
*-lowering-*.f64N/A
metadata-evalN/A
/-lowering-/.f64N/A
metadata-evalN/A
--lowering--.f64N/A
--lowering--.f64N/A
/-lowering-/.f64N/A
exp-lowering-exp.f64100.0%
Applied egg-rr100.0%
Final simplification99.8%
(FPCore (wj x)
:precision binary64
(let* ((t_0 (+ wj (/ (- wj (/ x (exp wj))) (- -1.0 wj)))))
(if (<= wj -6.5e-6)
t_0
(if (<= wj 1.3e-8)
(+
x
(*
wj
(+
(* x -2.0)
(*
wj
(+
(+
1.0
(*
wj
(- (- -1.0 (* x 0.6666666666666666)) (+ (* x -3.0) (* x 5.0)))))
(* x 2.5))))))
t_0))))
double code(double wj, double x) {
double t_0 = wj + ((wj - (x / exp(wj))) / (-1.0 - wj));
double tmp;
if (wj <= -6.5e-6) {
tmp = t_0;
} else if (wj <= 1.3e-8) {
tmp = x + (wj * ((x * -2.0) + (wj * ((1.0 + (wj * ((-1.0 - (x * 0.6666666666666666)) - ((x * -3.0) + (x * 5.0))))) + (x * 2.5)))));
} else {
tmp = t_0;
}
return tmp;
}
real(8) function code(wj, x)
real(8), intent (in) :: wj
real(8), intent (in) :: x
real(8) :: t_0
real(8) :: tmp
t_0 = wj + ((wj - (x / exp(wj))) / ((-1.0d0) - wj))
if (wj <= (-6.5d-6)) then
tmp = t_0
else if (wj <= 1.3d-8) then
tmp = x + (wj * ((x * (-2.0d0)) + (wj * ((1.0d0 + (wj * (((-1.0d0) - (x * 0.6666666666666666d0)) - ((x * (-3.0d0)) + (x * 5.0d0))))) + (x * 2.5d0)))))
else
tmp = t_0
end if
code = tmp
end function
public static double code(double wj, double x) {
double t_0 = wj + ((wj - (x / Math.exp(wj))) / (-1.0 - wj));
double tmp;
if (wj <= -6.5e-6) {
tmp = t_0;
} else if (wj <= 1.3e-8) {
tmp = x + (wj * ((x * -2.0) + (wj * ((1.0 + (wj * ((-1.0 - (x * 0.6666666666666666)) - ((x * -3.0) + (x * 5.0))))) + (x * 2.5)))));
} else {
tmp = t_0;
}
return tmp;
}
def code(wj, x): t_0 = wj + ((wj - (x / math.exp(wj))) / (-1.0 - wj)) tmp = 0 if wj <= -6.5e-6: tmp = t_0 elif wj <= 1.3e-8: tmp = x + (wj * ((x * -2.0) + (wj * ((1.0 + (wj * ((-1.0 - (x * 0.6666666666666666)) - ((x * -3.0) + (x * 5.0))))) + (x * 2.5))))) else: tmp = t_0 return tmp
function code(wj, x) t_0 = Float64(wj + Float64(Float64(wj - Float64(x / exp(wj))) / Float64(-1.0 - wj))) tmp = 0.0 if (wj <= -6.5e-6) tmp = t_0; elseif (wj <= 1.3e-8) tmp = Float64(x + Float64(wj * Float64(Float64(x * -2.0) + Float64(wj * Float64(Float64(1.0 + Float64(wj * Float64(Float64(-1.0 - Float64(x * 0.6666666666666666)) - Float64(Float64(x * -3.0) + Float64(x * 5.0))))) + Float64(x * 2.5)))))); else tmp = t_0; end return tmp end
function tmp_2 = code(wj, x) t_0 = wj + ((wj - (x / exp(wj))) / (-1.0 - wj)); tmp = 0.0; if (wj <= -6.5e-6) tmp = t_0; elseif (wj <= 1.3e-8) tmp = x + (wj * ((x * -2.0) + (wj * ((1.0 + (wj * ((-1.0 - (x * 0.6666666666666666)) - ((x * -3.0) + (x * 5.0))))) + (x * 2.5))))); else tmp = t_0; end tmp_2 = tmp; end
code[wj_, x_] := Block[{t$95$0 = N[(wj + N[(N[(wj - N[(x / N[Exp[wj], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(-1.0 - wj), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[wj, -6.5e-6], t$95$0, If[LessEqual[wj, 1.3e-8], N[(x + N[(wj * N[(N[(x * -2.0), $MachinePrecision] + N[(wj * N[(N[(1.0 + N[(wj * N[(N[(-1.0 - N[(x * 0.6666666666666666), $MachinePrecision]), $MachinePrecision] - N[(N[(x * -3.0), $MachinePrecision] + N[(x * 5.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(x * 2.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$0]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := wj + \frac{wj - \frac{x}{e^{wj}}}{-1 - wj}\\
\mathbf{if}\;wj \leq -6.5 \cdot 10^{-6}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;wj \leq 1.3 \cdot 10^{-8}:\\
\;\;\;\;x + wj \cdot \left(x \cdot -2 + wj \cdot \left(\left(1 + wj \cdot \left(\left(-1 - x \cdot 0.6666666666666666\right) - \left(x \cdot -3 + x \cdot 5\right)\right)\right) + x \cdot 2.5\right)\right)\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if wj < -6.4999999999999996e-6 or 1.3000000000000001e-8 < wj Initial program 50.7%
sub-negN/A
+-lowering-+.f64N/A
distribute-rgt1-inN/A
associate-/l/N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
div-subN/A
associate-/l*N/A
*-inversesN/A
*-rgt-identityN/A
--lowering--.f64N/A
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
+-commutativeN/A
*-inversesN/A
distribute-neg-inN/A
Simplified97.5%
if -6.4999999999999996e-6 < wj < 1.3000000000000001e-8Initial program 79.4%
Taylor expanded in wj around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
cancel-sign-sub-invN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
Simplified100.0%
Final simplification99.8%
(FPCore (wj x)
:precision binary64
(if (<= wj 1.3e-8)
(+
x
(*
wj
(+
(* x -2.0)
(*
wj
(+
(+
1.0
(* wj (- (- -1.0 (* x 0.6666666666666666)) (+ (* x -3.0) (* x 5.0)))))
(* x 2.5))))))
(+
wj
(/
(+
wj
(/ x (+ -1.0 (* wj (- -1.0 (* wj (+ 0.5 (* wj 0.16666666666666666))))))))
(- -1.0 wj)))))
double code(double wj, double x) {
double tmp;
if (wj <= 1.3e-8) {
tmp = x + (wj * ((x * -2.0) + (wj * ((1.0 + (wj * ((-1.0 - (x * 0.6666666666666666)) - ((x * -3.0) + (x * 5.0))))) + (x * 2.5)))));
} else {
tmp = wj + ((wj + (x / (-1.0 + (wj * (-1.0 - (wj * (0.5 + (wj * 0.16666666666666666)))))))) / (-1.0 - wj));
}
return tmp;
}
real(8) function code(wj, x)
real(8), intent (in) :: wj
real(8), intent (in) :: x
real(8) :: tmp
if (wj <= 1.3d-8) then
tmp = x + (wj * ((x * (-2.0d0)) + (wj * ((1.0d0 + (wj * (((-1.0d0) - (x * 0.6666666666666666d0)) - ((x * (-3.0d0)) + (x * 5.0d0))))) + (x * 2.5d0)))))
else
tmp = wj + ((wj + (x / ((-1.0d0) + (wj * ((-1.0d0) - (wj * (0.5d0 + (wj * 0.16666666666666666d0)))))))) / ((-1.0d0) - wj))
end if
code = tmp
end function
public static double code(double wj, double x) {
double tmp;
if (wj <= 1.3e-8) {
tmp = x + (wj * ((x * -2.0) + (wj * ((1.0 + (wj * ((-1.0 - (x * 0.6666666666666666)) - ((x * -3.0) + (x * 5.0))))) + (x * 2.5)))));
} else {
tmp = wj + ((wj + (x / (-1.0 + (wj * (-1.0 - (wj * (0.5 + (wj * 0.16666666666666666)))))))) / (-1.0 - wj));
}
return tmp;
}
def code(wj, x): tmp = 0 if wj <= 1.3e-8: tmp = x + (wj * ((x * -2.0) + (wj * ((1.0 + (wj * ((-1.0 - (x * 0.6666666666666666)) - ((x * -3.0) + (x * 5.0))))) + (x * 2.5))))) else: tmp = wj + ((wj + (x / (-1.0 + (wj * (-1.0 - (wj * (0.5 + (wj * 0.16666666666666666)))))))) / (-1.0 - wj)) return tmp
function code(wj, x) tmp = 0.0 if (wj <= 1.3e-8) tmp = Float64(x + Float64(wj * Float64(Float64(x * -2.0) + Float64(wj * Float64(Float64(1.0 + Float64(wj * Float64(Float64(-1.0 - Float64(x * 0.6666666666666666)) - Float64(Float64(x * -3.0) + Float64(x * 5.0))))) + Float64(x * 2.5)))))); else tmp = Float64(wj + Float64(Float64(wj + Float64(x / Float64(-1.0 + Float64(wj * Float64(-1.0 - Float64(wj * Float64(0.5 + Float64(wj * 0.16666666666666666)))))))) / Float64(-1.0 - wj))); end return tmp end
function tmp_2 = code(wj, x) tmp = 0.0; if (wj <= 1.3e-8) tmp = x + (wj * ((x * -2.0) + (wj * ((1.0 + (wj * ((-1.0 - (x * 0.6666666666666666)) - ((x * -3.0) + (x * 5.0))))) + (x * 2.5))))); else tmp = wj + ((wj + (x / (-1.0 + (wj * (-1.0 - (wj * (0.5 + (wj * 0.16666666666666666)))))))) / (-1.0 - wj)); end tmp_2 = tmp; end
code[wj_, x_] := If[LessEqual[wj, 1.3e-8], N[(x + N[(wj * N[(N[(x * -2.0), $MachinePrecision] + N[(wj * N[(N[(1.0 + N[(wj * N[(N[(-1.0 - N[(x * 0.6666666666666666), $MachinePrecision]), $MachinePrecision] - N[(N[(x * -3.0), $MachinePrecision] + N[(x * 5.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(x * 2.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(wj + N[(N[(wj + N[(x / N[(-1.0 + N[(wj * N[(-1.0 - N[(wj * N[(0.5 + N[(wj * 0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(-1.0 - wj), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;wj \leq 1.3 \cdot 10^{-8}:\\
\;\;\;\;x + wj \cdot \left(x \cdot -2 + wj \cdot \left(\left(1 + wj \cdot \left(\left(-1 - x \cdot 0.6666666666666666\right) - \left(x \cdot -3 + x \cdot 5\right)\right)\right) + x \cdot 2.5\right)\right)\\
\mathbf{else}:\\
\;\;\;\;wj + \frac{wj + \frac{x}{-1 + wj \cdot \left(-1 - wj \cdot \left(0.5 + wj \cdot 0.16666666666666666\right)\right)}}{-1 - wj}\\
\end{array}
\end{array}
if wj < 1.3000000000000001e-8Initial program 78.4%
Taylor expanded in wj around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
cancel-sign-sub-invN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
Simplified97.3%
if 1.3000000000000001e-8 < wj Initial program 49.7%
sub-negN/A
+-lowering-+.f64N/A
distribute-rgt1-inN/A
associate-/l/N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
div-subN/A
associate-/l*N/A
*-inversesN/A
*-rgt-identityN/A
--lowering--.f64N/A
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
+-commutativeN/A
*-inversesN/A
distribute-neg-inN/A
Simplified99.7%
Taylor expanded in wj around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6498.3%
Simplified98.3%
Final simplification97.3%
(FPCore (wj x)
:precision binary64
(if (<= wj 1.3e-8)
(-
x
(*
x
(-
(* (* wj wj) (/ (+ wj -1.0) x))
(* wj (+ -2.0 (* wj (+ 2.5 (* wj -2.6666666666666665))))))))
(+
wj
(/
(+
wj
(/ x (+ -1.0 (* wj (- -1.0 (* wj (+ 0.5 (* wj 0.16666666666666666))))))))
(- -1.0 wj)))))
double code(double wj, double x) {
double tmp;
if (wj <= 1.3e-8) {
tmp = x - (x * (((wj * wj) * ((wj + -1.0) / x)) - (wj * (-2.0 + (wj * (2.5 + (wj * -2.6666666666666665)))))));
} else {
tmp = wj + ((wj + (x / (-1.0 + (wj * (-1.0 - (wj * (0.5 + (wj * 0.16666666666666666)))))))) / (-1.0 - wj));
}
return tmp;
}
real(8) function code(wj, x)
real(8), intent (in) :: wj
real(8), intent (in) :: x
real(8) :: tmp
if (wj <= 1.3d-8) then
tmp = x - (x * (((wj * wj) * ((wj + (-1.0d0)) / x)) - (wj * ((-2.0d0) + (wj * (2.5d0 + (wj * (-2.6666666666666665d0))))))))
else
tmp = wj + ((wj + (x / ((-1.0d0) + (wj * ((-1.0d0) - (wj * (0.5d0 + (wj * 0.16666666666666666d0)))))))) / ((-1.0d0) - wj))
end if
code = tmp
end function
public static double code(double wj, double x) {
double tmp;
if (wj <= 1.3e-8) {
tmp = x - (x * (((wj * wj) * ((wj + -1.0) / x)) - (wj * (-2.0 + (wj * (2.5 + (wj * -2.6666666666666665)))))));
} else {
tmp = wj + ((wj + (x / (-1.0 + (wj * (-1.0 - (wj * (0.5 + (wj * 0.16666666666666666)))))))) / (-1.0 - wj));
}
return tmp;
}
def code(wj, x): tmp = 0 if wj <= 1.3e-8: tmp = x - (x * (((wj * wj) * ((wj + -1.0) / x)) - (wj * (-2.0 + (wj * (2.5 + (wj * -2.6666666666666665))))))) else: tmp = wj + ((wj + (x / (-1.0 + (wj * (-1.0 - (wj * (0.5 + (wj * 0.16666666666666666)))))))) / (-1.0 - wj)) return tmp
function code(wj, x) tmp = 0.0 if (wj <= 1.3e-8) tmp = Float64(x - Float64(x * Float64(Float64(Float64(wj * wj) * Float64(Float64(wj + -1.0) / x)) - Float64(wj * Float64(-2.0 + Float64(wj * Float64(2.5 + Float64(wj * -2.6666666666666665)))))))); else tmp = Float64(wj + Float64(Float64(wj + Float64(x / Float64(-1.0 + Float64(wj * Float64(-1.0 - Float64(wj * Float64(0.5 + Float64(wj * 0.16666666666666666)))))))) / Float64(-1.0 - wj))); end return tmp end
function tmp_2 = code(wj, x) tmp = 0.0; if (wj <= 1.3e-8) tmp = x - (x * (((wj * wj) * ((wj + -1.0) / x)) - (wj * (-2.0 + (wj * (2.5 + (wj * -2.6666666666666665))))))); else tmp = wj + ((wj + (x / (-1.0 + (wj * (-1.0 - (wj * (0.5 + (wj * 0.16666666666666666)))))))) / (-1.0 - wj)); end tmp_2 = tmp; end
code[wj_, x_] := If[LessEqual[wj, 1.3e-8], N[(x - N[(x * N[(N[(N[(wj * wj), $MachinePrecision] * N[(N[(wj + -1.0), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision] - N[(wj * N[(-2.0 + N[(wj * N[(2.5 + N[(wj * -2.6666666666666665), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(wj + N[(N[(wj + N[(x / N[(-1.0 + N[(wj * N[(-1.0 - N[(wj * N[(0.5 + N[(wj * 0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(-1.0 - wj), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;wj \leq 1.3 \cdot 10^{-8}:\\
\;\;\;\;x - x \cdot \left(\left(wj \cdot wj\right) \cdot \frac{wj + -1}{x} - wj \cdot \left(-2 + wj \cdot \left(2.5 + wj \cdot -2.6666666666666665\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;wj + \frac{wj + \frac{x}{-1 + wj \cdot \left(-1 - wj \cdot \left(0.5 + wj \cdot 0.16666666666666666\right)\right)}}{-1 - wj}\\
\end{array}
\end{array}
if wj < 1.3000000000000001e-8Initial program 78.4%
Taylor expanded in wj around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
cancel-sign-sub-invN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
Simplified97.3%
Taylor expanded in x around inf
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
cancel-sign-sub-invN/A
metadata-evalN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-/l*N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
--lowering--.f6497.2%
Simplified97.2%
if 1.3000000000000001e-8 < wj Initial program 49.7%
sub-negN/A
+-lowering-+.f64N/A
distribute-rgt1-inN/A
associate-/l/N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
div-subN/A
associate-/l*N/A
*-inversesN/A
*-rgt-identityN/A
--lowering--.f64N/A
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
+-commutativeN/A
*-inversesN/A
distribute-neg-inN/A
Simplified99.7%
Taylor expanded in wj around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6498.3%
Simplified98.3%
Final simplification97.3%
(FPCore (wj x)
:precision binary64
(if (<= wj -0.00032)
(+ wj (/ (/ (- (* wj wj) (* x x)) (+ wj x)) (- -1.0 wj)))
(if (<= wj 1.3e-8)
(+ x (* wj (+ (* x -2.0) (* wj (- 1.0 wj)))))
(+ wj (/ (+ wj (/ x (+ -1.0 (* wj (- -1.0 (* wj 0.5)))))) (- -1.0 wj))))))
double code(double wj, double x) {
double tmp;
if (wj <= -0.00032) {
tmp = wj + ((((wj * wj) - (x * x)) / (wj + x)) / (-1.0 - wj));
} else if (wj <= 1.3e-8) {
tmp = x + (wj * ((x * -2.0) + (wj * (1.0 - wj))));
} else {
tmp = wj + ((wj + (x / (-1.0 + (wj * (-1.0 - (wj * 0.5)))))) / (-1.0 - wj));
}
return tmp;
}
real(8) function code(wj, x)
real(8), intent (in) :: wj
real(8), intent (in) :: x
real(8) :: tmp
if (wj <= (-0.00032d0)) then
tmp = wj + ((((wj * wj) - (x * x)) / (wj + x)) / ((-1.0d0) - wj))
else if (wj <= 1.3d-8) then
tmp = x + (wj * ((x * (-2.0d0)) + (wj * (1.0d0 - wj))))
else
tmp = wj + ((wj + (x / ((-1.0d0) + (wj * ((-1.0d0) - (wj * 0.5d0)))))) / ((-1.0d0) - wj))
end if
code = tmp
end function
public static double code(double wj, double x) {
double tmp;
if (wj <= -0.00032) {
tmp = wj + ((((wj * wj) - (x * x)) / (wj + x)) / (-1.0 - wj));
} else if (wj <= 1.3e-8) {
tmp = x + (wj * ((x * -2.0) + (wj * (1.0 - wj))));
} else {
tmp = wj + ((wj + (x / (-1.0 + (wj * (-1.0 - (wj * 0.5)))))) / (-1.0 - wj));
}
return tmp;
}
def code(wj, x): tmp = 0 if wj <= -0.00032: tmp = wj + ((((wj * wj) - (x * x)) / (wj + x)) / (-1.0 - wj)) elif wj <= 1.3e-8: tmp = x + (wj * ((x * -2.0) + (wj * (1.0 - wj)))) else: tmp = wj + ((wj + (x / (-1.0 + (wj * (-1.0 - (wj * 0.5)))))) / (-1.0 - wj)) return tmp
function code(wj, x) tmp = 0.0 if (wj <= -0.00032) tmp = Float64(wj + Float64(Float64(Float64(Float64(wj * wj) - Float64(x * x)) / Float64(wj + x)) / Float64(-1.0 - wj))); elseif (wj <= 1.3e-8) tmp = Float64(x + Float64(wj * Float64(Float64(x * -2.0) + Float64(wj * Float64(1.0 - wj))))); else tmp = Float64(wj + Float64(Float64(wj + Float64(x / Float64(-1.0 + Float64(wj * Float64(-1.0 - Float64(wj * 0.5)))))) / Float64(-1.0 - wj))); end return tmp end
function tmp_2 = code(wj, x) tmp = 0.0; if (wj <= -0.00032) tmp = wj + ((((wj * wj) - (x * x)) / (wj + x)) / (-1.0 - wj)); elseif (wj <= 1.3e-8) tmp = x + (wj * ((x * -2.0) + (wj * (1.0 - wj)))); else tmp = wj + ((wj + (x / (-1.0 + (wj * (-1.0 - (wj * 0.5)))))) / (-1.0 - wj)); end tmp_2 = tmp; end
code[wj_, x_] := If[LessEqual[wj, -0.00032], N[(wj + N[(N[(N[(N[(wj * wj), $MachinePrecision] - N[(x * x), $MachinePrecision]), $MachinePrecision] / N[(wj + x), $MachinePrecision]), $MachinePrecision] / N[(-1.0 - wj), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[wj, 1.3e-8], N[(x + N[(wj * N[(N[(x * -2.0), $MachinePrecision] + N[(wj * N[(1.0 - wj), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(wj + N[(N[(wj + N[(x / N[(-1.0 + N[(wj * N[(-1.0 - N[(wj * 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(-1.0 - wj), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;wj \leq -0.00032:\\
\;\;\;\;wj + \frac{\frac{wj \cdot wj - x \cdot x}{wj + x}}{-1 - wj}\\
\mathbf{elif}\;wj \leq 1.3 \cdot 10^{-8}:\\
\;\;\;\;x + wj \cdot \left(x \cdot -2 + wj \cdot \left(1 - wj\right)\right)\\
\mathbf{else}:\\
\;\;\;\;wj + \frac{wj + \frac{x}{-1 + wj \cdot \left(-1 - wj \cdot 0.5\right)}}{-1 - wj}\\
\end{array}
\end{array}
if wj < -3.20000000000000026e-4Initial program 51.4%
sub-negN/A
+-lowering-+.f64N/A
distribute-rgt1-inN/A
associate-/l/N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
div-subN/A
associate-/l*N/A
*-inversesN/A
*-rgt-identityN/A
--lowering--.f64N/A
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
+-commutativeN/A
*-inversesN/A
distribute-neg-inN/A
Simplified96.1%
Taylor expanded in wj around 0
Simplified45.6%
flip--N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f6466.7%
Applied egg-rr66.7%
if -3.20000000000000026e-4 < wj < 1.3000000000000001e-8Initial program 79.4%
Taylor expanded in wj around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
cancel-sign-sub-invN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
Simplified100.0%
Taylor expanded in x around 0
*-lowering-*.f64N/A
--lowering--.f6499.8%
Simplified99.8%
if 1.3000000000000001e-8 < wj Initial program 49.7%
sub-negN/A
+-lowering-+.f64N/A
distribute-rgt1-inN/A
associate-/l/N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
div-subN/A
associate-/l*N/A
*-inversesN/A
*-rgt-identityN/A
--lowering--.f64N/A
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
+-commutativeN/A
*-inversesN/A
distribute-neg-inN/A
Simplified99.7%
Taylor expanded in wj around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6495.1%
Simplified95.1%
Final simplification98.5%
(FPCore (wj x)
:precision binary64
(if (<= wj 1.3e-8)
(+
x
(*
x
(*
wj
(+
-2.0
(*
wj
(+ (/ 1.0 x) (+ 2.5 (* wj (+ -2.6666666666666665 (/ -1.0 x))))))))))
(+
wj
(/
(+
wj
(/ x (+ -1.0 (* wj (- -1.0 (* wj (+ 0.5 (* wj 0.16666666666666666))))))))
(- -1.0 wj)))))
double code(double wj, double x) {
double tmp;
if (wj <= 1.3e-8) {
tmp = x + (x * (wj * (-2.0 + (wj * ((1.0 / x) + (2.5 + (wj * (-2.6666666666666665 + (-1.0 / x)))))))));
} else {
tmp = wj + ((wj + (x / (-1.0 + (wj * (-1.0 - (wj * (0.5 + (wj * 0.16666666666666666)))))))) / (-1.0 - wj));
}
return tmp;
}
real(8) function code(wj, x)
real(8), intent (in) :: wj
real(8), intent (in) :: x
real(8) :: tmp
if (wj <= 1.3d-8) then
tmp = x + (x * (wj * ((-2.0d0) + (wj * ((1.0d0 / x) + (2.5d0 + (wj * ((-2.6666666666666665d0) + ((-1.0d0) / x)))))))))
else
tmp = wj + ((wj + (x / ((-1.0d0) + (wj * ((-1.0d0) - (wj * (0.5d0 + (wj * 0.16666666666666666d0)))))))) / ((-1.0d0) - wj))
end if
code = tmp
end function
public static double code(double wj, double x) {
double tmp;
if (wj <= 1.3e-8) {
tmp = x + (x * (wj * (-2.0 + (wj * ((1.0 / x) + (2.5 + (wj * (-2.6666666666666665 + (-1.0 / x)))))))));
} else {
tmp = wj + ((wj + (x / (-1.0 + (wj * (-1.0 - (wj * (0.5 + (wj * 0.16666666666666666)))))))) / (-1.0 - wj));
}
return tmp;
}
def code(wj, x): tmp = 0 if wj <= 1.3e-8: tmp = x + (x * (wj * (-2.0 + (wj * ((1.0 / x) + (2.5 + (wj * (-2.6666666666666665 + (-1.0 / x))))))))) else: tmp = wj + ((wj + (x / (-1.0 + (wj * (-1.0 - (wj * (0.5 + (wj * 0.16666666666666666)))))))) / (-1.0 - wj)) return tmp
function code(wj, x) tmp = 0.0 if (wj <= 1.3e-8) tmp = Float64(x + Float64(x * Float64(wj * Float64(-2.0 + Float64(wj * Float64(Float64(1.0 / x) + Float64(2.5 + Float64(wj * Float64(-2.6666666666666665 + Float64(-1.0 / x)))))))))); else tmp = Float64(wj + Float64(Float64(wj + Float64(x / Float64(-1.0 + Float64(wj * Float64(-1.0 - Float64(wj * Float64(0.5 + Float64(wj * 0.16666666666666666)))))))) / Float64(-1.0 - wj))); end return tmp end
function tmp_2 = code(wj, x) tmp = 0.0; if (wj <= 1.3e-8) tmp = x + (x * (wj * (-2.0 + (wj * ((1.0 / x) + (2.5 + (wj * (-2.6666666666666665 + (-1.0 / x))))))))); else tmp = wj + ((wj + (x / (-1.0 + (wj * (-1.0 - (wj * (0.5 + (wj * 0.16666666666666666)))))))) / (-1.0 - wj)); end tmp_2 = tmp; end
code[wj_, x_] := If[LessEqual[wj, 1.3e-8], N[(x + N[(x * N[(wj * N[(-2.0 + N[(wj * N[(N[(1.0 / x), $MachinePrecision] + N[(2.5 + N[(wj * N[(-2.6666666666666665 + N[(-1.0 / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(wj + N[(N[(wj + N[(x / N[(-1.0 + N[(wj * N[(-1.0 - N[(wj * N[(0.5 + N[(wj * 0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(-1.0 - wj), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;wj \leq 1.3 \cdot 10^{-8}:\\
\;\;\;\;x + x \cdot \left(wj \cdot \left(-2 + wj \cdot \left(\frac{1}{x} + \left(2.5 + wj \cdot \left(-2.6666666666666665 + \frac{-1}{x}\right)\right)\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;wj + \frac{wj + \frac{x}{-1 + wj \cdot \left(-1 - wj \cdot \left(0.5 + wj \cdot 0.16666666666666666\right)\right)}}{-1 - wj}\\
\end{array}
\end{array}
if wj < 1.3000000000000001e-8Initial program 78.4%
Taylor expanded in wj around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
cancel-sign-sub-invN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
Simplified97.3%
Taylor expanded in x around inf
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
cancel-sign-sub-invN/A
metadata-evalN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-/l*N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
--lowering--.f6497.2%
Simplified97.2%
Taylor expanded in wj around 0
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
mul-1-negN/A
distribute-rgt-neg-inN/A
*-lowering-*.f64N/A
distribute-neg-inN/A
metadata-evalN/A
+-lowering-+.f64N/A
distribute-neg-fracN/A
metadata-evalN/A
/-lowering-/.f6497.2%
Simplified97.2%
if 1.3000000000000001e-8 < wj Initial program 49.7%
sub-negN/A
+-lowering-+.f64N/A
distribute-rgt1-inN/A
associate-/l/N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
div-subN/A
associate-/l*N/A
*-inversesN/A
*-rgt-identityN/A
--lowering--.f64N/A
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
+-commutativeN/A
*-inversesN/A
distribute-neg-inN/A
Simplified99.7%
Taylor expanded in wj around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6498.3%
Simplified98.3%
Final simplification97.2%
(FPCore (wj x)
:precision binary64
(if (<= wj 1.3e-8)
(+ x (* x (* wj (+ -2.0 (* wj (+ (/ 1.0 x) (- 2.5 (/ wj x))))))))
(+
wj
(/
(+
wj
(/ x (+ -1.0 (* wj (- -1.0 (* wj (+ 0.5 (* wj 0.16666666666666666))))))))
(- -1.0 wj)))))
double code(double wj, double x) {
double tmp;
if (wj <= 1.3e-8) {
tmp = x + (x * (wj * (-2.0 + (wj * ((1.0 / x) + (2.5 - (wj / x)))))));
} else {
tmp = wj + ((wj + (x / (-1.0 + (wj * (-1.0 - (wj * (0.5 + (wj * 0.16666666666666666)))))))) / (-1.0 - wj));
}
return tmp;
}
real(8) function code(wj, x)
real(8), intent (in) :: wj
real(8), intent (in) :: x
real(8) :: tmp
if (wj <= 1.3d-8) then
tmp = x + (x * (wj * ((-2.0d0) + (wj * ((1.0d0 / x) + (2.5d0 - (wj / x)))))))
else
tmp = wj + ((wj + (x / ((-1.0d0) + (wj * ((-1.0d0) - (wj * (0.5d0 + (wj * 0.16666666666666666d0)))))))) / ((-1.0d0) - wj))
end if
code = tmp
end function
public static double code(double wj, double x) {
double tmp;
if (wj <= 1.3e-8) {
tmp = x + (x * (wj * (-2.0 + (wj * ((1.0 / x) + (2.5 - (wj / x)))))));
} else {
tmp = wj + ((wj + (x / (-1.0 + (wj * (-1.0 - (wj * (0.5 + (wj * 0.16666666666666666)))))))) / (-1.0 - wj));
}
return tmp;
}
def code(wj, x): tmp = 0 if wj <= 1.3e-8: tmp = x + (x * (wj * (-2.0 + (wj * ((1.0 / x) + (2.5 - (wj / x))))))) else: tmp = wj + ((wj + (x / (-1.0 + (wj * (-1.0 - (wj * (0.5 + (wj * 0.16666666666666666)))))))) / (-1.0 - wj)) return tmp
function code(wj, x) tmp = 0.0 if (wj <= 1.3e-8) tmp = Float64(x + Float64(x * Float64(wj * Float64(-2.0 + Float64(wj * Float64(Float64(1.0 / x) + Float64(2.5 - Float64(wj / x)))))))); else tmp = Float64(wj + Float64(Float64(wj + Float64(x / Float64(-1.0 + Float64(wj * Float64(-1.0 - Float64(wj * Float64(0.5 + Float64(wj * 0.16666666666666666)))))))) / Float64(-1.0 - wj))); end return tmp end
function tmp_2 = code(wj, x) tmp = 0.0; if (wj <= 1.3e-8) tmp = x + (x * (wj * (-2.0 + (wj * ((1.0 / x) + (2.5 - (wj / x))))))); else tmp = wj + ((wj + (x / (-1.0 + (wj * (-1.0 - (wj * (0.5 + (wj * 0.16666666666666666)))))))) / (-1.0 - wj)); end tmp_2 = tmp; end
code[wj_, x_] := If[LessEqual[wj, 1.3e-8], N[(x + N[(x * N[(wj * N[(-2.0 + N[(wj * N[(N[(1.0 / x), $MachinePrecision] + N[(2.5 - N[(wj / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(wj + N[(N[(wj + N[(x / N[(-1.0 + N[(wj * N[(-1.0 - N[(wj * N[(0.5 + N[(wj * 0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(-1.0 - wj), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;wj \leq 1.3 \cdot 10^{-8}:\\
\;\;\;\;x + x \cdot \left(wj \cdot \left(-2 + wj \cdot \left(\frac{1}{x} + \left(2.5 - \frac{wj}{x}\right)\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;wj + \frac{wj + \frac{x}{-1 + wj \cdot \left(-1 - wj \cdot \left(0.5 + wj \cdot 0.16666666666666666\right)\right)}}{-1 - wj}\\
\end{array}
\end{array}
if wj < 1.3000000000000001e-8Initial program 78.4%
Taylor expanded in wj around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
cancel-sign-sub-invN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
Simplified97.3%
Taylor expanded in x around inf
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
cancel-sign-sub-invN/A
metadata-evalN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-/l*N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
--lowering--.f6497.2%
Simplified97.2%
Taylor expanded in wj around 0
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
mul-1-negN/A
distribute-rgt-neg-inN/A
*-lowering-*.f64N/A
distribute-neg-inN/A
metadata-evalN/A
+-lowering-+.f64N/A
distribute-neg-fracN/A
metadata-evalN/A
/-lowering-/.f6497.2%
Simplified97.2%
Taylor expanded in x around 0
mul-1-negN/A
distribute-neg-frac2N/A
mul-1-negN/A
/-lowering-/.f64N/A
mul-1-negN/A
neg-sub0N/A
--lowering--.f6497.1%
Simplified97.1%
if 1.3000000000000001e-8 < wj Initial program 49.7%
sub-negN/A
+-lowering-+.f64N/A
distribute-rgt1-inN/A
associate-/l/N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
div-subN/A
associate-/l*N/A
*-inversesN/A
*-rgt-identityN/A
--lowering--.f64N/A
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
+-commutativeN/A
*-inversesN/A
distribute-neg-inN/A
Simplified99.7%
Taylor expanded in wj around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6498.3%
Simplified98.3%
Final simplification97.2%
(FPCore (wj x) :precision binary64 (if (<= wj 1.3e-8) (+ x (* x (* wj (+ -2.0 (* wj (+ (/ 1.0 x) (- 2.5 (/ wj x)))))))) (+ wj (/ (+ wj (/ x (+ -1.0 (* wj (- -1.0 (* wj 0.5)))))) (- -1.0 wj)))))
double code(double wj, double x) {
double tmp;
if (wj <= 1.3e-8) {
tmp = x + (x * (wj * (-2.0 + (wj * ((1.0 / x) + (2.5 - (wj / x)))))));
} else {
tmp = wj + ((wj + (x / (-1.0 + (wj * (-1.0 - (wj * 0.5)))))) / (-1.0 - wj));
}
return tmp;
}
real(8) function code(wj, x)
real(8), intent (in) :: wj
real(8), intent (in) :: x
real(8) :: tmp
if (wj <= 1.3d-8) then
tmp = x + (x * (wj * ((-2.0d0) + (wj * ((1.0d0 / x) + (2.5d0 - (wj / x)))))))
else
tmp = wj + ((wj + (x / ((-1.0d0) + (wj * ((-1.0d0) - (wj * 0.5d0)))))) / ((-1.0d0) - wj))
end if
code = tmp
end function
public static double code(double wj, double x) {
double tmp;
if (wj <= 1.3e-8) {
tmp = x + (x * (wj * (-2.0 + (wj * ((1.0 / x) + (2.5 - (wj / x)))))));
} else {
tmp = wj + ((wj + (x / (-1.0 + (wj * (-1.0 - (wj * 0.5)))))) / (-1.0 - wj));
}
return tmp;
}
def code(wj, x): tmp = 0 if wj <= 1.3e-8: tmp = x + (x * (wj * (-2.0 + (wj * ((1.0 / x) + (2.5 - (wj / x))))))) else: tmp = wj + ((wj + (x / (-1.0 + (wj * (-1.0 - (wj * 0.5)))))) / (-1.0 - wj)) return tmp
function code(wj, x) tmp = 0.0 if (wj <= 1.3e-8) tmp = Float64(x + Float64(x * Float64(wj * Float64(-2.0 + Float64(wj * Float64(Float64(1.0 / x) + Float64(2.5 - Float64(wj / x)))))))); else tmp = Float64(wj + Float64(Float64(wj + Float64(x / Float64(-1.0 + Float64(wj * Float64(-1.0 - Float64(wj * 0.5)))))) / Float64(-1.0 - wj))); end return tmp end
function tmp_2 = code(wj, x) tmp = 0.0; if (wj <= 1.3e-8) tmp = x + (x * (wj * (-2.0 + (wj * ((1.0 / x) + (2.5 - (wj / x))))))); else tmp = wj + ((wj + (x / (-1.0 + (wj * (-1.0 - (wj * 0.5)))))) / (-1.0 - wj)); end tmp_2 = tmp; end
code[wj_, x_] := If[LessEqual[wj, 1.3e-8], N[(x + N[(x * N[(wj * N[(-2.0 + N[(wj * N[(N[(1.0 / x), $MachinePrecision] + N[(2.5 - N[(wj / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(wj + N[(N[(wj + N[(x / N[(-1.0 + N[(wj * N[(-1.0 - N[(wj * 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(-1.0 - wj), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;wj \leq 1.3 \cdot 10^{-8}:\\
\;\;\;\;x + x \cdot \left(wj \cdot \left(-2 + wj \cdot \left(\frac{1}{x} + \left(2.5 - \frac{wj}{x}\right)\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;wj + \frac{wj + \frac{x}{-1 + wj \cdot \left(-1 - wj \cdot 0.5\right)}}{-1 - wj}\\
\end{array}
\end{array}
if wj < 1.3000000000000001e-8Initial program 78.4%
Taylor expanded in wj around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
cancel-sign-sub-invN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
Simplified97.3%
Taylor expanded in x around inf
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
cancel-sign-sub-invN/A
metadata-evalN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-/l*N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
--lowering--.f6497.2%
Simplified97.2%
Taylor expanded in wj around 0
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
mul-1-negN/A
distribute-rgt-neg-inN/A
*-lowering-*.f64N/A
distribute-neg-inN/A
metadata-evalN/A
+-lowering-+.f64N/A
distribute-neg-fracN/A
metadata-evalN/A
/-lowering-/.f6497.2%
Simplified97.2%
Taylor expanded in x around 0
mul-1-negN/A
distribute-neg-frac2N/A
mul-1-negN/A
/-lowering-/.f64N/A
mul-1-negN/A
neg-sub0N/A
--lowering--.f6497.1%
Simplified97.1%
if 1.3000000000000001e-8 < wj Initial program 49.7%
sub-negN/A
+-lowering-+.f64N/A
distribute-rgt1-inN/A
associate-/l/N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
div-subN/A
associate-/l*N/A
*-inversesN/A
*-rgt-identityN/A
--lowering--.f64N/A
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
+-commutativeN/A
*-inversesN/A
distribute-neg-inN/A
Simplified99.7%
Taylor expanded in wj around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6495.1%
Simplified95.1%
Final simplification97.1%
(FPCore (wj x)
:precision binary64
(if (<= wj -0.000235)
(+ wj (/ (/ (- (* wj wj) (* x x)) (+ wj x)) (- -1.0 wj)))
(if (<= wj 1.3e-8)
(+ x (* wj (+ (* x -2.0) (* wj (- 1.0 wj)))))
(+ wj (/ (+ wj (/ x (- -1.0 wj))) (- -1.0 wj))))))
double code(double wj, double x) {
double tmp;
if (wj <= -0.000235) {
tmp = wj + ((((wj * wj) - (x * x)) / (wj + x)) / (-1.0 - wj));
} else if (wj <= 1.3e-8) {
tmp = x + (wj * ((x * -2.0) + (wj * (1.0 - wj))));
} else {
tmp = wj + ((wj + (x / (-1.0 - wj))) / (-1.0 - wj));
}
return tmp;
}
real(8) function code(wj, x)
real(8), intent (in) :: wj
real(8), intent (in) :: x
real(8) :: tmp
if (wj <= (-0.000235d0)) then
tmp = wj + ((((wj * wj) - (x * x)) / (wj + x)) / ((-1.0d0) - wj))
else if (wj <= 1.3d-8) then
tmp = x + (wj * ((x * (-2.0d0)) + (wj * (1.0d0 - wj))))
else
tmp = wj + ((wj + (x / ((-1.0d0) - wj))) / ((-1.0d0) - wj))
end if
code = tmp
end function
public static double code(double wj, double x) {
double tmp;
if (wj <= -0.000235) {
tmp = wj + ((((wj * wj) - (x * x)) / (wj + x)) / (-1.0 - wj));
} else if (wj <= 1.3e-8) {
tmp = x + (wj * ((x * -2.0) + (wj * (1.0 - wj))));
} else {
tmp = wj + ((wj + (x / (-1.0 - wj))) / (-1.0 - wj));
}
return tmp;
}
def code(wj, x): tmp = 0 if wj <= -0.000235: tmp = wj + ((((wj * wj) - (x * x)) / (wj + x)) / (-1.0 - wj)) elif wj <= 1.3e-8: tmp = x + (wj * ((x * -2.0) + (wj * (1.0 - wj)))) else: tmp = wj + ((wj + (x / (-1.0 - wj))) / (-1.0 - wj)) return tmp
function code(wj, x) tmp = 0.0 if (wj <= -0.000235) tmp = Float64(wj + Float64(Float64(Float64(Float64(wj * wj) - Float64(x * x)) / Float64(wj + x)) / Float64(-1.0 - wj))); elseif (wj <= 1.3e-8) tmp = Float64(x + Float64(wj * Float64(Float64(x * -2.0) + Float64(wj * Float64(1.0 - wj))))); else tmp = Float64(wj + Float64(Float64(wj + Float64(x / Float64(-1.0 - wj))) / Float64(-1.0 - wj))); end return tmp end
function tmp_2 = code(wj, x) tmp = 0.0; if (wj <= -0.000235) tmp = wj + ((((wj * wj) - (x * x)) / (wj + x)) / (-1.0 - wj)); elseif (wj <= 1.3e-8) tmp = x + (wj * ((x * -2.0) + (wj * (1.0 - wj)))); else tmp = wj + ((wj + (x / (-1.0 - wj))) / (-1.0 - wj)); end tmp_2 = tmp; end
code[wj_, x_] := If[LessEqual[wj, -0.000235], N[(wj + N[(N[(N[(N[(wj * wj), $MachinePrecision] - N[(x * x), $MachinePrecision]), $MachinePrecision] / N[(wj + x), $MachinePrecision]), $MachinePrecision] / N[(-1.0 - wj), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[wj, 1.3e-8], N[(x + N[(wj * N[(N[(x * -2.0), $MachinePrecision] + N[(wj * N[(1.0 - wj), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(wj + N[(N[(wj + N[(x / N[(-1.0 - wj), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(-1.0 - wj), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;wj \leq -0.000235:\\
\;\;\;\;wj + \frac{\frac{wj \cdot wj - x \cdot x}{wj + x}}{-1 - wj}\\
\mathbf{elif}\;wj \leq 1.3 \cdot 10^{-8}:\\
\;\;\;\;x + wj \cdot \left(x \cdot -2 + wj \cdot \left(1 - wj\right)\right)\\
\mathbf{else}:\\
\;\;\;\;wj + \frac{wj + \frac{x}{-1 - wj}}{-1 - wj}\\
\end{array}
\end{array}
if wj < -2.34999999999999993e-4Initial program 51.4%
sub-negN/A
+-lowering-+.f64N/A
distribute-rgt1-inN/A
associate-/l/N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
div-subN/A
associate-/l*N/A
*-inversesN/A
*-rgt-identityN/A
--lowering--.f64N/A
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
+-commutativeN/A
*-inversesN/A
distribute-neg-inN/A
Simplified96.1%
Taylor expanded in wj around 0
Simplified45.6%
flip--N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f6466.7%
Applied egg-rr66.7%
if -2.34999999999999993e-4 < wj < 1.3000000000000001e-8Initial program 79.4%
Taylor expanded in wj around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
cancel-sign-sub-invN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
Simplified100.0%
Taylor expanded in x around 0
*-lowering-*.f64N/A
--lowering--.f6499.8%
Simplified99.8%
if 1.3000000000000001e-8 < wj Initial program 49.7%
sub-negN/A
+-lowering-+.f64N/A
distribute-rgt1-inN/A
associate-/l/N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
div-subN/A
associate-/l*N/A
*-inversesN/A
*-rgt-identityN/A
--lowering--.f64N/A
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
+-commutativeN/A
*-inversesN/A
distribute-neg-inN/A
Simplified99.7%
Taylor expanded in wj around 0
+-commutativeN/A
+-lowering-+.f6492.0%
Simplified92.0%
Final simplification98.4%
(FPCore (wj x) :precision binary64 (let* ((t_0 (+ wj (/ (- wj x) (- -1.0 wj))))) (if (<= wj -1.15e-9) t_0 (if (<= wj 2.2e-14) (+ x (* wj (* x -2.0))) t_0))))
double code(double wj, double x) {
double t_0 = wj + ((wj - x) / (-1.0 - wj));
double tmp;
if (wj <= -1.15e-9) {
tmp = t_0;
} else if (wj <= 2.2e-14) {
tmp = x + (wj * (x * -2.0));
} else {
tmp = t_0;
}
return tmp;
}
real(8) function code(wj, x)
real(8), intent (in) :: wj
real(8), intent (in) :: x
real(8) :: t_0
real(8) :: tmp
t_0 = wj + ((wj - x) / ((-1.0d0) - wj))
if (wj <= (-1.15d-9)) then
tmp = t_0
else if (wj <= 2.2d-14) then
tmp = x + (wj * (x * (-2.0d0)))
else
tmp = t_0
end if
code = tmp
end function
public static double code(double wj, double x) {
double t_0 = wj + ((wj - x) / (-1.0 - wj));
double tmp;
if (wj <= -1.15e-9) {
tmp = t_0;
} else if (wj <= 2.2e-14) {
tmp = x + (wj * (x * -2.0));
} else {
tmp = t_0;
}
return tmp;
}
def code(wj, x): t_0 = wj + ((wj - x) / (-1.0 - wj)) tmp = 0 if wj <= -1.15e-9: tmp = t_0 elif wj <= 2.2e-14: tmp = x + (wj * (x * -2.0)) else: tmp = t_0 return tmp
function code(wj, x) t_0 = Float64(wj + Float64(Float64(wj - x) / Float64(-1.0 - wj))) tmp = 0.0 if (wj <= -1.15e-9) tmp = t_0; elseif (wj <= 2.2e-14) tmp = Float64(x + Float64(wj * Float64(x * -2.0))); else tmp = t_0; end return tmp end
function tmp_2 = code(wj, x) t_0 = wj + ((wj - x) / (-1.0 - wj)); tmp = 0.0; if (wj <= -1.15e-9) tmp = t_0; elseif (wj <= 2.2e-14) tmp = x + (wj * (x * -2.0)); else tmp = t_0; end tmp_2 = tmp; end
code[wj_, x_] := Block[{t$95$0 = N[(wj + N[(N[(wj - x), $MachinePrecision] / N[(-1.0 - wj), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[wj, -1.15e-9], t$95$0, If[LessEqual[wj, 2.2e-14], N[(x + N[(wj * N[(x * -2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$0]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := wj + \frac{wj - x}{-1 - wj}\\
\mathbf{if}\;wj \leq -1.15 \cdot 10^{-9}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;wj \leq 2.2 \cdot 10^{-14}:\\
\;\;\;\;x + wj \cdot \left(x \cdot -2\right)\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if wj < -1.15e-9 or 2.2000000000000001e-14 < wj Initial program 57.3%
sub-negN/A
+-lowering-+.f64N/A
distribute-rgt1-inN/A
associate-/l/N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
div-subN/A
associate-/l*N/A
*-inversesN/A
*-rgt-identityN/A
--lowering--.f64N/A
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
+-commutativeN/A
*-inversesN/A
distribute-neg-inN/A
Simplified90.5%
Taylor expanded in wj around 0
Simplified59.1%
if -1.15e-9 < wj < 2.2000000000000001e-14Initial program 79.6%
Taylor expanded in wj around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
cancel-sign-sub-invN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
Simplified100.0%
Taylor expanded in wj around 0
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6491.7%
Simplified91.7%
(FPCore (wj x) :precision binary64 (if (<= wj 1.3e-8) (+ x (* wj (+ (* x -2.0) (* wj (- 1.0 wj))))) (+ wj (/ (+ wj (/ x (- -1.0 wj))) (- -1.0 wj)))))
double code(double wj, double x) {
double tmp;
if (wj <= 1.3e-8) {
tmp = x + (wj * ((x * -2.0) + (wj * (1.0 - wj))));
} else {
tmp = wj + ((wj + (x / (-1.0 - wj))) / (-1.0 - wj));
}
return tmp;
}
real(8) function code(wj, x)
real(8), intent (in) :: wj
real(8), intent (in) :: x
real(8) :: tmp
if (wj <= 1.3d-8) then
tmp = x + (wj * ((x * (-2.0d0)) + (wj * (1.0d0 - wj))))
else
tmp = wj + ((wj + (x / ((-1.0d0) - wj))) / ((-1.0d0) - wj))
end if
code = tmp
end function
public static double code(double wj, double x) {
double tmp;
if (wj <= 1.3e-8) {
tmp = x + (wj * ((x * -2.0) + (wj * (1.0 - wj))));
} else {
tmp = wj + ((wj + (x / (-1.0 - wj))) / (-1.0 - wj));
}
return tmp;
}
def code(wj, x): tmp = 0 if wj <= 1.3e-8: tmp = x + (wj * ((x * -2.0) + (wj * (1.0 - wj)))) else: tmp = wj + ((wj + (x / (-1.0 - wj))) / (-1.0 - wj)) return tmp
function code(wj, x) tmp = 0.0 if (wj <= 1.3e-8) tmp = Float64(x + Float64(wj * Float64(Float64(x * -2.0) + Float64(wj * Float64(1.0 - wj))))); else tmp = Float64(wj + Float64(Float64(wj + Float64(x / Float64(-1.0 - wj))) / Float64(-1.0 - wj))); end return tmp end
function tmp_2 = code(wj, x) tmp = 0.0; if (wj <= 1.3e-8) tmp = x + (wj * ((x * -2.0) + (wj * (1.0 - wj)))); else tmp = wj + ((wj + (x / (-1.0 - wj))) / (-1.0 - wj)); end tmp_2 = tmp; end
code[wj_, x_] := If[LessEqual[wj, 1.3e-8], N[(x + N[(wj * N[(N[(x * -2.0), $MachinePrecision] + N[(wj * N[(1.0 - wj), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(wj + N[(N[(wj + N[(x / N[(-1.0 - wj), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(-1.0 - wj), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;wj \leq 1.3 \cdot 10^{-8}:\\
\;\;\;\;x + wj \cdot \left(x \cdot -2 + wj \cdot \left(1 - wj\right)\right)\\
\mathbf{else}:\\
\;\;\;\;wj + \frac{wj + \frac{x}{-1 - wj}}{-1 - wj}\\
\end{array}
\end{array}
if wj < 1.3000000000000001e-8Initial program 78.4%
Taylor expanded in wj around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
cancel-sign-sub-invN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
Simplified97.3%
Taylor expanded in x around 0
*-lowering-*.f64N/A
--lowering--.f6497.0%
Simplified97.0%
if 1.3000000000000001e-8 < wj Initial program 49.7%
sub-negN/A
+-lowering-+.f64N/A
distribute-rgt1-inN/A
associate-/l/N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
div-subN/A
associate-/l*N/A
*-inversesN/A
*-rgt-identityN/A
--lowering--.f64N/A
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
+-commutativeN/A
*-inversesN/A
distribute-neg-inN/A
Simplified99.7%
Taylor expanded in wj around 0
+-commutativeN/A
+-lowering-+.f6492.0%
Simplified92.0%
Final simplification96.9%
(FPCore (wj x) :precision binary64 (if (<= wj 5.7e-9) (+ x (* wj (* wj (- 1.0 wj)))) (+ wj (/ (+ wj (/ x (- -1.0 wj))) (- -1.0 wj)))))
double code(double wj, double x) {
double tmp;
if (wj <= 5.7e-9) {
tmp = x + (wj * (wj * (1.0 - wj)));
} else {
tmp = wj + ((wj + (x / (-1.0 - wj))) / (-1.0 - wj));
}
return tmp;
}
real(8) function code(wj, x)
real(8), intent (in) :: wj
real(8), intent (in) :: x
real(8) :: tmp
if (wj <= 5.7d-9) then
tmp = x + (wj * (wj * (1.0d0 - wj)))
else
tmp = wj + ((wj + (x / ((-1.0d0) - wj))) / ((-1.0d0) - wj))
end if
code = tmp
end function
public static double code(double wj, double x) {
double tmp;
if (wj <= 5.7e-9) {
tmp = x + (wj * (wj * (1.0 - wj)));
} else {
tmp = wj + ((wj + (x / (-1.0 - wj))) / (-1.0 - wj));
}
return tmp;
}
def code(wj, x): tmp = 0 if wj <= 5.7e-9: tmp = x + (wj * (wj * (1.0 - wj))) else: tmp = wj + ((wj + (x / (-1.0 - wj))) / (-1.0 - wj)) return tmp
function code(wj, x) tmp = 0.0 if (wj <= 5.7e-9) tmp = Float64(x + Float64(wj * Float64(wj * Float64(1.0 - wj)))); else tmp = Float64(wj + Float64(Float64(wj + Float64(x / Float64(-1.0 - wj))) / Float64(-1.0 - wj))); end return tmp end
function tmp_2 = code(wj, x) tmp = 0.0; if (wj <= 5.7e-9) tmp = x + (wj * (wj * (1.0 - wj))); else tmp = wj + ((wj + (x / (-1.0 - wj))) / (-1.0 - wj)); end tmp_2 = tmp; end
code[wj_, x_] := If[LessEqual[wj, 5.7e-9], N[(x + N[(wj * N[(wj * N[(1.0 - wj), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(wj + N[(N[(wj + N[(x / N[(-1.0 - wj), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(-1.0 - wj), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;wj \leq 5.7 \cdot 10^{-9}:\\
\;\;\;\;x + wj \cdot \left(wj \cdot \left(1 - wj\right)\right)\\
\mathbf{else}:\\
\;\;\;\;wj + \frac{wj + \frac{x}{-1 - wj}}{-1 - wj}\\
\end{array}
\end{array}
if wj < 5.6999999999999998e-9Initial program 78.4%
Taylor expanded in wj around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
cancel-sign-sub-invN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
Simplified97.3%
Taylor expanded in x around 0
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
--lowering--.f6496.3%
Simplified96.3%
if 5.6999999999999998e-9 < wj Initial program 49.7%
sub-negN/A
+-lowering-+.f64N/A
distribute-rgt1-inN/A
associate-/l/N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
div-subN/A
associate-/l*N/A
*-inversesN/A
*-rgt-identityN/A
--lowering--.f64N/A
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
+-commutativeN/A
*-inversesN/A
distribute-neg-inN/A
Simplified99.7%
Taylor expanded in wj around 0
+-commutativeN/A
+-lowering-+.f6492.0%
Simplified92.0%
Final simplification96.2%
(FPCore (wj x) :precision binary64 (+ x (* wj (* wj (- 1.0 wj)))))
double code(double wj, double x) {
return x + (wj * (wj * (1.0 - wj)));
}
real(8) function code(wj, x)
real(8), intent (in) :: wj
real(8), intent (in) :: x
code = x + (wj * (wj * (1.0d0 - wj)))
end function
public static double code(double wj, double x) {
return x + (wj * (wj * (1.0 - wj)));
}
def code(wj, x): return x + (wj * (wj * (1.0 - wj)))
function code(wj, x) return Float64(x + Float64(wj * Float64(wj * Float64(1.0 - wj)))) end
function tmp = code(wj, x) tmp = x + (wj * (wj * (1.0 - wj))); end
code[wj_, x_] := N[(x + N[(wj * N[(wj * N[(1.0 - wj), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x + wj \cdot \left(wj \cdot \left(1 - wj\right)\right)
\end{array}
Initial program 77.7%
Taylor expanded in wj around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
cancel-sign-sub-invN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
Simplified95.7%
Taylor expanded in x around 0
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
--lowering--.f6494.4%
Simplified94.4%
(FPCore (wj x) :precision binary64 (+ x (* wj (* x -2.0))))
double code(double wj, double x) {
return x + (wj * (x * -2.0));
}
real(8) function code(wj, x)
real(8), intent (in) :: wj
real(8), intent (in) :: x
code = x + (wj * (x * (-2.0d0)))
end function
public static double code(double wj, double x) {
return x + (wj * (x * -2.0));
}
def code(wj, x): return x + (wj * (x * -2.0))
function code(wj, x) return Float64(x + Float64(wj * Float64(x * -2.0))) end
function tmp = code(wj, x) tmp = x + (wj * (x * -2.0)); end
code[wj_, x_] := N[(x + N[(wj * N[(x * -2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x + wj \cdot \left(x \cdot -2\right)
\end{array}
Initial program 77.7%
Taylor expanded in wj around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
cancel-sign-sub-invN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
Simplified95.7%
Taylor expanded in wj around 0
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6485.8%
Simplified85.8%
(FPCore (wj x) :precision binary64 (* x (+ 1.0 (* wj -2.0))))
double code(double wj, double x) {
return x * (1.0 + (wj * -2.0));
}
real(8) function code(wj, x)
real(8), intent (in) :: wj
real(8), intent (in) :: x
code = x * (1.0d0 + (wj * (-2.0d0)))
end function
public static double code(double wj, double x) {
return x * (1.0 + (wj * -2.0));
}
def code(wj, x): return x * (1.0 + (wj * -2.0))
function code(wj, x) return Float64(x * Float64(1.0 + Float64(wj * -2.0))) end
function tmp = code(wj, x) tmp = x * (1.0 + (wj * -2.0)); end
code[wj_, x_] := N[(x * N[(1.0 + N[(wj * -2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x \cdot \left(1 + wj \cdot -2\right)
\end{array}
Initial program 77.7%
Taylor expanded in wj around 0
associate-*r*N/A
distribute-rgt1-inN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6485.8%
Simplified85.8%
Final simplification85.8%
(FPCore (wj x) :precision binary64 (/ x (+ wj 1.0)))
double code(double wj, double x) {
return x / (wj + 1.0);
}
real(8) function code(wj, x)
real(8), intent (in) :: wj
real(8), intent (in) :: x
code = x / (wj + 1.0d0)
end function
public static double code(double wj, double x) {
return x / (wj + 1.0);
}
def code(wj, x): return x / (wj + 1.0)
function code(wj, x) return Float64(x / Float64(wj + 1.0)) end
function tmp = code(wj, x) tmp = x / (wj + 1.0); end
code[wj_, x_] := N[(x / N[(wj + 1.0), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{x}{wj + 1}
\end{array}
Initial program 77.7%
sub-negN/A
+-lowering-+.f64N/A
distribute-rgt1-inN/A
associate-/l/N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
div-subN/A
associate-/l*N/A
*-inversesN/A
*-rgt-identityN/A
--lowering--.f64N/A
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
+-commutativeN/A
*-inversesN/A
distribute-neg-inN/A
Simplified80.5%
Taylor expanded in wj around 0
Simplified77.7%
Taylor expanded in x around inf
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f6485.1%
Simplified85.1%
(FPCore (wj x) :precision binary64 x)
double code(double wj, double x) {
return x;
}
real(8) function code(wj, x)
real(8), intent (in) :: wj
real(8), intent (in) :: x
code = x
end function
public static double code(double wj, double x) {
return x;
}
def code(wj, x): return x
function code(wj, x) return x end
function tmp = code(wj, x) tmp = x; end
code[wj_, x_] := x
\begin{array}{l}
\\
x
\end{array}
Initial program 77.7%
sub-negN/A
+-lowering-+.f64N/A
distribute-rgt1-inN/A
associate-/l/N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
div-subN/A
associate-/l*N/A
*-inversesN/A
*-rgt-identityN/A
--lowering--.f64N/A
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
+-commutativeN/A
*-inversesN/A
distribute-neg-inN/A
Simplified80.5%
Taylor expanded in wj around 0
Simplified85.0%
(FPCore (wj x) :precision binary64 wj)
double code(double wj, double x) {
return wj;
}
real(8) function code(wj, x)
real(8), intent (in) :: wj
real(8), intent (in) :: x
code = wj
end function
public static double code(double wj, double x) {
return wj;
}
def code(wj, x): return wj
function code(wj, x) return wj end
function tmp = code(wj, x) tmp = wj; end
code[wj_, x_] := wj
\begin{array}{l}
\\
wj
\end{array}
Initial program 77.7%
sub-negN/A
+-lowering-+.f64N/A
distribute-rgt1-inN/A
associate-/l/N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
div-subN/A
associate-/l*N/A
*-inversesN/A
*-rgt-identityN/A
--lowering--.f64N/A
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
+-commutativeN/A
*-inversesN/A
distribute-neg-inN/A
Simplified80.5%
Taylor expanded in wj around inf
Simplified4.3%
(FPCore (wj x) :precision binary64 -1.0)
double code(double wj, double x) {
return -1.0;
}
real(8) function code(wj, x)
real(8), intent (in) :: wj
real(8), intent (in) :: x
code = -1.0d0
end function
public static double code(double wj, double x) {
return -1.0;
}
def code(wj, x): return -1.0
function code(wj, x) return -1.0 end
function tmp = code(wj, x) tmp = -1.0; end
code[wj_, x_] := -1.0
\begin{array}{l}
\\
-1
\end{array}
Initial program 77.7%
Taylor expanded in wj around inf
Simplified4.5%
Taylor expanded in wj around 0
Simplified3.5%
(FPCore (wj x) :precision binary64 (- wj (- (/ wj (+ wj 1.0)) (/ x (+ (exp wj) (* wj (exp wj)))))))
double code(double wj, double x) {
return wj - ((wj / (wj + 1.0)) - (x / (exp(wj) + (wj * exp(wj)))));
}
real(8) function code(wj, x)
real(8), intent (in) :: wj
real(8), intent (in) :: x
code = wj - ((wj / (wj + 1.0d0)) - (x / (exp(wj) + (wj * exp(wj)))))
end function
public static double code(double wj, double x) {
return wj - ((wj / (wj + 1.0)) - (x / (Math.exp(wj) + (wj * Math.exp(wj)))));
}
def code(wj, x): return wj - ((wj / (wj + 1.0)) - (x / (math.exp(wj) + (wj * math.exp(wj)))))
function code(wj, x) return Float64(wj - Float64(Float64(wj / Float64(wj + 1.0)) - Float64(x / Float64(exp(wj) + Float64(wj * exp(wj)))))) end
function tmp = code(wj, x) tmp = wj - ((wj / (wj + 1.0)) - (x / (exp(wj) + (wj * exp(wj))))); end
code[wj_, x_] := N[(wj - N[(N[(wj / N[(wj + 1.0), $MachinePrecision]), $MachinePrecision] - N[(x / N[(N[Exp[wj], $MachinePrecision] + N[(wj * N[Exp[wj], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
wj - \left(\frac{wj}{wj + 1} - \frac{x}{e^{wj} + wj \cdot e^{wj}}\right)
\end{array}
herbie shell --seed 2024144
(FPCore (wj x)
:name "Jmat.Real.lambertw, newton loop step"
:precision binary64
:alt
(! :herbie-platform default (let ((ew (exp wj))) (- wj (- (/ wj (+ wj 1)) (/ x (+ ew (* wj ew)))))))
(- wj (/ (- (* wj (exp wj)) x) (+ (exp wj) (* wj (exp wj))))))