
(FPCore (a b) :precision binary64 (- (+ (pow (+ (* a a) (* b b)) 2.0) (* 4.0 (+ (* (* a a) (- 1.0 a)) (* (* b b) (+ 3.0 a))))) 1.0))
double code(double a, double b) {
return (pow(((a * a) + (b * b)), 2.0) + (4.0 * (((a * a) * (1.0 - a)) + ((b * b) * (3.0 + a))))) - 1.0;
}
real(8) function code(a, b)
real(8), intent (in) :: a
real(8), intent (in) :: b
code = ((((a * a) + (b * b)) ** 2.0d0) + (4.0d0 * (((a * a) * (1.0d0 - a)) + ((b * b) * (3.0d0 + a))))) - 1.0d0
end function
public static double code(double a, double b) {
return (Math.pow(((a * a) + (b * b)), 2.0) + (4.0 * (((a * a) * (1.0 - a)) + ((b * b) * (3.0 + a))))) - 1.0;
}
def code(a, b): return (math.pow(((a * a) + (b * b)), 2.0) + (4.0 * (((a * a) * (1.0 - a)) + ((b * b) * (3.0 + a))))) - 1.0
function code(a, b) return Float64(Float64((Float64(Float64(a * a) + Float64(b * b)) ^ 2.0) + Float64(4.0 * Float64(Float64(Float64(a * a) * Float64(1.0 - a)) + Float64(Float64(b * b) * Float64(3.0 + a))))) - 1.0) end
function tmp = code(a, b) tmp = ((((a * a) + (b * b)) ^ 2.0) + (4.0 * (((a * a) * (1.0 - a)) + ((b * b) * (3.0 + a))))) - 1.0; end
code[a_, b_] := N[(N[(N[Power[N[(N[(a * a), $MachinePrecision] + N[(b * b), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision] + N[(4.0 * N[(N[(N[(a * a), $MachinePrecision] * N[(1.0 - a), $MachinePrecision]), $MachinePrecision] + N[(N[(b * b), $MachinePrecision] * N[(3.0 + a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - 1.0), $MachinePrecision]
\begin{array}{l}
\\
\left({\left(a \cdot a + b \cdot b\right)}^{2} + 4 \cdot \left(\left(a \cdot a\right) \cdot \left(1 - a\right) + \left(b \cdot b\right) \cdot \left(3 + a\right)\right)\right) - 1
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 14 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (a b) :precision binary64 (- (+ (pow (+ (* a a) (* b b)) 2.0) (* 4.0 (+ (* (* a a) (- 1.0 a)) (* (* b b) (+ 3.0 a))))) 1.0))
double code(double a, double b) {
return (pow(((a * a) + (b * b)), 2.0) + (4.0 * (((a * a) * (1.0 - a)) + ((b * b) * (3.0 + a))))) - 1.0;
}
real(8) function code(a, b)
real(8), intent (in) :: a
real(8), intent (in) :: b
code = ((((a * a) + (b * b)) ** 2.0d0) + (4.0d0 * (((a * a) * (1.0d0 - a)) + ((b * b) * (3.0d0 + a))))) - 1.0d0
end function
public static double code(double a, double b) {
return (Math.pow(((a * a) + (b * b)), 2.0) + (4.0 * (((a * a) * (1.0 - a)) + ((b * b) * (3.0 + a))))) - 1.0;
}
def code(a, b): return (math.pow(((a * a) + (b * b)), 2.0) + (4.0 * (((a * a) * (1.0 - a)) + ((b * b) * (3.0 + a))))) - 1.0
function code(a, b) return Float64(Float64((Float64(Float64(a * a) + Float64(b * b)) ^ 2.0) + Float64(4.0 * Float64(Float64(Float64(a * a) * Float64(1.0 - a)) + Float64(Float64(b * b) * Float64(3.0 + a))))) - 1.0) end
function tmp = code(a, b) tmp = ((((a * a) + (b * b)) ^ 2.0) + (4.0 * (((a * a) * (1.0 - a)) + ((b * b) * (3.0 + a))))) - 1.0; end
code[a_, b_] := N[(N[(N[Power[N[(N[(a * a), $MachinePrecision] + N[(b * b), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision] + N[(4.0 * N[(N[(N[(a * a), $MachinePrecision] * N[(1.0 - a), $MachinePrecision]), $MachinePrecision] + N[(N[(b * b), $MachinePrecision] * N[(3.0 + a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - 1.0), $MachinePrecision]
\begin{array}{l}
\\
\left({\left(a \cdot a + b \cdot b\right)}^{2} + 4 \cdot \left(\left(a \cdot a\right) \cdot \left(1 - a\right) + \left(b \cdot b\right) \cdot \left(3 + a\right)\right)\right) - 1
\end{array}
(FPCore (a b)
:precision binary64
(let* ((t_0
(+
(pow (+ (* a a) (* b b)) 2.0)
(* 4.0 (+ (* (* a a) (- 1.0 a)) (* (* b b) (+ 3.0 a)))))))
(if (<= t_0 INFINITY)
(- t_0 1.0)
(fma (* (* (* a a) 2.0) b) b (fma (* a a) (fma (- a 4.0) a 4.0) -1.0)))))
double code(double a, double b) {
double t_0 = pow(((a * a) + (b * b)), 2.0) + (4.0 * (((a * a) * (1.0 - a)) + ((b * b) * (3.0 + a))));
double tmp;
if (t_0 <= ((double) INFINITY)) {
tmp = t_0 - 1.0;
} else {
tmp = fma((((a * a) * 2.0) * b), b, fma((a * a), fma((a - 4.0), a, 4.0), -1.0));
}
return tmp;
}
function code(a, b) t_0 = Float64((Float64(Float64(a * a) + Float64(b * b)) ^ 2.0) + Float64(4.0 * Float64(Float64(Float64(a * a) * Float64(1.0 - a)) + Float64(Float64(b * b) * Float64(3.0 + a))))) tmp = 0.0 if (t_0 <= Inf) tmp = Float64(t_0 - 1.0); else tmp = fma(Float64(Float64(Float64(a * a) * 2.0) * b), b, fma(Float64(a * a), fma(Float64(a - 4.0), a, 4.0), -1.0)); end return tmp end
code[a_, b_] := Block[{t$95$0 = N[(N[Power[N[(N[(a * a), $MachinePrecision] + N[(b * b), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision] + N[(4.0 * N[(N[(N[(a * a), $MachinePrecision] * N[(1.0 - a), $MachinePrecision]), $MachinePrecision] + N[(N[(b * b), $MachinePrecision] * N[(3.0 + a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, Infinity], N[(t$95$0 - 1.0), $MachinePrecision], N[(N[(N[(N[(a * a), $MachinePrecision] * 2.0), $MachinePrecision] * b), $MachinePrecision] * b + N[(N[(a * a), $MachinePrecision] * N[(N[(a - 4.0), $MachinePrecision] * a + 4.0), $MachinePrecision] + -1.0), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := {\left(a \cdot a + b \cdot b\right)}^{2} + 4 \cdot \left(\left(a \cdot a\right) \cdot \left(1 - a\right) + \left(b \cdot b\right) \cdot \left(3 + a\right)\right)\\
\mathbf{if}\;t\_0 \leq \infty:\\
\;\;\;\;t\_0 - 1\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\left(\left(a \cdot a\right) \cdot 2\right) \cdot b, b, \mathsf{fma}\left(a \cdot a, \mathsf{fma}\left(a - 4, a, 4\right), -1\right)\right)\\
\end{array}
\end{array}
if (+.f64 (pow.f64 (+.f64 (*.f64 a a) (*.f64 b b)) #s(literal 2 binary64)) (*.f64 #s(literal 4 binary64) (+.f64 (*.f64 (*.f64 a a) (-.f64 #s(literal 1 binary64) a)) (*.f64 (*.f64 b b) (+.f64 #s(literal 3 binary64) a))))) < +inf.0Initial program 99.8%
if +inf.0 < (+.f64 (pow.f64 (+.f64 (*.f64 a a) (*.f64 b b)) #s(literal 2 binary64)) (*.f64 #s(literal 4 binary64) (+.f64 (*.f64 (*.f64 a a) (-.f64 #s(literal 1 binary64) a)) (*.f64 (*.f64 b b) (+.f64 #s(literal 3 binary64) a))))) Initial program 0.0%
Taylor expanded in a around inf
unpow3N/A
unpow2N/A
associate-*l*N/A
sub-negN/A
remove-double-negN/A
distribute-neg-inN/A
+-commutativeN/A
sub-negN/A
distribute-rgt-neg-inN/A
*-commutativeN/A
sub-negN/A
+-commutativeN/A
distribute-lft1-inN/A
distribute-lft-neg-outN/A
lft-mult-inverseN/A
metadata-evalN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
Applied rewrites40.8%
Taylor expanded in b around 0
+-commutativeN/A
associate-+l+N/A
+-commutativeN/A
associate--l+N/A
*-commutativeN/A
distribute-lft-inN/A
metadata-evalN/A
unpow2N/A
associate-*r*N/A
lower-fma.f64N/A
Applied rewrites100.0%
Taylor expanded in a around 0
Applied rewrites100.0%
Taylor expanded in a around inf
Applied rewrites100.0%
(FPCore (a b) :precision binary64 (if (<= a 4.8e+49) (- (+ (pow (+ (* a a) (* b b)) 2.0) (* 4.0 (* (fma (- a) a a) a))) 1.0) (pow a 4.0)))
double code(double a, double b) {
double tmp;
if (a <= 4.8e+49) {
tmp = (pow(((a * a) + (b * b)), 2.0) + (4.0 * (fma(-a, a, a) * a))) - 1.0;
} else {
tmp = pow(a, 4.0);
}
return tmp;
}
function code(a, b) tmp = 0.0 if (a <= 4.8e+49) tmp = Float64(Float64((Float64(Float64(a * a) + Float64(b * b)) ^ 2.0) + Float64(4.0 * Float64(fma(Float64(-a), a, a) * a))) - 1.0); else tmp = a ^ 4.0; end return tmp end
code[a_, b_] := If[LessEqual[a, 4.8e+49], N[(N[(N[Power[N[(N[(a * a), $MachinePrecision] + N[(b * b), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision] + N[(4.0 * N[(N[((-a) * a + a), $MachinePrecision] * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - 1.0), $MachinePrecision], N[Power[a, 4.0], $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq 4.8 \cdot 10^{+49}:\\
\;\;\;\;\left({\left(a \cdot a + b \cdot b\right)}^{2} + 4 \cdot \left(\mathsf{fma}\left(-a, a, a\right) \cdot a\right)\right) - 1\\
\mathbf{else}:\\
\;\;\;\;{a}^{4}\\
\end{array}
\end{array}
if a < 4.8e49Initial program 85.6%
Taylor expanded in a around inf
unpow3N/A
unpow2N/A
associate-*l*N/A
sub-negN/A
remove-double-negN/A
distribute-neg-inN/A
+-commutativeN/A
sub-negN/A
distribute-rgt-neg-inN/A
*-commutativeN/A
sub-negN/A
+-commutativeN/A
distribute-lft1-inN/A
distribute-lft-neg-outN/A
lft-mult-inverseN/A
metadata-evalN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
Applied rewrites99.6%
if 4.8e49 < a Initial program 20.6%
Taylor expanded in a around inf
unpow3N/A
unpow2N/A
associate-*l*N/A
sub-negN/A
remove-double-negN/A
distribute-neg-inN/A
+-commutativeN/A
sub-negN/A
distribute-rgt-neg-inN/A
*-commutativeN/A
sub-negN/A
+-commutativeN/A
distribute-lft1-inN/A
distribute-lft-neg-outN/A
lft-mult-inverseN/A
metadata-evalN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
Applied rewrites20.6%
Taylor expanded in b around 0
+-commutativeN/A
associate-+l+N/A
+-commutativeN/A
associate--l+N/A
*-commutativeN/A
distribute-lft-inN/A
metadata-evalN/A
unpow2N/A
associate-*r*N/A
lower-fma.f64N/A
Applied rewrites99.9%
Taylor expanded in a around inf
lower-pow.f64100.0
Applied rewrites100.0%
(FPCore (a b)
:precision binary64
(let* ((t_0 (fma (* a a) (fma (- a 4.0) a 4.0) -1.0)))
(if (<= a -1600000000.0)
(fma (* (* (* a a) 2.0) b) b t_0)
(if (<= a 6.2e-30)
(- (fma (* b b) 12.0 (pow b 4.0)) 1.0)
(fma (* (fma a (fma 2.0 a 4.0) 12.0) b) b t_0)))))
double code(double a, double b) {
double t_0 = fma((a * a), fma((a - 4.0), a, 4.0), -1.0);
double tmp;
if (a <= -1600000000.0) {
tmp = fma((((a * a) * 2.0) * b), b, t_0);
} else if (a <= 6.2e-30) {
tmp = fma((b * b), 12.0, pow(b, 4.0)) - 1.0;
} else {
tmp = fma((fma(a, fma(2.0, a, 4.0), 12.0) * b), b, t_0);
}
return tmp;
}
function code(a, b) t_0 = fma(Float64(a * a), fma(Float64(a - 4.0), a, 4.0), -1.0) tmp = 0.0 if (a <= -1600000000.0) tmp = fma(Float64(Float64(Float64(a * a) * 2.0) * b), b, t_0); elseif (a <= 6.2e-30) tmp = Float64(fma(Float64(b * b), 12.0, (b ^ 4.0)) - 1.0); else tmp = fma(Float64(fma(a, fma(2.0, a, 4.0), 12.0) * b), b, t_0); end return tmp end
code[a_, b_] := Block[{t$95$0 = N[(N[(a * a), $MachinePrecision] * N[(N[(a - 4.0), $MachinePrecision] * a + 4.0), $MachinePrecision] + -1.0), $MachinePrecision]}, If[LessEqual[a, -1600000000.0], N[(N[(N[(N[(a * a), $MachinePrecision] * 2.0), $MachinePrecision] * b), $MachinePrecision] * b + t$95$0), $MachinePrecision], If[LessEqual[a, 6.2e-30], N[(N[(N[(b * b), $MachinePrecision] * 12.0 + N[Power[b, 4.0], $MachinePrecision]), $MachinePrecision] - 1.0), $MachinePrecision], N[(N[(N[(a * N[(2.0 * a + 4.0), $MachinePrecision] + 12.0), $MachinePrecision] * b), $MachinePrecision] * b + t$95$0), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \mathsf{fma}\left(a \cdot a, \mathsf{fma}\left(a - 4, a, 4\right), -1\right)\\
\mathbf{if}\;a \leq -1600000000:\\
\;\;\;\;\mathsf{fma}\left(\left(\left(a \cdot a\right) \cdot 2\right) \cdot b, b, t\_0\right)\\
\mathbf{elif}\;a \leq 6.2 \cdot 10^{-30}:\\
\;\;\;\;\mathsf{fma}\left(b \cdot b, 12, {b}^{4}\right) - 1\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(a, \mathsf{fma}\left(2, a, 4\right), 12\right) \cdot b, b, t\_0\right)\\
\end{array}
\end{array}
if a < -1.6e9Initial program 53.9%
Taylor expanded in a around inf
unpow3N/A
unpow2N/A
associate-*l*N/A
sub-negN/A
remove-double-negN/A
distribute-neg-inN/A
+-commutativeN/A
sub-negN/A
distribute-rgt-neg-inN/A
*-commutativeN/A
sub-negN/A
+-commutativeN/A
distribute-lft1-inN/A
distribute-lft-neg-outN/A
lft-mult-inverseN/A
metadata-evalN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
Applied rewrites100.0%
Taylor expanded in b around 0
+-commutativeN/A
associate-+l+N/A
+-commutativeN/A
associate--l+N/A
*-commutativeN/A
distribute-lft-inN/A
metadata-evalN/A
unpow2N/A
associate-*r*N/A
lower-fma.f64N/A
Applied rewrites98.5%
Taylor expanded in a around 0
Applied rewrites98.5%
Taylor expanded in a around inf
Applied rewrites98.5%
if -1.6e9 < a < 6.19999999999999982e-30Initial program 99.9%
Taylor expanded in a around 0
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
lower-pow.f6499.8
Applied rewrites99.8%
if 6.19999999999999982e-30 < a Initial program 33.2%
Taylor expanded in a around inf
unpow3N/A
unpow2N/A
associate-*l*N/A
sub-negN/A
remove-double-negN/A
distribute-neg-inN/A
+-commutativeN/A
sub-negN/A
distribute-rgt-neg-inN/A
*-commutativeN/A
sub-negN/A
+-commutativeN/A
distribute-lft1-inN/A
distribute-lft-neg-outN/A
lft-mult-inverseN/A
metadata-evalN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
Applied rewrites33.1%
Taylor expanded in b around 0
+-commutativeN/A
associate-+l+N/A
+-commutativeN/A
associate--l+N/A
*-commutativeN/A
distribute-lft-inN/A
metadata-evalN/A
unpow2N/A
associate-*r*N/A
lower-fma.f64N/A
Applied rewrites98.4%
Taylor expanded in a around 0
Applied rewrites98.4%
(FPCore (a b)
:precision binary64
(let* ((t_0 (fma (* a a) (fma (- a 4.0) a 4.0) -1.0)))
(if (<= a -1600000000.0)
(fma (* (* (* a a) 2.0) b) b t_0)
(if (<= a 6.2e-30)
(fma (* (fma b b 12.0) b) b -1.0)
(fma (* (fma a (fma 2.0 a 4.0) 12.0) b) b t_0)))))
double code(double a, double b) {
double t_0 = fma((a * a), fma((a - 4.0), a, 4.0), -1.0);
double tmp;
if (a <= -1600000000.0) {
tmp = fma((((a * a) * 2.0) * b), b, t_0);
} else if (a <= 6.2e-30) {
tmp = fma((fma(b, b, 12.0) * b), b, -1.0);
} else {
tmp = fma((fma(a, fma(2.0, a, 4.0), 12.0) * b), b, t_0);
}
return tmp;
}
function code(a, b) t_0 = fma(Float64(a * a), fma(Float64(a - 4.0), a, 4.0), -1.0) tmp = 0.0 if (a <= -1600000000.0) tmp = fma(Float64(Float64(Float64(a * a) * 2.0) * b), b, t_0); elseif (a <= 6.2e-30) tmp = fma(Float64(fma(b, b, 12.0) * b), b, -1.0); else tmp = fma(Float64(fma(a, fma(2.0, a, 4.0), 12.0) * b), b, t_0); end return tmp end
code[a_, b_] := Block[{t$95$0 = N[(N[(a * a), $MachinePrecision] * N[(N[(a - 4.0), $MachinePrecision] * a + 4.0), $MachinePrecision] + -1.0), $MachinePrecision]}, If[LessEqual[a, -1600000000.0], N[(N[(N[(N[(a * a), $MachinePrecision] * 2.0), $MachinePrecision] * b), $MachinePrecision] * b + t$95$0), $MachinePrecision], If[LessEqual[a, 6.2e-30], N[(N[(N[(b * b + 12.0), $MachinePrecision] * b), $MachinePrecision] * b + -1.0), $MachinePrecision], N[(N[(N[(a * N[(2.0 * a + 4.0), $MachinePrecision] + 12.0), $MachinePrecision] * b), $MachinePrecision] * b + t$95$0), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \mathsf{fma}\left(a \cdot a, \mathsf{fma}\left(a - 4, a, 4\right), -1\right)\\
\mathbf{if}\;a \leq -1600000000:\\
\;\;\;\;\mathsf{fma}\left(\left(\left(a \cdot a\right) \cdot 2\right) \cdot b, b, t\_0\right)\\
\mathbf{elif}\;a \leq 6.2 \cdot 10^{-30}:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(b, b, 12\right) \cdot b, b, -1\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(a, \mathsf{fma}\left(2, a, 4\right), 12\right) \cdot b, b, t\_0\right)\\
\end{array}
\end{array}
if a < -1.6e9Initial program 53.9%
Taylor expanded in a around inf
unpow3N/A
unpow2N/A
associate-*l*N/A
sub-negN/A
remove-double-negN/A
distribute-neg-inN/A
+-commutativeN/A
sub-negN/A
distribute-rgt-neg-inN/A
*-commutativeN/A
sub-negN/A
+-commutativeN/A
distribute-lft1-inN/A
distribute-lft-neg-outN/A
lft-mult-inverseN/A
metadata-evalN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
Applied rewrites100.0%
Taylor expanded in b around 0
+-commutativeN/A
associate-+l+N/A
+-commutativeN/A
associate--l+N/A
*-commutativeN/A
distribute-lft-inN/A
metadata-evalN/A
unpow2N/A
associate-*r*N/A
lower-fma.f64N/A
Applied rewrites98.5%
Taylor expanded in a around 0
Applied rewrites98.5%
Taylor expanded in a around inf
Applied rewrites98.5%
if -1.6e9 < a < 6.19999999999999982e-30Initial program 99.9%
Taylor expanded in a around inf
unpow3N/A
unpow2N/A
associate-*l*N/A
sub-negN/A
remove-double-negN/A
distribute-neg-inN/A
+-commutativeN/A
sub-negN/A
distribute-rgt-neg-inN/A
*-commutativeN/A
sub-negN/A
+-commutativeN/A
distribute-lft1-inN/A
distribute-lft-neg-outN/A
lft-mult-inverseN/A
metadata-evalN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
Applied rewrites99.4%
Taylor expanded in a around 0
sub-negN/A
+-commutativeN/A
metadata-evalN/A
metadata-evalN/A
pow-sqrN/A
distribute-rgt-inN/A
+-commutativeN/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
unpow2N/A
lower-fma.f6499.7
Applied rewrites99.7%
if 6.19999999999999982e-30 < a Initial program 33.2%
Taylor expanded in a around inf
unpow3N/A
unpow2N/A
associate-*l*N/A
sub-negN/A
remove-double-negN/A
distribute-neg-inN/A
+-commutativeN/A
sub-negN/A
distribute-rgt-neg-inN/A
*-commutativeN/A
sub-negN/A
+-commutativeN/A
distribute-lft1-inN/A
distribute-lft-neg-outN/A
lft-mult-inverseN/A
metadata-evalN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
Applied rewrites33.1%
Taylor expanded in b around 0
+-commutativeN/A
associate-+l+N/A
+-commutativeN/A
associate--l+N/A
*-commutativeN/A
distribute-lft-inN/A
metadata-evalN/A
unpow2N/A
associate-*r*N/A
lower-fma.f64N/A
Applied rewrites98.4%
Taylor expanded in a around 0
Applied rewrites98.4%
(FPCore (a b) :precision binary64 (if (or (<= a -1600000000.0) (not (<= a 6.2e-30))) (fma (* (* (* a a) 2.0) b) b (fma (* a a) (fma (- a 4.0) a 4.0) -1.0)) (fma (* (fma b b 12.0) b) b -1.0)))
double code(double a, double b) {
double tmp;
if ((a <= -1600000000.0) || !(a <= 6.2e-30)) {
tmp = fma((((a * a) * 2.0) * b), b, fma((a * a), fma((a - 4.0), a, 4.0), -1.0));
} else {
tmp = fma((fma(b, b, 12.0) * b), b, -1.0);
}
return tmp;
}
function code(a, b) tmp = 0.0 if ((a <= -1600000000.0) || !(a <= 6.2e-30)) tmp = fma(Float64(Float64(Float64(a * a) * 2.0) * b), b, fma(Float64(a * a), fma(Float64(a - 4.0), a, 4.0), -1.0)); else tmp = fma(Float64(fma(b, b, 12.0) * b), b, -1.0); end return tmp end
code[a_, b_] := If[Or[LessEqual[a, -1600000000.0], N[Not[LessEqual[a, 6.2e-30]], $MachinePrecision]], N[(N[(N[(N[(a * a), $MachinePrecision] * 2.0), $MachinePrecision] * b), $MachinePrecision] * b + N[(N[(a * a), $MachinePrecision] * N[(N[(a - 4.0), $MachinePrecision] * a + 4.0), $MachinePrecision] + -1.0), $MachinePrecision]), $MachinePrecision], N[(N[(N[(b * b + 12.0), $MachinePrecision] * b), $MachinePrecision] * b + -1.0), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq -1600000000 \lor \neg \left(a \leq 6.2 \cdot 10^{-30}\right):\\
\;\;\;\;\mathsf{fma}\left(\left(\left(a \cdot a\right) \cdot 2\right) \cdot b, b, \mathsf{fma}\left(a \cdot a, \mathsf{fma}\left(a - 4, a, 4\right), -1\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(b, b, 12\right) \cdot b, b, -1\right)\\
\end{array}
\end{array}
if a < -1.6e9 or 6.19999999999999982e-30 < a Initial program 43.6%
Taylor expanded in a around inf
unpow3N/A
unpow2N/A
associate-*l*N/A
sub-negN/A
remove-double-negN/A
distribute-neg-inN/A
+-commutativeN/A
sub-negN/A
distribute-rgt-neg-inN/A
*-commutativeN/A
sub-negN/A
+-commutativeN/A
distribute-lft1-inN/A
distribute-lft-neg-outN/A
lft-mult-inverseN/A
metadata-evalN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
Applied rewrites66.6%
Taylor expanded in b around 0
+-commutativeN/A
associate-+l+N/A
+-commutativeN/A
associate--l+N/A
*-commutativeN/A
distribute-lft-inN/A
metadata-evalN/A
unpow2N/A
associate-*r*N/A
lower-fma.f64N/A
Applied rewrites98.4%
Taylor expanded in a around 0
Applied rewrites98.4%
Taylor expanded in a around inf
Applied rewrites98.4%
if -1.6e9 < a < 6.19999999999999982e-30Initial program 99.9%
Taylor expanded in a around inf
unpow3N/A
unpow2N/A
associate-*l*N/A
sub-negN/A
remove-double-negN/A
distribute-neg-inN/A
+-commutativeN/A
sub-negN/A
distribute-rgt-neg-inN/A
*-commutativeN/A
sub-negN/A
+-commutativeN/A
distribute-lft1-inN/A
distribute-lft-neg-outN/A
lft-mult-inverseN/A
metadata-evalN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
Applied rewrites99.4%
Taylor expanded in a around 0
sub-negN/A
+-commutativeN/A
metadata-evalN/A
metadata-evalN/A
pow-sqrN/A
distribute-rgt-inN/A
+-commutativeN/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
unpow2N/A
lower-fma.f6499.7
Applied rewrites99.7%
Final simplification99.1%
(FPCore (a b)
:precision binary64
(if (<= a -9e+43)
(- (* (* a a) (* a a)) 1.0)
(if (<= a 2.5e+22)
(- (* (* (fma b b (fma 4.0 a 12.0)) b) b) 1.0)
(fma (* (fma 4.0 a 12.0) b) b (fma (* a a) (fma (- a 4.0) a 4.0) -1.0)))))
double code(double a, double b) {
double tmp;
if (a <= -9e+43) {
tmp = ((a * a) * (a * a)) - 1.0;
} else if (a <= 2.5e+22) {
tmp = ((fma(b, b, fma(4.0, a, 12.0)) * b) * b) - 1.0;
} else {
tmp = fma((fma(4.0, a, 12.0) * b), b, fma((a * a), fma((a - 4.0), a, 4.0), -1.0));
}
return tmp;
}
function code(a, b) tmp = 0.0 if (a <= -9e+43) tmp = Float64(Float64(Float64(a * a) * Float64(a * a)) - 1.0); elseif (a <= 2.5e+22) tmp = Float64(Float64(Float64(fma(b, b, fma(4.0, a, 12.0)) * b) * b) - 1.0); else tmp = fma(Float64(fma(4.0, a, 12.0) * b), b, fma(Float64(a * a), fma(Float64(a - 4.0), a, 4.0), -1.0)); end return tmp end
code[a_, b_] := If[LessEqual[a, -9e+43], N[(N[(N[(a * a), $MachinePrecision] * N[(a * a), $MachinePrecision]), $MachinePrecision] - 1.0), $MachinePrecision], If[LessEqual[a, 2.5e+22], N[(N[(N[(N[(b * b + N[(4.0 * a + 12.0), $MachinePrecision]), $MachinePrecision] * b), $MachinePrecision] * b), $MachinePrecision] - 1.0), $MachinePrecision], N[(N[(N[(4.0 * a + 12.0), $MachinePrecision] * b), $MachinePrecision] * b + N[(N[(a * a), $MachinePrecision] * N[(N[(a - 4.0), $MachinePrecision] * a + 4.0), $MachinePrecision] + -1.0), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq -9 \cdot 10^{+43}:\\
\;\;\;\;\left(a \cdot a\right) \cdot \left(a \cdot a\right) - 1\\
\mathbf{elif}\;a \leq 2.5 \cdot 10^{+22}:\\
\;\;\;\;\left(\mathsf{fma}\left(b, b, \mathsf{fma}\left(4, a, 12\right)\right) \cdot b\right) \cdot b - 1\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(4, a, 12\right) \cdot b, b, \mathsf{fma}\left(a \cdot a, \mathsf{fma}\left(a - 4, a, 4\right), -1\right)\right)\\
\end{array}
\end{array}
if a < -9e43Initial program 55.4%
Taylor expanded in a around inf
lower-pow.f6498.3
Applied rewrites98.3%
Applied rewrites98.3%
if -9e43 < a < 2.4999999999999998e22Initial program 97.1%
Taylor expanded in a around 0
associate-+r+N/A
associate-*r*N/A
distribute-rgt-outN/A
metadata-evalN/A
distribute-lft-inN/A
+-commutativeN/A
metadata-evalN/A
pow-sqrN/A
distribute-lft-inN/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
lower-*.f64N/A
Applied rewrites97.7%
if 2.4999999999999998e22 < a Initial program 23.5%
Taylor expanded in a around inf
unpow3N/A
unpow2N/A
associate-*l*N/A
sub-negN/A
remove-double-negN/A
distribute-neg-inN/A
+-commutativeN/A
sub-negN/A
distribute-rgt-neg-inN/A
*-commutativeN/A
sub-negN/A
+-commutativeN/A
distribute-lft1-inN/A
distribute-lft-neg-outN/A
lft-mult-inverseN/A
metadata-evalN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
Applied rewrites23.5%
Taylor expanded in b around 0
+-commutativeN/A
associate-+l+N/A
+-commutativeN/A
associate--l+N/A
*-commutativeN/A
distribute-lft-inN/A
metadata-evalN/A
unpow2N/A
associate-*r*N/A
lower-fma.f64N/A
Applied rewrites99.9%
Taylor expanded in a around 0
Applied rewrites99.9%
Taylor expanded in a around 0
Applied rewrites99.9%
(FPCore (a b) :precision binary64 (if (or (<= a -9e+43) (not (<= a 3.05e+29))) (- (* (* a a) (* a a)) 1.0) (- (* (* (fma b b (fma 4.0 a 12.0)) b) b) 1.0)))
double code(double a, double b) {
double tmp;
if ((a <= -9e+43) || !(a <= 3.05e+29)) {
tmp = ((a * a) * (a * a)) - 1.0;
} else {
tmp = ((fma(b, b, fma(4.0, a, 12.0)) * b) * b) - 1.0;
}
return tmp;
}
function code(a, b) tmp = 0.0 if ((a <= -9e+43) || !(a <= 3.05e+29)) tmp = Float64(Float64(Float64(a * a) * Float64(a * a)) - 1.0); else tmp = Float64(Float64(Float64(fma(b, b, fma(4.0, a, 12.0)) * b) * b) - 1.0); end return tmp end
code[a_, b_] := If[Or[LessEqual[a, -9e+43], N[Not[LessEqual[a, 3.05e+29]], $MachinePrecision]], N[(N[(N[(a * a), $MachinePrecision] * N[(a * a), $MachinePrecision]), $MachinePrecision] - 1.0), $MachinePrecision], N[(N[(N[(N[(b * b + N[(4.0 * a + 12.0), $MachinePrecision]), $MachinePrecision] * b), $MachinePrecision] * b), $MachinePrecision] - 1.0), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq -9 \cdot 10^{+43} \lor \neg \left(a \leq 3.05 \cdot 10^{+29}\right):\\
\;\;\;\;\left(a \cdot a\right) \cdot \left(a \cdot a\right) - 1\\
\mathbf{else}:\\
\;\;\;\;\left(\mathsf{fma}\left(b, b, \mathsf{fma}\left(4, a, 12\right)\right) \cdot b\right) \cdot b - 1\\
\end{array}
\end{array}
if a < -9e43 or 3.0499999999999999e29 < a Initial program 39.6%
Taylor expanded in a around inf
lower-pow.f6498.3
Applied rewrites98.3%
Applied rewrites98.2%
if -9e43 < a < 3.0499999999999999e29Initial program 97.1%
Taylor expanded in a around 0
associate-+r+N/A
associate-*r*N/A
distribute-rgt-outN/A
metadata-evalN/A
distribute-lft-inN/A
+-commutativeN/A
metadata-evalN/A
pow-sqrN/A
distribute-lft-inN/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
lower-*.f64N/A
Applied rewrites97.7%
Final simplification97.9%
(FPCore (a b) :precision binary64 (if (or (<= a -9e+43) (not (<= a 3.05e+29))) (- (* (* a a) (* a a)) 1.0) (fma (* (fma b b 12.0) b) b -1.0)))
double code(double a, double b) {
double tmp;
if ((a <= -9e+43) || !(a <= 3.05e+29)) {
tmp = ((a * a) * (a * a)) - 1.0;
} else {
tmp = fma((fma(b, b, 12.0) * b), b, -1.0);
}
return tmp;
}
function code(a, b) tmp = 0.0 if ((a <= -9e+43) || !(a <= 3.05e+29)) tmp = Float64(Float64(Float64(a * a) * Float64(a * a)) - 1.0); else tmp = fma(Float64(fma(b, b, 12.0) * b), b, -1.0); end return tmp end
code[a_, b_] := If[Or[LessEqual[a, -9e+43], N[Not[LessEqual[a, 3.05e+29]], $MachinePrecision]], N[(N[(N[(a * a), $MachinePrecision] * N[(a * a), $MachinePrecision]), $MachinePrecision] - 1.0), $MachinePrecision], N[(N[(N[(b * b + 12.0), $MachinePrecision] * b), $MachinePrecision] * b + -1.0), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq -9 \cdot 10^{+43} \lor \neg \left(a \leq 3.05 \cdot 10^{+29}\right):\\
\;\;\;\;\left(a \cdot a\right) \cdot \left(a \cdot a\right) - 1\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(b, b, 12\right) \cdot b, b, -1\right)\\
\end{array}
\end{array}
if a < -9e43 or 3.0499999999999999e29 < a Initial program 39.6%
Taylor expanded in a around inf
lower-pow.f6498.3
Applied rewrites98.3%
Applied rewrites98.2%
if -9e43 < a < 3.0499999999999999e29Initial program 97.1%
Taylor expanded in a around inf
unpow3N/A
unpow2N/A
associate-*l*N/A
sub-negN/A
remove-double-negN/A
distribute-neg-inN/A
+-commutativeN/A
sub-negN/A
distribute-rgt-neg-inN/A
*-commutativeN/A
sub-negN/A
+-commutativeN/A
distribute-lft1-inN/A
distribute-lft-neg-outN/A
lft-mult-inverseN/A
metadata-evalN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
Applied rewrites99.4%
Taylor expanded in a around 0
sub-negN/A
+-commutativeN/A
metadata-evalN/A
metadata-evalN/A
pow-sqrN/A
distribute-rgt-inN/A
+-commutativeN/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
unpow2N/A
lower-fma.f6497.7
Applied rewrites97.7%
Final simplification97.9%
(FPCore (a b)
:precision binary64
(if (<= a -4.3e+101)
(fma (* a a) (fma -4.0 a 4.0) -1.0)
(if (<= a 6.6e+153)
(fma (* (fma b b 12.0) b) b -1.0)
(fma (* a a) 4.0 -1.0))))
double code(double a, double b) {
double tmp;
if (a <= -4.3e+101) {
tmp = fma((a * a), fma(-4.0, a, 4.0), -1.0);
} else if (a <= 6.6e+153) {
tmp = fma((fma(b, b, 12.0) * b), b, -1.0);
} else {
tmp = fma((a * a), 4.0, -1.0);
}
return tmp;
}
function code(a, b) tmp = 0.0 if (a <= -4.3e+101) tmp = fma(Float64(a * a), fma(-4.0, a, 4.0), -1.0); elseif (a <= 6.6e+153) tmp = fma(Float64(fma(b, b, 12.0) * b), b, -1.0); else tmp = fma(Float64(a * a), 4.0, -1.0); end return tmp end
code[a_, b_] := If[LessEqual[a, -4.3e+101], N[(N[(a * a), $MachinePrecision] * N[(-4.0 * a + 4.0), $MachinePrecision] + -1.0), $MachinePrecision], If[LessEqual[a, 6.6e+153], N[(N[(N[(b * b + 12.0), $MachinePrecision] * b), $MachinePrecision] * b + -1.0), $MachinePrecision], N[(N[(a * a), $MachinePrecision] * 4.0 + -1.0), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq -4.3 \cdot 10^{+101}:\\
\;\;\;\;\mathsf{fma}\left(a \cdot a, \mathsf{fma}\left(-4, a, 4\right), -1\right)\\
\mathbf{elif}\;a \leq 6.6 \cdot 10^{+153}:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(b, b, 12\right) \cdot b, b, -1\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(a \cdot a, 4, -1\right)\\
\end{array}
\end{array}
if a < -4.3000000000000001e101Initial program 54.3%
Taylor expanded in a around inf
unpow3N/A
unpow2N/A
associate-*l*N/A
sub-negN/A
remove-double-negN/A
distribute-neg-inN/A
+-commutativeN/A
sub-negN/A
distribute-rgt-neg-inN/A
*-commutativeN/A
sub-negN/A
+-commutativeN/A
distribute-lft1-inN/A
distribute-lft-neg-outN/A
lft-mult-inverseN/A
metadata-evalN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
Applied rewrites100.0%
Taylor expanded in b around 0
sub-negN/A
*-commutativeN/A
associate-*r*N/A
metadata-evalN/A
pow-sqrN/A
distribute-rgt-outN/A
metadata-evalN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
lower-fma.f64N/A
lower--.f64N/A
unpow2N/A
lower-*.f64100.0
Applied rewrites100.0%
Taylor expanded in a around 0
Applied rewrites98.1%
if -4.3000000000000001e101 < a < 6.59999999999999989e153Initial program 88.7%
Taylor expanded in a around inf
unpow3N/A
unpow2N/A
associate-*l*N/A
sub-negN/A
remove-double-negN/A
distribute-neg-inN/A
+-commutativeN/A
sub-negN/A
distribute-rgt-neg-inN/A
*-commutativeN/A
sub-negN/A
+-commutativeN/A
distribute-lft1-inN/A
distribute-lft-neg-outN/A
lft-mult-inverseN/A
metadata-evalN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
Applied rewrites92.8%
Taylor expanded in a around 0
sub-negN/A
+-commutativeN/A
metadata-evalN/A
metadata-evalN/A
pow-sqrN/A
distribute-rgt-inN/A
+-commutativeN/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
unpow2N/A
lower-fma.f6485.1
Applied rewrites85.1%
if 6.59999999999999989e153 < a Initial program 0.0%
Taylor expanded in a around inf
unpow3N/A
unpow2N/A
associate-*l*N/A
sub-negN/A
remove-double-negN/A
distribute-neg-inN/A
+-commutativeN/A
sub-negN/A
distribute-rgt-neg-inN/A
*-commutativeN/A
sub-negN/A
+-commutativeN/A
distribute-lft1-inN/A
distribute-lft-neg-outN/A
lft-mult-inverseN/A
metadata-evalN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
Applied rewrites0.0%
Taylor expanded in b around 0
sub-negN/A
*-commutativeN/A
associate-*r*N/A
metadata-evalN/A
pow-sqrN/A
distribute-rgt-outN/A
metadata-evalN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
lower-fma.f64N/A
lower--.f64N/A
unpow2N/A
lower-*.f64100.0
Applied rewrites100.0%
Taylor expanded in a around 0
Applied rewrites100.0%
(FPCore (a b) :precision binary64 (if (or (<= a -5.6e+126) (not (<= a 6.6e+153))) (fma (* a a) 4.0 -1.0) (fma (* (* b b) b) b -1.0)))
double code(double a, double b) {
double tmp;
if ((a <= -5.6e+126) || !(a <= 6.6e+153)) {
tmp = fma((a * a), 4.0, -1.0);
} else {
tmp = fma(((b * b) * b), b, -1.0);
}
return tmp;
}
function code(a, b) tmp = 0.0 if ((a <= -5.6e+126) || !(a <= 6.6e+153)) tmp = fma(Float64(a * a), 4.0, -1.0); else tmp = fma(Float64(Float64(b * b) * b), b, -1.0); end return tmp end
code[a_, b_] := If[Or[LessEqual[a, -5.6e+126], N[Not[LessEqual[a, 6.6e+153]], $MachinePrecision]], N[(N[(a * a), $MachinePrecision] * 4.0 + -1.0), $MachinePrecision], N[(N[(N[(b * b), $MachinePrecision] * b), $MachinePrecision] * b + -1.0), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq -5.6 \cdot 10^{+126} \lor \neg \left(a \leq 6.6 \cdot 10^{+153}\right):\\
\;\;\;\;\mathsf{fma}\left(a \cdot a, 4, -1\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\left(b \cdot b\right) \cdot b, b, -1\right)\\
\end{array}
\end{array}
if a < -5.60000000000000018e126 or 6.59999999999999989e153 < a Initial program 29.4%
Taylor expanded in a around inf
unpow3N/A
unpow2N/A
associate-*l*N/A
sub-negN/A
remove-double-negN/A
distribute-neg-inN/A
+-commutativeN/A
sub-negN/A
distribute-rgt-neg-inN/A
*-commutativeN/A
sub-negN/A
+-commutativeN/A
distribute-lft1-inN/A
distribute-lft-neg-outN/A
lft-mult-inverseN/A
metadata-evalN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
Applied rewrites55.9%
Taylor expanded in b around 0
sub-negN/A
*-commutativeN/A
associate-*r*N/A
metadata-evalN/A
pow-sqrN/A
distribute-rgt-outN/A
metadata-evalN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
lower-fma.f64N/A
lower--.f64N/A
unpow2N/A
lower-*.f64100.0
Applied rewrites100.0%
Taylor expanded in a around 0
Applied rewrites92.0%
if -5.60000000000000018e126 < a < 6.59999999999999989e153Initial program 87.6%
Taylor expanded in a around inf
unpow3N/A
unpow2N/A
associate-*l*N/A
sub-negN/A
remove-double-negN/A
distribute-neg-inN/A
+-commutativeN/A
sub-negN/A
distribute-rgt-neg-inN/A
*-commutativeN/A
sub-negN/A
+-commutativeN/A
distribute-lft1-inN/A
distribute-lft-neg-outN/A
lft-mult-inverseN/A
metadata-evalN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
Applied rewrites93.1%
Taylor expanded in a around 0
sub-negN/A
+-commutativeN/A
metadata-evalN/A
metadata-evalN/A
pow-sqrN/A
distribute-rgt-inN/A
+-commutativeN/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
unpow2N/A
lower-fma.f6483.1
Applied rewrites83.1%
Taylor expanded in b around 0
Applied rewrites65.2%
Taylor expanded in b around inf
Applied rewrites82.7%
Final simplification85.2%
(FPCore (a b) :precision binary64 (if (<= a -4.3e+101) (fma (* a a) (fma -4.0 a 4.0) -1.0) (if (<= a 6.6e+153) (fma (* (* b b) b) b -1.0) (fma (* a a) 4.0 -1.0))))
double code(double a, double b) {
double tmp;
if (a <= -4.3e+101) {
tmp = fma((a * a), fma(-4.0, a, 4.0), -1.0);
} else if (a <= 6.6e+153) {
tmp = fma(((b * b) * b), b, -1.0);
} else {
tmp = fma((a * a), 4.0, -1.0);
}
return tmp;
}
function code(a, b) tmp = 0.0 if (a <= -4.3e+101) tmp = fma(Float64(a * a), fma(-4.0, a, 4.0), -1.0); elseif (a <= 6.6e+153) tmp = fma(Float64(Float64(b * b) * b), b, -1.0); else tmp = fma(Float64(a * a), 4.0, -1.0); end return tmp end
code[a_, b_] := If[LessEqual[a, -4.3e+101], N[(N[(a * a), $MachinePrecision] * N[(-4.0 * a + 4.0), $MachinePrecision] + -1.0), $MachinePrecision], If[LessEqual[a, 6.6e+153], N[(N[(N[(b * b), $MachinePrecision] * b), $MachinePrecision] * b + -1.0), $MachinePrecision], N[(N[(a * a), $MachinePrecision] * 4.0 + -1.0), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq -4.3 \cdot 10^{+101}:\\
\;\;\;\;\mathsf{fma}\left(a \cdot a, \mathsf{fma}\left(-4, a, 4\right), -1\right)\\
\mathbf{elif}\;a \leq 6.6 \cdot 10^{+153}:\\
\;\;\;\;\mathsf{fma}\left(\left(b \cdot b\right) \cdot b, b, -1\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(a \cdot a, 4, -1\right)\\
\end{array}
\end{array}
if a < -4.3000000000000001e101Initial program 54.3%
Taylor expanded in a around inf
unpow3N/A
unpow2N/A
associate-*l*N/A
sub-negN/A
remove-double-negN/A
distribute-neg-inN/A
+-commutativeN/A
sub-negN/A
distribute-rgt-neg-inN/A
*-commutativeN/A
sub-negN/A
+-commutativeN/A
distribute-lft1-inN/A
distribute-lft-neg-outN/A
lft-mult-inverseN/A
metadata-evalN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
Applied rewrites100.0%
Taylor expanded in b around 0
sub-negN/A
*-commutativeN/A
associate-*r*N/A
metadata-evalN/A
pow-sqrN/A
distribute-rgt-outN/A
metadata-evalN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
lower-fma.f64N/A
lower--.f64N/A
unpow2N/A
lower-*.f64100.0
Applied rewrites100.0%
Taylor expanded in a around 0
Applied rewrites98.1%
if -4.3000000000000001e101 < a < 6.59999999999999989e153Initial program 88.7%
Taylor expanded in a around inf
unpow3N/A
unpow2N/A
associate-*l*N/A
sub-negN/A
remove-double-negN/A
distribute-neg-inN/A
+-commutativeN/A
sub-negN/A
distribute-rgt-neg-inN/A
*-commutativeN/A
sub-negN/A
+-commutativeN/A
distribute-lft1-inN/A
distribute-lft-neg-outN/A
lft-mult-inverseN/A
metadata-evalN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
Applied rewrites92.8%
Taylor expanded in a around 0
sub-negN/A
+-commutativeN/A
metadata-evalN/A
metadata-evalN/A
pow-sqrN/A
distribute-rgt-inN/A
+-commutativeN/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
unpow2N/A
lower-fma.f6485.1
Applied rewrites85.1%
Taylor expanded in b around 0
Applied rewrites66.4%
Taylor expanded in b around inf
Applied rewrites84.7%
if 6.59999999999999989e153 < a Initial program 0.0%
Taylor expanded in a around inf
unpow3N/A
unpow2N/A
associate-*l*N/A
sub-negN/A
remove-double-negN/A
distribute-neg-inN/A
+-commutativeN/A
sub-negN/A
distribute-rgt-neg-inN/A
*-commutativeN/A
sub-negN/A
+-commutativeN/A
distribute-lft1-inN/A
distribute-lft-neg-outN/A
lft-mult-inverseN/A
metadata-evalN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
Applied rewrites0.0%
Taylor expanded in b around 0
sub-negN/A
*-commutativeN/A
associate-*r*N/A
metadata-evalN/A
pow-sqrN/A
distribute-rgt-outN/A
metadata-evalN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
lower-fma.f64N/A
lower--.f64N/A
unpow2N/A
lower-*.f64100.0
Applied rewrites100.0%
Taylor expanded in a around 0
Applied rewrites100.0%
(FPCore (a b) :precision binary64 (if (<= (* b b) 2e+296) (fma (* a a) 4.0 -1.0) (fma (* b b) 12.0 -1.0)))
double code(double a, double b) {
double tmp;
if ((b * b) <= 2e+296) {
tmp = fma((a * a), 4.0, -1.0);
} else {
tmp = fma((b * b), 12.0, -1.0);
}
return tmp;
}
function code(a, b) tmp = 0.0 if (Float64(b * b) <= 2e+296) tmp = fma(Float64(a * a), 4.0, -1.0); else tmp = fma(Float64(b * b), 12.0, -1.0); end return tmp end
code[a_, b_] := If[LessEqual[N[(b * b), $MachinePrecision], 2e+296], N[(N[(a * a), $MachinePrecision] * 4.0 + -1.0), $MachinePrecision], N[(N[(b * b), $MachinePrecision] * 12.0 + -1.0), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \cdot b \leq 2 \cdot 10^{+296}:\\
\;\;\;\;\mathsf{fma}\left(a \cdot a, 4, -1\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(b \cdot b, 12, -1\right)\\
\end{array}
\end{array}
if (*.f64 b b) < 1.99999999999999996e296Initial program 77.6%
Taylor expanded in a around inf
unpow3N/A
unpow2N/A
associate-*l*N/A
sub-negN/A
remove-double-negN/A
distribute-neg-inN/A
+-commutativeN/A
sub-negN/A
distribute-rgt-neg-inN/A
*-commutativeN/A
sub-negN/A
+-commutativeN/A
distribute-lft1-inN/A
distribute-lft-neg-outN/A
lft-mult-inverseN/A
metadata-evalN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
Applied rewrites82.0%
Taylor expanded in b around 0
sub-negN/A
*-commutativeN/A
associate-*r*N/A
metadata-evalN/A
pow-sqrN/A
distribute-rgt-outN/A
metadata-evalN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
lower-fma.f64N/A
lower--.f64N/A
unpow2N/A
lower-*.f6482.9
Applied rewrites82.9%
Taylor expanded in a around 0
Applied rewrites63.3%
if 1.99999999999999996e296 < (*.f64 b b) Initial program 56.7%
Taylor expanded in a around inf
unpow3N/A
unpow2N/A
associate-*l*N/A
sub-negN/A
remove-double-negN/A
distribute-neg-inN/A
+-commutativeN/A
sub-negN/A
distribute-rgt-neg-inN/A
*-commutativeN/A
sub-negN/A
+-commutativeN/A
distribute-lft1-inN/A
distribute-lft-neg-outN/A
lft-mult-inverseN/A
metadata-evalN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
Applied rewrites86.6%
Taylor expanded in b around 0
+-commutativeN/A
associate-+l+N/A
+-commutativeN/A
associate--l+N/A
*-commutativeN/A
distribute-lft-inN/A
metadata-evalN/A
unpow2N/A
associate-*r*N/A
lower-fma.f64N/A
Applied rewrites98.7%
Taylor expanded in a around 0
Applied rewrites97.4%
(FPCore (a b) :precision binary64 (fma (* b b) 12.0 -1.0))
double code(double a, double b) {
return fma((b * b), 12.0, -1.0);
}
function code(a, b) return fma(Float64(b * b), 12.0, -1.0) end
code[a_, b_] := N[(N[(b * b), $MachinePrecision] * 12.0 + -1.0), $MachinePrecision]
\begin{array}{l}
\\
\mathsf{fma}\left(b \cdot b, 12, -1\right)
\end{array}
Initial program 72.1%
Taylor expanded in a around inf
unpow3N/A
unpow2N/A
associate-*l*N/A
sub-negN/A
remove-double-negN/A
distribute-neg-inN/A
+-commutativeN/A
sub-negN/A
distribute-rgt-neg-inN/A
*-commutativeN/A
sub-negN/A
+-commutativeN/A
distribute-lft1-inN/A
distribute-lft-neg-outN/A
lft-mult-inverseN/A
metadata-evalN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
Applied rewrites83.2%
Taylor expanded in b around 0
+-commutativeN/A
associate-+l+N/A
+-commutativeN/A
associate--l+N/A
*-commutativeN/A
distribute-lft-inN/A
metadata-evalN/A
unpow2N/A
associate-*r*N/A
lower-fma.f64N/A
Applied rewrites87.9%
Taylor expanded in a around 0
Applied rewrites54.3%
(FPCore (a b) :precision binary64 -1.0)
double code(double a, double b) {
return -1.0;
}
real(8) function code(a, b)
real(8), intent (in) :: a
real(8), intent (in) :: b
code = -1.0d0
end function
public static double code(double a, double b) {
return -1.0;
}
def code(a, b): return -1.0
function code(a, b) return -1.0 end
function tmp = code(a, b) tmp = -1.0; end
code[a_, b_] := -1.0
\begin{array}{l}
\\
-1
\end{array}
Initial program 72.1%
Taylor expanded in a around inf
unpow3N/A
unpow2N/A
associate-*l*N/A
sub-negN/A
remove-double-negN/A
distribute-neg-inN/A
+-commutativeN/A
sub-negN/A
distribute-rgt-neg-inN/A
*-commutativeN/A
sub-negN/A
+-commutativeN/A
distribute-lft1-inN/A
distribute-lft-neg-outN/A
lft-mult-inverseN/A
metadata-evalN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
Applied rewrites83.2%
Taylor expanded in b around 0
sub-negN/A
*-commutativeN/A
associate-*r*N/A
metadata-evalN/A
pow-sqrN/A
distribute-rgt-outN/A
metadata-evalN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
lower-fma.f64N/A
lower--.f64N/A
unpow2N/A
lower-*.f6471.2
Applied rewrites71.2%
Taylor expanded in a around 0
Applied rewrites27.6%
herbie shell --seed 2024314
(FPCore (a b)
:name "Bouland and Aaronson, Equation (24)"
:precision binary64
(- (+ (pow (+ (* a a) (* b b)) 2.0) (* 4.0 (+ (* (* a a) (- 1.0 a)) (* (* b b) (+ 3.0 a))))) 1.0))