
(FPCore (x) :precision binary64 (- (/ PI 2.0) (* 2.0 (asin (sqrt (/ (- 1.0 x) 2.0))))))
double code(double x) {
return (((double) M_PI) / 2.0) - (2.0 * asin(sqrt(((1.0 - x) / 2.0))));
}
public static double code(double x) {
return (Math.PI / 2.0) - (2.0 * Math.asin(Math.sqrt(((1.0 - x) / 2.0))));
}
def code(x): return (math.pi / 2.0) - (2.0 * math.asin(math.sqrt(((1.0 - x) / 2.0))))
function code(x) return Float64(Float64(pi / 2.0) - Float64(2.0 * asin(sqrt(Float64(Float64(1.0 - x) / 2.0))))) end
function tmp = code(x) tmp = (pi / 2.0) - (2.0 * asin(sqrt(((1.0 - x) / 2.0)))); end
code[x_] := N[(N[(Pi / 2.0), $MachinePrecision] - N[(2.0 * N[ArcSin[N[Sqrt[N[(N[(1.0 - x), $MachinePrecision] / 2.0), $MachinePrecision]], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\pi}{2} - 2 \cdot \sin^{-1} \left(\sqrt{\frac{1 - x}{2}}\right)
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 11 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x) :precision binary64 (- (/ PI 2.0) (* 2.0 (asin (sqrt (/ (- 1.0 x) 2.0))))))
double code(double x) {
return (((double) M_PI) / 2.0) - (2.0 * asin(sqrt(((1.0 - x) / 2.0))));
}
public static double code(double x) {
return (Math.PI / 2.0) - (2.0 * Math.asin(Math.sqrt(((1.0 - x) / 2.0))));
}
def code(x): return (math.pi / 2.0) - (2.0 * math.asin(math.sqrt(((1.0 - x) / 2.0))))
function code(x) return Float64(Float64(pi / 2.0) - Float64(2.0 * asin(sqrt(Float64(Float64(1.0 - x) / 2.0))))) end
function tmp = code(x) tmp = (pi / 2.0) - (2.0 * asin(sqrt(((1.0 - x) / 2.0)))); end
code[x_] := N[(N[(Pi / 2.0), $MachinePrecision] - N[(2.0 * N[ArcSin[N[Sqrt[N[(N[(1.0 - x), $MachinePrecision] / 2.0), $MachinePrecision]], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\pi}{2} - 2 \cdot \sin^{-1} \left(\sqrt{\frac{1 - x}{2}}\right)
\end{array}
(FPCore (x)
:precision binary64
(/
1.0
(/
(fma 2.0 (asin (sqrt (- 0.5 (* 0.5 x)))) (* 0.5 (cbrt (pow PI 3.0))))
(log
(exp
(-
(* (pow PI 2.0) 0.25)
(* 4.0 (expm1 (log1p (pow (asin (sqrt (+ 0.5 (* x -0.5)))) 2.0))))))))))
double code(double x) {
return 1.0 / (fma(2.0, asin(sqrt((0.5 - (0.5 * x)))), (0.5 * cbrt(pow(((double) M_PI), 3.0)))) / log(exp(((pow(((double) M_PI), 2.0) * 0.25) - (4.0 * expm1(log1p(pow(asin(sqrt((0.5 + (x * -0.5)))), 2.0))))))));
}
function code(x) return Float64(1.0 / Float64(fma(2.0, asin(sqrt(Float64(0.5 - Float64(0.5 * x)))), Float64(0.5 * cbrt((pi ^ 3.0)))) / log(exp(Float64(Float64((pi ^ 2.0) * 0.25) - Float64(4.0 * expm1(log1p((asin(sqrt(Float64(0.5 + Float64(x * -0.5)))) ^ 2.0))))))))) end
code[x_] := N[(1.0 / N[(N[(2.0 * N[ArcSin[N[Sqrt[N[(0.5 - N[(0.5 * x), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], $MachinePrecision] + N[(0.5 * N[Power[N[Power[Pi, 3.0], $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[Log[N[Exp[N[(N[(N[Power[Pi, 2.0], $MachinePrecision] * 0.25), $MachinePrecision] - N[(4.0 * N[(Exp[N[Log[1 + N[Power[N[ArcSin[N[Sqrt[N[(0.5 + N[(x * -0.5), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], $MachinePrecision], 2.0], $MachinePrecision]], $MachinePrecision]] - 1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{1}{\frac{\mathsf{fma}\left(2, \sin^{-1} \left(\sqrt{0.5 - 0.5 \cdot x}\right), 0.5 \cdot \sqrt[3]{{\pi}^{3}}\right)}{\log \left(e^{{\pi}^{2} \cdot 0.25 - 4 \cdot \mathsf{expm1}\left(\mathsf{log1p}\left({\sin^{-1} \left(\sqrt{0.5 + x \cdot -0.5}\right)}^{2}\right)\right)}\right)}}
\end{array}
Initial program 6.6%
flip--6.5%
clear-num6.5%
Applied egg-rr6.5%
add-log-exp6.5%
unpow-prod-down6.5%
metadata-eval6.5%
unpow-prod-down6.5%
metadata-eval6.5%
*-commutative6.5%
Applied egg-rr6.5%
expm1-log1p-u8.3%
expm1-undefine8.3%
add-sqr-sqrt8.3%
add-sqr-sqrt8.3%
cancel-sign-sub-inv8.3%
metadata-eval8.3%
Applied egg-rr8.3%
expm1-define8.3%
*-commutative8.3%
Simplified8.3%
add-cbrt-cube8.3%
pow38.3%
Applied egg-rr8.3%
Final simplification8.3%
(FPCore (x)
:precision binary64
(/
1.0
(/
(fma 2.0 (- (* 0.5 PI) (acos (sqrt (- 0.5 (* 0.5 x))))) (* 0.5 PI))
(log
(exp
(-
(* (pow PI 2.0) 0.25)
(* 4.0 (expm1 (log1p (pow (asin (sqrt (+ 0.5 (* x -0.5)))) 2.0))))))))))
double code(double x) {
return 1.0 / (fma(2.0, ((0.5 * ((double) M_PI)) - acos(sqrt((0.5 - (0.5 * x))))), (0.5 * ((double) M_PI))) / log(exp(((pow(((double) M_PI), 2.0) * 0.25) - (4.0 * expm1(log1p(pow(asin(sqrt((0.5 + (x * -0.5)))), 2.0))))))));
}
function code(x) return Float64(1.0 / Float64(fma(2.0, Float64(Float64(0.5 * pi) - acos(sqrt(Float64(0.5 - Float64(0.5 * x))))), Float64(0.5 * pi)) / log(exp(Float64(Float64((pi ^ 2.0) * 0.25) - Float64(4.0 * expm1(log1p((asin(sqrt(Float64(0.5 + Float64(x * -0.5)))) ^ 2.0))))))))) end
code[x_] := N[(1.0 / N[(N[(2.0 * N[(N[(0.5 * Pi), $MachinePrecision] - N[ArcCos[N[Sqrt[N[(0.5 - N[(0.5 * x), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], $MachinePrecision]), $MachinePrecision] + N[(0.5 * Pi), $MachinePrecision]), $MachinePrecision] / N[Log[N[Exp[N[(N[(N[Power[Pi, 2.0], $MachinePrecision] * 0.25), $MachinePrecision] - N[(4.0 * N[(Exp[N[Log[1 + N[Power[N[ArcSin[N[Sqrt[N[(0.5 + N[(x * -0.5), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], $MachinePrecision], 2.0], $MachinePrecision]], $MachinePrecision]] - 1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{1}{\frac{\mathsf{fma}\left(2, 0.5 \cdot \pi - \cos^{-1} \left(\sqrt{0.5 - 0.5 \cdot x}\right), 0.5 \cdot \pi\right)}{\log \left(e^{{\pi}^{2} \cdot 0.25 - 4 \cdot \mathsf{expm1}\left(\mathsf{log1p}\left({\sin^{-1} \left(\sqrt{0.5 + x \cdot -0.5}\right)}^{2}\right)\right)}\right)}}
\end{array}
Initial program 6.6%
flip--6.5%
clear-num6.5%
Applied egg-rr6.5%
add-log-exp6.5%
unpow-prod-down6.5%
metadata-eval6.5%
unpow-prod-down6.5%
metadata-eval6.5%
*-commutative6.5%
Applied egg-rr6.5%
expm1-log1p-u8.3%
expm1-undefine8.3%
add-sqr-sqrt8.3%
add-sqr-sqrt8.3%
cancel-sign-sub-inv8.3%
metadata-eval8.3%
Applied egg-rr8.3%
expm1-define8.3%
*-commutative8.3%
Simplified8.3%
asin-acos8.2%
div-inv8.2%
metadata-eval8.2%
*-commutative8.2%
Applied egg-rr8.3%
Final simplification8.3%
(FPCore (x)
:precision binary64
(/
1.0
(/
(fma 2.0 (asin (sqrt (- 0.5 (* 0.5 x)))) (* 0.5 PI))
(log
(exp
(-
(* (pow PI 2.0) 0.25)
(* 4.0 (expm1 (log1p (pow (asin (sqrt (+ 0.5 (* x -0.5)))) 2.0))))))))))
double code(double x) {
return 1.0 / (fma(2.0, asin(sqrt((0.5 - (0.5 * x)))), (0.5 * ((double) M_PI))) / log(exp(((pow(((double) M_PI), 2.0) * 0.25) - (4.0 * expm1(log1p(pow(asin(sqrt((0.5 + (x * -0.5)))), 2.0))))))));
}
function code(x) return Float64(1.0 / Float64(fma(2.0, asin(sqrt(Float64(0.5 - Float64(0.5 * x)))), Float64(0.5 * pi)) / log(exp(Float64(Float64((pi ^ 2.0) * 0.25) - Float64(4.0 * expm1(log1p((asin(sqrt(Float64(0.5 + Float64(x * -0.5)))) ^ 2.0))))))))) end
code[x_] := N[(1.0 / N[(N[(2.0 * N[ArcSin[N[Sqrt[N[(0.5 - N[(0.5 * x), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], $MachinePrecision] + N[(0.5 * Pi), $MachinePrecision]), $MachinePrecision] / N[Log[N[Exp[N[(N[(N[Power[Pi, 2.0], $MachinePrecision] * 0.25), $MachinePrecision] - N[(4.0 * N[(Exp[N[Log[1 + N[Power[N[ArcSin[N[Sqrt[N[(0.5 + N[(x * -0.5), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], $MachinePrecision], 2.0], $MachinePrecision]], $MachinePrecision]] - 1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{1}{\frac{\mathsf{fma}\left(2, \sin^{-1} \left(\sqrt{0.5 - 0.5 \cdot x}\right), 0.5 \cdot \pi\right)}{\log \left(e^{{\pi}^{2} \cdot 0.25 - 4 \cdot \mathsf{expm1}\left(\mathsf{log1p}\left({\sin^{-1} \left(\sqrt{0.5 + x \cdot -0.5}\right)}^{2}\right)\right)}\right)}}
\end{array}
Initial program 6.6%
flip--6.5%
clear-num6.5%
Applied egg-rr6.5%
add-log-exp6.5%
unpow-prod-down6.5%
metadata-eval6.5%
unpow-prod-down6.5%
metadata-eval6.5%
*-commutative6.5%
Applied egg-rr6.5%
expm1-log1p-u8.3%
expm1-undefine8.3%
add-sqr-sqrt8.3%
add-sqr-sqrt8.3%
cancel-sign-sub-inv8.3%
metadata-eval8.3%
Applied egg-rr8.3%
expm1-define8.3%
*-commutative8.3%
Simplified8.3%
Final simplification8.3%
(FPCore (x) :precision binary64 (/ (- (* (pow PI 2.0) 0.25) (* 4.0 (expm1 (log1p (pow (asin (sqrt (+ 0.5 (* x -0.5)))) 2.0))))) (+ (* 0.5 PI) (* 2.0 (asin (sqrt (- 0.5 (* 0.5 x))))))))
double code(double x) {
return ((pow(((double) M_PI), 2.0) * 0.25) - (4.0 * expm1(log1p(pow(asin(sqrt((0.5 + (x * -0.5)))), 2.0))))) / ((0.5 * ((double) M_PI)) + (2.0 * asin(sqrt((0.5 - (0.5 * x))))));
}
public static double code(double x) {
return ((Math.pow(Math.PI, 2.0) * 0.25) - (4.0 * Math.expm1(Math.log1p(Math.pow(Math.asin(Math.sqrt((0.5 + (x * -0.5)))), 2.0))))) / ((0.5 * Math.PI) + (2.0 * Math.asin(Math.sqrt((0.5 - (0.5 * x))))));
}
def code(x): return ((math.pow(math.pi, 2.0) * 0.25) - (4.0 * math.expm1(math.log1p(math.pow(math.asin(math.sqrt((0.5 + (x * -0.5)))), 2.0))))) / ((0.5 * math.pi) + (2.0 * math.asin(math.sqrt((0.5 - (0.5 * x))))))
function code(x) return Float64(Float64(Float64((pi ^ 2.0) * 0.25) - Float64(4.0 * expm1(log1p((asin(sqrt(Float64(0.5 + Float64(x * -0.5)))) ^ 2.0))))) / Float64(Float64(0.5 * pi) + Float64(2.0 * asin(sqrt(Float64(0.5 - Float64(0.5 * x))))))) end
code[x_] := N[(N[(N[(N[Power[Pi, 2.0], $MachinePrecision] * 0.25), $MachinePrecision] - N[(4.0 * N[(Exp[N[Log[1 + N[Power[N[ArcSin[N[Sqrt[N[(0.5 + N[(x * -0.5), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], $MachinePrecision], 2.0], $MachinePrecision]], $MachinePrecision]] - 1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(N[(0.5 * Pi), $MachinePrecision] + N[(2.0 * N[ArcSin[N[Sqrt[N[(0.5 - N[(0.5 * x), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{{\pi}^{2} \cdot 0.25 - 4 \cdot \mathsf{expm1}\left(\mathsf{log1p}\left({\sin^{-1} \left(\sqrt{0.5 + x \cdot -0.5}\right)}^{2}\right)\right)}{0.5 \cdot \pi + 2 \cdot \sin^{-1} \left(\sqrt{0.5 - 0.5 \cdot x}\right)}
\end{array}
Initial program 6.6%
flip--6.5%
clear-num6.5%
Applied egg-rr6.5%
Taylor expanded in x around 0 6.5%
expm1-log1p-u8.3%
expm1-undefine8.3%
add-sqr-sqrt8.3%
add-sqr-sqrt8.3%
cancel-sign-sub-inv8.3%
metadata-eval8.3%
Applied egg-rr8.2%
expm1-define8.3%
*-commutative8.3%
Simplified8.2%
Final simplification8.2%
(FPCore (x)
:precision binary64
(let* ((t_0 (sqrt (- 0.5 (* 0.5 x)))))
(/
1.0
(/
(fma 2.0 (asin t_0) (* 0.5 PI))
(- (pow (* 0.5 PI) 2.0) (pow (* 2.0 (- (* 0.5 PI) (acos t_0))) 2.0))))))
double code(double x) {
double t_0 = sqrt((0.5 - (0.5 * x)));
return 1.0 / (fma(2.0, asin(t_0), (0.5 * ((double) M_PI))) / (pow((0.5 * ((double) M_PI)), 2.0) - pow((2.0 * ((0.5 * ((double) M_PI)) - acos(t_0))), 2.0)));
}
function code(x) t_0 = sqrt(Float64(0.5 - Float64(0.5 * x))) return Float64(1.0 / Float64(fma(2.0, asin(t_0), Float64(0.5 * pi)) / Float64((Float64(0.5 * pi) ^ 2.0) - (Float64(2.0 * Float64(Float64(0.5 * pi) - acos(t_0))) ^ 2.0)))) end
code[x_] := Block[{t$95$0 = N[Sqrt[N[(0.5 - N[(0.5 * x), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]}, N[(1.0 / N[(N[(2.0 * N[ArcSin[t$95$0], $MachinePrecision] + N[(0.5 * Pi), $MachinePrecision]), $MachinePrecision] / N[(N[Power[N[(0.5 * Pi), $MachinePrecision], 2.0], $MachinePrecision] - N[Power[N[(2.0 * N[(N[(0.5 * Pi), $MachinePrecision] - N[ArcCos[t$95$0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sqrt{0.5 - 0.5 \cdot x}\\
\frac{1}{\frac{\mathsf{fma}\left(2, \sin^{-1} t\_0, 0.5 \cdot \pi\right)}{{\left(0.5 \cdot \pi\right)}^{2} - {\left(2 \cdot \left(0.5 \cdot \pi - \cos^{-1} t\_0\right)\right)}^{2}}}
\end{array}
\end{array}
Initial program 6.6%
flip--6.5%
clear-num6.5%
Applied egg-rr6.5%
asin-acos8.2%
div-inv8.2%
metadata-eval8.2%
*-commutative8.2%
Applied egg-rr8.2%
Final simplification8.2%
(FPCore (x) :precision binary64 (/ (- (* (pow PI 2.0) 0.25) (* 4.0 (pow (- (* 0.5 PI) (acos (sqrt (+ 0.5 (* x -0.5))))) 2.0))) (+ (* 0.5 PI) (* 2.0 (asin (sqrt (- 0.5 (* 0.5 x))))))))
double code(double x) {
return ((pow(((double) M_PI), 2.0) * 0.25) - (4.0 * pow(((0.5 * ((double) M_PI)) - acos(sqrt((0.5 + (x * -0.5))))), 2.0))) / ((0.5 * ((double) M_PI)) + (2.0 * asin(sqrt((0.5 - (0.5 * x))))));
}
public static double code(double x) {
return ((Math.pow(Math.PI, 2.0) * 0.25) - (4.0 * Math.pow(((0.5 * Math.PI) - Math.acos(Math.sqrt((0.5 + (x * -0.5))))), 2.0))) / ((0.5 * Math.PI) + (2.0 * Math.asin(Math.sqrt((0.5 - (0.5 * x))))));
}
def code(x): return ((math.pow(math.pi, 2.0) * 0.25) - (4.0 * math.pow(((0.5 * math.pi) - math.acos(math.sqrt((0.5 + (x * -0.5))))), 2.0))) / ((0.5 * math.pi) + (2.0 * math.asin(math.sqrt((0.5 - (0.5 * x))))))
function code(x) return Float64(Float64(Float64((pi ^ 2.0) * 0.25) - Float64(4.0 * (Float64(Float64(0.5 * pi) - acos(sqrt(Float64(0.5 + Float64(x * -0.5))))) ^ 2.0))) / Float64(Float64(0.5 * pi) + Float64(2.0 * asin(sqrt(Float64(0.5 - Float64(0.5 * x))))))) end
function tmp = code(x) tmp = (((pi ^ 2.0) * 0.25) - (4.0 * (((0.5 * pi) - acos(sqrt((0.5 + (x * -0.5))))) ^ 2.0))) / ((0.5 * pi) + (2.0 * asin(sqrt((0.5 - (0.5 * x)))))); end
code[x_] := N[(N[(N[(N[Power[Pi, 2.0], $MachinePrecision] * 0.25), $MachinePrecision] - N[(4.0 * N[Power[N[(N[(0.5 * Pi), $MachinePrecision] - N[ArcCos[N[Sqrt[N[(0.5 + N[(x * -0.5), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(N[(0.5 * Pi), $MachinePrecision] + N[(2.0 * N[ArcSin[N[Sqrt[N[(0.5 - N[(0.5 * x), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{{\pi}^{2} \cdot 0.25 - 4 \cdot {\left(0.5 \cdot \pi - \cos^{-1} \left(\sqrt{0.5 + x \cdot -0.5}\right)\right)}^{2}}{0.5 \cdot \pi + 2 \cdot \sin^{-1} \left(\sqrt{0.5 - 0.5 \cdot x}\right)}
\end{array}
Initial program 6.6%
flip--6.5%
clear-num6.5%
Applied egg-rr6.5%
Taylor expanded in x around 0 6.5%
asin-acos8.2%
div-inv8.2%
metadata-eval8.2%
sub-neg8.2%
cancel-sign-sub-inv8.2%
metadata-eval8.2%
Applied egg-rr8.2%
sub-neg8.2%
*-commutative8.2%
Simplified8.2%
Final simplification8.2%
(FPCore (x) :precision binary64 (+ (/ PI 2.0) (* 2.0 (- (acos (sqrt (+ 0.5 (* x -0.5)))) (/ PI 2.0)))))
double code(double x) {
return (((double) M_PI) / 2.0) + (2.0 * (acos(sqrt((0.5 + (x * -0.5)))) - (((double) M_PI) / 2.0)));
}
public static double code(double x) {
return (Math.PI / 2.0) + (2.0 * (Math.acos(Math.sqrt((0.5 + (x * -0.5)))) - (Math.PI / 2.0)));
}
def code(x): return (math.pi / 2.0) + (2.0 * (math.acos(math.sqrt((0.5 + (x * -0.5)))) - (math.pi / 2.0)))
function code(x) return Float64(Float64(pi / 2.0) + Float64(2.0 * Float64(acos(sqrt(Float64(0.5 + Float64(x * -0.5)))) - Float64(pi / 2.0)))) end
function tmp = code(x) tmp = (pi / 2.0) + (2.0 * (acos(sqrt((0.5 + (x * -0.5)))) - (pi / 2.0))); end
code[x_] := N[(N[(Pi / 2.0), $MachinePrecision] + N[(2.0 * N[(N[ArcCos[N[Sqrt[N[(0.5 + N[(x * -0.5), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], $MachinePrecision] - N[(Pi / 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\pi}{2} + 2 \cdot \left(\cos^{-1} \left(\sqrt{0.5 + x \cdot -0.5}\right) - \frac{\pi}{2}\right)
\end{array}
Initial program 6.6%
asin-acos8.2%
add-cube-cbrt6.3%
associate-/l*6.3%
fmm-def6.3%
pow26.3%
div-sub6.3%
metadata-eval6.3%
div-inv6.3%
metadata-eval6.3%
Applied egg-rr6.3%
fmm-undef6.3%
associate-*r/6.3%
unpow26.3%
rem-3cbrt-lft8.2%
sub-neg8.2%
distribute-rgt-neg-in8.2%
metadata-eval8.2%
Simplified8.2%
Final simplification8.2%
(FPCore (x) :precision binary64 (- (/ PI 2.0) (* 2.0 (asin (/ 1.0 (sqrt (/ 2.0 (- 1.0 x))))))))
double code(double x) {
return (((double) M_PI) / 2.0) - (2.0 * asin((1.0 / sqrt((2.0 / (1.0 - x))))));
}
public static double code(double x) {
return (Math.PI / 2.0) - (2.0 * Math.asin((1.0 / Math.sqrt((2.0 / (1.0 - x))))));
}
def code(x): return (math.pi / 2.0) - (2.0 * math.asin((1.0 / math.sqrt((2.0 / (1.0 - x))))))
function code(x) return Float64(Float64(pi / 2.0) - Float64(2.0 * asin(Float64(1.0 / sqrt(Float64(2.0 / Float64(1.0 - x))))))) end
function tmp = code(x) tmp = (pi / 2.0) - (2.0 * asin((1.0 / sqrt((2.0 / (1.0 - x)))))); end
code[x_] := N[(N[(Pi / 2.0), $MachinePrecision] - N[(2.0 * N[ArcSin[N[(1.0 / N[Sqrt[N[(2.0 / N[(1.0 - x), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\pi}{2} - 2 \cdot \sin^{-1} \left(\frac{1}{\sqrt{\frac{2}{1 - x}}}\right)
\end{array}
Initial program 6.6%
clear-num6.5%
sqrt-div6.7%
metadata-eval6.7%
Applied egg-rr6.7%
Final simplification6.7%
(FPCore (x) :precision binary64 (- (/ PI 2.0) (* 2.0 (asin (sqrt (/ (- 1.0 x) 2.0))))))
double code(double x) {
return (((double) M_PI) / 2.0) - (2.0 * asin(sqrt(((1.0 - x) / 2.0))));
}
public static double code(double x) {
return (Math.PI / 2.0) - (2.0 * Math.asin(Math.sqrt(((1.0 - x) / 2.0))));
}
def code(x): return (math.pi / 2.0) - (2.0 * math.asin(math.sqrt(((1.0 - x) / 2.0))))
function code(x) return Float64(Float64(pi / 2.0) - Float64(2.0 * asin(sqrt(Float64(Float64(1.0 - x) / 2.0))))) end
function tmp = code(x) tmp = (pi / 2.0) - (2.0 * asin(sqrt(((1.0 - x) / 2.0)))); end
code[x_] := N[(N[(Pi / 2.0), $MachinePrecision] - N[(2.0 * N[ArcSin[N[Sqrt[N[(N[(1.0 - x), $MachinePrecision] / 2.0), $MachinePrecision]], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\pi}{2} - 2 \cdot \sin^{-1} \left(\sqrt{\frac{1 - x}{2}}\right)
\end{array}
Initial program 6.6%
Final simplification6.6%
(FPCore (x) :precision binary64 (+ (* 0.5 PI) (* 2.0 (asin (sqrt 0.5)))))
double code(double x) {
return (0.5 * ((double) M_PI)) + (2.0 * asin(sqrt(0.5)));
}
public static double code(double x) {
return (0.5 * Math.PI) + (2.0 * Math.asin(Math.sqrt(0.5)));
}
def code(x): return (0.5 * math.pi) + (2.0 * math.asin(math.sqrt(0.5)))
function code(x) return Float64(Float64(0.5 * pi) + Float64(2.0 * asin(sqrt(0.5)))) end
function tmp = code(x) tmp = (0.5 * pi) + (2.0 * asin(sqrt(0.5))); end
code[x_] := N[(N[(0.5 * Pi), $MachinePrecision] + N[(2.0 * N[ArcSin[N[Sqrt[0.5], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
0.5 \cdot \pi + 2 \cdot \sin^{-1} \left(\sqrt{0.5}\right)
\end{array}
Initial program 6.6%
add-cube-cbrt6.6%
pow36.6%
Applied egg-rr6.6%
rem-cube-cbrt6.6%
fma-undefine6.6%
add-sqr-sqrt0.0%
sqrt-unprod3.7%
swap-sqr3.7%
pow23.7%
*-commutative3.7%
metadata-eval3.7%
sqrt-prod3.7%
Applied egg-rr3.7%
Taylor expanded in x around 0 3.7%
Final simplification3.7%
(FPCore (x) :precision binary64 (- (/ PI 2.0) (* 2.0 (asin (sqrt 0.5)))))
double code(double x) {
return (((double) M_PI) / 2.0) - (2.0 * asin(sqrt(0.5)));
}
public static double code(double x) {
return (Math.PI / 2.0) - (2.0 * Math.asin(Math.sqrt(0.5)));
}
def code(x): return (math.pi / 2.0) - (2.0 * math.asin(math.sqrt(0.5)))
function code(x) return Float64(Float64(pi / 2.0) - Float64(2.0 * asin(sqrt(0.5)))) end
function tmp = code(x) tmp = (pi / 2.0) - (2.0 * asin(sqrt(0.5))); end
code[x_] := N[(N[(Pi / 2.0), $MachinePrecision] - N[(2.0 * N[ArcSin[N[Sqrt[0.5], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\pi}{2} - 2 \cdot \sin^{-1} \left(\sqrt{0.5}\right)
\end{array}
Initial program 6.6%
Taylor expanded in x around 0 3.9%
Final simplification3.9%
(FPCore (x) :precision binary64 (asin x))
double code(double x) {
return asin(x);
}
real(8) function code(x)
real(8), intent (in) :: x
code = asin(x)
end function
public static double code(double x) {
return Math.asin(x);
}
def code(x): return math.asin(x)
function code(x) return asin(x) end
function tmp = code(x) tmp = asin(x); end
code[x_] := N[ArcSin[x], $MachinePrecision]
\begin{array}{l}
\\
\sin^{-1} x
\end{array}
herbie shell --seed 2024112
(FPCore (x)
:name "Ian Simplification"
:precision binary64
:alt
(asin x)
(- (/ PI 2.0) (* 2.0 (asin (sqrt (/ (- 1.0 x) 2.0))))))