ab-angle->ABCF C

Percentage Accurate: 79.3% → 79.3%

Time: 9.3s

Alternatives: 9

Speedup: N/A×

• Could not converge on a ground truth

  1. a = 7.262699913972246e-84
  2. b = 1.05994409102362e-10
  3. angle = 4.0274830554216396e+229

• 36.7% of points produce a very large (infinite) output. You may want to add a precondition.

Specification

Math

\[\begin{array}{l} \\ \begin{array}{l} t_0 := \mathsf{PI}\left(\right) \cdot \frac{angle}{180}\\ {\left(a \cdot \cos t\_0\right)}^{2} + {\left(b \cdot \sin t\_0\right)}^{2} \end{array} \end{array} \]

FPCore

(FPCore (a b angle)
 :precision binary64
 (let* ((t_0 (* (PI) (/ angle 180.0))))
   (+ (pow (* a (cos t_0)) 2.0) (pow (* b (sin t_0)) 2.0))))

Initial Program: 79.3% accurate, 1.0× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} t_0 := \mathsf{PI}\left(\right) \cdot \frac{angle}{180}\\ {\left(a \cdot \cos t\_0\right)}^{2} + {\left(b \cdot \sin t\_0\right)}^{2} \end{array} \end{array} \]

FPCore

(FPCore (a b angle)
 :precision binary64
 (let* ((t_0 (* (PI) (/ angle 180.0))))
   (+ (pow (* a (cos t_0)) 2.0) (pow (* b (sin t_0)) 2.0))))

Alternative 1: 79.3% accurate, 2.0× speedup

Math

