
(FPCore (x y) :precision binary64 (+ (- 1.0 x) (* y (sqrt x))))
double code(double x, double y) {
return (1.0 - x) + (y * sqrt(x));
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = (1.0d0 - x) + (y * sqrt(x))
end function
public static double code(double x, double y) {
return (1.0 - x) + (y * Math.sqrt(x));
}
def code(x, y): return (1.0 - x) + (y * math.sqrt(x))
function code(x, y) return Float64(Float64(1.0 - x) + Float64(y * sqrt(x))) end
function tmp = code(x, y) tmp = (1.0 - x) + (y * sqrt(x)); end
code[x_, y_] := N[(N[(1.0 - x), $MachinePrecision] + N[(y * N[Sqrt[x], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(1 - x\right) + y \cdot \sqrt{x}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 10 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x y) :precision binary64 (+ (- 1.0 x) (* y (sqrt x))))
double code(double x, double y) {
return (1.0 - x) + (y * sqrt(x));
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = (1.0d0 - x) + (y * sqrt(x))
end function
public static double code(double x, double y) {
return (1.0 - x) + (y * Math.sqrt(x));
}
def code(x, y): return (1.0 - x) + (y * math.sqrt(x))
function code(x, y) return Float64(Float64(1.0 - x) + Float64(y * sqrt(x))) end
function tmp = code(x, y) tmp = (1.0 - x) + (y * sqrt(x)); end
code[x_, y_] := N[(N[(1.0 - x), $MachinePrecision] + N[(y * N[Sqrt[x], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(1 - x\right) + y \cdot \sqrt{x}
\end{array}
(FPCore (x y) :precision binary64 (+ (- 1.0 x) (* y (sqrt x))))
double code(double x, double y) {
return (1.0 - x) + (y * sqrt(x));
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = (1.0d0 - x) + (y * sqrt(x))
end function
public static double code(double x, double y) {
return (1.0 - x) + (y * Math.sqrt(x));
}
def code(x, y): return (1.0 - x) + (y * math.sqrt(x))
function code(x, y) return Float64(Float64(1.0 - x) + Float64(y * sqrt(x))) end
function tmp = code(x, y) tmp = (1.0 - x) + (y * sqrt(x)); end
code[x_, y_] := N[(N[(1.0 - x), $MachinePrecision] + N[(y * N[Sqrt[x], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(1 - x\right) + y \cdot \sqrt{x}
\end{array}
Initial program 99.9%
(FPCore (x y) :precision binary64 (let* ((t_0 (* y (sqrt x)))) (if (<= y -1.8e+53) t_0 (if (<= y 3100000000000.0) (- 1.0 x) (+ 1.0 t_0)))))
double code(double x, double y) {
double t_0 = y * sqrt(x);
double tmp;
if (y <= -1.8e+53) {
tmp = t_0;
} else if (y <= 3100000000000.0) {
tmp = 1.0 - x;
} else {
tmp = 1.0 + t_0;
}
return tmp;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8) :: t_0
real(8) :: tmp
t_0 = y * sqrt(x)
if (y <= (-1.8d+53)) then
tmp = t_0
else if (y <= 3100000000000.0d0) then
tmp = 1.0d0 - x
else
tmp = 1.0d0 + t_0
end if
code = tmp
end function
public static double code(double x, double y) {
double t_0 = y * Math.sqrt(x);
double tmp;
if (y <= -1.8e+53) {
tmp = t_0;
} else if (y <= 3100000000000.0) {
tmp = 1.0 - x;
} else {
tmp = 1.0 + t_0;
}
return tmp;
}
def code(x, y): t_0 = y * math.sqrt(x) tmp = 0 if y <= -1.8e+53: tmp = t_0 elif y <= 3100000000000.0: tmp = 1.0 - x else: tmp = 1.0 + t_0 return tmp
function code(x, y) t_0 = Float64(y * sqrt(x)) tmp = 0.0 if (y <= -1.8e+53) tmp = t_0; elseif (y <= 3100000000000.0) tmp = Float64(1.0 - x); else tmp = Float64(1.0 + t_0); end return tmp end
function tmp_2 = code(x, y) t_0 = y * sqrt(x); tmp = 0.0; if (y <= -1.8e+53) tmp = t_0; elseif (y <= 3100000000000.0) tmp = 1.0 - x; else tmp = 1.0 + t_0; end tmp_2 = tmp; end
code[x_, y_] := Block[{t$95$0 = N[(y * N[Sqrt[x], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -1.8e+53], t$95$0, If[LessEqual[y, 3100000000000.0], N[(1.0 - x), $MachinePrecision], N[(1.0 + t$95$0), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := y \cdot \sqrt{x}\\
\mathbf{if}\;y \leq -1.8 \cdot 10^{+53}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;y \leq 3100000000000:\\
\;\;\;\;1 - x\\
\mathbf{else}:\\
\;\;\;\;1 + t\_0\\
\end{array}
\end{array}
if y < -1.8e53Initial program 99.5%
+-commutative99.5%
*-commutative99.5%
add-sqr-sqrt99.3%
associate-*l*99.4%
fma-define99.4%
pow1/299.4%
sqrt-pow199.4%
metadata-eval99.4%
pow1/299.4%
sqrt-pow199.4%
metadata-eval99.4%
Applied egg-rr99.4%
Taylor expanded in y around inf 91.6%
if -1.8e53 < y < 3.1e12Initial program 100.0%
Taylor expanded in y around 0 98.9%
if 3.1e12 < y Initial program 99.8%
Taylor expanded in x around 0 89.4%
Final simplification95.2%
(FPCore (x y) :precision binary64 (if (or (<= y -9.2e+55) (not (<= y 2.45e+50))) (* y (sqrt x)) (- 1.0 x)))
double code(double x, double y) {
double tmp;
if ((y <= -9.2e+55) || !(y <= 2.45e+50)) {
tmp = y * sqrt(x);
} else {
tmp = 1.0 - x;
}
return tmp;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8) :: tmp
if ((y <= (-9.2d+55)) .or. (.not. (y <= 2.45d+50))) then
tmp = y * sqrt(x)
else
tmp = 1.0d0 - x
end if
code = tmp
end function
public static double code(double x, double y) {
double tmp;
if ((y <= -9.2e+55) || !(y <= 2.45e+50)) {
tmp = y * Math.sqrt(x);
} else {
tmp = 1.0 - x;
}
return tmp;
}
def code(x, y): tmp = 0 if (y <= -9.2e+55) or not (y <= 2.45e+50): tmp = y * math.sqrt(x) else: tmp = 1.0 - x return tmp
function code(x, y) tmp = 0.0 if ((y <= -9.2e+55) || !(y <= 2.45e+50)) tmp = Float64(y * sqrt(x)); else tmp = Float64(1.0 - x); end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if ((y <= -9.2e+55) || ~((y <= 2.45e+50))) tmp = y * sqrt(x); else tmp = 1.0 - x; end tmp_2 = tmp; end
code[x_, y_] := If[Or[LessEqual[y, -9.2e+55], N[Not[LessEqual[y, 2.45e+50]], $MachinePrecision]], N[(y * N[Sqrt[x], $MachinePrecision]), $MachinePrecision], N[(1.0 - x), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -9.2 \cdot 10^{+55} \lor \neg \left(y \leq 2.45 \cdot 10^{+50}\right):\\
\;\;\;\;y \cdot \sqrt{x}\\
\mathbf{else}:\\
\;\;\;\;1 - x\\
\end{array}
\end{array}
if y < -9.1999999999999995e55 or 2.4500000000000001e50 < y Initial program 99.7%
+-commutative99.7%
*-commutative99.7%
add-sqr-sqrt99.3%
associate-*l*99.4%
fma-define99.4%
pow1/299.4%
sqrt-pow199.6%
metadata-eval99.6%
pow1/299.6%
sqrt-pow199.4%
metadata-eval99.4%
Applied egg-rr99.4%
Taylor expanded in y around inf 91.5%
if -9.1999999999999995e55 < y < 2.4500000000000001e50Initial program 100.0%
Taylor expanded in y around 0 96.8%
Final simplification94.6%
(FPCore (x y) :precision binary64 (let* ((t_0 (* y (sqrt x)))) (if (<= x 0.09) (+ 1.0 t_0) (- t_0 x))))
double code(double x, double y) {
double t_0 = y * sqrt(x);
double tmp;
if (x <= 0.09) {
tmp = 1.0 + t_0;
} else {
tmp = t_0 - x;
}
return tmp;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8) :: t_0
real(8) :: tmp
t_0 = y * sqrt(x)
if (x <= 0.09d0) then
tmp = 1.0d0 + t_0
else
tmp = t_0 - x
end if
code = tmp
end function
public static double code(double x, double y) {
double t_0 = y * Math.sqrt(x);
double tmp;
if (x <= 0.09) {
tmp = 1.0 + t_0;
} else {
tmp = t_0 - x;
}
return tmp;
}
def code(x, y): t_0 = y * math.sqrt(x) tmp = 0 if x <= 0.09: tmp = 1.0 + t_0 else: tmp = t_0 - x return tmp
function code(x, y) t_0 = Float64(y * sqrt(x)) tmp = 0.0 if (x <= 0.09) tmp = Float64(1.0 + t_0); else tmp = Float64(t_0 - x); end return tmp end
function tmp_2 = code(x, y) t_0 = y * sqrt(x); tmp = 0.0; if (x <= 0.09) tmp = 1.0 + t_0; else tmp = t_0 - x; end tmp_2 = tmp; end
code[x_, y_] := Block[{t$95$0 = N[(y * N[Sqrt[x], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[x, 0.09], N[(1.0 + t$95$0), $MachinePrecision], N[(t$95$0 - x), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := y \cdot \sqrt{x}\\
\mathbf{if}\;x \leq 0.09:\\
\;\;\;\;1 + t\_0\\
\mathbf{else}:\\
\;\;\;\;t\_0 - x\\
\end{array}
\end{array}
if x < 0.089999999999999997Initial program 99.8%
Taylor expanded in x around 0 98.0%
if 0.089999999999999997 < x Initial program 99.9%
Taylor expanded in x around inf 99.7%
neg-mul-199.7%
Simplified99.7%
Taylor expanded in y around 0 99.7%
Final simplification98.8%
(FPCore (x y) :precision binary64 (if (<= y -7e+114) (/ (* y (- -1.0 x)) y) (if (<= y 1.15e+104) (- 1.0 x) (/ (* y (+ 1.0 x)) y))))
double code(double x, double y) {
double tmp;
if (y <= -7e+114) {
tmp = (y * (-1.0 - x)) / y;
} else if (y <= 1.15e+104) {
tmp = 1.0 - x;
} else {
tmp = (y * (1.0 + x)) / y;
}
return tmp;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8) :: tmp
if (y <= (-7d+114)) then
tmp = (y * ((-1.0d0) - x)) / y
else if (y <= 1.15d+104) then
tmp = 1.0d0 - x
else
tmp = (y * (1.0d0 + x)) / y
end if
code = tmp
end function
public static double code(double x, double y) {
double tmp;
if (y <= -7e+114) {
tmp = (y * (-1.0 - x)) / y;
} else if (y <= 1.15e+104) {
tmp = 1.0 - x;
} else {
tmp = (y * (1.0 + x)) / y;
}
return tmp;
}
def code(x, y): tmp = 0 if y <= -7e+114: tmp = (y * (-1.0 - x)) / y elif y <= 1.15e+104: tmp = 1.0 - x else: tmp = (y * (1.0 + x)) / y return tmp
function code(x, y) tmp = 0.0 if (y <= -7e+114) tmp = Float64(Float64(y * Float64(-1.0 - x)) / y); elseif (y <= 1.15e+104) tmp = Float64(1.0 - x); else tmp = Float64(Float64(y * Float64(1.0 + x)) / y); end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if (y <= -7e+114) tmp = (y * (-1.0 - x)) / y; elseif (y <= 1.15e+104) tmp = 1.0 - x; else tmp = (y * (1.0 + x)) / y; end tmp_2 = tmp; end
code[x_, y_] := If[LessEqual[y, -7e+114], N[(N[(y * N[(-1.0 - x), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision], If[LessEqual[y, 1.15e+104], N[(1.0 - x), $MachinePrecision], N[(N[(y * N[(1.0 + x), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -7 \cdot 10^{+114}:\\
\;\;\;\;\frac{y \cdot \left(-1 - x\right)}{y}\\
\mathbf{elif}\;y \leq 1.15 \cdot 10^{+104}:\\
\;\;\;\;1 - x\\
\mathbf{else}:\\
\;\;\;\;\frac{y \cdot \left(1 + x\right)}{y}\\
\end{array}
\end{array}
if y < -7.0000000000000001e114Initial program 99.5%
Taylor expanded in y around inf 99.5%
associate--l+99.5%
div-sub99.5%
Simplified99.5%
Taylor expanded in y around 0 2.9%
*-commutative2.9%
frac-2neg2.9%
associate-*l/27.3%
Applied egg-rr29.8%
if -7.0000000000000001e114 < y < 1.14999999999999992e104Initial program 100.0%
Taylor expanded in y around 0 88.2%
if 1.14999999999999992e104 < y Initial program 99.7%
Taylor expanded in y around inf 99.7%
associate--l+99.7%
div-sub99.7%
Simplified99.7%
Taylor expanded in y around 0 3.5%
associate-*r/3.5%
sub-neg3.5%
pow13.5%
metadata-eval3.5%
sqrt-pow116.3%
pow216.3%
sqr-neg16.3%
sqrt-prod16.3%
add-sqr-sqrt16.3%
+-commutative16.3%
Applied egg-rr16.3%
Final simplification67.8%
(FPCore (x y) :precision binary64 (if (<= y -1.7e+114) (* (/ -1.0 y) (* x y)) (if (<= y 1.15e+104) (- 1.0 x) (/ (* y (+ 1.0 x)) y))))
double code(double x, double y) {
double tmp;
if (y <= -1.7e+114) {
tmp = (-1.0 / y) * (x * y);
} else if (y <= 1.15e+104) {
tmp = 1.0 - x;
} else {
tmp = (y * (1.0 + x)) / y;
}
return tmp;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8) :: tmp
if (y <= (-1.7d+114)) then
tmp = ((-1.0d0) / y) * (x * y)
else if (y <= 1.15d+104) then
tmp = 1.0d0 - x
else
tmp = (y * (1.0d0 + x)) / y
end if
code = tmp
end function
public static double code(double x, double y) {
double tmp;
if (y <= -1.7e+114) {
tmp = (-1.0 / y) * (x * y);
} else if (y <= 1.15e+104) {
tmp = 1.0 - x;
} else {
tmp = (y * (1.0 + x)) / y;
}
return tmp;
}
def code(x, y): tmp = 0 if y <= -1.7e+114: tmp = (-1.0 / y) * (x * y) elif y <= 1.15e+104: tmp = 1.0 - x else: tmp = (y * (1.0 + x)) / y return tmp
function code(x, y) tmp = 0.0 if (y <= -1.7e+114) tmp = Float64(Float64(-1.0 / y) * Float64(x * y)); elseif (y <= 1.15e+104) tmp = Float64(1.0 - x); else tmp = Float64(Float64(y * Float64(1.0 + x)) / y); end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if (y <= -1.7e+114) tmp = (-1.0 / y) * (x * y); elseif (y <= 1.15e+104) tmp = 1.0 - x; else tmp = (y * (1.0 + x)) / y; end tmp_2 = tmp; end
code[x_, y_] := If[LessEqual[y, -1.7e+114], N[(N[(-1.0 / y), $MachinePrecision] * N[(x * y), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 1.15e+104], N[(1.0 - x), $MachinePrecision], N[(N[(y * N[(1.0 + x), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -1.7 \cdot 10^{+114}:\\
\;\;\;\;\frac{-1}{y} \cdot \left(x \cdot y\right)\\
\mathbf{elif}\;y \leq 1.15 \cdot 10^{+104}:\\
\;\;\;\;1 - x\\
\mathbf{else}:\\
\;\;\;\;\frac{y \cdot \left(1 + x\right)}{y}\\
\end{array}
\end{array}
if y < -1.7e114Initial program 99.5%
Taylor expanded in y around inf 99.5%
associate--l+99.5%
div-sub99.5%
Simplified99.5%
Taylor expanded in x around inf 75.2%
associate-*r*97.0%
*-commutative97.0%
sub-neg97.0%
distribute-neg-frac97.0%
metadata-eval97.0%
Simplified97.0%
Taylor expanded in y around 0 28.6%
if -1.7e114 < y < 1.14999999999999992e104Initial program 100.0%
Taylor expanded in y around 0 88.2%
if 1.14999999999999992e104 < y Initial program 99.7%
Taylor expanded in y around inf 99.7%
associate--l+99.7%
div-sub99.7%
Simplified99.7%
Taylor expanded in y around 0 3.5%
associate-*r/3.5%
sub-neg3.5%
pow13.5%
metadata-eval3.5%
sqrt-pow116.3%
pow216.3%
sqr-neg16.3%
sqrt-prod16.3%
add-sqr-sqrt16.3%
+-commutative16.3%
Applied egg-rr16.3%
Final simplification67.7%
(FPCore (x y) :precision binary64 (if (<= y -7.6e+114) (* (/ -1.0 y) (* x y)) (- 1.0 x)))
double code(double x, double y) {
double tmp;
if (y <= -7.6e+114) {
tmp = (-1.0 / y) * (x * y);
} else {
tmp = 1.0 - x;
}
return tmp;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8) :: tmp
if (y <= (-7.6d+114)) then
tmp = ((-1.0d0) / y) * (x * y)
else
tmp = 1.0d0 - x
end if
code = tmp
end function
public static double code(double x, double y) {
double tmp;
if (y <= -7.6e+114) {
tmp = (-1.0 / y) * (x * y);
} else {
tmp = 1.0 - x;
}
return tmp;
}
def code(x, y): tmp = 0 if y <= -7.6e+114: tmp = (-1.0 / y) * (x * y) else: tmp = 1.0 - x return tmp
function code(x, y) tmp = 0.0 if (y <= -7.6e+114) tmp = Float64(Float64(-1.0 / y) * Float64(x * y)); else tmp = Float64(1.0 - x); end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if (y <= -7.6e+114) tmp = (-1.0 / y) * (x * y); else tmp = 1.0 - x; end tmp_2 = tmp; end
code[x_, y_] := If[LessEqual[y, -7.6e+114], N[(N[(-1.0 / y), $MachinePrecision] * N[(x * y), $MachinePrecision]), $MachinePrecision], N[(1.0 - x), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -7.6 \cdot 10^{+114}:\\
\;\;\;\;\frac{-1}{y} \cdot \left(x \cdot y\right)\\
\mathbf{else}:\\
\;\;\;\;1 - x\\
\end{array}
\end{array}
if y < -7.6000000000000001e114Initial program 99.5%
Taylor expanded in y around inf 99.5%
associate--l+99.5%
div-sub99.5%
Simplified99.5%
Taylor expanded in x around inf 75.2%
associate-*r*97.0%
*-commutative97.0%
sub-neg97.0%
distribute-neg-frac97.0%
metadata-eval97.0%
Simplified97.0%
Taylor expanded in y around 0 28.6%
if -7.6000000000000001e114 < y Initial program 99.9%
Taylor expanded in y around 0 71.3%
(FPCore (x y) :precision binary64 (if (<= x 65000000000.0) 1.0 (- x)))
double code(double x, double y) {
double tmp;
if (x <= 65000000000.0) {
tmp = 1.0;
} else {
tmp = -x;
}
return tmp;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8) :: tmp
if (x <= 65000000000.0d0) then
tmp = 1.0d0
else
tmp = -x
end if
code = tmp
end function
public static double code(double x, double y) {
double tmp;
if (x <= 65000000000.0) {
tmp = 1.0;
} else {
tmp = -x;
}
return tmp;
}
def code(x, y): tmp = 0 if x <= 65000000000.0: tmp = 1.0 else: tmp = -x return tmp
function code(x, y) tmp = 0.0 if (x <= 65000000000.0) tmp = 1.0; else tmp = Float64(-x); end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if (x <= 65000000000.0) tmp = 1.0; else tmp = -x; end tmp_2 = tmp; end
code[x_, y_] := If[LessEqual[x, 65000000000.0], 1.0, (-x)]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 65000000000:\\
\;\;\;\;1\\
\mathbf{else}:\\
\;\;\;\;-x\\
\end{array}
\end{array}
if x < 6.5e10Initial program 99.8%
Taylor expanded in y around 0 58.9%
Taylor expanded in x around 0 57.3%
if 6.5e10 < x Initial program 99.9%
Taylor expanded in x around inf 99.7%
neg-mul-199.7%
Simplified99.7%
Taylor expanded in y around 0 65.3%
mul-1-neg65.3%
Simplified65.3%
(FPCore (x y) :precision binary64 (- 1.0 x))
double code(double x, double y) {
return 1.0 - x;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = 1.0d0 - x
end function
public static double code(double x, double y) {
return 1.0 - x;
}
def code(x, y): return 1.0 - x
function code(x, y) return Float64(1.0 - x) end
function tmp = code(x, y) tmp = 1.0 - x; end
code[x_, y_] := N[(1.0 - x), $MachinePrecision]
\begin{array}{l}
\\
1 - x
\end{array}
Initial program 99.9%
Taylor expanded in y around 0 62.0%
(FPCore (x y) :precision binary64 1.0)
double code(double x, double y) {
return 1.0;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = 1.0d0
end function
public static double code(double x, double y) {
return 1.0;
}
def code(x, y): return 1.0
function code(x, y) return 1.0 end
function tmp = code(x, y) tmp = 1.0; end
code[x_, y_] := 1.0
\begin{array}{l}
\\
1
\end{array}
Initial program 99.9%
Taylor expanded in y around 0 62.0%
Taylor expanded in x around 0 31.4%
herbie shell --seed 2024137
(FPCore (x y)
:name "Numeric.SpecFunctions:invIncompleteBetaWorker from math-functions-0.1.5.2, E"
:precision binary64
(+ (- 1.0 x) (* y (sqrt x))))