
(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 (/ (* (exp (- (+ (* (log a) (- t 1.0)) (* (log z) y)) b)) x) y))
double code(double x, double y, double z, double t, double a, double b) {
return (exp((((log(a) * (t - 1.0)) + (log(z) * y)) - b)) * x) / 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 = (exp((((log(a) * (t - 1.0d0)) + (log(z) * y)) - b)) * x) / y
end function
public static double code(double x, double y, double z, double t, double a, double b) {
return (Math.exp((((Math.log(a) * (t - 1.0)) + (Math.log(z) * y)) - b)) * x) / y;
}
def code(x, y, z, t, a, b): return (math.exp((((math.log(a) * (t - 1.0)) + (math.log(z) * y)) - b)) * x) / y
function code(x, y, z, t, a, b) return Float64(Float64(exp(Float64(Float64(Float64(log(a) * Float64(t - 1.0)) + Float64(log(z) * y)) - b)) * x) / y) end
function tmp = code(x, y, z, t, a, b) tmp = (exp((((log(a) * (t - 1.0)) + (log(z) * y)) - b)) * x) / y; end
code[x_, y_, z_, t_, a_, b_] := N[(N[(N[Exp[N[(N[(N[(N[Log[a], $MachinePrecision] * N[(t - 1.0), $MachinePrecision]), $MachinePrecision] + N[(N[Log[z], $MachinePrecision] * y), $MachinePrecision]), $MachinePrecision] - b), $MachinePrecision]], $MachinePrecision] * x), $MachinePrecision] / y), $MachinePrecision]
\begin{array}{l}
\\
\frac{e^{\left(\log a \cdot \left(t - 1\right) + \log z \cdot y\right) - b} \cdot x}{y}
\end{array}
Initial program 97.4%
Final simplification97.4%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (/ (* (exp (- (+ (* (log a) (- t 1.0)) (* (log z) y)) b)) x) y)))
(if (<= t_1 -5e+305)
(/ (/ x a) y)
(if (<= t_1 1e-205)
(/ (/ x (fma (fma (* (fma 0.16666666666666666 b 0.5) y) b y) b y)) a)
(* (/ (* (/ 1.0 a) 1.0) (* 1.0 y)) x)))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (exp((((log(a) * (t - 1.0)) + (log(z) * y)) - b)) * x) / y;
double tmp;
if (t_1 <= -5e+305) {
tmp = (x / a) / y;
} else if (t_1 <= 1e-205) {
tmp = (x / fma(fma((fma(0.16666666666666666, b, 0.5) * y), b, y), b, y)) / a;
} else {
tmp = (((1.0 / a) * 1.0) / (1.0 * y)) * x;
}
return tmp;
}
function code(x, y, z, t, a, b) t_1 = Float64(Float64(exp(Float64(Float64(Float64(log(a) * Float64(t - 1.0)) + Float64(log(z) * y)) - b)) * x) / y) tmp = 0.0 if (t_1 <= -5e+305) tmp = Float64(Float64(x / a) / y); elseif (t_1 <= 1e-205) tmp = Float64(Float64(x / fma(fma(Float64(fma(0.16666666666666666, b, 0.5) * y), b, y), b, y)) / a); else tmp = Float64(Float64(Float64(Float64(1.0 / a) * 1.0) / Float64(1.0 * y)) * x); end return tmp end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(N[Exp[N[(N[(N[(N[Log[a], $MachinePrecision] * N[(t - 1.0), $MachinePrecision]), $MachinePrecision] + N[(N[Log[z], $MachinePrecision] * y), $MachinePrecision]), $MachinePrecision] - b), $MachinePrecision]], $MachinePrecision] * x), $MachinePrecision] / y), $MachinePrecision]}, If[LessEqual[t$95$1, -5e+305], N[(N[(x / a), $MachinePrecision] / y), $MachinePrecision], If[LessEqual[t$95$1, 1e-205], N[(N[(x / N[(N[(N[(N[(0.16666666666666666 * b + 0.5), $MachinePrecision] * y), $MachinePrecision] * b + y), $MachinePrecision] * b + y), $MachinePrecision]), $MachinePrecision] / a), $MachinePrecision], N[(N[(N[(N[(1.0 / a), $MachinePrecision] * 1.0), $MachinePrecision] / N[(1.0 * y), $MachinePrecision]), $MachinePrecision] * x), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{e^{\left(\log a \cdot \left(t - 1\right) + \log z \cdot y\right) - b} \cdot x}{y}\\
\mathbf{if}\;t\_1 \leq -5 \cdot 10^{+305}:\\
\;\;\;\;\frac{\frac{x}{a}}{y}\\
\mathbf{elif}\;t\_1 \leq 10^{-205}:\\
\;\;\;\;\frac{\frac{x}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(0.16666666666666666, b, 0.5\right) \cdot y, b, y\right), b, y\right)}}{a}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{1}{a} \cdot 1}{1 \cdot y} \cdot x\\
\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) < -5.00000000000000009e305Initial program 99.8%
Taylor expanded in b around 0
exp-sumN/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
*-commutativeN/A
exp-to-powN/A
lower-pow.f64N/A
exp-prodN/A
lower-pow.f64N/A
rem-exp-logN/A
lower--.f6468.2
Applied rewrites68.2%
Taylor expanded in t around 0
Applied rewrites68.9%
Taylor expanded in y around 0
Applied rewrites21.1%
if -5.00000000000000009e305 < (/.f64 (*.f64 x (exp.f64 (-.f64 (+.f64 (*.f64 y (log.f64 z)) (*.f64 (-.f64 t #s(literal 1 binary64)) (log.f64 a))) b))) y) < 1e-205Initial program 95.6%
lift-/.f64N/A
lift-*.f64N/A
*-commutativeN/A
associate-/l*N/A
lift-exp.f64N/A
lift--.f64N/A
exp-diffN/A
clear-numN/A
frac-timesN/A
lower-/.f64N/A
lower-*.f64N/A
lower-*.f64N/A
Applied rewrites80.6%
Taylor expanded in t around 0
*-commutativeN/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-exp.f64N/A
lower-/.f64N/A
lower-pow.f6476.3
Applied rewrites76.3%
Taylor expanded in y around 0
Applied rewrites64.7%
Taylor expanded in b around 0
Applied rewrites56.5%
if 1e-205 < (/.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%
lift-/.f64N/A
lift-*.f64N/A
associate-/l*N/A
*-commutativeN/A
lower-*.f64N/A
Applied rewrites78.5%
Taylor expanded in y around 0
Applied rewrites78.7%
Taylor expanded in b around 0
Applied rewrites76.4%
Taylor expanded in t around 0
lower-/.f6440.2
Applied rewrites40.2%
Final simplification44.5%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (/ (* (exp (- (+ (* (log a) (- t 1.0)) (* (log z) y)) b)) x) y)))
(if (<= t_1 -5e+305)
(/ (/ x a) y)
(if (<= t_1 1e-205)
(/ (/ x (fma (* (fma 0.5 b 1.0) y) b y)) a)
(* (/ (* (/ 1.0 a) 1.0) (* 1.0 y)) x)))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (exp((((log(a) * (t - 1.0)) + (log(z) * y)) - b)) * x) / y;
double tmp;
if (t_1 <= -5e+305) {
tmp = (x / a) / y;
} else if (t_1 <= 1e-205) {
tmp = (x / fma((fma(0.5, b, 1.0) * y), b, y)) / a;
} else {
tmp = (((1.0 / a) * 1.0) / (1.0 * y)) * x;
}
return tmp;
}
function code(x, y, z, t, a, b) t_1 = Float64(Float64(exp(Float64(Float64(Float64(log(a) * Float64(t - 1.0)) + Float64(log(z) * y)) - b)) * x) / y) tmp = 0.0 if (t_1 <= -5e+305) tmp = Float64(Float64(x / a) / y); elseif (t_1 <= 1e-205) tmp = Float64(Float64(x / fma(Float64(fma(0.5, b, 1.0) * y), b, y)) / a); else tmp = Float64(Float64(Float64(Float64(1.0 / a) * 1.0) / Float64(1.0 * y)) * x); end return tmp end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(N[Exp[N[(N[(N[(N[Log[a], $MachinePrecision] * N[(t - 1.0), $MachinePrecision]), $MachinePrecision] + N[(N[Log[z], $MachinePrecision] * y), $MachinePrecision]), $MachinePrecision] - b), $MachinePrecision]], $MachinePrecision] * x), $MachinePrecision] / y), $MachinePrecision]}, If[LessEqual[t$95$1, -5e+305], N[(N[(x / a), $MachinePrecision] / y), $MachinePrecision], If[LessEqual[t$95$1, 1e-205], N[(N[(x / N[(N[(N[(0.5 * b + 1.0), $MachinePrecision] * y), $MachinePrecision] * b + y), $MachinePrecision]), $MachinePrecision] / a), $MachinePrecision], N[(N[(N[(N[(1.0 / a), $MachinePrecision] * 1.0), $MachinePrecision] / N[(1.0 * y), $MachinePrecision]), $MachinePrecision] * x), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{e^{\left(\log a \cdot \left(t - 1\right) + \log z \cdot y\right) - b} \cdot x}{y}\\
\mathbf{if}\;t\_1 \leq -5 \cdot 10^{+305}:\\
\;\;\;\;\frac{\frac{x}{a}}{y}\\
\mathbf{elif}\;t\_1 \leq 10^{-205}:\\
\;\;\;\;\frac{\frac{x}{\mathsf{fma}\left(\mathsf{fma}\left(0.5, b, 1\right) \cdot y, b, y\right)}}{a}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{1}{a} \cdot 1}{1 \cdot y} \cdot x\\
\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) < -5.00000000000000009e305Initial program 99.8%
Taylor expanded in b around 0
exp-sumN/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
*-commutativeN/A
exp-to-powN/A
lower-pow.f64N/A
exp-prodN/A
lower-pow.f64N/A
rem-exp-logN/A
lower--.f6468.2
Applied rewrites68.2%
Taylor expanded in t around 0
Applied rewrites68.9%
Taylor expanded in y around 0
Applied rewrites21.1%
if -5.00000000000000009e305 < (/.f64 (*.f64 x (exp.f64 (-.f64 (+.f64 (*.f64 y (log.f64 z)) (*.f64 (-.f64 t #s(literal 1 binary64)) (log.f64 a))) b))) y) < 1e-205Initial program 95.6%
lift-/.f64N/A
lift-*.f64N/A
*-commutativeN/A
associate-/l*N/A
lift-exp.f64N/A
lift--.f64N/A
exp-diffN/A
clear-numN/A
frac-timesN/A
lower-/.f64N/A
lower-*.f64N/A
lower-*.f64N/A
Applied rewrites80.6%
Taylor expanded in t around 0
*-commutativeN/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-exp.f64N/A
lower-/.f64N/A
lower-pow.f6476.3
Applied rewrites76.3%
Taylor expanded in y around 0
Applied rewrites64.7%
Taylor expanded in b around 0
Applied rewrites53.3%
if 1e-205 < (/.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%
lift-/.f64N/A
lift-*.f64N/A
associate-/l*N/A
*-commutativeN/A
lower-*.f64N/A
Applied rewrites78.5%
Taylor expanded in y around 0
Applied rewrites78.7%
Taylor expanded in b around 0
Applied rewrites76.4%
Taylor expanded in t around 0
lower-/.f6440.2
Applied rewrites40.2%
Final simplification42.8%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (/ (* (exp (- (+ (* (log a) (- t 1.0)) (* (log z) y)) b)) x) y)))
(if (<= t_1 -5e+305)
(/ (/ x a) y)
(if (<= t_1 1e-205)
(/ (/ x (fma b y y)) a)
(* (/ (* (/ 1.0 a) 1.0) (* 1.0 y)) x)))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (exp((((log(a) * (t - 1.0)) + (log(z) * y)) - b)) * x) / y;
double tmp;
if (t_1 <= -5e+305) {
tmp = (x / a) / y;
} else if (t_1 <= 1e-205) {
tmp = (x / fma(b, y, y)) / a;
} else {
tmp = (((1.0 / a) * 1.0) / (1.0 * y)) * x;
}
return tmp;
}
function code(x, y, z, t, a, b) t_1 = Float64(Float64(exp(Float64(Float64(Float64(log(a) * Float64(t - 1.0)) + Float64(log(z) * y)) - b)) * x) / y) tmp = 0.0 if (t_1 <= -5e+305) tmp = Float64(Float64(x / a) / y); elseif (t_1 <= 1e-205) tmp = Float64(Float64(x / fma(b, y, y)) / a); else tmp = Float64(Float64(Float64(Float64(1.0 / a) * 1.0) / Float64(1.0 * y)) * x); end return tmp end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(N[Exp[N[(N[(N[(N[Log[a], $MachinePrecision] * N[(t - 1.0), $MachinePrecision]), $MachinePrecision] + N[(N[Log[z], $MachinePrecision] * y), $MachinePrecision]), $MachinePrecision] - b), $MachinePrecision]], $MachinePrecision] * x), $MachinePrecision] / y), $MachinePrecision]}, If[LessEqual[t$95$1, -5e+305], N[(N[(x / a), $MachinePrecision] / y), $MachinePrecision], If[LessEqual[t$95$1, 1e-205], N[(N[(x / N[(b * y + y), $MachinePrecision]), $MachinePrecision] / a), $MachinePrecision], N[(N[(N[(N[(1.0 / a), $MachinePrecision] * 1.0), $MachinePrecision] / N[(1.0 * y), $MachinePrecision]), $MachinePrecision] * x), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{e^{\left(\log a \cdot \left(t - 1\right) + \log z \cdot y\right) - b} \cdot x}{y}\\
\mathbf{if}\;t\_1 \leq -5 \cdot 10^{+305}:\\
\;\;\;\;\frac{\frac{x}{a}}{y}\\
\mathbf{elif}\;t\_1 \leq 10^{-205}:\\
\;\;\;\;\frac{\frac{x}{\mathsf{fma}\left(b, y, y\right)}}{a}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{1}{a} \cdot 1}{1 \cdot y} \cdot x\\
\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) < -5.00000000000000009e305Initial program 99.8%
Taylor expanded in b around 0
exp-sumN/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
*-commutativeN/A
exp-to-powN/A
lower-pow.f64N/A
exp-prodN/A
lower-pow.f64N/A
rem-exp-logN/A
lower--.f6468.2
Applied rewrites68.2%
Taylor expanded in t around 0
Applied rewrites68.9%
Taylor expanded in y around 0
Applied rewrites21.1%
if -5.00000000000000009e305 < (/.f64 (*.f64 x (exp.f64 (-.f64 (+.f64 (*.f64 y (log.f64 z)) (*.f64 (-.f64 t #s(literal 1 binary64)) (log.f64 a))) b))) y) < 1e-205Initial program 95.6%
lift-/.f64N/A
lift-*.f64N/A
*-commutativeN/A
associate-/l*N/A
lift-exp.f64N/A
lift--.f64N/A
exp-diffN/A
clear-numN/A
frac-timesN/A
lower-/.f64N/A
lower-*.f64N/A
lower-*.f64N/A
Applied rewrites80.6%
Taylor expanded in t around 0
*-commutativeN/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-exp.f64N/A
lower-/.f64N/A
lower-pow.f6476.3
Applied rewrites76.3%
Taylor expanded in y around 0
Applied rewrites64.7%
Taylor expanded in b around 0
Applied rewrites46.0%
if 1e-205 < (/.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%
lift-/.f64N/A
lift-*.f64N/A
associate-/l*N/A
*-commutativeN/A
lower-*.f64N/A
Applied rewrites78.5%
Taylor expanded in y around 0
Applied rewrites78.7%
Taylor expanded in b around 0
Applied rewrites76.4%
Taylor expanded in t around 0
lower-/.f6440.2
Applied rewrites40.2%
Final simplification39.2%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (/ (* (pow a (- t 1.0)) x) y)))
(if (<= t -55000000000000.0)
t_1
(if (<= t 1.9e+82) (* (/ (pow z y) a) (/ x (* (exp b) y))) t_1))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (pow(a, (t - 1.0)) * x) / y;
double tmp;
if (t <= -55000000000000.0) {
tmp = t_1;
} else if (t <= 1.9e+82) {
tmp = (pow(z, y) / a) * (x / (exp(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 = ((a ** (t - 1.0d0)) * x) / y
if (t <= (-55000000000000.0d0)) then
tmp = t_1
else if (t <= 1.9d+82) then
tmp = ((z ** y) / a) * (x / (exp(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 = (Math.pow(a, (t - 1.0)) * x) / y;
double tmp;
if (t <= -55000000000000.0) {
tmp = t_1;
} else if (t <= 1.9e+82) {
tmp = (Math.pow(z, y) / a) * (x / (Math.exp(b) * y));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = (math.pow(a, (t - 1.0)) * x) / y tmp = 0 if t <= -55000000000000.0: tmp = t_1 elif t <= 1.9e+82: tmp = (math.pow(z, y) / a) * (x / (math.exp(b) * y)) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(Float64((a ^ Float64(t - 1.0)) * x) / y) tmp = 0.0 if (t <= -55000000000000.0) tmp = t_1; elseif (t <= 1.9e+82) tmp = Float64(Float64((z ^ y) / a) * Float64(x / Float64(exp(b) * y))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = ((a ^ (t - 1.0)) * x) / y; tmp = 0.0; if (t <= -55000000000000.0) tmp = t_1; elseif (t <= 1.9e+82) tmp = ((z ^ y) / a) * (x / (exp(b) * y)); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(N[Power[a, N[(t - 1.0), $MachinePrecision]], $MachinePrecision] * x), $MachinePrecision] / y), $MachinePrecision]}, If[LessEqual[t, -55000000000000.0], t$95$1, If[LessEqual[t, 1.9e+82], N[(N[(N[Power[z, y], $MachinePrecision] / a), $MachinePrecision] * N[(x / N[(N[Exp[b], $MachinePrecision] * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{{a}^{\left(t - 1\right)} \cdot x}{y}\\
\mathbf{if}\;t \leq -55000000000000:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t \leq 1.9 \cdot 10^{+82}:\\
\;\;\;\;\frac{{z}^{y}}{a} \cdot \frac{x}{e^{b} \cdot y}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if t < -5.5e13 or 1.90000000000000017e82 < t Initial program 100.0%
Taylor expanded in b around 0
exp-sumN/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
*-commutativeN/A
exp-to-powN/A
lower-pow.f64N/A
exp-prodN/A
lower-pow.f64N/A
rem-exp-logN/A
lower--.f6473.2
Applied rewrites73.2%
Taylor expanded in y around 0
Applied rewrites85.4%
if -5.5e13 < t < 1.90000000000000017e82Initial program 95.6%
lift-/.f64N/A
lift-*.f64N/A
*-commutativeN/A
associate-/l*N/A
lift-exp.f64N/A
lift--.f64N/A
exp-diffN/A
clear-numN/A
frac-timesN/A
lower-/.f64N/A
lower-*.f64N/A
lower-*.f64N/A
Applied rewrites87.9%
Taylor expanded in t around 0
*-commutativeN/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-exp.f64N/A
lower-/.f64N/A
lower-pow.f6491.3
Applied rewrites91.3%
Final simplification88.8%
(FPCore (x y z t a b)
:precision binary64
(if (<= y -6.5e+51)
(/ (/ (* (pow z y) x) a) y)
(if (<= y -1.42e-260)
(/ (/ x (* (exp b) y)) a)
(if (<= y 6.5e+35)
(* (/ (* 1.0 (pow a (- t 1.0))) (* 1.0 y)) x)
(/ (* (exp (* (log z) y)) x) y)))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (y <= -6.5e+51) {
tmp = ((pow(z, y) * x) / a) / y;
} else if (y <= -1.42e-260) {
tmp = (x / (exp(b) * y)) / a;
} else if (y <= 6.5e+35) {
tmp = ((1.0 * pow(a, (t - 1.0))) / (1.0 * y)) * x;
} else {
tmp = (exp((log(z) * y)) * x) / 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) :: tmp
if (y <= (-6.5d+51)) then
tmp = (((z ** y) * x) / a) / y
else if (y <= (-1.42d-260)) then
tmp = (x / (exp(b) * y)) / a
else if (y <= 6.5d+35) then
tmp = ((1.0d0 * (a ** (t - 1.0d0))) / (1.0d0 * y)) * x
else
tmp = (exp((log(z) * y)) * x) / y
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 (y <= -6.5e+51) {
tmp = ((Math.pow(z, y) * x) / a) / y;
} else if (y <= -1.42e-260) {
tmp = (x / (Math.exp(b) * y)) / a;
} else if (y <= 6.5e+35) {
tmp = ((1.0 * Math.pow(a, (t - 1.0))) / (1.0 * y)) * x;
} else {
tmp = (Math.exp((Math.log(z) * y)) * x) / y;
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if y <= -6.5e+51: tmp = ((math.pow(z, y) * x) / a) / y elif y <= -1.42e-260: tmp = (x / (math.exp(b) * y)) / a elif y <= 6.5e+35: tmp = ((1.0 * math.pow(a, (t - 1.0))) / (1.0 * y)) * x else: tmp = (math.exp((math.log(z) * y)) * x) / y return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (y <= -6.5e+51) tmp = Float64(Float64(Float64((z ^ y) * x) / a) / y); elseif (y <= -1.42e-260) tmp = Float64(Float64(x / Float64(exp(b) * y)) / a); elseif (y <= 6.5e+35) tmp = Float64(Float64(Float64(1.0 * (a ^ Float64(t - 1.0))) / Float64(1.0 * y)) * x); else tmp = Float64(Float64(exp(Float64(log(z) * y)) * x) / y); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (y <= -6.5e+51) tmp = (((z ^ y) * x) / a) / y; elseif (y <= -1.42e-260) tmp = (x / (exp(b) * y)) / a; elseif (y <= 6.5e+35) tmp = ((1.0 * (a ^ (t - 1.0))) / (1.0 * y)) * x; else tmp = (exp((log(z) * y)) * x) / y; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[y, -6.5e+51], N[(N[(N[(N[Power[z, y], $MachinePrecision] * x), $MachinePrecision] / a), $MachinePrecision] / y), $MachinePrecision], If[LessEqual[y, -1.42e-260], N[(N[(x / N[(N[Exp[b], $MachinePrecision] * y), $MachinePrecision]), $MachinePrecision] / a), $MachinePrecision], If[LessEqual[y, 6.5e+35], N[(N[(N[(1.0 * N[Power[a, N[(t - 1.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / N[(1.0 * y), $MachinePrecision]), $MachinePrecision] * x), $MachinePrecision], N[(N[(N[Exp[N[(N[Log[z], $MachinePrecision] * y), $MachinePrecision]], $MachinePrecision] * x), $MachinePrecision] / y), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -6.5 \cdot 10^{+51}:\\
\;\;\;\;\frac{\frac{{z}^{y} \cdot x}{a}}{y}\\
\mathbf{elif}\;y \leq -1.42 \cdot 10^{-260}:\\
\;\;\;\;\frac{\frac{x}{e^{b} \cdot y}}{a}\\
\mathbf{elif}\;y \leq 6.5 \cdot 10^{+35}:\\
\;\;\;\;\frac{1 \cdot {a}^{\left(t - 1\right)}}{1 \cdot y} \cdot x\\
\mathbf{else}:\\
\;\;\;\;\frac{e^{\log z \cdot y} \cdot x}{y}\\
\end{array}
\end{array}
if y < -6.5e51Initial program 100.0%
Taylor expanded in b around 0
exp-sumN/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
*-commutativeN/A
exp-to-powN/A
lower-pow.f64N/A
exp-prodN/A
lower-pow.f64N/A
rem-exp-logN/A
lower--.f6476.7
Applied rewrites76.7%
Taylor expanded in t around 0
Applied rewrites91.6%
if -6.5e51 < y < -1.42000000000000004e-260Initial program 95.7%
lift-/.f64N/A
lift-*.f64N/A
*-commutativeN/A
associate-/l*N/A
lift-exp.f64N/A
lift--.f64N/A
exp-diffN/A
clear-numN/A
frac-timesN/A
lower-/.f64N/A
lower-*.f64N/A
lower-*.f64N/A
Applied rewrites82.0%
Taylor expanded in t around 0
*-commutativeN/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-exp.f64N/A
lower-/.f64N/A
lower-pow.f6479.0
Applied rewrites79.0%
Taylor expanded in y around 0
Applied rewrites83.0%
if -1.42000000000000004e-260 < y < 6.5000000000000003e35Initial program 95.5%
lift-/.f64N/A
lift-*.f64N/A
associate-/l*N/A
*-commutativeN/A
lower-*.f64N/A
Applied rewrites83.7%
Taylor expanded in y around 0
Applied rewrites84.8%
Taylor expanded in b around 0
Applied rewrites75.5%
if 6.5000000000000003e35 < y Initial program 100.0%
Taylor expanded in y around inf
*-commutativeN/A
lower-*.f64N/A
lower-log.f6487.1
Applied rewrites87.1%
Final simplification82.7%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (/ (* (pow a (- t 1.0)) x) y)))
(if (<= t -55000000000000.0)
t_1
(if (<= t 8e+82) (/ (* (pow z y) x) (* (* (exp b) y) a)) t_1))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (pow(a, (t - 1.0)) * x) / y;
double tmp;
if (t <= -55000000000000.0) {
tmp = t_1;
} else if (t <= 8e+82) {
tmp = (pow(z, y) * x) / ((exp(b) * y) * 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 = ((a ** (t - 1.0d0)) * x) / y
if (t <= (-55000000000000.0d0)) then
tmp = t_1
else if (t <= 8d+82) then
tmp = ((z ** y) * x) / ((exp(b) * y) * 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 = (Math.pow(a, (t - 1.0)) * x) / y;
double tmp;
if (t <= -55000000000000.0) {
tmp = t_1;
} else if (t <= 8e+82) {
tmp = (Math.pow(z, y) * x) / ((Math.exp(b) * y) * a);
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = (math.pow(a, (t - 1.0)) * x) / y tmp = 0 if t <= -55000000000000.0: tmp = t_1 elif t <= 8e+82: tmp = (math.pow(z, y) * x) / ((math.exp(b) * y) * a) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(Float64((a ^ Float64(t - 1.0)) * x) / y) tmp = 0.0 if (t <= -55000000000000.0) tmp = t_1; elseif (t <= 8e+82) tmp = Float64(Float64((z ^ y) * x) / Float64(Float64(exp(b) * y) * a)); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = ((a ^ (t - 1.0)) * x) / y; tmp = 0.0; if (t <= -55000000000000.0) tmp = t_1; elseif (t <= 8e+82) tmp = ((z ^ y) * x) / ((exp(b) * y) * a); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(N[Power[a, N[(t - 1.0), $MachinePrecision]], $MachinePrecision] * x), $MachinePrecision] / y), $MachinePrecision]}, If[LessEqual[t, -55000000000000.0], t$95$1, If[LessEqual[t, 8e+82], N[(N[(N[Power[z, y], $MachinePrecision] * x), $MachinePrecision] / N[(N[(N[Exp[b], $MachinePrecision] * y), $MachinePrecision] * a), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{{a}^{\left(t - 1\right)} \cdot x}{y}\\
\mathbf{if}\;t \leq -55000000000000:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t \leq 8 \cdot 10^{+82}:\\
\;\;\;\;\frac{{z}^{y} \cdot x}{\left(e^{b} \cdot y\right) \cdot a}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if t < -5.5e13 or 7.9999999999999997e82 < t Initial program 100.0%
Taylor expanded in b around 0
exp-sumN/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
*-commutativeN/A
exp-to-powN/A
lower-pow.f64N/A
exp-prodN/A
lower-pow.f64N/A
rem-exp-logN/A
lower--.f6473.2
Applied rewrites73.2%
Taylor expanded in y around 0
Applied rewrites85.4%
if -5.5e13 < t < 7.9999999999999997e82Initial program 95.6%
lift-/.f64N/A
lift-*.f64N/A
*-commutativeN/A
associate-/l*N/A
lift-exp.f64N/A
lift--.f64N/A
exp-diffN/A
clear-numN/A
frac-timesN/A
lower-/.f64N/A
lower-*.f64N/A
lower-*.f64N/A
Applied rewrites87.9%
Taylor expanded in t around 0
*-commutativeN/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-exp.f64N/A
lower-/.f64N/A
lower-pow.f6491.3
Applied rewrites91.3%
Applied rewrites90.6%
(FPCore (x y z t a b)
:precision binary64
(if (<= b -3.6e-10)
(/ (/ x (* (exp b) y)) a)
(if (<= b 6500000000.0)
(/ (* (* (pow z y) x) (pow a (- t 1.0))) y)
(* (/ (exp (- b)) y) x))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (b <= -3.6e-10) {
tmp = (x / (exp(b) * y)) / a;
} else if (b <= 6500000000.0) {
tmp = ((pow(z, y) * x) * pow(a, (t - 1.0))) / y;
} else {
tmp = (exp(-b) / y) * x;
}
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 <= (-3.6d-10)) then
tmp = (x / (exp(b) * y)) / a
else if (b <= 6500000000.0d0) then
tmp = (((z ** y) * x) * (a ** (t - 1.0d0))) / y
else
tmp = (exp(-b) / y) * x
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 <= -3.6e-10) {
tmp = (x / (Math.exp(b) * y)) / a;
} else if (b <= 6500000000.0) {
tmp = ((Math.pow(z, y) * x) * Math.pow(a, (t - 1.0))) / y;
} else {
tmp = (Math.exp(-b) / y) * x;
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if b <= -3.6e-10: tmp = (x / (math.exp(b) * y)) / a elif b <= 6500000000.0: tmp = ((math.pow(z, y) * x) * math.pow(a, (t - 1.0))) / y else: tmp = (math.exp(-b) / y) * x return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (b <= -3.6e-10) tmp = Float64(Float64(x / Float64(exp(b) * y)) / a); elseif (b <= 6500000000.0) tmp = Float64(Float64(Float64((z ^ y) * x) * (a ^ Float64(t - 1.0))) / y); else tmp = Float64(Float64(exp(Float64(-b)) / y) * x); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (b <= -3.6e-10) tmp = (x / (exp(b) * y)) / a; elseif (b <= 6500000000.0) tmp = (((z ^ y) * x) * (a ^ (t - 1.0))) / y; else tmp = (exp(-b) / y) * x; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[b, -3.6e-10], N[(N[(x / N[(N[Exp[b], $MachinePrecision] * y), $MachinePrecision]), $MachinePrecision] / a), $MachinePrecision], If[LessEqual[b, 6500000000.0], N[(N[(N[(N[Power[z, y], $MachinePrecision] * x), $MachinePrecision] * N[Power[a, N[(t - 1.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision], N[(N[(N[Exp[(-b)], $MachinePrecision] / y), $MachinePrecision] * x), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq -3.6 \cdot 10^{-10}:\\
\;\;\;\;\frac{\frac{x}{e^{b} \cdot y}}{a}\\
\mathbf{elif}\;b \leq 6500000000:\\
\;\;\;\;\frac{\left({z}^{y} \cdot x\right) \cdot {a}^{\left(t - 1\right)}}{y}\\
\mathbf{else}:\\
\;\;\;\;\frac{e^{-b}}{y} \cdot x\\
\end{array}
\end{array}
if b < -3.6e-10Initial program 96.6%
lift-/.f64N/A
lift-*.f64N/A
*-commutativeN/A
associate-/l*N/A
lift-exp.f64N/A
lift--.f64N/A
exp-diffN/A
clear-numN/A
frac-timesN/A
lower-/.f64N/A
lower-*.f64N/A
lower-*.f64N/A
Applied rewrites58.2%
Taylor expanded in t around 0
*-commutativeN/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-exp.f64N/A
lower-/.f64N/A
lower-pow.f6467.5
Applied rewrites67.5%
Taylor expanded in y around 0
Applied rewrites78.5%
if -3.6e-10 < b < 6.5e9Initial program 96.6%
Taylor expanded in b around 0
exp-sumN/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
*-commutativeN/A
exp-to-powN/A
lower-pow.f64N/A
exp-prodN/A
lower-pow.f64N/A
rem-exp-logN/A
lower--.f6488.9
Applied rewrites88.9%
if 6.5e9 < b Initial program 100.0%
Taylor expanded in b around inf
mul-1-negN/A
lower-neg.f6481.8
Applied rewrites81.8%
lift-/.f64N/A
lift-*.f64N/A
associate-/l*N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f6481.8
Applied rewrites81.8%
Final simplification84.8%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (* (/ (exp (- b)) y) x)))
(if (<= b -880000000.0)
t_1
(if (<= b 6500000000.0)
(* (/ (pow a (- t 1.0)) y) (* (pow z y) x))
t_1))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (exp(-b) / y) * x;
double tmp;
if (b <= -880000000.0) {
tmp = t_1;
} else if (b <= 6500000000.0) {
tmp = (pow(a, (t - 1.0)) / y) * (pow(z, y) * x);
} 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 = (exp(-b) / y) * x
if (b <= (-880000000.0d0)) then
tmp = t_1
else if (b <= 6500000000.0d0) then
tmp = ((a ** (t - 1.0d0)) / y) * ((z ** y) * x)
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 = (Math.exp(-b) / y) * x;
double tmp;
if (b <= -880000000.0) {
tmp = t_1;
} else if (b <= 6500000000.0) {
tmp = (Math.pow(a, (t - 1.0)) / y) * (Math.pow(z, y) * x);
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = (math.exp(-b) / y) * x tmp = 0 if b <= -880000000.0: tmp = t_1 elif b <= 6500000000.0: tmp = (math.pow(a, (t - 1.0)) / y) * (math.pow(z, y) * x) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(Float64(exp(Float64(-b)) / y) * x) tmp = 0.0 if (b <= -880000000.0) tmp = t_1; elseif (b <= 6500000000.0) tmp = Float64(Float64((a ^ Float64(t - 1.0)) / y) * Float64((z ^ y) * x)); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = (exp(-b) / y) * x; tmp = 0.0; if (b <= -880000000.0) tmp = t_1; elseif (b <= 6500000000.0) tmp = ((a ^ (t - 1.0)) / y) * ((z ^ y) * x); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(N[Exp[(-b)], $MachinePrecision] / y), $MachinePrecision] * x), $MachinePrecision]}, If[LessEqual[b, -880000000.0], t$95$1, If[LessEqual[b, 6500000000.0], N[(N[(N[Power[a, N[(t - 1.0), $MachinePrecision]], $MachinePrecision] / y), $MachinePrecision] * N[(N[Power[z, y], $MachinePrecision] * x), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{e^{-b}}{y} \cdot x\\
\mathbf{if}\;b \leq -880000000:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;b \leq 6500000000:\\
\;\;\;\;\frac{{a}^{\left(t - 1\right)}}{y} \cdot \left({z}^{y} \cdot x\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if b < -8.8e8 or 6.5e9 < b Initial program 100.0%
Taylor expanded in b around inf
mul-1-negN/A
lower-neg.f6480.5
Applied rewrites80.5%
lift-/.f64N/A
lift-*.f64N/A
associate-/l*N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f6480.5
Applied rewrites80.5%
if -8.8e8 < b < 6.5e9Initial program 95.3%
Taylor expanded in b around 0
exp-sumN/A
associate-*r*N/A
associate-/l*N/A
lower-*.f64N/A
lower-*.f64N/A
*-commutativeN/A
exp-to-powN/A
lower-pow.f64N/A
lower-/.f64N/A
exp-prodN/A
lower-pow.f64N/A
rem-exp-logN/A
lower--.f6487.9
Applied rewrites87.9%
Final simplification84.5%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (/ (/ (* (pow z y) x) a) y)))
(if (<= y -6.5e+51)
t_1
(if (<= y -1.42e-260)
(/ (/ x (* (exp b) y)) a)
(if (<= y 6.5e+35)
(* (/ (* 1.0 (pow a (- t 1.0))) (* 1.0 y)) x)
t_1)))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = ((pow(z, y) * x) / a) / y;
double tmp;
if (y <= -6.5e+51) {
tmp = t_1;
} else if (y <= -1.42e-260) {
tmp = (x / (exp(b) * y)) / a;
} else if (y <= 6.5e+35) {
tmp = ((1.0 * pow(a, (t - 1.0))) / (1.0 * y)) * x;
} 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 = (((z ** y) * x) / a) / y
if (y <= (-6.5d+51)) then
tmp = t_1
else if (y <= (-1.42d-260)) then
tmp = (x / (exp(b) * y)) / a
else if (y <= 6.5d+35) then
tmp = ((1.0d0 * (a ** (t - 1.0d0))) / (1.0d0 * y)) * x
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 = ((Math.pow(z, y) * x) / a) / y;
double tmp;
if (y <= -6.5e+51) {
tmp = t_1;
} else if (y <= -1.42e-260) {
tmp = (x / (Math.exp(b) * y)) / a;
} else if (y <= 6.5e+35) {
tmp = ((1.0 * Math.pow(a, (t - 1.0))) / (1.0 * y)) * x;
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = ((math.pow(z, y) * x) / a) / y tmp = 0 if y <= -6.5e+51: tmp = t_1 elif y <= -1.42e-260: tmp = (x / (math.exp(b) * y)) / a elif y <= 6.5e+35: tmp = ((1.0 * math.pow(a, (t - 1.0))) / (1.0 * y)) * x else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(Float64(Float64((z ^ y) * x) / a) / y) tmp = 0.0 if (y <= -6.5e+51) tmp = t_1; elseif (y <= -1.42e-260) tmp = Float64(Float64(x / Float64(exp(b) * y)) / a); elseif (y <= 6.5e+35) tmp = Float64(Float64(Float64(1.0 * (a ^ Float64(t - 1.0))) / Float64(1.0 * y)) * x); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = (((z ^ y) * x) / a) / y; tmp = 0.0; if (y <= -6.5e+51) tmp = t_1; elseif (y <= -1.42e-260) tmp = (x / (exp(b) * y)) / a; elseif (y <= 6.5e+35) tmp = ((1.0 * (a ^ (t - 1.0))) / (1.0 * y)) * x; else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(N[(N[Power[z, y], $MachinePrecision] * x), $MachinePrecision] / a), $MachinePrecision] / y), $MachinePrecision]}, If[LessEqual[y, -6.5e+51], t$95$1, If[LessEqual[y, -1.42e-260], N[(N[(x / N[(N[Exp[b], $MachinePrecision] * y), $MachinePrecision]), $MachinePrecision] / a), $MachinePrecision], If[LessEqual[y, 6.5e+35], N[(N[(N[(1.0 * N[Power[a, N[(t - 1.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / N[(1.0 * y), $MachinePrecision]), $MachinePrecision] * x), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{\frac{{z}^{y} \cdot x}{a}}{y}\\
\mathbf{if}\;y \leq -6.5 \cdot 10^{+51}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq -1.42 \cdot 10^{-260}:\\
\;\;\;\;\frac{\frac{x}{e^{b} \cdot y}}{a}\\
\mathbf{elif}\;y \leq 6.5 \cdot 10^{+35}:\\
\;\;\;\;\frac{1 \cdot {a}^{\left(t - 1\right)}}{1 \cdot y} \cdot x\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -6.5e51 or 6.5000000000000003e35 < y Initial program 100.0%
Taylor expanded in b around 0
exp-sumN/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
*-commutativeN/A
exp-to-powN/A
lower-pow.f64N/A
exp-prodN/A
lower-pow.f64N/A
rem-exp-logN/A
lower--.f6474.1
Applied rewrites74.1%
Taylor expanded in t around 0
Applied rewrites89.1%
if -6.5e51 < y < -1.42000000000000004e-260Initial program 95.7%
lift-/.f64N/A
lift-*.f64N/A
*-commutativeN/A
associate-/l*N/A
lift-exp.f64N/A
lift--.f64N/A
exp-diffN/A
clear-numN/A
frac-timesN/A
lower-/.f64N/A
lower-*.f64N/A
lower-*.f64N/A
Applied rewrites82.0%
Taylor expanded in t around 0
*-commutativeN/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-exp.f64N/A
lower-/.f64N/A
lower-pow.f6479.0
Applied rewrites79.0%
Taylor expanded in y around 0
Applied rewrites83.0%
if -1.42000000000000004e-260 < y < 6.5000000000000003e35Initial program 95.5%
lift-/.f64N/A
lift-*.f64N/A
associate-/l*N/A
*-commutativeN/A
lower-*.f64N/A
Applied rewrites83.7%
Taylor expanded in y around 0
Applied rewrites84.8%
Taylor expanded in b around 0
Applied rewrites75.5%
Final simplification82.7%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (/ (* (pow z y) x) (* a y))) (t_2 (* (/ (exp (- b)) y) x)))
(if (<= b -6.6e-7)
t_2
(if (<= b -1.6e-308)
t_1
(if (<= b 4e-55)
(/ (* (pow a (- t 1.0)) x) y)
(if (<= b 6200000000.0) t_1 t_2))))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (pow(z, y) * x) / (a * y);
double t_2 = (exp(-b) / y) * x;
double tmp;
if (b <= -6.6e-7) {
tmp = t_2;
} else if (b <= -1.6e-308) {
tmp = t_1;
} else if (b <= 4e-55) {
tmp = (pow(a, (t - 1.0)) * x) / y;
} else if (b <= 6200000000.0) {
tmp = 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 = ((z ** y) * x) / (a * y)
t_2 = (exp(-b) / y) * x
if (b <= (-6.6d-7)) then
tmp = t_2
else if (b <= (-1.6d-308)) then
tmp = t_1
else if (b <= 4d-55) then
tmp = ((a ** (t - 1.0d0)) * x) / y
else if (b <= 6200000000.0d0) then
tmp = 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 = (Math.pow(z, y) * x) / (a * y);
double t_2 = (Math.exp(-b) / y) * x;
double tmp;
if (b <= -6.6e-7) {
tmp = t_2;
} else if (b <= -1.6e-308) {
tmp = t_1;
} else if (b <= 4e-55) {
tmp = (Math.pow(a, (t - 1.0)) * x) / y;
} else if (b <= 6200000000.0) {
tmp = t_1;
} else {
tmp = t_2;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = (math.pow(z, y) * x) / (a * y) t_2 = (math.exp(-b) / y) * x tmp = 0 if b <= -6.6e-7: tmp = t_2 elif b <= -1.6e-308: tmp = t_1 elif b <= 4e-55: tmp = (math.pow(a, (t - 1.0)) * x) / y elif b <= 6200000000.0: tmp = t_1 else: tmp = t_2 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(Float64((z ^ y) * x) / Float64(a * y)) t_2 = Float64(Float64(exp(Float64(-b)) / y) * x) tmp = 0.0 if (b <= -6.6e-7) tmp = t_2; elseif (b <= -1.6e-308) tmp = t_1; elseif (b <= 4e-55) tmp = Float64(Float64((a ^ Float64(t - 1.0)) * x) / y); elseif (b <= 6200000000.0) tmp = t_1; else tmp = t_2; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = ((z ^ y) * x) / (a * y); t_2 = (exp(-b) / y) * x; tmp = 0.0; if (b <= -6.6e-7) tmp = t_2; elseif (b <= -1.6e-308) tmp = t_1; elseif (b <= 4e-55) tmp = ((a ^ (t - 1.0)) * x) / y; elseif (b <= 6200000000.0) tmp = t_1; else tmp = t_2; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(N[Power[z, y], $MachinePrecision] * x), $MachinePrecision] / N[(a * y), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(N[(N[Exp[(-b)], $MachinePrecision] / y), $MachinePrecision] * x), $MachinePrecision]}, If[LessEqual[b, -6.6e-7], t$95$2, If[LessEqual[b, -1.6e-308], t$95$1, If[LessEqual[b, 4e-55], N[(N[(N[Power[a, N[(t - 1.0), $MachinePrecision]], $MachinePrecision] * x), $MachinePrecision] / y), $MachinePrecision], If[LessEqual[b, 6200000000.0], t$95$1, t$95$2]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{{z}^{y} \cdot x}{a \cdot y}\\
t_2 := \frac{e^{-b}}{y} \cdot x\\
\mathbf{if}\;b \leq -6.6 \cdot 10^{-7}:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;b \leq -1.6 \cdot 10^{-308}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;b \leq 4 \cdot 10^{-55}:\\
\;\;\;\;\frac{{a}^{\left(t - 1\right)} \cdot x}{y}\\
\mathbf{elif}\;b \leq 6200000000:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;t\_2\\
\end{array}
\end{array}
if b < -6.6000000000000003e-7 or 6.2e9 < b Initial program 100.0%
Taylor expanded in b around inf
mul-1-negN/A
lower-neg.f6480.0
Applied rewrites80.0%
lift-/.f64N/A
lift-*.f64N/A
associate-/l*N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f6480.0
Applied rewrites80.0%
if -6.6000000000000003e-7 < b < -1.6000000000000001e-308 or 3.99999999999999998e-55 < b < 6.2e9Initial program 94.1%
lift-/.f64N/A
lift-*.f64N/A
*-commutativeN/A
associate-/l*N/A
lift-exp.f64N/A
lift--.f64N/A
exp-diffN/A
clear-numN/A
frac-timesN/A
lower-/.f64N/A
lower-*.f64N/A
lower-*.f64N/A
Applied rewrites87.9%
Taylor expanded in t around 0
*-commutativeN/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-exp.f64N/A
lower-/.f64N/A
lower-pow.f6477.6
Applied rewrites77.6%
Applied rewrites77.6%
Taylor expanded in b around 0
Applied rewrites77.6%
if -1.6000000000000001e-308 < b < 3.99999999999999998e-55Initial program 96.9%
Taylor expanded in b around 0
exp-sumN/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
*-commutativeN/A
exp-to-powN/A
lower-pow.f64N/A
exp-prodN/A
lower-pow.f64N/A
rem-exp-logN/A
lower--.f6490.7
Applied rewrites90.7%
Taylor expanded in y around 0
Applied rewrites84.0%
Final simplification80.1%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (/ (/ (* (pow z y) x) a) y)))
(if (<= y -6.5e+51)
t_1
(if (<= y -2.35e-260)
(/ (/ x (* (exp b) y)) a)
(if (<= y 5.2e+35) (/ (* (pow a (- t 1.0)) x) y) t_1)))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = ((pow(z, y) * x) / a) / y;
double tmp;
if (y <= -6.5e+51) {
tmp = t_1;
} else if (y <= -2.35e-260) {
tmp = (x / (exp(b) * y)) / a;
} else if (y <= 5.2e+35) {
tmp = (pow(a, (t - 1.0)) * x) / 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 = (((z ** y) * x) / a) / y
if (y <= (-6.5d+51)) then
tmp = t_1
else if (y <= (-2.35d-260)) then
tmp = (x / (exp(b) * y)) / a
else if (y <= 5.2d+35) then
tmp = ((a ** (t - 1.0d0)) * x) / 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 = ((Math.pow(z, y) * x) / a) / y;
double tmp;
if (y <= -6.5e+51) {
tmp = t_1;
} else if (y <= -2.35e-260) {
tmp = (x / (Math.exp(b) * y)) / a;
} else if (y <= 5.2e+35) {
tmp = (Math.pow(a, (t - 1.0)) * x) / y;
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = ((math.pow(z, y) * x) / a) / y tmp = 0 if y <= -6.5e+51: tmp = t_1 elif y <= -2.35e-260: tmp = (x / (math.exp(b) * y)) / a elif y <= 5.2e+35: tmp = (math.pow(a, (t - 1.0)) * x) / y else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(Float64(Float64((z ^ y) * x) / a) / y) tmp = 0.0 if (y <= -6.5e+51) tmp = t_1; elseif (y <= -2.35e-260) tmp = Float64(Float64(x / Float64(exp(b) * y)) / a); elseif (y <= 5.2e+35) tmp = Float64(Float64((a ^ Float64(t - 1.0)) * x) / y); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = (((z ^ y) * x) / a) / y; tmp = 0.0; if (y <= -6.5e+51) tmp = t_1; elseif (y <= -2.35e-260) tmp = (x / (exp(b) * y)) / a; elseif (y <= 5.2e+35) tmp = ((a ^ (t - 1.0)) * x) / y; else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(N[(N[Power[z, y], $MachinePrecision] * x), $MachinePrecision] / a), $MachinePrecision] / y), $MachinePrecision]}, If[LessEqual[y, -6.5e+51], t$95$1, If[LessEqual[y, -2.35e-260], N[(N[(x / N[(N[Exp[b], $MachinePrecision] * y), $MachinePrecision]), $MachinePrecision] / a), $MachinePrecision], If[LessEqual[y, 5.2e+35], N[(N[(N[Power[a, N[(t - 1.0), $MachinePrecision]], $MachinePrecision] * x), $MachinePrecision] / y), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{\frac{{z}^{y} \cdot x}{a}}{y}\\
\mathbf{if}\;y \leq -6.5 \cdot 10^{+51}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq -2.35 \cdot 10^{-260}:\\
\;\;\;\;\frac{\frac{x}{e^{b} \cdot y}}{a}\\
\mathbf{elif}\;y \leq 5.2 \cdot 10^{+35}:\\
\;\;\;\;\frac{{a}^{\left(t - 1\right)} \cdot x}{y}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -6.5e51 or 5.20000000000000013e35 < y Initial program 100.0%
Taylor expanded in b around 0
exp-sumN/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
*-commutativeN/A
exp-to-powN/A
lower-pow.f64N/A
exp-prodN/A
lower-pow.f64N/A
rem-exp-logN/A
lower--.f6474.1
Applied rewrites74.1%
Taylor expanded in t around 0
Applied rewrites89.1%
if -6.5e51 < y < -2.35e-260Initial program 95.7%
lift-/.f64N/A
lift-*.f64N/A
*-commutativeN/A
associate-/l*N/A
lift-exp.f64N/A
lift--.f64N/A
exp-diffN/A
clear-numN/A
frac-timesN/A
lower-/.f64N/A
lower-*.f64N/A
lower-*.f64N/A
Applied rewrites82.0%
Taylor expanded in t around 0
*-commutativeN/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-exp.f64N/A
lower-/.f64N/A
lower-pow.f6479.0
Applied rewrites79.0%
Taylor expanded in y around 0
Applied rewrites83.0%
if -2.35e-260 < y < 5.20000000000000013e35Initial program 95.5%
Taylor expanded in b around 0
exp-sumN/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
*-commutativeN/A
exp-to-powN/A
lower-pow.f64N/A
exp-prodN/A
lower-pow.f64N/A
rem-exp-logN/A
lower--.f6474.5
Applied rewrites74.5%
Taylor expanded in y around 0
Applied rewrites74.5%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (/ (* (pow a (- t 1.0)) x) y)))
(if (<= t -105000.0)
t_1
(if (<= t 3.9e+74) (/ (/ x (* (exp b) y)) a) t_1))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (pow(a, (t - 1.0)) * x) / y;
double tmp;
if (t <= -105000.0) {
tmp = t_1;
} else if (t <= 3.9e+74) {
tmp = (x / (exp(b) * y)) / 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 = ((a ** (t - 1.0d0)) * x) / y
if (t <= (-105000.0d0)) then
tmp = t_1
else if (t <= 3.9d+74) then
tmp = (x / (exp(b) * y)) / 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 = (Math.pow(a, (t - 1.0)) * x) / y;
double tmp;
if (t <= -105000.0) {
tmp = t_1;
} else if (t <= 3.9e+74) {
tmp = (x / (Math.exp(b) * y)) / a;
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = (math.pow(a, (t - 1.0)) * x) / y tmp = 0 if t <= -105000.0: tmp = t_1 elif t <= 3.9e+74: tmp = (x / (math.exp(b) * y)) / a else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(Float64((a ^ Float64(t - 1.0)) * x) / y) tmp = 0.0 if (t <= -105000.0) tmp = t_1; elseif (t <= 3.9e+74) tmp = Float64(Float64(x / Float64(exp(b) * y)) / a); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = ((a ^ (t - 1.0)) * x) / y; tmp = 0.0; if (t <= -105000.0) tmp = t_1; elseif (t <= 3.9e+74) tmp = (x / (exp(b) * y)) / a; else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(N[Power[a, N[(t - 1.0), $MachinePrecision]], $MachinePrecision] * x), $MachinePrecision] / y), $MachinePrecision]}, If[LessEqual[t, -105000.0], t$95$1, If[LessEqual[t, 3.9e+74], N[(N[(x / N[(N[Exp[b], $MachinePrecision] * y), $MachinePrecision]), $MachinePrecision] / a), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{{a}^{\left(t - 1\right)} \cdot x}{y}\\
\mathbf{if}\;t \leq -105000:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t \leq 3.9 \cdot 10^{+74}:\\
\;\;\;\;\frac{\frac{x}{e^{b} \cdot y}}{a}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if t < -105000 or 3.90000000000000008e74 < t Initial program 100.0%
Taylor expanded in b around 0
exp-sumN/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
*-commutativeN/A
exp-to-powN/A
lower-pow.f64N/A
exp-prodN/A
lower-pow.f64N/A
rem-exp-logN/A
lower--.f6473.5
Applied rewrites73.5%
Taylor expanded in y around 0
Applied rewrites85.2%
if -105000 < t < 3.90000000000000008e74Initial program 95.4%
lift-/.f64N/A
lift-*.f64N/A
*-commutativeN/A
associate-/l*N/A
lift-exp.f64N/A
lift--.f64N/A
exp-diffN/A
clear-numN/A
frac-timesN/A
lower-/.f64N/A
lower-*.f64N/A
lower-*.f64N/A
Applied rewrites88.2%
Taylor expanded in t around 0
*-commutativeN/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-exp.f64N/A
lower-/.f64N/A
lower-pow.f6491.7
Applied rewrites91.7%
Taylor expanded in y around 0
Applied rewrites75.7%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (* (/ (exp (- b)) y) x)))
(if (<= b -1.3e+21)
t_1
(if (<= b 8000000000.0) (/ (* (pow a (- t 1.0)) x) y) t_1))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (exp(-b) / y) * x;
double tmp;
if (b <= -1.3e+21) {
tmp = t_1;
} else if (b <= 8000000000.0) {
tmp = (pow(a, (t - 1.0)) * x) / 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 = (exp(-b) / y) * x
if (b <= (-1.3d+21)) then
tmp = t_1
else if (b <= 8000000000.0d0) then
tmp = ((a ** (t - 1.0d0)) * x) / 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 = (Math.exp(-b) / y) * x;
double tmp;
if (b <= -1.3e+21) {
tmp = t_1;
} else if (b <= 8000000000.0) {
tmp = (Math.pow(a, (t - 1.0)) * x) / y;
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = (math.exp(-b) / y) * x tmp = 0 if b <= -1.3e+21: tmp = t_1 elif b <= 8000000000.0: tmp = (math.pow(a, (t - 1.0)) * x) / y else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(Float64(exp(Float64(-b)) / y) * x) tmp = 0.0 if (b <= -1.3e+21) tmp = t_1; elseif (b <= 8000000000.0) tmp = Float64(Float64((a ^ Float64(t - 1.0)) * x) / y); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = (exp(-b) / y) * x; tmp = 0.0; if (b <= -1.3e+21) tmp = t_1; elseif (b <= 8000000000.0) tmp = ((a ^ (t - 1.0)) * x) / y; else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(N[Exp[(-b)], $MachinePrecision] / y), $MachinePrecision] * x), $MachinePrecision]}, If[LessEqual[b, -1.3e+21], t$95$1, If[LessEqual[b, 8000000000.0], N[(N[(N[Power[a, N[(t - 1.0), $MachinePrecision]], $MachinePrecision] * x), $MachinePrecision] / y), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{e^{-b}}{y} \cdot x\\
\mathbf{if}\;b \leq -1.3 \cdot 10^{+21}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;b \leq 8000000000:\\
\;\;\;\;\frac{{a}^{\left(t - 1\right)} \cdot x}{y}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if b < -1.3e21 or 8e9 < b Initial program 100.0%
Taylor expanded in b around inf
mul-1-negN/A
lower-neg.f6481.0
Applied rewrites81.0%
lift-/.f64N/A
lift-*.f64N/A
associate-/l*N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f6481.0
Applied rewrites81.0%
if -1.3e21 < b < 8e9Initial program 95.4%
Taylor expanded in b around 0
exp-sumN/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
*-commutativeN/A
exp-to-powN/A
lower-pow.f64N/A
exp-prodN/A
lower-pow.f64N/A
rem-exp-logN/A
lower--.f6486.6
Applied rewrites86.6%
Taylor expanded in y around 0
Applied rewrites72.6%
(FPCore (x y z t a b) :precision binary64 (let* ((t_1 (* (/ (exp (- b)) y) x))) (if (<= b -6.6e-7) t_1 (if (<= b 4.3e-26) (/ 1.0 (/ y (/ x a))) t_1))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (exp(-b) / y) * x;
double tmp;
if (b <= -6.6e-7) {
tmp = t_1;
} else if (b <= 4.3e-26) {
tmp = 1.0 / (y / (x / 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 = (exp(-b) / y) * x
if (b <= (-6.6d-7)) then
tmp = t_1
else if (b <= 4.3d-26) then
tmp = 1.0d0 / (y / (x / 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 = (Math.exp(-b) / y) * x;
double tmp;
if (b <= -6.6e-7) {
tmp = t_1;
} else if (b <= 4.3e-26) {
tmp = 1.0 / (y / (x / a));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = (math.exp(-b) / y) * x tmp = 0 if b <= -6.6e-7: tmp = t_1 elif b <= 4.3e-26: tmp = 1.0 / (y / (x / a)) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(Float64(exp(Float64(-b)) / y) * x) tmp = 0.0 if (b <= -6.6e-7) tmp = t_1; elseif (b <= 4.3e-26) tmp = Float64(1.0 / Float64(y / Float64(x / a))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = (exp(-b) / y) * x; tmp = 0.0; if (b <= -6.6e-7) tmp = t_1; elseif (b <= 4.3e-26) tmp = 1.0 / (y / (x / a)); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(N[Exp[(-b)], $MachinePrecision] / y), $MachinePrecision] * x), $MachinePrecision]}, If[LessEqual[b, -6.6e-7], t$95$1, If[LessEqual[b, 4.3e-26], N[(1.0 / N[(y / N[(x / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{e^{-b}}{y} \cdot x\\
\mathbf{if}\;b \leq -6.6 \cdot 10^{-7}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;b \leq 4.3 \cdot 10^{-26}:\\
\;\;\;\;\frac{1}{\frac{y}{\frac{x}{a}}}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if b < -6.6000000000000003e-7 or 4.29999999999999988e-26 < b Initial program 99.3%
Taylor expanded in b around inf
mul-1-negN/A
lower-neg.f6476.8
Applied rewrites76.8%
lift-/.f64N/A
lift-*.f64N/A
associate-/l*N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f6476.8
Applied rewrites76.8%
if -6.6000000000000003e-7 < b < 4.29999999999999988e-26Initial program 95.6%
Taylor expanded in b around 0
exp-sumN/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
*-commutativeN/A
exp-to-powN/A
lower-pow.f64N/A
exp-prodN/A
lower-pow.f64N/A
rem-exp-logN/A
lower--.f6489.1
Applied rewrites89.1%
Taylor expanded in t around 0
Applied rewrites76.5%
Taylor expanded in y around 0
Applied rewrites44.4%
lift-/.f64N/A
clear-numN/A
lower-/.f64N/A
lower-/.f6444.4
Applied rewrites44.4%
(FPCore (x y z t a b) :precision binary64 (let* ((t_1 (* (/ x y) (exp (- b))))) (if (<= b -6.6e-7) t_1 (if (<= b 4.3e-26) (/ 1.0 (/ y (/ x a))) t_1))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (x / y) * exp(-b);
double tmp;
if (b <= -6.6e-7) {
tmp = t_1;
} else if (b <= 4.3e-26) {
tmp = 1.0 / (y / (x / 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 / y) * exp(-b)
if (b <= (-6.6d-7)) then
tmp = t_1
else if (b <= 4.3d-26) then
tmp = 1.0d0 / (y / (x / 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 / y) * Math.exp(-b);
double tmp;
if (b <= -6.6e-7) {
tmp = t_1;
} else if (b <= 4.3e-26) {
tmp = 1.0 / (y / (x / a));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = (x / y) * math.exp(-b) tmp = 0 if b <= -6.6e-7: tmp = t_1 elif b <= 4.3e-26: tmp = 1.0 / (y / (x / a)) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(Float64(x / y) * exp(Float64(-b))) tmp = 0.0 if (b <= -6.6e-7) tmp = t_1; elseif (b <= 4.3e-26) tmp = Float64(1.0 / Float64(y / Float64(x / a))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = (x / y) * exp(-b); tmp = 0.0; if (b <= -6.6e-7) tmp = t_1; elseif (b <= 4.3e-26) tmp = 1.0 / (y / (x / a)); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(x / y), $MachinePrecision] * N[Exp[(-b)], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[b, -6.6e-7], t$95$1, If[LessEqual[b, 4.3e-26], N[(1.0 / N[(y / N[(x / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{x}{y} \cdot e^{-b}\\
\mathbf{if}\;b \leq -6.6 \cdot 10^{-7}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;b \leq 4.3 \cdot 10^{-26}:\\
\;\;\;\;\frac{1}{\frac{y}{\frac{x}{a}}}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if b < -6.6000000000000003e-7 or 4.29999999999999988e-26 < b Initial program 99.3%
Taylor expanded in b around inf
mul-1-negN/A
lower-neg.f6476.8
Applied rewrites76.8%
lift-/.f64N/A
lift-*.f64N/A
*-commutativeN/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f6466.6
Applied rewrites66.6%
if -6.6000000000000003e-7 < b < 4.29999999999999988e-26Initial program 95.6%
Taylor expanded in b around 0
exp-sumN/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
*-commutativeN/A
exp-to-powN/A
lower-pow.f64N/A
exp-prodN/A
lower-pow.f64N/A
rem-exp-logN/A
lower--.f6489.1
Applied rewrites89.1%
Taylor expanded in t around 0
Applied rewrites76.5%
Taylor expanded in y around 0
Applied rewrites44.4%
lift-/.f64N/A
clear-numN/A
lower-/.f64N/A
lower-/.f6444.4
Applied rewrites44.4%
Final simplification55.4%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (/ x (* a y))))
(if (<= b -3.5e+178)
(fma (- (fma (* -0.5 t_1) b t_1)) b t_1)
(/ (/ x (fma (fma (* (fma 0.16666666666666666 b 0.5) y) b y) b y)) a))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = x / (a * y);
double tmp;
if (b <= -3.5e+178) {
tmp = fma(-fma((-0.5 * t_1), b, t_1), b, t_1);
} else {
tmp = (x / fma(fma((fma(0.16666666666666666, b, 0.5) * y), b, y), b, y)) / a;
}
return tmp;
}
function code(x, y, z, t, a, b) t_1 = Float64(x / Float64(a * y)) tmp = 0.0 if (b <= -3.5e+178) tmp = fma(Float64(-fma(Float64(-0.5 * t_1), b, t_1)), b, t_1); else tmp = Float64(Float64(x / fma(fma(Float64(fma(0.16666666666666666, b, 0.5) * y), b, y), b, y)) / a); end return tmp end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(x / N[(a * y), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[b, -3.5e+178], N[((-N[(N[(-0.5 * t$95$1), $MachinePrecision] * b + t$95$1), $MachinePrecision]) * b + t$95$1), $MachinePrecision], N[(N[(x / N[(N[(N[(N[(0.16666666666666666 * b + 0.5), $MachinePrecision] * y), $MachinePrecision] * b + y), $MachinePrecision] * b + y), $MachinePrecision]), $MachinePrecision] / a), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{x}{a \cdot y}\\
\mathbf{if}\;b \leq -3.5 \cdot 10^{+178}:\\
\;\;\;\;\mathsf{fma}\left(-\mathsf{fma}\left(-0.5 \cdot t\_1, b, t\_1\right), b, t\_1\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{x}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(0.16666666666666666, b, 0.5\right) \cdot y, b, y\right), b, y\right)}}{a}\\
\end{array}
\end{array}
if b < -3.5e178Initial program 100.0%
lift-/.f64N/A
lift-*.f64N/A
*-commutativeN/A
associate-/l*N/A
lift-exp.f64N/A
lift--.f64N/A
exp-diffN/A
clear-numN/A
frac-timesN/A
lower-/.f64N/A
lower-*.f64N/A
lower-*.f64N/A
Applied rewrites66.7%
Taylor expanded in t around 0
*-commutativeN/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-exp.f64N/A
lower-/.f64N/A
lower-pow.f6480.1
Applied rewrites80.1%
Taylor expanded in y around 0
Applied rewrites90.2%
Taylor expanded in b around 0
Applied rewrites77.5%
if -3.5e178 < b Initial program 97.1%
lift-/.f64N/A
lift-*.f64N/A
*-commutativeN/A
associate-/l*N/A
lift-exp.f64N/A
lift--.f64N/A
exp-diffN/A
clear-numN/A
frac-timesN/A
lower-/.f64N/A
lower-*.f64N/A
lower-*.f64N/A
Applied rewrites77.9%
Taylor expanded in t around 0
*-commutativeN/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-exp.f64N/A
lower-/.f64N/A
lower-pow.f6471.3
Applied rewrites71.3%
Taylor expanded in y around 0
Applied rewrites55.9%
Taylor expanded in b around 0
Applied rewrites44.7%
Final simplification48.6%
(FPCore (x y z t a b) :precision binary64 (if (<= b -3.5e+178) (/ (fma (- (fma (* -0.5 (/ x y)) b (/ x y))) b (/ x y)) a) (/ (/ x (fma (fma (* (fma 0.16666666666666666 b 0.5) y) b y) b y)) a)))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (b <= -3.5e+178) {
tmp = fma(-fma((-0.5 * (x / y)), b, (x / y)), b, (x / y)) / a;
} else {
tmp = (x / fma(fma((fma(0.16666666666666666, b, 0.5) * y), b, y), b, y)) / a;
}
return tmp;
}
function code(x, y, z, t, a, b) tmp = 0.0 if (b <= -3.5e+178) tmp = Float64(fma(Float64(-fma(Float64(-0.5 * Float64(x / y)), b, Float64(x / y))), b, Float64(x / y)) / a); else tmp = Float64(Float64(x / fma(fma(Float64(fma(0.16666666666666666, b, 0.5) * y), b, y), b, y)) / a); end return tmp end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[b, -3.5e+178], N[(N[((-N[(N[(-0.5 * N[(x / y), $MachinePrecision]), $MachinePrecision] * b + N[(x / y), $MachinePrecision]), $MachinePrecision]) * b + N[(x / y), $MachinePrecision]), $MachinePrecision] / a), $MachinePrecision], N[(N[(x / N[(N[(N[(N[(0.16666666666666666 * b + 0.5), $MachinePrecision] * y), $MachinePrecision] * b + y), $MachinePrecision] * b + y), $MachinePrecision]), $MachinePrecision] / a), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq -3.5 \cdot 10^{+178}:\\
\;\;\;\;\frac{\mathsf{fma}\left(-\mathsf{fma}\left(-0.5 \cdot \frac{x}{y}, b, \frac{x}{y}\right), b, \frac{x}{y}\right)}{a}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{x}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(0.16666666666666666, b, 0.5\right) \cdot y, b, y\right), b, y\right)}}{a}\\
\end{array}
\end{array}
if b < -3.5e178Initial program 100.0%
lift-/.f64N/A
lift-*.f64N/A
*-commutativeN/A
associate-/l*N/A
lift-exp.f64N/A
lift--.f64N/A
exp-diffN/A
clear-numN/A
frac-timesN/A
lower-/.f64N/A
lower-*.f64N/A
lower-*.f64N/A
Applied rewrites66.7%
Taylor expanded in t around 0
*-commutativeN/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-exp.f64N/A
lower-/.f64N/A
lower-pow.f6480.1
Applied rewrites80.1%
Taylor expanded in y around 0
Applied rewrites90.2%
Taylor expanded in b around 0
Applied rewrites77.3%
if -3.5e178 < b Initial program 97.1%
lift-/.f64N/A
lift-*.f64N/A
*-commutativeN/A
associate-/l*N/A
lift-exp.f64N/A
lift--.f64N/A
exp-diffN/A
clear-numN/A
frac-timesN/A
lower-/.f64N/A
lower-*.f64N/A
lower-*.f64N/A
Applied rewrites77.9%
Taylor expanded in t around 0
*-commutativeN/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-exp.f64N/A
lower-/.f64N/A
lower-pow.f6471.3
Applied rewrites71.3%
Taylor expanded in y around 0
Applied rewrites55.9%
Taylor expanded in b around 0
Applied rewrites44.7%
Final simplification48.5%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (/ x (* a y))))
(if (<= b -3.5e+178)
(fma (- b) t_1 t_1)
(/ (/ x (fma (* (fma 0.5 b 1.0) y) b y)) a))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = x / (a * y);
double tmp;
if (b <= -3.5e+178) {
tmp = fma(-b, t_1, t_1);
} else {
tmp = (x / fma((fma(0.5, b, 1.0) * y), b, y)) / a;
}
return tmp;
}
function code(x, y, z, t, a, b) t_1 = Float64(x / Float64(a * y)) tmp = 0.0 if (b <= -3.5e+178) tmp = fma(Float64(-b), t_1, t_1); else tmp = Float64(Float64(x / fma(Float64(fma(0.5, b, 1.0) * y), b, y)) / a); end return tmp end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(x / N[(a * y), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[b, -3.5e+178], N[((-b) * t$95$1 + t$95$1), $MachinePrecision], N[(N[(x / N[(N[(N[(0.5 * b + 1.0), $MachinePrecision] * y), $MachinePrecision] * b + y), $MachinePrecision]), $MachinePrecision] / a), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{x}{a \cdot y}\\
\mathbf{if}\;b \leq -3.5 \cdot 10^{+178}:\\
\;\;\;\;\mathsf{fma}\left(-b, t\_1, t\_1\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{x}{\mathsf{fma}\left(\mathsf{fma}\left(0.5, b, 1\right) \cdot y, b, y\right)}}{a}\\
\end{array}
\end{array}
if b < -3.5e178Initial program 100.0%
lift-/.f64N/A
lift-*.f64N/A
*-commutativeN/A
associate-/l*N/A
lift-exp.f64N/A
lift--.f64N/A
exp-diffN/A
clear-numN/A
frac-timesN/A
lower-/.f64N/A
lower-*.f64N/A
lower-*.f64N/A
Applied rewrites66.7%
Taylor expanded in t around 0
*-commutativeN/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-exp.f64N/A
lower-/.f64N/A
lower-pow.f6480.1
Applied rewrites80.1%
Taylor expanded in y around 0
Applied rewrites90.2%
Taylor expanded in b around 0
Applied rewrites55.0%
if -3.5e178 < b Initial program 97.1%
lift-/.f64N/A
lift-*.f64N/A
*-commutativeN/A
associate-/l*N/A
lift-exp.f64N/A
lift--.f64N/A
exp-diffN/A
clear-numN/A
frac-timesN/A
lower-/.f64N/A
lower-*.f64N/A
lower-*.f64N/A
Applied rewrites77.9%
Taylor expanded in t around 0
*-commutativeN/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-exp.f64N/A
lower-/.f64N/A
lower-pow.f6471.3
Applied rewrites71.3%
Taylor expanded in y around 0
Applied rewrites55.9%
Taylor expanded in b around 0
Applied rewrites43.4%
Final simplification44.7%
(FPCore (x y z t a b) :precision binary64 (if (<= b 1.2e-257) (/ 1.0 (/ y (/ x a))) (/ (/ x (fma b y y)) a)))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (b <= 1.2e-257) {
tmp = 1.0 / (y / (x / a));
} else {
tmp = (x / fma(b, y, y)) / a;
}
return tmp;
}
function code(x, y, z, t, a, b) tmp = 0.0 if (b <= 1.2e-257) tmp = Float64(1.0 / Float64(y / Float64(x / a))); else tmp = Float64(Float64(x / fma(b, y, y)) / a); end return tmp end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[b, 1.2e-257], N[(1.0 / N[(y / N[(x / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(x / N[(b * y + y), $MachinePrecision]), $MachinePrecision] / a), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq 1.2 \cdot 10^{-257}:\\
\;\;\;\;\frac{1}{\frac{y}{\frac{x}{a}}}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{x}{\mathsf{fma}\left(b, y, y\right)}}{a}\\
\end{array}
\end{array}
if b < 1.20000000000000008e-257Initial program 97.1%
Taylor expanded in b around 0
exp-sumN/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
*-commutativeN/A
exp-to-powN/A
lower-pow.f64N/A
exp-prodN/A
lower-pow.f64N/A
rem-exp-logN/A
lower--.f6474.6
Applied rewrites74.6%
Taylor expanded in t around 0
Applied rewrites68.5%
Taylor expanded in y around 0
Applied rewrites40.1%
lift-/.f64N/A
clear-numN/A
lower-/.f64N/A
lower-/.f6440.1
Applied rewrites40.1%
if 1.20000000000000008e-257 < b Initial program 97.8%
lift-/.f64N/A
lift-*.f64N/A
*-commutativeN/A
associate-/l*N/A
lift-exp.f64N/A
lift--.f64N/A
exp-diffN/A
clear-numN/A
frac-timesN/A
lower-/.f64N/A
lower-*.f64N/A
lower-*.f64N/A
Applied rewrites75.9%
Taylor expanded in t around 0
*-commutativeN/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-exp.f64N/A
lower-/.f64N/A
lower-pow.f6472.0
Applied rewrites72.0%
Taylor expanded in y around 0
Applied rewrites64.4%
Taylor expanded in b around 0
Applied rewrites41.1%
(FPCore (x y z t a b) :precision binary64 (if (<= b 1e-265) (/ (/ x a) y) (/ (/ x (fma b y y)) a)))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (b <= 1e-265) {
tmp = (x / a) / y;
} else {
tmp = (x / fma(b, y, y)) / a;
}
return tmp;
}
function code(x, y, z, t, a, b) tmp = 0.0 if (b <= 1e-265) tmp = Float64(Float64(x / a) / y); else tmp = Float64(Float64(x / fma(b, y, y)) / a); end return tmp end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[b, 1e-265], N[(N[(x / a), $MachinePrecision] / y), $MachinePrecision], N[(N[(x / N[(b * y + y), $MachinePrecision]), $MachinePrecision] / a), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq 10^{-265}:\\
\;\;\;\;\frac{\frac{x}{a}}{y}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{x}{\mathsf{fma}\left(b, y, y\right)}}{a}\\
\end{array}
\end{array}
if b < 9.99999999999999985e-266Initial program 97.1%
Taylor expanded in b around 0
exp-sumN/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
*-commutativeN/A
exp-to-powN/A
lower-pow.f64N/A
exp-prodN/A
lower-pow.f64N/A
rem-exp-logN/A
lower--.f6474.2
Applied rewrites74.2%
Taylor expanded in t around 0
Applied rewrites68.0%
Taylor expanded in y around 0
Applied rewrites39.9%
if 9.99999999999999985e-266 < b Initial program 97.8%
lift-/.f64N/A
lift-*.f64N/A
*-commutativeN/A
associate-/l*N/A
lift-exp.f64N/A
lift--.f64N/A
exp-diffN/A
clear-numN/A
frac-timesN/A
lower-/.f64N/A
lower-*.f64N/A
lower-*.f64N/A
Applied rewrites76.2%
Taylor expanded in t around 0
*-commutativeN/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-exp.f64N/A
lower-/.f64N/A
lower-pow.f6472.5
Applied rewrites72.5%
Taylor expanded in y around 0
Applied rewrites64.2%
Taylor expanded in b around 0
Applied rewrites41.3%
(FPCore (x y z t a b) :precision binary64 (if (<= y -6.4e-171) (/ (/ x y) a) (/ x (* a y))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (y <= -6.4e-171) {
tmp = (x / y) / a;
} else {
tmp = x / (a * 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) :: tmp
if (y <= (-6.4d-171)) then
tmp = (x / y) / a
else
tmp = x / (a * y)
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 (y <= -6.4e-171) {
tmp = (x / y) / a;
} else {
tmp = x / (a * y);
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if y <= -6.4e-171: tmp = (x / y) / a else: tmp = x / (a * y) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (y <= -6.4e-171) tmp = Float64(Float64(x / y) / a); else tmp = Float64(x / Float64(a * y)); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (y <= -6.4e-171) tmp = (x / y) / a; else tmp = x / (a * y); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[y, -6.4e-171], N[(N[(x / y), $MachinePrecision] / a), $MachinePrecision], N[(x / N[(a * y), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -6.4 \cdot 10^{-171}:\\
\;\;\;\;\frac{\frac{x}{y}}{a}\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{a \cdot y}\\
\end{array}
\end{array}
if y < -6.4000000000000003e-171Initial program 97.7%
lift-/.f64N/A
lift-*.f64N/A
*-commutativeN/A
associate-/l*N/A
lift-exp.f64N/A
lift--.f64N/A
exp-diffN/A
clear-numN/A
frac-timesN/A
lower-/.f64N/A
lower-*.f64N/A
lower-*.f64N/A
Applied rewrites74.7%
Taylor expanded in t around 0
*-commutativeN/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-exp.f64N/A
lower-/.f64N/A
lower-pow.f6477.4
Applied rewrites77.4%
Taylor expanded in y around 0
Applied rewrites56.9%
Taylor expanded in b around 0
Applied rewrites28.7%
if -6.4000000000000003e-171 < y Initial program 97.3%
lift-/.f64N/A
lift-*.f64N/A
*-commutativeN/A
associate-/l*N/A
lift-exp.f64N/A
lift--.f64N/A
exp-diffN/A
clear-numN/A
frac-timesN/A
lower-/.f64N/A
lower-*.f64N/A
lower-*.f64N/A
Applied rewrites77.6%
Taylor expanded in t around 0
*-commutativeN/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-exp.f64N/A
lower-/.f64N/A
lower-pow.f6469.8
Applied rewrites69.8%
Taylor expanded in y around 0
Applied rewrites61.5%
Taylor expanded in b around 0
Applied rewrites38.8%
Final simplification35.3%
(FPCore (x y z t a b) :precision binary64 (/ (/ x a) y))
double code(double x, double y, double z, double t, double a, double b) {
return (x / a) / 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 / a) / y
end function
public static double code(double x, double y, double z, double t, double a, double b) {
return (x / a) / y;
}
def code(x, y, z, t, a, b): return (x / a) / y
function code(x, y, z, t, a, b) return Float64(Float64(x / a) / y) end
function tmp = code(x, y, z, t, a, b) tmp = (x / a) / y; end
code[x_, y_, z_, t_, a_, b_] := N[(N[(x / a), $MachinePrecision] / y), $MachinePrecision]
\begin{array}{l}
\\
\frac{\frac{x}{a}}{y}
\end{array}
Initial program 97.4%
Taylor expanded in b around 0
exp-sumN/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
*-commutativeN/A
exp-to-powN/A
lower-pow.f64N/A
exp-prodN/A
lower-pow.f64N/A
rem-exp-logN/A
lower--.f6471.8
Applied rewrites71.8%
Taylor expanded in t around 0
Applied rewrites64.2%
Taylor expanded in y around 0
Applied rewrites34.5%
(FPCore (x y z t a b) :precision binary64 (/ x (* a y)))
double code(double x, double y, double z, double t, double a, double b) {
return x / (a * 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 / (a * y)
end function
public static double code(double x, double y, double z, double t, double a, double b) {
return x / (a * y);
}
def code(x, y, z, t, a, b): return x / (a * y)
function code(x, y, z, t, a, b) return Float64(x / Float64(a * y)) end
function tmp = code(x, y, z, t, a, b) tmp = x / (a * y); end
code[x_, y_, z_, t_, a_, b_] := N[(x / N[(a * y), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{x}{a \cdot y}
\end{array}
Initial program 97.4%
lift-/.f64N/A
lift-*.f64N/A
*-commutativeN/A
associate-/l*N/A
lift-exp.f64N/A
lift--.f64N/A
exp-diffN/A
clear-numN/A
frac-timesN/A
lower-/.f64N/A
lower-*.f64N/A
lower-*.f64N/A
Applied rewrites76.6%
Taylor expanded in t around 0
*-commutativeN/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-exp.f64N/A
lower-/.f64N/A
lower-pow.f6472.4
Applied rewrites72.4%
Taylor expanded in y around 0
Applied rewrites59.9%
Taylor expanded in b around 0
Applied rewrites33.1%
Final simplification33.1%
(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 2024268
(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))