
(FPCore (x y z t) :precision binary64 (- x (/ (log (+ (- 1.0 y) (* y (exp z)))) t)))
double code(double x, double y, double z, double t) {
return x - (log(((1.0 - y) + (y * exp(z)))) / t);
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
code = x - (log(((1.0d0 - y) + (y * exp(z)))) / t)
end function
public static double code(double x, double y, double z, double t) {
return x - (Math.log(((1.0 - y) + (y * Math.exp(z)))) / t);
}
def code(x, y, z, t): return x - (math.log(((1.0 - y) + (y * math.exp(z)))) / t)
function code(x, y, z, t) return Float64(x - Float64(log(Float64(Float64(1.0 - y) + Float64(y * exp(z)))) / t)) end
function tmp = code(x, y, z, t) tmp = x - (log(((1.0 - y) + (y * exp(z)))) / t); end
code[x_, y_, z_, t_] := N[(x - N[(N[Log[N[(N[(1.0 - y), $MachinePrecision] + N[(y * N[Exp[z], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x - \frac{\log \left(\left(1 - y\right) + y \cdot e^{z}\right)}{t}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 14 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x y z t) :precision binary64 (- x (/ (log (+ (- 1.0 y) (* y (exp z)))) t)))
double code(double x, double y, double z, double t) {
return x - (log(((1.0 - y) + (y * exp(z)))) / t);
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
code = x - (log(((1.0d0 - y) + (y * exp(z)))) / t)
end function
public static double code(double x, double y, double z, double t) {
return x - (Math.log(((1.0 - y) + (y * Math.exp(z)))) / t);
}
def code(x, y, z, t): return x - (math.log(((1.0 - y) + (y * math.exp(z)))) / t)
function code(x, y, z, t) return Float64(x - Float64(log(Float64(Float64(1.0 - y) + Float64(y * exp(z)))) / t)) end
function tmp = code(x, y, z, t) tmp = x - (log(((1.0 - y) + (y * exp(z)))) / t); end
code[x_, y_, z_, t_] := N[(x - N[(N[Log[N[(N[(1.0 - y), $MachinePrecision] + N[(y * N[Exp[z], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x - \frac{\log \left(\left(1 - y\right) + y \cdot e^{z}\right)}{t}
\end{array}
(FPCore (x y z t) :precision binary64 (- x (/ (log1p (* y (expm1 z))) t)))
double code(double x, double y, double z, double t) {
return x - (log1p((y * expm1(z))) / t);
}
public static double code(double x, double y, double z, double t) {
return x - (Math.log1p((y * Math.expm1(z))) / t);
}
def code(x, y, z, t): return x - (math.log1p((y * math.expm1(z))) / t)
function code(x, y, z, t) return Float64(x - Float64(log1p(Float64(y * expm1(z))) / t)) end
code[x_, y_, z_, t_] := N[(x - N[(N[Log[1 + N[(y * N[(Exp[z] - 1), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x - \frac{\mathsf{log1p}\left(y \cdot \mathsf{expm1}\left(z\right)\right)}{t}
\end{array}
Initial program 60.4%
--lowering--.f64N/A
/-lowering-/.f64N/A
sub-negN/A
associate-+l+N/A
log1p-defineN/A
log1p-lowering-log1p.f64N/A
neg-mul-1N/A
*-commutativeN/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
+-commutativeN/A
metadata-evalN/A
sub-negN/A
expm1-defineN/A
expm1-lowering-expm1.f6498.8%
Simplified98.8%
(FPCore (x y z t)
:precision binary64
(if (<= z -0.027)
(+ x (/ -1.0 (/ (+ (* 0.5 (* y t)) (/ t (expm1 z))) y)))
(+
x
(/
-1.0
(/
t
(log1p
(*
y
(*
z
(+
1.0
(*
z
(+
0.5
(* z (+ 0.16666666666666666 (* z 0.041666666666666664))))))))))))))
double code(double x, double y, double z, double t) {
double tmp;
if (z <= -0.027) {
tmp = x + (-1.0 / (((0.5 * (y * t)) + (t / expm1(z))) / y));
} else {
tmp = x + (-1.0 / (t / log1p((y * (z * (1.0 + (z * (0.5 + (z * (0.16666666666666666 + (z * 0.041666666666666664)))))))))));
}
return tmp;
}
public static double code(double x, double y, double z, double t) {
double tmp;
if (z <= -0.027) {
tmp = x + (-1.0 / (((0.5 * (y * t)) + (t / Math.expm1(z))) / y));
} else {
tmp = x + (-1.0 / (t / Math.log1p((y * (z * (1.0 + (z * (0.5 + (z * (0.16666666666666666 + (z * 0.041666666666666664)))))))))));
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if z <= -0.027: tmp = x + (-1.0 / (((0.5 * (y * t)) + (t / math.expm1(z))) / y)) else: tmp = x + (-1.0 / (t / math.log1p((y * (z * (1.0 + (z * (0.5 + (z * (0.16666666666666666 + (z * 0.041666666666666664))))))))))) return tmp
function code(x, y, z, t) tmp = 0.0 if (z <= -0.027) tmp = Float64(x + Float64(-1.0 / Float64(Float64(Float64(0.5 * Float64(y * t)) + Float64(t / expm1(z))) / y))); else tmp = Float64(x + Float64(-1.0 / Float64(t / log1p(Float64(y * Float64(z * Float64(1.0 + Float64(z * Float64(0.5 + Float64(z * Float64(0.16666666666666666 + Float64(z * 0.041666666666666664)))))))))))); end return tmp end
code[x_, y_, z_, t_] := If[LessEqual[z, -0.027], N[(x + N[(-1.0 / N[(N[(N[(0.5 * N[(y * t), $MachinePrecision]), $MachinePrecision] + N[(t / N[(Exp[z] - 1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x + N[(-1.0 / N[(t / N[Log[1 + N[(y * N[(z * N[(1.0 + N[(z * N[(0.5 + N[(z * N[(0.16666666666666666 + N[(z * 0.041666666666666664), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -0.027:\\
\;\;\;\;x + \frac{-1}{\frac{0.5 \cdot \left(y \cdot t\right) + \frac{t}{\mathsf{expm1}\left(z\right)}}{y}}\\
\mathbf{else}:\\
\;\;\;\;x + \frac{-1}{\frac{t}{\mathsf{log1p}\left(y \cdot \left(z \cdot \left(1 + z \cdot \left(0.5 + z \cdot \left(0.16666666666666666 + z \cdot 0.041666666666666664\right)\right)\right)\right)\right)}}\\
\end{array}
\end{array}
if z < -0.0269999999999999997Initial program 82.5%
--lowering--.f64N/A
/-lowering-/.f64N/A
sub-negN/A
associate-+l+N/A
log1p-defineN/A
log1p-lowering-log1p.f64N/A
neg-mul-1N/A
*-commutativeN/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
+-commutativeN/A
metadata-evalN/A
sub-negN/A
expm1-defineN/A
expm1-lowering-expm1.f6499.9%
Simplified99.9%
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
log1p-defineN/A
log1p-lowering-log1p.f64N/A
*-lowering-*.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f6499.9%
Applied egg-rr99.9%
Taylor expanded in y around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f6483.7%
Simplified83.7%
if -0.0269999999999999997 < z Initial program 49.5%
--lowering--.f64N/A
/-lowering-/.f64N/A
sub-negN/A
associate-+l+N/A
log1p-defineN/A
log1p-lowering-log1p.f64N/A
neg-mul-1N/A
*-commutativeN/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
+-commutativeN/A
metadata-evalN/A
sub-negN/A
expm1-defineN/A
expm1-lowering-expm1.f6498.2%
Simplified98.2%
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
log1p-defineN/A
log1p-lowering-log1p.f64N/A
*-lowering-*.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f6498.2%
Applied egg-rr98.2%
Taylor expanded in z around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6498.3%
Simplified98.3%
Final simplification93.4%
(FPCore (x y z t)
:precision binary64
(if (<= z -0.0075)
(+ x (/ -1.0 (/ (+ (* 0.5 (* y t)) (/ t (expm1 z))) y)))
(-
x
(/
(log1p (* z (+ y (* z (+ (* y (* z 0.16666666666666666)) (* y 0.5))))))
t))))
double code(double x, double y, double z, double t) {
double tmp;
if (z <= -0.0075) {
tmp = x + (-1.0 / (((0.5 * (y * t)) + (t / expm1(z))) / y));
} else {
tmp = x - (log1p((z * (y + (z * ((y * (z * 0.16666666666666666)) + (y * 0.5)))))) / t);
}
return tmp;
}
public static double code(double x, double y, double z, double t) {
double tmp;
if (z <= -0.0075) {
tmp = x + (-1.0 / (((0.5 * (y * t)) + (t / Math.expm1(z))) / y));
} else {
tmp = x - (Math.log1p((z * (y + (z * ((y * (z * 0.16666666666666666)) + (y * 0.5)))))) / t);
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if z <= -0.0075: tmp = x + (-1.0 / (((0.5 * (y * t)) + (t / math.expm1(z))) / y)) else: tmp = x - (math.log1p((z * (y + (z * ((y * (z * 0.16666666666666666)) + (y * 0.5)))))) / t) return tmp
function code(x, y, z, t) tmp = 0.0 if (z <= -0.0075) tmp = Float64(x + Float64(-1.0 / Float64(Float64(Float64(0.5 * Float64(y * t)) + Float64(t / expm1(z))) / y))); else tmp = Float64(x - Float64(log1p(Float64(z * Float64(y + Float64(z * Float64(Float64(y * Float64(z * 0.16666666666666666)) + Float64(y * 0.5)))))) / t)); end return tmp end
code[x_, y_, z_, t_] := If[LessEqual[z, -0.0075], N[(x + N[(-1.0 / N[(N[(N[(0.5 * N[(y * t), $MachinePrecision]), $MachinePrecision] + N[(t / N[(Exp[z] - 1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x - N[(N[Log[1 + N[(z * N[(y + N[(z * N[(N[(y * N[(z * 0.16666666666666666), $MachinePrecision]), $MachinePrecision] + N[(y * 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -0.0075:\\
\;\;\;\;x + \frac{-1}{\frac{0.5 \cdot \left(y \cdot t\right) + \frac{t}{\mathsf{expm1}\left(z\right)}}{y}}\\
\mathbf{else}:\\
\;\;\;\;x - \frac{\mathsf{log1p}\left(z \cdot \left(y + z \cdot \left(y \cdot \left(z \cdot 0.16666666666666666\right) + y \cdot 0.5\right)\right)\right)}{t}\\
\end{array}
\end{array}
if z < -0.0074999999999999997Initial program 82.5%
--lowering--.f64N/A
/-lowering-/.f64N/A
sub-negN/A
associate-+l+N/A
log1p-defineN/A
log1p-lowering-log1p.f64N/A
neg-mul-1N/A
*-commutativeN/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
+-commutativeN/A
metadata-evalN/A
sub-negN/A
expm1-defineN/A
expm1-lowering-expm1.f6499.9%
Simplified99.9%
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
log1p-defineN/A
log1p-lowering-log1p.f64N/A
*-lowering-*.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f6499.9%
Applied egg-rr99.9%
Taylor expanded in y around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f6483.7%
Simplified83.7%
if -0.0074999999999999997 < z Initial program 49.5%
--lowering--.f64N/A
/-lowering-/.f64N/A
sub-negN/A
associate-+l+N/A
log1p-defineN/A
log1p-lowering-log1p.f64N/A
neg-mul-1N/A
*-commutativeN/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
+-commutativeN/A
metadata-evalN/A
sub-negN/A
expm1-defineN/A
expm1-lowering-expm1.f6498.2%
Simplified98.2%
Taylor expanded in z around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6498.2%
Simplified98.2%
Final simplification93.4%
(FPCore (x y z t)
:precision binary64
(if (<= z -1.5e+158)
(- x (/ (* y (expm1 z)) t))
(if (<= z -0.03)
(+ x (/ 1.0 (/ (- (* z (* 0.5 (- (/ t y) t))) (/ t y)) z)))
(- x (/ (log1p (* z (+ y (* y (* z 0.5))))) t)))))
double code(double x, double y, double z, double t) {
double tmp;
if (z <= -1.5e+158) {
tmp = x - ((y * expm1(z)) / t);
} else if (z <= -0.03) {
tmp = x + (1.0 / (((z * (0.5 * ((t / y) - t))) - (t / y)) / z));
} else {
tmp = x - (log1p((z * (y + (y * (z * 0.5))))) / t);
}
return tmp;
}
public static double code(double x, double y, double z, double t) {
double tmp;
if (z <= -1.5e+158) {
tmp = x - ((y * Math.expm1(z)) / t);
} else if (z <= -0.03) {
tmp = x + (1.0 / (((z * (0.5 * ((t / y) - t))) - (t / y)) / z));
} else {
tmp = x - (Math.log1p((z * (y + (y * (z * 0.5))))) / t);
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if z <= -1.5e+158: tmp = x - ((y * math.expm1(z)) / t) elif z <= -0.03: tmp = x + (1.0 / (((z * (0.5 * ((t / y) - t))) - (t / y)) / z)) else: tmp = x - (math.log1p((z * (y + (y * (z * 0.5))))) / t) return tmp
function code(x, y, z, t) tmp = 0.0 if (z <= -1.5e+158) tmp = Float64(x - Float64(Float64(y * expm1(z)) / t)); elseif (z <= -0.03) tmp = Float64(x + Float64(1.0 / Float64(Float64(Float64(z * Float64(0.5 * Float64(Float64(t / y) - t))) - Float64(t / y)) / z))); else tmp = Float64(x - Float64(log1p(Float64(z * Float64(y + Float64(y * Float64(z * 0.5))))) / t)); end return tmp end
code[x_, y_, z_, t_] := If[LessEqual[z, -1.5e+158], N[(x - N[(N[(y * N[(Exp[z] - 1), $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, -0.03], N[(x + N[(1.0 / N[(N[(N[(z * N[(0.5 * N[(N[(t / y), $MachinePrecision] - t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(t / y), $MachinePrecision]), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x - N[(N[Log[1 + N[(z * N[(y + N[(y * N[(z * 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -1.5 \cdot 10^{+158}:\\
\;\;\;\;x - \frac{y \cdot \mathsf{expm1}\left(z\right)}{t}\\
\mathbf{elif}\;z \leq -0.03:\\
\;\;\;\;x + \frac{1}{\frac{z \cdot \left(0.5 \cdot \left(\frac{t}{y} - t\right)\right) - \frac{t}{y}}{z}}\\
\mathbf{else}:\\
\;\;\;\;x - \frac{\mathsf{log1p}\left(z \cdot \left(y + y \cdot \left(z \cdot 0.5\right)\right)\right)}{t}\\
\end{array}
\end{array}
if z < -1.5e158Initial program 70.7%
--lowering--.f64N/A
/-lowering-/.f64N/A
sub-negN/A
associate-+l+N/A
log1p-defineN/A
log1p-lowering-log1p.f64N/A
neg-mul-1N/A
*-commutativeN/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
+-commutativeN/A
metadata-evalN/A
sub-negN/A
expm1-defineN/A
expm1-lowering-expm1.f6499.9%
Simplified99.9%
Taylor expanded in y around 0
*-lowering-*.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f6482.7%
Simplified82.7%
if -1.5e158 < z < -0.029999999999999999Initial program 92.4%
--lowering--.f64N/A
/-lowering-/.f64N/A
sub-negN/A
associate-+l+N/A
log1p-defineN/A
log1p-lowering-log1p.f64N/A
neg-mul-1N/A
*-commutativeN/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
+-commutativeN/A
metadata-evalN/A
sub-negN/A
expm1-defineN/A
expm1-lowering-expm1.f6499.9%
Simplified99.9%
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
log1p-defineN/A
log1p-lowering-log1p.f64N/A
*-lowering-*.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f6499.9%
Applied egg-rr99.9%
Taylor expanded in y around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f6481.2%
Simplified81.2%
Taylor expanded in z around 0
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
distribute-lft-out--N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
/-lowering-/.f6475.4%
Simplified75.4%
if -0.029999999999999999 < z Initial program 49.5%
--lowering--.f64N/A
/-lowering-/.f64N/A
sub-negN/A
associate-+l+N/A
log1p-defineN/A
log1p-lowering-log1p.f64N/A
neg-mul-1N/A
*-commutativeN/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
+-commutativeN/A
metadata-evalN/A
sub-negN/A
expm1-defineN/A
expm1-lowering-expm1.f6498.2%
Simplified98.2%
Taylor expanded in z around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6498.1%
Simplified98.1%
Final simplification91.7%
(FPCore (x y z t) :precision binary64 (if (<= z -0.005) (+ x (/ -1.0 (/ (+ (* 0.5 (* y t)) (/ t (expm1 z))) y))) (- x (/ (log1p (* z (+ y (* y (* z 0.5))))) t))))
double code(double x, double y, double z, double t) {
double tmp;
if (z <= -0.005) {
tmp = x + (-1.0 / (((0.5 * (y * t)) + (t / expm1(z))) / y));
} else {
tmp = x - (log1p((z * (y + (y * (z * 0.5))))) / t);
}
return tmp;
}
public static double code(double x, double y, double z, double t) {
double tmp;
if (z <= -0.005) {
tmp = x + (-1.0 / (((0.5 * (y * t)) + (t / Math.expm1(z))) / y));
} else {
tmp = x - (Math.log1p((z * (y + (y * (z * 0.5))))) / t);
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if z <= -0.005: tmp = x + (-1.0 / (((0.5 * (y * t)) + (t / math.expm1(z))) / y)) else: tmp = x - (math.log1p((z * (y + (y * (z * 0.5))))) / t) return tmp
function code(x, y, z, t) tmp = 0.0 if (z <= -0.005) tmp = Float64(x + Float64(-1.0 / Float64(Float64(Float64(0.5 * Float64(y * t)) + Float64(t / expm1(z))) / y))); else tmp = Float64(x - Float64(log1p(Float64(z * Float64(y + Float64(y * Float64(z * 0.5))))) / t)); end return tmp end
code[x_, y_, z_, t_] := If[LessEqual[z, -0.005], N[(x + N[(-1.0 / N[(N[(N[(0.5 * N[(y * t), $MachinePrecision]), $MachinePrecision] + N[(t / N[(Exp[z] - 1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x - N[(N[Log[1 + N[(z * N[(y + N[(y * N[(z * 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -0.005:\\
\;\;\;\;x + \frac{-1}{\frac{0.5 \cdot \left(y \cdot t\right) + \frac{t}{\mathsf{expm1}\left(z\right)}}{y}}\\
\mathbf{else}:\\
\;\;\;\;x - \frac{\mathsf{log1p}\left(z \cdot \left(y + y \cdot \left(z \cdot 0.5\right)\right)\right)}{t}\\
\end{array}
\end{array}
if z < -0.0050000000000000001Initial program 82.5%
--lowering--.f64N/A
/-lowering-/.f64N/A
sub-negN/A
associate-+l+N/A
log1p-defineN/A
log1p-lowering-log1p.f64N/A
neg-mul-1N/A
*-commutativeN/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
+-commutativeN/A
metadata-evalN/A
sub-negN/A
expm1-defineN/A
expm1-lowering-expm1.f6499.9%
Simplified99.9%
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
log1p-defineN/A
log1p-lowering-log1p.f64N/A
*-lowering-*.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f6499.9%
Applied egg-rr99.9%
Taylor expanded in y around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f6483.7%
Simplified83.7%
if -0.0050000000000000001 < z Initial program 49.5%
--lowering--.f64N/A
/-lowering-/.f64N/A
sub-negN/A
associate-+l+N/A
log1p-defineN/A
log1p-lowering-log1p.f64N/A
neg-mul-1N/A
*-commutativeN/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
+-commutativeN/A
metadata-evalN/A
sub-negN/A
expm1-defineN/A
expm1-lowering-expm1.f6498.2%
Simplified98.2%
Taylor expanded in z around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6498.1%
Simplified98.1%
Final simplification93.3%
(FPCore (x y z t)
:precision binary64
(if (<= z -2.45e+156)
(- x (/ (* y (expm1 z)) t))
(if (<= z -4e+54)
(+ x (/ 1.0 (/ (- (* z (* 0.5 (- (/ t y) t))) (/ t y)) z)))
(+ x (/ -1.0 (/ t (log1p (* y z))))))))
double code(double x, double y, double z, double t) {
double tmp;
if (z <= -2.45e+156) {
tmp = x - ((y * expm1(z)) / t);
} else if (z <= -4e+54) {
tmp = x + (1.0 / (((z * (0.5 * ((t / y) - t))) - (t / y)) / z));
} else {
tmp = x + (-1.0 / (t / log1p((y * z))));
}
return tmp;
}
public static double code(double x, double y, double z, double t) {
double tmp;
if (z <= -2.45e+156) {
tmp = x - ((y * Math.expm1(z)) / t);
} else if (z <= -4e+54) {
tmp = x + (1.0 / (((z * (0.5 * ((t / y) - t))) - (t / y)) / z));
} else {
tmp = x + (-1.0 / (t / Math.log1p((y * z))));
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if z <= -2.45e+156: tmp = x - ((y * math.expm1(z)) / t) elif z <= -4e+54: tmp = x + (1.0 / (((z * (0.5 * ((t / y) - t))) - (t / y)) / z)) else: tmp = x + (-1.0 / (t / math.log1p((y * z)))) return tmp
function code(x, y, z, t) tmp = 0.0 if (z <= -2.45e+156) tmp = Float64(x - Float64(Float64(y * expm1(z)) / t)); elseif (z <= -4e+54) tmp = Float64(x + Float64(1.0 / Float64(Float64(Float64(z * Float64(0.5 * Float64(Float64(t / y) - t))) - Float64(t / y)) / z))); else tmp = Float64(x + Float64(-1.0 / Float64(t / log1p(Float64(y * z))))); end return tmp end
code[x_, y_, z_, t_] := If[LessEqual[z, -2.45e+156], N[(x - N[(N[(y * N[(Exp[z] - 1), $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, -4e+54], N[(x + N[(1.0 / N[(N[(N[(z * N[(0.5 * N[(N[(t / y), $MachinePrecision] - t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(t / y), $MachinePrecision]), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x + N[(-1.0 / N[(t / N[Log[1 + N[(y * z), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -2.45 \cdot 10^{+156}:\\
\;\;\;\;x - \frac{y \cdot \mathsf{expm1}\left(z\right)}{t}\\
\mathbf{elif}\;z \leq -4 \cdot 10^{+54}:\\
\;\;\;\;x + \frac{1}{\frac{z \cdot \left(0.5 \cdot \left(\frac{t}{y} - t\right)\right) - \frac{t}{y}}{z}}\\
\mathbf{else}:\\
\;\;\;\;x + \frac{-1}{\frac{t}{\mathsf{log1p}\left(y \cdot z\right)}}\\
\end{array}
\end{array}
if z < -2.44999999999999984e156Initial program 70.7%
--lowering--.f64N/A
/-lowering-/.f64N/A
sub-negN/A
associate-+l+N/A
log1p-defineN/A
log1p-lowering-log1p.f64N/A
neg-mul-1N/A
*-commutativeN/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
+-commutativeN/A
metadata-evalN/A
sub-negN/A
expm1-defineN/A
expm1-lowering-expm1.f6499.9%
Simplified99.9%
Taylor expanded in y around 0
*-lowering-*.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f6482.7%
Simplified82.7%
if -2.44999999999999984e156 < z < -4.0000000000000003e54Initial program 95.6%
--lowering--.f64N/A
/-lowering-/.f64N/A
sub-negN/A
associate-+l+N/A
log1p-defineN/A
log1p-lowering-log1p.f64N/A
neg-mul-1N/A
*-commutativeN/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
+-commutativeN/A
metadata-evalN/A
sub-negN/A
expm1-defineN/A
expm1-lowering-expm1.f64100.0%
Simplified100.0%
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
log1p-defineN/A
log1p-lowering-log1p.f64N/A
*-lowering-*.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f64100.0%
Applied egg-rr100.0%
Taylor expanded in y around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f6487.1%
Simplified87.1%
Taylor expanded in z around 0
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
distribute-lft-out--N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
/-lowering-/.f6483.5%
Simplified83.5%
if -4.0000000000000003e54 < z Initial program 51.7%
--lowering--.f64N/A
/-lowering-/.f64N/A
sub-negN/A
associate-+l+N/A
log1p-defineN/A
log1p-lowering-log1p.f64N/A
neg-mul-1N/A
*-commutativeN/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
+-commutativeN/A
metadata-evalN/A
sub-negN/A
expm1-defineN/A
expm1-lowering-expm1.f6498.3%
Simplified98.3%
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
log1p-defineN/A
log1p-lowering-log1p.f64N/A
*-lowering-*.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f6498.3%
Applied egg-rr98.3%
Taylor expanded in z around 0
*-lowering-*.f6494.8%
Simplified94.8%
Final simplification91.5%
(FPCore (x y z t)
:precision binary64
(if (<= z -2.35e+156)
(- x (/ (* y (expm1 z)) t))
(if (<= z -4e+54)
(+ x (/ 1.0 (/ (- (* z (* 0.5 (- (/ t y) t))) (/ t y)) z)))
(- x (/ (log1p (* y z)) t)))))
double code(double x, double y, double z, double t) {
double tmp;
if (z <= -2.35e+156) {
tmp = x - ((y * expm1(z)) / t);
} else if (z <= -4e+54) {
tmp = x + (1.0 / (((z * (0.5 * ((t / y) - t))) - (t / y)) / z));
} else {
tmp = x - (log1p((y * z)) / t);
}
return tmp;
}
public static double code(double x, double y, double z, double t) {
double tmp;
if (z <= -2.35e+156) {
tmp = x - ((y * Math.expm1(z)) / t);
} else if (z <= -4e+54) {
tmp = x + (1.0 / (((z * (0.5 * ((t / y) - t))) - (t / y)) / z));
} else {
tmp = x - (Math.log1p((y * z)) / t);
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if z <= -2.35e+156: tmp = x - ((y * math.expm1(z)) / t) elif z <= -4e+54: tmp = x + (1.0 / (((z * (0.5 * ((t / y) - t))) - (t / y)) / z)) else: tmp = x - (math.log1p((y * z)) / t) return tmp
function code(x, y, z, t) tmp = 0.0 if (z <= -2.35e+156) tmp = Float64(x - Float64(Float64(y * expm1(z)) / t)); elseif (z <= -4e+54) tmp = Float64(x + Float64(1.0 / Float64(Float64(Float64(z * Float64(0.5 * Float64(Float64(t / y) - t))) - Float64(t / y)) / z))); else tmp = Float64(x - Float64(log1p(Float64(y * z)) / t)); end return tmp end
code[x_, y_, z_, t_] := If[LessEqual[z, -2.35e+156], N[(x - N[(N[(y * N[(Exp[z] - 1), $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, -4e+54], N[(x + N[(1.0 / N[(N[(N[(z * N[(0.5 * N[(N[(t / y), $MachinePrecision] - t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(t / y), $MachinePrecision]), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x - N[(N[Log[1 + N[(y * z), $MachinePrecision]], $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -2.35 \cdot 10^{+156}:\\
\;\;\;\;x - \frac{y \cdot \mathsf{expm1}\left(z\right)}{t}\\
\mathbf{elif}\;z \leq -4 \cdot 10^{+54}:\\
\;\;\;\;x + \frac{1}{\frac{z \cdot \left(0.5 \cdot \left(\frac{t}{y} - t\right)\right) - \frac{t}{y}}{z}}\\
\mathbf{else}:\\
\;\;\;\;x - \frac{\mathsf{log1p}\left(y \cdot z\right)}{t}\\
\end{array}
\end{array}
if z < -2.35e156Initial program 70.7%
--lowering--.f64N/A
/-lowering-/.f64N/A
sub-negN/A
associate-+l+N/A
log1p-defineN/A
log1p-lowering-log1p.f64N/A
neg-mul-1N/A
*-commutativeN/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
+-commutativeN/A
metadata-evalN/A
sub-negN/A
expm1-defineN/A
expm1-lowering-expm1.f6499.9%
Simplified99.9%
Taylor expanded in y around 0
*-lowering-*.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f6482.7%
Simplified82.7%
if -2.35e156 < z < -4.0000000000000003e54Initial program 95.6%
--lowering--.f64N/A
/-lowering-/.f64N/A
sub-negN/A
associate-+l+N/A
log1p-defineN/A
log1p-lowering-log1p.f64N/A
neg-mul-1N/A
*-commutativeN/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
+-commutativeN/A
metadata-evalN/A
sub-negN/A
expm1-defineN/A
expm1-lowering-expm1.f64100.0%
Simplified100.0%
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
log1p-defineN/A
log1p-lowering-log1p.f64N/A
*-lowering-*.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f64100.0%
Applied egg-rr100.0%
Taylor expanded in y around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f6487.1%
Simplified87.1%
Taylor expanded in z around 0
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
distribute-lft-out--N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
/-lowering-/.f6483.5%
Simplified83.5%
if -4.0000000000000003e54 < z Initial program 51.7%
--lowering--.f64N/A
/-lowering-/.f64N/A
sub-negN/A
associate-+l+N/A
log1p-defineN/A
log1p-lowering-log1p.f64N/A
neg-mul-1N/A
*-commutativeN/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
+-commutativeN/A
metadata-evalN/A
sub-negN/A
expm1-defineN/A
expm1-lowering-expm1.f6498.3%
Simplified98.3%
Taylor expanded in z around 0
*-lowering-*.f6494.8%
Simplified94.8%
Final simplification91.5%
(FPCore (x y z t)
:precision binary64
(if (<= z -4.5e+156)
(- x (* y (/ (expm1 z) t)))
(if (<= z -4e+54)
(+ x (/ 1.0 (/ (- (* z (* 0.5 (- (/ t y) t))) (/ t y)) z)))
(- x (/ (log1p (* y z)) t)))))
double code(double x, double y, double z, double t) {
double tmp;
if (z <= -4.5e+156) {
tmp = x - (y * (expm1(z) / t));
} else if (z <= -4e+54) {
tmp = x + (1.0 / (((z * (0.5 * ((t / y) - t))) - (t / y)) / z));
} else {
tmp = x - (log1p((y * z)) / t);
}
return tmp;
}
public static double code(double x, double y, double z, double t) {
double tmp;
if (z <= -4.5e+156) {
tmp = x - (y * (Math.expm1(z) / t));
} else if (z <= -4e+54) {
tmp = x + (1.0 / (((z * (0.5 * ((t / y) - t))) - (t / y)) / z));
} else {
tmp = x - (Math.log1p((y * z)) / t);
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if z <= -4.5e+156: tmp = x - (y * (math.expm1(z) / t)) elif z <= -4e+54: tmp = x + (1.0 / (((z * (0.5 * ((t / y) - t))) - (t / y)) / z)) else: tmp = x - (math.log1p((y * z)) / t) return tmp
function code(x, y, z, t) tmp = 0.0 if (z <= -4.5e+156) tmp = Float64(x - Float64(y * Float64(expm1(z) / t))); elseif (z <= -4e+54) tmp = Float64(x + Float64(1.0 / Float64(Float64(Float64(z * Float64(0.5 * Float64(Float64(t / y) - t))) - Float64(t / y)) / z))); else tmp = Float64(x - Float64(log1p(Float64(y * z)) / t)); end return tmp end
code[x_, y_, z_, t_] := If[LessEqual[z, -4.5e+156], N[(x - N[(y * N[(N[(Exp[z] - 1), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, -4e+54], N[(x + N[(1.0 / N[(N[(N[(z * N[(0.5 * N[(N[(t / y), $MachinePrecision] - t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(t / y), $MachinePrecision]), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x - N[(N[Log[1 + N[(y * z), $MachinePrecision]], $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -4.5 \cdot 10^{+156}:\\
\;\;\;\;x - y \cdot \frac{\mathsf{expm1}\left(z\right)}{t}\\
\mathbf{elif}\;z \leq -4 \cdot 10^{+54}:\\
\;\;\;\;x + \frac{1}{\frac{z \cdot \left(0.5 \cdot \left(\frac{t}{y} - t\right)\right) - \frac{t}{y}}{z}}\\
\mathbf{else}:\\
\;\;\;\;x - \frac{\mathsf{log1p}\left(y \cdot z\right)}{t}\\
\end{array}
\end{array}
if z < -4.50000000000000031e156Initial program 70.7%
--lowering--.f64N/A
/-lowering-/.f64N/A
sub-negN/A
associate-+l+N/A
log1p-defineN/A
log1p-lowering-log1p.f64N/A
neg-mul-1N/A
*-commutativeN/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
+-commutativeN/A
metadata-evalN/A
sub-negN/A
expm1-defineN/A
expm1-lowering-expm1.f6499.9%
Simplified99.9%
Taylor expanded in y around 0
associate-/l*N/A
div-subN/A
*-lowering-*.f64N/A
div-subN/A
/-lowering-/.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f6482.6%
Simplified82.6%
if -4.50000000000000031e156 < z < -4.0000000000000003e54Initial program 95.6%
--lowering--.f64N/A
/-lowering-/.f64N/A
sub-negN/A
associate-+l+N/A
log1p-defineN/A
log1p-lowering-log1p.f64N/A
neg-mul-1N/A
*-commutativeN/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
+-commutativeN/A
metadata-evalN/A
sub-negN/A
expm1-defineN/A
expm1-lowering-expm1.f64100.0%
Simplified100.0%
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
log1p-defineN/A
log1p-lowering-log1p.f64N/A
*-lowering-*.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f64100.0%
Applied egg-rr100.0%
Taylor expanded in y around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f6487.1%
Simplified87.1%
Taylor expanded in z around 0
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
distribute-lft-out--N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
/-lowering-/.f6483.5%
Simplified83.5%
if -4.0000000000000003e54 < z Initial program 51.7%
--lowering--.f64N/A
/-lowering-/.f64N/A
sub-negN/A
associate-+l+N/A
log1p-defineN/A
log1p-lowering-log1p.f64N/A
neg-mul-1N/A
*-commutativeN/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
+-commutativeN/A
metadata-evalN/A
sub-negN/A
expm1-defineN/A
expm1-lowering-expm1.f6498.3%
Simplified98.3%
Taylor expanded in z around 0
*-lowering-*.f6494.8%
Simplified94.8%
Final simplification91.4%
(FPCore (x y z t) :precision binary64 (if (<= y -950000.0) (+ x (/ 1.0 (- (/ (- (* t 0.5) (/ t z)) y) (* t 0.5)))) (- x (* y (/ (expm1 z) t)))))
double code(double x, double y, double z, double t) {
double tmp;
if (y <= -950000.0) {
tmp = x + (1.0 / ((((t * 0.5) - (t / z)) / y) - (t * 0.5)));
} else {
tmp = x - (y * (expm1(z) / t));
}
return tmp;
}
public static double code(double x, double y, double z, double t) {
double tmp;
if (y <= -950000.0) {
tmp = x + (1.0 / ((((t * 0.5) - (t / z)) / y) - (t * 0.5)));
} else {
tmp = x - (y * (Math.expm1(z) / t));
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if y <= -950000.0: tmp = x + (1.0 / ((((t * 0.5) - (t / z)) / y) - (t * 0.5))) else: tmp = x - (y * (math.expm1(z) / t)) return tmp
function code(x, y, z, t) tmp = 0.0 if (y <= -950000.0) tmp = Float64(x + Float64(1.0 / Float64(Float64(Float64(Float64(t * 0.5) - Float64(t / z)) / y) - Float64(t * 0.5)))); else tmp = Float64(x - Float64(y * Float64(expm1(z) / t))); end return tmp end
code[x_, y_, z_, t_] := If[LessEqual[y, -950000.0], N[(x + N[(1.0 / N[(N[(N[(N[(t * 0.5), $MachinePrecision] - N[(t / z), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision] - N[(t * 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x - N[(y * N[(N[(Exp[z] - 1), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -950000:\\
\;\;\;\;x + \frac{1}{\frac{t \cdot 0.5 - \frac{t}{z}}{y} - t \cdot 0.5}\\
\mathbf{else}:\\
\;\;\;\;x - y \cdot \frac{\mathsf{expm1}\left(z\right)}{t}\\
\end{array}
\end{array}
if y < -9.5e5Initial program 50.4%
--lowering--.f64N/A
/-lowering-/.f64N/A
sub-negN/A
associate-+l+N/A
log1p-defineN/A
log1p-lowering-log1p.f64N/A
neg-mul-1N/A
*-commutativeN/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
+-commutativeN/A
metadata-evalN/A
sub-negN/A
expm1-defineN/A
expm1-lowering-expm1.f6499.8%
Simplified99.8%
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
log1p-defineN/A
log1p-lowering-log1p.f64N/A
*-lowering-*.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f6499.8%
Applied egg-rr99.8%
clear-numN/A
associate-/r/N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
log1p-defineN/A
log1p-lowering-log1p.f64N/A
*-lowering-*.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f6499.7%
Applied egg-rr99.7%
Taylor expanded in z around 0
/-lowering-/.f64N/A
Simplified58.9%
Taylor expanded in y around -inf
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f6467.5%
Simplified67.5%
if -9.5e5 < y Initial program 64.3%
--lowering--.f64N/A
/-lowering-/.f64N/A
sub-negN/A
associate-+l+N/A
log1p-defineN/A
log1p-lowering-log1p.f64N/A
neg-mul-1N/A
*-commutativeN/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
+-commutativeN/A
metadata-evalN/A
sub-negN/A
expm1-defineN/A
expm1-lowering-expm1.f6498.4%
Simplified98.4%
Taylor expanded in y around 0
associate-/l*N/A
div-subN/A
*-lowering-*.f64N/A
div-subN/A
/-lowering-/.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f6494.0%
Simplified94.0%
Final simplification86.6%
(FPCore (x y z t) :precision binary64 (if (<= z -650000000.0) (+ x (* (/ 2.0 t) (/ y (- 1.0 y)))) (+ x (* y (/ (* z (- -1.0 (* z (+ 0.5 (* z 0.16666666666666666))))) t)))))
double code(double x, double y, double z, double t) {
double tmp;
if (z <= -650000000.0) {
tmp = x + ((2.0 / t) * (y / (1.0 - y)));
} else {
tmp = x + (y * ((z * (-1.0 - (z * (0.5 + (z * 0.16666666666666666))))) / t));
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: tmp
if (z <= (-650000000.0d0)) then
tmp = x + ((2.0d0 / t) * (y / (1.0d0 - y)))
else
tmp = x + (y * ((z * ((-1.0d0) - (z * (0.5d0 + (z * 0.16666666666666666d0))))) / t))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double tmp;
if (z <= -650000000.0) {
tmp = x + ((2.0 / t) * (y / (1.0 - y)));
} else {
tmp = x + (y * ((z * (-1.0 - (z * (0.5 + (z * 0.16666666666666666))))) / t));
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if z <= -650000000.0: tmp = x + ((2.0 / t) * (y / (1.0 - y))) else: tmp = x + (y * ((z * (-1.0 - (z * (0.5 + (z * 0.16666666666666666))))) / t)) return tmp
function code(x, y, z, t) tmp = 0.0 if (z <= -650000000.0) tmp = Float64(x + Float64(Float64(2.0 / t) * Float64(y / Float64(1.0 - y)))); else tmp = Float64(x + Float64(y * Float64(Float64(z * Float64(-1.0 - Float64(z * Float64(0.5 + Float64(z * 0.16666666666666666))))) / t))); end return tmp end
function tmp_2 = code(x, y, z, t) tmp = 0.0; if (z <= -650000000.0) tmp = x + ((2.0 / t) * (y / (1.0 - y))); else tmp = x + (y * ((z * (-1.0 - (z * (0.5 + (z * 0.16666666666666666))))) / t)); end tmp_2 = tmp; end
code[x_, y_, z_, t_] := If[LessEqual[z, -650000000.0], N[(x + N[(N[(2.0 / t), $MachinePrecision] * N[(y / N[(1.0 - y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x + N[(y * N[(N[(z * N[(-1.0 - N[(z * N[(0.5 + N[(z * 0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -650000000:\\
\;\;\;\;x + \frac{2}{t} \cdot \frac{y}{1 - y}\\
\mathbf{else}:\\
\;\;\;\;x + y \cdot \frac{z \cdot \left(-1 - z \cdot \left(0.5 + z \cdot 0.16666666666666666\right)\right)}{t}\\
\end{array}
\end{array}
if z < -6.5e8Initial program 83.2%
--lowering--.f64N/A
/-lowering-/.f64N/A
sub-negN/A
associate-+l+N/A
log1p-defineN/A
log1p-lowering-log1p.f64N/A
neg-mul-1N/A
*-commutativeN/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
+-commutativeN/A
metadata-evalN/A
sub-negN/A
expm1-defineN/A
expm1-lowering-expm1.f6499.9%
Simplified99.9%
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
log1p-defineN/A
log1p-lowering-log1p.f64N/A
*-lowering-*.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f6499.8%
Applied egg-rr99.8%
clear-numN/A
associate-/r/N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
log1p-defineN/A
log1p-lowering-log1p.f64N/A
*-lowering-*.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f6499.8%
Applied egg-rr99.8%
Taylor expanded in z around 0
/-lowering-/.f64N/A
Simplified66.5%
Taylor expanded in z around inf
cancel-sign-sub-invN/A
metadata-evalN/A
+-lowering-+.f64N/A
associate-*r/N/A
times-fracN/A
metadata-evalN/A
associate-*r/N/A
*-lowering-*.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
--lowering--.f6468.8%
Simplified68.8%
if -6.5e8 < z Initial program 49.5%
--lowering--.f64N/A
/-lowering-/.f64N/A
sub-negN/A
associate-+l+N/A
log1p-defineN/A
log1p-lowering-log1p.f64N/A
neg-mul-1N/A
*-commutativeN/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
+-commutativeN/A
metadata-evalN/A
sub-negN/A
expm1-defineN/A
expm1-lowering-expm1.f6498.2%
Simplified98.2%
Taylor expanded in z around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified68.7%
Taylor expanded in y around 0
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
associate-/l*N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6487.1%
Simplified87.1%
Final simplification81.2%
(FPCore (x y z t) :precision binary64 (+ x (/ 1.0 (- (/ (- (* t 0.5) (/ t z)) y) (* t 0.5)))))
double code(double x, double y, double z, double t) {
return x + (1.0 / ((((t * 0.5) - (t / z)) / y) - (t * 0.5)));
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
code = x + (1.0d0 / ((((t * 0.5d0) - (t / z)) / y) - (t * 0.5d0)))
end function
public static double code(double x, double y, double z, double t) {
return x + (1.0 / ((((t * 0.5) - (t / z)) / y) - (t * 0.5)));
}
def code(x, y, z, t): return x + (1.0 / ((((t * 0.5) - (t / z)) / y) - (t * 0.5)))
function code(x, y, z, t) return Float64(x + Float64(1.0 / Float64(Float64(Float64(Float64(t * 0.5) - Float64(t / z)) / y) - Float64(t * 0.5)))) end
function tmp = code(x, y, z, t) tmp = x + (1.0 / ((((t * 0.5) - (t / z)) / y) - (t * 0.5))); end
code[x_, y_, z_, t_] := N[(x + N[(1.0 / N[(N[(N[(N[(t * 0.5), $MachinePrecision] - N[(t / z), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision] - N[(t * 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x + \frac{1}{\frac{t \cdot 0.5 - \frac{t}{z}}{y} - t \cdot 0.5}
\end{array}
Initial program 60.4%
--lowering--.f64N/A
/-lowering-/.f64N/A
sub-negN/A
associate-+l+N/A
log1p-defineN/A
log1p-lowering-log1p.f64N/A
neg-mul-1N/A
*-commutativeN/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
+-commutativeN/A
metadata-evalN/A
sub-negN/A
expm1-defineN/A
expm1-lowering-expm1.f6498.8%
Simplified98.8%
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
log1p-defineN/A
log1p-lowering-log1p.f64N/A
*-lowering-*.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f6498.8%
Applied egg-rr98.8%
clear-numN/A
associate-/r/N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
log1p-defineN/A
log1p-lowering-log1p.f64N/A
*-lowering-*.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f6498.7%
Applied egg-rr98.7%
Taylor expanded in z around 0
/-lowering-/.f64N/A
Simplified77.8%
Taylor expanded in y around -inf
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f6483.0%
Simplified83.0%
Final simplification83.0%
(FPCore (x y z t) :precision binary64 (if (<= z -2.1e-17) (+ x (* (/ 2.0 t) (/ y (- 1.0 y)))) (- x (* y (/ z t)))))
double code(double x, double y, double z, double t) {
double tmp;
if (z <= -2.1e-17) {
tmp = x + ((2.0 / t) * (y / (1.0 - y)));
} else {
tmp = x - (y * (z / t));
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: tmp
if (z <= (-2.1d-17)) then
tmp = x + ((2.0d0 / t) * (y / (1.0d0 - y)))
else
tmp = x - (y * (z / t))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double tmp;
if (z <= -2.1e-17) {
tmp = x + ((2.0 / t) * (y / (1.0 - y)));
} else {
tmp = x - (y * (z / t));
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if z <= -2.1e-17: tmp = x + ((2.0 / t) * (y / (1.0 - y))) else: tmp = x - (y * (z / t)) return tmp
function code(x, y, z, t) tmp = 0.0 if (z <= -2.1e-17) tmp = Float64(x + Float64(Float64(2.0 / t) * Float64(y / Float64(1.0 - y)))); else tmp = Float64(x - Float64(y * Float64(z / t))); end return tmp end
function tmp_2 = code(x, y, z, t) tmp = 0.0; if (z <= -2.1e-17) tmp = x + ((2.0 / t) * (y / (1.0 - y))); else tmp = x - (y * (z / t)); end tmp_2 = tmp; end
code[x_, y_, z_, t_] := If[LessEqual[z, -2.1e-17], N[(x + N[(N[(2.0 / t), $MachinePrecision] * N[(y / N[(1.0 - y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x - N[(y * N[(z / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -2.1 \cdot 10^{-17}:\\
\;\;\;\;x + \frac{2}{t} \cdot \frac{y}{1 - y}\\
\mathbf{else}:\\
\;\;\;\;x - y \cdot \frac{z}{t}\\
\end{array}
\end{array}
if z < -2.09999999999999992e-17Initial program 81.2%
--lowering--.f64N/A
/-lowering-/.f64N/A
sub-negN/A
associate-+l+N/A
log1p-defineN/A
log1p-lowering-log1p.f64N/A
neg-mul-1N/A
*-commutativeN/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
+-commutativeN/A
metadata-evalN/A
sub-negN/A
expm1-defineN/A
expm1-lowering-expm1.f6499.9%
Simplified99.9%
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
log1p-defineN/A
log1p-lowering-log1p.f64N/A
*-lowering-*.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f6499.8%
Applied egg-rr99.8%
clear-numN/A
associate-/r/N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
log1p-defineN/A
log1p-lowering-log1p.f64N/A
*-lowering-*.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f6499.8%
Applied egg-rr99.8%
Taylor expanded in z around 0
/-lowering-/.f64N/A
Simplified64.8%
Taylor expanded in z around inf
cancel-sign-sub-invN/A
metadata-evalN/A
+-lowering-+.f64N/A
associate-*r/N/A
times-fracN/A
metadata-evalN/A
associate-*r/N/A
*-lowering-*.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
--lowering--.f6467.0%
Simplified67.0%
if -2.09999999999999992e-17 < z Initial program 49.4%
--lowering--.f64N/A
/-lowering-/.f64N/A
sub-negN/A
associate-+l+N/A
log1p-defineN/A
log1p-lowering-log1p.f64N/A
neg-mul-1N/A
*-commutativeN/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
+-commutativeN/A
metadata-evalN/A
sub-negN/A
expm1-defineN/A
expm1-lowering-expm1.f6498.2%
Simplified98.2%
Taylor expanded in z around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified69.7%
Taylor expanded in z around 0
Simplified87.2%
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
/-lowering-/.f6488.4%
Applied egg-rr88.4%
(FPCore (x y z t) :precision binary64 (if (<= z -3e+16) x (- x (* y (/ z t)))))
double code(double x, double y, double z, double t) {
double tmp;
if (z <= -3e+16) {
tmp = x;
} else {
tmp = x - (y * (z / t));
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: tmp
if (z <= (-3d+16)) then
tmp = x
else
tmp = x - (y * (z / t))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double tmp;
if (z <= -3e+16) {
tmp = x;
} else {
tmp = x - (y * (z / t));
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if z <= -3e+16: tmp = x else: tmp = x - (y * (z / t)) return tmp
function code(x, y, z, t) tmp = 0.0 if (z <= -3e+16) tmp = x; else tmp = Float64(x - Float64(y * Float64(z / t))); end return tmp end
function tmp_2 = code(x, y, z, t) tmp = 0.0; if (z <= -3e+16) tmp = x; else tmp = x - (y * (z / t)); end tmp_2 = tmp; end
code[x_, y_, z_, t_] := If[LessEqual[z, -3e+16], x, N[(x - N[(y * N[(z / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -3 \cdot 10^{+16}:\\
\;\;\;\;x\\
\mathbf{else}:\\
\;\;\;\;x - y \cdot \frac{z}{t}\\
\end{array}
\end{array}
if z < -3e16Initial program 84.0%
--lowering--.f64N/A
/-lowering-/.f64N/A
sub-negN/A
associate-+l+N/A
log1p-defineN/A
log1p-lowering-log1p.f64N/A
neg-mul-1N/A
*-commutativeN/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
+-commutativeN/A
metadata-evalN/A
sub-negN/A
expm1-defineN/A
expm1-lowering-expm1.f6499.9%
Simplified99.9%
Taylor expanded in x around inf
Simplified65.9%
if -3e16 < z Initial program 49.5%
--lowering--.f64N/A
/-lowering-/.f64N/A
sub-negN/A
associate-+l+N/A
log1p-defineN/A
log1p-lowering-log1p.f64N/A
neg-mul-1N/A
*-commutativeN/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
+-commutativeN/A
metadata-evalN/A
sub-negN/A
expm1-defineN/A
expm1-lowering-expm1.f6498.2%
Simplified98.2%
Taylor expanded in z around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified68.0%
Taylor expanded in z around 0
Simplified84.7%
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
/-lowering-/.f6485.8%
Applied egg-rr85.8%
(FPCore (x y z t) :precision binary64 x)
double code(double x, double y, double z, double t) {
return x;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
code = x
end function
public static double code(double x, double y, double z, double t) {
return x;
}
def code(x, y, z, t): return x
function code(x, y, z, t) return x end
function tmp = code(x, y, z, t) tmp = x; end
code[x_, y_, z_, t_] := x
\begin{array}{l}
\\
x
\end{array}
Initial program 60.4%
--lowering--.f64N/A
/-lowering-/.f64N/A
sub-negN/A
associate-+l+N/A
log1p-defineN/A
log1p-lowering-log1p.f64N/A
neg-mul-1N/A
*-commutativeN/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
+-commutativeN/A
metadata-evalN/A
sub-negN/A
expm1-defineN/A
expm1-lowering-expm1.f6498.8%
Simplified98.8%
Taylor expanded in x around inf
Simplified69.2%
(FPCore (x y z t)
:precision binary64
(let* ((t_1 (/ (- 0.5) (* y t))))
(if (< z -2.8874623088207947e+119)
(- (- x (/ t_1 (* z z))) (* t_1 (/ (/ 2.0 z) (* z z))))
(- x (/ (log (+ 1.0 (* z y))) t)))))
double code(double x, double y, double z, double t) {
double t_1 = -0.5 / (y * t);
double tmp;
if (z < -2.8874623088207947e+119) {
tmp = (x - (t_1 / (z * z))) - (t_1 * ((2.0 / z) / (z * z)));
} else {
tmp = x - (log((1.0 + (z * y))) / t);
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: t_1
real(8) :: tmp
t_1 = -0.5d0 / (y * t)
if (z < (-2.8874623088207947d+119)) then
tmp = (x - (t_1 / (z * z))) - (t_1 * ((2.0d0 / z) / (z * z)))
else
tmp = x - (log((1.0d0 + (z * y))) / t)
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double t_1 = -0.5 / (y * t);
double tmp;
if (z < -2.8874623088207947e+119) {
tmp = (x - (t_1 / (z * z))) - (t_1 * ((2.0 / z) / (z * z)));
} else {
tmp = x - (Math.log((1.0 + (z * y))) / t);
}
return tmp;
}
def code(x, y, z, t): t_1 = -0.5 / (y * t) tmp = 0 if z < -2.8874623088207947e+119: tmp = (x - (t_1 / (z * z))) - (t_1 * ((2.0 / z) / (z * z))) else: tmp = x - (math.log((1.0 + (z * y))) / t) return tmp
function code(x, y, z, t) t_1 = Float64(Float64(-0.5) / Float64(y * t)) tmp = 0.0 if (z < -2.8874623088207947e+119) tmp = Float64(Float64(x - Float64(t_1 / Float64(z * z))) - Float64(t_1 * Float64(Float64(2.0 / z) / Float64(z * z)))); else tmp = Float64(x - Float64(log(Float64(1.0 + Float64(z * y))) / t)); end return tmp end
function tmp_2 = code(x, y, z, t) t_1 = -0.5 / (y * t); tmp = 0.0; if (z < -2.8874623088207947e+119) tmp = (x - (t_1 / (z * z))) - (t_1 * ((2.0 / z) / (z * z))); else tmp = x - (log((1.0 + (z * y))) / t); end tmp_2 = tmp; end
code[x_, y_, z_, t_] := Block[{t$95$1 = N[((-0.5) / N[(y * t), $MachinePrecision]), $MachinePrecision]}, If[Less[z, -2.8874623088207947e+119], N[(N[(x - N[(t$95$1 / N[(z * z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(t$95$1 * N[(N[(2.0 / z), $MachinePrecision] / N[(z * z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x - N[(N[Log[N[(1.0 + N[(z * y), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{-0.5}{y \cdot t}\\
\mathbf{if}\;z < -2.8874623088207947 \cdot 10^{+119}:\\
\;\;\;\;\left(x - \frac{t\_1}{z \cdot z}\right) - t\_1 \cdot \frac{\frac{2}{z}}{z \cdot z}\\
\mathbf{else}:\\
\;\;\;\;x - \frac{\log \left(1 + z \cdot y\right)}{t}\\
\end{array}
\end{array}
herbie shell --seed 2024160
(FPCore (x y z t)
:name "System.Random.MWC.Distributions:truncatedExp from mwc-random-0.13.3.2"
:precision binary64
:alt
(! :herbie-platform default (if (< z -288746230882079470000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000) (- (- x (/ (/ (- 1/2) (* y t)) (* z z))) (* (/ (- 1/2) (* y t)) (/ (/ 2 z) (* z z)))) (- x (/ (log (+ 1 (* z y))) t))))
(- x (/ (log (+ (- 1.0 y) (* y (exp z)))) t)))