
(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}
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 (let* ((t_0 (acos (sqrt (fma x -0.5 0.5)))) (t_1 (- (/ PI 2.0) t_0))) (/ (- (* t_0 t_0) (* t_1 t_1)) (+ t_0 t_1))))
double code(double x) {
double t_0 = acos(sqrt(fma(x, -0.5, 0.5)));
double t_1 = (((double) M_PI) / 2.0) - t_0;
return ((t_0 * t_0) - (t_1 * t_1)) / (t_0 + t_1);
}
function code(x) t_0 = acos(sqrt(fma(x, -0.5, 0.5))) t_1 = Float64(Float64(pi / 2.0) - t_0) return Float64(Float64(Float64(t_0 * t_0) - Float64(t_1 * t_1)) / Float64(t_0 + t_1)) end
code[x_] := Block[{t$95$0 = N[ArcCos[N[Sqrt[N[(x * -0.5 + 0.5), $MachinePrecision]], $MachinePrecision]], $MachinePrecision]}, Block[{t$95$1 = N[(N[(Pi / 2.0), $MachinePrecision] - t$95$0), $MachinePrecision]}, N[(N[(N[(t$95$0 * t$95$0), $MachinePrecision] - N[(t$95$1 * t$95$1), $MachinePrecision]), $MachinePrecision] / N[(t$95$0 + t$95$1), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \cos^{-1} \left(\sqrt{\mathsf{fma}\left(x, -0.5, 0.5\right)}\right)\\
t_1 := \frac{\pi}{2} - t\_0\\
\frac{t\_0 \cdot t\_0 - t\_1 \cdot t\_1}{t\_0 + t\_1}
\end{array}
\end{array}
Initial program 6.8%
lift--.f64N/A
lift-*.f64N/A
lift-asin.f64N/A
lift-sqrt.f64N/A
lift--.f64N/A
lift-/.f64N/A
count-2-revN/A
associate--r+N/A
lift-PI.f64N/A
lift-/.f64N/A
acos-asinN/A
lower--.f64N/A
Applied rewrites6.8%
lift--.f64N/A
lift-acos.f64N/A
lift-fma.f64N/A
lift-sqrt.f64N/A
lift-asin.f64N/A
lift-fma.f64N/A
lift-sqrt.f64N/A
flip--N/A
lower-/.f64N/A
Applied rewrites6.8%
lift-asin.f64N/A
lift-fma.f64N/A
lift-sqrt.f64N/A
*-commutativeN/A
flip-+N/A
flip-+N/A
+-commutativeN/A
metadata-evalN/A
fp-cancel-sub-sign-invN/A
asin-acosN/A
lower--.f64N/A
lift-/.f64N/A
lift-PI.f64N/A
fp-cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
*-commutativeN/A
lift-sqrt.f64N/A
lift-fma.f64N/A
Applied rewrites6.8%
lift-asin.f64N/A
lift-fma.f64N/A
lift-sqrt.f64N/A
*-commutativeN/A
flip-+N/A
flip-+N/A
+-commutativeN/A
metadata-evalN/A
fp-cancel-sub-sign-invN/A
asin-acosN/A
lower--.f64N/A
lift-/.f64N/A
lift-PI.f64N/A
fp-cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
*-commutativeN/A
lift-sqrt.f64N/A
lift-fma.f64N/A
Applied rewrites8.2%
lift-asin.f64N/A
lift-fma.f64N/A
lift-sqrt.f64N/A
*-commutativeN/A
flip-+N/A
flip-+N/A
+-commutativeN/A
metadata-evalN/A
fp-cancel-sub-sign-invN/A
asin-acosN/A
lower--.f64N/A
lift-/.f64N/A
lift-PI.f64N/A
fp-cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
*-commutativeN/A
lift-sqrt.f64N/A
lift-fma.f64N/A
Applied rewrites8.2%
(FPCore (x) :precision binary64 (let* ((t_0 (- (* 0.5 PI) (acos (sqrt (fma x -0.5 0.5)))))) (/ (fma (* PI PI) 0.25 (* -4.0 (pow t_0 2.0))) (fma 0.5 PI (* t_0 2.0)))))
double code(double x) {
double t_0 = (0.5 * ((double) M_PI)) - acos(sqrt(fma(x, -0.5, 0.5)));
return fma((((double) M_PI) * ((double) M_PI)), 0.25, (-4.0 * pow(t_0, 2.0))) / fma(0.5, ((double) M_PI), (t_0 * 2.0));
}
function code(x) t_0 = Float64(Float64(0.5 * pi) - acos(sqrt(fma(x, -0.5, 0.5)))) return Float64(fma(Float64(pi * pi), 0.25, Float64(-4.0 * (t_0 ^ 2.0))) / fma(0.5, pi, Float64(t_0 * 2.0))) end
code[x_] := Block[{t$95$0 = N[(N[(0.5 * Pi), $MachinePrecision] - N[ArcCos[N[Sqrt[N[(x * -0.5 + 0.5), $MachinePrecision]], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, N[(N[(N[(Pi * Pi), $MachinePrecision] * 0.25 + N[(-4.0 * N[Power[t$95$0, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(0.5 * Pi + N[(t$95$0 * 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := 0.5 \cdot \pi - \cos^{-1} \left(\sqrt{\mathsf{fma}\left(x, -0.5, 0.5\right)}\right)\\
\frac{\mathsf{fma}\left(\pi \cdot \pi, 0.25, -4 \cdot {t\_0}^{2}\right)}{\mathsf{fma}\left(0.5, \pi, t\_0 \cdot 2\right)}
\end{array}
\end{array}
Initial program 6.8%
lift-asin.f64N/A
lift-sqrt.f64N/A
lift--.f64N/A
lift-/.f64N/A
asin-acosN/A
lift-/.f64N/A
lift-PI.f64N/A
lower--.f64N/A
lower-acos.f64N/A
lift-/.f64N/A
lift--.f64N/A
lift-sqrt.f648.2
lift--.f64N/A
lift-/.f64N/A
div-subN/A
metadata-evalN/A
*-lft-identityN/A
associate-*l/N/A
metadata-evalN/A
metadata-evalN/A
fp-cancel-sign-sub-invN/A
+-commutativeN/A
lower-fma.f648.2
Applied rewrites8.2%
Taylor expanded in x around 0
fp-cancel-sub-sign-invN/A
lift-PI.f64N/A
metadata-evalN/A
fp-cancel-sub-sign-invN/A
lower-fma.f64N/A
metadata-evalN/A
lower-*.f64N/A
Applied rewrites8.2%
Applied rewrites8.2%
Applied rewrites8.2%
(FPCore (x) :precision binary64 (fma 0.5 PI (* -2.0 (- (* 0.5 PI) (acos (sqrt (fma x -0.5 0.5)))))))
double code(double x) {
return fma(0.5, ((double) M_PI), (-2.0 * ((0.5 * ((double) M_PI)) - acos(sqrt(fma(x, -0.5, 0.5))))));
}
function code(x) return fma(0.5, pi, Float64(-2.0 * Float64(Float64(0.5 * pi) - acos(sqrt(fma(x, -0.5, 0.5)))))) end
code[x_] := N[(0.5 * Pi + N[(-2.0 * N[(N[(0.5 * Pi), $MachinePrecision] - N[ArcCos[N[Sqrt[N[(x * -0.5 + 0.5), $MachinePrecision]], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\mathsf{fma}\left(0.5, \pi, -2 \cdot \left(0.5 \cdot \pi - \cos^{-1} \left(\sqrt{\mathsf{fma}\left(x, -0.5, 0.5\right)}\right)\right)\right)
\end{array}
Initial program 6.8%
lift-asin.f64N/A
lift-sqrt.f64N/A
lift--.f64N/A
lift-/.f64N/A
asin-acosN/A
lift-/.f64N/A
lift-PI.f64N/A
lower--.f64N/A
lower-acos.f64N/A
lift-/.f64N/A
lift--.f64N/A
lift-sqrt.f648.2
lift--.f64N/A
lift-/.f64N/A
div-subN/A
metadata-evalN/A
*-lft-identityN/A
associate-*l/N/A
metadata-evalN/A
metadata-evalN/A
fp-cancel-sign-sub-invN/A
+-commutativeN/A
lower-fma.f648.2
Applied rewrites8.2%
Taylor expanded in x around 0
fp-cancel-sub-sign-invN/A
lift-PI.f64N/A
metadata-evalN/A
fp-cancel-sub-sign-invN/A
lower-fma.f64N/A
metadata-evalN/A
lower-*.f64N/A
Applied rewrites8.2%
(FPCore (x) :precision binary64 (fma -2.0 (asin (sqrt (fma -0.5 x 0.5))) (* 0.5 PI)))
double code(double x) {
return fma(-2.0, asin(sqrt(fma(-0.5, x, 0.5))), (0.5 * ((double) M_PI)));
}
function code(x) return fma(-2.0, asin(sqrt(fma(-0.5, x, 0.5))), Float64(0.5 * pi)) end
code[x_] := N[(-2.0 * N[ArcSin[N[Sqrt[N[(-0.5 * x + 0.5), $MachinePrecision]], $MachinePrecision]], $MachinePrecision] + N[(0.5 * Pi), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\mathsf{fma}\left(-2, \sin^{-1} \left(\sqrt{\mathsf{fma}\left(-0.5, x, 0.5\right)}\right), 0.5 \cdot \pi\right)
\end{array}
Initial program 6.8%
Taylor expanded in x around 0
fp-cancel-sub-sign-invN/A
+-commutativeN/A
distribute-lft-out--N/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
associate-*l/N/A
*-lft-identityN/A
div-subN/A
lower-fma.f64N/A
Applied rewrites6.8%
(FPCore (x) :precision binary64 (fma -2.0 (asin (sqrt 0.5)) (* 0.5 PI)))
double code(double x) {
return fma(-2.0, asin(sqrt(0.5)), (0.5 * ((double) M_PI)));
}
function code(x) return fma(-2.0, asin(sqrt(0.5)), Float64(0.5 * pi)) end
code[x_] := N[(-2.0 * N[ArcSin[N[Sqrt[0.5], $MachinePrecision]], $MachinePrecision] + N[(0.5 * Pi), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\mathsf{fma}\left(-2, \sin^{-1} \left(\sqrt{0.5}\right), 0.5 \cdot \pi\right)
\end{array}
Initial program 6.8%
Taylor expanded in x around 0
fp-cancel-sub-sign-invN/A
+-commutativeN/A
distribute-lft-out--N/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
associate-*l/N/A
*-lft-identityN/A
div-subN/A
lower-fma.f64N/A
Applied rewrites6.8%
Taylor expanded in x around 0
Applied rewrites4.1%
herbie shell --seed 2025135
(FPCore (x)
:name "Ian Simplification"
:precision binary64
(- (/ PI 2.0) (* 2.0 (asin (sqrt (/ (- 1.0 x) 2.0))))))