
(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 16 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 (/ 1.0 (expm1 z)))) t)))
double code(double x, double y, double z, double t) {
return x - (log1p((y / (1.0 / expm1(z)))) / t);
}
public static double code(double x, double y, double z, double t) {
return x - (Math.log1p((y / (1.0 / Math.expm1(z)))) / t);
}
def code(x, y, z, t): return x - (math.log1p((y / (1.0 / math.expm1(z)))) / t)
function code(x, y, z, t) return Float64(x - Float64(log1p(Float64(y / Float64(1.0 / expm1(z)))) / t)) end
code[x_, y_, z_, t_] := N[(x - N[(N[Log[1 + N[(y / N[(1.0 / N[(Exp[z] - 1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x - \frac{\mathsf{log1p}\left(\frac{y}{\frac{1}{\mathsf{expm1}\left(z\right)}}\right)}{t}
\end{array}
Initial program 57.9%
--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.f6497.8%
Simplified97.8%
flip--N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
clear-numN/A
flip--N/A
/-lowering-/.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f6497.8%
Applied egg-rr97.8%
(FPCore (x y z t) :precision binary64 (if (<= (exp z) 0.0) (+ 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 (exp(z) <= 0.0) {
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 (Math.exp(z) <= 0.0) {
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 math.exp(z) <= 0.0: 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 (exp(z) <= 0.0) 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[N[Exp[z], $MachinePrecision], 0.0], 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}\;e^{z} \leq 0:\\
\;\;\;\;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 (exp.f64 z) < 0.0Initial program 81.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.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.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.f6484.2%
Simplified84.2%
if 0.0 < (exp.f64 z) Initial program 49.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.f6497.0%
Simplified97.0%
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 simplification94.3%
(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 57.9%
--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.f6497.8%
Simplified97.8%
(FPCore (x y z t)
:precision binary64
(if (<= z -0.088)
(+ x (/ -1.0 (/ (+ (* 0.5 (* y t)) (/ t (expm1 z))) y)))
(-
x
(/
(log1p
(/
y
(/
1.0
(*
z
(+
1.0
(*
z
(+
0.5
(* z (+ 0.16666666666666666 (* z 0.041666666666666664))))))))))
t))))
double code(double x, double y, double z, double t) {
double tmp;
if (z <= -0.088) {
tmp = x + (-1.0 / (((0.5 * (y * t)) + (t / expm1(z))) / y));
} else {
tmp = x - (log1p((y / (1.0 / (z * (1.0 + (z * (0.5 + (z * (0.16666666666666666 + (z * 0.041666666666666664)))))))))) / t);
}
return tmp;
}
public static double code(double x, double y, double z, double t) {
double tmp;
if (z <= -0.088) {
tmp = x + (-1.0 / (((0.5 * (y * t)) + (t / Math.expm1(z))) / y));
} else {
tmp = x - (Math.log1p((y / (1.0 / (z * (1.0 + (z * (0.5 + (z * (0.16666666666666666 + (z * 0.041666666666666664)))))))))) / t);
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if z <= -0.088: tmp = x + (-1.0 / (((0.5 * (y * t)) + (t / math.expm1(z))) / y)) else: tmp = x - (math.log1p((y / (1.0 / (z * (1.0 + (z * (0.5 + (z * (0.16666666666666666 + (z * 0.041666666666666664)))))))))) / t) return tmp
function code(x, y, z, t) tmp = 0.0 if (z <= -0.088) 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(y / Float64(1.0 / Float64(z * Float64(1.0 + Float64(z * Float64(0.5 + Float64(z * Float64(0.16666666666666666 + Float64(z * 0.041666666666666664)))))))))) / t)); end return tmp end
code[x_, y_, z_, t_] := If[LessEqual[z, -0.088], 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[(y / N[(1.0 / 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] / t), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -0.088:\\
\;\;\;\;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(\frac{y}{\frac{1}{z \cdot \left(1 + z \cdot \left(0.5 + z \cdot \left(0.16666666666666666 + z \cdot 0.041666666666666664\right)\right)\right)}}\right)}{t}\\
\end{array}
\end{array}
if z < -0.087999999999999995Initial program 81.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.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.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.f6484.2%
Simplified84.2%
if -0.087999999999999995 < z Initial program 49.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.f6497.0%
Simplified97.0%
flip--N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
clear-numN/A
flip--N/A
/-lowering-/.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f6497.0%
Applied egg-rr97.0%
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.1%
Simplified98.1%
Final simplification94.3%
(FPCore (x y z t)
:precision binary64
(if (<= y -7.6e+38)
(- x (/ (log1p (* y z)) t))
(if (<= y 2e+53)
(- x (* y (/ (expm1 z) t)))
(- x (/ (log1p (* z (+ y (* y (* z 0.5))))) t)))))
double code(double x, double y, double z, double t) {
double tmp;
if (y <= -7.6e+38) {
tmp = x - (log1p((y * z)) / t);
} else if (y <= 2e+53) {
tmp = x - (y * (expm1(z) / t));
} 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 (y <= -7.6e+38) {
tmp = x - (Math.log1p((y * z)) / t);
} else if (y <= 2e+53) {
tmp = x - (y * (Math.expm1(z) / t));
} else {
tmp = x - (Math.log1p((z * (y + (y * (z * 0.5))))) / t);
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if y <= -7.6e+38: tmp = x - (math.log1p((y * z)) / t) elif y <= 2e+53: tmp = x - (y * (math.expm1(z) / t)) else: tmp = x - (math.log1p((z * (y + (y * (z * 0.5))))) / t) return tmp
function code(x, y, z, t) tmp = 0.0 if (y <= -7.6e+38) tmp = Float64(x - Float64(log1p(Float64(y * z)) / t)); elseif (y <= 2e+53) tmp = Float64(x - Float64(y * Float64(expm1(z) / t))); 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[y, -7.6e+38], N[(x - N[(N[Log[1 + N[(y * z), $MachinePrecision]], $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 2e+53], N[(x - N[(y * N[(N[(Exp[z] - 1), $MachinePrecision] / t), $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}\;y \leq -7.6 \cdot 10^{+38}:\\
\;\;\;\;x - \frac{\mathsf{log1p}\left(y \cdot z\right)}{t}\\
\mathbf{elif}\;y \leq 2 \cdot 10^{+53}:\\
\;\;\;\;x - y \cdot \frac{\mathsf{expm1}\left(z\right)}{t}\\
\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 y < -7.5999999999999996e38Initial program 44.8%
--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 z around 0
*-lowering-*.f6475.2%
Simplified75.2%
if -7.5999999999999996e38 < y < 2e53Initial program 72.9%
--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
/-lowering-/.f64N/A
*-lowering-*.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f6498.4%
Simplified98.4%
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f6499.9%
Applied egg-rr99.9%
if 2e53 < y Initial program 3.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.f6490.8%
Simplified90.8%
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-*.f6497.3%
Simplified97.3%
Final simplification94.1%
(FPCore (x y z t)
:precision binary64
(if (<= z -0.0295)
(+ x (/ -1.0 (/ (+ (* 0.5 (* y t)) (/ t (expm1 z))) y)))
(-
x
(/ (log1p (* z (+ y (* z (* y (+ 0.5 (* z 0.16666666666666666))))))) t))))
double code(double x, double y, double z, double t) {
double tmp;
if (z <= -0.0295) {
tmp = x + (-1.0 / (((0.5 * (y * t)) + (t / expm1(z))) / y));
} else {
tmp = x - (log1p((z * (y + (z * (y * (0.5 + (z * 0.16666666666666666))))))) / t);
}
return tmp;
}
public static double code(double x, double y, double z, double t) {
double tmp;
if (z <= -0.0295) {
tmp = x + (-1.0 / (((0.5 * (y * t)) + (t / Math.expm1(z))) / y));
} else {
tmp = x - (Math.log1p((z * (y + (z * (y * (0.5 + (z * 0.16666666666666666))))))) / t);
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if z <= -0.0295: tmp = x + (-1.0 / (((0.5 * (y * t)) + (t / math.expm1(z))) / y)) else: tmp = x - (math.log1p((z * (y + (z * (y * (0.5 + (z * 0.16666666666666666))))))) / t) return tmp
function code(x, y, z, t) tmp = 0.0 if (z <= -0.0295) 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(y * Float64(0.5 + Float64(z * 0.16666666666666666))))))) / t)); end return tmp end
code[x_, y_, z_, t_] := If[LessEqual[z, -0.0295], 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[(y * N[(0.5 + N[(z * 0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -0.0295:\\
\;\;\;\;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(0.5 + z \cdot 0.16666666666666666\right)\right)\right)\right)}{t}\\
\end{array}
\end{array}
if z < -0.029499999999999998Initial program 81.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.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.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.f6484.2%
Simplified84.2%
if -0.029499999999999998 < z Initial program 49.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.f6497.0%
Simplified97.0%
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.1%
Simplified98.1%
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
distribute-lft-outN/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f6498.1%
Applied egg-rr98.1%
Final simplification94.3%
(FPCore (x y z t) :precision binary64 (let* ((t_1 (- x (/ (log1p (* y z)) t)))) (if (<= y -7.5e+38) t_1 (if (<= y 2e+49) (- x (* y (/ (expm1 z) t))) t_1))))
double code(double x, double y, double z, double t) {
double t_1 = x - (log1p((y * z)) / t);
double tmp;
if (y <= -7.5e+38) {
tmp = t_1;
} else if (y <= 2e+49) {
tmp = x - (y * (expm1(z) / t));
} else {
tmp = t_1;
}
return tmp;
}
public static double code(double x, double y, double z, double t) {
double t_1 = x - (Math.log1p((y * z)) / t);
double tmp;
if (y <= -7.5e+38) {
tmp = t_1;
} else if (y <= 2e+49) {
tmp = x - (y * (Math.expm1(z) / t));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t): t_1 = x - (math.log1p((y * z)) / t) tmp = 0 if y <= -7.5e+38: tmp = t_1 elif y <= 2e+49: tmp = x - (y * (math.expm1(z) / t)) else: tmp = t_1 return tmp
function code(x, y, z, t) t_1 = Float64(x - Float64(log1p(Float64(y * z)) / t)) tmp = 0.0 if (y <= -7.5e+38) tmp = t_1; elseif (y <= 2e+49) tmp = Float64(x - Float64(y * Float64(expm1(z) / t))); else tmp = t_1; end return tmp end
code[x_, y_, z_, t_] := Block[{t$95$1 = N[(x - N[(N[Log[1 + N[(y * z), $MachinePrecision]], $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -7.5e+38], t$95$1, If[LessEqual[y, 2e+49], N[(x - N[(y * N[(N[(Exp[z] - 1), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x - \frac{\mathsf{log1p}\left(y \cdot z\right)}{t}\\
\mathbf{if}\;y \leq -7.5 \cdot 10^{+38}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 2 \cdot 10^{+49}:\\
\;\;\;\;x - y \cdot \frac{\mathsf{expm1}\left(z\right)}{t}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -7.4999999999999999e38 or 1.99999999999999989e49 < y Initial program 29.8%
--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.f6496.6%
Simplified96.6%
Taylor expanded in z around 0
*-lowering-*.f6483.1%
Simplified83.1%
if -7.4999999999999999e38 < y < 1.99999999999999989e49Initial program 72.9%
--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
/-lowering-/.f64N/A
*-lowering-*.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f6498.4%
Simplified98.4%
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f6499.9%
Applied egg-rr99.9%
Final simplification94.1%
(FPCore (x y z t) :precision binary64 (let* ((t_1 (+ x (/ 1.0 (/ (/ (- (* z (* 0.5 (- t (* y t)))) t) z) y))))) (if (<= y -7e+38) t_1 (if (<= y 9.5e+60) (- x (* y (/ (expm1 z) t))) t_1))))
double code(double x, double y, double z, double t) {
double t_1 = x + (1.0 / ((((z * (0.5 * (t - (y * t)))) - t) / z) / y));
double tmp;
if (y <= -7e+38) {
tmp = t_1;
} else if (y <= 9.5e+60) {
tmp = x - (y * (expm1(z) / t));
} else {
tmp = t_1;
}
return tmp;
}
public static double code(double x, double y, double z, double t) {
double t_1 = x + (1.0 / ((((z * (0.5 * (t - (y * t)))) - t) / z) / y));
double tmp;
if (y <= -7e+38) {
tmp = t_1;
} else if (y <= 9.5e+60) {
tmp = x - (y * (Math.expm1(z) / t));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t): t_1 = x + (1.0 / ((((z * (0.5 * (t - (y * t)))) - t) / z) / y)) tmp = 0 if y <= -7e+38: tmp = t_1 elif y <= 9.5e+60: tmp = x - (y * (math.expm1(z) / t)) else: tmp = t_1 return tmp
function code(x, y, z, t) t_1 = Float64(x + Float64(1.0 / Float64(Float64(Float64(Float64(z * Float64(0.5 * Float64(t - Float64(y * t)))) - t) / z) / y))) tmp = 0.0 if (y <= -7e+38) tmp = t_1; elseif (y <= 9.5e+60) tmp = Float64(x - Float64(y * Float64(expm1(z) / t))); else tmp = t_1; end return tmp end
code[x_, y_, z_, t_] := Block[{t$95$1 = N[(x + N[(1.0 / N[(N[(N[(N[(z * N[(0.5 * N[(t - N[(y * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - t), $MachinePrecision] / z), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -7e+38], t$95$1, If[LessEqual[y, 9.5e+60], N[(x - N[(y * N[(N[(Exp[z] - 1), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x + \frac{1}{\frac{\frac{z \cdot \left(0.5 \cdot \left(t - y \cdot t\right)\right) - t}{z}}{y}}\\
\mathbf{if}\;y \leq -7 \cdot 10^{+38}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 9.5 \cdot 10^{+60}:\\
\;\;\;\;x - y \cdot \frac{\mathsf{expm1}\left(z\right)}{t}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -7.00000000000000003e38 or 9.49999999999999988e60 < y Initial program 30.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.f6496.5%
Simplified96.5%
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.f6496.5%
Applied egg-rr96.5%
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.f6459.6%
Simplified59.6%
Taylor expanded in z around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
distribute-lft-out--N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6467.3%
Simplified67.3%
if -7.00000000000000003e38 < y < 9.49999999999999988e60Initial program 72.1%
--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
/-lowering-/.f64N/A
*-lowering-*.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f6498.4%
Simplified98.4%
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f6499.9%
Applied egg-rr99.9%
Final simplification88.8%
(FPCore (x y z t)
:precision binary64
(if (<= z -6.5e-147)
(+ x (/ 1.0 (/ (- (* z (* 0.5 (- (/ t y) t))) (/ t y)) z)))
(if (<= z 5e-224)
(- x (* y (/ z t)))
(+ x (/ -1.0 (/ (+ (* 0.5 (* y t)) (/ t z)) y))))))
double code(double x, double y, double z, double t) {
double tmp;
if (z <= -6.5e-147) {
tmp = x + (1.0 / (((z * (0.5 * ((t / y) - t))) - (t / y)) / z));
} else if (z <= 5e-224) {
tmp = x - (y * (z / t));
} else {
tmp = x + (-1.0 / (((0.5 * (y * t)) + (t / z)) / y));
}
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 <= (-6.5d-147)) then
tmp = x + (1.0d0 / (((z * (0.5d0 * ((t / y) - t))) - (t / y)) / z))
else if (z <= 5d-224) then
tmp = x - (y * (z / t))
else
tmp = x + ((-1.0d0) / (((0.5d0 * (y * t)) + (t / z)) / y))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double tmp;
if (z <= -6.5e-147) {
tmp = x + (1.0 / (((z * (0.5 * ((t / y) - t))) - (t / y)) / z));
} else if (z <= 5e-224) {
tmp = x - (y * (z / t));
} else {
tmp = x + (-1.0 / (((0.5 * (y * t)) + (t / z)) / y));
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if z <= -6.5e-147: tmp = x + (1.0 / (((z * (0.5 * ((t / y) - t))) - (t / y)) / z)) elif z <= 5e-224: tmp = x - (y * (z / t)) else: tmp = x + (-1.0 / (((0.5 * (y * t)) + (t / z)) / y)) return tmp
function code(x, y, z, t) tmp = 0.0 if (z <= -6.5e-147) tmp = Float64(x + Float64(1.0 / Float64(Float64(Float64(z * Float64(0.5 * Float64(Float64(t / y) - t))) - Float64(t / y)) / z))); elseif (z <= 5e-224) tmp = Float64(x - Float64(y * Float64(z / t))); else tmp = Float64(x + Float64(-1.0 / Float64(Float64(Float64(0.5 * Float64(y * t)) + Float64(t / z)) / y))); end return tmp end
function tmp_2 = code(x, y, z, t) tmp = 0.0; if (z <= -6.5e-147) tmp = x + (1.0 / (((z * (0.5 * ((t / y) - t))) - (t / y)) / z)); elseif (z <= 5e-224) tmp = x - (y * (z / t)); else tmp = x + (-1.0 / (((0.5 * (y * t)) + (t / z)) / y)); end tmp_2 = tmp; end
code[x_, y_, z_, t_] := If[LessEqual[z, -6.5e-147], 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], If[LessEqual[z, 5e-224], N[(x - N[(y * N[(z / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x + N[(-1.0 / N[(N[(N[(0.5 * N[(y * t), $MachinePrecision]), $MachinePrecision] + N[(t / z), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -6.5 \cdot 10^{-147}:\\
\;\;\;\;x + \frac{1}{\frac{z \cdot \left(0.5 \cdot \left(\frac{t}{y} - t\right)\right) - \frac{t}{y}}{z}}\\
\mathbf{elif}\;z \leq 5 \cdot 10^{-224}:\\
\;\;\;\;x - y \cdot \frac{z}{t}\\
\mathbf{else}:\\
\;\;\;\;x + \frac{-1}{\frac{0.5 \cdot \left(y \cdot t\right) + \frac{t}{z}}{y}}\\
\end{array}
\end{array}
if z < -6.49999999999999967e-147Initial program 69.9%
--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%
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.f6484.0%
Simplified84.0%
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-/.f6474.2%
Simplified74.2%
if -6.49999999999999967e-147 < z < 4.9999999999999999e-224Initial program 48.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.0%
Simplified99.0%
Taylor expanded in y around 0
/-lowering-/.f64N/A
*-lowering-*.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f6497.8%
Simplified97.8%
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f6498.6%
Applied egg-rr98.6%
Taylor expanded in z around 0
/-lowering-/.f6498.6%
Simplified98.6%
if 4.9999999999999999e-224 < z Initial program 49.1%
--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.f6493.2%
Simplified93.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.f6493.2%
Applied egg-rr93.2%
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.f6486.6%
Simplified86.6%
Taylor expanded in z around 0
/-lowering-/.f6486.6%
Simplified86.6%
Final simplification84.7%
(FPCore (x y z t)
:precision binary64
(if (<= z -4e+15)
x
(if (<= z 8.2e-227)
(+ x (/ -1.0 (/ (/ t z) y)))
(+ x (/ -1.0 (/ (+ (* 0.5 (* y t)) (/ t z)) y))))))
double code(double x, double y, double z, double t) {
double tmp;
if (z <= -4e+15) {
tmp = x;
} else if (z <= 8.2e-227) {
tmp = x + (-1.0 / ((t / z) / y));
} else {
tmp = x + (-1.0 / (((0.5 * (y * t)) + (t / z)) / y));
}
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 <= (-4d+15)) then
tmp = x
else if (z <= 8.2d-227) then
tmp = x + ((-1.0d0) / ((t / z) / y))
else
tmp = x + ((-1.0d0) / (((0.5d0 * (y * t)) + (t / z)) / y))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double tmp;
if (z <= -4e+15) {
tmp = x;
} else if (z <= 8.2e-227) {
tmp = x + (-1.0 / ((t / z) / y));
} else {
tmp = x + (-1.0 / (((0.5 * (y * t)) + (t / z)) / y));
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if z <= -4e+15: tmp = x elif z <= 8.2e-227: tmp = x + (-1.0 / ((t / z) / y)) else: tmp = x + (-1.0 / (((0.5 * (y * t)) + (t / z)) / y)) return tmp
function code(x, y, z, t) tmp = 0.0 if (z <= -4e+15) tmp = x; elseif (z <= 8.2e-227) tmp = Float64(x + Float64(-1.0 / Float64(Float64(t / z) / y))); else tmp = Float64(x + Float64(-1.0 / Float64(Float64(Float64(0.5 * Float64(y * t)) + Float64(t / z)) / y))); end return tmp end
function tmp_2 = code(x, y, z, t) tmp = 0.0; if (z <= -4e+15) tmp = x; elseif (z <= 8.2e-227) tmp = x + (-1.0 / ((t / z) / y)); else tmp = x + (-1.0 / (((0.5 * (y * t)) + (t / z)) / y)); end tmp_2 = tmp; end
code[x_, y_, z_, t_] := If[LessEqual[z, -4e+15], x, If[LessEqual[z, 8.2e-227], N[(x + N[(-1.0 / N[(N[(t / z), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x + N[(-1.0 / N[(N[(N[(0.5 * N[(y * t), $MachinePrecision]), $MachinePrecision] + N[(t / z), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -4 \cdot 10^{+15}:\\
\;\;\;\;x\\
\mathbf{elif}\;z \leq 8.2 \cdot 10^{-227}:\\
\;\;\;\;x + \frac{-1}{\frac{\frac{t}{z}}{y}}\\
\mathbf{else}:\\
\;\;\;\;x + \frac{-1}{\frac{0.5 \cdot \left(y \cdot t\right) + \frac{t}{z}}{y}}\\
\end{array}
\end{array}
if z < -4e15Initial 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.f64100.0%
Simplified100.0%
Taylor expanded in x around inf
Simplified66.1%
if -4e15 < z < 8.20000000000000018e-227Initial program 48.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.3%
Simplified99.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.f6499.2%
Applied egg-rr99.2%
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.2%
Simplified87.2%
Taylor expanded in z around 0
/-lowering-/.f6491.1%
Simplified91.1%
if 8.20000000000000018e-227 < z Initial program 48.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.f6493.3%
Simplified93.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.f6493.3%
Applied egg-rr93.3%
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.f6486.8%
Simplified86.8%
Taylor expanded in z around 0
/-lowering-/.f6486.8%
Simplified86.8%
Final simplification83.2%
(FPCore (x y z t) :precision binary64 (+ x (/ 1.0 (/ (/ (- (* z (* 0.5 (- t (* y t)))) t) z) y))))
double code(double x, double y, double z, double t) {
return x + (1.0 / ((((z * (0.5 * (t - (y * t)))) - t) / z) / y));
}
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 / ((((z * (0.5d0 * (t - (y * t)))) - t) / z) / y))
end function
public static double code(double x, double y, double z, double t) {
return x + (1.0 / ((((z * (0.5 * (t - (y * t)))) - t) / z) / y));
}
def code(x, y, z, t): return x + (1.0 / ((((z * (0.5 * (t - (y * t)))) - t) / z) / y))
function code(x, y, z, t) return Float64(x + Float64(1.0 / Float64(Float64(Float64(Float64(z * Float64(0.5 * Float64(t - Float64(y * t)))) - t) / z) / y))) end
function tmp = code(x, y, z, t) tmp = x + (1.0 / ((((z * (0.5 * (t - (y * t)))) - t) / z) / y)); end
code[x_, y_, z_, t_] := N[(x + N[(1.0 / N[(N[(N[(N[(z * N[(0.5 * N[(t - N[(y * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - t), $MachinePrecision] / z), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x + \frac{1}{\frac{\frac{z \cdot \left(0.5 \cdot \left(t - y \cdot t\right)\right) - t}{z}}{y}}
\end{array}
Initial program 57.9%
--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.f6497.8%
Simplified97.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.f6497.7%
Applied egg-rr97.7%
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.f6486.2%
Simplified86.2%
Taylor expanded in z around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
distribute-lft-out--N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6484.5%
Simplified84.5%
Final simplification84.5%
(FPCore (x y z t) :precision binary64 (if (<= t -4.8e-244) x (if (<= t 5.8e-207) (/ (* y z) (- 0.0 t)) x)))
double code(double x, double y, double z, double t) {
double tmp;
if (t <= -4.8e-244) {
tmp = x;
} else if (t <= 5.8e-207) {
tmp = (y * z) / (0.0 - t);
} else {
tmp = x;
}
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 (t <= (-4.8d-244)) then
tmp = x
else if (t <= 5.8d-207) then
tmp = (y * z) / (0.0d0 - t)
else
tmp = x
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double tmp;
if (t <= -4.8e-244) {
tmp = x;
} else if (t <= 5.8e-207) {
tmp = (y * z) / (0.0 - t);
} else {
tmp = x;
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if t <= -4.8e-244: tmp = x elif t <= 5.8e-207: tmp = (y * z) / (0.0 - t) else: tmp = x return tmp
function code(x, y, z, t) tmp = 0.0 if (t <= -4.8e-244) tmp = x; elseif (t <= 5.8e-207) tmp = Float64(Float64(y * z) / Float64(0.0 - t)); else tmp = x; end return tmp end
function tmp_2 = code(x, y, z, t) tmp = 0.0; if (t <= -4.8e-244) tmp = x; elseif (t <= 5.8e-207) tmp = (y * z) / (0.0 - t); else tmp = x; end tmp_2 = tmp; end
code[x_, y_, z_, t_] := If[LessEqual[t, -4.8e-244], x, If[LessEqual[t, 5.8e-207], N[(N[(y * z), $MachinePrecision] / N[(0.0 - t), $MachinePrecision]), $MachinePrecision], x]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -4.8 \cdot 10^{-244}:\\
\;\;\;\;x\\
\mathbf{elif}\;t \leq 5.8 \cdot 10^{-207}:\\
\;\;\;\;\frac{y \cdot z}{0 - t}\\
\mathbf{else}:\\
\;\;\;\;x\\
\end{array}
\end{array}
if t < -4.80000000000000032e-244 or 5.80000000000000022e-207 < t Initial program 64.1%
--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 x around inf
Simplified75.4%
if -4.80000000000000032e-244 < t < 5.80000000000000022e-207Initial program 17.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.f6495.0%
Simplified95.0%
Taylor expanded in z around 0
distribute-lft-inN/A
fma-defineN/A
*-commutativeN/A
associate-/l*N/A
associate-*r*N/A
*-commutativeN/A
fma-undefineN/A
distribute-lft-inN/A
+-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
associate-*r/N/A
associate-*r/N/A
Simplified53.6%
Taylor expanded in t around 0
mul-1-negN/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
Simplified48.3%
Taylor expanded in z around 0
Simplified49.6%
Final simplification72.0%
(FPCore (x y z t) :precision binary64 (if (<= z -960000000.0) x (+ x (/ -1.0 (/ (/ t z) y)))))
double code(double x, double y, double z, double t) {
double tmp;
if (z <= -960000000.0) {
tmp = x;
} else {
tmp = x + (-1.0 / ((t / z) / y));
}
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 <= (-960000000.0d0)) then
tmp = x
else
tmp = x + ((-1.0d0) / ((t / z) / y))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double tmp;
if (z <= -960000000.0) {
tmp = x;
} else {
tmp = x + (-1.0 / ((t / z) / y));
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if z <= -960000000.0: tmp = x else: tmp = x + (-1.0 / ((t / z) / y)) return tmp
function code(x, y, z, t) tmp = 0.0 if (z <= -960000000.0) tmp = x; else tmp = Float64(x + Float64(-1.0 / Float64(Float64(t / z) / y))); end return tmp end
function tmp_2 = code(x, y, z, t) tmp = 0.0; if (z <= -960000000.0) tmp = x; else tmp = x + (-1.0 / ((t / z) / y)); end tmp_2 = tmp; end
code[x_, y_, z_, t_] := If[LessEqual[z, -960000000.0], x, N[(x + N[(-1.0 / N[(N[(t / z), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -960000000:\\
\;\;\;\;x\\
\mathbf{else}:\\
\;\;\;\;x + \frac{-1}{\frac{\frac{t}{z}}{y}}\\
\end{array}
\end{array}
if z < -9.6e8Initial 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.f64100.0%
Simplified100.0%
Taylor expanded in x around inf
Simplified66.1%
if -9.6e8 < z Initial program 48.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.f6497.0%
Simplified97.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.f6496.9%
Applied egg-rr96.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.f6487.0%
Simplified87.0%
Taylor expanded in z around 0
/-lowering-/.f6486.3%
Simplified86.3%
Final simplification80.9%
(FPCore (x y z t) :precision binary64 (if (<= z -5.3e-32) x (- x (* y (/ z t)))))
double code(double x, double y, double z, double t) {
double tmp;
if (z <= -5.3e-32) {
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 <= (-5.3d-32)) 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 <= -5.3e-32) {
tmp = x;
} else {
tmp = x - (y * (z / t));
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if z <= -5.3e-32: tmp = x else: tmp = x - (y * (z / t)) return tmp
function code(x, y, z, t) tmp = 0.0 if (z <= -5.3e-32) 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 <= -5.3e-32) tmp = x; else tmp = x - (y * (z / t)); end tmp_2 = tmp; end
code[x_, y_, z_, t_] := If[LessEqual[z, -5.3e-32], x, N[(x - N[(y * N[(z / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -5.3 \cdot 10^{-32}:\\
\;\;\;\;x\\
\mathbf{else}:\\
\;\;\;\;x - y \cdot \frac{z}{t}\\
\end{array}
\end{array}
if z < -5.3000000000000001e-32Initial program 78.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%
Taylor expanded in x around inf
Simplified64.3%
if -5.3000000000000001e-32 < z Initial program 49.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.f6496.8%
Simplified96.8%
Taylor expanded in y around 0
/-lowering-/.f64N/A
*-lowering-*.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f6485.8%
Simplified85.8%
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f6487.7%
Applied egg-rr87.7%
Taylor expanded in z around 0
/-lowering-/.f6487.8%
Simplified87.8%
Final simplification80.6%
(FPCore (x y z t) :precision binary64 (if (<= z -2.25e-31) x (- x (* z (/ y t)))))
double code(double x, double y, double z, double t) {
double tmp;
if (z <= -2.25e-31) {
tmp = x;
} else {
tmp = x - (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) :: tmp
if (z <= (-2.25d-31)) then
tmp = x
else
tmp = x - (z * (y / t))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double tmp;
if (z <= -2.25e-31) {
tmp = x;
} else {
tmp = x - (z * (y / t));
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if z <= -2.25e-31: tmp = x else: tmp = x - (z * (y / t)) return tmp
function code(x, y, z, t) tmp = 0.0 if (z <= -2.25e-31) tmp = x; else tmp = Float64(x - Float64(z * Float64(y / t))); end return tmp end
function tmp_2 = code(x, y, z, t) tmp = 0.0; if (z <= -2.25e-31) tmp = x; else tmp = x - (z * (y / t)); end tmp_2 = tmp; end
code[x_, y_, z_, t_] := If[LessEqual[z, -2.25e-31], x, N[(x - N[(z * N[(y / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -2.25 \cdot 10^{-31}:\\
\;\;\;\;x\\
\mathbf{else}:\\
\;\;\;\;x - z \cdot \frac{y}{t}\\
\end{array}
\end{array}
if z < -2.2500000000000002e-31Initial program 78.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%
Taylor expanded in x around inf
Simplified64.3%
if -2.2500000000000002e-31 < z Initial program 49.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.f6496.8%
Simplified96.8%
Taylor expanded in z around 0
distribute-lft-inN/A
fma-defineN/A
*-commutativeN/A
associate-/l*N/A
associate-*r*N/A
*-commutativeN/A
fma-undefineN/A
distribute-lft-inN/A
+-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
associate-*r/N/A
associate-*r/N/A
Simplified74.1%
Taylor expanded in z around 0
/-lowering-/.f6481.3%
Simplified81.3%
(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 57.9%
--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.f6497.8%
Simplified97.8%
Taylor expanded in x around inf
Simplified68.5%
(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 2024138
(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)))