\[\begin{array}{l} \\ a \cdot a + {\left(b \cdot \sin \left(\left(0.005555555555555556 \cdot angle\right) \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} \end{array} \]

FPCore

(FPCore (a b angle)
 :precision binary64
 (+ (* a a) (pow (* b (sin (* (* 0.005555555555555556 angle) (PI)))) 2.0)))

Derivation

1. Initial program 85.2%

   \[{\left(a \cdot \cos \left(\mathsf{PI}\left(\right) \cdot \frac{angle}{180}\right)\right)}^{2} + {\left(b \cdot \sin \left(\mathsf{PI}\left(\right) \cdot \frac{angle}{180}\right)\right)}^{2} \]
2. Add Preprocessing

3. Taylor expanded in angle around 0

   \[\leadsto {\left(a \cdot \color{blue}{1}\right)}^{2} + {\left(b \cdot \sin \left(\mathsf{PI}\left(\right) \cdot \frac{angle}{180}\right)\right)}^{2} \]
4. Step-by-step derivation

5. Applied rewrites 85.3%

   \[\leadsto {\left(a \cdot \color{blue}{1}\right)}^{2} + {\left(b \cdot \sin \left(\mathsf{PI}\left(\right) \cdot \frac{angle}{180}\right)\right)}^{2} \]
6. Step-by-step derivation

   1. lift-*.f64 N/A

      \[\leadsto {\left(a \cdot 1\right)}^{2} + {\left(b \cdot \sin \color{blue}{\left(\mathsf{PI}\left(\right) \cdot \frac{angle}{180}\right)}\right)}^{2} \]

   2. lift-/.f64 N/A

      \[\leadsto {\left(a \cdot 1\right)}^{2} + {\left(b \cdot \sin \left(\mathsf{PI}\left(\right) \cdot \color{blue}{\frac{angle}{180}}\right)\right)}^{2} \]

   3. div-inv N/A

      \[\leadsto {\left(a \cdot 1\right)}^{2} + {\left(b \cdot \sin \left(\mathsf{PI}\left(\right) \cdot \color{blue}{\left(angle \cdot \frac{1}{180}\right)}\right)\right)}^{2} \]

   4. metadata-eval N/A

      \[\leadsto {\left(a \cdot 1\right)}^{2} + {\left(b \cdot \sin \left(\mathsf{PI}\left(\right) \cdot \left(angle \cdot \color{blue}{\frac{1}{180}}\right)\right)\right)}^{2} \]

   5. *-commutative N/A

      \[\leadsto {\left(a \cdot 1\right)}^{2} + {\left(b \cdot \sin \left(\mathsf{PI}\left(\right) \cdot \color{blue}{\left(\frac{1}{180} \cdot angle\right)}\right)\right)}^{2} \]

   6. lift-*.f64 N/A

      \[\leadsto {\left(a \cdot 1\right)}^{2} + {\left(b \cdot \sin \left(\mathsf{PI}\left(\right) \cdot \color{blue}{\left(\frac{1}{180} \cdot angle\right)}\right)\right)}^{2} \]

   7. *-commutative N/A

      \[\leadsto {\left(a \cdot 1\right)}^{2} + {\left(b \cdot \sin \color{blue}{\left(\left(\frac{1}{180} \cdot angle\right) \cdot \mathsf{PI}\left(\right)\right)}\right)}^{2} \]

   8. lift-*.f64 85.4

      \[\leadsto {\left(a \cdot 1\right)}^{2} + {\left(b \cdot \sin \color{blue}{\left(\left(0.005555555555555556 \cdot angle\right) \cdot \mathsf{PI}\left(\right)\right)}\right)}^{2} \]

7. Applied rewrites 85.4%

   \[\leadsto {\left(a \cdot 1\right)}^{2} + {\left(b \cdot \sin \color{blue}{\left(\left(0.005555555555555556 \cdot angle\right) \cdot \mathsf{PI}\left(\right)\right)}\right)}^{2} \]
8. Step-by-step derivation

   1. lift-pow.f64 N/A

      \[\leadsto \color{blue}{{\left(a \cdot 1\right)}^{2}} + {\left(b \cdot \sin \left(\left(\frac{1}{180} \cdot angle\right) \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} \]

   2. unpow2 N/A

      \[\leadsto \color{blue}{\left(a \cdot 1\right) \cdot \left(a \cdot 1\right)} + {\left(b \cdot \sin \left(\left(\frac{1}{180} \cdot angle\right) \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} \]

   3. lift-*.f64 N/A

      \[\leadsto \left(a \cdot 1\right) \cdot \color{blue}{\left(a \cdot 1\right)} + {\left(b \cdot \sin \left(\left(\frac{1}{180} \cdot angle\right) \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} \]

   4. *-commutative N/A

      \[\leadsto \left(a \cdot 1\right) \cdot \color{blue}{\left(1 \cdot a\right)} + {\left(b \cdot \sin \left(\left(\frac{1}{180} \cdot angle\right) \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} \]

   5. associate-*r* N/A

      \[\leadsto \color{blue}{\left(\left(a \cdot 1\right) \cdot 1\right) \cdot a} + {\left(b \cdot \sin \left(\left(\frac{1}{180} \cdot angle\right) \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} \]

   6. lower-*.f64 N/A

      \[\leadsto \color{blue}{\left(\left(a \cdot 1\right) \cdot 1\right) \cdot a} + {\left(b \cdot \sin \left(\left(\frac{1}{180} \cdot angle\right) \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} \]

   7. lower-*.f64 85.4

      \[\leadsto \color{blue}{\left(\left(a \cdot 1\right) \cdot 1\right)} \cdot a + {\left(b \cdot \sin \left(\left(0.005555555555555556 \cdot angle\right) \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} \]

   8. lift-*.f64 N/A

      \[\leadsto \left(\color{blue}{\left(a \cdot 1\right)} \cdot 1\right) \cdot a + {\left(b \cdot \sin \left(\left(\frac{1}{180} \cdot angle\right) \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} \]

   9. *-commutative N/A

      \[\leadsto \left(\color{blue}{\left(1 \cdot a\right)} \cdot 1\right) \cdot a + {\left(b \cdot \sin \left(\left(\frac{1}{180} \cdot angle\right) \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} \]

   10. lower-*.f64 85.4

       \[\leadsto \left(\color{blue}{\left(1 \cdot a\right)} \cdot 1\right) \cdot a + {\left(b \cdot \sin \left(\left(0.005555555555555556 \cdot angle\right) \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} \]

9. Applied rewrites 85.4%

   \[\leadsto \color{blue}{\left(\left(1 \cdot a\right) \cdot 1\right) \cdot a} + {\left(b \cdot \sin \left(\left(0.005555555555555556 \cdot angle\right) \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} \]

10. Taylor expanded in angle around 0

    \[\leadsto \color{blue}{{a}^{2}} + {\left(b \cdot \sin \left(\left(\frac{1}{180} \cdot angle\right) \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} \]
11. Step-by-step derivation

    1. unpow2 N/A

       \[\leadsto \color{blue}{a \cdot a} + {\left(b \cdot \sin \left(\left(\frac{1}{180} \cdot angle\right) \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} \]

    2. lower-*.f64 85.4

       \[\leadsto \color{blue}{a \cdot a} + {\left(b \cdot \sin \left(\left(0.005555555555555556 \cdot angle\right) \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} \]

12. Applied rewrites 85.4%

    \[\leadsto \color{blue}{a \cdot a} + {\left(b \cdot \sin \left(\left(0.005555555555555556 \cdot angle\right) \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} \]
13. Add Preprocessing

Alternative 2: 66.8% accurate, 2.8× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;b \leq 3.4 \cdot 10^{-45}:\\ \;\;\;\;a \cdot a\\ \mathbf{else}:\\ \;\;\;\;\left(\left(1 \cdot a\right) \cdot 1\right) \cdot a + {\left(b \cdot \left(\left(\mathsf{PI}\left(\right) \cdot \mathsf{fma}\left(-2.8577960676726107 \cdot 10^{-8} \cdot \left(angle \cdot angle\right), \mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right), 0.005555555555555556\right)\right) \cdot angle\right)\right)}^{2}\\ \end{array} \end{array} \]

FPCore

(FPCore (a b angle)
 :precision binary64
 (if (<= b 3.4e-45)
   (* a a)
   (+
    (* (* (* 1.0 a) 1.0) a)
    (pow
     (*
      b
      (*
       (*
        (PI)
        (fma
         (* -2.8577960676726107e-8 (* angle angle))
         (* (PI) (PI))
         0.005555555555555556))
       angle))
     2.0))))

Derivation

1. Split input into 2 regimes

2. if b < 3.40000000000000004e-45

   1. Initial program 83.3%

      \[{\left(a \cdot \cos \left(\mathsf{PI}\left(\right) \cdot \frac{angle}{180}\right)\right)}^{2} + {\left(b \cdot \sin \left(\mathsf{PI}\left(\right) \cdot \frac{angle}{180}\right)\right)}^{2} \]
   2. Add Preprocessing

   3. Taylor expanded in angle around 0

      \[\leadsto \color{blue}{{a}^{2}} \]
   4. Step-by-step derivation

      1. unpow2 N/A

         \[\leadsto \color{blue}{a \cdot a} \]

      2. lower-*.f64 65.9

         \[\leadsto \color{blue}{a \cdot a} \]

   5. Applied rewrites 65.9%

      \[\leadsto \color{blue}{a \cdot a} \]

   if 3.40000000000000004e-45 < b

   1. Initial program 89.0%

      \[{\left(a \cdot \cos \left(\mathsf{PI}\left(\right) \cdot \frac{angle}{180}\right)\right)}^{2} + {\left(b \cdot \sin \left(\mathsf{PI}\left(\right) \cdot \frac{angle}{180}\right)\right)}^{2} \]
   2. Add Preprocessing

   3. Taylor expanded in angle around 0

      \[\leadsto {\left(a \cdot \color{blue}{1}\right)}^{2} + {\left(b \cdot \sin \left(\mathsf{PI}\left(\right) \cdot \frac{angle}{180}\right)\right)}^{2} \]
   4. Step-by-step derivation

   5. Applied rewrites 88.9%

      \[\leadsto {\left(a \cdot \color{blue}{1}\right)}^{2} + {\left(b \cdot \sin \left(\mathsf{PI}\left(\right) \cdot \frac{angle}{180}\right)\right)}^{2} \]
   6. Step-by-step derivation

      1. lift-*.f64 N/A

         \[\leadsto {\left(a \cdot 1\right)}^{2} + {\left(b \cdot \sin \color{blue}{\left(\mathsf{PI}\left(\right) \cdot \frac{angle}{180}\right)}\right)}^{2} \]

      2. lift-/.f64 N/A

         \[\leadsto {\left(a \cdot 1\right)}^{2} + {\left(b \cdot \sin \left(\mathsf{PI}\left(\right) \cdot \color{blue}{\frac{angle}{180}}\right)\right)}^{2} \]

      3. div-inv N/A

         \[\leadsto {\left(a \cdot 1\right)}^{2} + {\left(b \cdot \sin \left(\mathsf{PI}\left(\right) \cdot \color{blue}{\left(angle \cdot \frac{1}{180}\right)}\right)\right)}^{2} \]

      4. metadata-eval N/A

         \[\leadsto {\left(a \cdot 1\right)}^{2} + {\left(b \cdot \sin \left(\mathsf{PI}\left(\right) \cdot \left(angle \cdot \color{blue}{\frac{1}{180}}\right)\right)\right)}^{2} \]

      5. *-commutative N/A

         \[\leadsto {\left(a \cdot 1\right)}^{2} + {\left(b \cdot \sin \left(\mathsf{PI}\left(\right) \cdot \color{blue}{\left(\frac{1}{180} \cdot angle\right)}\right)\right)}^{2} \]

      6. lift-*.f64 N/A

         \[\leadsto {\left(a \cdot 1\right)}^{2} + {\left(b \cdot \sin \left(\mathsf{PI}\left(\right) \cdot \color{blue}{\left(\frac{1}{180} \cdot angle\right)}\right)\right)}^{2} \]

      7. *-commutative N/A

         \[\leadsto {\left(a \cdot 1\right)}^{2} + {\left(b \cdot \sin \color{blue}{\left(\left(\frac{1}{180} \cdot angle\right) \cdot \mathsf{PI}\left(\right)\right)}\right)}^{2} \]

      8. lift-*.f64 89.0

         \[\leadsto {\left(a \cdot 1\right)}^{2} + {\left(b \cdot \sin \color{blue}{\left(\left(0.005555555555555556 \cdot angle\right) \cdot \mathsf{PI}\left(\right)\right)}\right)}^{2} \]

   7. Applied rewrites 89.0%

      \[\leadsto {\left(a \cdot 1\right)}^{2} + {\left(b \cdot \sin \color{blue}{\left(\left(0.005555555555555556 \cdot angle\right) \cdot \mathsf{PI}\left(\right)\right)}\right)}^{2} \]
   8. Step-by-step derivation

      1. lift-pow.f64 N/A

         \[\leadsto \color{blue}{{\left(a \cdot 1\right)}^{2}} + {\left(b \cdot \sin \left(\left(\frac{1}{180} \cdot angle\right) \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} \]

      2. unpow2 N/A

         \[\leadsto \color{blue}{\left(a \cdot 1\right) \cdot \left(a \cdot 1\right)} + {\left(b \cdot \sin \left(\left(\frac{1}{180} \cdot angle\right) \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} \]

      3. lift-*.f64 N/A

         \[\leadsto \left(a \cdot 1\right) \cdot \color{blue}{\left(a \cdot 1\right)} + {\left(b \cdot \sin \left(\left(\frac{1}{180} \cdot angle\right) \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} \]

      4. *-commutative N/A

         \[\leadsto \left(a \cdot 1\right) \cdot \color{blue}{\left(1 \cdot a\right)} + {\left(b \cdot \sin \left(\left(\frac{1}{180} \cdot angle\right) \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} \]

      5. associate-*r* N/A

         \[\leadsto \color{blue}{\left(\left(a \cdot 1\right) \cdot 1\right) \cdot a} + {\left(b \cdot \sin \left(\left(\frac{1}{180} \cdot angle\right) \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} \]

      6. lower-*.f64 N/A

         \[\leadsto \color{blue}{\left(\left(a \cdot 1\right) \cdot 1\right) \cdot a} + {\left(b \cdot \sin \left(\left(\frac{1}{180} \cdot angle\right) \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} \]

      7. lower-*.f64 89.0

         \[\leadsto \color{blue}{\left(\left(a \cdot 1\right) \cdot 1\right)} \cdot a + {\left(b \cdot \sin \left(\left(0.005555555555555556 \cdot angle\right) \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} \]

      8. lift-*.f64 N/A

         \[\leadsto \left(\color{blue}{\left(a \cdot 1\right)} \cdot 1\right) \cdot a + {\left(b \cdot \sin \left(\left(\frac{1}{180} \cdot angle\right) \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} \]

      9. *-commutative N/A

         \[\leadsto \left(\color{blue}{\left(1 \cdot a\right)} \cdot 1\right) \cdot a + {\left(b \cdot \sin \left(\left(\frac{1}{180} \cdot angle\right) \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} \]

      10. lower-*.f64 89.0

          \[\leadsto \left(\color{blue}{\left(1 \cdot a\right)} \cdot 1\right) \cdot a + {\left(b \cdot \sin \left(\left(0.005555555555555556 \cdot angle\right) \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} \]

   9. Applied rewrites 89.0%

      \[\leadsto \color{blue}{\left(\left(1 \cdot a\right) \cdot 1\right) \cdot a} + {\left(b \cdot \sin \left(\left(0.005555555555555556 \cdot angle\right) \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} \]

   10. Taylor expanded in angle around 0

       \[\leadsto \left(\left(1 \cdot a\right) \cdot 1\right) \cdot a + {\left(b \cdot \color{blue}{\left(angle \cdot \left(\frac{-1}{34992000} \cdot \left({angle}^{2} \cdot {\mathsf{PI}\left(\right)}^{3}\right) + \frac{1}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}\right)}^{2} \]
   11. Step-by-step derivation

       1. *-commutative N/A

          \[\leadsto \left(\left(1 \cdot a\right) \cdot 1\right) \cdot a + {\left(b \cdot \color{blue}{\left(\left(\frac{-1}{34992000} \cdot \left({angle}^{2} \cdot {\mathsf{PI}\left(\right)}^{3}\right) + \frac{1}{180} \cdot \mathsf{PI}\left(\right)\right) \cdot angle\right)}\right)}^{2} \]

       2. +-commutative N/A

          \[\leadsto \left(\left(1 \cdot a\right) \cdot 1\right) \cdot a + {\left(b \cdot \left(\color{blue}{\left(\frac{1}{180} \cdot \mathsf{PI}\left(\right) + \frac{-1}{34992000} \cdot \left({angle}^{2} \cdot {\mathsf{PI}\left(\right)}^{3}\right)\right)} \cdot angle\right)\right)}^{2} \]

       3. *-commutative N/A

          \[\leadsto \left(\left(1 \cdot a\right) \cdot 1\right) \cdot a + {\left(b \cdot \left(\left(\frac{1}{180} \cdot \mathsf{PI}\left(\right) + \color{blue}{\left({angle}^{2} \cdot {\mathsf{PI}\left(\right)}^{3}\right) \cdot \frac{-1}{34992000}}\right) \cdot angle\right)\right)}^{2} \]

       4. associate-*r* N/A

          \[\leadsto \left(\left(1 \cdot a\right) \cdot 1\right) \cdot a + {\left(b \cdot \left(\left(\frac{1}{180} \cdot \mathsf{PI}\left(\right) + \color{blue}{{angle}^{2} \cdot \left({\mathsf{PI}\left(\right)}^{3} \cdot \frac{-1}{34992000}\right)}\right) \cdot angle\right)\right)}^{2} \]

       5. *-commutative N/A

          \[\leadsto \left(\left(1 \cdot a\right) \cdot 1\right) \cdot a + {\left(b \cdot \left(\left(\frac{1}{180} \cdot \mathsf{PI}\left(\right) + {angle}^{2} \cdot \color{blue}{\left(\frac{-1}{34992000} \cdot {\mathsf{PI}\left(\right)}^{3}\right)}\right) \cdot angle\right)\right)}^{2} \]

       6. lower-*.f64 N/A

          \[\leadsto \left(\left(1 \cdot a\right) \cdot 1\right) \cdot a + {\left(b \cdot \color{blue}{\left(\left(\frac{1}{180} \cdot \mathsf{PI}\left(\right) + {angle}^{2} \cdot \left(\frac{-1}{34992000} \cdot {\mathsf{PI}\left(\right)}^{3}\right)\right) \cdot angle\right)}\right)}^{2} \]

   12. Applied rewrites 87.5%

       \[\leadsto \left(\left(1 \cdot a\right) \cdot 1\right) \cdot a + {\left(b \cdot \color{blue}{\left(\left(\mathsf{PI}\left(\right) \cdot \mathsf{fma}\left(-2.8577960676726107 \cdot 10^{-8} \cdot \left(angle \cdot angle\right), \mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right), 0.005555555555555556\right)\right) \cdot angle\right)}\right)}^{2} \]
3. Recombined 2 regimes into one program.
4. Add Preprocessing

Alternative 3: 64.4% accurate, 3.5× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;b \leq 3.4 \cdot 10^{-45}:\\ \;\;\;\;a \cdot a\\ \mathbf{elif}\;b \leq 1.9 \cdot 10^{+158}:\\ \;\;\;\;\mathsf{fma}\left(\left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right) \cdot \left(3.08641975308642 \cdot 10^{-5} \cdot \left(b \cdot b\right)\right), angle \cdot angle, a \cdot a\right)\\ \mathbf{else}:\\ \;\;\;\;{\left(\left(\mathsf{PI}\left(\right) \cdot b\right) \cdot angle\right)}^{2} \cdot 3.08641975308642 \cdot 10^{-5}\\ \end{array} \end{array} \]

FPCore

(FPCore (a b angle)
 :precision binary64
 (if (<= b 3.4e-45)
   (* a a)
   (if (<= b 1.9e+158)
     (fma
      (* (* (PI) (PI)) (* 3.08641975308642e-5 (* b b)))
      (* angle angle)
      (* a a))
     (* (pow (* (* (PI) b) angle) 2.0) 3.08641975308642e-5))))

Derivation

1. Split input into 3 regimes

2. if b < 3.40000000000000004e-45

   1. Initial program 83.3%

      \[{\left(a \cdot \cos \left(\mathsf{PI}\left(\right) \cdot \frac{angle}{180}\right)\right)}^{2} + {\left(b \cdot \sin \left(\mathsf{PI}\left(\right) \cdot \frac{angle}{180}\right)\right)}^{2} \]
   2. Add Preprocessing

   3. Taylor expanded in angle around 0

      \[\leadsto \color{blue}{{a}^{2}} \]
   4. Step-by-step derivation

      1. unpow2 N/A

         \[\leadsto \color{blue}{a \cdot a} \]

      2. lower-*.f64 65.9

         \[\leadsto \color{blue}{a \cdot a} \]

   5. Applied rewrites 65.9%

      \[\leadsto \color{blue}{a \cdot a} \]

   if 3.40000000000000004e-45 < b < 1.8999999999999999e158

   1. Initial program 82.3%

      \[{\left(a \cdot \cos \left(\mathsf{PI}\left(\right) \cdot \frac{angle}{180}\right)\right)}^{2} + {\left(b \cdot \sin \left(\mathsf{PI}\left(\right) \cdot \frac{angle}{180}\right)\right)}^{2} \]
   2. Add Preprocessing

   3. Taylor expanded in angle around 0

      \[\leadsto \color{blue}{{angle}^{2} \cdot \left(\frac{-1}{32400} \cdot \left({a}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right) + \frac{1}{32400} \cdot \left({b}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right)\right) + {a}^{2}} \]
   4. Step-by-step derivation

      1. *-commutative N/A

         \[\leadsto \color{blue}{\left(\frac{-1}{32400} \cdot \left({a}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right) + \frac{1}{32400} \cdot \left({b}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right)\right) \cdot {angle}^{2}} + {a}^{2} \]

      2. lower-fma.f64 N/A

         \[\leadsto \color{blue}{\mathsf{fma}\left(\frac{-1}{32400} \cdot \left({a}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right) + \frac{1}{32400} \cdot \left({b}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right), {angle}^{2}, {a}^{2}\right)} \]

   5. Applied rewrites 53.3%

      \[\leadsto \color{blue}{\mathsf{fma}\left(\left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right) \cdot \mathsf{fma}\left(3.08641975308642 \cdot 10^{-5}, b \cdot b, -3.08641975308642 \cdot 10^{-5} \cdot \left(a \cdot a\right)\right), angle \cdot angle, a \cdot a\right)} \]

   6. Taylor expanded in a around 0

      \[\leadsto \mathsf{fma}\left(\left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right) \cdot \left(\frac{1}{32400} \cdot {b}^{2}\right), angle \cdot angle, a \cdot a\right) \]
   7. Step-by-step derivation

   8. Applied rewrites 77.0%

      \[\leadsto \mathsf{fma}\left(\left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right) \cdot \left(3.08641975308642 \cdot 10^{-5} \cdot \left(b \cdot b\right)\right), angle \cdot angle, a \cdot a\right) \]

   if 1.8999999999999999e158 < b

   1. Initial program 97.2%

      \[{\left(a \cdot \cos \left(\mathsf{PI}\left(\right) \cdot \frac{angle}{180}\right)\right)}^{2} + {\left(b \cdot \sin \left(\mathsf{PI}\left(\right) \cdot \frac{angle}{180}\right)\right)}^{2} \]
   2. Add Preprocessing

   3. Taylor expanded in angle around 0

      \[\leadsto \color{blue}{{angle}^{2} \cdot \left(\frac{-1}{32400} \cdot \left({a}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right) + \frac{1}{32400} \cdot \left({b}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right)\right) + {a}^{2}} \]
   4. Step-by-step derivation

      1. *-commutative N/A

         \[\leadsto \color{blue}{\left(\frac{-1}{32400} \cdot \left({a}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right) + \frac{1}{32400} \cdot \left({b}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right)\right) \cdot {angle}^{2}} + {a}^{2} \]

      2. lower-fma.f64 N/A

         \[\leadsto \color{blue}{\mathsf{fma}\left(\frac{-1}{32400} \cdot \left({a}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right) + \frac{1}{32400} \cdot \left({b}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right), {angle}^{2}, {a}^{2}\right)} \]

   5. Applied rewrites 43.0%

      \[\leadsto \color{blue}{\mathsf{fma}\left(\left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right) \cdot \mathsf{fma}\left(3.08641975308642 \cdot 10^{-5}, b \cdot b, -3.08641975308642 \cdot 10^{-5} \cdot \left(a \cdot a\right)\right), angle \cdot angle, a \cdot a\right)} \]

   6. Taylor expanded in a around 0

      \[\leadsto \frac{1}{32400} \cdot \color{blue}{\left({angle}^{2} \cdot \left({b}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right)\right)} \]
   7. Step-by-step derivation

   8. Applied rewrites 64.0%

      \[\leadsto \left(\left(angle \cdot angle\right) \cdot 3.08641975308642 \cdot 10^{-5}\right) \cdot \color{blue}{\left(\left(\left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right) \cdot b\right) \cdot b\right)} \]
   9. Step-by-step derivation

   10. Applied rewrites 75.5%

       \[\leadsto \color{blue}{{\left(\left(\mathsf{PI}\left(\right) \cdot b\right) \cdot angle\right)}^{2} \cdot 3.08641975308642 \cdot 10^{-5}} \]
3. Recombined 3 regimes into one program.
4. Add Preprocessing

Alternative 4: 62.3% accurate, 6.3× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;\frac{angle}{180} \leq 2 \cdot 10^{-160}:\\ \;\;\;\;a \cdot a\\ \mathbf{elif}\;\frac{angle}{180} \leq 2 \cdot 10^{+145}:\\ \;\;\;\;\mathsf{fma}\left(\left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right) \cdot \left(3.08641975308642 \cdot 10^{-5} \cdot \left(b \cdot b\right)\right), angle \cdot angle, a \cdot a\right)\\ \mathbf{else}:\\ \;\;\;\;\left(\left(\left(\left(angle \cdot angle\right) \cdot 3.08641975308642 \cdot 10^{-5}\right) \cdot \mathsf{PI}\left(\right)\right) \cdot \left(\mathsf{PI}\left(\right) \cdot b\right)\right) \cdot b\\ \end{array} \end{array} \]

FPCore

(FPCore (a b angle)
 :precision binary64
 (if (<= (/ angle 180.0) 2e-160)
   (* a a)
   (if (<= (/ angle 180.0) 2e+145)
     (fma
      (* (* (PI) (PI)) (* 3.08641975308642e-5 (* b b)))
      (* angle angle)
      (* a a))
     (* (* (* (* (* angle angle) 3.08641975308642e-5) (PI)) (* (PI) b)) b))))

Derivation

1. Split input into 3 regimes

2. if (/.f64 angle #s(literal 180 binary64)) < 2e-160

   1. Initial program 86.5%

      \[{\left(a \cdot \cos \left(\mathsf{PI}\left(\right) \cdot \frac{angle}{180}\right)\right)}^{2} + {\left(b \cdot \sin \left(\mathsf{PI}\left(\right) \cdot \frac{angle}{180}\right)\right)}^{2} \]
   2. Add Preprocessing

   3. Taylor expanded in angle around 0

      \[\leadsto \color{blue}{{a}^{2}} \]
   4. Step-by-step derivation

      1. unpow2 N/A

         \[\leadsto \color{blue}{a \cdot a} \]

      2. lower-*.f64 63.1

         \[\leadsto \color{blue}{a \cdot a} \]

   5. Applied rewrites 63.1%

      \[\leadsto \color{blue}{a \cdot a} \]

   if 2e-160 < (/.f64 angle #s(literal 180 binary64)) < 2e145

   1. Initial program 85.4%

      \[{\left(a \cdot \cos \left(\mathsf{PI}\left(\right) \cdot \frac{angle}{180}\right)\right)}^{2} + {\left(b \cdot \sin \left(\mathsf{PI}\left(\right) \cdot \frac{angle}{180}\right)\right)}^{2} \]
   2. Add Preprocessing

   3. Taylor expanded in angle around 0

      \[\leadsto \color{blue}{{angle}^{2} \cdot \left(\frac{-1}{32400} \cdot \left({a}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right) + \frac{1}{32400} \cdot \left({b}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right)\right) + {a}^{2}} \]
   4. Step-by-step derivation

      1. *-commutative N/A

         \[\leadsto \color{blue}{\left(\frac{-1}{32400} \cdot \left({a}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right) + \frac{1}{32400} \cdot \left({b}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right)\right) \cdot {angle}^{2}} + {a}^{2} \]

      2. lower-fma.f64 N/A

         \[\leadsto \color{blue}{\mathsf{fma}\left(\frac{-1}{32400} \cdot \left({a}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right) + \frac{1}{32400} \cdot \left({b}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right), {angle}^{2}, {a}^{2}\right)} \]

   5. Applied rewrites 47.2%

      \[\leadsto \color{blue}{\mathsf{fma}\left(\left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right) \cdot \mathsf{fma}\left(3.08641975308642 \cdot 10^{-5}, b \cdot b, -3.08641975308642 \cdot 10^{-5} \cdot \left(a \cdot a\right)\right), angle \cdot angle, a \cdot a\right)} \]

   6. Taylor expanded in a around 0

      \[\leadsto \mathsf{fma}\left(\left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right) \cdot \left(\frac{1}{32400} \cdot {b}^{2}\right), angle \cdot angle, a \cdot a\right) \]
   7. Step-by-step derivation

   8. Applied rewrites 78.6%

      \[\leadsto \mathsf{fma}\left(\left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right) \cdot \left(3.08641975308642 \cdot 10^{-5} \cdot \left(b \cdot b\right)\right), angle \cdot angle, a \cdot a\right) \]

   if 2e145 < (/.f64 angle #s(literal 180 binary64))

   1. Initial program 76.2%

      \[{\left(a \cdot \cos \left(\mathsf{PI}\left(\right) \cdot \frac{angle}{180}\right)\right)}^{2} + {\left(b \cdot \sin \left(\mathsf{PI}\left(\right) \cdot \frac{angle}{180}\right)\right)}^{2} \]
   2. Add Preprocessing

   3. Taylor expanded in angle around 0

      \[\leadsto \color{blue}{{angle}^{2} \cdot \left(\frac{-1}{32400} \cdot \left({a}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right) + \frac{1}{32400} \cdot \left({b}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right)\right) + {a}^{2}} \]
   4. Step-by-step derivation

      1. *-commutative N/A

         \[\leadsto \color{blue}{\left(\frac{-1}{32400} \cdot \left({a}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right) + \frac{1}{32400} \cdot \left({b}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right)\right) \cdot {angle}^{2}} + {a}^{2} \]

      2. lower-fma.f64 N/A

         \[\leadsto \color{blue}{\mathsf{fma}\left(\frac{-1}{32400} \cdot \left({a}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right) + \frac{1}{32400} \cdot \left({b}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right), {angle}^{2}, {a}^{2}\right)} \]

   5. Applied rewrites 15.4%

      \[\leadsto \color{blue}{\mathsf{fma}\left(\left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right) \cdot \mathsf{fma}\left(3.08641975308642 \cdot 10^{-5}, b \cdot b, -3.08641975308642 \cdot 10^{-5} \cdot \left(a \cdot a\right)\right), angle \cdot angle, a \cdot a\right)} \]

   6. Taylor expanded in a around 0

      \[\leadsto \frac{1}{32400} \cdot \color{blue}{\left({angle}^{2} \cdot \left({b}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right)\right)} \]
   7. Step-by-step derivation

   8. Applied rewrites 49.1%

      \[\leadsto \left(\left(angle \cdot angle\right) \cdot 3.08641975308642 \cdot 10^{-5}\right) \cdot \color{blue}{\left(\left(\left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right) \cdot b\right) \cdot b\right)} \]
   9. Step-by-step derivation

   10. Applied rewrites 67.9%

       \[\leadsto \left(\left(\left(\left(angle \cdot angle\right) \cdot 3.08641975308642 \cdot 10^{-5}\right) \cdot \mathsf{PI}\left(\right)\right) \cdot \left(\mathsf{PI}\left(\right) \cdot b\right)\right) \cdot b \]
3. Recombined 3 regimes into one program.
4. Add Preprocessing

Alternative 5: 62.3% accurate, 6.3× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;\frac{angle}{180} \leq 2 \cdot 10^{-160}:\\ \;\;\;\;a \cdot a\\ \mathbf{elif}\;\frac{angle}{180} \leq 2 \cdot 10^{+145}:\\ \;\;\;\;\mathsf{fma}\left(\left(\left(\left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right) \cdot 3.08641975308642 \cdot 10^{-5}\right) \cdot b\right) \cdot b, angle \cdot angle, a \cdot a\right)\\ \mathbf{else}:\\ \;\;\;\;\left(\left(\left(\left(angle \cdot angle\right) \cdot 3.08641975308642 \cdot 10^{-5}\right) \cdot \mathsf{PI}\left(\right)\right) \cdot \left(\mathsf{PI}\left(\right) \cdot b\right)\right) \cdot b\\ \end{array} \end{array} \]

FPCore

(FPCore (a b angle)
 :precision binary64
 (if (<= (/ angle 180.0) 2e-160)
   (* a a)
   (if (<= (/ angle 180.0) 2e+145)
     (fma
      (* (* (* (* (PI) (PI)) 3.08641975308642e-5) b) b)
      (* angle angle)
      (* a a))
     (* (* (* (* (* angle angle) 3.08641975308642e-5) (PI)) (* (PI) b)) b))))

Derivation

1. Split input into 3 regimes

2. if (/.f64 angle #s(literal 180 binary64)) < 2e-160

   1. Initial program 86.5%

      \[{\left(a \cdot \cos \left(\mathsf{PI}\left(\right) \cdot \frac{angle}{180}\right)\right)}^{2} + {\left(b \cdot \sin \left(\mathsf{PI}\left(\right) \cdot \frac{angle}{180}\right)\right)}^{2} \]
   2. Add Preprocessing

   3. Taylor expanded in angle around 0

      \[\leadsto \color{blue}{{a}^{2}} \]
   4. Step-by-step derivation

      1. unpow2 N/A

         \[\leadsto \color{blue}{a \cdot a} \]

      2. lower-*.f64 63.1

         \[\leadsto \color{blue}{a \cdot a} \]

   5. Applied rewrites 63.1%

      \[\leadsto \color{blue}{a \cdot a} \]

   if 2e-160 < (/.f64 angle #s(literal 180 binary64)) < 2e145

   1. Initial program 85.4%

      \[{\left(a \cdot \cos \left(\mathsf{PI}\left(\right) \cdot \frac{angle}{180}\right)\right)}^{2} + {\left(b \cdot \sin \left(\mathsf{PI}\left(\right) \cdot \frac{angle}{180}\right)\right)}^{2} \]
   2. Add Preprocessing

   3. Taylor expanded in angle around 0

      \[\leadsto \color{blue}{{angle}^{2} \cdot \left(\frac{-1}{32400} \cdot \left({a}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right) + \frac{1}{32400} \cdot \left({b}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right)\right) + {a}^{2}} \]
   4. Step-by-step derivation

      1. *-commutative N/A

         \[\leadsto \color{blue}{\left(\frac{-1}{32400} \cdot \left({a}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right) + \frac{1}{32400} \cdot \left({b}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right)\right) \cdot {angle}^{2}} + {a}^{2} \]

      2. lower-fma.f64 N/A

         \[\leadsto \color{blue}{\mathsf{fma}\left(\frac{-1}{32400} \cdot \left({a}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right) + \frac{1}{32400} \cdot \left({b}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right), {angle}^{2}, {a}^{2}\right)} \]

   5. Applied rewrites 47.2%

      \[\leadsto \color{blue}{\mathsf{fma}\left(\left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right) \cdot \mathsf{fma}\left(3.08641975308642 \cdot 10^{-5}, b \cdot b, -3.08641975308642 \cdot 10^{-5} \cdot \left(a \cdot a\right)\right), angle \cdot angle, a \cdot a\right)} \]

   6. Taylor expanded in a around 0

      \[\leadsto \mathsf{fma}\left(\frac{1}{32400} \cdot \left({b}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right), \color{blue}{angle} \cdot angle, a \cdot a\right) \]
   7. Step-by-step derivation

   8. Applied rewrites 78.5%

      \[\leadsto \mathsf{fma}\left(\left(\left(\left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right) \cdot 3.08641975308642 \cdot 10^{-5}\right) \cdot b\right) \cdot b, \color{blue}{angle} \cdot angle, a \cdot a\right) \]

   if 2e145 < (/.f64 angle #s(literal 180 binary64))

   1. Initial program 76.2%

      \[{\left(a \cdot \cos \left(\mathsf{PI}\left(\right) \cdot \frac{angle}{180}\right)\right)}^{2} + {\left(b \cdot \sin \left(\mathsf{PI}\left(\right) \cdot \frac{angle}{180}\right)\right)}^{2} \]
   2. Add Preprocessing

   3. Taylor expanded in angle around 0

      \[\leadsto \color{blue}{{angle}^{2} \cdot \left(\frac{-1}{32400} \cdot \left({a}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right) + \frac{1}{32400} \cdot \left({b}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right)\right) + {a}^{2}} \]
   4. Step-by-step derivation

      1. *-commutative N/A

         \[\leadsto \color{blue}{\left(\frac{-1}{32400} \cdot \left({a}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right) + \frac{1}{32400} \cdot \left({b}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right)\right) \cdot {angle}^{2}} + {a}^{2} \]

      2. lower-fma.f64 N/A

         \[\leadsto \color{blue}{\mathsf{fma}\left(\frac{-1}{32400} \cdot \left({a}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right) + \frac{1}{32400} \cdot \left({b}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right), {angle}^{2}, {a}^{2}\right)} \]

   5. Applied rewrites 15.4%

      \[\leadsto \color{blue}{\mathsf{fma}\left(\left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right) \cdot \mathsf{fma}\left(3.08641975308642 \cdot 10^{-5}, b \cdot b, -3.08641975308642 \cdot 10^{-5} \cdot \left(a \cdot a\right)\right), angle \cdot angle, a \cdot a\right)} \]

   6. Taylor expanded in a around 0

      \[\leadsto \frac{1}{32400} \cdot \color{blue}{\left({angle}^{2} \cdot \left({b}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right)\right)} \]
   7. Step-by-step derivation

   8. Applied rewrites 49.1%

      \[\leadsto \left(\left(angle \cdot angle\right) \cdot 3.08641975308642 \cdot 10^{-5}\right) \cdot \color{blue}{\left(\left(\left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right) \cdot b\right) \cdot b\right)} \]
   9. Step-by-step derivation

   10. Applied rewrites 67.9%

       \[\leadsto \left(\left(\left(\left(angle \cdot angle\right) \cdot 3.08641975308642 \cdot 10^{-5}\right) \cdot \mathsf{PI}\left(\right)\right) \cdot \left(\mathsf{PI}\left(\right) \cdot b\right)\right) \cdot b \]
3. Recombined 3 regimes into one program.
4. Add Preprocessing

Alternative 6: 55.5% accurate, 8.3× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;a \leq 2.8 \cdot 10^{+90}:\\ \;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(b \cdot b, 3.08641975308642 \cdot 10^{-5}, \left(a \cdot a\right) \cdot -3.08641975308642 \cdot 10^{-5}\right) \cdot \left(\left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right) \cdot angle\right), angle, a \cdot a\right)\\ \mathbf{else}:\\ \;\;\;\;a \cdot a\\ \end{array} \end{array} \]

FPCore

(FPCore (a b angle)
 :precision binary64
 (if (<= a 2.8e+90)
   (fma
    (*
     (fma (* b b) 3.08641975308642e-5 (* (* a a) -3.08641975308642e-5))
     (* (* (PI) (PI)) angle))
    angle
    (* a a))
   (* a a)))

Derivation

1. Split input into 2 regimes

2. if a < 2.8e90

   1. Initial program 82.5%

      \[{\left(a \cdot \cos \left(\mathsf{PI}\left(\right) \cdot \frac{angle}{180}\right)\right)}^{2} + {\left(b \cdot \sin \left(\mathsf{PI}\left(\right) \cdot \frac{angle}{180}\right)\right)}^{2} \]
   2. Add Preprocessing

   3. Taylor expanded in angle around 0

      \[\leadsto \color{blue}{{angle}^{2} \cdot \left(\frac{-1}{32400} \cdot \left({a}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right) + \frac{1}{32400} \cdot \left({b}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right)\right) + {a}^{2}} \]
   4. Step-by-step derivation

      1. *-commutative N/A

         \[\leadsto \color{blue}{\left(\frac{-1}{32400} \cdot \left({a}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right) + \frac{1}{32400} \cdot \left({b}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right)\right) \cdot {angle}^{2}} + {a}^{2} \]

      2. lower-fma.f64 N/A

         \[\leadsto \color{blue}{\mathsf{fma}\left(\frac{-1}{32400} \cdot \left({a}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right) + \frac{1}{32400} \cdot \left({b}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right), {angle}^{2}, {a}^{2}\right)} \]

   5. Applied rewrites 47.7%

      \[\leadsto \color{blue}{\mathsf{fma}\left(\left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right) \cdot \mathsf{fma}\left(3.08641975308642 \cdot 10^{-5}, b \cdot b, -3.08641975308642 \cdot 10^{-5} \cdot \left(a \cdot a\right)\right), angle \cdot angle, a \cdot a\right)} \]
   6. Step-by-step derivation

   7. Applied rewrites 50.1%

      \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(b \cdot b, 3.08641975308642 \cdot 10^{-5}, \left(a \cdot a\right) \cdot -3.08641975308642 \cdot 10^{-5}\right) \cdot \left(\left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right) \cdot angle\right), \color{blue}{angle}, a \cdot a\right) \]

   if 2.8e90 < a

   1. Initial program 95.4%

      \[{\left(a \cdot \cos \left(\mathsf{PI}\left(\right) \cdot \frac{angle}{180}\right)\right)}^{2} + {\left(b \cdot \sin \left(\mathsf{PI}\left(\right) \cdot \frac{angle}{180}\right)\right)}^{2} \]
   2. Add Preprocessing

   3. Taylor expanded in angle around 0

      \[\leadsto \color{blue}{{a}^{2}} \]
   4. Step-by-step derivation

      1. unpow2 N/A

         \[\leadsto \color{blue}{a \cdot a} \]

      2. lower-*.f64 90.5

         \[\leadsto \color{blue}{a \cdot a} \]

   5. Applied rewrites 90.5%

      \[\leadsto \color{blue}{a \cdot a} \]
3. Recombined 2 regimes into one program.
4. Add Preprocessing

Alternative 7: 55.5% accurate, 8.3× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;a \leq 2.8 \cdot 10^{+90}:\\ \;\;\;\;\mathsf{fma}\left(\mathsf{PI}\left(\right), \mathsf{PI}\left(\right) \cdot \left(\left(\mathsf{fma}\left(b \cdot b, 3.08641975308642 \cdot 10^{-5}, \left(a \cdot a\right) \cdot -3.08641975308642 \cdot 10^{-5}\right) \cdot angle\right) \cdot angle\right), a \cdot a\right)\\ \mathbf{else}:\\ \;\;\;\;a \cdot a\\ \end{array} \end{array} \]

FPCore

(FPCore (a b angle)
 :precision binary64
 (if (<= a 2.8e+90)
   (fma
    (PI)
    (*
     (PI)
     (*
      (*
       (fma (* b b) 3.08641975308642e-5 (* (* a a) -3.08641975308642e-5))
       angle)
      angle))
    (* a a))
   (* a a)))

Derivation

1. Split input into 2 regimes

2. if a < 2.8e90

   1. Initial program 82.5%

      \[{\left(a \cdot \cos \left(\mathsf{PI}\left(\right) \cdot \frac{angle}{180}\right)\right)}^{2} + {\left(b \cdot \sin \left(\mathsf{PI}\left(\right) \cdot \frac{angle}{180}\right)\right)}^{2} \]
   2. Add Preprocessing

   3. Taylor expanded in angle around 0

      \[\leadsto \color{blue}{{angle}^{2} \cdot \left(\frac{-1}{32400} \cdot \left({a}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right) + \frac{1}{32400} \cdot \left({b}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right)\right) + {a}^{2}} \]
   4. Step-by-step derivation

      1. *-commutative N/A

         \[\leadsto \color{blue}{\left(\frac{-1}{32400} \cdot \left({a}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right) + \frac{1}{32400} \cdot \left({b}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right)\right) \cdot {angle}^{2}} + {a}^{2} \]

      2. lower-fma.f64 N/A

         \[\leadsto \color{blue}{\mathsf{fma}\left(\frac{-1}{32400} \cdot \left({a}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right) + \frac{1}{32400} \cdot \left({b}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right), {angle}^{2}, {a}^{2}\right)} \]

   5. Applied rewrites 47.7%

      \[\leadsto \color{blue}{\mathsf{fma}\left(\left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right) \cdot \mathsf{fma}\left(3.08641975308642 \cdot 10^{-5}, b \cdot b, -3.08641975308642 \cdot 10^{-5} \cdot \left(a \cdot a\right)\right), angle \cdot angle, a \cdot a\right)} \]
   6. Step-by-step derivation

   7. Applied rewrites 50.1%

      \[\leadsto \mathsf{fma}\left(\mathsf{PI}\left(\right), \color{blue}{\mathsf{PI}\left(\right) \cdot \left(\left(\mathsf{fma}\left(b \cdot b, 3.08641975308642 \cdot 10^{-5}, \left(a \cdot a\right) \cdot -3.08641975308642 \cdot 10^{-5}\right) \cdot angle\right) \cdot angle\right)}, a \cdot a\right) \]

   if 2.8e90 < a

   1. Initial program 95.4%

      \[{\left(a \cdot \cos \left(\mathsf{PI}\left(\right) \cdot \frac{angle}{180}\right)\right)}^{2} + {\left(b \cdot \sin \left(\mathsf{PI}\left(\right) \cdot \frac{angle}{180}\right)\right)}^{2} \]
   2. Add Preprocessing

   3. Taylor expanded in angle around 0

      \[\leadsto \color{blue}{{a}^{2}} \]
   4. Step-by-step derivation

      1. unpow2 N/A

         \[\leadsto \color{blue}{a \cdot a} \]

      2. lower-*.f64 90.5

         \[\leadsto \color{blue}{a \cdot a} \]

   5. Applied rewrites 90.5%

      \[\leadsto \color{blue}{a \cdot a} \]
3. Recombined 2 regimes into one program.
4. Add Preprocessing

Alternative 8: 60.5% accurate, 12.1× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;b \leq 5.9 \cdot 10^{+190}:\\ \;\;\;\;a \cdot a\\ \mathbf{else}:\\ \;\;\;\;\left(3.08641975308642 \cdot 10^{-5} \cdot \left(\left(\left(angle \cdot angle\right) \cdot b\right) \cdot b\right)\right) \cdot \left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right)\\ \end{array} \end{array} \]

FPCore

(FPCore (a b angle)
 :precision binary64
 (if (<= b 5.9e+190)
   (* a a)
   (* (* 3.08641975308642e-5 (* (* (* angle angle) b) b)) (* (PI) (PI)))))

Derivation

1. Split input into 2 regimes

2. if b < 5.89999999999999972e190

   1. Initial program 82.9%

      \[{\left(a \cdot \cos \left(\mathsf{PI}\left(\right) \cdot \frac{angle}{180}\right)\right)}^{2} + {\left(b \cdot \sin \left(\mathsf{PI}\left(\right) \cdot \frac{angle}{180}\right)\right)}^{2} \]
   2. Add Preprocessing

   3. Taylor expanded in angle around 0

      \[\leadsto \color{blue}{{a}^{2}} \]
   4. Step-by-step derivation

      1. unpow2 N/A

         \[\leadsto \color{blue}{a \cdot a} \]

      2. lower-*.f64 63.3

         \[\leadsto \color{blue}{a \cdot a} \]

   5. Applied rewrites 63.3%

      \[\leadsto \color{blue}{a \cdot a} \]

   if 5.89999999999999972e190 < b

   1. Initial program 99.9%

      \[{\left(a \cdot \cos \left(\mathsf{PI}\left(\right) \cdot \frac{angle}{180}\right)\right)}^{2} + {\left(b \cdot \sin \left(\mathsf{PI}\left(\right) \cdot \frac{angle}{180}\right)\right)}^{2} \]
   2. Add Preprocessing

   3. Taylor expanded in angle around 0

      \[\leadsto \color{blue}{{angle}^{2} \cdot \left(\frac{-1}{32400} \cdot \left({a}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right) + \frac{1}{32400} \cdot \left({b}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right)\right) + {a}^{2}} \]
   4. Step-by-step derivation

      1. *-commutative N/A

         \[\leadsto \color{blue}{\left(\frac{-1}{32400} \cdot \left({a}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right) + \frac{1}{32400} \cdot \left({b}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right)\right) \cdot {angle}^{2}} + {a}^{2} \]

      2. lower-fma.f64 N/A

         \[\leadsto \color{blue}{\mathsf{fma}\left(\frac{-1}{32400} \cdot \left({a}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right) + \frac{1}{32400} \cdot \left({b}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right), {angle}^{2}, {a}^{2}\right)} \]

   5. Applied rewrites 47.8%

      \[\leadsto \color{blue}{\mathsf{fma}\left(\left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right) \cdot \mathsf{fma}\left(3.08641975308642 \cdot 10^{-5}, b \cdot b, -3.08641975308642 \cdot 10^{-5} \cdot \left(a \cdot a\right)\right), angle \cdot angle, a \cdot a\right)} \]

   6. Taylor expanded in a around 0

      \[\leadsto \frac{1}{32400} \cdot \color{blue}{\left({angle}^{2} \cdot \left({b}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right)\right)} \]
   7. Step-by-step derivation

   8. Applied rewrites 71.3%

      \[\leadsto \left(\left(angle \cdot angle\right) \cdot 3.08641975308642 \cdot 10^{-5}\right) \cdot \color{blue}{\left(\left(\left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right) \cdot b\right) \cdot b\right)} \]

   9. Taylor expanded in a around 0

      \[\leadsto \frac{1}{32400} \cdot \color{blue}{\left({angle}^{2} \cdot \left({b}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right)\right)} \]
   10. Step-by-step derivation

   11. Applied rewrites 77.1%

       \[\leadsto \left(3.08641975308642 \cdot 10^{-5} \cdot \left(\left(\left(angle \cdot angle\right) \cdot b\right) \cdot b\right)\right) \cdot \color{blue}{\left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right)} \]
3. Recombined 2 regimes into one program.
4. Add Preprocessing

Alternative 9: 56.9% accurate, 74.7× speedup

Math

\[\begin{array}{l} \\ a \cdot a \end{array} \]

FPCore

(FPCore (a b angle) :precision binary64 (* a a))

C

double code(double a, double b, double angle) {
	return a * a;
}

Fortran

real(8) function code(a, b, angle)
    real(8), intent (in) :: a
    real(8), intent (in) :: b
    real(8), intent (in) :: angle
    code = a * a
end function

Java

public static double code(double a, double b, double angle) {
	return a * a;
}

Python

def code(a, b, angle):
	return a * a

Julia

function code(a, b, angle)
	return Float64(a * a)
end

MATLAB

function tmp = code(a, b, angle)
	tmp = a * a;
end

Wolfram

code[a_, b_, angle_] := N[(a * a), $MachinePrecision]

Derivation

1. Initial program 85.2%

   \[{\left(a \cdot \cos \left(\mathsf{PI}\left(\right) \cdot \frac{angle}{180}\right)\right)}^{2} + {\left(b \cdot \sin \left(\mathsf{PI}\left(\right) \cdot \frac{angle}{180}\right)\right)}^{2} \]
2. Add Preprocessing

3. Taylor expanded in angle around 0

   \[\leadsto \color{blue}{{a}^{2}} \]
4. Step-by-step derivation

   1. unpow2 N/A

      \[\leadsto \color{blue}{a \cdot a} \]

   2. lower-*.f64 59.9

      \[\leadsto \color{blue}{a \cdot a} \]

5. Applied rewrites 59.9%

   \[\leadsto \color{blue}{a \cdot a} \]
6. Add Preprocessing

Reproduce

herbie shell --seed 2024327 
(FPCore (a b angle)
  :name "ab-angle->ABCF C"
  :precision binary64
  (+ (pow (* a (cos (* (PI) (/ angle 180.0)))) 2.0) (pow (* b (sin (* (PI) (/ angle 180.0)))) 2.0)))