
(FPCore (a b) :precision binary64 (- (+ (pow (+ (* a a) (* b b)) 2.0) (* 4.0 (+ (* (* a a) (+ 1.0 a)) (* (* b b) (- 1.0 (* 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) * (1.0 - (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) * (1.0d0 - (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) * (1.0 - (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) * (1.0 - (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(1.0 - 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) * (1.0 - (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[(1.0 - N[(3.0 * a), $MachinePrecision]), $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(1 - 3 \cdot 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) (- 1.0 (* 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) * (1.0 - (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) * (1.0d0 - (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) * (1.0 - (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) * (1.0 - (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(1.0 - 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) * (1.0 - (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[(1.0 - N[(3.0 * a), $MachinePrecision]), $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(1 - 3 \cdot a\right)\right)\right) - 1
\end{array}
(FPCore (a b) :precision binary64 (let* ((t_0 (+ (* a a) (* b b)))) (+ (/ t_0 (/ 1.0 t_0)) (+ (* (* b b) 4.0) -1.0))))
double code(double a, double b) {
double t_0 = (a * a) + (b * b);
return (t_0 / (1.0 / t_0)) + (((b * b) * 4.0) + -1.0);
}
real(8) function code(a, b)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: t_0
t_0 = (a * a) + (b * b)
code = (t_0 / (1.0d0 / t_0)) + (((b * b) * 4.0d0) + (-1.0d0))
end function
public static double code(double a, double b) {
double t_0 = (a * a) + (b * b);
return (t_0 / (1.0 / t_0)) + (((b * b) * 4.0) + -1.0);
}
def code(a, b): t_0 = (a * a) + (b * b) return (t_0 / (1.0 / t_0)) + (((b * b) * 4.0) + -1.0)
function code(a, b) t_0 = Float64(Float64(a * a) + Float64(b * b)) return Float64(Float64(t_0 / Float64(1.0 / t_0)) + Float64(Float64(Float64(b * b) * 4.0) + -1.0)) end
function tmp = code(a, b) t_0 = (a * a) + (b * b); tmp = (t_0 / (1.0 / t_0)) + (((b * b) * 4.0) + -1.0); end
code[a_, b_] := Block[{t$95$0 = N[(N[(a * a), $MachinePrecision] + N[(b * b), $MachinePrecision]), $MachinePrecision]}, N[(N[(t$95$0 / N[(1.0 / t$95$0), $MachinePrecision]), $MachinePrecision] + N[(N[(N[(b * b), $MachinePrecision] * 4.0), $MachinePrecision] + -1.0), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := a \cdot a + b \cdot b\\
\frac{t\_0}{\frac{1}{t\_0}} + \left(\left(b \cdot b\right) \cdot 4 + -1\right)
\end{array}
\end{array}
Initial program 73.4%
associate--l+N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
Simplified73.4%
Taylor expanded in a around 0
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6498.8%
Simplified98.8%
+-commutativeN/A
flip-+N/A
associate-*r*N/A
associate-*r*N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
clear-numN/A
metadata-evalN/A
Applied egg-rr98.9%
Final simplification98.9%
(FPCore (a b)
:precision binary64
(let* ((t_0 (* a (* a (* a a)))))
(if (<= b 8e-282)
t_0
(if (<= b 2.65e-235)
-1.0
(if (<= b 3.05e-130)
(* (* a a) (* a a))
(if (<= b 0.000195)
-1.0
(if (<= b 4e+73) t_0 (* b (* b (* b b))))))))))
double code(double a, double b) {
double t_0 = a * (a * (a * a));
double tmp;
if (b <= 8e-282) {
tmp = t_0;
} else if (b <= 2.65e-235) {
tmp = -1.0;
} else if (b <= 3.05e-130) {
tmp = (a * a) * (a * a);
} else if (b <= 0.000195) {
tmp = -1.0;
} else if (b <= 4e+73) {
tmp = t_0;
} else {
tmp = b * (b * (b * b));
}
return tmp;
}
real(8) function code(a, b)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: t_0
real(8) :: tmp
t_0 = a * (a * (a * a))
if (b <= 8d-282) then
tmp = t_0
else if (b <= 2.65d-235) then
tmp = -1.0d0
else if (b <= 3.05d-130) then
tmp = (a * a) * (a * a)
else if (b <= 0.000195d0) then
tmp = -1.0d0
else if (b <= 4d+73) then
tmp = t_0
else
tmp = b * (b * (b * b))
end if
code = tmp
end function
public static double code(double a, double b) {
double t_0 = a * (a * (a * a));
double tmp;
if (b <= 8e-282) {
tmp = t_0;
} else if (b <= 2.65e-235) {
tmp = -1.0;
} else if (b <= 3.05e-130) {
tmp = (a * a) * (a * a);
} else if (b <= 0.000195) {
tmp = -1.0;
} else if (b <= 4e+73) {
tmp = t_0;
} else {
tmp = b * (b * (b * b));
}
return tmp;
}
def code(a, b): t_0 = a * (a * (a * a)) tmp = 0 if b <= 8e-282: tmp = t_0 elif b <= 2.65e-235: tmp = -1.0 elif b <= 3.05e-130: tmp = (a * a) * (a * a) elif b <= 0.000195: tmp = -1.0 elif b <= 4e+73: tmp = t_0 else: tmp = b * (b * (b * b)) return tmp
function code(a, b) t_0 = Float64(a * Float64(a * Float64(a * a))) tmp = 0.0 if (b <= 8e-282) tmp = t_0; elseif (b <= 2.65e-235) tmp = -1.0; elseif (b <= 3.05e-130) tmp = Float64(Float64(a * a) * Float64(a * a)); elseif (b <= 0.000195) tmp = -1.0; elseif (b <= 4e+73) tmp = t_0; else tmp = Float64(b * Float64(b * Float64(b * b))); end return tmp end
function tmp_2 = code(a, b) t_0 = a * (a * (a * a)); tmp = 0.0; if (b <= 8e-282) tmp = t_0; elseif (b <= 2.65e-235) tmp = -1.0; elseif (b <= 3.05e-130) tmp = (a * a) * (a * a); elseif (b <= 0.000195) tmp = -1.0; elseif (b <= 4e+73) tmp = t_0; else tmp = b * (b * (b * b)); end tmp_2 = tmp; end
code[a_, b_] := Block[{t$95$0 = N[(a * N[(a * N[(a * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[b, 8e-282], t$95$0, If[LessEqual[b, 2.65e-235], -1.0, If[LessEqual[b, 3.05e-130], N[(N[(a * a), $MachinePrecision] * N[(a * a), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 0.000195], -1.0, If[LessEqual[b, 4e+73], t$95$0, N[(b * N[(b * N[(b * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := a \cdot \left(a \cdot \left(a \cdot a\right)\right)\\
\mathbf{if}\;b \leq 8 \cdot 10^{-282}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;b \leq 2.65 \cdot 10^{-235}:\\
\;\;\;\;-1\\
\mathbf{elif}\;b \leq 3.05 \cdot 10^{-130}:\\
\;\;\;\;\left(a \cdot a\right) \cdot \left(a \cdot a\right)\\
\mathbf{elif}\;b \leq 0.000195:\\
\;\;\;\;-1\\
\mathbf{elif}\;b \leq 4 \cdot 10^{+73}:\\
\;\;\;\;t\_0\\
\mathbf{else}:\\
\;\;\;\;b \cdot \left(b \cdot \left(b \cdot b\right)\right)\\
\end{array}
\end{array}
if b < 8.0000000000000001e-282 or 1.94999999999999996e-4 < b < 3.99999999999999993e73Initial program 70.9%
associate--l+N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
Simplified70.9%
Taylor expanded in a around 0
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6498.8%
Simplified98.8%
Taylor expanded in a around inf
metadata-evalN/A
pow-sqrN/A
unpow2N/A
associate-*l*N/A
unpow2N/A
cube-multN/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6446.9%
Simplified46.9%
if 8.0000000000000001e-282 < b < 2.6500000000000001e-235 or 3.04999999999999998e-130 < b < 1.94999999999999996e-4Initial program 99.9%
associate--l+N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
Simplified99.9%
Taylor expanded in a around 0
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
+-commutativeN/A
metadata-evalN/A
pow-sqrN/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f6475.2%
Simplified75.2%
Taylor expanded in b around 0
Simplified73.7%
if 2.6500000000000001e-235 < b < 3.04999999999999998e-130Initial program 78.3%
associate--l+N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
Simplified78.3%
Taylor expanded in a around inf
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6470.4%
Simplified70.4%
if 3.99999999999999993e73 < b Initial program 59.6%
associate--l+N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
Simplified59.6%
Taylor expanded in b around inf
metadata-evalN/A
pow-sqrN/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64100.0%
Simplified100.0%
(FPCore (a b)
:precision binary64
(let* ((t_0 (* a (* a (* a a)))))
(if (<= b 9e-281)
t_0
(if (<= b 4.5e-235)
-1.0
(if (<= b 2.1e-130)
t_0
(if (<= b 4.3e-5)
-1.0
(if (<= b 7.5e+73) t_0 (* b (* b (* b b))))))))))
double code(double a, double b) {
double t_0 = a * (a * (a * a));
double tmp;
if (b <= 9e-281) {
tmp = t_0;
} else if (b <= 4.5e-235) {
tmp = -1.0;
} else if (b <= 2.1e-130) {
tmp = t_0;
} else if (b <= 4.3e-5) {
tmp = -1.0;
} else if (b <= 7.5e+73) {
tmp = t_0;
} else {
tmp = b * (b * (b * b));
}
return tmp;
}
real(8) function code(a, b)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: t_0
real(8) :: tmp
t_0 = a * (a * (a * a))
if (b <= 9d-281) then
tmp = t_0
else if (b <= 4.5d-235) then
tmp = -1.0d0
else if (b <= 2.1d-130) then
tmp = t_0
else if (b <= 4.3d-5) then
tmp = -1.0d0
else if (b <= 7.5d+73) then
tmp = t_0
else
tmp = b * (b * (b * b))
end if
code = tmp
end function
public static double code(double a, double b) {
double t_0 = a * (a * (a * a));
double tmp;
if (b <= 9e-281) {
tmp = t_0;
} else if (b <= 4.5e-235) {
tmp = -1.0;
} else if (b <= 2.1e-130) {
tmp = t_0;
} else if (b <= 4.3e-5) {
tmp = -1.0;
} else if (b <= 7.5e+73) {
tmp = t_0;
} else {
tmp = b * (b * (b * b));
}
return tmp;
}
def code(a, b): t_0 = a * (a * (a * a)) tmp = 0 if b <= 9e-281: tmp = t_0 elif b <= 4.5e-235: tmp = -1.0 elif b <= 2.1e-130: tmp = t_0 elif b <= 4.3e-5: tmp = -1.0 elif b <= 7.5e+73: tmp = t_0 else: tmp = b * (b * (b * b)) return tmp
function code(a, b) t_0 = Float64(a * Float64(a * Float64(a * a))) tmp = 0.0 if (b <= 9e-281) tmp = t_0; elseif (b <= 4.5e-235) tmp = -1.0; elseif (b <= 2.1e-130) tmp = t_0; elseif (b <= 4.3e-5) tmp = -1.0; elseif (b <= 7.5e+73) tmp = t_0; else tmp = Float64(b * Float64(b * Float64(b * b))); end return tmp end
function tmp_2 = code(a, b) t_0 = a * (a * (a * a)); tmp = 0.0; if (b <= 9e-281) tmp = t_0; elseif (b <= 4.5e-235) tmp = -1.0; elseif (b <= 2.1e-130) tmp = t_0; elseif (b <= 4.3e-5) tmp = -1.0; elseif (b <= 7.5e+73) tmp = t_0; else tmp = b * (b * (b * b)); end tmp_2 = tmp; end
code[a_, b_] := Block[{t$95$0 = N[(a * N[(a * N[(a * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[b, 9e-281], t$95$0, If[LessEqual[b, 4.5e-235], -1.0, If[LessEqual[b, 2.1e-130], t$95$0, If[LessEqual[b, 4.3e-5], -1.0, If[LessEqual[b, 7.5e+73], t$95$0, N[(b * N[(b * N[(b * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := a \cdot \left(a \cdot \left(a \cdot a\right)\right)\\
\mathbf{if}\;b \leq 9 \cdot 10^{-281}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;b \leq 4.5 \cdot 10^{-235}:\\
\;\;\;\;-1\\
\mathbf{elif}\;b \leq 2.1 \cdot 10^{-130}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;b \leq 4.3 \cdot 10^{-5}:\\
\;\;\;\;-1\\
\mathbf{elif}\;b \leq 7.5 \cdot 10^{+73}:\\
\;\;\;\;t\_0\\
\mathbf{else}:\\
\;\;\;\;b \cdot \left(b \cdot \left(b \cdot b\right)\right)\\
\end{array}
\end{array}
if b < 8.99999999999999986e-281 or 4.4999999999999998e-235 < b < 2.10000000000000002e-130 or 4.3000000000000002e-5 < b < 7.5e73Initial program 71.9%
associate--l+N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
Simplified71.9%
Taylor expanded in a around 0
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6499.0%
Simplified99.0%
Taylor expanded in a around inf
metadata-evalN/A
pow-sqrN/A
unpow2N/A
associate-*l*N/A
unpow2N/A
cube-multN/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6450.1%
Simplified50.1%
if 8.99999999999999986e-281 < b < 4.4999999999999998e-235 or 2.10000000000000002e-130 < b < 4.3000000000000002e-5Initial program 99.9%
associate--l+N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
Simplified99.9%
Taylor expanded in a around 0
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
+-commutativeN/A
metadata-evalN/A
pow-sqrN/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f6475.2%
Simplified75.2%
Taylor expanded in b around 0
Simplified73.7%
if 7.5e73 < b Initial program 59.6%
associate--l+N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
Simplified59.6%
Taylor expanded in b around inf
metadata-evalN/A
pow-sqrN/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64100.0%
Simplified100.0%
(FPCore (a b) :precision binary64 (if (<= (* b b) 0.0005) (+ -1.0 (* a (* a (+ (* a a) (+ 4.0 (* a 4.0)))))) (+ (+ (* (* b b) 4.0) -1.0) (/ (+ (* a a) (* b b)) (/ 1.0 (* b b))))))
double code(double a, double b) {
double tmp;
if ((b * b) <= 0.0005) {
tmp = -1.0 + (a * (a * ((a * a) + (4.0 + (a * 4.0)))));
} else {
tmp = (((b * b) * 4.0) + -1.0) + (((a * a) + (b * b)) / (1.0 / (b * b)));
}
return tmp;
}
real(8) function code(a, b)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: tmp
if ((b * b) <= 0.0005d0) then
tmp = (-1.0d0) + (a * (a * ((a * a) + (4.0d0 + (a * 4.0d0)))))
else
tmp = (((b * b) * 4.0d0) + (-1.0d0)) + (((a * a) + (b * b)) / (1.0d0 / (b * b)))
end if
code = tmp
end function
public static double code(double a, double b) {
double tmp;
if ((b * b) <= 0.0005) {
tmp = -1.0 + (a * (a * ((a * a) + (4.0 + (a * 4.0)))));
} else {
tmp = (((b * b) * 4.0) + -1.0) + (((a * a) + (b * b)) / (1.0 / (b * b)));
}
return tmp;
}
def code(a, b): tmp = 0 if (b * b) <= 0.0005: tmp = -1.0 + (a * (a * ((a * a) + (4.0 + (a * 4.0))))) else: tmp = (((b * b) * 4.0) + -1.0) + (((a * a) + (b * b)) / (1.0 / (b * b))) return tmp
function code(a, b) tmp = 0.0 if (Float64(b * b) <= 0.0005) tmp = Float64(-1.0 + Float64(a * Float64(a * Float64(Float64(a * a) + Float64(4.0 + Float64(a * 4.0)))))); else tmp = Float64(Float64(Float64(Float64(b * b) * 4.0) + -1.0) + Float64(Float64(Float64(a * a) + Float64(b * b)) / Float64(1.0 / Float64(b * b)))); end return tmp end
function tmp_2 = code(a, b) tmp = 0.0; if ((b * b) <= 0.0005) tmp = -1.0 + (a * (a * ((a * a) + (4.0 + (a * 4.0))))); else tmp = (((b * b) * 4.0) + -1.0) + (((a * a) + (b * b)) / (1.0 / (b * b))); end tmp_2 = tmp; end
code[a_, b_] := If[LessEqual[N[(b * b), $MachinePrecision], 0.0005], N[(-1.0 + N[(a * N[(a * N[(N[(a * a), $MachinePrecision] + N[(4.0 + N[(a * 4.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(N[(b * b), $MachinePrecision] * 4.0), $MachinePrecision] + -1.0), $MachinePrecision] + N[(N[(N[(a * a), $MachinePrecision] + N[(b * b), $MachinePrecision]), $MachinePrecision] / N[(1.0 / N[(b * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \cdot b \leq 0.0005:\\
\;\;\;\;-1 + a \cdot \left(a \cdot \left(a \cdot a + \left(4 + a \cdot 4\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\left(\left(b \cdot b\right) \cdot 4 + -1\right) + \frac{a \cdot a + b \cdot b}{\frac{1}{b \cdot b}}\\
\end{array}
\end{array}
if (*.f64 b b) < 5.0000000000000001e-4Initial program 86.3%
associate--l+N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
Simplified86.3%
Taylor expanded in b around 0
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
+-commutativeN/A
metadata-evalN/A
pow-sqrN/A
distribute-lft-outN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f6499.4%
Simplified99.4%
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f6499.5%
Applied egg-rr99.5%
if 5.0000000000000001e-4 < (*.f64 b b) Initial program 62.3%
associate--l+N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
Simplified62.3%
Taylor expanded in a around 0
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6499.9%
Simplified99.9%
+-commutativeN/A
flip-+N/A
associate-*r*N/A
associate-*r*N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
clear-numN/A
metadata-evalN/A
Applied egg-rr99.9%
Taylor expanded in a around 0
unpow2N/A
*-lowering-*.f6496.0%
Simplified96.0%
Final simplification97.6%
(FPCore (a b) :precision binary64 (let* ((t_0 (+ (* a a) (* b b)))) (+ (+ (* (* b b) 4.0) -1.0) (* t_0 t_0))))
double code(double a, double b) {
double t_0 = (a * a) + (b * b);
return (((b * b) * 4.0) + -1.0) + (t_0 * t_0);
}
real(8) function code(a, b)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: t_0
t_0 = (a * a) + (b * b)
code = (((b * b) * 4.0d0) + (-1.0d0)) + (t_0 * t_0)
end function
public static double code(double a, double b) {
double t_0 = (a * a) + (b * b);
return (((b * b) * 4.0) + -1.0) + (t_0 * t_0);
}
def code(a, b): t_0 = (a * a) + (b * b) return (((b * b) * 4.0) + -1.0) + (t_0 * t_0)
function code(a, b) t_0 = Float64(Float64(a * a) + Float64(b * b)) return Float64(Float64(Float64(Float64(b * b) * 4.0) + -1.0) + Float64(t_0 * t_0)) end
function tmp = code(a, b) t_0 = (a * a) + (b * b); tmp = (((b * b) * 4.0) + -1.0) + (t_0 * t_0); end
code[a_, b_] := Block[{t$95$0 = N[(N[(a * a), $MachinePrecision] + N[(b * b), $MachinePrecision]), $MachinePrecision]}, N[(N[(N[(N[(b * b), $MachinePrecision] * 4.0), $MachinePrecision] + -1.0), $MachinePrecision] + N[(t$95$0 * t$95$0), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := a \cdot a + b \cdot b\\
\left(\left(b \cdot b\right) \cdot 4 + -1\right) + t\_0 \cdot t\_0
\end{array}
\end{array}
Initial program 73.4%
associate--l+N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
Simplified73.4%
Taylor expanded in a around 0
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6498.8%
Simplified98.8%
Final simplification98.8%
(FPCore (a b) :precision binary64 (if (<= (* b b) 2e+147) (+ -1.0 (* a (* a (+ (* a a) (+ 4.0 (* a 4.0)))))) (* b (* b (* b b)))))
double code(double a, double b) {
double tmp;
if ((b * b) <= 2e+147) {
tmp = -1.0 + (a * (a * ((a * a) + (4.0 + (a * 4.0)))));
} else {
tmp = b * (b * (b * b));
}
return tmp;
}
real(8) function code(a, b)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: tmp
if ((b * b) <= 2d+147) then
tmp = (-1.0d0) + (a * (a * ((a * a) + (4.0d0 + (a * 4.0d0)))))
else
tmp = b * (b * (b * b))
end if
code = tmp
end function
public static double code(double a, double b) {
double tmp;
if ((b * b) <= 2e+147) {
tmp = -1.0 + (a * (a * ((a * a) + (4.0 + (a * 4.0)))));
} else {
tmp = b * (b * (b * b));
}
return tmp;
}
def code(a, b): tmp = 0 if (b * b) <= 2e+147: tmp = -1.0 + (a * (a * ((a * a) + (4.0 + (a * 4.0))))) else: tmp = b * (b * (b * b)) return tmp
function code(a, b) tmp = 0.0 if (Float64(b * b) <= 2e+147) tmp = Float64(-1.0 + Float64(a * Float64(a * Float64(Float64(a * a) + Float64(4.0 + Float64(a * 4.0)))))); else tmp = Float64(b * Float64(b * Float64(b * b))); end return tmp end
function tmp_2 = code(a, b) tmp = 0.0; if ((b * b) <= 2e+147) tmp = -1.0 + (a * (a * ((a * a) + (4.0 + (a * 4.0))))); else tmp = b * (b * (b * b)); end tmp_2 = tmp; end
code[a_, b_] := If[LessEqual[N[(b * b), $MachinePrecision], 2e+147], N[(-1.0 + N[(a * N[(a * N[(N[(a * a), $MachinePrecision] + N[(4.0 + N[(a * 4.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(b * N[(b * N[(b * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \cdot b \leq 2 \cdot 10^{+147}:\\
\;\;\;\;-1 + a \cdot \left(a \cdot \left(a \cdot a + \left(4 + a \cdot 4\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;b \cdot \left(b \cdot \left(b \cdot b\right)\right)\\
\end{array}
\end{array}
if (*.f64 b b) < 2e147Initial program 80.8%
associate--l+N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
Simplified80.8%
Taylor expanded in b around 0
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
+-commutativeN/A
metadata-evalN/A
pow-sqrN/A
distribute-lft-outN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f6492.1%
Simplified92.1%
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f6492.2%
Applied egg-rr92.2%
if 2e147 < (*.f64 b b) Initial program 63.3%
associate--l+N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
Simplified63.3%
Taylor expanded in b around inf
metadata-evalN/A
pow-sqrN/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64100.0%
Simplified100.0%
Final simplification95.5%
(FPCore (a b) :precision binary64 (if (<= (* b b) 2e-9) (+ -1.0 (* (* a a) 4.0)) (if (<= (* b b) 2e+147) (* a (* a (* a a))) (* b (* b (* b b))))))
double code(double a, double b) {
double tmp;
if ((b * b) <= 2e-9) {
tmp = -1.0 + ((a * a) * 4.0);
} else if ((b * b) <= 2e+147) {
tmp = a * (a * (a * a));
} else {
tmp = b * (b * (b * b));
}
return tmp;
}
real(8) function code(a, b)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: tmp
if ((b * b) <= 2d-9) then
tmp = (-1.0d0) + ((a * a) * 4.0d0)
else if ((b * b) <= 2d+147) then
tmp = a * (a * (a * a))
else
tmp = b * (b * (b * b))
end if
code = tmp
end function
public static double code(double a, double b) {
double tmp;
if ((b * b) <= 2e-9) {
tmp = -1.0 + ((a * a) * 4.0);
} else if ((b * b) <= 2e+147) {
tmp = a * (a * (a * a));
} else {
tmp = b * (b * (b * b));
}
return tmp;
}
def code(a, b): tmp = 0 if (b * b) <= 2e-9: tmp = -1.0 + ((a * a) * 4.0) elif (b * b) <= 2e+147: tmp = a * (a * (a * a)) else: tmp = b * (b * (b * b)) return tmp
function code(a, b) tmp = 0.0 if (Float64(b * b) <= 2e-9) tmp = Float64(-1.0 + Float64(Float64(a * a) * 4.0)); elseif (Float64(b * b) <= 2e+147) tmp = Float64(a * Float64(a * Float64(a * a))); else tmp = Float64(b * Float64(b * Float64(b * b))); end return tmp end
function tmp_2 = code(a, b) tmp = 0.0; if ((b * b) <= 2e-9) tmp = -1.0 + ((a * a) * 4.0); elseif ((b * b) <= 2e+147) tmp = a * (a * (a * a)); else tmp = b * (b * (b * b)); end tmp_2 = tmp; end
code[a_, b_] := If[LessEqual[N[(b * b), $MachinePrecision], 2e-9], N[(-1.0 + N[(N[(a * a), $MachinePrecision] * 4.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[N[(b * b), $MachinePrecision], 2e+147], N[(a * N[(a * N[(a * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(b * N[(b * N[(b * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \cdot b \leq 2 \cdot 10^{-9}:\\
\;\;\;\;-1 + \left(a \cdot a\right) \cdot 4\\
\mathbf{elif}\;b \cdot b \leq 2 \cdot 10^{+147}:\\
\;\;\;\;a \cdot \left(a \cdot \left(a \cdot a\right)\right)\\
\mathbf{else}:\\
\;\;\;\;b \cdot \left(b \cdot \left(b \cdot b\right)\right)\\
\end{array}
\end{array}
if (*.f64 b b) < 2.00000000000000012e-9Initial program 86.2%
associate--l+N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
Simplified86.2%
Taylor expanded in b around 0
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
+-commutativeN/A
metadata-evalN/A
pow-sqrN/A
distribute-lft-outN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f6499.5%
Simplified99.5%
Taylor expanded in a around 0
Simplified75.3%
if 2.00000000000000012e-9 < (*.f64 b b) < 2e147Initial program 59.6%
associate--l+N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
Simplified59.6%
Taylor expanded in a around 0
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6499.6%
Simplified99.6%
Taylor expanded in a around inf
metadata-evalN/A
pow-sqrN/A
unpow2N/A
associate-*l*N/A
unpow2N/A
cube-multN/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6464.2%
Simplified64.2%
if 2e147 < (*.f64 b b) Initial program 63.3%
associate--l+N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
Simplified63.3%
Taylor expanded in b around inf
metadata-evalN/A
pow-sqrN/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64100.0%
Simplified100.0%
(FPCore (a b) :precision binary64 (if (<= (* b b) 2e+147) (+ -1.0 (* a (* a (+ 4.0 (* a (+ a 4.0)))))) (* b (* b (* b b)))))
double code(double a, double b) {
double tmp;
if ((b * b) <= 2e+147) {
tmp = -1.0 + (a * (a * (4.0 + (a * (a + 4.0)))));
} else {
tmp = b * (b * (b * b));
}
return tmp;
}
real(8) function code(a, b)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: tmp
if ((b * b) <= 2d+147) then
tmp = (-1.0d0) + (a * (a * (4.0d0 + (a * (a + 4.0d0)))))
else
tmp = b * (b * (b * b))
end if
code = tmp
end function
public static double code(double a, double b) {
double tmp;
if ((b * b) <= 2e+147) {
tmp = -1.0 + (a * (a * (4.0 + (a * (a + 4.0)))));
} else {
tmp = b * (b * (b * b));
}
return tmp;
}
def code(a, b): tmp = 0 if (b * b) <= 2e+147: tmp = -1.0 + (a * (a * (4.0 + (a * (a + 4.0))))) else: tmp = b * (b * (b * b)) return tmp
function code(a, b) tmp = 0.0 if (Float64(b * b) <= 2e+147) tmp = Float64(-1.0 + Float64(a * Float64(a * Float64(4.0 + Float64(a * Float64(a + 4.0)))))); else tmp = Float64(b * Float64(b * Float64(b * b))); end return tmp end
function tmp_2 = code(a, b) tmp = 0.0; if ((b * b) <= 2e+147) tmp = -1.0 + (a * (a * (4.0 + (a * (a + 4.0))))); else tmp = b * (b * (b * b)); end tmp_2 = tmp; end
code[a_, b_] := If[LessEqual[N[(b * b), $MachinePrecision], 2e+147], N[(-1.0 + N[(a * N[(a * N[(4.0 + N[(a * N[(a + 4.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(b * N[(b * N[(b * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \cdot b \leq 2 \cdot 10^{+147}:\\
\;\;\;\;-1 + a \cdot \left(a \cdot \left(4 + a \cdot \left(a + 4\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;b \cdot \left(b \cdot \left(b \cdot b\right)\right)\\
\end{array}
\end{array}
if (*.f64 b b) < 2e147Initial program 80.8%
associate--l+N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
Simplified80.8%
Taylor expanded in b around 0
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
+-commutativeN/A
metadata-evalN/A
pow-sqrN/A
distribute-lft-outN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f6492.1%
Simplified92.1%
Taylor expanded in a around 0
sub-negN/A
distribute-lft-inN/A
*-commutativeN/A
distribute-rgt-inN/A
unpow2N/A
distribute-rgt-inN/A
pow-sqrN/A
metadata-evalN/A
associate-+l+N/A
distribute-rgt-inN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
distribute-rgt-inN/A
Simplified92.1%
+-commutativeN/A
+-lowering-+.f64N/A
associate-*l*N/A
distribute-lft-inN/A
associate-+r+N/A
+-commutativeN/A
*-lowering-*.f64N/A
+-commutativeN/A
associate-+r+N/A
distribute-lft-inN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f6492.2%
Applied egg-rr92.2%
if 2e147 < (*.f64 b b) Initial program 63.3%
associate--l+N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
Simplified63.3%
Taylor expanded in b around inf
metadata-evalN/A
pow-sqrN/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64100.0%
Simplified100.0%
Final simplification95.5%
(FPCore (a b) :precision binary64 (if (<= (* b b) 2e+147) (+ -1.0 (* (* a a) (+ 4.0 (* a (+ a 4.0))))) (* b (* b (* b b)))))
double code(double a, double b) {
double tmp;
if ((b * b) <= 2e+147) {
tmp = -1.0 + ((a * a) * (4.0 + (a * (a + 4.0))));
} else {
tmp = b * (b * (b * b));
}
return tmp;
}
real(8) function code(a, b)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: tmp
if ((b * b) <= 2d+147) then
tmp = (-1.0d0) + ((a * a) * (4.0d0 + (a * (a + 4.0d0))))
else
tmp = b * (b * (b * b))
end if
code = tmp
end function
public static double code(double a, double b) {
double tmp;
if ((b * b) <= 2e+147) {
tmp = -1.0 + ((a * a) * (4.0 + (a * (a + 4.0))));
} else {
tmp = b * (b * (b * b));
}
return tmp;
}
def code(a, b): tmp = 0 if (b * b) <= 2e+147: tmp = -1.0 + ((a * a) * (4.0 + (a * (a + 4.0)))) else: tmp = b * (b * (b * b)) return tmp
function code(a, b) tmp = 0.0 if (Float64(b * b) <= 2e+147) tmp = Float64(-1.0 + Float64(Float64(a * a) * Float64(4.0 + Float64(a * Float64(a + 4.0))))); else tmp = Float64(b * Float64(b * Float64(b * b))); end return tmp end
function tmp_2 = code(a, b) tmp = 0.0; if ((b * b) <= 2e+147) tmp = -1.0 + ((a * a) * (4.0 + (a * (a + 4.0)))); else tmp = b * (b * (b * b)); end tmp_2 = tmp; end
code[a_, b_] := If[LessEqual[N[(b * b), $MachinePrecision], 2e+147], N[(-1.0 + N[(N[(a * a), $MachinePrecision] * N[(4.0 + N[(a * N[(a + 4.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(b * N[(b * N[(b * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \cdot b \leq 2 \cdot 10^{+147}:\\
\;\;\;\;-1 + \left(a \cdot a\right) \cdot \left(4 + a \cdot \left(a + 4\right)\right)\\
\mathbf{else}:\\
\;\;\;\;b \cdot \left(b \cdot \left(b \cdot b\right)\right)\\
\end{array}
\end{array}
if (*.f64 b b) < 2e147Initial program 80.8%
associate--l+N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
Simplified80.8%
Taylor expanded in b around 0
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
+-commutativeN/A
metadata-evalN/A
pow-sqrN/A
distribute-lft-outN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f6492.1%
Simplified92.1%
Taylor expanded in a around 0
sub-negN/A
distribute-lft-inN/A
*-commutativeN/A
distribute-rgt-inN/A
unpow2N/A
distribute-rgt-inN/A
pow-sqrN/A
metadata-evalN/A
associate-+l+N/A
distribute-rgt-inN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
distribute-rgt-inN/A
Simplified92.1%
if 2e147 < (*.f64 b b) Initial program 63.3%
associate--l+N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
Simplified63.3%
Taylor expanded in b around inf
metadata-evalN/A
pow-sqrN/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64100.0%
Simplified100.0%
Final simplification95.5%
(FPCore (a b) :precision binary64 (if (<= (* b b) 2e+147) (+ -1.0 (* a (* a (* a (+ a 4.0))))) (* b (* b (* b b)))))
double code(double a, double b) {
double tmp;
if ((b * b) <= 2e+147) {
tmp = -1.0 + (a * (a * (a * (a + 4.0))));
} else {
tmp = b * (b * (b * b));
}
return tmp;
}
real(8) function code(a, b)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: tmp
if ((b * b) <= 2d+147) then
tmp = (-1.0d0) + (a * (a * (a * (a + 4.0d0))))
else
tmp = b * (b * (b * b))
end if
code = tmp
end function
public static double code(double a, double b) {
double tmp;
if ((b * b) <= 2e+147) {
tmp = -1.0 + (a * (a * (a * (a + 4.0))));
} else {
tmp = b * (b * (b * b));
}
return tmp;
}
def code(a, b): tmp = 0 if (b * b) <= 2e+147: tmp = -1.0 + (a * (a * (a * (a + 4.0)))) else: tmp = b * (b * (b * b)) return tmp
function code(a, b) tmp = 0.0 if (Float64(b * b) <= 2e+147) tmp = Float64(-1.0 + Float64(a * Float64(a * Float64(a * Float64(a + 4.0))))); else tmp = Float64(b * Float64(b * Float64(b * b))); end return tmp end
function tmp_2 = code(a, b) tmp = 0.0; if ((b * b) <= 2e+147) tmp = -1.0 + (a * (a * (a * (a + 4.0)))); else tmp = b * (b * (b * b)); end tmp_2 = tmp; end
code[a_, b_] := If[LessEqual[N[(b * b), $MachinePrecision], 2e+147], N[(-1.0 + N[(a * N[(a * N[(a * N[(a + 4.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(b * N[(b * N[(b * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \cdot b \leq 2 \cdot 10^{+147}:\\
\;\;\;\;-1 + a \cdot \left(a \cdot \left(a \cdot \left(a + 4\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;b \cdot \left(b \cdot \left(b \cdot b\right)\right)\\
\end{array}
\end{array}
if (*.f64 b b) < 2e147Initial program 80.8%
associate--l+N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
Simplified80.8%
Taylor expanded in b around 0
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
+-commutativeN/A
metadata-evalN/A
pow-sqrN/A
distribute-lft-outN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f6492.1%
Simplified92.1%
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f6492.2%
Applied egg-rr92.2%
Taylor expanded in a around inf
+-commutativeN/A
distribute-rgt-inN/A
cube-multN/A
unpow2N/A
associate-*r*N/A
*-lft-identityN/A
fma-defineN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
fma-defineN/A
cube-multN/A
unpow2N/A
distribute-rgt-inN/A
unpow2N/A
associate-*r*N/A
*-lowering-*.f64N/A
distribute-lft-inN/A
unpow2N/A
fma-defineN/A
metadata-evalN/A
lft-mult-inverseN/A
associate-*l*N/A
fma-undefineN/A
unpow2N/A
Simplified90.9%
if 2e147 < (*.f64 b b) Initial program 63.3%
associate--l+N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
Simplified63.3%
Taylor expanded in b around inf
metadata-evalN/A
pow-sqrN/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64100.0%
Simplified100.0%
Final simplification94.8%
(FPCore (a b) :precision binary64 (if (<= (* b b) 2e+147) (+ -1.0 (* (* a a) (* a (+ a 4.0)))) (* b (* b (* b b)))))
double code(double a, double b) {
double tmp;
if ((b * b) <= 2e+147) {
tmp = -1.0 + ((a * a) * (a * (a + 4.0)));
} else {
tmp = b * (b * (b * b));
}
return tmp;
}
real(8) function code(a, b)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: tmp
if ((b * b) <= 2d+147) then
tmp = (-1.0d0) + ((a * a) * (a * (a + 4.0d0)))
else
tmp = b * (b * (b * b))
end if
code = tmp
end function
public static double code(double a, double b) {
double tmp;
if ((b * b) <= 2e+147) {
tmp = -1.0 + ((a * a) * (a * (a + 4.0)));
} else {
tmp = b * (b * (b * b));
}
return tmp;
}
def code(a, b): tmp = 0 if (b * b) <= 2e+147: tmp = -1.0 + ((a * a) * (a * (a + 4.0))) else: tmp = b * (b * (b * b)) return tmp
function code(a, b) tmp = 0.0 if (Float64(b * b) <= 2e+147) tmp = Float64(-1.0 + Float64(Float64(a * a) * Float64(a * Float64(a + 4.0)))); else tmp = Float64(b * Float64(b * Float64(b * b))); end return tmp end
function tmp_2 = code(a, b) tmp = 0.0; if ((b * b) <= 2e+147) tmp = -1.0 + ((a * a) * (a * (a + 4.0))); else tmp = b * (b * (b * b)); end tmp_2 = tmp; end
code[a_, b_] := If[LessEqual[N[(b * b), $MachinePrecision], 2e+147], N[(-1.0 + N[(N[(a * a), $MachinePrecision] * N[(a * N[(a + 4.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(b * N[(b * N[(b * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \cdot b \leq 2 \cdot 10^{+147}:\\
\;\;\;\;-1 + \left(a \cdot a\right) \cdot \left(a \cdot \left(a + 4\right)\right)\\
\mathbf{else}:\\
\;\;\;\;b \cdot \left(b \cdot \left(b \cdot b\right)\right)\\
\end{array}
\end{array}
if (*.f64 b b) < 2e147Initial program 80.8%
associate--l+N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
Simplified80.8%
Taylor expanded in b around 0
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
+-commutativeN/A
metadata-evalN/A
pow-sqrN/A
distribute-lft-outN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f6492.1%
Simplified92.1%
Taylor expanded in a around 0
sub-negN/A
distribute-lft-inN/A
*-commutativeN/A
distribute-rgt-inN/A
unpow2N/A
distribute-rgt-inN/A
pow-sqrN/A
metadata-evalN/A
associate-+l+N/A
distribute-rgt-inN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
distribute-rgt-inN/A
Simplified92.1%
Taylor expanded in a around inf
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-commutativeN/A
distribute-lft1-inN/A
associate-*l*N/A
lft-mult-inverseN/A
metadata-evalN/A
+-lowering-+.f6490.9%
Simplified90.9%
if 2e147 < (*.f64 b b) Initial program 63.3%
associate--l+N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
Simplified63.3%
Taylor expanded in b around inf
metadata-evalN/A
pow-sqrN/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64100.0%
Simplified100.0%
Final simplification94.8%
(FPCore (a b) :precision binary64 (let* ((t_0 (* a (* a (* a a))))) (if (<= a -2.4) t_0 (if (<= a 1.25e-6) -1.0 t_0))))
double code(double a, double b) {
double t_0 = a * (a * (a * a));
double tmp;
if (a <= -2.4) {
tmp = t_0;
} else if (a <= 1.25e-6) {
tmp = -1.0;
} else {
tmp = t_0;
}
return tmp;
}
real(8) function code(a, b)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: t_0
real(8) :: tmp
t_0 = a * (a * (a * a))
if (a <= (-2.4d0)) then
tmp = t_0
else if (a <= 1.25d-6) then
tmp = -1.0d0
else
tmp = t_0
end if
code = tmp
end function
public static double code(double a, double b) {
double t_0 = a * (a * (a * a));
double tmp;
if (a <= -2.4) {
tmp = t_0;
} else if (a <= 1.25e-6) {
tmp = -1.0;
} else {
tmp = t_0;
}
return tmp;
}
def code(a, b): t_0 = a * (a * (a * a)) tmp = 0 if a <= -2.4: tmp = t_0 elif a <= 1.25e-6: tmp = -1.0 else: tmp = t_0 return tmp
function code(a, b) t_0 = Float64(a * Float64(a * Float64(a * a))) tmp = 0.0 if (a <= -2.4) tmp = t_0; elseif (a <= 1.25e-6) tmp = -1.0; else tmp = t_0; end return tmp end
function tmp_2 = code(a, b) t_0 = a * (a * (a * a)); tmp = 0.0; if (a <= -2.4) tmp = t_0; elseif (a <= 1.25e-6) tmp = -1.0; else tmp = t_0; end tmp_2 = tmp; end
code[a_, b_] := Block[{t$95$0 = N[(a * N[(a * N[(a * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[a, -2.4], t$95$0, If[LessEqual[a, 1.25e-6], -1.0, t$95$0]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := a \cdot \left(a \cdot \left(a \cdot a\right)\right)\\
\mathbf{if}\;a \leq -2.4:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;a \leq 1.25 \cdot 10^{-6}:\\
\;\;\;\;-1\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if a < -2.39999999999999991 or 1.2500000000000001e-6 < a Initial program 45.5%
associate--l+N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
Simplified45.5%
Taylor expanded in a around 0
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6498.1%
Simplified98.1%
Taylor expanded in a around inf
metadata-evalN/A
pow-sqrN/A
unpow2N/A
associate-*l*N/A
unpow2N/A
cube-multN/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6488.2%
Simplified88.2%
if -2.39999999999999991 < a < 1.2500000000000001e-6Initial program 100.0%
associate--l+N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
Simplified100.0%
Taylor expanded in a around 0
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
+-commutativeN/A
metadata-evalN/A
pow-sqrN/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f6499.6%
Simplified99.6%
Taylor expanded in b around 0
Simplified46.9%
(FPCore (a b) :precision binary64 (let* ((t_0 (* (* a a) 4.0))) (if (<= a -2.4) t_0 (if (<= a 1.25e-6) -1.0 t_0))))
double code(double a, double b) {
double t_0 = (a * a) * 4.0;
double tmp;
if (a <= -2.4) {
tmp = t_0;
} else if (a <= 1.25e-6) {
tmp = -1.0;
} else {
tmp = t_0;
}
return tmp;
}
real(8) function code(a, b)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: t_0
real(8) :: tmp
t_0 = (a * a) * 4.0d0
if (a <= (-2.4d0)) then
tmp = t_0
else if (a <= 1.25d-6) then
tmp = -1.0d0
else
tmp = t_0
end if
code = tmp
end function
public static double code(double a, double b) {
double t_0 = (a * a) * 4.0;
double tmp;
if (a <= -2.4) {
tmp = t_0;
} else if (a <= 1.25e-6) {
tmp = -1.0;
} else {
tmp = t_0;
}
return tmp;
}
def code(a, b): t_0 = (a * a) * 4.0 tmp = 0 if a <= -2.4: tmp = t_0 elif a <= 1.25e-6: tmp = -1.0 else: tmp = t_0 return tmp
function code(a, b) t_0 = Float64(Float64(a * a) * 4.0) tmp = 0.0 if (a <= -2.4) tmp = t_0; elseif (a <= 1.25e-6) tmp = -1.0; else tmp = t_0; end return tmp end
function tmp_2 = code(a, b) t_0 = (a * a) * 4.0; tmp = 0.0; if (a <= -2.4) tmp = t_0; elseif (a <= 1.25e-6) tmp = -1.0; else tmp = t_0; end tmp_2 = tmp; end
code[a_, b_] := Block[{t$95$0 = N[(N[(a * a), $MachinePrecision] * 4.0), $MachinePrecision]}, If[LessEqual[a, -2.4], t$95$0, If[LessEqual[a, 1.25e-6], -1.0, t$95$0]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(a \cdot a\right) \cdot 4\\
\mathbf{if}\;a \leq -2.4:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;a \leq 1.25 \cdot 10^{-6}:\\
\;\;\;\;-1\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if a < -2.39999999999999991 or 1.2500000000000001e-6 < a Initial program 45.5%
associate--l+N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
Simplified45.5%
Taylor expanded in b around 0
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
+-commutativeN/A
metadata-evalN/A
pow-sqrN/A
distribute-lft-outN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f6490.0%
Simplified90.0%
Taylor expanded in a around 0
Simplified51.1%
Taylor expanded in a around inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6451.1%
Simplified51.1%
if -2.39999999999999991 < a < 1.2500000000000001e-6Initial program 100.0%
associate--l+N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
Simplified100.0%
Taylor expanded in a around 0
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
+-commutativeN/A
metadata-evalN/A
pow-sqrN/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f6499.6%
Simplified99.6%
Taylor expanded in b around 0
Simplified46.9%
Final simplification49.0%
(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 73.4%
associate--l+N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
Simplified73.4%
Taylor expanded in a around 0
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
+-commutativeN/A
metadata-evalN/A
pow-sqrN/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f6471.4%
Simplified71.4%
Taylor expanded in b around 0
Simplified24.3%
herbie shell --seed 2024144
(FPCore (a b)
:name "Bouland and Aaronson, Equation (25)"
:precision binary64
(- (+ (pow (+ (* a a) (* b b)) 2.0) (* 4.0 (+ (* (* a a) (+ 1.0 a)) (* (* b b) (- 1.0 (* 3.0 a)))))) 1.0))