
(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 5 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 (+ (/ PI 2.0) (* 2.0 (- (acos (* (sqrt (- 1.0 x)) (sqrt 0.5))) (/ PI 2.0)))))
double code(double x) {
return (((double) M_PI) / 2.0) + (2.0 * (acos((sqrt((1.0 - x)) * sqrt(0.5))) - (((double) M_PI) / 2.0)));
}
public static double code(double x) {
return (Math.PI / 2.0) + (2.0 * (Math.acos((Math.sqrt((1.0 - x)) * Math.sqrt(0.5))) - (Math.PI / 2.0)));
}
def code(x): return (math.pi / 2.0) + (2.0 * (math.acos((math.sqrt((1.0 - x)) * math.sqrt(0.5))) - (math.pi / 2.0)))
function code(x) return Float64(Float64(pi / 2.0) + Float64(2.0 * Float64(acos(Float64(sqrt(Float64(1.0 - x)) * sqrt(0.5))) - Float64(pi / 2.0)))) end
function tmp = code(x) tmp = (pi / 2.0) + (2.0 * (acos((sqrt((1.0 - x)) * sqrt(0.5))) - (pi / 2.0))); end
code[x_] := N[(N[(Pi / 2.0), $MachinePrecision] + N[(2.0 * N[(N[ArcCos[N[(N[Sqrt[N[(1.0 - x), $MachinePrecision]], $MachinePrecision] * N[Sqrt[0.5], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] - N[(Pi / 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\pi}{2} + 2 \cdot \left(\cos^{-1} \left(\sqrt{1 - x} \cdot \sqrt{0.5}\right) - \frac{\pi}{2}\right)
\end{array}
Initial program 7.9%
asin-acos9.1%
add-sqr-sqrt7.4%
associate-/l*7.4%
fma-neg7.4%
div-inv7.4%
metadata-eval7.4%
Applied egg-rr7.4%
fma-neg7.4%
associate-*r/7.4%
rem-square-sqrt9.1%
*-commutative9.1%
Simplified9.1%
*-commutative9.1%
sqrt-prod9.1%
Applied egg-rr9.1%
Final simplification9.1%
(FPCore (x) :precision binary64 (+ (/ PI 2.0) (* 2.0 (- (acos (sqrt (* (- 1.0 x) 0.5))) (/ PI 2.0)))))
double code(double x) {
return (((double) M_PI) / 2.0) + (2.0 * (acos(sqrt(((1.0 - 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(((1.0 - x) * 0.5))) - (Math.PI / 2.0)));
}
def code(x): return (math.pi / 2.0) + (2.0 * (math.acos(math.sqrt(((1.0 - x) * 0.5))) - (math.pi / 2.0)))
function code(x) return Float64(Float64(pi / 2.0) + Float64(2.0 * Float64(acos(sqrt(Float64(Float64(1.0 - x) * 0.5))) - Float64(pi / 2.0)))) end
function tmp = code(x) tmp = (pi / 2.0) + (2.0 * (acos(sqrt(((1.0 - x) * 0.5))) - (pi / 2.0))); end
code[x_] := N[(N[(Pi / 2.0), $MachinePrecision] + N[(2.0 * N[(N[ArcCos[N[Sqrt[N[(N[(1.0 - x), $MachinePrecision] * 0.5), $MachinePrecision]], $MachinePrecision]], $MachinePrecision] - N[(Pi / 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\pi}{2} + 2 \cdot \left(\cos^{-1} \left(\sqrt{\left(1 - x\right) \cdot 0.5}\right) - \frac{\pi}{2}\right)
\end{array}
Initial program 7.9%
asin-acos9.1%
add-sqr-sqrt7.4%
associate-/l*7.4%
fma-neg7.4%
div-inv7.4%
metadata-eval7.4%
Applied egg-rr7.4%
fma-neg7.4%
associate-*r/7.4%
rem-square-sqrt9.1%
*-commutative9.1%
Simplified9.1%
Final simplification9.1%
(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 7.9%
(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 7.9%
Taylor expanded in x around 0 4.2%
(FPCore (x) :precision binary64 (+ (* 2.0 (asin (sqrt 0.5))) (* PI 0.5)))
double code(double x) {
return (2.0 * asin(sqrt(0.5))) + (((double) M_PI) * 0.5);
}
public static double code(double x) {
return (2.0 * Math.asin(Math.sqrt(0.5))) + (Math.PI * 0.5);
}
def code(x): return (2.0 * math.asin(math.sqrt(0.5))) + (math.pi * 0.5)
function code(x) return Float64(Float64(2.0 * asin(sqrt(0.5))) + Float64(pi * 0.5)) end
function tmp = code(x) tmp = (2.0 * asin(sqrt(0.5))) + (pi * 0.5); end
code[x_] := N[(N[(2.0 * N[ArcSin[N[Sqrt[0.5], $MachinePrecision]], $MachinePrecision]), $MachinePrecision] + N[(Pi * 0.5), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
2 \cdot \sin^{-1} \left(\sqrt{0.5}\right) + \pi \cdot 0.5
\end{array}
Initial program 7.9%
asin-acos9.1%
add-sqr-sqrt7.4%
associate-/l*7.4%
fma-neg7.4%
div-inv7.4%
metadata-eval7.4%
Applied egg-rr7.4%
fma-neg7.4%
associate-*r/7.4%
rem-square-sqrt9.1%
*-commutative9.1%
Simplified9.1%
cancel-sign-sub-inv9.1%
div-inv9.1%
metadata-eval9.1%
metadata-eval9.1%
*-commutative9.1%
asin-acos7.9%
*-commutative7.9%
add-sqr-sqrt0.0%
sqrt-unprod3.8%
swap-sqr3.8%
pow23.8%
metadata-eval3.8%
metadata-eval3.8%
Applied egg-rr3.8%
Taylor expanded in x around 0 3.8%
Final simplification3.8%
(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 2024096
(FPCore (x)
:name "Ian Simplification"
:precision binary64
:alt
(asin x)
(- (/ PI 2.0) (* 2.0 (asin (sqrt (/ (- 1.0 x) 2.0))))))