
(FPCore (x y z t a b) :precision binary64 (/ (* x (exp (- (+ (* y (log z)) (* (- t 1.0) (log a))) b))) y))
double code(double x, double y, double z, double t, double a, double b) {
return (x * exp((((y * log(z)) + ((t - 1.0) * log(a))) - b))) / y;
}
real(8) function code(x, y, z, t, a, b)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
code = (x * exp((((y * log(z)) + ((t - 1.0d0) * log(a))) - b))) / y
end function
public static double code(double x, double y, double z, double t, double a, double b) {
return (x * Math.exp((((y * Math.log(z)) + ((t - 1.0) * Math.log(a))) - b))) / y;
}
def code(x, y, z, t, a, b): return (x * math.exp((((y * math.log(z)) + ((t - 1.0) * math.log(a))) - b))) / y
function code(x, y, z, t, a, b) return Float64(Float64(x * exp(Float64(Float64(Float64(y * log(z)) + Float64(Float64(t - 1.0) * log(a))) - b))) / y) end
function tmp = code(x, y, z, t, a, b) tmp = (x * exp((((y * log(z)) + ((t - 1.0) * log(a))) - b))) / y; end
code[x_, y_, z_, t_, a_, b_] := N[(N[(x * N[Exp[N[(N[(N[(y * N[Log[z], $MachinePrecision]), $MachinePrecision] + N[(N[(t - 1.0), $MachinePrecision] * N[Log[a], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - b), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]
\begin{array}{l}
\\
\frac{x \cdot e^{\left(y \cdot \log z + \left(t - 1\right) \cdot \log a\right) - b}}{y}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 24 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x y z t a b) :precision binary64 (/ (* x (exp (- (+ (* y (log z)) (* (- t 1.0) (log a))) b))) y))
double code(double x, double y, double z, double t, double a, double b) {
return (x * exp((((y * log(z)) + ((t - 1.0) * log(a))) - b))) / y;
}
real(8) function code(x, y, z, t, a, b)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
code = (x * exp((((y * log(z)) + ((t - 1.0d0) * log(a))) - b))) / y
end function
public static double code(double x, double y, double z, double t, double a, double b) {
return (x * Math.exp((((y * Math.log(z)) + ((t - 1.0) * Math.log(a))) - b))) / y;
}
def code(x, y, z, t, a, b): return (x * math.exp((((y * math.log(z)) + ((t - 1.0) * math.log(a))) - b))) / y
function code(x, y, z, t, a, b) return Float64(Float64(x * exp(Float64(Float64(Float64(y * log(z)) + Float64(Float64(t - 1.0) * log(a))) - b))) / y) end
function tmp = code(x, y, z, t, a, b) tmp = (x * exp((((y * log(z)) + ((t - 1.0) * log(a))) - b))) / y; end
code[x_, y_, z_, t_, a_, b_] := N[(N[(x * N[Exp[N[(N[(N[(y * N[Log[z], $MachinePrecision]), $MachinePrecision] + N[(N[(t - 1.0), $MachinePrecision] * N[Log[a], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - b), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]
\begin{array}{l}
\\
\frac{x \cdot e^{\left(y \cdot \log z + \left(t - 1\right) \cdot \log a\right) - b}}{y}
\end{array}
(FPCore (x y z t a b) :precision binary64 (/ (* x (exp (- (+ (* y (log z)) (* (log a) (+ t -1.0))) b))) y))
double code(double x, double y, double z, double t, double a, double b) {
return (x * exp((((y * log(z)) + (log(a) * (t + -1.0))) - b))) / y;
}
real(8) function code(x, y, z, t, a, b)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
code = (x * exp((((y * log(z)) + (log(a) * (t + (-1.0d0)))) - b))) / y
end function
public static double code(double x, double y, double z, double t, double a, double b) {
return (x * Math.exp((((y * Math.log(z)) + (Math.log(a) * (t + -1.0))) - b))) / y;
}
def code(x, y, z, t, a, b): return (x * math.exp((((y * math.log(z)) + (math.log(a) * (t + -1.0))) - b))) / y
function code(x, y, z, t, a, b) return Float64(Float64(x * exp(Float64(Float64(Float64(y * log(z)) + Float64(log(a) * Float64(t + -1.0))) - b))) / y) end
function tmp = code(x, y, z, t, a, b) tmp = (x * exp((((y * log(z)) + (log(a) * (t + -1.0))) - b))) / y; end
code[x_, y_, z_, t_, a_, b_] := N[(N[(x * N[Exp[N[(N[(N[(y * N[Log[z], $MachinePrecision]), $MachinePrecision] + N[(N[Log[a], $MachinePrecision] * N[(t + -1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - b), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]
\begin{array}{l}
\\
\frac{x \cdot e^{\left(y \cdot \log z + \log a \cdot \left(t + -1\right)\right) - b}}{y}
\end{array}
Initial program 98.7%
Final simplification98.7%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (/ (* x (exp (- (+ (* y (log z)) (* (log a) (+ t -1.0))) b))) y)))
(if (<= t_1 (- INFINITY))
(* x (/ (fma b (fma 0.5 (/ b y) (/ -1.0 y)) (/ 1.0 y)) a))
(if (<= t_1 4e-108)
(/
x
(* y (fma b (fma b (fma a 0.5 (* 0.16666666666666666 (* a b))) a) a)))
(* x (/ (fma b (fma 0.5 (/ b a) (/ -1.0 a)) (/ 1.0 a)) y))))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (x * exp((((y * log(z)) + (log(a) * (t + -1.0))) - b))) / y;
double tmp;
if (t_1 <= -((double) INFINITY)) {
tmp = x * (fma(b, fma(0.5, (b / y), (-1.0 / y)), (1.0 / y)) / a);
} else if (t_1 <= 4e-108) {
tmp = x / (y * fma(b, fma(b, fma(a, 0.5, (0.16666666666666666 * (a * b))), a), a));
} else {
tmp = x * (fma(b, fma(0.5, (b / a), (-1.0 / a)), (1.0 / a)) / y);
}
return tmp;
}
function code(x, y, z, t, a, b) t_1 = Float64(Float64(x * exp(Float64(Float64(Float64(y * log(z)) + Float64(log(a) * Float64(t + -1.0))) - b))) / y) tmp = 0.0 if (t_1 <= Float64(-Inf)) tmp = Float64(x * Float64(fma(b, fma(0.5, Float64(b / y), Float64(-1.0 / y)), Float64(1.0 / y)) / a)); elseif (t_1 <= 4e-108) tmp = Float64(x / Float64(y * fma(b, fma(b, fma(a, 0.5, Float64(0.16666666666666666 * Float64(a * b))), a), a))); else tmp = Float64(x * Float64(fma(b, fma(0.5, Float64(b / a), Float64(-1.0 / a)), Float64(1.0 / a)) / y)); end return tmp end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(x * N[Exp[N[(N[(N[(y * N[Log[z], $MachinePrecision]), $MachinePrecision] + N[(N[Log[a], $MachinePrecision] * N[(t + -1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - b), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]}, If[LessEqual[t$95$1, (-Infinity)], N[(x * N[(N[(b * N[(0.5 * N[(b / y), $MachinePrecision] + N[(-1.0 / y), $MachinePrecision]), $MachinePrecision] + N[(1.0 / y), $MachinePrecision]), $MachinePrecision] / a), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$1, 4e-108], N[(x / N[(y * N[(b * N[(b * N[(a * 0.5 + N[(0.16666666666666666 * N[(a * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + a), $MachinePrecision] + a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x * N[(N[(b * N[(0.5 * N[(b / a), $MachinePrecision] + N[(-1.0 / a), $MachinePrecision]), $MachinePrecision] + N[(1.0 / a), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{x \cdot e^{\left(y \cdot \log z + \log a \cdot \left(t + -1\right)\right) - b}}{y}\\
\mathbf{if}\;t\_1 \leq -\infty:\\
\;\;\;\;x \cdot \frac{\mathsf{fma}\left(b, \mathsf{fma}\left(0.5, \frac{b}{y}, \frac{-1}{y}\right), \frac{1}{y}\right)}{a}\\
\mathbf{elif}\;t\_1 \leq 4 \cdot 10^{-108}:\\
\;\;\;\;\frac{x}{y \cdot \mathsf{fma}\left(b, \mathsf{fma}\left(b, \mathsf{fma}\left(a, 0.5, 0.16666666666666666 \cdot \left(a \cdot b\right)\right), a\right), a\right)}\\
\mathbf{else}:\\
\;\;\;\;x \cdot \frac{\mathsf{fma}\left(b, \mathsf{fma}\left(0.5, \frac{b}{a}, \frac{-1}{a}\right), \frac{1}{a}\right)}{y}\\
\end{array}
\end{array}
if (/.f64 (*.f64 x (exp.f64 (-.f64 (+.f64 (*.f64 y (log.f64 z)) (*.f64 (-.f64 t #s(literal 1 binary64)) (log.f64 a))) b))) y) < -inf.0Initial program 100.0%
Taylor expanded in t around 0
associate-/l*N/A
lower-*.f64N/A
exp-diffN/A
associate-/l/N/A
lower-/.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
exp-diffN/A
lower-/.f64N/A
*-commutativeN/A
exp-to-powN/A
lower-pow.f64N/A
rem-exp-logN/A
lower-*.f64N/A
lower-exp.f6472.7
Applied rewrites72.7%
Taylor expanded in y around 0
Applied rewrites55.9%
Taylor expanded in b around 0
Applied rewrites47.8%
Taylor expanded in a around 0
Applied rewrites53.5%
if -inf.0 < (/.f64 (*.f64 x (exp.f64 (-.f64 (+.f64 (*.f64 y (log.f64 z)) (*.f64 (-.f64 t #s(literal 1 binary64)) (log.f64 a))) b))) y) < 4.00000000000000016e-108Initial program 97.8%
Taylor expanded in y around 0
*-commutativeN/A
exp-diffN/A
associate-*l/N/A
associate-/l/N/A
lower-/.f64N/A
lower-*.f64N/A
exp-prodN/A
lower-pow.f64N/A
rem-exp-logN/A
sub-negN/A
metadata-evalN/A
lower-+.f64N/A
lower-*.f64N/A
lower-exp.f6474.7
Applied rewrites74.7%
Taylor expanded in t around 0
Applied rewrites66.7%
Taylor expanded in b around 0
Applied rewrites63.7%
if 4.00000000000000016e-108 < (/.f64 (*.f64 x (exp.f64 (-.f64 (+.f64 (*.f64 y (log.f64 z)) (*.f64 (-.f64 t #s(literal 1 binary64)) (log.f64 a))) b))) y) Initial program 99.0%
Taylor expanded in t around 0
associate-/l*N/A
lower-*.f64N/A
exp-diffN/A
associate-/l/N/A
lower-/.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
exp-diffN/A
lower-/.f64N/A
*-commutativeN/A
exp-to-powN/A
lower-pow.f64N/A
rem-exp-logN/A
lower-*.f64N/A
lower-exp.f6473.8
Applied rewrites73.8%
Taylor expanded in y around 0
Applied rewrites63.9%
Taylor expanded in b around 0
Applied rewrites48.4%
Taylor expanded in y around 0
Applied rewrites55.7%
Final simplification58.8%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (/ (* x (exp (- (+ (* y (log z)) (* (log a) (+ t -1.0))) b))) y))
(t_2 (* x (/ (fma b (fma 0.5 (/ b a) (/ -1.0 a)) (/ 1.0 a)) y))))
(if (<= t_1 (- INFINITY))
t_2
(if (<= t_1 4e-108)
(/
x
(* y (fma b (fma b (fma a 0.5 (* 0.16666666666666666 (* a b))) a) a)))
t_2))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (x * exp((((y * log(z)) + (log(a) * (t + -1.0))) - b))) / y;
double t_2 = x * (fma(b, fma(0.5, (b / a), (-1.0 / a)), (1.0 / a)) / y);
double tmp;
if (t_1 <= -((double) INFINITY)) {
tmp = t_2;
} else if (t_1 <= 4e-108) {
tmp = x / (y * fma(b, fma(b, fma(a, 0.5, (0.16666666666666666 * (a * b))), a), a));
} else {
tmp = t_2;
}
return tmp;
}
function code(x, y, z, t, a, b) t_1 = Float64(Float64(x * exp(Float64(Float64(Float64(y * log(z)) + Float64(log(a) * Float64(t + -1.0))) - b))) / y) t_2 = Float64(x * Float64(fma(b, fma(0.5, Float64(b / a), Float64(-1.0 / a)), Float64(1.0 / a)) / y)) tmp = 0.0 if (t_1 <= Float64(-Inf)) tmp = t_2; elseif (t_1 <= 4e-108) tmp = Float64(x / Float64(y * fma(b, fma(b, fma(a, 0.5, Float64(0.16666666666666666 * Float64(a * b))), a), a))); else tmp = t_2; end return tmp end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(x * N[Exp[N[(N[(N[(y * N[Log[z], $MachinePrecision]), $MachinePrecision] + N[(N[Log[a], $MachinePrecision] * N[(t + -1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - b), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]}, Block[{t$95$2 = N[(x * N[(N[(b * N[(0.5 * N[(b / a), $MachinePrecision] + N[(-1.0 / a), $MachinePrecision]), $MachinePrecision] + N[(1.0 / a), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, (-Infinity)], t$95$2, If[LessEqual[t$95$1, 4e-108], N[(x / N[(y * N[(b * N[(b * N[(a * 0.5 + N[(0.16666666666666666 * N[(a * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + a), $MachinePrecision] + a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$2]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{x \cdot e^{\left(y \cdot \log z + \log a \cdot \left(t + -1\right)\right) - b}}{y}\\
t_2 := x \cdot \frac{\mathsf{fma}\left(b, \mathsf{fma}\left(0.5, \frac{b}{a}, \frac{-1}{a}\right), \frac{1}{a}\right)}{y}\\
\mathbf{if}\;t\_1 \leq -\infty:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;t\_1 \leq 4 \cdot 10^{-108}:\\
\;\;\;\;\frac{x}{y \cdot \mathsf{fma}\left(b, \mathsf{fma}\left(b, \mathsf{fma}\left(a, 0.5, 0.16666666666666666 \cdot \left(a \cdot b\right)\right), a\right), a\right)}\\
\mathbf{else}:\\
\;\;\;\;t\_2\\
\end{array}
\end{array}
if (/.f64 (*.f64 x (exp.f64 (-.f64 (+.f64 (*.f64 y (log.f64 z)) (*.f64 (-.f64 t #s(literal 1 binary64)) (log.f64 a))) b))) y) < -inf.0 or 4.00000000000000016e-108 < (/.f64 (*.f64 x (exp.f64 (-.f64 (+.f64 (*.f64 y (log.f64 z)) (*.f64 (-.f64 t #s(literal 1 binary64)) (log.f64 a))) b))) y) Initial program 99.5%
Taylor expanded in t around 0
associate-/l*N/A
lower-*.f64N/A
exp-diffN/A
associate-/l/N/A
lower-/.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
exp-diffN/A
lower-/.f64N/A
*-commutativeN/A
exp-to-powN/A
lower-pow.f64N/A
rem-exp-logN/A
lower-*.f64N/A
lower-exp.f6473.2
Applied rewrites73.2%
Taylor expanded in y around 0
Applied rewrites59.9%
Taylor expanded in b around 0
Applied rewrites48.1%
Taylor expanded in y around 0
Applied rewrites56.7%
if -inf.0 < (/.f64 (*.f64 x (exp.f64 (-.f64 (+.f64 (*.f64 y (log.f64 z)) (*.f64 (-.f64 t #s(literal 1 binary64)) (log.f64 a))) b))) y) < 4.00000000000000016e-108Initial program 97.8%
Taylor expanded in y around 0
*-commutativeN/A
exp-diffN/A
associate-*l/N/A
associate-/l/N/A
lower-/.f64N/A
lower-*.f64N/A
exp-prodN/A
lower-pow.f64N/A
rem-exp-logN/A
sub-negN/A
metadata-evalN/A
lower-+.f64N/A
lower-*.f64N/A
lower-exp.f6474.7
Applied rewrites74.7%
Taylor expanded in t around 0
Applied rewrites66.7%
Taylor expanded in b around 0
Applied rewrites63.7%
Final simplification59.9%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (/ (* x (exp (- (+ (* y (log z)) (* (log a) (+ t -1.0))) b))) y))
(t_2 (* x (* b (* b (/ 0.5 (* y a)))))))
(if (<= t_1 (- INFINITY))
t_2
(if (<= t_1 1e+251)
(/
x
(* y (fma b (fma b (fma a 0.5 (* 0.16666666666666666 (* a b))) a) a)))
t_2))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (x * exp((((y * log(z)) + (log(a) * (t + -1.0))) - b))) / y;
double t_2 = x * (b * (b * (0.5 / (y * a))));
double tmp;
if (t_1 <= -((double) INFINITY)) {
tmp = t_2;
} else if (t_1 <= 1e+251) {
tmp = x / (y * fma(b, fma(b, fma(a, 0.5, (0.16666666666666666 * (a * b))), a), a));
} else {
tmp = t_2;
}
return tmp;
}
function code(x, y, z, t, a, b) t_1 = Float64(Float64(x * exp(Float64(Float64(Float64(y * log(z)) + Float64(log(a) * Float64(t + -1.0))) - b))) / y) t_2 = Float64(x * Float64(b * Float64(b * Float64(0.5 / Float64(y * a))))) tmp = 0.0 if (t_1 <= Float64(-Inf)) tmp = t_2; elseif (t_1 <= 1e+251) tmp = Float64(x / Float64(y * fma(b, fma(b, fma(a, 0.5, Float64(0.16666666666666666 * Float64(a * b))), a), a))); else tmp = t_2; end return tmp end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(x * N[Exp[N[(N[(N[(y * N[Log[z], $MachinePrecision]), $MachinePrecision] + N[(N[Log[a], $MachinePrecision] * N[(t + -1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - b), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]}, Block[{t$95$2 = N[(x * N[(b * N[(b * N[(0.5 / N[(y * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, (-Infinity)], t$95$2, If[LessEqual[t$95$1, 1e+251], N[(x / N[(y * N[(b * N[(b * N[(a * 0.5 + N[(0.16666666666666666 * N[(a * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + a), $MachinePrecision] + a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$2]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{x \cdot e^{\left(y \cdot \log z + \log a \cdot \left(t + -1\right)\right) - b}}{y}\\
t_2 := x \cdot \left(b \cdot \left(b \cdot \frac{0.5}{y \cdot a}\right)\right)\\
\mathbf{if}\;t\_1 \leq -\infty:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;t\_1 \leq 10^{+251}:\\
\;\;\;\;\frac{x}{y \cdot \mathsf{fma}\left(b, \mathsf{fma}\left(b, \mathsf{fma}\left(a, 0.5, 0.16666666666666666 \cdot \left(a \cdot b\right)\right), a\right), a\right)}\\
\mathbf{else}:\\
\;\;\;\;t\_2\\
\end{array}
\end{array}
if (/.f64 (*.f64 x (exp.f64 (-.f64 (+.f64 (*.f64 y (log.f64 z)) (*.f64 (-.f64 t #s(literal 1 binary64)) (log.f64 a))) b))) y) < -inf.0 or 1e251 < (/.f64 (*.f64 x (exp.f64 (-.f64 (+.f64 (*.f64 y (log.f64 z)) (*.f64 (-.f64 t #s(literal 1 binary64)) (log.f64 a))) b))) y) Initial program 100.0%
Taylor expanded in t around 0
associate-/l*N/A
lower-*.f64N/A
exp-diffN/A
associate-/l/N/A
lower-/.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
exp-diffN/A
lower-/.f64N/A
*-commutativeN/A
exp-to-powN/A
lower-pow.f64N/A
rem-exp-logN/A
lower-*.f64N/A
lower-exp.f6473.1
Applied rewrites73.1%
Taylor expanded in y around 0
Applied rewrites59.0%
Taylor expanded in b around 0
Applied rewrites46.5%
Taylor expanded in b around inf
Applied rewrites45.0%
if -inf.0 < (/.f64 (*.f64 x (exp.f64 (-.f64 (+.f64 (*.f64 y (log.f64 z)) (*.f64 (-.f64 t #s(literal 1 binary64)) (log.f64 a))) b))) y) < 1e251Initial program 97.4%
Taylor expanded in y around 0
*-commutativeN/A
exp-diffN/A
associate-*l/N/A
associate-/l/N/A
lower-/.f64N/A
lower-*.f64N/A
exp-prodN/A
lower-pow.f64N/A
rem-exp-logN/A
sub-negN/A
metadata-evalN/A
lower-+.f64N/A
lower-*.f64N/A
lower-exp.f6476.2
Applied rewrites76.2%
Taylor expanded in t around 0
Applied rewrites67.1%
Taylor expanded in b around 0
Applied rewrites64.3%
Final simplification54.6%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (* y (exp b))) (t_2 (/ (* x (exp (- (* y (log z)) b))) y)))
(if (<= y -1.1e+22)
t_2
(if (<= y -5.5e-183)
(/ (/ x t_1) a)
(if (<= y 4.4e-212)
(/ (* x (exp (- (* t (log a)) b))) y)
(if (<= y 4.8e-11) (/ x (* a t_1)) t_2))))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = y * exp(b);
double t_2 = (x * exp(((y * log(z)) - b))) / y;
double tmp;
if (y <= -1.1e+22) {
tmp = t_2;
} else if (y <= -5.5e-183) {
tmp = (x / t_1) / a;
} else if (y <= 4.4e-212) {
tmp = (x * exp(((t * log(a)) - b))) / y;
} else if (y <= 4.8e-11) {
tmp = x / (a * t_1);
} else {
tmp = t_2;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: t_1
real(8) :: t_2
real(8) :: tmp
t_1 = y * exp(b)
t_2 = (x * exp(((y * log(z)) - b))) / y
if (y <= (-1.1d+22)) then
tmp = t_2
else if (y <= (-5.5d-183)) then
tmp = (x / t_1) / a
else if (y <= 4.4d-212) then
tmp = (x * exp(((t * log(a)) - b))) / y
else if (y <= 4.8d-11) then
tmp = x / (a * t_1)
else
tmp = t_2
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double t_1 = y * Math.exp(b);
double t_2 = (x * Math.exp(((y * Math.log(z)) - b))) / y;
double tmp;
if (y <= -1.1e+22) {
tmp = t_2;
} else if (y <= -5.5e-183) {
tmp = (x / t_1) / a;
} else if (y <= 4.4e-212) {
tmp = (x * Math.exp(((t * Math.log(a)) - b))) / y;
} else if (y <= 4.8e-11) {
tmp = x / (a * t_1);
} else {
tmp = t_2;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = y * math.exp(b) t_2 = (x * math.exp(((y * math.log(z)) - b))) / y tmp = 0 if y <= -1.1e+22: tmp = t_2 elif y <= -5.5e-183: tmp = (x / t_1) / a elif y <= 4.4e-212: tmp = (x * math.exp(((t * math.log(a)) - b))) / y elif y <= 4.8e-11: tmp = x / (a * t_1) else: tmp = t_2 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(y * exp(b)) t_2 = Float64(Float64(x * exp(Float64(Float64(y * log(z)) - b))) / y) tmp = 0.0 if (y <= -1.1e+22) tmp = t_2; elseif (y <= -5.5e-183) tmp = Float64(Float64(x / t_1) / a); elseif (y <= 4.4e-212) tmp = Float64(Float64(x * exp(Float64(Float64(t * log(a)) - b))) / y); elseif (y <= 4.8e-11) tmp = Float64(x / Float64(a * t_1)); else tmp = t_2; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = y * exp(b); t_2 = (x * exp(((y * log(z)) - b))) / y; tmp = 0.0; if (y <= -1.1e+22) tmp = t_2; elseif (y <= -5.5e-183) tmp = (x / t_1) / a; elseif (y <= 4.4e-212) tmp = (x * exp(((t * log(a)) - b))) / y; elseif (y <= 4.8e-11) tmp = x / (a * t_1); else tmp = t_2; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(y * N[Exp[b], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(N[(x * N[Exp[N[(N[(y * N[Log[z], $MachinePrecision]), $MachinePrecision] - b), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]}, If[LessEqual[y, -1.1e+22], t$95$2, If[LessEqual[y, -5.5e-183], N[(N[(x / t$95$1), $MachinePrecision] / a), $MachinePrecision], If[LessEqual[y, 4.4e-212], N[(N[(x * N[Exp[N[(N[(t * N[Log[a], $MachinePrecision]), $MachinePrecision] - b), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision], If[LessEqual[y, 4.8e-11], N[(x / N[(a * t$95$1), $MachinePrecision]), $MachinePrecision], t$95$2]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := y \cdot e^{b}\\
t_2 := \frac{x \cdot e^{y \cdot \log z - b}}{y}\\
\mathbf{if}\;y \leq -1.1 \cdot 10^{+22}:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;y \leq -5.5 \cdot 10^{-183}:\\
\;\;\;\;\frac{\frac{x}{t\_1}}{a}\\
\mathbf{elif}\;y \leq 4.4 \cdot 10^{-212}:\\
\;\;\;\;\frac{x \cdot e^{t \cdot \log a - b}}{y}\\
\mathbf{elif}\;y \leq 4.8 \cdot 10^{-11}:\\
\;\;\;\;\frac{x}{a \cdot t\_1}\\
\mathbf{else}:\\
\;\;\;\;t\_2\\
\end{array}
\end{array}
if y < -1.1e22 or 4.8000000000000002e-11 < y Initial program 100.0%
Taylor expanded in y around inf
lower-*.f64N/A
lower-log.f6493.8
Applied rewrites93.8%
if -1.1e22 < y < -5.4999999999999999e-183Initial program 98.2%
Taylor expanded in y around 0
*-commutativeN/A
exp-diffN/A
associate-*l/N/A
associate-/l/N/A
lower-/.f64N/A
lower-*.f64N/A
exp-prodN/A
lower-pow.f64N/A
rem-exp-logN/A
sub-negN/A
metadata-evalN/A
lower-+.f64N/A
lower-*.f64N/A
lower-exp.f6491.1
Applied rewrites91.1%
Taylor expanded in t around 0
Applied rewrites84.9%
Applied rewrites87.3%
if -5.4999999999999999e-183 < y < 4.40000000000000006e-212Initial program 98.9%
Taylor expanded in t around inf
*-commutativeN/A
lower-*.f64N/A
rem-exp-logN/A
lower-log.f64N/A
rem-exp-log90.0
Applied rewrites90.0%
if 4.40000000000000006e-212 < y < 4.8000000000000002e-11Initial program 94.8%
Taylor expanded in y around 0
*-commutativeN/A
exp-diffN/A
associate-*l/N/A
associate-/l/N/A
lower-/.f64N/A
lower-*.f64N/A
exp-prodN/A
lower-pow.f64N/A
rem-exp-logN/A
sub-negN/A
metadata-evalN/A
lower-+.f64N/A
lower-*.f64N/A
lower-exp.f6488.1
Applied rewrites88.1%
Taylor expanded in t around 0
Applied rewrites81.9%
Applied rewrites81.9%
Final simplification90.3%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (/ (* x (exp (- (* y (log z)) b))) y)))
(if (<= y -6e+26)
t_1
(if (<= y 1.32e+63) (/ (* x (exp (- (* (log a) (+ t -1.0)) b))) y) t_1))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (x * exp(((y * log(z)) - b))) / y;
double tmp;
if (y <= -6e+26) {
tmp = t_1;
} else if (y <= 1.32e+63) {
tmp = (x * exp(((log(a) * (t + -1.0)) - b))) / y;
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: t_1
real(8) :: tmp
t_1 = (x * exp(((y * log(z)) - b))) / y
if (y <= (-6d+26)) then
tmp = t_1
else if (y <= 1.32d+63) then
tmp = (x * exp(((log(a) * (t + (-1.0d0))) - b))) / y
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (x * Math.exp(((y * Math.log(z)) - b))) / y;
double tmp;
if (y <= -6e+26) {
tmp = t_1;
} else if (y <= 1.32e+63) {
tmp = (x * Math.exp(((Math.log(a) * (t + -1.0)) - b))) / y;
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = (x * math.exp(((y * math.log(z)) - b))) / y tmp = 0 if y <= -6e+26: tmp = t_1 elif y <= 1.32e+63: tmp = (x * math.exp(((math.log(a) * (t + -1.0)) - b))) / y else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(Float64(x * exp(Float64(Float64(y * log(z)) - b))) / y) tmp = 0.0 if (y <= -6e+26) tmp = t_1; elseif (y <= 1.32e+63) tmp = Float64(Float64(x * exp(Float64(Float64(log(a) * Float64(t + -1.0)) - b))) / y); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = (x * exp(((y * log(z)) - b))) / y; tmp = 0.0; if (y <= -6e+26) tmp = t_1; elseif (y <= 1.32e+63) tmp = (x * exp(((log(a) * (t + -1.0)) - b))) / y; else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(x * N[Exp[N[(N[(y * N[Log[z], $MachinePrecision]), $MachinePrecision] - b), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]}, If[LessEqual[y, -6e+26], t$95$1, If[LessEqual[y, 1.32e+63], N[(N[(x * N[Exp[N[(N[(N[Log[a], $MachinePrecision] * N[(t + -1.0), $MachinePrecision]), $MachinePrecision] - b), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{x \cdot e^{y \cdot \log z - b}}{y}\\
\mathbf{if}\;y \leq -6 \cdot 10^{+26}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 1.32 \cdot 10^{+63}:\\
\;\;\;\;\frac{x \cdot e^{\log a \cdot \left(t + -1\right) - b}}{y}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -5.99999999999999994e26 or 1.32e63 < y Initial program 100.0%
Taylor expanded in y around inf
lower-*.f64N/A
lower-log.f6495.4
Applied rewrites95.4%
if -5.99999999999999994e26 < y < 1.32e63Initial program 97.8%
Taylor expanded in y around 0
lower-*.f64N/A
rem-exp-logN/A
lower-log.f64N/A
rem-exp-logN/A
sub-negN/A
metadata-evalN/A
lower-+.f6496.2
Applied rewrites96.2%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (/ (* x (exp (- (* y (log z)) b))) y)))
(if (<= y -3.5e-5)
t_1
(if (<= y 4.8e-11) (/ (* x (pow a (+ t -1.0))) (* y (exp b))) t_1))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (x * exp(((y * log(z)) - b))) / y;
double tmp;
if (y <= -3.5e-5) {
tmp = t_1;
} else if (y <= 4.8e-11) {
tmp = (x * pow(a, (t + -1.0))) / (y * exp(b));
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: t_1
real(8) :: tmp
t_1 = (x * exp(((y * log(z)) - b))) / y
if (y <= (-3.5d-5)) then
tmp = t_1
else if (y <= 4.8d-11) then
tmp = (x * (a ** (t + (-1.0d0)))) / (y * exp(b))
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (x * Math.exp(((y * Math.log(z)) - b))) / y;
double tmp;
if (y <= -3.5e-5) {
tmp = t_1;
} else if (y <= 4.8e-11) {
tmp = (x * Math.pow(a, (t + -1.0))) / (y * Math.exp(b));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = (x * math.exp(((y * math.log(z)) - b))) / y tmp = 0 if y <= -3.5e-5: tmp = t_1 elif y <= 4.8e-11: tmp = (x * math.pow(a, (t + -1.0))) / (y * math.exp(b)) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(Float64(x * exp(Float64(Float64(y * log(z)) - b))) / y) tmp = 0.0 if (y <= -3.5e-5) tmp = t_1; elseif (y <= 4.8e-11) tmp = Float64(Float64(x * (a ^ Float64(t + -1.0))) / Float64(y * exp(b))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = (x * exp(((y * log(z)) - b))) / y; tmp = 0.0; if (y <= -3.5e-5) tmp = t_1; elseif (y <= 4.8e-11) tmp = (x * (a ^ (t + -1.0))) / (y * exp(b)); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(x * N[Exp[N[(N[(y * N[Log[z], $MachinePrecision]), $MachinePrecision] - b), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]}, If[LessEqual[y, -3.5e-5], t$95$1, If[LessEqual[y, 4.8e-11], N[(N[(x * N[Power[a, N[(t + -1.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / N[(y * N[Exp[b], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{x \cdot e^{y \cdot \log z - b}}{y}\\
\mathbf{if}\;y \leq -3.5 \cdot 10^{-5}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 4.8 \cdot 10^{-11}:\\
\;\;\;\;\frac{x \cdot {a}^{\left(t + -1\right)}}{y \cdot e^{b}}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -3.4999999999999997e-5 or 4.8000000000000002e-11 < y Initial program 100.0%
Taylor expanded in y around inf
lower-*.f64N/A
lower-log.f6493.3
Applied rewrites93.3%
if -3.4999999999999997e-5 < y < 4.8000000000000002e-11Initial program 97.3%
Taylor expanded in y around 0
*-commutativeN/A
exp-diffN/A
associate-*l/N/A
associate-/l/N/A
lower-/.f64N/A
lower-*.f64N/A
exp-prodN/A
lower-pow.f64N/A
rem-exp-logN/A
sub-negN/A
metadata-evalN/A
lower-+.f64N/A
lower-*.f64N/A
lower-exp.f6490.7
Applied rewrites90.7%
Final simplification92.1%
(FPCore (x y z t a b) :precision binary64 (let* ((t_1 (/ (* x (exp (- (* y (log z)) b))) y))) (if (<= y -1.1e+22) t_1 (if (<= y 4.8e-11) (/ (/ x (* y (exp b))) a) t_1))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (x * exp(((y * log(z)) - b))) / y;
double tmp;
if (y <= -1.1e+22) {
tmp = t_1;
} else if (y <= 4.8e-11) {
tmp = (x / (y * exp(b))) / a;
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: t_1
real(8) :: tmp
t_1 = (x * exp(((y * log(z)) - b))) / y
if (y <= (-1.1d+22)) then
tmp = t_1
else if (y <= 4.8d-11) then
tmp = (x / (y * exp(b))) / a
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (x * Math.exp(((y * Math.log(z)) - b))) / y;
double tmp;
if (y <= -1.1e+22) {
tmp = t_1;
} else if (y <= 4.8e-11) {
tmp = (x / (y * Math.exp(b))) / a;
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = (x * math.exp(((y * math.log(z)) - b))) / y tmp = 0 if y <= -1.1e+22: tmp = t_1 elif y <= 4.8e-11: tmp = (x / (y * math.exp(b))) / a else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(Float64(x * exp(Float64(Float64(y * log(z)) - b))) / y) tmp = 0.0 if (y <= -1.1e+22) tmp = t_1; elseif (y <= 4.8e-11) tmp = Float64(Float64(x / Float64(y * exp(b))) / a); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = (x * exp(((y * log(z)) - b))) / y; tmp = 0.0; if (y <= -1.1e+22) tmp = t_1; elseif (y <= 4.8e-11) tmp = (x / (y * exp(b))) / a; else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(x * N[Exp[N[(N[(y * N[Log[z], $MachinePrecision]), $MachinePrecision] - b), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]}, If[LessEqual[y, -1.1e+22], t$95$1, If[LessEqual[y, 4.8e-11], N[(N[(x / N[(y * N[Exp[b], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / a), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{x \cdot e^{y \cdot \log z - b}}{y}\\
\mathbf{if}\;y \leq -1.1 \cdot 10^{+22}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 4.8 \cdot 10^{-11}:\\
\;\;\;\;\frac{\frac{x}{y \cdot e^{b}}}{a}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -1.1e22 or 4.8000000000000002e-11 < y Initial program 100.0%
Taylor expanded in y around inf
lower-*.f64N/A
lower-log.f6493.8
Applied rewrites93.8%
if -1.1e22 < y < 4.8000000000000002e-11Initial program 97.4%
Taylor expanded in y around 0
*-commutativeN/A
exp-diffN/A
associate-*l/N/A
associate-/l/N/A
lower-/.f64N/A
lower-*.f64N/A
exp-prodN/A
lower-pow.f64N/A
rem-exp-logN/A
sub-negN/A
metadata-evalN/A
lower-+.f64N/A
lower-*.f64N/A
lower-exp.f6489.5
Applied rewrites89.5%
Taylor expanded in t around 0
Applied rewrites77.4%
Applied rewrites79.0%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (* x (/ (/ (pow z y) y) a))))
(if (<= y -8.2e+26)
t_1
(if (<= y 1.65e+33) (/ (/ x (* y (exp b))) a) t_1))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = x * ((pow(z, y) / y) / a);
double tmp;
if (y <= -8.2e+26) {
tmp = t_1;
} else if (y <= 1.65e+33) {
tmp = (x / (y * exp(b))) / a;
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: t_1
real(8) :: tmp
t_1 = x * (((z ** y) / y) / a)
if (y <= (-8.2d+26)) then
tmp = t_1
else if (y <= 1.65d+33) then
tmp = (x / (y * exp(b))) / a
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double t_1 = x * ((Math.pow(z, y) / y) / a);
double tmp;
if (y <= -8.2e+26) {
tmp = t_1;
} else if (y <= 1.65e+33) {
tmp = (x / (y * Math.exp(b))) / a;
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = x * ((math.pow(z, y) / y) / a) tmp = 0 if y <= -8.2e+26: tmp = t_1 elif y <= 1.65e+33: tmp = (x / (y * math.exp(b))) / a else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(x * Float64(Float64((z ^ y) / y) / a)) tmp = 0.0 if (y <= -8.2e+26) tmp = t_1; elseif (y <= 1.65e+33) tmp = Float64(Float64(x / Float64(y * exp(b))) / a); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = x * (((z ^ y) / y) / a); tmp = 0.0; if (y <= -8.2e+26) tmp = t_1; elseif (y <= 1.65e+33) tmp = (x / (y * exp(b))) / a; else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(x * N[(N[(N[Power[z, y], $MachinePrecision] / y), $MachinePrecision] / a), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -8.2e+26], t$95$1, If[LessEqual[y, 1.65e+33], N[(N[(x / N[(y * N[Exp[b], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / a), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x \cdot \frac{\frac{{z}^{y}}{y}}{a}\\
\mathbf{if}\;y \leq -8.2 \cdot 10^{+26}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 1.65 \cdot 10^{+33}:\\
\;\;\;\;\frac{\frac{x}{y \cdot e^{b}}}{a}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -8.19999999999999967e26 or 1.64999999999999988e33 < y Initial program 100.0%
Taylor expanded in t around 0
associate-/l*N/A
lower-*.f64N/A
exp-diffN/A
associate-/l/N/A
lower-/.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
exp-diffN/A
lower-/.f64N/A
*-commutativeN/A
exp-to-powN/A
lower-pow.f64N/A
rem-exp-logN/A
lower-*.f64N/A
lower-exp.f6466.7
Applied rewrites66.7%
Applied rewrites66.7%
Taylor expanded in b around 0
Applied rewrites84.0%
if -8.19999999999999967e26 < y < 1.64999999999999988e33Initial program 97.6%
Taylor expanded in y around 0
*-commutativeN/A
exp-diffN/A
associate-*l/N/A
associate-/l/N/A
lower-/.f64N/A
lower-*.f64N/A
exp-prodN/A
lower-pow.f64N/A
rem-exp-logN/A
sub-negN/A
metadata-evalN/A
lower-+.f64N/A
lower-*.f64N/A
lower-exp.f6489.0
Applied rewrites89.0%
Taylor expanded in t around 0
Applied rewrites77.1%
Applied rewrites78.6%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (* y (exp b))))
(if (<= b -4.6e-67)
(/ (/ x t_1) a)
(if (<= b 9e-19) (/ (* x (pow a (+ t -1.0))) y) (/ x (* a t_1))))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = y * exp(b);
double tmp;
if (b <= -4.6e-67) {
tmp = (x / t_1) / a;
} else if (b <= 9e-19) {
tmp = (x * pow(a, (t + -1.0))) / y;
} else {
tmp = x / (a * t_1);
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: t_1
real(8) :: tmp
t_1 = y * exp(b)
if (b <= (-4.6d-67)) then
tmp = (x / t_1) / a
else if (b <= 9d-19) then
tmp = (x * (a ** (t + (-1.0d0)))) / y
else
tmp = x / (a * t_1)
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double t_1 = y * Math.exp(b);
double tmp;
if (b <= -4.6e-67) {
tmp = (x / t_1) / a;
} else if (b <= 9e-19) {
tmp = (x * Math.pow(a, (t + -1.0))) / y;
} else {
tmp = x / (a * t_1);
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = y * math.exp(b) tmp = 0 if b <= -4.6e-67: tmp = (x / t_1) / a elif b <= 9e-19: tmp = (x * math.pow(a, (t + -1.0))) / y else: tmp = x / (a * t_1) return tmp
function code(x, y, z, t, a, b) t_1 = Float64(y * exp(b)) tmp = 0.0 if (b <= -4.6e-67) tmp = Float64(Float64(x / t_1) / a); elseif (b <= 9e-19) tmp = Float64(Float64(x * (a ^ Float64(t + -1.0))) / y); else tmp = Float64(x / Float64(a * t_1)); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = y * exp(b); tmp = 0.0; if (b <= -4.6e-67) tmp = (x / t_1) / a; elseif (b <= 9e-19) tmp = (x * (a ^ (t + -1.0))) / y; else tmp = x / (a * t_1); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(y * N[Exp[b], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[b, -4.6e-67], N[(N[(x / t$95$1), $MachinePrecision] / a), $MachinePrecision], If[LessEqual[b, 9e-19], N[(N[(x * N[Power[a, N[(t + -1.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision], N[(x / N[(a * t$95$1), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := y \cdot e^{b}\\
\mathbf{if}\;b \leq -4.6 \cdot 10^{-67}:\\
\;\;\;\;\frac{\frac{x}{t\_1}}{a}\\
\mathbf{elif}\;b \leq 9 \cdot 10^{-19}:\\
\;\;\;\;\frac{x \cdot {a}^{\left(t + -1\right)}}{y}\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{a \cdot t\_1}\\
\end{array}
\end{array}
if b < -4.6000000000000001e-67Initial program 98.6%
Taylor expanded in y around 0
*-commutativeN/A
exp-diffN/A
associate-*l/N/A
associate-/l/N/A
lower-/.f64N/A
lower-*.f64N/A
exp-prodN/A
lower-pow.f64N/A
rem-exp-logN/A
sub-negN/A
metadata-evalN/A
lower-+.f64N/A
lower-*.f64N/A
lower-exp.f6473.3
Applied rewrites73.3%
Taylor expanded in t around 0
Applied rewrites84.2%
Applied rewrites87.7%
if -4.6000000000000001e-67 < b < 9.00000000000000026e-19Initial program 98.0%
Taylor expanded in y around 0
*-commutativeN/A
exp-diffN/A
associate-*l/N/A
associate-/l/N/A
lower-/.f64N/A
lower-*.f64N/A
exp-prodN/A
lower-pow.f64N/A
rem-exp-logN/A
sub-negN/A
metadata-evalN/A
lower-+.f64N/A
lower-*.f64N/A
lower-exp.f6469.7
Applied rewrites69.7%
Taylor expanded in t around 0
Applied rewrites37.9%
Taylor expanded in b around 0
Applied rewrites69.7%
if 9.00000000000000026e-19 < b Initial program 99.8%
Taylor expanded in y around 0
*-commutativeN/A
exp-diffN/A
associate-*l/N/A
associate-/l/N/A
lower-/.f64N/A
lower-*.f64N/A
exp-prodN/A
lower-pow.f64N/A
rem-exp-logN/A
sub-negN/A
metadata-evalN/A
lower-+.f64N/A
lower-*.f64N/A
lower-exp.f6460.2
Applied rewrites60.2%
Taylor expanded in t around 0
Applied rewrites75.4%
Applied rewrites75.4%
Final simplification77.0%
(FPCore (x y z t a b) :precision binary64 (if (<= b -1.5e-15) (/ x (* y (* a (exp b)))) (if (<= b 9e-19) (/ (* x (pow a (+ t -1.0))) y) (/ x (* a (* y (exp b)))))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (b <= -1.5e-15) {
tmp = x / (y * (a * exp(b)));
} else if (b <= 9e-19) {
tmp = (x * pow(a, (t + -1.0))) / y;
} else {
tmp = x / (a * (y * exp(b)));
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: tmp
if (b <= (-1.5d-15)) then
tmp = x / (y * (a * exp(b)))
else if (b <= 9d-19) then
tmp = (x * (a ** (t + (-1.0d0)))) / y
else
tmp = x / (a * (y * exp(b)))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (b <= -1.5e-15) {
tmp = x / (y * (a * Math.exp(b)));
} else if (b <= 9e-19) {
tmp = (x * Math.pow(a, (t + -1.0))) / y;
} else {
tmp = x / (a * (y * Math.exp(b)));
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if b <= -1.5e-15: tmp = x / (y * (a * math.exp(b))) elif b <= 9e-19: tmp = (x * math.pow(a, (t + -1.0))) / y else: tmp = x / (a * (y * math.exp(b))) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (b <= -1.5e-15) tmp = Float64(x / Float64(y * Float64(a * exp(b)))); elseif (b <= 9e-19) tmp = Float64(Float64(x * (a ^ Float64(t + -1.0))) / y); else tmp = Float64(x / Float64(a * Float64(y * exp(b)))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (b <= -1.5e-15) tmp = x / (y * (a * exp(b))); elseif (b <= 9e-19) tmp = (x * (a ^ (t + -1.0))) / y; else tmp = x / (a * (y * exp(b))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[b, -1.5e-15], N[(x / N[(y * N[(a * N[Exp[b], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 9e-19], N[(N[(x * N[Power[a, N[(t + -1.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision], N[(x / N[(a * N[(y * N[Exp[b], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq -1.5 \cdot 10^{-15}:\\
\;\;\;\;\frac{x}{y \cdot \left(a \cdot e^{b}\right)}\\
\mathbf{elif}\;b \leq 9 \cdot 10^{-19}:\\
\;\;\;\;\frac{x \cdot {a}^{\left(t + -1\right)}}{y}\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{a \cdot \left(y \cdot e^{b}\right)}\\
\end{array}
\end{array}
if b < -1.5e-15Initial program 99.7%
Taylor expanded in y around 0
*-commutativeN/A
exp-diffN/A
associate-*l/N/A
associate-/l/N/A
lower-/.f64N/A
lower-*.f64N/A
exp-prodN/A
lower-pow.f64N/A
rem-exp-logN/A
sub-negN/A
metadata-evalN/A
lower-+.f64N/A
lower-*.f64N/A
lower-exp.f6472.1
Applied rewrites72.1%
Taylor expanded in t around 0
Applied rewrites88.4%
if -1.5e-15 < b < 9.00000000000000026e-19Initial program 97.4%
Taylor expanded in y around 0
*-commutativeN/A
exp-diffN/A
associate-*l/N/A
associate-/l/N/A
lower-/.f64N/A
lower-*.f64N/A
exp-prodN/A
lower-pow.f64N/A
rem-exp-logN/A
sub-negN/A
metadata-evalN/A
lower-+.f64N/A
lower-*.f64N/A
lower-exp.f6470.8
Applied rewrites70.8%
Taylor expanded in t around 0
Applied rewrites40.2%
Taylor expanded in b around 0
Applied rewrites70.8%
if 9.00000000000000026e-19 < b Initial program 99.8%
Taylor expanded in y around 0
*-commutativeN/A
exp-diffN/A
associate-*l/N/A
associate-/l/N/A
lower-/.f64N/A
lower-*.f64N/A
exp-prodN/A
lower-pow.f64N/A
rem-exp-logN/A
sub-negN/A
metadata-evalN/A
lower-+.f64N/A
lower-*.f64N/A
lower-exp.f6460.2
Applied rewrites60.2%
Taylor expanded in t around 0
Applied rewrites75.4%
Applied rewrites75.4%
Final simplification76.8%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (/ x (* y (* a (exp b))))))
(if (<= b -1.5e-15)
t_1
(if (<= b 9e-19) (/ (* x (pow a (+ t -1.0))) y) t_1))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = x / (y * (a * exp(b)));
double tmp;
if (b <= -1.5e-15) {
tmp = t_1;
} else if (b <= 9e-19) {
tmp = (x * pow(a, (t + -1.0))) / y;
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: t_1
real(8) :: tmp
t_1 = x / (y * (a * exp(b)))
if (b <= (-1.5d-15)) then
tmp = t_1
else if (b <= 9d-19) then
tmp = (x * (a ** (t + (-1.0d0)))) / y
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double t_1 = x / (y * (a * Math.exp(b)));
double tmp;
if (b <= -1.5e-15) {
tmp = t_1;
} else if (b <= 9e-19) {
tmp = (x * Math.pow(a, (t + -1.0))) / y;
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = x / (y * (a * math.exp(b))) tmp = 0 if b <= -1.5e-15: tmp = t_1 elif b <= 9e-19: tmp = (x * math.pow(a, (t + -1.0))) / y else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(x / Float64(y * Float64(a * exp(b)))) tmp = 0.0 if (b <= -1.5e-15) tmp = t_1; elseif (b <= 9e-19) tmp = Float64(Float64(x * (a ^ Float64(t + -1.0))) / y); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = x / (y * (a * exp(b))); tmp = 0.0; if (b <= -1.5e-15) tmp = t_1; elseif (b <= 9e-19) tmp = (x * (a ^ (t + -1.0))) / y; else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(x / N[(y * N[(a * N[Exp[b], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[b, -1.5e-15], t$95$1, If[LessEqual[b, 9e-19], N[(N[(x * N[Power[a, N[(t + -1.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{x}{y \cdot \left(a \cdot e^{b}\right)}\\
\mathbf{if}\;b \leq -1.5 \cdot 10^{-15}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;b \leq 9 \cdot 10^{-19}:\\
\;\;\;\;\frac{x \cdot {a}^{\left(t + -1\right)}}{y}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if b < -1.5e-15 or 9.00000000000000026e-19 < b Initial program 99.7%
Taylor expanded in y around 0
*-commutativeN/A
exp-diffN/A
associate-*l/N/A
associate-/l/N/A
lower-/.f64N/A
lower-*.f64N/A
exp-prodN/A
lower-pow.f64N/A
rem-exp-logN/A
sub-negN/A
metadata-evalN/A
lower-+.f64N/A
lower-*.f64N/A
lower-exp.f6465.9
Applied rewrites65.9%
Taylor expanded in t around 0
Applied rewrites81.7%
if -1.5e-15 < b < 9.00000000000000026e-19Initial program 97.4%
Taylor expanded in y around 0
*-commutativeN/A
exp-diffN/A
associate-*l/N/A
associate-/l/N/A
lower-/.f64N/A
lower-*.f64N/A
exp-prodN/A
lower-pow.f64N/A
rem-exp-logN/A
sub-negN/A
metadata-evalN/A
lower-+.f64N/A
lower-*.f64N/A
lower-exp.f6470.8
Applied rewrites70.8%
Taylor expanded in t around 0
Applied rewrites40.2%
Taylor expanded in b around 0
Applied rewrites70.8%
Final simplification76.8%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (exp (- b))))
(if (<= b -5800000000.0)
(* x (/ t_1 y))
(if (<= b 0.78) (/ (* x (pow a (+ t -1.0))) y) (/ (* x t_1) y)))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = exp(-b);
double tmp;
if (b <= -5800000000.0) {
tmp = x * (t_1 / y);
} else if (b <= 0.78) {
tmp = (x * pow(a, (t + -1.0))) / y;
} else {
tmp = (x * t_1) / y;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: t_1
real(8) :: tmp
t_1 = exp(-b)
if (b <= (-5800000000.0d0)) then
tmp = x * (t_1 / y)
else if (b <= 0.78d0) then
tmp = (x * (a ** (t + (-1.0d0)))) / y
else
tmp = (x * t_1) / y
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double t_1 = Math.exp(-b);
double tmp;
if (b <= -5800000000.0) {
tmp = x * (t_1 / y);
} else if (b <= 0.78) {
tmp = (x * Math.pow(a, (t + -1.0))) / y;
} else {
tmp = (x * t_1) / y;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = math.exp(-b) tmp = 0 if b <= -5800000000.0: tmp = x * (t_1 / y) elif b <= 0.78: tmp = (x * math.pow(a, (t + -1.0))) / y else: tmp = (x * t_1) / y return tmp
function code(x, y, z, t, a, b) t_1 = exp(Float64(-b)) tmp = 0.0 if (b <= -5800000000.0) tmp = Float64(x * Float64(t_1 / y)); elseif (b <= 0.78) tmp = Float64(Float64(x * (a ^ Float64(t + -1.0))) / y); else tmp = Float64(Float64(x * t_1) / y); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = exp(-b); tmp = 0.0; if (b <= -5800000000.0) tmp = x * (t_1 / y); elseif (b <= 0.78) tmp = (x * (a ^ (t + -1.0))) / y; else tmp = (x * t_1) / y; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[Exp[(-b)], $MachinePrecision]}, If[LessEqual[b, -5800000000.0], N[(x * N[(t$95$1 / y), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 0.78], N[(N[(x * N[Power[a, N[(t + -1.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision], N[(N[(x * t$95$1), $MachinePrecision] / y), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := e^{-b}\\
\mathbf{if}\;b \leq -5800000000:\\
\;\;\;\;x \cdot \frac{t\_1}{y}\\
\mathbf{elif}\;b \leq 0.78:\\
\;\;\;\;\frac{x \cdot {a}^{\left(t + -1\right)}}{y}\\
\mathbf{else}:\\
\;\;\;\;\frac{x \cdot t\_1}{y}\\
\end{array}
\end{array}
if b < -5.8e9Initial program 100.0%
Taylor expanded in y around 0
lower-*.f64N/A
rem-exp-logN/A
lower-log.f64N/A
rem-exp-logN/A
sub-negN/A
metadata-evalN/A
lower-+.f6493.8
Applied rewrites93.8%
Taylor expanded in b around inf
neg-mul-1N/A
lower-neg.f6487.7
Applied rewrites87.7%
lift-/.f64N/A
lift-*.f64N/A
associate-/l*N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f6487.7
Applied rewrites87.7%
if -5.8e9 < b < 0.78000000000000003Initial program 97.3%
Taylor expanded in y around 0
*-commutativeN/A
exp-diffN/A
associate-*l/N/A
associate-/l/N/A
lower-/.f64N/A
lower-*.f64N/A
exp-prodN/A
lower-pow.f64N/A
rem-exp-logN/A
sub-negN/A
metadata-evalN/A
lower-+.f64N/A
lower-*.f64N/A
lower-exp.f6472.5
Applied rewrites72.5%
Taylor expanded in t around 0
Applied rewrites43.9%
Taylor expanded in b around 0
Applied rewrites72.2%
if 0.78000000000000003 < b Initial program 100.0%
Taylor expanded in y around 0
lower-*.f64N/A
rem-exp-logN/A
lower-log.f64N/A
rem-exp-logN/A
sub-negN/A
metadata-evalN/A
lower-+.f6480.0
Applied rewrites80.0%
Taylor expanded in b around inf
neg-mul-1N/A
lower-neg.f6474.4
Applied rewrites74.4%
Final simplification76.6%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (exp (- b))))
(if (<= b -160000000.0)
(* x (/ t_1 y))
(if (<= b 0.78) (* (pow a (+ t -1.0)) (/ x y)) (/ (* x t_1) y)))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = exp(-b);
double tmp;
if (b <= -160000000.0) {
tmp = x * (t_1 / y);
} else if (b <= 0.78) {
tmp = pow(a, (t + -1.0)) * (x / y);
} else {
tmp = (x * t_1) / y;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: t_1
real(8) :: tmp
t_1 = exp(-b)
if (b <= (-160000000.0d0)) then
tmp = x * (t_1 / y)
else if (b <= 0.78d0) then
tmp = (a ** (t + (-1.0d0))) * (x / y)
else
tmp = (x * t_1) / y
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double t_1 = Math.exp(-b);
double tmp;
if (b <= -160000000.0) {
tmp = x * (t_1 / y);
} else if (b <= 0.78) {
tmp = Math.pow(a, (t + -1.0)) * (x / y);
} else {
tmp = (x * t_1) / y;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = math.exp(-b) tmp = 0 if b <= -160000000.0: tmp = x * (t_1 / y) elif b <= 0.78: tmp = math.pow(a, (t + -1.0)) * (x / y) else: tmp = (x * t_1) / y return tmp
function code(x, y, z, t, a, b) t_1 = exp(Float64(-b)) tmp = 0.0 if (b <= -160000000.0) tmp = Float64(x * Float64(t_1 / y)); elseif (b <= 0.78) tmp = Float64((a ^ Float64(t + -1.0)) * Float64(x / y)); else tmp = Float64(Float64(x * t_1) / y); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = exp(-b); tmp = 0.0; if (b <= -160000000.0) tmp = x * (t_1 / y); elseif (b <= 0.78) tmp = (a ^ (t + -1.0)) * (x / y); else tmp = (x * t_1) / y; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[Exp[(-b)], $MachinePrecision]}, If[LessEqual[b, -160000000.0], N[(x * N[(t$95$1 / y), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 0.78], N[(N[Power[a, N[(t + -1.0), $MachinePrecision]], $MachinePrecision] * N[(x / y), $MachinePrecision]), $MachinePrecision], N[(N[(x * t$95$1), $MachinePrecision] / y), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := e^{-b}\\
\mathbf{if}\;b \leq -160000000:\\
\;\;\;\;x \cdot \frac{t\_1}{y}\\
\mathbf{elif}\;b \leq 0.78:\\
\;\;\;\;{a}^{\left(t + -1\right)} \cdot \frac{x}{y}\\
\mathbf{else}:\\
\;\;\;\;\frac{x \cdot t\_1}{y}\\
\end{array}
\end{array}
if b < -1.6e8Initial program 100.0%
Taylor expanded in y around 0
lower-*.f64N/A
rem-exp-logN/A
lower-log.f64N/A
rem-exp-logN/A
sub-negN/A
metadata-evalN/A
lower-+.f6493.8
Applied rewrites93.8%
Taylor expanded in b around inf
neg-mul-1N/A
lower-neg.f6487.7
Applied rewrites87.7%
lift-/.f64N/A
lift-*.f64N/A
associate-/l*N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f6487.7
Applied rewrites87.7%
if -1.6e8 < b < 0.78000000000000003Initial program 97.3%
Taylor expanded in y around 0
*-commutativeN/A
exp-diffN/A
associate-*l/N/A
associate-/l/N/A
lower-/.f64N/A
lower-*.f64N/A
exp-prodN/A
lower-pow.f64N/A
rem-exp-logN/A
sub-negN/A
metadata-evalN/A
lower-+.f64N/A
lower-*.f64N/A
lower-exp.f6472.5
Applied rewrites72.5%
Taylor expanded in b around 0
Applied rewrites67.8%
if 0.78000000000000003 < b Initial program 100.0%
Taylor expanded in y around 0
lower-*.f64N/A
rem-exp-logN/A
lower-log.f64N/A
rem-exp-logN/A
sub-negN/A
metadata-evalN/A
lower-+.f6480.0
Applied rewrites80.0%
Taylor expanded in b around inf
neg-mul-1N/A
lower-neg.f6474.4
Applied rewrites74.4%
Final simplification74.6%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (exp (- b))))
(if (<= b -16000.0)
(* x (/ t_1 y))
(if (<= b 0.78) (/ (* (/ x y) (+ b -1.0)) (- a)) (/ (* x t_1) y)))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = exp(-b);
double tmp;
if (b <= -16000.0) {
tmp = x * (t_1 / y);
} else if (b <= 0.78) {
tmp = ((x / y) * (b + -1.0)) / -a;
} else {
tmp = (x * t_1) / y;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: t_1
real(8) :: tmp
t_1 = exp(-b)
if (b <= (-16000.0d0)) then
tmp = x * (t_1 / y)
else if (b <= 0.78d0) then
tmp = ((x / y) * (b + (-1.0d0))) / -a
else
tmp = (x * t_1) / y
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double t_1 = Math.exp(-b);
double tmp;
if (b <= -16000.0) {
tmp = x * (t_1 / y);
} else if (b <= 0.78) {
tmp = ((x / y) * (b + -1.0)) / -a;
} else {
tmp = (x * t_1) / y;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = math.exp(-b) tmp = 0 if b <= -16000.0: tmp = x * (t_1 / y) elif b <= 0.78: tmp = ((x / y) * (b + -1.0)) / -a else: tmp = (x * t_1) / y return tmp
function code(x, y, z, t, a, b) t_1 = exp(Float64(-b)) tmp = 0.0 if (b <= -16000.0) tmp = Float64(x * Float64(t_1 / y)); elseif (b <= 0.78) tmp = Float64(Float64(Float64(x / y) * Float64(b + -1.0)) / Float64(-a)); else tmp = Float64(Float64(x * t_1) / y); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = exp(-b); tmp = 0.0; if (b <= -16000.0) tmp = x * (t_1 / y); elseif (b <= 0.78) tmp = ((x / y) * (b + -1.0)) / -a; else tmp = (x * t_1) / y; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[Exp[(-b)], $MachinePrecision]}, If[LessEqual[b, -16000.0], N[(x * N[(t$95$1 / y), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 0.78], N[(N[(N[(x / y), $MachinePrecision] * N[(b + -1.0), $MachinePrecision]), $MachinePrecision] / (-a)), $MachinePrecision], N[(N[(x * t$95$1), $MachinePrecision] / y), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := e^{-b}\\
\mathbf{if}\;b \leq -16000:\\
\;\;\;\;x \cdot \frac{t\_1}{y}\\
\mathbf{elif}\;b \leq 0.78:\\
\;\;\;\;\frac{\frac{x}{y} \cdot \left(b + -1\right)}{-a}\\
\mathbf{else}:\\
\;\;\;\;\frac{x \cdot t\_1}{y}\\
\end{array}
\end{array}
if b < -16000Initial program 100.0%
Taylor expanded in y around 0
lower-*.f64N/A
rem-exp-logN/A
lower-log.f64N/A
rem-exp-logN/A
sub-negN/A
metadata-evalN/A
lower-+.f6493.8
Applied rewrites93.8%
Taylor expanded in b around inf
neg-mul-1N/A
lower-neg.f6487.7
Applied rewrites87.7%
lift-/.f64N/A
lift-*.f64N/A
associate-/l*N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f6487.7
Applied rewrites87.7%
if -16000 < b < 0.78000000000000003Initial program 97.3%
Taylor expanded in y around 0
*-commutativeN/A
exp-diffN/A
associate-*l/N/A
associate-/l/N/A
lower-/.f64N/A
lower-*.f64N/A
exp-prodN/A
lower-pow.f64N/A
rem-exp-logN/A
sub-negN/A
metadata-evalN/A
lower-+.f64N/A
lower-*.f64N/A
lower-exp.f6472.5
Applied rewrites72.5%
Taylor expanded in t around 0
Applied rewrites43.9%
Taylor expanded in b around 0
Applied rewrites38.1%
Taylor expanded in a around -inf
Applied rewrites46.3%
if 0.78000000000000003 < b Initial program 100.0%
Taylor expanded in y around 0
lower-*.f64N/A
rem-exp-logN/A
lower-log.f64N/A
rem-exp-logN/A
sub-negN/A
metadata-evalN/A
lower-+.f6480.0
Applied rewrites80.0%
Taylor expanded in b around inf
neg-mul-1N/A
lower-neg.f6474.4
Applied rewrites74.4%
Final simplification64.2%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (* x (/ (exp (- b)) y))))
(if (<= b -16000.0)
t_1
(if (<= b 0.78) (/ (* (/ x y) (+ b -1.0)) (- a)) t_1))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = x * (exp(-b) / y);
double tmp;
if (b <= -16000.0) {
tmp = t_1;
} else if (b <= 0.78) {
tmp = ((x / y) * (b + -1.0)) / -a;
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: t_1
real(8) :: tmp
t_1 = x * (exp(-b) / y)
if (b <= (-16000.0d0)) then
tmp = t_1
else if (b <= 0.78d0) then
tmp = ((x / y) * (b + (-1.0d0))) / -a
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double t_1 = x * (Math.exp(-b) / y);
double tmp;
if (b <= -16000.0) {
tmp = t_1;
} else if (b <= 0.78) {
tmp = ((x / y) * (b + -1.0)) / -a;
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = x * (math.exp(-b) / y) tmp = 0 if b <= -16000.0: tmp = t_1 elif b <= 0.78: tmp = ((x / y) * (b + -1.0)) / -a else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(x * Float64(exp(Float64(-b)) / y)) tmp = 0.0 if (b <= -16000.0) tmp = t_1; elseif (b <= 0.78) tmp = Float64(Float64(Float64(x / y) * Float64(b + -1.0)) / Float64(-a)); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = x * (exp(-b) / y); tmp = 0.0; if (b <= -16000.0) tmp = t_1; elseif (b <= 0.78) tmp = ((x / y) * (b + -1.0)) / -a; else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(x * N[(N[Exp[(-b)], $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[b, -16000.0], t$95$1, If[LessEqual[b, 0.78], N[(N[(N[(x / y), $MachinePrecision] * N[(b + -1.0), $MachinePrecision]), $MachinePrecision] / (-a)), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x \cdot \frac{e^{-b}}{y}\\
\mathbf{if}\;b \leq -16000:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;b \leq 0.78:\\
\;\;\;\;\frac{\frac{x}{y} \cdot \left(b + -1\right)}{-a}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if b < -16000 or 0.78000000000000003 < b Initial program 100.0%
Taylor expanded in y around 0
lower-*.f64N/A
rem-exp-logN/A
lower-log.f64N/A
rem-exp-logN/A
sub-negN/A
metadata-evalN/A
lower-+.f6486.7
Applied rewrites86.7%
Taylor expanded in b around inf
neg-mul-1N/A
lower-neg.f6480.8
Applied rewrites80.8%
lift-/.f64N/A
lift-*.f64N/A
associate-/l*N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f6480.8
Applied rewrites80.8%
if -16000 < b < 0.78000000000000003Initial program 97.3%
Taylor expanded in y around 0
*-commutativeN/A
exp-diffN/A
associate-*l/N/A
associate-/l/N/A
lower-/.f64N/A
lower-*.f64N/A
exp-prodN/A
lower-pow.f64N/A
rem-exp-logN/A
sub-negN/A
metadata-evalN/A
lower-+.f64N/A
lower-*.f64N/A
lower-exp.f6472.5
Applied rewrites72.5%
Taylor expanded in t around 0
Applied rewrites43.9%
Taylor expanded in b around 0
Applied rewrites38.1%
Taylor expanded in a around -inf
Applied rewrites46.3%
Final simplification64.2%
(FPCore (x y z t a b) :precision binary64 (if (<= (log a) -280.0) (* x (/ 1.0 (* y a))) (/ x (* y (fma a b a)))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (log(a) <= -280.0) {
tmp = x * (1.0 / (y * a));
} else {
tmp = x / (y * fma(a, b, a));
}
return tmp;
}
function code(x, y, z, t, a, b) tmp = 0.0 if (log(a) <= -280.0) tmp = Float64(x * Float64(1.0 / Float64(y * a))); else tmp = Float64(x / Float64(y * fma(a, b, a))); end return tmp end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[N[Log[a], $MachinePrecision], -280.0], N[(x * N[(1.0 / N[(y * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x / N[(y * N[(a * b + a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\log a \leq -280:\\
\;\;\;\;x \cdot \frac{1}{y \cdot a}\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{y \cdot \mathsf{fma}\left(a, b, a\right)}\\
\end{array}
\end{array}
if (log.f64 a) < -280Initial program 99.4%
Taylor expanded in t around 0
associate-/l*N/A
lower-*.f64N/A
exp-diffN/A
associate-/l/N/A
lower-/.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
exp-diffN/A
lower-/.f64N/A
*-commutativeN/A
exp-to-powN/A
lower-pow.f64N/A
rem-exp-logN/A
lower-*.f64N/A
lower-exp.f6469.3
Applied rewrites69.3%
Taylor expanded in y around 0
Applied rewrites58.4%
Taylor expanded in b around 0
Applied rewrites29.6%
if -280 < (log.f64 a) Initial program 98.3%
Taylor expanded in y around 0
*-commutativeN/A
exp-diffN/A
associate-*l/N/A
associate-/l/N/A
lower-/.f64N/A
lower-*.f64N/A
exp-prodN/A
lower-pow.f64N/A
rem-exp-logN/A
sub-negN/A
metadata-evalN/A
lower-+.f64N/A
lower-*.f64N/A
lower-exp.f6469.5
Applied rewrites69.5%
Taylor expanded in t around 0
Applied rewrites65.5%
Taylor expanded in b around 0
Applied rewrites38.3%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (* x (* b (* b (/ 0.5 (* y a)))))))
(if (<= b -100000000.0)
t_1
(if (<= b -3.8e-105)
(/ (* (/ x y) (+ b -1.0)) (- a))
(if (<= b -6e-290) t_1 (/ x (* y (fma b (fma a (* b 0.5) a) a))))))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = x * (b * (b * (0.5 / (y * a))));
double tmp;
if (b <= -100000000.0) {
tmp = t_1;
} else if (b <= -3.8e-105) {
tmp = ((x / y) * (b + -1.0)) / -a;
} else if (b <= -6e-290) {
tmp = t_1;
} else {
tmp = x / (y * fma(b, fma(a, (b * 0.5), a), a));
}
return tmp;
}
function code(x, y, z, t, a, b) t_1 = Float64(x * Float64(b * Float64(b * Float64(0.5 / Float64(y * a))))) tmp = 0.0 if (b <= -100000000.0) tmp = t_1; elseif (b <= -3.8e-105) tmp = Float64(Float64(Float64(x / y) * Float64(b + -1.0)) / Float64(-a)); elseif (b <= -6e-290) tmp = t_1; else tmp = Float64(x / Float64(y * fma(b, fma(a, Float64(b * 0.5), a), a))); end return tmp end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(x * N[(b * N[(b * N[(0.5 / N[(y * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[b, -100000000.0], t$95$1, If[LessEqual[b, -3.8e-105], N[(N[(N[(x / y), $MachinePrecision] * N[(b + -1.0), $MachinePrecision]), $MachinePrecision] / (-a)), $MachinePrecision], If[LessEqual[b, -6e-290], t$95$1, N[(x / N[(y * N[(b * N[(a * N[(b * 0.5), $MachinePrecision] + a), $MachinePrecision] + a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x \cdot \left(b \cdot \left(b \cdot \frac{0.5}{y \cdot a}\right)\right)\\
\mathbf{if}\;b \leq -100000000:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;b \leq -3.8 \cdot 10^{-105}:\\
\;\;\;\;\frac{\frac{x}{y} \cdot \left(b + -1\right)}{-a}\\
\mathbf{elif}\;b \leq -6 \cdot 10^{-290}:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{y \cdot \mathsf{fma}\left(b, \mathsf{fma}\left(a, b \cdot 0.5, a\right), a\right)}\\
\end{array}
\end{array}
if b < -1e8 or -3.7999999999999998e-105 < b < -5.99999999999999985e-290Initial program 99.9%
Taylor expanded in t around 0
associate-/l*N/A
lower-*.f64N/A
exp-diffN/A
associate-/l/N/A
lower-/.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
exp-diffN/A
lower-/.f64N/A
*-commutativeN/A
exp-to-powN/A
lower-pow.f64N/A
rem-exp-logN/A
lower-*.f64N/A
lower-exp.f6468.4
Applied rewrites68.4%
Taylor expanded in y around 0
Applied rewrites67.9%
Taylor expanded in b around 0
Applied rewrites47.8%
Taylor expanded in b around inf
Applied rewrites54.6%
if -1e8 < b < -3.7999999999999998e-105Initial program 95.6%
Taylor expanded in y around 0
*-commutativeN/A
exp-diffN/A
associate-*l/N/A
associate-/l/N/A
lower-/.f64N/A
lower-*.f64N/A
exp-prodN/A
lower-pow.f64N/A
rem-exp-logN/A
sub-negN/A
metadata-evalN/A
lower-+.f64N/A
lower-*.f64N/A
lower-exp.f6474.5
Applied rewrites74.5%
Taylor expanded in t around 0
Applied rewrites57.4%
Taylor expanded in b around 0
Applied rewrites57.4%
Taylor expanded in a around -inf
Applied rewrites72.5%
if -5.99999999999999985e-290 < b Initial program 98.4%
Taylor expanded in y around 0
*-commutativeN/A
exp-diffN/A
associate-*l/N/A
associate-/l/N/A
lower-/.f64N/A
lower-*.f64N/A
exp-prodN/A
lower-pow.f64N/A
rem-exp-logN/A
sub-negN/A
metadata-evalN/A
lower-+.f64N/A
lower-*.f64N/A
lower-exp.f6463.7
Applied rewrites63.7%
Taylor expanded in t around 0
Applied rewrites60.4%
Taylor expanded in b around 0
Applied rewrites48.8%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (* x (* b (* b (/ 0.5 (* y a)))))))
(if (<= b -100000000.0)
t_1
(if (<= b -3.8e-105)
(/ (* (/ x y) (+ b -1.0)) (- a))
(if (<= b -6e-290) t_1 (/ x (* y (fma a b a))))))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = x * (b * (b * (0.5 / (y * a))));
double tmp;
if (b <= -100000000.0) {
tmp = t_1;
} else if (b <= -3.8e-105) {
tmp = ((x / y) * (b + -1.0)) / -a;
} else if (b <= -6e-290) {
tmp = t_1;
} else {
tmp = x / (y * fma(a, b, a));
}
return tmp;
}
function code(x, y, z, t, a, b) t_1 = Float64(x * Float64(b * Float64(b * Float64(0.5 / Float64(y * a))))) tmp = 0.0 if (b <= -100000000.0) tmp = t_1; elseif (b <= -3.8e-105) tmp = Float64(Float64(Float64(x / y) * Float64(b + -1.0)) / Float64(-a)); elseif (b <= -6e-290) tmp = t_1; else tmp = Float64(x / Float64(y * fma(a, b, a))); end return tmp end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(x * N[(b * N[(b * N[(0.5 / N[(y * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[b, -100000000.0], t$95$1, If[LessEqual[b, -3.8e-105], N[(N[(N[(x / y), $MachinePrecision] * N[(b + -1.0), $MachinePrecision]), $MachinePrecision] / (-a)), $MachinePrecision], If[LessEqual[b, -6e-290], t$95$1, N[(x / N[(y * N[(a * b + a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x \cdot \left(b \cdot \left(b \cdot \frac{0.5}{y \cdot a}\right)\right)\\
\mathbf{if}\;b \leq -100000000:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;b \leq -3.8 \cdot 10^{-105}:\\
\;\;\;\;\frac{\frac{x}{y} \cdot \left(b + -1\right)}{-a}\\
\mathbf{elif}\;b \leq -6 \cdot 10^{-290}:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{y \cdot \mathsf{fma}\left(a, b, a\right)}\\
\end{array}
\end{array}
if b < -1e8 or -3.7999999999999998e-105 < b < -5.99999999999999985e-290Initial program 99.9%
Taylor expanded in t around 0
associate-/l*N/A
lower-*.f64N/A
exp-diffN/A
associate-/l/N/A
lower-/.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
exp-diffN/A
lower-/.f64N/A
*-commutativeN/A
exp-to-powN/A
lower-pow.f64N/A
rem-exp-logN/A
lower-*.f64N/A
lower-exp.f6468.4
Applied rewrites68.4%
Taylor expanded in y around 0
Applied rewrites67.9%
Taylor expanded in b around 0
Applied rewrites47.8%
Taylor expanded in b around inf
Applied rewrites54.6%
if -1e8 < b < -3.7999999999999998e-105Initial program 95.6%
Taylor expanded in y around 0
*-commutativeN/A
exp-diffN/A
associate-*l/N/A
associate-/l/N/A
lower-/.f64N/A
lower-*.f64N/A
exp-prodN/A
lower-pow.f64N/A
rem-exp-logN/A
sub-negN/A
metadata-evalN/A
lower-+.f64N/A
lower-*.f64N/A
lower-exp.f6474.5
Applied rewrites74.5%
Taylor expanded in t around 0
Applied rewrites57.4%
Taylor expanded in b around 0
Applied rewrites57.4%
Taylor expanded in a around -inf
Applied rewrites72.5%
if -5.99999999999999985e-290 < b Initial program 98.4%
Taylor expanded in y around 0
*-commutativeN/A
exp-diffN/A
associate-*l/N/A
associate-/l/N/A
lower-/.f64N/A
lower-*.f64N/A
exp-prodN/A
lower-pow.f64N/A
rem-exp-logN/A
sub-negN/A
metadata-evalN/A
lower-+.f64N/A
lower-*.f64N/A
lower-exp.f6463.7
Applied rewrites63.7%
Taylor expanded in t around 0
Applied rewrites60.4%
Taylor expanded in b around 0
Applied rewrites39.0%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (* x (* b (* b (/ 0.5 (* y a)))))))
(if (<= b -410000.0)
t_1
(if (<= b -3.8e-105)
(* (/ x y) (/ 1.0 a))
(if (<= b -6e-290) t_1 (/ x (* y (fma a b a))))))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = x * (b * (b * (0.5 / (y * a))));
double tmp;
if (b <= -410000.0) {
tmp = t_1;
} else if (b <= -3.8e-105) {
tmp = (x / y) * (1.0 / a);
} else if (b <= -6e-290) {
tmp = t_1;
} else {
tmp = x / (y * fma(a, b, a));
}
return tmp;
}
function code(x, y, z, t, a, b) t_1 = Float64(x * Float64(b * Float64(b * Float64(0.5 / Float64(y * a))))) tmp = 0.0 if (b <= -410000.0) tmp = t_1; elseif (b <= -3.8e-105) tmp = Float64(Float64(x / y) * Float64(1.0 / a)); elseif (b <= -6e-290) tmp = t_1; else tmp = Float64(x / Float64(y * fma(a, b, a))); end return tmp end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(x * N[(b * N[(b * N[(0.5 / N[(y * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[b, -410000.0], t$95$1, If[LessEqual[b, -3.8e-105], N[(N[(x / y), $MachinePrecision] * N[(1.0 / a), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, -6e-290], t$95$1, N[(x / N[(y * N[(a * b + a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x \cdot \left(b \cdot \left(b \cdot \frac{0.5}{y \cdot a}\right)\right)\\
\mathbf{if}\;b \leq -410000:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;b \leq -3.8 \cdot 10^{-105}:\\
\;\;\;\;\frac{x}{y} \cdot \frac{1}{a}\\
\mathbf{elif}\;b \leq -6 \cdot 10^{-290}:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{y \cdot \mathsf{fma}\left(a, b, a\right)}\\
\end{array}
\end{array}
if b < -4.1e5 or -3.7999999999999998e-105 < b < -5.99999999999999985e-290Initial program 99.9%
Taylor expanded in t around 0
associate-/l*N/A
lower-*.f64N/A
exp-diffN/A
associate-/l/N/A
lower-/.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
exp-diffN/A
lower-/.f64N/A
*-commutativeN/A
exp-to-powN/A
lower-pow.f64N/A
rem-exp-logN/A
lower-*.f64N/A
lower-exp.f6468.4
Applied rewrites68.4%
Taylor expanded in y around 0
Applied rewrites67.9%
Taylor expanded in b around 0
Applied rewrites47.8%
Taylor expanded in b around inf
Applied rewrites54.6%
if -4.1e5 < b < -3.7999999999999998e-105Initial program 95.6%
Taylor expanded in y around 0
*-commutativeN/A
exp-diffN/A
associate-*l/N/A
associate-/l/N/A
lower-/.f64N/A
lower-*.f64N/A
exp-prodN/A
lower-pow.f64N/A
rem-exp-logN/A
sub-negN/A
metadata-evalN/A
lower-+.f64N/A
lower-*.f64N/A
lower-exp.f6474.5
Applied rewrites74.5%
Taylor expanded in b around 0
Applied rewrites75.1%
Taylor expanded in t around 0
Applied rewrites71.4%
if -5.99999999999999985e-290 < b Initial program 98.4%
Taylor expanded in y around 0
*-commutativeN/A
exp-diffN/A
associate-*l/N/A
associate-/l/N/A
lower-/.f64N/A
lower-*.f64N/A
exp-prodN/A
lower-pow.f64N/A
rem-exp-logN/A
sub-negN/A
metadata-evalN/A
lower-+.f64N/A
lower-*.f64N/A
lower-exp.f6463.7
Applied rewrites63.7%
Taylor expanded in t around 0
Applied rewrites60.4%
Taylor expanded in b around 0
Applied rewrites39.0%
Final simplification48.2%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (* (- b) (/ x (* y a)))))
(if (<= b -6.4e+49)
t_1
(if (<= b -1.78e-102)
(* (/ x y) (/ 1.0 a))
(if (<= b -5.6e-301) t_1 (/ x (* y (fma a b a))))))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = -b * (x / (y * a));
double tmp;
if (b <= -6.4e+49) {
tmp = t_1;
} else if (b <= -1.78e-102) {
tmp = (x / y) * (1.0 / a);
} else if (b <= -5.6e-301) {
tmp = t_1;
} else {
tmp = x / (y * fma(a, b, a));
}
return tmp;
}
function code(x, y, z, t, a, b) t_1 = Float64(Float64(-b) * Float64(x / Float64(y * a))) tmp = 0.0 if (b <= -6.4e+49) tmp = t_1; elseif (b <= -1.78e-102) tmp = Float64(Float64(x / y) * Float64(1.0 / a)); elseif (b <= -5.6e-301) tmp = t_1; else tmp = Float64(x / Float64(y * fma(a, b, a))); end return tmp end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[((-b) * N[(x / N[(y * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[b, -6.4e+49], t$95$1, If[LessEqual[b, -1.78e-102], N[(N[(x / y), $MachinePrecision] * N[(1.0 / a), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, -5.6e-301], t$95$1, N[(x / N[(y * N[(a * b + a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(-b\right) \cdot \frac{x}{y \cdot a}\\
\mathbf{if}\;b \leq -6.4 \cdot 10^{+49}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;b \leq -1.78 \cdot 10^{-102}:\\
\;\;\;\;\frac{x}{y} \cdot \frac{1}{a}\\
\mathbf{elif}\;b \leq -5.6 \cdot 10^{-301}:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{y \cdot \mathsf{fma}\left(a, b, a\right)}\\
\end{array}
\end{array}
if b < -6.40000000000000028e49 or -1.78e-102 < b < -5.6000000000000002e-301Initial program 99.9%
Taylor expanded in y around 0
*-commutativeN/A
exp-diffN/A
associate-*l/N/A
associate-/l/N/A
lower-/.f64N/A
lower-*.f64N/A
exp-prodN/A
lower-pow.f64N/A
rem-exp-logN/A
sub-negN/A
metadata-evalN/A
lower-+.f64N/A
lower-*.f64N/A
lower-exp.f6472.1
Applied rewrites72.1%
Taylor expanded in t around 0
Applied rewrites65.6%
Taylor expanded in b around 0
Applied rewrites34.5%
Taylor expanded in b around inf
Applied rewrites39.5%
if -6.40000000000000028e49 < b < -1.78e-102Initial program 96.6%
Taylor expanded in y around 0
*-commutativeN/A
exp-diffN/A
associate-*l/N/A
associate-/l/N/A
lower-/.f64N/A
lower-*.f64N/A
exp-prodN/A
lower-pow.f64N/A
rem-exp-logN/A
sub-negN/A
metadata-evalN/A
lower-+.f64N/A
lower-*.f64N/A
lower-exp.f6474.7
Applied rewrites74.7%
Taylor expanded in b around 0
Applied rewrites66.6%
Taylor expanded in t around 0
Applied rewrites61.0%
if -5.6000000000000002e-301 < b Initial program 98.3%
Taylor expanded in y around 0
*-commutativeN/A
exp-diffN/A
associate-*l/N/A
associate-/l/N/A
lower-/.f64N/A
lower-*.f64N/A
exp-prodN/A
lower-pow.f64N/A
rem-exp-logN/A
sub-negN/A
metadata-evalN/A
lower-+.f64N/A
lower-*.f64N/A
lower-exp.f6463.3
Applied rewrites63.3%
Taylor expanded in t around 0
Applied rewrites60.7%
Taylor expanded in b around 0
Applied rewrites38.5%
Final simplification41.8%
(FPCore (x y z t a b) :precision binary64 (if (<= b -0.82) (* (- b) (/ x (* y a))) (/ x (* y (fma a b a)))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (b <= -0.82) {
tmp = -b * (x / (y * a));
} else {
tmp = x / (y * fma(a, b, a));
}
return tmp;
}
function code(x, y, z, t, a, b) tmp = 0.0 if (b <= -0.82) tmp = Float64(Float64(-b) * Float64(x / Float64(y * a))); else tmp = Float64(x / Float64(y * fma(a, b, a))); end return tmp end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[b, -0.82], N[((-b) * N[(x / N[(y * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x / N[(y * N[(a * b + a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq -0.82:\\
\;\;\;\;\left(-b\right) \cdot \frac{x}{y \cdot a}\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{y \cdot \mathsf{fma}\left(a, b, a\right)}\\
\end{array}
\end{array}
if b < -0.819999999999999951Initial program 100.0%
Taylor expanded in y around 0
*-commutativeN/A
exp-diffN/A
associate-*l/N/A
associate-/l/N/A
lower-/.f64N/A
lower-*.f64N/A
exp-prodN/A
lower-pow.f64N/A
rem-exp-logN/A
sub-negN/A
metadata-evalN/A
lower-+.f64N/A
lower-*.f64N/A
lower-exp.f6470.8
Applied rewrites70.8%
Taylor expanded in t around 0
Applied rewrites87.9%
Taylor expanded in b around 0
Applied rewrites33.4%
Taylor expanded in b around inf
Applied rewrites33.4%
if -0.819999999999999951 < b Initial program 98.3%
Taylor expanded in y around 0
*-commutativeN/A
exp-diffN/A
associate-*l/N/A
associate-/l/N/A
lower-/.f64N/A
lower-*.f64N/A
exp-prodN/A
lower-pow.f64N/A
rem-exp-logN/A
sub-negN/A
metadata-evalN/A
lower-+.f64N/A
lower-*.f64N/A
lower-exp.f6467.2
Applied rewrites67.2%
Taylor expanded in t around 0
Applied rewrites54.6%
Taylor expanded in b around 0
Applied rewrites39.9%
Final simplification38.2%
(FPCore (x y z t a b) :precision binary64 (* x (/ 1.0 (* y a))))
double code(double x, double y, double z, double t, double a, double b) {
return x * (1.0 / (y * a));
}
real(8) function code(x, y, z, t, a, b)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
code = x * (1.0d0 / (y * a))
end function
public static double code(double x, double y, double z, double t, double a, double b) {
return x * (1.0 / (y * a));
}
def code(x, y, z, t, a, b): return x * (1.0 / (y * a))
function code(x, y, z, t, a, b) return Float64(x * Float64(1.0 / Float64(y * a))) end
function tmp = code(x, y, z, t, a, b) tmp = x * (1.0 / (y * a)); end
code[x_, y_, z_, t_, a_, b_] := N[(x * N[(1.0 / N[(y * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x \cdot \frac{1}{y \cdot a}
\end{array}
Initial program 98.7%
Taylor expanded in t around 0
associate-/l*N/A
lower-*.f64N/A
exp-diffN/A
associate-/l/N/A
lower-/.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
exp-diffN/A
lower-/.f64N/A
*-commutativeN/A
exp-to-powN/A
lower-pow.f64N/A
rem-exp-logN/A
lower-*.f64N/A
lower-exp.f6471.0
Applied rewrites71.0%
Taylor expanded in y around 0
Applied rewrites63.0%
Taylor expanded in b around 0
Applied rewrites28.9%
(FPCore (x y z t a b) :precision binary64 (/ x (* y a)))
double code(double x, double y, double z, double t, double a, double b) {
return x / (y * a);
}
real(8) function code(x, y, z, t, a, b)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
code = x / (y * a)
end function
public static double code(double x, double y, double z, double t, double a, double b) {
return x / (y * a);
}
def code(x, y, z, t, a, b): return x / (y * a)
function code(x, y, z, t, a, b) return Float64(x / Float64(y * a)) end
function tmp = code(x, y, z, t, a, b) tmp = x / (y * a); end
code[x_, y_, z_, t_, a_, b_] := N[(x / N[(y * a), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{x}{y \cdot a}
\end{array}
Initial program 98.7%
Taylor expanded in y around 0
*-commutativeN/A
exp-diffN/A
associate-*l/N/A
associate-/l/N/A
lower-/.f64N/A
lower-*.f64N/A
exp-prodN/A
lower-pow.f64N/A
rem-exp-logN/A
sub-negN/A
metadata-evalN/A
lower-+.f64N/A
lower-*.f64N/A
lower-exp.f6468.1
Applied rewrites68.1%
Taylor expanded in t around 0
Applied rewrites63.0%
Taylor expanded in b around 0
Applied rewrites28.6%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (pow a (- t 1.0)))
(t_2 (/ (* x (/ t_1 y)) (- (+ b 1.0) (* y (log z))))))
(if (< t -0.8845848504127471)
t_2
(if (< t 852031.2288374073)
(/ (* (/ x y) t_1) (exp (- b (* (log z) y))))
t_2))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = pow(a, (t - 1.0));
double t_2 = (x * (t_1 / y)) / ((b + 1.0) - (y * log(z)));
double tmp;
if (t < -0.8845848504127471) {
tmp = t_2;
} else if (t < 852031.2288374073) {
tmp = ((x / y) * t_1) / exp((b - (log(z) * y)));
} else {
tmp = t_2;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: t_1
real(8) :: t_2
real(8) :: tmp
t_1 = a ** (t - 1.0d0)
t_2 = (x * (t_1 / y)) / ((b + 1.0d0) - (y * log(z)))
if (t < (-0.8845848504127471d0)) then
tmp = t_2
else if (t < 852031.2288374073d0) then
tmp = ((x / y) * t_1) / exp((b - (log(z) * y)))
else
tmp = t_2
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double t_1 = Math.pow(a, (t - 1.0));
double t_2 = (x * (t_1 / y)) / ((b + 1.0) - (y * Math.log(z)));
double tmp;
if (t < -0.8845848504127471) {
tmp = t_2;
} else if (t < 852031.2288374073) {
tmp = ((x / y) * t_1) / Math.exp((b - (Math.log(z) * y)));
} else {
tmp = t_2;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = math.pow(a, (t - 1.0)) t_2 = (x * (t_1 / y)) / ((b + 1.0) - (y * math.log(z))) tmp = 0 if t < -0.8845848504127471: tmp = t_2 elif t < 852031.2288374073: tmp = ((x / y) * t_1) / math.exp((b - (math.log(z) * y))) else: tmp = t_2 return tmp
function code(x, y, z, t, a, b) t_1 = a ^ Float64(t - 1.0) t_2 = Float64(Float64(x * Float64(t_1 / y)) / Float64(Float64(b + 1.0) - Float64(y * log(z)))) tmp = 0.0 if (t < -0.8845848504127471) tmp = t_2; elseif (t < 852031.2288374073) tmp = Float64(Float64(Float64(x / y) * t_1) / exp(Float64(b - Float64(log(z) * y)))); else tmp = t_2; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = a ^ (t - 1.0); t_2 = (x * (t_1 / y)) / ((b + 1.0) - (y * log(z))); tmp = 0.0; if (t < -0.8845848504127471) tmp = t_2; elseif (t < 852031.2288374073) tmp = ((x / y) * t_1) / exp((b - (log(z) * y))); else tmp = t_2; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[Power[a, N[(t - 1.0), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$2 = N[(N[(x * N[(t$95$1 / y), $MachinePrecision]), $MachinePrecision] / N[(N[(b + 1.0), $MachinePrecision] - N[(y * N[Log[z], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[Less[t, -0.8845848504127471], t$95$2, If[Less[t, 852031.2288374073], N[(N[(N[(x / y), $MachinePrecision] * t$95$1), $MachinePrecision] / N[Exp[N[(b - N[(N[Log[z], $MachinePrecision] * y), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], t$95$2]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := {a}^{\left(t - 1\right)}\\
t_2 := \frac{x \cdot \frac{t\_1}{y}}{\left(b + 1\right) - y \cdot \log z}\\
\mathbf{if}\;t < -0.8845848504127471:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;t < 852031.2288374073:\\
\;\;\;\;\frac{\frac{x}{y} \cdot t\_1}{e^{b - \log z \cdot y}}\\
\mathbf{else}:\\
\;\;\;\;t\_2\\
\end{array}
\end{array}
herbie shell --seed 2024233
(FPCore (x y z t a b)
:name "Numeric.SpecFunctions:incompleteBetaWorker from math-functions-0.1.5.2, A"
:precision binary64
:alt
(! :herbie-platform default (if (< t -8845848504127471/10000000000000000) (/ (* x (/ (pow a (- t 1)) y)) (- (+ b 1) (* y (log z)))) (if (< t 8520312288374073/10000000000) (/ (* (/ x y) (pow a (- t 1))) (exp (- b (* (log z) y)))) (/ (* x (/ (pow a (- t 1)) y)) (- (+ b 1) (* y (log z)))))))
(/ (* x (exp (- (+ (* y (log z)) (* (- t 1.0) (log a))) b))) y))