
(FPCore (x y z t a) :precision binary64 (+ x (* (- y z) (/ (- t x) (- a z)))))
double code(double x, double y, double z, double t, double a) {
return x + ((y - z) * ((t - x) / (a - z)));
}
real(8) function code(x, y, z, t, a)
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
code = x + ((y - z) * ((t - x) / (a - z)))
end function
public static double code(double x, double y, double z, double t, double a) {
return x + ((y - z) * ((t - x) / (a - z)));
}
def code(x, y, z, t, a): return x + ((y - z) * ((t - x) / (a - z)))
function code(x, y, z, t, a) return Float64(x + Float64(Float64(y - z) * Float64(Float64(t - x) / Float64(a - z)))) end
function tmp = code(x, y, z, t, a) tmp = x + ((y - z) * ((t - x) / (a - z))); end
code[x_, y_, z_, t_, a_] := N[(x + N[(N[(y - z), $MachinePrecision] * N[(N[(t - x), $MachinePrecision] / N[(a - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x + \left(y - z\right) \cdot \frac{t - x}{a - z}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 16 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x y z t a) :precision binary64 (+ x (* (- y z) (/ (- t x) (- a z)))))
double code(double x, double y, double z, double t, double a) {
return x + ((y - z) * ((t - x) / (a - z)));
}
real(8) function code(x, y, z, t, a)
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
code = x + ((y - z) * ((t - x) / (a - z)))
end function
public static double code(double x, double y, double z, double t, double a) {
return x + ((y - z) * ((t - x) / (a - z)));
}
def code(x, y, z, t, a): return x + ((y - z) * ((t - x) / (a - z)))
function code(x, y, z, t, a) return Float64(x + Float64(Float64(y - z) * Float64(Float64(t - x) / Float64(a - z)))) end
function tmp = code(x, y, z, t, a) tmp = x + ((y - z) * ((t - x) / (a - z))); end
code[x_, y_, z_, t_, a_] := N[(x + N[(N[(y - z), $MachinePrecision] * N[(N[(t - x), $MachinePrecision] / N[(a - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x + \left(y - z\right) \cdot \frac{t - x}{a - z}
\end{array}
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (+ x (* (- y z) (/ (- t x) (- a z))))))
(if (<= t_1 -2e-259)
(fma (/ (- y z) (- a z)) (- t x) x)
(if (<= t_1 0.0)
(fma (- t x) (/ (- a y) z) t)
(fma (+ (/ y (- a z)) (/ z (- z a))) (- t x) x)))))
double code(double x, double y, double z, double t, double a) {
double t_1 = x + ((y - z) * ((t - x) / (a - z)));
double tmp;
if (t_1 <= -2e-259) {
tmp = fma(((y - z) / (a - z)), (t - x), x);
} else if (t_1 <= 0.0) {
tmp = fma((t - x), ((a - y) / z), t);
} else {
tmp = fma(((y / (a - z)) + (z / (z - a))), (t - x), x);
}
return tmp;
}
function code(x, y, z, t, a) t_1 = Float64(x + Float64(Float64(y - z) * Float64(Float64(t - x) / Float64(a - z)))) tmp = 0.0 if (t_1 <= -2e-259) tmp = fma(Float64(Float64(y - z) / Float64(a - z)), Float64(t - x), x); elseif (t_1 <= 0.0) tmp = fma(Float64(t - x), Float64(Float64(a - y) / z), t); else tmp = fma(Float64(Float64(y / Float64(a - z)) + Float64(z / Float64(z - a))), Float64(t - x), x); end return tmp end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(x + N[(N[(y - z), $MachinePrecision] * N[(N[(t - x), $MachinePrecision] / N[(a - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, -2e-259], N[(N[(N[(y - z), $MachinePrecision] / N[(a - z), $MachinePrecision]), $MachinePrecision] * N[(t - x), $MachinePrecision] + x), $MachinePrecision], If[LessEqual[t$95$1, 0.0], N[(N[(t - x), $MachinePrecision] * N[(N[(a - y), $MachinePrecision] / z), $MachinePrecision] + t), $MachinePrecision], N[(N[(N[(y / N[(a - z), $MachinePrecision]), $MachinePrecision] + N[(z / N[(z - a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(t - x), $MachinePrecision] + x), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x + \left(y - z\right) \cdot \frac{t - x}{a - z}\\
\mathbf{if}\;t\_1 \leq -2 \cdot 10^{-259}:\\
\;\;\;\;\mathsf{fma}\left(\frac{y - z}{a - z}, t - x, x\right)\\
\mathbf{elif}\;t\_1 \leq 0:\\
\;\;\;\;\mathsf{fma}\left(t - x, \frac{a - y}{z}, t\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\frac{y}{a - z} + \frac{z}{z - a}, t - x, x\right)\\
\end{array}
\end{array}
if (+.f64 x (*.f64 (-.f64 y z) (/.f64 (-.f64 t x) (-.f64 a z)))) < -2.0000000000000001e-259Initial program 90.5%
+-commutativeN/A
clear-numN/A
associate-*r/N/A
div-invN/A
times-fracN/A
flip3--N/A
clear-numN/A
clear-numN/A
flip3--N/A
accelerator-lowering-fma.f64N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
--lowering--.f64N/A
--lowering--.f6496.2
Applied egg-rr96.2%
if -2.0000000000000001e-259 < (+.f64 x (*.f64 (-.f64 y z) (/.f64 (-.f64 t x) (-.f64 a z)))) < 0.0Initial program 2.9%
+-commutativeN/A
clear-numN/A
associate-*r/N/A
div-invN/A
times-fracN/A
flip3--N/A
clear-numN/A
clear-numN/A
flip3--N/A
accelerator-lowering-fma.f64N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
--lowering--.f64N/A
--lowering--.f642.9
Applied egg-rr2.9%
div-subN/A
--lowering--.f64N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
/-lowering-/.f64N/A
--lowering--.f642.9
Applied egg-rr2.9%
Taylor expanded in z around inf
+-commutativeN/A
*-commutativeN/A
+-commutativeN/A
mul-1-negN/A
remove-double-negN/A
distribute-neg-inN/A
sub-negN/A
distribute-rgt-neg-inN/A
distribute-neg-fracN/A
associate-/l*N/A
distribute-rgt-neg-inN/A
mul-1-negN/A
accelerator-lowering-fma.f64N/A
Simplified99.8%
if 0.0 < (+.f64 x (*.f64 (-.f64 y z) (/.f64 (-.f64 t x) (-.f64 a z)))) Initial program 90.8%
+-commutativeN/A
clear-numN/A
associate-*r/N/A
div-invN/A
times-fracN/A
flip3--N/A
clear-numN/A
clear-numN/A
flip3--N/A
accelerator-lowering-fma.f64N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
--lowering--.f64N/A
--lowering--.f6494.6
Applied egg-rr94.6%
div-subN/A
--lowering--.f64N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
/-lowering-/.f64N/A
--lowering--.f6494.6
Applied egg-rr94.6%
Final simplification95.9%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (fma (/ (- y z) (- a z)) (- t x) x))
(t_2 (+ x (* (- y z) (/ (- t x) (- a z))))))
(if (<= t_2 -2e-259)
t_1
(if (<= t_2 0.0) (fma (- t x) (/ (- a y) z) t) t_1))))
double code(double x, double y, double z, double t, double a) {
double t_1 = fma(((y - z) / (a - z)), (t - x), x);
double t_2 = x + ((y - z) * ((t - x) / (a - z)));
double tmp;
if (t_2 <= -2e-259) {
tmp = t_1;
} else if (t_2 <= 0.0) {
tmp = fma((t - x), ((a - y) / z), t);
} else {
tmp = t_1;
}
return tmp;
}
function code(x, y, z, t, a) t_1 = fma(Float64(Float64(y - z) / Float64(a - z)), Float64(t - x), x) t_2 = Float64(x + Float64(Float64(y - z) * Float64(Float64(t - x) / Float64(a - z)))) tmp = 0.0 if (t_2 <= -2e-259) tmp = t_1; elseif (t_2 <= 0.0) tmp = fma(Float64(t - x), Float64(Float64(a - y) / z), t); else tmp = t_1; end return tmp end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(N[(N[(y - z), $MachinePrecision] / N[(a - z), $MachinePrecision]), $MachinePrecision] * N[(t - x), $MachinePrecision] + x), $MachinePrecision]}, Block[{t$95$2 = N[(x + N[(N[(y - z), $MachinePrecision] * N[(N[(t - x), $MachinePrecision] / N[(a - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$2, -2e-259], t$95$1, If[LessEqual[t$95$2, 0.0], N[(N[(t - x), $MachinePrecision] * N[(N[(a - y), $MachinePrecision] / z), $MachinePrecision] + t), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \mathsf{fma}\left(\frac{y - z}{a - z}, t - x, x\right)\\
t_2 := x + \left(y - z\right) \cdot \frac{t - x}{a - z}\\
\mathbf{if}\;t\_2 \leq -2 \cdot 10^{-259}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t\_2 \leq 0:\\
\;\;\;\;\mathsf{fma}\left(t - x, \frac{a - y}{z}, t\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if (+.f64 x (*.f64 (-.f64 y z) (/.f64 (-.f64 t x) (-.f64 a z)))) < -2.0000000000000001e-259 or 0.0 < (+.f64 x (*.f64 (-.f64 y z) (/.f64 (-.f64 t x) (-.f64 a z)))) Initial program 90.7%
+-commutativeN/A
clear-numN/A
associate-*r/N/A
div-invN/A
times-fracN/A
flip3--N/A
clear-numN/A
clear-numN/A
flip3--N/A
accelerator-lowering-fma.f64N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
--lowering--.f64N/A
--lowering--.f6495.4
Applied egg-rr95.4%
if -2.0000000000000001e-259 < (+.f64 x (*.f64 (-.f64 y z) (/.f64 (-.f64 t x) (-.f64 a z)))) < 0.0Initial program 2.9%
+-commutativeN/A
clear-numN/A
associate-*r/N/A
div-invN/A
times-fracN/A
flip3--N/A
clear-numN/A
clear-numN/A
flip3--N/A
accelerator-lowering-fma.f64N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
--lowering--.f64N/A
--lowering--.f642.9
Applied egg-rr2.9%
div-subN/A
--lowering--.f64N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
/-lowering-/.f64N/A
--lowering--.f642.9
Applied egg-rr2.9%
Taylor expanded in z around inf
+-commutativeN/A
*-commutativeN/A
+-commutativeN/A
mul-1-negN/A
remove-double-negN/A
distribute-neg-inN/A
sub-negN/A
distribute-rgt-neg-inN/A
distribute-neg-fracN/A
associate-/l*N/A
distribute-rgt-neg-inN/A
mul-1-negN/A
accelerator-lowering-fma.f64N/A
Simplified99.8%
(FPCore (x y z t a)
:precision binary64
(if (<= a -1.9e+47)
(fma (- y z) (/ (- t x) a) x)
(if (<= a 2.4e-119)
(fma (- t x) (/ (- a y) z) t)
(+ x (* (- y z) (/ t (- a z)))))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (a <= -1.9e+47) {
tmp = fma((y - z), ((t - x) / a), x);
} else if (a <= 2.4e-119) {
tmp = fma((t - x), ((a - y) / z), t);
} else {
tmp = x + ((y - z) * (t / (a - z)));
}
return tmp;
}
function code(x, y, z, t, a) tmp = 0.0 if (a <= -1.9e+47) tmp = fma(Float64(y - z), Float64(Float64(t - x) / a), x); elseif (a <= 2.4e-119) tmp = fma(Float64(t - x), Float64(Float64(a - y) / z), t); else tmp = Float64(x + Float64(Float64(y - z) * Float64(t / Float64(a - z)))); end return tmp end
code[x_, y_, z_, t_, a_] := If[LessEqual[a, -1.9e+47], N[(N[(y - z), $MachinePrecision] * N[(N[(t - x), $MachinePrecision] / a), $MachinePrecision] + x), $MachinePrecision], If[LessEqual[a, 2.4e-119], N[(N[(t - x), $MachinePrecision] * N[(N[(a - y), $MachinePrecision] / z), $MachinePrecision] + t), $MachinePrecision], N[(x + N[(N[(y - z), $MachinePrecision] * N[(t / N[(a - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq -1.9 \cdot 10^{+47}:\\
\;\;\;\;\mathsf{fma}\left(y - z, \frac{t - x}{a}, x\right)\\
\mathbf{elif}\;a \leq 2.4 \cdot 10^{-119}:\\
\;\;\;\;\mathsf{fma}\left(t - x, \frac{a - y}{z}, t\right)\\
\mathbf{else}:\\
\;\;\;\;x + \left(y - z\right) \cdot \frac{t}{a - z}\\
\end{array}
\end{array}
if a < -1.9000000000000002e47Initial program 93.8%
Taylor expanded in a around inf
+-commutativeN/A
*-commutativeN/A
associate-/l*N/A
accelerator-lowering-fma.f64N/A
--lowering--.f64N/A
/-lowering-/.f64N/A
--lowering--.f6486.7
Simplified86.7%
if -1.9000000000000002e47 < a < 2.40000000000000009e-119Initial program 72.2%
+-commutativeN/A
clear-numN/A
associate-*r/N/A
div-invN/A
times-fracN/A
flip3--N/A
clear-numN/A
clear-numN/A
flip3--N/A
accelerator-lowering-fma.f64N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
--lowering--.f64N/A
--lowering--.f6478.1
Applied egg-rr78.1%
div-subN/A
--lowering--.f64N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
/-lowering-/.f64N/A
--lowering--.f6478.1
Applied egg-rr78.1%
Taylor expanded in z around inf
+-commutativeN/A
*-commutativeN/A
+-commutativeN/A
mul-1-negN/A
remove-double-negN/A
distribute-neg-inN/A
sub-negN/A
distribute-rgt-neg-inN/A
distribute-neg-fracN/A
associate-/l*N/A
distribute-rgt-neg-inN/A
mul-1-negN/A
accelerator-lowering-fma.f64N/A
Simplified85.9%
if 2.40000000000000009e-119 < a Initial program 84.1%
Taylor expanded in t around inf
/-lowering-/.f64N/A
--lowering--.f6475.0
Simplified75.0%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (fma (- y z) (/ (- t x) a) x)))
(if (<= a -1.1e+54)
t_1
(if (<= a 4.8e-21) (fma (- t x) (/ (- a y) z) t) t_1))))
double code(double x, double y, double z, double t, double a) {
double t_1 = fma((y - z), ((t - x) / a), x);
double tmp;
if (a <= -1.1e+54) {
tmp = t_1;
} else if (a <= 4.8e-21) {
tmp = fma((t - x), ((a - y) / z), t);
} else {
tmp = t_1;
}
return tmp;
}
function code(x, y, z, t, a) t_1 = fma(Float64(y - z), Float64(Float64(t - x) / a), x) tmp = 0.0 if (a <= -1.1e+54) tmp = t_1; elseif (a <= 4.8e-21) tmp = fma(Float64(t - x), Float64(Float64(a - y) / z), t); else tmp = t_1; end return tmp end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(N[(y - z), $MachinePrecision] * N[(N[(t - x), $MachinePrecision] / a), $MachinePrecision] + x), $MachinePrecision]}, If[LessEqual[a, -1.1e+54], t$95$1, If[LessEqual[a, 4.8e-21], N[(N[(t - x), $MachinePrecision] * N[(N[(a - y), $MachinePrecision] / z), $MachinePrecision] + t), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \mathsf{fma}\left(y - z, \frac{t - x}{a}, x\right)\\
\mathbf{if}\;a \leq -1.1 \cdot 10^{+54}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;a \leq 4.8 \cdot 10^{-21}:\\
\;\;\;\;\mathsf{fma}\left(t - x, \frac{a - y}{z}, t\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if a < -1.09999999999999995e54 or 4.7999999999999999e-21 < a Initial program 90.0%
Taylor expanded in a around inf
+-commutativeN/A
*-commutativeN/A
associate-/l*N/A
accelerator-lowering-fma.f64N/A
--lowering--.f64N/A
/-lowering-/.f64N/A
--lowering--.f6479.0
Simplified79.0%
if -1.09999999999999995e54 < a < 4.7999999999999999e-21Initial program 72.6%
+-commutativeN/A
clear-numN/A
associate-*r/N/A
div-invN/A
times-fracN/A
flip3--N/A
clear-numN/A
clear-numN/A
flip3--N/A
accelerator-lowering-fma.f64N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
--lowering--.f64N/A
--lowering--.f6479.6
Applied egg-rr79.6%
div-subN/A
--lowering--.f64N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
/-lowering-/.f64N/A
--lowering--.f6479.6
Applied egg-rr79.6%
Taylor expanded in z around inf
+-commutativeN/A
*-commutativeN/A
+-commutativeN/A
mul-1-negN/A
remove-double-negN/A
distribute-neg-inN/A
sub-negN/A
distribute-rgt-neg-inN/A
distribute-neg-fracN/A
associate-/l*N/A
distribute-rgt-neg-inN/A
mul-1-negN/A
accelerator-lowering-fma.f64N/A
Simplified80.1%
(FPCore (x y z t a) :precision binary64 (let* ((t_1 (fma (- t x) (/ (- a y) z) t))) (if (<= z -3.4e-64) t_1 (if (<= z 2.4e+15) (fma (/ y a) (- t x) x) t_1))))
double code(double x, double y, double z, double t, double a) {
double t_1 = fma((t - x), ((a - y) / z), t);
double tmp;
if (z <= -3.4e-64) {
tmp = t_1;
} else if (z <= 2.4e+15) {
tmp = fma((y / a), (t - x), x);
} else {
tmp = t_1;
}
return tmp;
}
function code(x, y, z, t, a) t_1 = fma(Float64(t - x), Float64(Float64(a - y) / z), t) tmp = 0.0 if (z <= -3.4e-64) tmp = t_1; elseif (z <= 2.4e+15) tmp = fma(Float64(y / a), Float64(t - x), x); else tmp = t_1; end return tmp end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(N[(t - x), $MachinePrecision] * N[(N[(a - y), $MachinePrecision] / z), $MachinePrecision] + t), $MachinePrecision]}, If[LessEqual[z, -3.4e-64], t$95$1, If[LessEqual[z, 2.4e+15], N[(N[(y / a), $MachinePrecision] * N[(t - x), $MachinePrecision] + x), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \mathsf{fma}\left(t - x, \frac{a - y}{z}, t\right)\\
\mathbf{if}\;z \leq -3.4 \cdot 10^{-64}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;z \leq 2.4 \cdot 10^{+15}:\\
\;\;\;\;\mathsf{fma}\left(\frac{y}{a}, t - x, x\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if z < -3.40000000000000012e-64 or 2.4e15 < z Initial program 73.1%
+-commutativeN/A
clear-numN/A
associate-*r/N/A
div-invN/A
times-fracN/A
flip3--N/A
clear-numN/A
clear-numN/A
flip3--N/A
accelerator-lowering-fma.f64N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
--lowering--.f64N/A
--lowering--.f6475.5
Applied egg-rr75.5%
div-subN/A
--lowering--.f64N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
/-lowering-/.f64N/A
--lowering--.f6475.5
Applied egg-rr75.5%
Taylor expanded in z around inf
+-commutativeN/A
*-commutativeN/A
+-commutativeN/A
mul-1-negN/A
remove-double-negN/A
distribute-neg-inN/A
sub-negN/A
distribute-rgt-neg-inN/A
distribute-neg-fracN/A
associate-/l*N/A
distribute-rgt-neg-inN/A
mul-1-negN/A
accelerator-lowering-fma.f64N/A
Simplified74.3%
if -3.40000000000000012e-64 < z < 2.4e15Initial program 90.2%
+-commutativeN/A
clear-numN/A
associate-*r/N/A
div-invN/A
times-fracN/A
flip3--N/A
clear-numN/A
clear-numN/A
flip3--N/A
accelerator-lowering-fma.f64N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
--lowering--.f64N/A
--lowering--.f6496.6
Applied egg-rr96.6%
Taylor expanded in z around 0
/-lowering-/.f6475.5
Simplified75.5%
(FPCore (x y z t a) :precision binary64 (if (<= z -4.2e-65) (* (- y z) (/ t (- a z))) (if (<= z 1.35e+48) (fma (/ y a) (- t x) x) (fma t (/ y (- z)) t))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (z <= -4.2e-65) {
tmp = (y - z) * (t / (a - z));
} else if (z <= 1.35e+48) {
tmp = fma((y / a), (t - x), x);
} else {
tmp = fma(t, (y / -z), t);
}
return tmp;
}
function code(x, y, z, t, a) tmp = 0.0 if (z <= -4.2e-65) tmp = Float64(Float64(y - z) * Float64(t / Float64(a - z))); elseif (z <= 1.35e+48) tmp = fma(Float64(y / a), Float64(t - x), x); else tmp = fma(t, Float64(y / Float64(-z)), t); end return tmp end
code[x_, y_, z_, t_, a_] := If[LessEqual[z, -4.2e-65], N[(N[(y - z), $MachinePrecision] * N[(t / N[(a - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, 1.35e+48], N[(N[(y / a), $MachinePrecision] * N[(t - x), $MachinePrecision] + x), $MachinePrecision], N[(t * N[(y / (-z)), $MachinePrecision] + t), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -4.2 \cdot 10^{-65}:\\
\;\;\;\;\left(y - z\right) \cdot \frac{t}{a - z}\\
\mathbf{elif}\;z \leq 1.35 \cdot 10^{+48}:\\
\;\;\;\;\mathsf{fma}\left(\frac{y}{a}, t - x, x\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(t, \frac{y}{-z}, t\right)\\
\end{array}
\end{array}
if z < -4.20000000000000006e-65Initial program 72.9%
Taylor expanded in x around 0
/-lowering-/.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
--lowering--.f6442.9
Simplified42.9%
associate-/l*N/A
clear-numN/A
associate-*r/N/A
div-invN/A
times-fracN/A
flip--N/A
clear-numN/A
clear-numN/A
flip--N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
--lowering--.f6453.6
Applied egg-rr53.6%
if -4.20000000000000006e-65 < z < 1.35000000000000002e48Initial program 90.6%
+-commutativeN/A
clear-numN/A
associate-*r/N/A
div-invN/A
times-fracN/A
flip3--N/A
clear-numN/A
clear-numN/A
flip3--N/A
accelerator-lowering-fma.f64N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
--lowering--.f64N/A
--lowering--.f6496.7
Applied egg-rr96.7%
Taylor expanded in z around 0
/-lowering-/.f6474.9
Simplified74.9%
if 1.35000000000000002e48 < z Initial program 70.7%
+-commutativeN/A
clear-numN/A
associate-*r/N/A
div-invN/A
times-fracN/A
flip3--N/A
clear-numN/A
clear-numN/A
flip3--N/A
accelerator-lowering-fma.f64N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
--lowering--.f64N/A
--lowering--.f6475.7
Applied egg-rr75.7%
Taylor expanded in z around inf
associate--l+N/A
associate-*r/N/A
associate-*r/N/A
mul-1-negN/A
div-subN/A
mul-1-negN/A
distribute-lft-out--N/A
associate-*r/N/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
/-lowering-/.f64N/A
Simplified64.6%
Taylor expanded in t around -inf
+-commutativeN/A
distribute-lft-inN/A
*-rgt-identityN/A
accelerator-lowering-fma.f64N/A
associate-*r/N/A
/-lowering-/.f64N/A
sub-negN/A
distribute-lft-inN/A
neg-mul-1N/A
remove-double-negN/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
--lowering--.f6459.9
Simplified59.9%
Taylor expanded in a around 0
mul-1-negN/A
neg-lowering-neg.f6459.9
Simplified59.9%
Final simplification65.2%
(FPCore (x y z t a) :precision binary64 (let* ((t_1 (fma t (/ y (- z)) t))) (if (<= z -3.25e+66) t_1 (if (<= z 8.5e+39) (fma (/ y a) (- t x) x) t_1))))
double code(double x, double y, double z, double t, double a) {
double t_1 = fma(t, (y / -z), t);
double tmp;
if (z <= -3.25e+66) {
tmp = t_1;
} else if (z <= 8.5e+39) {
tmp = fma((y / a), (t - x), x);
} else {
tmp = t_1;
}
return tmp;
}
function code(x, y, z, t, a) t_1 = fma(t, Float64(y / Float64(-z)), t) tmp = 0.0 if (z <= -3.25e+66) tmp = t_1; elseif (z <= 8.5e+39) tmp = fma(Float64(y / a), Float64(t - x), x); else tmp = t_1; end return tmp end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(t * N[(y / (-z)), $MachinePrecision] + t), $MachinePrecision]}, If[LessEqual[z, -3.25e+66], t$95$1, If[LessEqual[z, 8.5e+39], N[(N[(y / a), $MachinePrecision] * N[(t - x), $MachinePrecision] + x), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \mathsf{fma}\left(t, \frac{y}{-z}, t\right)\\
\mathbf{if}\;z \leq -3.25 \cdot 10^{+66}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;z \leq 8.5 \cdot 10^{+39}:\\
\;\;\;\;\mathsf{fma}\left(\frac{y}{a}, t - x, x\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if z < -3.2500000000000001e66 or 8.49999999999999971e39 < z Initial program 64.8%
+-commutativeN/A
clear-numN/A
associate-*r/N/A
div-invN/A
times-fracN/A
flip3--N/A
clear-numN/A
clear-numN/A
flip3--N/A
accelerator-lowering-fma.f64N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
--lowering--.f64N/A
--lowering--.f6468.1
Applied egg-rr68.1%
Taylor expanded in z around inf
associate--l+N/A
associate-*r/N/A
associate-*r/N/A
mul-1-negN/A
div-subN/A
mul-1-negN/A
distribute-lft-out--N/A
associate-*r/N/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
/-lowering-/.f64N/A
Simplified60.3%
Taylor expanded in t around -inf
+-commutativeN/A
distribute-lft-inN/A
*-rgt-identityN/A
accelerator-lowering-fma.f64N/A
associate-*r/N/A
/-lowering-/.f64N/A
sub-negN/A
distribute-lft-inN/A
neg-mul-1N/A
remove-double-negN/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
--lowering--.f6455.9
Simplified55.9%
Taylor expanded in a around 0
mul-1-negN/A
neg-lowering-neg.f6456.0
Simplified56.0%
if -3.2500000000000001e66 < z < 8.49999999999999971e39Initial program 91.3%
+-commutativeN/A
clear-numN/A
associate-*r/N/A
div-invN/A
times-fracN/A
flip3--N/A
clear-numN/A
clear-numN/A
flip3--N/A
accelerator-lowering-fma.f64N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
--lowering--.f64N/A
--lowering--.f6496.1
Applied egg-rr96.1%
Taylor expanded in z around 0
/-lowering-/.f6468.0
Simplified68.0%
Final simplification63.3%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (fma y (/ (- t x) a) x)))
(if (<= a -1.1e+44)
t_1
(if (<= a 190000000000.0) (fma t (/ (- a y) z) t) t_1))))
double code(double x, double y, double z, double t, double a) {
double t_1 = fma(y, ((t - x) / a), x);
double tmp;
if (a <= -1.1e+44) {
tmp = t_1;
} else if (a <= 190000000000.0) {
tmp = fma(t, ((a - y) / z), t);
} else {
tmp = t_1;
}
return tmp;
}
function code(x, y, z, t, a) t_1 = fma(y, Float64(Float64(t - x) / a), x) tmp = 0.0 if (a <= -1.1e+44) tmp = t_1; elseif (a <= 190000000000.0) tmp = fma(t, Float64(Float64(a - y) / z), t); else tmp = t_1; end return tmp end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(y * N[(N[(t - x), $MachinePrecision] / a), $MachinePrecision] + x), $MachinePrecision]}, If[LessEqual[a, -1.1e+44], t$95$1, If[LessEqual[a, 190000000000.0], N[(t * N[(N[(a - y), $MachinePrecision] / z), $MachinePrecision] + t), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \mathsf{fma}\left(y, \frac{t - x}{a}, x\right)\\
\mathbf{if}\;a \leq -1.1 \cdot 10^{+44}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;a \leq 190000000000:\\
\;\;\;\;\mathsf{fma}\left(t, \frac{a - y}{z}, t\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if a < -1.09999999999999998e44 or 1.9e11 < a Initial program 90.4%
Taylor expanded in z around 0
+-commutativeN/A
associate-/l*N/A
accelerator-lowering-fma.f64N/A
/-lowering-/.f64N/A
--lowering--.f6470.5
Simplified70.5%
if -1.09999999999999998e44 < a < 1.9e11Initial program 73.6%
+-commutativeN/A
clear-numN/A
associate-*r/N/A
div-invN/A
times-fracN/A
flip3--N/A
clear-numN/A
clear-numN/A
flip3--N/A
accelerator-lowering-fma.f64N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
--lowering--.f64N/A
--lowering--.f6480.5
Applied egg-rr80.5%
Taylor expanded in z around inf
associate--l+N/A
associate-*r/N/A
associate-*r/N/A
mul-1-negN/A
div-subN/A
mul-1-negN/A
distribute-lft-out--N/A
associate-*r/N/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
/-lowering-/.f64N/A
Simplified68.1%
Taylor expanded in t around -inf
+-commutativeN/A
distribute-lft-inN/A
*-rgt-identityN/A
accelerator-lowering-fma.f64N/A
associate-*r/N/A
/-lowering-/.f64N/A
sub-negN/A
distribute-lft-inN/A
neg-mul-1N/A
remove-double-negN/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
--lowering--.f6455.0
Simplified55.0%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (fma (/ y a) t x)))
(if (<= a -1.56e+66)
t_1
(if (<= a 16000000000.0) (fma t (/ (- a y) z) t) t_1))))
double code(double x, double y, double z, double t, double a) {
double t_1 = fma((y / a), t, x);
double tmp;
if (a <= -1.56e+66) {
tmp = t_1;
} else if (a <= 16000000000.0) {
tmp = fma(t, ((a - y) / z), t);
} else {
tmp = t_1;
}
return tmp;
}
function code(x, y, z, t, a) t_1 = fma(Float64(y / a), t, x) tmp = 0.0 if (a <= -1.56e+66) tmp = t_1; elseif (a <= 16000000000.0) tmp = fma(t, Float64(Float64(a - y) / z), t); else tmp = t_1; end return tmp end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(N[(y / a), $MachinePrecision] * t + x), $MachinePrecision]}, If[LessEqual[a, -1.56e+66], t$95$1, If[LessEqual[a, 16000000000.0], N[(t * N[(N[(a - y), $MachinePrecision] / z), $MachinePrecision] + t), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \mathsf{fma}\left(\frac{y}{a}, t, x\right)\\
\mathbf{if}\;a \leq -1.56 \cdot 10^{+66}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;a \leq 16000000000:\\
\;\;\;\;\mathsf{fma}\left(t, \frac{a - y}{z}, t\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if a < -1.5599999999999999e66 or 1.6e10 < a Initial program 90.3%
+-commutativeN/A
clear-numN/A
associate-*r/N/A
div-invN/A
times-fracN/A
flip3--N/A
clear-numN/A
clear-numN/A
flip3--N/A
accelerator-lowering-fma.f64N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
--lowering--.f64N/A
--lowering--.f6491.2
Applied egg-rr91.2%
Taylor expanded in z around 0
/-lowering-/.f6470.2
Simplified70.2%
Taylor expanded in t around inf
Simplified64.5%
if -1.5599999999999999e66 < a < 1.6e10Initial program 73.8%
+-commutativeN/A
clear-numN/A
associate-*r/N/A
div-invN/A
times-fracN/A
flip3--N/A
clear-numN/A
clear-numN/A
flip3--N/A
accelerator-lowering-fma.f64N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
--lowering--.f64N/A
--lowering--.f6480.6
Applied egg-rr80.6%
Taylor expanded in z around inf
associate--l+N/A
associate-*r/N/A
associate-*r/N/A
mul-1-negN/A
div-subN/A
mul-1-negN/A
distribute-lft-out--N/A
associate-*r/N/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
/-lowering-/.f64N/A
Simplified67.7%
Taylor expanded in t around -inf
+-commutativeN/A
distribute-lft-inN/A
*-rgt-identityN/A
accelerator-lowering-fma.f64N/A
associate-*r/N/A
/-lowering-/.f64N/A
sub-negN/A
distribute-lft-inN/A
neg-mul-1N/A
remove-double-negN/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
--lowering--.f6454.6
Simplified54.6%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (fma (/ y a) t x)))
(if (<= a -4.5e+62)
t_1
(if (<= a 90000000000.0) (fma t (/ y (- z)) t) t_1))))
double code(double x, double y, double z, double t, double a) {
double t_1 = fma((y / a), t, x);
double tmp;
if (a <= -4.5e+62) {
tmp = t_1;
} else if (a <= 90000000000.0) {
tmp = fma(t, (y / -z), t);
} else {
tmp = t_1;
}
return tmp;
}
function code(x, y, z, t, a) t_1 = fma(Float64(y / a), t, x) tmp = 0.0 if (a <= -4.5e+62) tmp = t_1; elseif (a <= 90000000000.0) tmp = fma(t, Float64(y / Float64(-z)), t); else tmp = t_1; end return tmp end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(N[(y / a), $MachinePrecision] * t + x), $MachinePrecision]}, If[LessEqual[a, -4.5e+62], t$95$1, If[LessEqual[a, 90000000000.0], N[(t * N[(y / (-z)), $MachinePrecision] + t), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \mathsf{fma}\left(\frac{y}{a}, t, x\right)\\
\mathbf{if}\;a \leq -4.5 \cdot 10^{+62}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;a \leq 90000000000:\\
\;\;\;\;\mathsf{fma}\left(t, \frac{y}{-z}, t\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if a < -4.49999999999999999e62 or 9e10 < a Initial program 90.3%
+-commutativeN/A
clear-numN/A
associate-*r/N/A
div-invN/A
times-fracN/A
flip3--N/A
clear-numN/A
clear-numN/A
flip3--N/A
accelerator-lowering-fma.f64N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
--lowering--.f64N/A
--lowering--.f6491.2
Applied egg-rr91.2%
Taylor expanded in z around 0
/-lowering-/.f6470.2
Simplified70.2%
Taylor expanded in t around inf
Simplified64.5%
if -4.49999999999999999e62 < a < 9e10Initial program 73.8%
+-commutativeN/A
clear-numN/A
associate-*r/N/A
div-invN/A
times-fracN/A
flip3--N/A
clear-numN/A
clear-numN/A
flip3--N/A
accelerator-lowering-fma.f64N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
--lowering--.f64N/A
--lowering--.f6480.6
Applied egg-rr80.6%
Taylor expanded in z around inf
associate--l+N/A
associate-*r/N/A
associate-*r/N/A
mul-1-negN/A
div-subN/A
mul-1-negN/A
distribute-lft-out--N/A
associate-*r/N/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
/-lowering-/.f64N/A
Simplified67.7%
Taylor expanded in t around -inf
+-commutativeN/A
distribute-lft-inN/A
*-rgt-identityN/A
accelerator-lowering-fma.f64N/A
associate-*r/N/A
/-lowering-/.f64N/A
sub-negN/A
distribute-lft-inN/A
neg-mul-1N/A
remove-double-negN/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
--lowering--.f6454.6
Simplified54.6%
Taylor expanded in a around 0
mul-1-negN/A
neg-lowering-neg.f6454.3
Simplified54.3%
Final simplification58.8%
(FPCore (x y z t a) :precision binary64 (if (<= z -5.6e+74) t (if (<= z 1e+113) (fma (/ y a) t x) (fma t (/ a z) t))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (z <= -5.6e+74) {
tmp = t;
} else if (z <= 1e+113) {
tmp = fma((y / a), t, x);
} else {
tmp = fma(t, (a / z), t);
}
return tmp;
}
function code(x, y, z, t, a) tmp = 0.0 if (z <= -5.6e+74) tmp = t; elseif (z <= 1e+113) tmp = fma(Float64(y / a), t, x); else tmp = fma(t, Float64(a / z), t); end return tmp end
code[x_, y_, z_, t_, a_] := If[LessEqual[z, -5.6e+74], t, If[LessEqual[z, 1e+113], N[(N[(y / a), $MachinePrecision] * t + x), $MachinePrecision], N[(t * N[(a / z), $MachinePrecision] + t), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -5.6 \cdot 10^{+74}:\\
\;\;\;\;t\\
\mathbf{elif}\;z \leq 10^{+113}:\\
\;\;\;\;\mathsf{fma}\left(\frac{y}{a}, t, x\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(t, \frac{a}{z}, t\right)\\
\end{array}
\end{array}
if z < -5.60000000000000003e74Initial program 58.5%
Taylor expanded in z around inf
Simplified45.4%
if -5.60000000000000003e74 < z < 1e113Initial program 90.6%
+-commutativeN/A
clear-numN/A
associate-*r/N/A
div-invN/A
times-fracN/A
flip3--N/A
clear-numN/A
clear-numN/A
flip3--N/A
accelerator-lowering-fma.f64N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
--lowering--.f64N/A
--lowering--.f6495.6
Applied egg-rr95.6%
Taylor expanded in z around 0
/-lowering-/.f6465.2
Simplified65.2%
Taylor expanded in t around inf
Simplified57.7%
if 1e113 < z Initial program 67.6%
+-commutativeN/A
clear-numN/A
associate-*r/N/A
div-invN/A
times-fracN/A
flip3--N/A
clear-numN/A
clear-numN/A
flip3--N/A
accelerator-lowering-fma.f64N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
--lowering--.f64N/A
--lowering--.f6471.8
Applied egg-rr71.8%
Taylor expanded in z around inf
associate--l+N/A
associate-*r/N/A
associate-*r/N/A
mul-1-negN/A
div-subN/A
mul-1-negN/A
distribute-lft-out--N/A
associate-*r/N/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
/-lowering-/.f64N/A
Simplified68.2%
Taylor expanded in t around -inf
+-commutativeN/A
distribute-lft-inN/A
*-rgt-identityN/A
accelerator-lowering-fma.f64N/A
associate-*r/N/A
/-lowering-/.f64N/A
sub-negN/A
distribute-lft-inN/A
neg-mul-1N/A
remove-double-negN/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
--lowering--.f6467.5
Simplified67.5%
Taylor expanded in a around inf
/-lowering-/.f6462.5
Simplified62.5%
(FPCore (x y z t a) :precision binary64 (if (<= a -9.2e+73) x (if (<= a 1800000000000.0) (fma t (/ a z) t) x)))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (a <= -9.2e+73) {
tmp = x;
} else if (a <= 1800000000000.0) {
tmp = fma(t, (a / z), t);
} else {
tmp = x;
}
return tmp;
}
function code(x, y, z, t, a) tmp = 0.0 if (a <= -9.2e+73) tmp = x; elseif (a <= 1800000000000.0) tmp = fma(t, Float64(a / z), t); else tmp = x; end return tmp end
code[x_, y_, z_, t_, a_] := If[LessEqual[a, -9.2e+73], x, If[LessEqual[a, 1800000000000.0], N[(t * N[(a / z), $MachinePrecision] + t), $MachinePrecision], x]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq -9.2 \cdot 10^{+73}:\\
\;\;\;\;x\\
\mathbf{elif}\;a \leq 1800000000000:\\
\;\;\;\;\mathsf{fma}\left(t, \frac{a}{z}, t\right)\\
\mathbf{else}:\\
\;\;\;\;x\\
\end{array}
\end{array}
if a < -9.199999999999999e73 or 1.8e12 < a Initial program 90.2%
Taylor expanded in a around inf
Simplified45.2%
if -9.199999999999999e73 < a < 1.8e12Initial program 73.9%
+-commutativeN/A
clear-numN/A
associate-*r/N/A
div-invN/A
times-fracN/A
flip3--N/A
clear-numN/A
clear-numN/A
flip3--N/A
accelerator-lowering-fma.f64N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
--lowering--.f64N/A
--lowering--.f6480.8
Applied egg-rr80.8%
Taylor expanded in z around inf
associate--l+N/A
associate-*r/N/A
associate-*r/N/A
mul-1-negN/A
div-subN/A
mul-1-negN/A
distribute-lft-out--N/A
associate-*r/N/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
/-lowering-/.f64N/A
Simplified67.2%
Taylor expanded in t around -inf
+-commutativeN/A
distribute-lft-inN/A
*-rgt-identityN/A
accelerator-lowering-fma.f64N/A
associate-*r/N/A
/-lowering-/.f64N/A
sub-negN/A
distribute-lft-inN/A
neg-mul-1N/A
remove-double-negN/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
--lowering--.f6454.2
Simplified54.2%
Taylor expanded in a around inf
/-lowering-/.f6438.4
Simplified38.4%
(FPCore (x y z t a) :precision binary64 (if (<= a -1.6e+70) x (if (<= a 4000000000000.0) (fma a (/ t z) t) x)))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (a <= -1.6e+70) {
tmp = x;
} else if (a <= 4000000000000.0) {
tmp = fma(a, (t / z), t);
} else {
tmp = x;
}
return tmp;
}
function code(x, y, z, t, a) tmp = 0.0 if (a <= -1.6e+70) tmp = x; elseif (a <= 4000000000000.0) tmp = fma(a, Float64(t / z), t); else tmp = x; end return tmp end
code[x_, y_, z_, t_, a_] := If[LessEqual[a, -1.6e+70], x, If[LessEqual[a, 4000000000000.0], N[(a * N[(t / z), $MachinePrecision] + t), $MachinePrecision], x]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq -1.6 \cdot 10^{+70}:\\
\;\;\;\;x\\
\mathbf{elif}\;a \leq 4000000000000:\\
\;\;\;\;\mathsf{fma}\left(a, \frac{t}{z}, t\right)\\
\mathbf{else}:\\
\;\;\;\;x\\
\end{array}
\end{array}
if a < -1.6000000000000001e70 or 4e12 < a Initial program 90.2%
Taylor expanded in a around inf
Simplified45.2%
if -1.6000000000000001e70 < a < 4e12Initial program 73.9%
+-commutativeN/A
clear-numN/A
associate-*r/N/A
div-invN/A
times-fracN/A
flip3--N/A
clear-numN/A
clear-numN/A
flip3--N/A
accelerator-lowering-fma.f64N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
--lowering--.f64N/A
--lowering--.f6480.8
Applied egg-rr80.8%
Taylor expanded in z around inf
associate--l+N/A
associate-*r/N/A
associate-*r/N/A
mul-1-negN/A
div-subN/A
mul-1-negN/A
distribute-lft-out--N/A
associate-*r/N/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
/-lowering-/.f64N/A
Simplified67.2%
Taylor expanded in t around -inf
+-commutativeN/A
distribute-lft-inN/A
*-rgt-identityN/A
accelerator-lowering-fma.f64N/A
associate-*r/N/A
/-lowering-/.f64N/A
sub-negN/A
distribute-lft-inN/A
neg-mul-1N/A
remove-double-negN/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
--lowering--.f6454.2
Simplified54.2%
Taylor expanded in y around 0
+-commutativeN/A
associate-/l*N/A
accelerator-lowering-fma.f64N/A
/-lowering-/.f6438.3
Simplified38.3%
(FPCore (x y z t a) :precision binary64 (if (<= a -6e+69) x (if (<= a 480000000000.0) t x)))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (a <= -6e+69) {
tmp = x;
} else if (a <= 480000000000.0) {
tmp = t;
} else {
tmp = x;
}
return tmp;
}
real(8) function code(x, y, z, t, a)
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) :: tmp
if (a <= (-6d+69)) then
tmp = x
else if (a <= 480000000000.0d0) then
tmp = t
else
tmp = x
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if (a <= -6e+69) {
tmp = x;
} else if (a <= 480000000000.0) {
tmp = t;
} else {
tmp = x;
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if a <= -6e+69: tmp = x elif a <= 480000000000.0: tmp = t else: tmp = x return tmp
function code(x, y, z, t, a) tmp = 0.0 if (a <= -6e+69) tmp = x; elseif (a <= 480000000000.0) tmp = t; else tmp = x; end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (a <= -6e+69) tmp = x; elseif (a <= 480000000000.0) tmp = t; else tmp = x; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[a, -6e+69], x, If[LessEqual[a, 480000000000.0], t, x]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq -6 \cdot 10^{+69}:\\
\;\;\;\;x\\
\mathbf{elif}\;a \leq 480000000000:\\
\;\;\;\;t\\
\mathbf{else}:\\
\;\;\;\;x\\
\end{array}
\end{array}
if a < -5.99999999999999967e69 or 4.8e11 < a Initial program 90.2%
Taylor expanded in a around inf
Simplified45.2%
if -5.99999999999999967e69 < a < 4.8e11Initial program 73.9%
Taylor expanded in z around inf
Simplified36.9%
(FPCore (x y z t a) :precision binary64 t)
double code(double x, double y, double z, double t, double a) {
return t;
}
real(8) function code(x, y, z, t, a)
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
code = t
end function
public static double code(double x, double y, double z, double t, double a) {
return t;
}
def code(x, y, z, t, a): return t
function code(x, y, z, t, a) return t end
function tmp = code(x, y, z, t, a) tmp = t; end
code[x_, y_, z_, t_, a_] := t
\begin{array}{l}
\\
t
\end{array}
Initial program 81.1%
Taylor expanded in z around inf
Simplified24.4%
(FPCore (x y z t a) :precision binary64 0.0)
double code(double x, double y, double z, double t, double a) {
return 0.0;
}
real(8) function code(x, y, z, t, a)
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
code = 0.0d0
end function
public static double code(double x, double y, double z, double t, double a) {
return 0.0;
}
def code(x, y, z, t, a): return 0.0
function code(x, y, z, t, a) return 0.0 end
function tmp = code(x, y, z, t, a) tmp = 0.0; end
code[x_, y_, z_, t_, a_] := 0.0
\begin{array}{l}
\\
0
\end{array}
Initial program 81.1%
Taylor expanded in x around inf
+-commutativeN/A
distribute-rgt-inN/A
mul-1-negN/A
distribute-lft-neg-outN/A
*-commutativeN/A
associate-/l*N/A
*-commutativeN/A
associate-/l*N/A
distribute-rgt-neg-inN/A
*-lft-identityN/A
accelerator-lowering-fma.f64N/A
--lowering--.f64N/A
neg-lowering-neg.f64N/A
/-lowering-/.f64N/A
--lowering--.f6441.9
Simplified41.9%
Taylor expanded in z around inf
distribute-rgt1-inN/A
metadata-evalN/A
mul0-lft2.7
Simplified2.7%
herbie shell --seed 2024204
(FPCore (x y z t a)
:name "Numeric.Signal:interpolate from hsignal-0.2.7.1"
:precision binary64
(+ x (* (- y z) (/ (- t x) (- a z)))))