
(FPCore (a b) :precision binary64 (/ (exp a) (+ (exp a) (exp b))))
double code(double a, double b) {
return exp(a) / (exp(a) + exp(b));
}
real(8) function code(a, b)
real(8), intent (in) :: a
real(8), intent (in) :: b
code = exp(a) / (exp(a) + exp(b))
end function
public static double code(double a, double b) {
return Math.exp(a) / (Math.exp(a) + Math.exp(b));
}
def code(a, b): return math.exp(a) / (math.exp(a) + math.exp(b))
function code(a, b) return Float64(exp(a) / Float64(exp(a) + exp(b))) end
function tmp = code(a, b) tmp = exp(a) / (exp(a) + exp(b)); end
code[a_, b_] := N[(N[Exp[a], $MachinePrecision] / N[(N[Exp[a], $MachinePrecision] + N[Exp[b], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{e^{a}}{e^{a} + e^{b}}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 17 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (a b) :precision binary64 (/ (exp a) (+ (exp a) (exp b))))
double code(double a, double b) {
return exp(a) / (exp(a) + exp(b));
}
real(8) function code(a, b)
real(8), intent (in) :: a
real(8), intent (in) :: b
code = exp(a) / (exp(a) + exp(b))
end function
public static double code(double a, double b) {
return Math.exp(a) / (Math.exp(a) + Math.exp(b));
}
def code(a, b): return math.exp(a) / (math.exp(a) + math.exp(b))
function code(a, b) return Float64(exp(a) / Float64(exp(a) + exp(b))) end
function tmp = code(a, b) tmp = exp(a) / (exp(a) + exp(b)); end
code[a_, b_] := N[(N[Exp[a], $MachinePrecision] / N[(N[Exp[a], $MachinePrecision] + N[Exp[b], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{e^{a}}{e^{a} + e^{b}}
\end{array}
(FPCore (a b)
:precision binary64
(let* ((t_0 (* a (+ 1.0 (* a (+ 0.5 (* a 0.16666666666666666)))))))
(if (<= (exp a) 0.9999)
(/ (exp a) (+ (exp a) 1.0))
(/ (+ 1.0 t_0) (+ t_0 (+ (exp b) 1.0))))))
double code(double a, double b) {
double t_0 = a * (1.0 + (a * (0.5 + (a * 0.16666666666666666))));
double tmp;
if (exp(a) <= 0.9999) {
tmp = exp(a) / (exp(a) + 1.0);
} else {
tmp = (1.0 + t_0) / (t_0 + (exp(b) + 1.0));
}
return tmp;
}
real(8) function code(a, b)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: t_0
real(8) :: tmp
t_0 = a * (1.0d0 + (a * (0.5d0 + (a * 0.16666666666666666d0))))
if (exp(a) <= 0.9999d0) then
tmp = exp(a) / (exp(a) + 1.0d0)
else
tmp = (1.0d0 + t_0) / (t_0 + (exp(b) + 1.0d0))
end if
code = tmp
end function
public static double code(double a, double b) {
double t_0 = a * (1.0 + (a * (0.5 + (a * 0.16666666666666666))));
double tmp;
if (Math.exp(a) <= 0.9999) {
tmp = Math.exp(a) / (Math.exp(a) + 1.0);
} else {
tmp = (1.0 + t_0) / (t_0 + (Math.exp(b) + 1.0));
}
return tmp;
}
def code(a, b): t_0 = a * (1.0 + (a * (0.5 + (a * 0.16666666666666666)))) tmp = 0 if math.exp(a) <= 0.9999: tmp = math.exp(a) / (math.exp(a) + 1.0) else: tmp = (1.0 + t_0) / (t_0 + (math.exp(b) + 1.0)) return tmp
function code(a, b) t_0 = Float64(a * Float64(1.0 + Float64(a * Float64(0.5 + Float64(a * 0.16666666666666666))))) tmp = 0.0 if (exp(a) <= 0.9999) tmp = Float64(exp(a) / Float64(exp(a) + 1.0)); else tmp = Float64(Float64(1.0 + t_0) / Float64(t_0 + Float64(exp(b) + 1.0))); end return tmp end
function tmp_2 = code(a, b) t_0 = a * (1.0 + (a * (0.5 + (a * 0.16666666666666666)))); tmp = 0.0; if (exp(a) <= 0.9999) tmp = exp(a) / (exp(a) + 1.0); else tmp = (1.0 + t_0) / (t_0 + (exp(b) + 1.0)); end tmp_2 = tmp; end
code[a_, b_] := Block[{t$95$0 = N[(a * N[(1.0 + N[(a * N[(0.5 + N[(a * 0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[N[Exp[a], $MachinePrecision], 0.9999], N[(N[Exp[a], $MachinePrecision] / N[(N[Exp[a], $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision], N[(N[(1.0 + t$95$0), $MachinePrecision] / N[(t$95$0 + N[(N[Exp[b], $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := a \cdot \left(1 + a \cdot \left(0.5 + a \cdot 0.16666666666666666\right)\right)\\
\mathbf{if}\;e^{a} \leq 0.9999:\\
\;\;\;\;\frac{e^{a}}{e^{a} + 1}\\
\mathbf{else}:\\
\;\;\;\;\frac{1 + t\_0}{t\_0 + \left(e^{b} + 1\right)}\\
\end{array}
\end{array}
if (exp.f64 a) < 0.99990000000000001Initial program 100.0%
Taylor expanded in b around 0
Simplified100.0%
if 0.99990000000000001 < (exp.f64 a) Initial program 98.9%
Taylor expanded in a around 0
associate-+r+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6497.9%
Simplified97.9%
Taylor expanded in a around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f6499.3%
Simplified99.3%
Final simplification99.5%
(FPCore (a b) :precision binary64 (/ (exp a) (+ (exp a) (exp b))))
double code(double a, double b) {
return exp(a) / (exp(a) + exp(b));
}
real(8) function code(a, b)
real(8), intent (in) :: a
real(8), intent (in) :: b
code = exp(a) / (exp(a) + exp(b))
end function
public static double code(double a, double b) {
return Math.exp(a) / (Math.exp(a) + Math.exp(b));
}
def code(a, b): return math.exp(a) / (math.exp(a) + math.exp(b))
function code(a, b) return Float64(exp(a) / Float64(exp(a) + exp(b))) end
function tmp = code(a, b) tmp = exp(a) / (exp(a) + exp(b)); end
code[a_, b_] := N[(N[Exp[a], $MachinePrecision] / N[(N[Exp[a], $MachinePrecision] + N[Exp[b], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{e^{a}}{e^{a} + e^{b}}
\end{array}
Initial program 99.2%
(FPCore (a b)
:precision binary64
(let* ((t_0 (* a (+ 1.0 (* a (+ 0.5 (* a 0.16666666666666666)))))))
(if (<= (exp a) 0.0)
(/ (exp a) 2.0)
(/ (+ 1.0 t_0) (+ t_0 (+ (exp b) 1.0))))))
double code(double a, double b) {
double t_0 = a * (1.0 + (a * (0.5 + (a * 0.16666666666666666))));
double tmp;
if (exp(a) <= 0.0) {
tmp = exp(a) / 2.0;
} else {
tmp = (1.0 + t_0) / (t_0 + (exp(b) + 1.0));
}
return tmp;
}
real(8) function code(a, b)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: t_0
real(8) :: tmp
t_0 = a * (1.0d0 + (a * (0.5d0 + (a * 0.16666666666666666d0))))
if (exp(a) <= 0.0d0) then
tmp = exp(a) / 2.0d0
else
tmp = (1.0d0 + t_0) / (t_0 + (exp(b) + 1.0d0))
end if
code = tmp
end function
public static double code(double a, double b) {
double t_0 = a * (1.0 + (a * (0.5 + (a * 0.16666666666666666))));
double tmp;
if (Math.exp(a) <= 0.0) {
tmp = Math.exp(a) / 2.0;
} else {
tmp = (1.0 + t_0) / (t_0 + (Math.exp(b) + 1.0));
}
return tmp;
}
def code(a, b): t_0 = a * (1.0 + (a * (0.5 + (a * 0.16666666666666666)))) tmp = 0 if math.exp(a) <= 0.0: tmp = math.exp(a) / 2.0 else: tmp = (1.0 + t_0) / (t_0 + (math.exp(b) + 1.0)) return tmp
function code(a, b) t_0 = Float64(a * Float64(1.0 + Float64(a * Float64(0.5 + Float64(a * 0.16666666666666666))))) tmp = 0.0 if (exp(a) <= 0.0) tmp = Float64(exp(a) / 2.0); else tmp = Float64(Float64(1.0 + t_0) / Float64(t_0 + Float64(exp(b) + 1.0))); end return tmp end
function tmp_2 = code(a, b) t_0 = a * (1.0 + (a * (0.5 + (a * 0.16666666666666666)))); tmp = 0.0; if (exp(a) <= 0.0) tmp = exp(a) / 2.0; else tmp = (1.0 + t_0) / (t_0 + (exp(b) + 1.0)); end tmp_2 = tmp; end
code[a_, b_] := Block[{t$95$0 = N[(a * N[(1.0 + N[(a * N[(0.5 + N[(a * 0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[N[Exp[a], $MachinePrecision], 0.0], N[(N[Exp[a], $MachinePrecision] / 2.0), $MachinePrecision], N[(N[(1.0 + t$95$0), $MachinePrecision] / N[(t$95$0 + N[(N[Exp[b], $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := a \cdot \left(1 + a \cdot \left(0.5 + a \cdot 0.16666666666666666\right)\right)\\
\mathbf{if}\;e^{a} \leq 0:\\
\;\;\;\;\frac{e^{a}}{2}\\
\mathbf{else}:\\
\;\;\;\;\frac{1 + t\_0}{t\_0 + \left(e^{b} + 1\right)}\\
\end{array}
\end{array}
if (exp.f64 a) < 0.0Initial program 100.0%
Taylor expanded in b around 0
Simplified100.0%
Taylor expanded in a around 0
Simplified100.0%
if 0.0 < (exp.f64 a) Initial program 98.9%
Taylor expanded in a around 0
associate-+r+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6497.9%
Simplified97.9%
Taylor expanded in a around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f6499.3%
Simplified99.3%
Final simplification99.5%
(FPCore (a b) :precision binary64 (if (<= (exp a) 0.99996) (/ (exp a) 2.0) (/ 1.0 (+ (exp b) 1.0))))
double code(double a, double b) {
double tmp;
if (exp(a) <= 0.99996) {
tmp = exp(a) / 2.0;
} else {
tmp = 1.0 / (exp(b) + 1.0);
}
return tmp;
}
real(8) function code(a, b)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: tmp
if (exp(a) <= 0.99996d0) then
tmp = exp(a) / 2.0d0
else
tmp = 1.0d0 / (exp(b) + 1.0d0)
end if
code = tmp
end function
public static double code(double a, double b) {
double tmp;
if (Math.exp(a) <= 0.99996) {
tmp = Math.exp(a) / 2.0;
} else {
tmp = 1.0 / (Math.exp(b) + 1.0);
}
return tmp;
}
def code(a, b): tmp = 0 if math.exp(a) <= 0.99996: tmp = math.exp(a) / 2.0 else: tmp = 1.0 / (math.exp(b) + 1.0) return tmp
function code(a, b) tmp = 0.0 if (exp(a) <= 0.99996) tmp = Float64(exp(a) / 2.0); else tmp = Float64(1.0 / Float64(exp(b) + 1.0)); end return tmp end
function tmp_2 = code(a, b) tmp = 0.0; if (exp(a) <= 0.99996) tmp = exp(a) / 2.0; else tmp = 1.0 / (exp(b) + 1.0); end tmp_2 = tmp; end
code[a_, b_] := If[LessEqual[N[Exp[a], $MachinePrecision], 0.99996], N[(N[Exp[a], $MachinePrecision] / 2.0), $MachinePrecision], N[(1.0 / N[(N[Exp[b], $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;e^{a} \leq 0.99996:\\
\;\;\;\;\frac{e^{a}}{2}\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{e^{b} + 1}\\
\end{array}
\end{array}
if (exp.f64 a) < 0.99995999999999996Initial program 100.0%
Taylor expanded in b around 0
Simplified100.0%
Taylor expanded in a around 0
Simplified98.3%
if 0.99995999999999996 < (exp.f64 a) Initial program 98.9%
Taylor expanded in a around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
exp-lowering-exp.f6497.7%
Simplified97.7%
Final simplification97.9%
(FPCore (a b)
:precision binary64
(let* ((t_0 (* b (+ 0.5 (* b 0.16666666666666666)))))
(if (<= b -17.5)
1.0
(if (<= b 4.8e+45)
(/ (exp a) 2.0)
(if (<= b 5e+102)
(/
(+ 2.0 (* b (+ -1.0 (* b (+ -0.5 (* b -0.16666666666666666))))))
(+ 4.0 (* (* b (+ 1.0 t_0)) (* b (- -1.0 t_0)))))
(/ -4.0 (* b (* b b))))))))
double code(double a, double b) {
double t_0 = b * (0.5 + (b * 0.16666666666666666));
double tmp;
if (b <= -17.5) {
tmp = 1.0;
} else if (b <= 4.8e+45) {
tmp = exp(a) / 2.0;
} else if (b <= 5e+102) {
tmp = (2.0 + (b * (-1.0 + (b * (-0.5 + (b * -0.16666666666666666)))))) / (4.0 + ((b * (1.0 + t_0)) * (b * (-1.0 - t_0))));
} else {
tmp = -4.0 / (b * (b * b));
}
return tmp;
}
real(8) function code(a, b)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: t_0
real(8) :: tmp
t_0 = b * (0.5d0 + (b * 0.16666666666666666d0))
if (b <= (-17.5d0)) then
tmp = 1.0d0
else if (b <= 4.8d+45) then
tmp = exp(a) / 2.0d0
else if (b <= 5d+102) then
tmp = (2.0d0 + (b * ((-1.0d0) + (b * ((-0.5d0) + (b * (-0.16666666666666666d0))))))) / (4.0d0 + ((b * (1.0d0 + t_0)) * (b * ((-1.0d0) - t_0))))
else
tmp = (-4.0d0) / (b * (b * b))
end if
code = tmp
end function
public static double code(double a, double b) {
double t_0 = b * (0.5 + (b * 0.16666666666666666));
double tmp;
if (b <= -17.5) {
tmp = 1.0;
} else if (b <= 4.8e+45) {
tmp = Math.exp(a) / 2.0;
} else if (b <= 5e+102) {
tmp = (2.0 + (b * (-1.0 + (b * (-0.5 + (b * -0.16666666666666666)))))) / (4.0 + ((b * (1.0 + t_0)) * (b * (-1.0 - t_0))));
} else {
tmp = -4.0 / (b * (b * b));
}
return tmp;
}
def code(a, b): t_0 = b * (0.5 + (b * 0.16666666666666666)) tmp = 0 if b <= -17.5: tmp = 1.0 elif b <= 4.8e+45: tmp = math.exp(a) / 2.0 elif b <= 5e+102: tmp = (2.0 + (b * (-1.0 + (b * (-0.5 + (b * -0.16666666666666666)))))) / (4.0 + ((b * (1.0 + t_0)) * (b * (-1.0 - t_0)))) else: tmp = -4.0 / (b * (b * b)) return tmp
function code(a, b) t_0 = Float64(b * Float64(0.5 + Float64(b * 0.16666666666666666))) tmp = 0.0 if (b <= -17.5) tmp = 1.0; elseif (b <= 4.8e+45) tmp = Float64(exp(a) / 2.0); elseif (b <= 5e+102) tmp = Float64(Float64(2.0 + Float64(b * Float64(-1.0 + Float64(b * Float64(-0.5 + Float64(b * -0.16666666666666666)))))) / Float64(4.0 + Float64(Float64(b * Float64(1.0 + t_0)) * Float64(b * Float64(-1.0 - t_0))))); else tmp = Float64(-4.0 / Float64(b * Float64(b * b))); end return tmp end
function tmp_2 = code(a, b) t_0 = b * (0.5 + (b * 0.16666666666666666)); tmp = 0.0; if (b <= -17.5) tmp = 1.0; elseif (b <= 4.8e+45) tmp = exp(a) / 2.0; elseif (b <= 5e+102) tmp = (2.0 + (b * (-1.0 + (b * (-0.5 + (b * -0.16666666666666666)))))) / (4.0 + ((b * (1.0 + t_0)) * (b * (-1.0 - t_0)))); else tmp = -4.0 / (b * (b * b)); end tmp_2 = tmp; end
code[a_, b_] := Block[{t$95$0 = N[(b * N[(0.5 + N[(b * 0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[b, -17.5], 1.0, If[LessEqual[b, 4.8e+45], N[(N[Exp[a], $MachinePrecision] / 2.0), $MachinePrecision], If[LessEqual[b, 5e+102], N[(N[(2.0 + N[(b * N[(-1.0 + N[(b * N[(-0.5 + N[(b * -0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(4.0 + N[(N[(b * N[(1.0 + t$95$0), $MachinePrecision]), $MachinePrecision] * N[(b * N[(-1.0 - t$95$0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(-4.0 / N[(b * N[(b * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := b \cdot \left(0.5 + b \cdot 0.16666666666666666\right)\\
\mathbf{if}\;b \leq -17.5:\\
\;\;\;\;1\\
\mathbf{elif}\;b \leq 4.8 \cdot 10^{+45}:\\
\;\;\;\;\frac{e^{a}}{2}\\
\mathbf{elif}\;b \leq 5 \cdot 10^{+102}:\\
\;\;\;\;\frac{2 + b \cdot \left(-1 + b \cdot \left(-0.5 + b \cdot -0.16666666666666666\right)\right)}{4 + \left(b \cdot \left(1 + t\_0\right)\right) \cdot \left(b \cdot \left(-1 - t\_0\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{-4}{b \cdot \left(b \cdot b\right)}\\
\end{array}
\end{array}
if b < -17.5Initial program 100.0%
Taylor expanded in a around 0
associate-+r+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6499.8%
Simplified99.8%
Taylor expanded in a around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64100.0%
Simplified100.0%
Taylor expanded in a around inf
Simplified100.0%
if -17.5 < b < 4.79999999999999979e45Initial program 99.3%
Taylor expanded in b around 0
Simplified95.2%
Taylor expanded in a around 0
Simplified92.2%
if 4.79999999999999979e45 < b < 5e102Initial program 100.0%
Taylor expanded in a around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
exp-lowering-exp.f64100.0%
Simplified100.0%
Taylor expanded in b around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f645.6%
Simplified5.6%
flip-+N/A
clear-numN/A
/-lowering-/.f64N/A
Applied egg-rr89.4%
if 5e102 < b Initial program 97.6%
Taylor expanded in a around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
exp-lowering-exp.f64100.0%
Simplified100.0%
Taylor expanded in b around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f6468.6%
Simplified68.6%
Taylor expanded in b around inf
/-lowering-/.f64N/A
sub-negN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
distribute-neg-fracN/A
/-lowering-/.f64N/A
metadata-evalN/A
unpow2N/A
*-lowering-*.f6468.6%
Simplified68.6%
Taylor expanded in b around 0
/-lowering-/.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64100.0%
Simplified100.0%
Final simplification94.8%
(FPCore (a b)
:precision binary64
(let* ((t_0 (+ 0.5 (* b 0.16666666666666666))) (t_1 (* b t_0)))
(if (<= b -1.55)
1.0
(if (<= b 2.8e+77)
(/
1.0
(+
2.0
(/
(* b (+ 1.0 (* (* b t_1) (* b (* t_0 t_0)))))
(+ 1.0 (* t_1 (+ -1.0 t_1))))))
(/ -4.0 (* b (* b b)))))))
double code(double a, double b) {
double t_0 = 0.5 + (b * 0.16666666666666666);
double t_1 = b * t_0;
double tmp;
if (b <= -1.55) {
tmp = 1.0;
} else if (b <= 2.8e+77) {
tmp = 1.0 / (2.0 + ((b * (1.0 + ((b * t_1) * (b * (t_0 * t_0))))) / (1.0 + (t_1 * (-1.0 + t_1)))));
} else {
tmp = -4.0 / (b * (b * b));
}
return tmp;
}
real(8) function code(a, b)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: t_0
real(8) :: t_1
real(8) :: tmp
t_0 = 0.5d0 + (b * 0.16666666666666666d0)
t_1 = b * t_0
if (b <= (-1.55d0)) then
tmp = 1.0d0
else if (b <= 2.8d+77) then
tmp = 1.0d0 / (2.0d0 + ((b * (1.0d0 + ((b * t_1) * (b * (t_0 * t_0))))) / (1.0d0 + (t_1 * ((-1.0d0) + t_1)))))
else
tmp = (-4.0d0) / (b * (b * b))
end if
code = tmp
end function
public static double code(double a, double b) {
double t_0 = 0.5 + (b * 0.16666666666666666);
double t_1 = b * t_0;
double tmp;
if (b <= -1.55) {
tmp = 1.0;
} else if (b <= 2.8e+77) {
tmp = 1.0 / (2.0 + ((b * (1.0 + ((b * t_1) * (b * (t_0 * t_0))))) / (1.0 + (t_1 * (-1.0 + t_1)))));
} else {
tmp = -4.0 / (b * (b * b));
}
return tmp;
}
def code(a, b): t_0 = 0.5 + (b * 0.16666666666666666) t_1 = b * t_0 tmp = 0 if b <= -1.55: tmp = 1.0 elif b <= 2.8e+77: tmp = 1.0 / (2.0 + ((b * (1.0 + ((b * t_1) * (b * (t_0 * t_0))))) / (1.0 + (t_1 * (-1.0 + t_1))))) else: tmp = -4.0 / (b * (b * b)) return tmp
function code(a, b) t_0 = Float64(0.5 + Float64(b * 0.16666666666666666)) t_1 = Float64(b * t_0) tmp = 0.0 if (b <= -1.55) tmp = 1.0; elseif (b <= 2.8e+77) tmp = Float64(1.0 / Float64(2.0 + Float64(Float64(b * Float64(1.0 + Float64(Float64(b * t_1) * Float64(b * Float64(t_0 * t_0))))) / Float64(1.0 + Float64(t_1 * Float64(-1.0 + t_1)))))); else tmp = Float64(-4.0 / Float64(b * Float64(b * b))); end return tmp end
function tmp_2 = code(a, b) t_0 = 0.5 + (b * 0.16666666666666666); t_1 = b * t_0; tmp = 0.0; if (b <= -1.55) tmp = 1.0; elseif (b <= 2.8e+77) tmp = 1.0 / (2.0 + ((b * (1.0 + ((b * t_1) * (b * (t_0 * t_0))))) / (1.0 + (t_1 * (-1.0 + t_1))))); else tmp = -4.0 / (b * (b * b)); end tmp_2 = tmp; end
code[a_, b_] := Block[{t$95$0 = N[(0.5 + N[(b * 0.16666666666666666), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(b * t$95$0), $MachinePrecision]}, If[LessEqual[b, -1.55], 1.0, If[LessEqual[b, 2.8e+77], N[(1.0 / N[(2.0 + N[(N[(b * N[(1.0 + N[(N[(b * t$95$1), $MachinePrecision] * N[(b * N[(t$95$0 * t$95$0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(1.0 + N[(t$95$1 * N[(-1.0 + t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(-4.0 / N[(b * N[(b * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := 0.5 + b \cdot 0.16666666666666666\\
t_1 := b \cdot t\_0\\
\mathbf{if}\;b \leq -1.55:\\
\;\;\;\;1\\
\mathbf{elif}\;b \leq 2.8 \cdot 10^{+77}:\\
\;\;\;\;\frac{1}{2 + \frac{b \cdot \left(1 + \left(b \cdot t\_1\right) \cdot \left(b \cdot \left(t\_0 \cdot t\_0\right)\right)\right)}{1 + t\_1 \cdot \left(-1 + t\_1\right)}}\\
\mathbf{else}:\\
\;\;\;\;\frac{-4}{b \cdot \left(b \cdot b\right)}\\
\end{array}
\end{array}
if b < -1.55000000000000004Initial program 100.0%
Taylor expanded in a around 0
associate-+r+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6499.8%
Simplified99.8%
Taylor expanded in a around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64100.0%
Simplified100.0%
Taylor expanded in a around inf
Simplified100.0%
if -1.55000000000000004 < b < 2.8e77Initial program 99.4%
Taylor expanded in a around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
exp-lowering-exp.f6466.6%
Simplified66.6%
Taylor expanded in b around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6457.9%
Simplified57.9%
*-commutativeN/A
flip3-+N/A
associate-*l/N/A
/-lowering-/.f64N/A
Applied egg-rr62.5%
if 2.8e77 < b Initial program 97.7%
Taylor expanded in a around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
exp-lowering-exp.f64100.0%
Simplified100.0%
Taylor expanded in b around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f6464.2%
Simplified64.2%
Taylor expanded in b around inf
/-lowering-/.f64N/A
sub-negN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
distribute-neg-fracN/A
/-lowering-/.f64N/A
metadata-evalN/A
unpow2N/A
*-lowering-*.f6464.2%
Simplified64.2%
Taylor expanded in b around 0
/-lowering-/.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6493.6%
Simplified93.6%
Final simplification74.9%
(FPCore (a b)
:precision binary64
(let* ((t_0 (* b (+ 0.5 (* b 0.16666666666666666)))))
(if (<= b -1.25)
1.0
(if (<= b 5e+102)
(/
(+ 2.0 (* b (+ -1.0 (* b (+ -0.5 (* b -0.16666666666666666))))))
(+ 4.0 (* (* b (+ 1.0 t_0)) (* b (- -1.0 t_0)))))
(/ -4.0 (* b (* b b)))))))
double code(double a, double b) {
double t_0 = b * (0.5 + (b * 0.16666666666666666));
double tmp;
if (b <= -1.25) {
tmp = 1.0;
} else if (b <= 5e+102) {
tmp = (2.0 + (b * (-1.0 + (b * (-0.5 + (b * -0.16666666666666666)))))) / (4.0 + ((b * (1.0 + t_0)) * (b * (-1.0 - t_0))));
} else {
tmp = -4.0 / (b * (b * b));
}
return tmp;
}
real(8) function code(a, b)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: t_0
real(8) :: tmp
t_0 = b * (0.5d0 + (b * 0.16666666666666666d0))
if (b <= (-1.25d0)) then
tmp = 1.0d0
else if (b <= 5d+102) then
tmp = (2.0d0 + (b * ((-1.0d0) + (b * ((-0.5d0) + (b * (-0.16666666666666666d0))))))) / (4.0d0 + ((b * (1.0d0 + t_0)) * (b * ((-1.0d0) - t_0))))
else
tmp = (-4.0d0) / (b * (b * b))
end if
code = tmp
end function
public static double code(double a, double b) {
double t_0 = b * (0.5 + (b * 0.16666666666666666));
double tmp;
if (b <= -1.25) {
tmp = 1.0;
} else if (b <= 5e+102) {
tmp = (2.0 + (b * (-1.0 + (b * (-0.5 + (b * -0.16666666666666666)))))) / (4.0 + ((b * (1.0 + t_0)) * (b * (-1.0 - t_0))));
} else {
tmp = -4.0 / (b * (b * b));
}
return tmp;
}
def code(a, b): t_0 = b * (0.5 + (b * 0.16666666666666666)) tmp = 0 if b <= -1.25: tmp = 1.0 elif b <= 5e+102: tmp = (2.0 + (b * (-1.0 + (b * (-0.5 + (b * -0.16666666666666666)))))) / (4.0 + ((b * (1.0 + t_0)) * (b * (-1.0 - t_0)))) else: tmp = -4.0 / (b * (b * b)) return tmp
function code(a, b) t_0 = Float64(b * Float64(0.5 + Float64(b * 0.16666666666666666))) tmp = 0.0 if (b <= -1.25) tmp = 1.0; elseif (b <= 5e+102) tmp = Float64(Float64(2.0 + Float64(b * Float64(-1.0 + Float64(b * Float64(-0.5 + Float64(b * -0.16666666666666666)))))) / Float64(4.0 + Float64(Float64(b * Float64(1.0 + t_0)) * Float64(b * Float64(-1.0 - t_0))))); else tmp = Float64(-4.0 / Float64(b * Float64(b * b))); end return tmp end
function tmp_2 = code(a, b) t_0 = b * (0.5 + (b * 0.16666666666666666)); tmp = 0.0; if (b <= -1.25) tmp = 1.0; elseif (b <= 5e+102) tmp = (2.0 + (b * (-1.0 + (b * (-0.5 + (b * -0.16666666666666666)))))) / (4.0 + ((b * (1.0 + t_0)) * (b * (-1.0 - t_0)))); else tmp = -4.0 / (b * (b * b)); end tmp_2 = tmp; end
code[a_, b_] := Block[{t$95$0 = N[(b * N[(0.5 + N[(b * 0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[b, -1.25], 1.0, If[LessEqual[b, 5e+102], N[(N[(2.0 + N[(b * N[(-1.0 + N[(b * N[(-0.5 + N[(b * -0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(4.0 + N[(N[(b * N[(1.0 + t$95$0), $MachinePrecision]), $MachinePrecision] * N[(b * N[(-1.0 - t$95$0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(-4.0 / N[(b * N[(b * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := b \cdot \left(0.5 + b \cdot 0.16666666666666666\right)\\
\mathbf{if}\;b \leq -1.25:\\
\;\;\;\;1\\
\mathbf{elif}\;b \leq 5 \cdot 10^{+102}:\\
\;\;\;\;\frac{2 + b \cdot \left(-1 + b \cdot \left(-0.5 + b \cdot -0.16666666666666666\right)\right)}{4 + \left(b \cdot \left(1 + t\_0\right)\right) \cdot \left(b \cdot \left(-1 - t\_0\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{-4}{b \cdot \left(b \cdot b\right)}\\
\end{array}
\end{array}
if b < -1.25Initial program 100.0%
Taylor expanded in a around 0
associate-+r+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6499.8%
Simplified99.8%
Taylor expanded in a around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64100.0%
Simplified100.0%
Taylor expanded in a around inf
Simplified100.0%
if -1.25 < b < 5e102Initial program 99.4%
Taylor expanded in a around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
exp-lowering-exp.f6467.2%
Simplified67.2%
Taylor expanded in b around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6457.0%
Simplified57.0%
flip-+N/A
clear-numN/A
/-lowering-/.f64N/A
Applied egg-rr61.5%
if 5e102 < b Initial program 97.6%
Taylor expanded in a around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
exp-lowering-exp.f64100.0%
Simplified100.0%
Taylor expanded in b around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f6468.6%
Simplified68.6%
Taylor expanded in b around inf
/-lowering-/.f64N/A
sub-negN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
distribute-neg-fracN/A
/-lowering-/.f64N/A
metadata-evalN/A
unpow2N/A
*-lowering-*.f6468.6%
Simplified68.6%
Taylor expanded in b around 0
/-lowering-/.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64100.0%
Simplified100.0%
Final simplification74.9%
(FPCore (a b)
:precision binary64
(if (<= b -2.4)
1.0
(if (<= b 420.0)
(+ 0.5 (* b (+ -0.25 (* b (* b 0.020833333333333332)))))
(if (<= b 1.08e+80)
(* a (* (* a a) -0.020833333333333332))
(/ -4.0 (* b (* b b)))))))
double code(double a, double b) {
double tmp;
if (b <= -2.4) {
tmp = 1.0;
} else if (b <= 420.0) {
tmp = 0.5 + (b * (-0.25 + (b * (b * 0.020833333333333332))));
} else if (b <= 1.08e+80) {
tmp = a * ((a * a) * -0.020833333333333332);
} else {
tmp = -4.0 / (b * (b * b));
}
return tmp;
}
real(8) function code(a, b)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: tmp
if (b <= (-2.4d0)) then
tmp = 1.0d0
else if (b <= 420.0d0) then
tmp = 0.5d0 + (b * ((-0.25d0) + (b * (b * 0.020833333333333332d0))))
else if (b <= 1.08d+80) then
tmp = a * ((a * a) * (-0.020833333333333332d0))
else
tmp = (-4.0d0) / (b * (b * b))
end if
code = tmp
end function
public static double code(double a, double b) {
double tmp;
if (b <= -2.4) {
tmp = 1.0;
} else if (b <= 420.0) {
tmp = 0.5 + (b * (-0.25 + (b * (b * 0.020833333333333332))));
} else if (b <= 1.08e+80) {
tmp = a * ((a * a) * -0.020833333333333332);
} else {
tmp = -4.0 / (b * (b * b));
}
return tmp;
}
def code(a, b): tmp = 0 if b <= -2.4: tmp = 1.0 elif b <= 420.0: tmp = 0.5 + (b * (-0.25 + (b * (b * 0.020833333333333332)))) elif b <= 1.08e+80: tmp = a * ((a * a) * -0.020833333333333332) else: tmp = -4.0 / (b * (b * b)) return tmp
function code(a, b) tmp = 0.0 if (b <= -2.4) tmp = 1.0; elseif (b <= 420.0) tmp = Float64(0.5 + Float64(b * Float64(-0.25 + Float64(b * Float64(b * 0.020833333333333332))))); elseif (b <= 1.08e+80) tmp = Float64(a * Float64(Float64(a * a) * -0.020833333333333332)); else tmp = Float64(-4.0 / Float64(b * Float64(b * b))); end return tmp end
function tmp_2 = code(a, b) tmp = 0.0; if (b <= -2.4) tmp = 1.0; elseif (b <= 420.0) tmp = 0.5 + (b * (-0.25 + (b * (b * 0.020833333333333332)))); elseif (b <= 1.08e+80) tmp = a * ((a * a) * -0.020833333333333332); else tmp = -4.0 / (b * (b * b)); end tmp_2 = tmp; end
code[a_, b_] := If[LessEqual[b, -2.4], 1.0, If[LessEqual[b, 420.0], N[(0.5 + N[(b * N[(-0.25 + N[(b * N[(b * 0.020833333333333332), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 1.08e+80], N[(a * N[(N[(a * a), $MachinePrecision] * -0.020833333333333332), $MachinePrecision]), $MachinePrecision], N[(-4.0 / N[(b * N[(b * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq -2.4:\\
\;\;\;\;1\\
\mathbf{elif}\;b \leq 420:\\
\;\;\;\;0.5 + b \cdot \left(-0.25 + b \cdot \left(b \cdot 0.020833333333333332\right)\right)\\
\mathbf{elif}\;b \leq 1.08 \cdot 10^{+80}:\\
\;\;\;\;a \cdot \left(\left(a \cdot a\right) \cdot -0.020833333333333332\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{-4}{b \cdot \left(b \cdot b\right)}\\
\end{array}
\end{array}
if b < -2.39999999999999991Initial program 100.0%
Taylor expanded in a around 0
associate-+r+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6499.8%
Simplified99.8%
Taylor expanded in a around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64100.0%
Simplified100.0%
Taylor expanded in a around inf
Simplified100.0%
if -2.39999999999999991 < b < 420Initial program 99.3%
Taylor expanded in a around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
exp-lowering-exp.f6463.5%
Simplified63.5%
Taylor expanded in b around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f6463.1%
Simplified63.1%
if 420 < b < 1.08e80Initial program 100.0%
Taylor expanded in b around 0
Simplified35.5%
Taylor expanded in a around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f642.7%
Simplified2.7%
Taylor expanded in a around inf
unpow3N/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6435.7%
Simplified35.7%
if 1.08e80 < b Initial program 97.7%
Taylor expanded in a around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
exp-lowering-exp.f64100.0%
Simplified100.0%
Taylor expanded in b around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f6465.6%
Simplified65.6%
Taylor expanded in b around inf
/-lowering-/.f64N/A
sub-negN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
distribute-neg-fracN/A
/-lowering-/.f64N/A
metadata-evalN/A
unpow2N/A
*-lowering-*.f6465.6%
Simplified65.6%
Taylor expanded in b around 0
/-lowering-/.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6495.7%
Simplified95.7%
(FPCore (a b)
:precision binary64
(if (<= b -0.061)
1.0
(if (<= b 350.0)
(+ 0.5 (* a 0.25))
(if (<= b 1.08e+80)
(* a (* (* a a) -0.020833333333333332))
(/ -4.0 (* b (* b b)))))))
double code(double a, double b) {
double tmp;
if (b <= -0.061) {
tmp = 1.0;
} else if (b <= 350.0) {
tmp = 0.5 + (a * 0.25);
} else if (b <= 1.08e+80) {
tmp = a * ((a * a) * -0.020833333333333332);
} else {
tmp = -4.0 / (b * (b * b));
}
return tmp;
}
real(8) function code(a, b)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: tmp
if (b <= (-0.061d0)) then
tmp = 1.0d0
else if (b <= 350.0d0) then
tmp = 0.5d0 + (a * 0.25d0)
else if (b <= 1.08d+80) then
tmp = a * ((a * a) * (-0.020833333333333332d0))
else
tmp = (-4.0d0) / (b * (b * b))
end if
code = tmp
end function
public static double code(double a, double b) {
double tmp;
if (b <= -0.061) {
tmp = 1.0;
} else if (b <= 350.0) {
tmp = 0.5 + (a * 0.25);
} else if (b <= 1.08e+80) {
tmp = a * ((a * a) * -0.020833333333333332);
} else {
tmp = -4.0 / (b * (b * b));
}
return tmp;
}
def code(a, b): tmp = 0 if b <= -0.061: tmp = 1.0 elif b <= 350.0: tmp = 0.5 + (a * 0.25) elif b <= 1.08e+80: tmp = a * ((a * a) * -0.020833333333333332) else: tmp = -4.0 / (b * (b * b)) return tmp
function code(a, b) tmp = 0.0 if (b <= -0.061) tmp = 1.0; elseif (b <= 350.0) tmp = Float64(0.5 + Float64(a * 0.25)); elseif (b <= 1.08e+80) tmp = Float64(a * Float64(Float64(a * a) * -0.020833333333333332)); else tmp = Float64(-4.0 / Float64(b * Float64(b * b))); end return tmp end
function tmp_2 = code(a, b) tmp = 0.0; if (b <= -0.061) tmp = 1.0; elseif (b <= 350.0) tmp = 0.5 + (a * 0.25); elseif (b <= 1.08e+80) tmp = a * ((a * a) * -0.020833333333333332); else tmp = -4.0 / (b * (b * b)); end tmp_2 = tmp; end
code[a_, b_] := If[LessEqual[b, -0.061], 1.0, If[LessEqual[b, 350.0], N[(0.5 + N[(a * 0.25), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 1.08e+80], N[(a * N[(N[(a * a), $MachinePrecision] * -0.020833333333333332), $MachinePrecision]), $MachinePrecision], N[(-4.0 / N[(b * N[(b * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq -0.061:\\
\;\;\;\;1\\
\mathbf{elif}\;b \leq 350:\\
\;\;\;\;0.5 + a \cdot 0.25\\
\mathbf{elif}\;b \leq 1.08 \cdot 10^{+80}:\\
\;\;\;\;a \cdot \left(\left(a \cdot a\right) \cdot -0.020833333333333332\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{-4}{b \cdot \left(b \cdot b\right)}\\
\end{array}
\end{array}
if b < -0.060999999999999999Initial program 100.0%
Taylor expanded in a around 0
associate-+r+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6499.8%
Simplified99.8%
Taylor expanded in a around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f6498.0%
Simplified98.0%
Taylor expanded in a around inf
Simplified98.0%
if -0.060999999999999999 < b < 350Initial program 99.3%
Taylor expanded in b around 0
Simplified97.5%
Taylor expanded in a around 0
+-lowering-+.f64N/A
*-lowering-*.f6463.4%
Simplified63.4%
if 350 < b < 1.08e80Initial program 100.0%
Taylor expanded in b around 0
Simplified35.5%
Taylor expanded in a around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f642.7%
Simplified2.7%
Taylor expanded in a around inf
unpow3N/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6435.7%
Simplified35.7%
if 1.08e80 < b Initial program 97.7%
Taylor expanded in a around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
exp-lowering-exp.f64100.0%
Simplified100.0%
Taylor expanded in b around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f6465.6%
Simplified65.6%
Taylor expanded in b around inf
/-lowering-/.f64N/A
sub-negN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
distribute-neg-fracN/A
/-lowering-/.f64N/A
metadata-evalN/A
unpow2N/A
*-lowering-*.f6465.6%
Simplified65.6%
Taylor expanded in b around 0
/-lowering-/.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6495.7%
Simplified95.7%
Final simplification73.8%
(FPCore (a b)
:precision binary64
(if (<= b -0.075)
1.0
(if (<= b 410.0)
(+ 0.5 (* a 0.25))
(if (<= b 1.85e+151)
(* a (* (* a a) -0.020833333333333332))
(/ 2.0 (* b b))))))
double code(double a, double b) {
double tmp;
if (b <= -0.075) {
tmp = 1.0;
} else if (b <= 410.0) {
tmp = 0.5 + (a * 0.25);
} else if (b <= 1.85e+151) {
tmp = a * ((a * a) * -0.020833333333333332);
} else {
tmp = 2.0 / (b * b);
}
return tmp;
}
real(8) function code(a, b)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: tmp
if (b <= (-0.075d0)) then
tmp = 1.0d0
else if (b <= 410.0d0) then
tmp = 0.5d0 + (a * 0.25d0)
else if (b <= 1.85d+151) then
tmp = a * ((a * a) * (-0.020833333333333332d0))
else
tmp = 2.0d0 / (b * b)
end if
code = tmp
end function
public static double code(double a, double b) {
double tmp;
if (b <= -0.075) {
tmp = 1.0;
} else if (b <= 410.0) {
tmp = 0.5 + (a * 0.25);
} else if (b <= 1.85e+151) {
tmp = a * ((a * a) * -0.020833333333333332);
} else {
tmp = 2.0 / (b * b);
}
return tmp;
}
def code(a, b): tmp = 0 if b <= -0.075: tmp = 1.0 elif b <= 410.0: tmp = 0.5 + (a * 0.25) elif b <= 1.85e+151: tmp = a * ((a * a) * -0.020833333333333332) else: tmp = 2.0 / (b * b) return tmp
function code(a, b) tmp = 0.0 if (b <= -0.075) tmp = 1.0; elseif (b <= 410.0) tmp = Float64(0.5 + Float64(a * 0.25)); elseif (b <= 1.85e+151) tmp = Float64(a * Float64(Float64(a * a) * -0.020833333333333332)); else tmp = Float64(2.0 / Float64(b * b)); end return tmp end
function tmp_2 = code(a, b) tmp = 0.0; if (b <= -0.075) tmp = 1.0; elseif (b <= 410.0) tmp = 0.5 + (a * 0.25); elseif (b <= 1.85e+151) tmp = a * ((a * a) * -0.020833333333333332); else tmp = 2.0 / (b * b); end tmp_2 = tmp; end
code[a_, b_] := If[LessEqual[b, -0.075], 1.0, If[LessEqual[b, 410.0], N[(0.5 + N[(a * 0.25), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 1.85e+151], N[(a * N[(N[(a * a), $MachinePrecision] * -0.020833333333333332), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(b * b), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq -0.075:\\
\;\;\;\;1\\
\mathbf{elif}\;b \leq 410:\\
\;\;\;\;0.5 + a \cdot 0.25\\
\mathbf{elif}\;b \leq 1.85 \cdot 10^{+151}:\\
\;\;\;\;a \cdot \left(\left(a \cdot a\right) \cdot -0.020833333333333332\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{b \cdot b}\\
\end{array}
\end{array}
if b < -0.0749999999999999972Initial program 100.0%
Taylor expanded in a around 0
associate-+r+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6499.8%
Simplified99.8%
Taylor expanded in a around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f6498.0%
Simplified98.0%
Taylor expanded in a around inf
Simplified98.0%
if -0.0749999999999999972 < b < 410Initial program 99.3%
Taylor expanded in b around 0
Simplified97.5%
Taylor expanded in a around 0
+-lowering-+.f64N/A
*-lowering-*.f6463.4%
Simplified63.4%
if 410 < b < 1.8499999999999999e151Initial program 100.0%
Taylor expanded in b around 0
Simplified33.2%
Taylor expanded in a around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f642.7%
Simplified2.7%
Taylor expanded in a around inf
unpow3N/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6438.9%
Simplified38.9%
if 1.8499999999999999e151 < b Initial program 96.6%
Taylor expanded in a around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
exp-lowering-exp.f64100.0%
Simplified100.0%
Taylor expanded in b around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f6494.0%
Simplified94.0%
Taylor expanded in b around inf
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6494.0%
Simplified94.0%
Final simplification70.7%
(FPCore (a b) :precision binary64 (if (<= b -1.55) 1.0 (/ 1.0 (+ 2.0 (* b (+ 1.0 (* b (+ 0.5 (* b 0.16666666666666666)))))))))
double code(double a, double b) {
double tmp;
if (b <= -1.55) {
tmp = 1.0;
} else {
tmp = 1.0 / (2.0 + (b * (1.0 + (b * (0.5 + (b * 0.16666666666666666))))));
}
return tmp;
}
real(8) function code(a, b)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: tmp
if (b <= (-1.55d0)) then
tmp = 1.0d0
else
tmp = 1.0d0 / (2.0d0 + (b * (1.0d0 + (b * (0.5d0 + (b * 0.16666666666666666d0))))))
end if
code = tmp
end function
public static double code(double a, double b) {
double tmp;
if (b <= -1.55) {
tmp = 1.0;
} else {
tmp = 1.0 / (2.0 + (b * (1.0 + (b * (0.5 + (b * 0.16666666666666666))))));
}
return tmp;
}
def code(a, b): tmp = 0 if b <= -1.55: tmp = 1.0 else: tmp = 1.0 / (2.0 + (b * (1.0 + (b * (0.5 + (b * 0.16666666666666666)))))) return tmp
function code(a, b) tmp = 0.0 if (b <= -1.55) tmp = 1.0; else tmp = Float64(1.0 / Float64(2.0 + Float64(b * Float64(1.0 + Float64(b * Float64(0.5 + Float64(b * 0.16666666666666666))))))); end return tmp end
function tmp_2 = code(a, b) tmp = 0.0; if (b <= -1.55) tmp = 1.0; else tmp = 1.0 / (2.0 + (b * (1.0 + (b * (0.5 + (b * 0.16666666666666666)))))); end tmp_2 = tmp; end
code[a_, b_] := If[LessEqual[b, -1.55], 1.0, N[(1.0 / N[(2.0 + N[(b * N[(1.0 + N[(b * N[(0.5 + N[(b * 0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq -1.55:\\
\;\;\;\;1\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{2 + b \cdot \left(1 + b \cdot \left(0.5 + b \cdot 0.16666666666666666\right)\right)}\\
\end{array}
\end{array}
if b < -1.55000000000000004Initial program 100.0%
Taylor expanded in a around 0
associate-+r+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6499.8%
Simplified99.8%
Taylor expanded in a around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64100.0%
Simplified100.0%
Taylor expanded in a around inf
Simplified100.0%
if -1.55000000000000004 < b Initial program 99.0%
Taylor expanded in a around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
exp-lowering-exp.f6473.7%
Simplified73.7%
Taylor expanded in b around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6465.4%
Simplified65.4%
(FPCore (a b) :precision binary64 (if (<= b -0.061) 1.0 (if (<= b 1.75) (+ 0.5 (* a 0.25)) (/ 2.0 (* b b)))))
double code(double a, double b) {
double tmp;
if (b <= -0.061) {
tmp = 1.0;
} else if (b <= 1.75) {
tmp = 0.5 + (a * 0.25);
} else {
tmp = 2.0 / (b * b);
}
return tmp;
}
real(8) function code(a, b)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: tmp
if (b <= (-0.061d0)) then
tmp = 1.0d0
else if (b <= 1.75d0) then
tmp = 0.5d0 + (a * 0.25d0)
else
tmp = 2.0d0 / (b * b)
end if
code = tmp
end function
public static double code(double a, double b) {
double tmp;
if (b <= -0.061) {
tmp = 1.0;
} else if (b <= 1.75) {
tmp = 0.5 + (a * 0.25);
} else {
tmp = 2.0 / (b * b);
}
return tmp;
}
def code(a, b): tmp = 0 if b <= -0.061: tmp = 1.0 elif b <= 1.75: tmp = 0.5 + (a * 0.25) else: tmp = 2.0 / (b * b) return tmp
function code(a, b) tmp = 0.0 if (b <= -0.061) tmp = 1.0; elseif (b <= 1.75) tmp = Float64(0.5 + Float64(a * 0.25)); else tmp = Float64(2.0 / Float64(b * b)); end return tmp end
function tmp_2 = code(a, b) tmp = 0.0; if (b <= -0.061) tmp = 1.0; elseif (b <= 1.75) tmp = 0.5 + (a * 0.25); else tmp = 2.0 / (b * b); end tmp_2 = tmp; end
code[a_, b_] := If[LessEqual[b, -0.061], 1.0, If[LessEqual[b, 1.75], N[(0.5 + N[(a * 0.25), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(b * b), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq -0.061:\\
\;\;\;\;1\\
\mathbf{elif}\;b \leq 1.75:\\
\;\;\;\;0.5 + a \cdot 0.25\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{b \cdot b}\\
\end{array}
\end{array}
if b < -0.060999999999999999Initial program 100.0%
Taylor expanded in a around 0
associate-+r+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6499.8%
Simplified99.8%
Taylor expanded in a around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f6498.0%
Simplified98.0%
Taylor expanded in a around inf
Simplified98.0%
if -0.060999999999999999 < b < 1.75Initial program 99.3%
Taylor expanded in b around 0
Simplified97.5%
Taylor expanded in a around 0
+-lowering-+.f64N/A
*-lowering-*.f6463.4%
Simplified63.4%
if 1.75 < b Initial program 98.3%
Taylor expanded in a around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
exp-lowering-exp.f64100.0%
Simplified100.0%
Taylor expanded in b around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f6449.6%
Simplified49.6%
Taylor expanded in b around inf
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6449.6%
Simplified49.6%
Final simplification66.9%
(FPCore (a b) :precision binary64 (if (<= b -0.95) 1.0 (/ 1.0 (+ b 2.0))))
double code(double a, double b) {
double tmp;
if (b <= -0.95) {
tmp = 1.0;
} else {
tmp = 1.0 / (b + 2.0);
}
return tmp;
}
real(8) function code(a, b)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: tmp
if (b <= (-0.95d0)) then
tmp = 1.0d0
else
tmp = 1.0d0 / (b + 2.0d0)
end if
code = tmp
end function
public static double code(double a, double b) {
double tmp;
if (b <= -0.95) {
tmp = 1.0;
} else {
tmp = 1.0 / (b + 2.0);
}
return tmp;
}
def code(a, b): tmp = 0 if b <= -0.95: tmp = 1.0 else: tmp = 1.0 / (b + 2.0) return tmp
function code(a, b) tmp = 0.0 if (b <= -0.95) tmp = 1.0; else tmp = Float64(1.0 / Float64(b + 2.0)); end return tmp end
function tmp_2 = code(a, b) tmp = 0.0; if (b <= -0.95) tmp = 1.0; else tmp = 1.0 / (b + 2.0); end tmp_2 = tmp; end
code[a_, b_] := If[LessEqual[b, -0.95], 1.0, N[(1.0 / N[(b + 2.0), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq -0.95:\\
\;\;\;\;1\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{b + 2}\\
\end{array}
\end{array}
if b < -0.94999999999999996Initial program 100.0%
Taylor expanded in a around 0
associate-+r+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6499.8%
Simplified99.8%
Taylor expanded in a around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64100.0%
Simplified100.0%
Taylor expanded in a around inf
Simplified100.0%
if -0.94999999999999996 < b Initial program 99.0%
Taylor expanded in a around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
exp-lowering-exp.f6473.7%
Simplified73.7%
Taylor expanded in b around 0
+-commutativeN/A
+-lowering-+.f6446.4%
Simplified46.4%
(FPCore (a b) :precision binary64 (if (<= b -0.075) 1.0 (+ 0.5 (* a 0.25))))
double code(double a, double b) {
double tmp;
if (b <= -0.075) {
tmp = 1.0;
} else {
tmp = 0.5 + (a * 0.25);
}
return tmp;
}
real(8) function code(a, b)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: tmp
if (b <= (-0.075d0)) then
tmp = 1.0d0
else
tmp = 0.5d0 + (a * 0.25d0)
end if
code = tmp
end function
public static double code(double a, double b) {
double tmp;
if (b <= -0.075) {
tmp = 1.0;
} else {
tmp = 0.5 + (a * 0.25);
}
return tmp;
}
def code(a, b): tmp = 0 if b <= -0.075: tmp = 1.0 else: tmp = 0.5 + (a * 0.25) return tmp
function code(a, b) tmp = 0.0 if (b <= -0.075) tmp = 1.0; else tmp = Float64(0.5 + Float64(a * 0.25)); end return tmp end
function tmp_2 = code(a, b) tmp = 0.0; if (b <= -0.075) tmp = 1.0; else tmp = 0.5 + (a * 0.25); end tmp_2 = tmp; end
code[a_, b_] := If[LessEqual[b, -0.075], 1.0, N[(0.5 + N[(a * 0.25), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq -0.075:\\
\;\;\;\;1\\
\mathbf{else}:\\
\;\;\;\;0.5 + a \cdot 0.25\\
\end{array}
\end{array}
if b < -0.0749999999999999972Initial program 100.0%
Taylor expanded in a around 0
associate-+r+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6499.8%
Simplified99.8%
Taylor expanded in a around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f6498.0%
Simplified98.0%
Taylor expanded in a around inf
Simplified98.0%
if -0.0749999999999999972 < b Initial program 99.0%
Taylor expanded in b around 0
Simplified79.4%
Taylor expanded in a around 0
+-lowering-+.f64N/A
*-lowering-*.f6446.4%
Simplified46.4%
Final simplification56.3%
(FPCore (a b) :precision binary64 (if (<= b -1.8) 1.0 (+ 0.5 (* b -0.25))))
double code(double a, double b) {
double tmp;
if (b <= -1.8) {
tmp = 1.0;
} else {
tmp = 0.5 + (b * -0.25);
}
return tmp;
}
real(8) function code(a, b)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: tmp
if (b <= (-1.8d0)) then
tmp = 1.0d0
else
tmp = 0.5d0 + (b * (-0.25d0))
end if
code = tmp
end function
public static double code(double a, double b) {
double tmp;
if (b <= -1.8) {
tmp = 1.0;
} else {
tmp = 0.5 + (b * -0.25);
}
return tmp;
}
def code(a, b): tmp = 0 if b <= -1.8: tmp = 1.0 else: tmp = 0.5 + (b * -0.25) return tmp
function code(a, b) tmp = 0.0 if (b <= -1.8) tmp = 1.0; else tmp = Float64(0.5 + Float64(b * -0.25)); end return tmp end
function tmp_2 = code(a, b) tmp = 0.0; if (b <= -1.8) tmp = 1.0; else tmp = 0.5 + (b * -0.25); end tmp_2 = tmp; end
code[a_, b_] := If[LessEqual[b, -1.8], 1.0, N[(0.5 + N[(b * -0.25), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq -1.8:\\
\;\;\;\;1\\
\mathbf{else}:\\
\;\;\;\;0.5 + b \cdot -0.25\\
\end{array}
\end{array}
if b < -1.80000000000000004Initial program 100.0%
Taylor expanded in a around 0
associate-+r+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6499.8%
Simplified99.8%
Taylor expanded in a around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64100.0%
Simplified100.0%
Taylor expanded in a around inf
Simplified100.0%
if -1.80000000000000004 < b Initial program 99.0%
Taylor expanded in a around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
exp-lowering-exp.f6473.7%
Simplified73.7%
Taylor expanded in b around 0
+-lowering-+.f64N/A
*-lowering-*.f6445.8%
Simplified45.8%
Final simplification56.0%
(FPCore (a b) :precision binary64 (if (<= b -0.94) 1.0 0.5))
double code(double a, double b) {
double tmp;
if (b <= -0.94) {
tmp = 1.0;
} else {
tmp = 0.5;
}
return tmp;
}
real(8) function code(a, b)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: tmp
if (b <= (-0.94d0)) then
tmp = 1.0d0
else
tmp = 0.5d0
end if
code = tmp
end function
public static double code(double a, double b) {
double tmp;
if (b <= -0.94) {
tmp = 1.0;
} else {
tmp = 0.5;
}
return tmp;
}
def code(a, b): tmp = 0 if b <= -0.94: tmp = 1.0 else: tmp = 0.5 return tmp
function code(a, b) tmp = 0.0 if (b <= -0.94) tmp = 1.0; else tmp = 0.5; end return tmp end
function tmp_2 = code(a, b) tmp = 0.0; if (b <= -0.94) tmp = 1.0; else tmp = 0.5; end tmp_2 = tmp; end
code[a_, b_] := If[LessEqual[b, -0.94], 1.0, 0.5]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq -0.94:\\
\;\;\;\;1\\
\mathbf{else}:\\
\;\;\;\;0.5\\
\end{array}
\end{array}
if b < -0.93999999999999995Initial program 100.0%
Taylor expanded in a around 0
associate-+r+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6499.8%
Simplified99.8%
Taylor expanded in a around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64100.0%
Simplified100.0%
Taylor expanded in a around inf
Simplified100.0%
if -0.93999999999999995 < b Initial program 99.0%
Taylor expanded in a around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
exp-lowering-exp.f6473.7%
Simplified73.7%
Taylor expanded in b around 0
Simplified45.3%
(FPCore (a b) :precision binary64 0.5)
double code(double a, double b) {
return 0.5;
}
real(8) function code(a, b)
real(8), intent (in) :: a
real(8), intent (in) :: b
code = 0.5d0
end function
public static double code(double a, double b) {
return 0.5;
}
def code(a, b): return 0.5
function code(a, b) return 0.5 end
function tmp = code(a, b) tmp = 0.5; end
code[a_, b_] := 0.5
\begin{array}{l}
\\
0.5
\end{array}
Initial program 99.2%
Taylor expanded in a around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
exp-lowering-exp.f6478.6%
Simplified78.6%
Taylor expanded in b around 0
Simplified40.3%
(FPCore (a b) :precision binary64 (/ 1.0 (+ 1.0 (exp (- b a)))))
double code(double a, double b) {
return 1.0 / (1.0 + exp((b - a)));
}
real(8) function code(a, b)
real(8), intent (in) :: a
real(8), intent (in) :: b
code = 1.0d0 / (1.0d0 + exp((b - a)))
end function
public static double code(double a, double b) {
return 1.0 / (1.0 + Math.exp((b - a)));
}
def code(a, b): return 1.0 / (1.0 + math.exp((b - a)))
function code(a, b) return Float64(1.0 / Float64(1.0 + exp(Float64(b - a)))) end
function tmp = code(a, b) tmp = 1.0 / (1.0 + exp((b - a))); end
code[a_, b_] := N[(1.0 / N[(1.0 + N[Exp[N[(b - a), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{1}{1 + e^{b - a}}
\end{array}
herbie shell --seed 2024192
(FPCore (a b)
:name "Quotient of sum of exps"
:precision binary64
:alt
(! :herbie-platform default (/ 1 (+ 1 (exp (- b a)))))
(/ (exp a) (+ (exp a) (exp b))))