
(FPCore (a b angle) :precision binary64 (let* ((t_0 (* (/ angle 180.0) (PI)))) (+ (pow (* a (sin t_0)) 2.0) (pow (* b (cos t_0)) 2.0))))
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{angle}{180} \cdot \mathsf{PI}\left(\right)\\
{\left(a \cdot \sin t\_0\right)}^{2} + {\left(b \cdot \cos t\_0\right)}^{2}
\end{array}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 13 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (a b angle) :precision binary64 (let* ((t_0 (* (/ angle 180.0) (PI)))) (+ (pow (* a (sin t_0)) 2.0) (pow (* b (cos t_0)) 2.0))))
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{angle}{180} \cdot \mathsf{PI}\left(\right)\\
{\left(a \cdot \sin t\_0\right)}^{2} + {\left(b \cdot \cos t\_0\right)}^{2}
\end{array}
\end{array}
angle_m = (fabs.f64 angle) (FPCore (a b angle_m) :precision binary64 (+ (pow (* a (sin (* (/ angle_m 180.0) (PI)))) 2.0) (pow (* b (cos (* 0.005555555555555556 (exp (log (* (PI) angle_m)))))) 2.0)))
\begin{array}{l}
angle_m = \left|angle\right|
\\
{\left(a \cdot \sin \left(\frac{angle\_m}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} + {\left(b \cdot \cos \left(0.005555555555555556 \cdot e^{\log \left(\mathsf{PI}\left(\right) \cdot angle\_m\right)}\right)\right)}^{2}
\end{array}
Initial program 80.4%
lift-*.f64N/A
lift-/.f64N/A
clear-numN/A
associate-*l/N/A
div-invN/A
times-fracN/A
lower-*.f64N/A
metadata-evalN/A
lower-/.f64N/A
inv-powN/A
lower-pow.f6480.4
Applied rewrites80.4%
lift-/.f64N/A
clear-numN/A
inv-powN/A
pow-to-expN/A
lower-exp.f64N/A
lower-*.f64N/A
lower-log.f64N/A
lift-pow.f64N/A
unpow-1N/A
associate-/l/N/A
lift-*.f64N/A
inv-powN/A
lower-pow.f6441.0
lift-*.f64N/A
*-commutativeN/A
lower-*.f6441.0
Applied rewrites41.0%
lift-*.f64N/A
*-commutativeN/A
mul-1-negN/A
lift-log.f64N/A
neg-logN/A
lift-pow.f64N/A
lift-*.f64N/A
*-commutativeN/A
lift-*.f64N/A
pow-flipN/A
metadata-evalN/A
unpow1N/A
lower-log.f6441.0
Applied rewrites41.0%
angle_m = (fabs.f64 angle) (FPCore (a b angle_m) :precision binary64 (+ (pow (* (sin (* (* 0.005555555555555556 (PI)) angle_m)) a) 2.0) (pow (* (cos (* -0.005555555555555556 (* (PI) angle_m))) b) 2.0)))
\begin{array}{l}
angle_m = \left|angle\right|
\\
{\left(\sin \left(\left(0.005555555555555556 \cdot \mathsf{PI}\left(\right)\right) \cdot angle\_m\right) \cdot a\right)}^{2} + {\left(\cos \left(-0.005555555555555556 \cdot \left(\mathsf{PI}\left(\right) \cdot angle\_m\right)\right) \cdot b\right)}^{2}
\end{array}
Initial program 80.4%
lift-*.f64N/A
lift-/.f64N/A
clear-numN/A
associate-*l/N/A
div-invN/A
times-fracN/A
lower-*.f64N/A
metadata-evalN/A
lower-/.f64N/A
inv-powN/A
lower-pow.f6480.4
Applied rewrites80.4%
lift-/.f64N/A
clear-numN/A
inv-powN/A
pow-to-expN/A
lower-exp.f64N/A
lower-*.f64N/A
lower-log.f64N/A
lift-pow.f64N/A
unpow-1N/A
associate-/l/N/A
lift-*.f64N/A
inv-powN/A
lower-pow.f6441.0
lift-*.f64N/A
*-commutativeN/A
lower-*.f6441.0
Applied rewrites41.0%
Applied rewrites80.5%
angle_m = (fabs.f64 angle) (FPCore (a b angle_m) :precision binary64 (+ (pow (* (sin (* (* 0.005555555555555556 angle_m) (PI))) a) 2.0) (pow (* (cos (* -0.005555555555555556 (* angle_m (PI)))) b) 2.0)))
\begin{array}{l}
angle_m = \left|angle\right|
\\
{\left(\sin \left(\left(0.005555555555555556 \cdot angle\_m\right) \cdot \mathsf{PI}\left(\right)\right) \cdot a\right)}^{2} + {\left(\cos \left(-0.005555555555555556 \cdot \left(angle\_m \cdot \mathsf{PI}\left(\right)\right)\right) \cdot b\right)}^{2}
\end{array}
Initial program 80.4%
lift-*.f64N/A
lift-/.f64N/A
clear-numN/A
associate-*l/N/A
div-invN/A
times-fracN/A
lower-*.f64N/A
metadata-evalN/A
lower-/.f64N/A
inv-powN/A
lower-pow.f6480.4
Applied rewrites80.4%
Applied rewrites80.4%
angle_m = (fabs.f64 angle) (FPCore (a b angle_m) :precision binary64 (+ (pow (* a (sin (* (/ angle_m 180.0) (PI)))) 2.0) (pow (* b 1.0) 2.0)))
\begin{array}{l}
angle_m = \left|angle\right|
\\
{\left(a \cdot \sin \left(\frac{angle\_m}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} + {\left(b \cdot 1\right)}^{2}
\end{array}
Initial program 80.4%
Taylor expanded in angle around 0
Applied rewrites79.9%
angle_m = (fabs.f64 angle)
(FPCore (a b angle_m)
:precision binary64
(if (<= a 1.55e-62)
(* (pow (cos (* -0.005555555555555556 (* (PI) angle_m))) 2.0) (* b b))
(if (<= a 8.4e+263)
(+
(* (* 3.08641975308642e-5 (* a (* (* angle_m angle_m) a))) (* (PI) (PI)))
(pow (* b (cos (* (/ angle_m 180.0) (PI)))) 2.0))
(* (pow (* (* a (PI)) angle_m) 2.0) 3.08641975308642e-5))))\begin{array}{l}
angle_m = \left|angle\right|
\\
\begin{array}{l}
\mathbf{if}\;a \leq 1.55 \cdot 10^{-62}:\\
\;\;\;\;{\cos \left(-0.005555555555555556 \cdot \left(\mathsf{PI}\left(\right) \cdot angle\_m\right)\right)}^{2} \cdot \left(b \cdot b\right)\\
\mathbf{elif}\;a \leq 8.4 \cdot 10^{+263}:\\
\;\;\;\;\left(3.08641975308642 \cdot 10^{-5} \cdot \left(a \cdot \left(\left(angle\_m \cdot angle\_m\right) \cdot a\right)\right)\right) \cdot \left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right) + {\left(b \cdot \cos \left(\frac{angle\_m}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}\\
\mathbf{else}:\\
\;\;\;\;{\left(\left(a \cdot \mathsf{PI}\left(\right)\right) \cdot angle\_m\right)}^{2} \cdot 3.08641975308642 \cdot 10^{-5}\\
\end{array}
\end{array}
if a < 1.55e-62Initial program 80.3%
Applied rewrites13.2%
Taylor expanded in a around 0
*-commutativeN/A
lower-*.f64N/A
lower-pow.f64N/A
lower-cos.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-PI.f64N/A
unpow2N/A
lower-*.f6465.4
Applied rewrites65.4%
if 1.55e-62 < a < 8.4000000000000004e263Initial program 79.0%
Taylor expanded in angle around 0
associate-*r*N/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
unpow2N/A
associate-*l*N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
lower-PI.f64N/A
lower-PI.f6471.6
Applied rewrites71.6%
if 8.4000000000000004e263 < a Initial program 98.7%
Taylor expanded in angle around 0
*-commutativeN/A
lower-fma.f64N/A
Applied rewrites20.0%
Taylor expanded in a around inf
Applied rewrites20.0%
Applied rewrites98.9%
angle_m = (fabs.f64 angle)
(FPCore (a b angle_m)
:precision binary64
(let* ((t_0 (* (PI) angle_m)))
(if (<= a 1.55e-62)
(* (pow (cos (* -0.005555555555555556 t_0)) 2.0) (* b b))
(if (<= a 2.05e+138)
(+
(*
(* (* angle_m angle_m) 3.08641975308642e-5)
(* (* (* a a) (PI)) (PI)))
(pow (* b (cos (* t_0 0.005555555555555556))) 2.0))
(* (pow (* (* a (PI)) angle_m) 2.0) 3.08641975308642e-5)))))\begin{array}{l}
angle_m = \left|angle\right|
\\
\begin{array}{l}
t_0 := \mathsf{PI}\left(\right) \cdot angle\_m\\
\mathbf{if}\;a \leq 1.55 \cdot 10^{-62}:\\
\;\;\;\;{\cos \left(-0.005555555555555556 \cdot t\_0\right)}^{2} \cdot \left(b \cdot b\right)\\
\mathbf{elif}\;a \leq 2.05 \cdot 10^{+138}:\\
\;\;\;\;\left(\left(angle\_m \cdot angle\_m\right) \cdot 3.08641975308642 \cdot 10^{-5}\right) \cdot \left(\left(\left(a \cdot a\right) \cdot \mathsf{PI}\left(\right)\right) \cdot \mathsf{PI}\left(\right)\right) + {\left(b \cdot \cos \left(t\_0 \cdot 0.005555555555555556\right)\right)}^{2}\\
\mathbf{else}:\\
\;\;\;\;{\left(\left(a \cdot \mathsf{PI}\left(\right)\right) \cdot angle\_m\right)}^{2} \cdot 3.08641975308642 \cdot 10^{-5}\\
\end{array}
\end{array}
if a < 1.55e-62Initial program 80.3%
Applied rewrites13.2%
Taylor expanded in a around 0
*-commutativeN/A
lower-*.f64N/A
lower-pow.f64N/A
lower-cos.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-PI.f64N/A
unpow2N/A
lower-*.f6465.4
Applied rewrites65.4%
if 1.55e-62 < a < 2.0499999999999999e138Initial program 68.5%
lift-*.f64N/A
lift-/.f64N/A
associate-*l/N/A
div-invN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
metadata-eval68.4
Applied rewrites68.4%
Taylor expanded in angle around 0
*-commutativeN/A
associate-*r*N/A
associate-*l*N/A
*-commutativeN/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
lower-PI.f64N/A
lower-PI.f6463.7
Applied rewrites63.7%
if 2.0499999999999999e138 < a Initial program 96.6%
Taylor expanded in angle around 0
*-commutativeN/A
lower-fma.f64N/A
Applied rewrites38.9%
Taylor expanded in a around inf
Applied rewrites59.6%
Applied rewrites86.6%
angle_m = (fabs.f64 angle)
(FPCore (a b angle_m)
:precision binary64
(if (<= a 1.55e-62)
(* (pow (cos (* -0.005555555555555556 (* (PI) angle_m))) 2.0) (* b b))
(if (<= a 2.05e+138)
(fma
(* (* (PI) (PI)) (* 3.08641975308642e-5 (* a a)))
(* angle_m angle_m)
(* b b))
(* (pow (* (* a (PI)) angle_m) 2.0) 3.08641975308642e-5))))\begin{array}{l}
angle_m = \left|angle\right|
\\
\begin{array}{l}
\mathbf{if}\;a \leq 1.55 \cdot 10^{-62}:\\
\;\;\;\;{\cos \left(-0.005555555555555556 \cdot \left(\mathsf{PI}\left(\right) \cdot angle\_m\right)\right)}^{2} \cdot \left(b \cdot b\right)\\
\mathbf{elif}\;a \leq 2.05 \cdot 10^{+138}:\\
\;\;\;\;\mathsf{fma}\left(\left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right) \cdot \left(3.08641975308642 \cdot 10^{-5} \cdot \left(a \cdot a\right)\right), angle\_m \cdot angle\_m, b \cdot b\right)\\
\mathbf{else}:\\
\;\;\;\;{\left(\left(a \cdot \mathsf{PI}\left(\right)\right) \cdot angle\_m\right)}^{2} \cdot 3.08641975308642 \cdot 10^{-5}\\
\end{array}
\end{array}
if a < 1.55e-62Initial program 80.3%
Applied rewrites13.2%
Taylor expanded in a around 0
*-commutativeN/A
lower-*.f64N/A
lower-pow.f64N/A
lower-cos.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-PI.f64N/A
unpow2N/A
lower-*.f6465.4
Applied rewrites65.4%
if 1.55e-62 < a < 2.0499999999999999e138Initial program 68.5%
Taylor expanded in angle around 0
*-commutativeN/A
lower-fma.f64N/A
Applied rewrites44.7%
Taylor expanded in a around inf
Applied rewrites63.7%
if 2.0499999999999999e138 < a Initial program 96.6%
Taylor expanded in angle around 0
*-commutativeN/A
lower-fma.f64N/A
Applied rewrites38.9%
Taylor expanded in a around inf
Applied rewrites59.6%
Applied rewrites86.6%
angle_m = (fabs.f64 angle)
(FPCore (a b angle_m)
:precision binary64
(if (<= a 1.55e-62)
(* b b)
(if (<= a 2.05e+138)
(fma
(* (* (PI) (PI)) (* 3.08641975308642e-5 (* a a)))
(* angle_m angle_m)
(* b b))
(* (pow (* (* a (PI)) angle_m) 2.0) 3.08641975308642e-5))))\begin{array}{l}
angle_m = \left|angle\right|
\\
\begin{array}{l}
\mathbf{if}\;a \leq 1.55 \cdot 10^{-62}:\\
\;\;\;\;b \cdot b\\
\mathbf{elif}\;a \leq 2.05 \cdot 10^{+138}:\\
\;\;\;\;\mathsf{fma}\left(\left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right) \cdot \left(3.08641975308642 \cdot 10^{-5} \cdot \left(a \cdot a\right)\right), angle\_m \cdot angle\_m, b \cdot b\right)\\
\mathbf{else}:\\
\;\;\;\;{\left(\left(a \cdot \mathsf{PI}\left(\right)\right) \cdot angle\_m\right)}^{2} \cdot 3.08641975308642 \cdot 10^{-5}\\
\end{array}
\end{array}
if a < 1.55e-62Initial program 80.3%
Taylor expanded in angle around 0
unpow2N/A
lower-*.f6464.7
Applied rewrites64.7%
if 1.55e-62 < a < 2.0499999999999999e138Initial program 68.5%
Taylor expanded in angle around 0
*-commutativeN/A
lower-fma.f64N/A
Applied rewrites44.7%
Taylor expanded in a around inf
Applied rewrites63.7%
if 2.0499999999999999e138 < a Initial program 96.6%
Taylor expanded in angle around 0
*-commutativeN/A
lower-fma.f64N/A
Applied rewrites38.9%
Taylor expanded in a around inf
Applied rewrites59.6%
Applied rewrites86.6%
angle_m = (fabs.f64 angle)
(FPCore (a b angle_m)
:precision binary64
(if (<= a 1.55e-62)
(* b b)
(fma
(* (* (PI) (PI)) (* 3.08641975308642e-5 (* a a)))
(* angle_m angle_m)
(* b b))))\begin{array}{l}
angle_m = \left|angle\right|
\\
\begin{array}{l}
\mathbf{if}\;a \leq 1.55 \cdot 10^{-62}:\\
\;\;\;\;b \cdot b\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right) \cdot \left(3.08641975308642 \cdot 10^{-5} \cdot \left(a \cdot a\right)\right), angle\_m \cdot angle\_m, b \cdot b\right)\\
\end{array}
\end{array}
if a < 1.55e-62Initial program 80.3%
Taylor expanded in angle around 0
unpow2N/A
lower-*.f6464.7
Applied rewrites64.7%
if 1.55e-62 < a Initial program 80.5%
Taylor expanded in angle around 0
*-commutativeN/A
lower-fma.f64N/A
Applied rewrites42.2%
Taylor expanded in a around inf
Applied rewrites61.9%
angle_m = (fabs.f64 angle)
(FPCore (a b angle_m)
:precision binary64
(if (<= a 1.55e-62)
(* b b)
(fma
(* (* (* 3.08641975308642e-5 (* a a)) (PI)) (PI))
(* angle_m angle_m)
(* b b))))\begin{array}{l}
angle_m = \left|angle\right|
\\
\begin{array}{l}
\mathbf{if}\;a \leq 1.55 \cdot 10^{-62}:\\
\;\;\;\;b \cdot b\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\left(\left(3.08641975308642 \cdot 10^{-5} \cdot \left(a \cdot a\right)\right) \cdot \mathsf{PI}\left(\right)\right) \cdot \mathsf{PI}\left(\right), angle\_m \cdot angle\_m, b \cdot b\right)\\
\end{array}
\end{array}
if a < 1.55e-62Initial program 80.3%
Taylor expanded in angle around 0
unpow2N/A
lower-*.f6464.7
Applied rewrites64.7%
if 1.55e-62 < a Initial program 80.5%
Taylor expanded in angle around 0
*-commutativeN/A
lower-fma.f64N/A
Applied rewrites42.2%
Taylor expanded in a around inf
Applied rewrites61.9%
angle_m = (fabs.f64 angle) (FPCore (a b angle_m) :precision binary64 (if (<= a 2.1e+166) (* b b) (* (* (* (* (* angle_m angle_m) 3.08641975308642e-5) (PI)) (* a (PI))) a)))
\begin{array}{l}
angle_m = \left|angle\right|
\\
\begin{array}{l}
\mathbf{if}\;a \leq 2.1 \cdot 10^{+166}:\\
\;\;\;\;b \cdot b\\
\mathbf{else}:\\
\;\;\;\;\left(\left(\left(\left(angle\_m \cdot angle\_m\right) \cdot 3.08641975308642 \cdot 10^{-5}\right) \cdot \mathsf{PI}\left(\right)\right) \cdot \left(a \cdot \mathsf{PI}\left(\right)\right)\right) \cdot a\\
\end{array}
\end{array}
if a < 2.1000000000000001e166Initial program 78.5%
Taylor expanded in angle around 0
unpow2N/A
lower-*.f6462.9
Applied rewrites62.9%
if 2.1000000000000001e166 < a Initial program 99.5%
Taylor expanded in angle around 0
*-commutativeN/A
lower-fma.f64N/A
Applied rewrites48.6%
Taylor expanded in a around inf
Applied rewrites66.0%
Applied rewrites74.8%
angle_m = (fabs.f64 angle) (FPCore (a b angle_m) :precision binary64 (if (<= a 1.35e+165) (* b b) (* (* 3.08641975308642e-5 (* a (* (* angle_m angle_m) a))) (* (PI) (PI)))))
\begin{array}{l}
angle_m = \left|angle\right|
\\
\begin{array}{l}
\mathbf{if}\;a \leq 1.35 \cdot 10^{+165}:\\
\;\;\;\;b \cdot b\\
\mathbf{else}:\\
\;\;\;\;\left(3.08641975308642 \cdot 10^{-5} \cdot \left(a \cdot \left(\left(angle\_m \cdot angle\_m\right) \cdot a\right)\right)\right) \cdot \left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right)\\
\end{array}
\end{array}
if a < 1.35e165Initial program 78.4%
Taylor expanded in angle around 0
unpow2N/A
lower-*.f6463.1
Applied rewrites63.1%
if 1.35e165 < a Initial program 99.5%
Taylor expanded in angle around 0
*-commutativeN/A
lower-fma.f64N/A
Applied rewrites46.6%
Taylor expanded in a around inf
Applied rewrites63.3%
Taylor expanded in a around 0
Applied rewrites71.7%
angle_m = (fabs.f64 angle) (FPCore (a b angle_m) :precision binary64 (* b b))
angle_m = fabs(angle);
double code(double a, double b, double angle_m) {
return b * b;
}
angle_m = abs(angle)
real(8) function code(a, b, angle_m)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: angle_m
code = b * b
end function
angle_m = Math.abs(angle);
public static double code(double a, double b, double angle_m) {
return b * b;
}
angle_m = math.fabs(angle) def code(a, b, angle_m): return b * b
angle_m = abs(angle) function code(a, b, angle_m) return Float64(b * b) end
angle_m = abs(angle); function tmp = code(a, b, angle_m) tmp = b * b; end
angle_m = N[Abs[angle], $MachinePrecision] code[a_, b_, angle$95$m_] := N[(b * b), $MachinePrecision]
\begin{array}{l}
angle_m = \left|angle\right|
\\
b \cdot b
\end{array}
Initial program 80.4%
Taylor expanded in angle around 0
unpow2N/A
lower-*.f6459.4
Applied rewrites59.4%
herbie shell --seed 2024324
(FPCore (a b angle)
:name "ab-angle->ABCF A"
:precision binary64
(+ (pow (* a (sin (* (/ angle 180.0) (PI)))) 2.0) (pow (* b (cos (* (/ angle 180.0) (PI)))) 2.0)))