
(FPCore (x) :precision binary64 (- (/ (PI) 2.0) (* 2.0 (asin (sqrt (/ (- 1.0 x) 2.0))))))
\begin{array}{l}
\\
\frac{\mathsf{PI}\left(\right)}{2} - 2 \cdot \sin^{-1} \left(\sqrt{\frac{1 - x}{2}}\right)
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 6 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x) :precision binary64 (- (/ (PI) 2.0) (* 2.0 (asin (sqrt (/ (- 1.0 x) 2.0))))))
\begin{array}{l}
\\
\frac{\mathsf{PI}\left(\right)}{2} - 2 \cdot \sin^{-1} \left(\sqrt{\frac{1 - x}{2}}\right)
\end{array}
(FPCore (x) :precision binary64 (let* ((t_0 (sqrt (/ (- 1.0 x) 2.0)))) (+ (- (/ (PI) 2.0) (- (asin t_0) (/ (PI) -2.0))) (acos t_0))))
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sqrt{\frac{1 - x}{2}}\\
\left(\frac{\mathsf{PI}\left(\right)}{2} - \left(\sin^{-1} t\_0 - \frac{\mathsf{PI}\left(\right)}{-2}\right)\right) + \cos^{-1} t\_0
\end{array}
\end{array}
Initial program 7.2%
lift--.f64N/A
lift-*.f64N/A
count-2-revN/A
lift-asin.f64N/A
asin-acosN/A
lift-PI.f64N/A
lift-/.f64N/A
associate-+r-N/A
associate--r-N/A
lower-+.f64N/A
Applied rewrites9.0%
(FPCore (x) :precision binary64 (fma -2.0 (- (/ (PI) 2.0) (acos (sqrt (* 0.5 (- 1.0 x))))) (* (PI) 0.5)))
\begin{array}{l}
\\
\mathsf{fma}\left(-2, \frac{\mathsf{PI}\left(\right)}{2} - \cos^{-1} \left(\sqrt{0.5 \cdot \left(1 - x\right)}\right), \mathsf{PI}\left(\right) \cdot 0.5\right)
\end{array}
Initial program 7.2%
Taylor expanded in x around 0
fp-cancel-sub-sign-invN/A
+-commutativeN/A
lower-fma.f64N/A
metadata-evalN/A
lower-asin.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-sqrt.f64N/A
lower--.f64N/A
lower-sqrt.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-PI.f647.1
Applied rewrites7.1%
Applied rewrites9.0%
(FPCore (x) :precision binary64 (fma (asin (sqrt (* 0.5 (- 1.0 x)))) -2.0 (* 0.5 (PI))))
\begin{array}{l}
\\
\mathsf{fma}\left(\sin^{-1} \left(\sqrt{0.5 \cdot \left(1 - x\right)}\right), -2, 0.5 \cdot \mathsf{PI}\left(\right)\right)
\end{array}
Initial program 7.2%
Taylor expanded in x around 0
fp-cancel-sub-sign-invN/A
+-commutativeN/A
lower-fma.f64N/A
metadata-evalN/A
lower-asin.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-sqrt.f64N/A
lower--.f64N/A
lower-sqrt.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-PI.f647.1
Applied rewrites7.1%
Applied rewrites7.2%
(FPCore (x) :precision binary64 (fma (asin (sqrt 0.5)) 2.0 (/ (PI) -2.0)))
\begin{array}{l}
\\
\mathsf{fma}\left(\sin^{-1} \left(\sqrt{0.5}\right), 2, \frac{\mathsf{PI}\left(\right)}{-2}\right)
\end{array}
Initial program 7.2%
Applied rewrites3.8%
Taylor expanded in x around 0
Applied rewrites4.3%
(FPCore (x) :precision binary64 (fma -2.0 (asin (sqrt 0.5)) (* (PI) 0.5)))
\begin{array}{l}
\\
\mathsf{fma}\left(-2, \sin^{-1} \left(\sqrt{0.5}\right), \mathsf{PI}\left(\right) \cdot 0.5\right)
\end{array}
Initial program 7.2%
Taylor expanded in x around 0
fp-cancel-sub-sign-invN/A
+-commutativeN/A
lower-fma.f64N/A
metadata-evalN/A
lower-asin.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-sqrt.f64N/A
lower--.f64N/A
lower-sqrt.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-PI.f647.1
Applied rewrites7.1%
Taylor expanded in x around 0
Applied rewrites3.8%
Taylor expanded in x around 0
Applied rewrites3.8%
(FPCore (x) :precision binary64 (/ 0.0 0.0))
double code(double x) {
return 0.0 / 0.0;
}
real(8) function code(x)
real(8), intent (in) :: x
code = 0.0d0 / 0.0d0
end function
public static double code(double x) {
return 0.0 / 0.0;
}
def code(x): return 0.0 / 0.0
function code(x) return Float64(0.0 / 0.0) end
function tmp = code(x) tmp = 0.0 / 0.0; end
code[x_] := N[(0.0 / 0.0), $MachinePrecision]
\begin{array}{l}
\\
\frac{0}{0}
\end{array}
Initial program 7.2%
lift--.f64N/A
lift-*.f64N/A
fp-cancel-sub-sign-invN/A
+-commutativeN/A
Applied rewrites0.0%
(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 2024337
(FPCore (x)
:name "Ian Simplification"
:precision binary64
:alt
(! :herbie-platform default (asin x))
(- (/ (PI) 2.0) (* 2.0 (asin (sqrt (/ (- 1.0 x) 2.0))))))