ab-angle->ABCF A

Percentage Accurate: 79.5% → 79.4%
Time: 7.5s
Alternatives: 11
Speedup: N/A×

Specification

?
\[\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} \]
(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:

Local Percentage Accuracy vs ?

The average percentage accuracy by input value. Horizontal axis shows value of an input variable; the variable is choosen in the title. Vertical axis is accuracy; higher is better. Red represent the original program, while blue represents Herbie's suggestion. These can be toggled with buttons below the plot. The line is an average while dots represent individual samples.

Accuracy vs Speed?

Herbie found 11 alternatives:

AlternativeAccuracySpeedup
The accuracy (vertical axis) and speed (horizontal axis) of each alternatives. Up and to the right is better. The red square shows the initial program, and each blue circle shows an alternative.The line shows the best available speed-accuracy tradeoffs.

Initial Program: 79.5% accurate, 1.0× speedup?

\[\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} \]
(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}

Alternative 1: 79.4% accurate, 1.3× speedup?

\[\begin{array}{l} \\ {\left(1 \cdot b\right)}^{2} + {\left(\sin \left(\mathsf{PI}\left(\right) \cdot \frac{angle}{180}\right) \cdot a\right)}^{2} \end{array} \]
(FPCore (a b angle)
 :precision binary64
 (+ (pow (* 1.0 b) 2.0) (pow (* (sin (* (PI) (/ angle 180.0))) a) 2.0)))
\begin{array}{l}

\\
{\left(1 \cdot b\right)}^{2} + {\left(\sin \left(\mathsf{PI}\left(\right) \cdot \frac{angle}{180}\right) \cdot a\right)}^{2}
\end{array}
Derivation
  1. Initial program 80.0%

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

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

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

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

    Alternative 2: 79.4% accurate, 1.9× speedup?

    \[\begin{array}{l} \\ \mathsf{fma}\left(\left(1 \cdot b\right) \cdot 1, b, {\left(\sin \left(\left(0.005555555555555556 \cdot angle\right) \cdot \mathsf{PI}\left(\right)\right) \cdot a\right)}^{2}\right) \end{array} \]
    (FPCore (a b angle)
     :precision binary64
     (fma
      (* (* 1.0 b) 1.0)
      b
      (pow (* (sin (* (* 0.005555555555555556 angle) (PI))) a) 2.0)))
    \begin{array}{l}
    
    \\
    \mathsf{fma}\left(\left(1 \cdot b\right) \cdot 1, b, {\left(\sin \left(\left(0.005555555555555556 \cdot angle\right) \cdot \mathsf{PI}\left(\right)\right) \cdot a\right)}^{2}\right)
    \end{array}
    
    Derivation
    1. Initial program 80.0%

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

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

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

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

          \[\leadsto \color{blue}{{\left(b \cdot 1\right)}^{2} + {\left(a \cdot \sin \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}} \]
        3. lift-pow.f64N/A

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

          \[\leadsto \color{blue}{\left(b \cdot 1\right) \cdot \left(b \cdot 1\right)} + {\left(a \cdot \sin \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} \]
        5. lift-*.f64N/A

          \[\leadsto \left(b \cdot 1\right) \cdot \color{blue}{\left(b \cdot 1\right)} + {\left(a \cdot \sin \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} \]
        6. *-commutativeN/A

          \[\leadsto \left(b \cdot 1\right) \cdot \color{blue}{\left(1 \cdot b\right)} + {\left(a \cdot \sin \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} \]
        7. associate-*r*N/A

          \[\leadsto \color{blue}{\left(\left(b \cdot 1\right) \cdot 1\right) \cdot b} + {\left(a \cdot \sin \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} \]
        8. lower-fma.f64N/A

          \[\leadsto \color{blue}{\mathsf{fma}\left(\left(b \cdot 1\right) \cdot 1, b, {\left(a \cdot \sin \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}\right)} \]
      3. Applied rewrites80.4%

        \[\leadsto \color{blue}{\mathsf{fma}\left(\left(1 \cdot b\right) \cdot 1, b, {\left(\sin \left(\mathsf{PI}\left(\right) \cdot \left(0.005555555555555556 \cdot angle\right)\right) \cdot a\right)}^{2}\right)} \]
      4. Final simplification80.4%

        \[\leadsto \mathsf{fma}\left(\left(1 \cdot b\right) \cdot 1, b, {\left(\sin \left(\left(0.005555555555555556 \cdot angle\right) \cdot \mathsf{PI}\left(\right)\right) \cdot a\right)}^{2}\right) \]
      5. Add Preprocessing

      Alternative 3: 79.3% accurate, 1.9× speedup?

      \[\begin{array}{l} \\ \mathsf{fma}\left(\left(1 \cdot b\right) \cdot 1, b, {\left(\sin \left(\left(\mathsf{PI}\left(\right) \cdot angle\right) \cdot 0.005555555555555556\right) \cdot a\right)}^{2}\right) \end{array} \]
      (FPCore (a b angle)
       :precision binary64
       (fma
        (* (* 1.0 b) 1.0)
        b
        (pow (* (sin (* (* (PI) angle) 0.005555555555555556)) a) 2.0)))
      \begin{array}{l}
      
      \\
      \mathsf{fma}\left(\left(1 \cdot b\right) \cdot 1, b, {\left(\sin \left(\left(\mathsf{PI}\left(\right) \cdot angle\right) \cdot 0.005555555555555556\right) \cdot a\right)}^{2}\right)
      \end{array}
      
      Derivation
      1. Initial program 80.0%

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

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

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

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

            \[\leadsto \color{blue}{{\left(b \cdot 1\right)}^{2} + {\left(a \cdot \sin \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}} \]
          3. lift-pow.f64N/A

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

            \[\leadsto \color{blue}{\left(b \cdot 1\right) \cdot \left(b \cdot 1\right)} + {\left(a \cdot \sin \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} \]
          5. lift-*.f64N/A

            \[\leadsto \left(b \cdot 1\right) \cdot \color{blue}{\left(b \cdot 1\right)} + {\left(a \cdot \sin \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} \]
          6. *-commutativeN/A

            \[\leadsto \left(b \cdot 1\right) \cdot \color{blue}{\left(1 \cdot b\right)} + {\left(a \cdot \sin \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} \]
          7. associate-*r*N/A

            \[\leadsto \color{blue}{\left(\left(b \cdot 1\right) \cdot 1\right) \cdot b} + {\left(a \cdot \sin \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} \]
          8. lower-fma.f64N/A

            \[\leadsto \color{blue}{\mathsf{fma}\left(\left(b \cdot 1\right) \cdot 1, b, {\left(a \cdot \sin \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}\right)} \]
        3. Applied rewrites80.4%

          \[\leadsto \color{blue}{\mathsf{fma}\left(\left(1 \cdot b\right) \cdot 1, b, {\left(\sin \left(\mathsf{PI}\left(\right) \cdot \left(0.005555555555555556 \cdot angle\right)\right) \cdot a\right)}^{2}\right)} \]
        4. Step-by-step derivation
          1. lift-*.f64N/A

            \[\leadsto \mathsf{fma}\left(\left(1 \cdot b\right) \cdot 1, b, {\color{blue}{\left(\sin \left(\mathsf{PI}\left(\right) \cdot \left(\frac{1}{180} \cdot angle\right)\right) \cdot a\right)}}^{2}\right) \]
          2. lift-*.f64N/A

            \[\leadsto \mathsf{fma}\left(\left(1 \cdot b\right) \cdot 1, b, {\left(\sin \color{blue}{\left(\mathsf{PI}\left(\right) \cdot \left(\frac{1}{180} \cdot angle\right)\right)} \cdot a\right)}^{2}\right) \]
          3. lift-*.f64N/A

            \[\leadsto \mathsf{fma}\left(\left(1 \cdot b\right) \cdot 1, b, {\left(\sin \left(\mathsf{PI}\left(\right) \cdot \color{blue}{\left(\frac{1}{180} \cdot angle\right)}\right) \cdot a\right)}^{2}\right) \]
          4. *-commutativeN/A

            \[\leadsto \mathsf{fma}\left(\left(1 \cdot b\right) \cdot 1, b, {\left(\sin \left(\mathsf{PI}\left(\right) \cdot \color{blue}{\left(angle \cdot \frac{1}{180}\right)}\right) \cdot a\right)}^{2}\right) \]
          5. metadata-evalN/A

            \[\leadsto \mathsf{fma}\left(\left(1 \cdot b\right) \cdot 1, b, {\left(\sin \left(\mathsf{PI}\left(\right) \cdot \left(angle \cdot \color{blue}{\frac{1}{180}}\right)\right) \cdot a\right)}^{2}\right) \]
          6. div-invN/A

            \[\leadsto \mathsf{fma}\left(\left(1 \cdot b\right) \cdot 1, b, {\left(\sin \left(\mathsf{PI}\left(\right) \cdot \color{blue}{\frac{angle}{180}}\right) \cdot a\right)}^{2}\right) \]
          7. lift-/.f64N/A

            \[\leadsto \mathsf{fma}\left(\left(1 \cdot b\right) \cdot 1, b, {\left(\sin \left(\mathsf{PI}\left(\right) \cdot \color{blue}{\frac{angle}{180}}\right) \cdot a\right)}^{2}\right) \]
          8. *-commutativeN/A

            \[\leadsto \mathsf{fma}\left(\left(1 \cdot b\right) \cdot 1, b, {\left(\sin \color{blue}{\left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)} \cdot a\right)}^{2}\right) \]
          9. lift-*.f64N/A

            \[\leadsto \mathsf{fma}\left(\left(1 \cdot b\right) \cdot 1, b, {\left(\sin \color{blue}{\left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)} \cdot a\right)}^{2}\right) \]
          10. *-commutativeN/A

            \[\leadsto \mathsf{fma}\left(\left(1 \cdot b\right) \cdot 1, b, {\color{blue}{\left(a \cdot \sin \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}}^{2}\right) \]
          11. lift-*.f6480.4

            \[\leadsto \mathsf{fma}\left(\left(1 \cdot b\right) \cdot 1, b, {\color{blue}{\left(a \cdot \sin \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}}^{2}\right) \]
        5. Applied rewrites80.3%

          \[\leadsto \mathsf{fma}\left(\left(1 \cdot b\right) \cdot 1, b, \color{blue}{{\left(a \cdot \sin \left(\left(\mathsf{PI}\left(\right) \cdot angle\right) \cdot 0.005555555555555556\right)\right)}^{2}}\right) \]
        6. Final simplification80.3%

          \[\leadsto \mathsf{fma}\left(\left(1 \cdot b\right) \cdot 1, b, {\left(\sin \left(\left(\mathsf{PI}\left(\right) \cdot angle\right) \cdot 0.005555555555555556\right) \cdot a\right)}^{2}\right) \]
        7. Add Preprocessing

        Alternative 4: 66.2% accurate, 1.9× speedup?

        \[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;a \leq 2.1 \cdot 10^{+30}:\\ \;\;\;\;b \cdot b\\ \mathbf{else}:\\ \;\;\;\;{\left(\left(\left(0.005555555555555556 \cdot \mathsf{PI}\left(\right)\right) \cdot angle\right) \cdot a\right)}^{2} + {\left(1 \cdot b\right)}^{2}\\ \end{array} \end{array} \]
        (FPCore (a b angle)
         :precision binary64
         (if (<= a 2.1e+30)
           (* b b)
           (+
            (pow (* (* (* 0.005555555555555556 (PI)) angle) a) 2.0)
            (pow (* 1.0 b) 2.0))))
        \begin{array}{l}
        
        \\
        \begin{array}{l}
        \mathbf{if}\;a \leq 2.1 \cdot 10^{+30}:\\
        \;\;\;\;b \cdot b\\
        
        \mathbf{else}:\\
        \;\;\;\;{\left(\left(\left(0.005555555555555556 \cdot \mathsf{PI}\left(\right)\right) \cdot angle\right) \cdot a\right)}^{2} + {\left(1 \cdot b\right)}^{2}\\
        
        
        \end{array}
        \end{array}
        
        Derivation
        1. Split input into 2 regimes
        2. if a < 2.1e30

          1. Initial program 75.5%

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

            \[\leadsto \color{blue}{{b}^{2}} \]
          4. Step-by-step derivation
            1. unpow2N/A

              \[\leadsto \color{blue}{b \cdot b} \]
            2. lower-*.f6462.7

              \[\leadsto \color{blue}{b \cdot b} \]
          5. Applied rewrites62.7%

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

          if 2.1e30 < a

          1. Initial program 95.6%

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

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

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

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

                \[\leadsto {\left(\frac{1}{180} \cdot \color{blue}{\left(\left(angle \cdot \mathsf{PI}\left(\right)\right) \cdot a\right)}\right)}^{2} + {\left(b \cdot 1\right)}^{2} \]
              2. associate-*r*N/A

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

                \[\leadsto {\left(\color{blue}{\left(\left(angle \cdot \mathsf{PI}\left(\right)\right) \cdot \frac{1}{180}\right)} \cdot a\right)}^{2} + {\left(b \cdot 1\right)}^{2} \]
              4. associate-*r*N/A

                \[\leadsto {\left(\color{blue}{\left(angle \cdot \left(\mathsf{PI}\left(\right) \cdot \frac{1}{180}\right)\right)} \cdot a\right)}^{2} + {\left(b \cdot 1\right)}^{2} \]
              5. *-commutativeN/A

                \[\leadsto {\left(\left(angle \cdot \color{blue}{\left(\frac{1}{180} \cdot \mathsf{PI}\left(\right)\right)}\right) \cdot a\right)}^{2} + {\left(b \cdot 1\right)}^{2} \]
              6. lower-*.f64N/A

                \[\leadsto {\color{blue}{\left(\left(angle \cdot \left(\frac{1}{180} \cdot \mathsf{PI}\left(\right)\right)\right) \cdot a\right)}}^{2} + {\left(b \cdot 1\right)}^{2} \]
              7. *-commutativeN/A

                \[\leadsto {\left(\color{blue}{\left(\left(\frac{1}{180} \cdot \mathsf{PI}\left(\right)\right) \cdot angle\right)} \cdot a\right)}^{2} + {\left(b \cdot 1\right)}^{2} \]
              8. lower-*.f64N/A

                \[\leadsto {\left(\color{blue}{\left(\left(\frac{1}{180} \cdot \mathsf{PI}\left(\right)\right) \cdot angle\right)} \cdot a\right)}^{2} + {\left(b \cdot 1\right)}^{2} \]
              9. lower-*.f64N/A

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

                \[\leadsto {\left(\left(\left(0.005555555555555556 \cdot \color{blue}{\mathsf{PI}\left(\right)}\right) \cdot angle\right) \cdot a\right)}^{2} + {\left(b \cdot 1\right)}^{2} \]
            4. Applied rewrites93.2%

              \[\leadsto {\color{blue}{\left(\left(\left(0.005555555555555556 \cdot \mathsf{PI}\left(\right)\right) \cdot angle\right) \cdot a\right)}}^{2} + {\left(b \cdot 1\right)}^{2} \]
          5. Recombined 2 regimes into one program.
          6. Final simplification69.6%

            \[\leadsto \begin{array}{l} \mathbf{if}\;a \leq 2.1 \cdot 10^{+30}:\\ \;\;\;\;b \cdot b\\ \mathbf{else}:\\ \;\;\;\;{\left(\left(\left(0.005555555555555556 \cdot \mathsf{PI}\left(\right)\right) \cdot angle\right) \cdot a\right)}^{2} + {\left(1 \cdot b\right)}^{2}\\ \end{array} \]
          7. Add Preprocessing

          Alternative 5: 65.8% accurate, 3.1× speedup?

          \[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;a \leq 2.1 \cdot 10^{+30}:\\ \;\;\;\;b \cdot b\\ \mathbf{else}:\\ \;\;\;\;\mathsf{fma}\left(\left(\left(3.08641975308642 \cdot 10^{-5} \cdot \left(\left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right) \cdot a\right)\right) \cdot angle\right) \cdot angle, a, {\left(1 \cdot b\right)}^{2}\right)\\ \end{array} \end{array} \]
          (FPCore (a b angle)
           :precision binary64
           (if (<= a 2.1e+30)
             (* b b)
             (fma
              (* (* (* 3.08641975308642e-5 (* (* (PI) (PI)) a)) angle) angle)
              a
              (pow (* 1.0 b) 2.0))))
          \begin{array}{l}
          
          \\
          \begin{array}{l}
          \mathbf{if}\;a \leq 2.1 \cdot 10^{+30}:\\
          \;\;\;\;b \cdot b\\
          
          \mathbf{else}:\\
          \;\;\;\;\mathsf{fma}\left(\left(\left(3.08641975308642 \cdot 10^{-5} \cdot \left(\left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right) \cdot a\right)\right) \cdot angle\right) \cdot angle, a, {\left(1 \cdot b\right)}^{2}\right)\\
          
          
          \end{array}
          \end{array}
          
          Derivation
          1. Split input into 2 regimes
          2. if a < 2.1e30

            1. Initial program 75.5%

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

              \[\leadsto \color{blue}{{b}^{2}} \]
            4. Step-by-step derivation
              1. unpow2N/A

                \[\leadsto \color{blue}{b \cdot b} \]
              2. lower-*.f6462.7

                \[\leadsto \color{blue}{b \cdot b} \]
            5. Applied rewrites62.7%

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

            if 2.1e30 < a

            1. Initial program 95.6%

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

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

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

                  \[\leadsto \color{blue}{{\left(a \cdot \sin \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} + {\left(b \cdot 1\right)}^{2}} \]
                2. lift-pow.f64N/A

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

                  \[\leadsto \color{blue}{\left(a \cdot \sin \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \left(a \cdot \sin \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)} + {\left(b \cdot 1\right)}^{2} \]
                4. lift-*.f64N/A

                  \[\leadsto \left(a \cdot \sin \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \color{blue}{\left(a \cdot \sin \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)} + {\left(b \cdot 1\right)}^{2} \]
                5. *-commutativeN/A

                  \[\leadsto \left(a \cdot \sin \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \color{blue}{\left(\sin \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right) \cdot a\right)} + {\left(b \cdot 1\right)}^{2} \]
                6. associate-*r*N/A

                  \[\leadsto \color{blue}{\left(\left(a \cdot \sin \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \sin \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right) \cdot a} + {\left(b \cdot 1\right)}^{2} \]
                7. lower-fma.f64N/A

                  \[\leadsto \color{blue}{\mathsf{fma}\left(\left(a \cdot \sin \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \sin \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right), a, {\left(b \cdot 1\right)}^{2}\right)} \]
              3. Applied rewrites95.6%

                \[\leadsto \color{blue}{\mathsf{fma}\left(\left(\sin \left(\mathsf{PI}\left(\right) \cdot \left(0.005555555555555556 \cdot angle\right)\right) \cdot a\right) \cdot \sin \left(\mathsf{PI}\left(\right) \cdot \left(0.005555555555555556 \cdot angle\right)\right), a, {\left(1 \cdot b\right)}^{2}\right)} \]
              4. Taylor expanded in angle around 0

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

                  \[\leadsto \mathsf{fma}\left(\frac{1}{32400} \cdot \left(a \cdot \color{blue}{\left({\mathsf{PI}\left(\right)}^{2} \cdot {angle}^{2}\right)}\right), a, {\left(1 \cdot b\right)}^{2}\right) \]
                2. associate-*r*N/A

                  \[\leadsto \mathsf{fma}\left(\frac{1}{32400} \cdot \color{blue}{\left(\left(a \cdot {\mathsf{PI}\left(\right)}^{2}\right) \cdot {angle}^{2}\right)}, a, {\left(1 \cdot b\right)}^{2}\right) \]
                3. associate-*l*N/A

                  \[\leadsto \mathsf{fma}\left(\color{blue}{\left(\frac{1}{32400} \cdot \left(a \cdot {\mathsf{PI}\left(\right)}^{2}\right)\right) \cdot {angle}^{2}}, a, {\left(1 \cdot b\right)}^{2}\right) \]
                4. unpow2N/A

                  \[\leadsto \mathsf{fma}\left(\left(\frac{1}{32400} \cdot \left(a \cdot {\mathsf{PI}\left(\right)}^{2}\right)\right) \cdot \color{blue}{\left(angle \cdot angle\right)}, a, {\left(1 \cdot b\right)}^{2}\right) \]
                5. associate-*r*N/A

                  \[\leadsto \mathsf{fma}\left(\color{blue}{\left(\left(\frac{1}{32400} \cdot \left(a \cdot {\mathsf{PI}\left(\right)}^{2}\right)\right) \cdot angle\right) \cdot angle}, a, {\left(1 \cdot b\right)}^{2}\right) \]
                6. lower-*.f64N/A

                  \[\leadsto \mathsf{fma}\left(\color{blue}{\left(\left(\frac{1}{32400} \cdot \left(a \cdot {\mathsf{PI}\left(\right)}^{2}\right)\right) \cdot angle\right) \cdot angle}, a, {\left(1 \cdot b\right)}^{2}\right) \]
                7. lower-*.f64N/A

                  \[\leadsto \mathsf{fma}\left(\color{blue}{\left(\left(\frac{1}{32400} \cdot \left(a \cdot {\mathsf{PI}\left(\right)}^{2}\right)\right) \cdot angle\right)} \cdot angle, a, {\left(1 \cdot b\right)}^{2}\right) \]
                8. *-commutativeN/A

                  \[\leadsto \mathsf{fma}\left(\left(\color{blue}{\left(\left(a \cdot {\mathsf{PI}\left(\right)}^{2}\right) \cdot \frac{1}{32400}\right)} \cdot angle\right) \cdot angle, a, {\left(1 \cdot b\right)}^{2}\right) \]
                9. lower-*.f64N/A

                  \[\leadsto \mathsf{fma}\left(\left(\color{blue}{\left(\left(a \cdot {\mathsf{PI}\left(\right)}^{2}\right) \cdot \frac{1}{32400}\right)} \cdot angle\right) \cdot angle, a, {\left(1 \cdot b\right)}^{2}\right) \]
                10. *-commutativeN/A

                  \[\leadsto \mathsf{fma}\left(\left(\left(\color{blue}{\left({\mathsf{PI}\left(\right)}^{2} \cdot a\right)} \cdot \frac{1}{32400}\right) \cdot angle\right) \cdot angle, a, {\left(1 \cdot b\right)}^{2}\right) \]
                11. lower-*.f64N/A

                  \[\leadsto \mathsf{fma}\left(\left(\left(\color{blue}{\left({\mathsf{PI}\left(\right)}^{2} \cdot a\right)} \cdot \frac{1}{32400}\right) \cdot angle\right) \cdot angle, a, {\left(1 \cdot b\right)}^{2}\right) \]
                12. unpow2N/A

                  \[\leadsto \mathsf{fma}\left(\left(\left(\left(\color{blue}{\left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right)} \cdot a\right) \cdot \frac{1}{32400}\right) \cdot angle\right) \cdot angle, a, {\left(1 \cdot b\right)}^{2}\right) \]
                13. lower-*.f64N/A

                  \[\leadsto \mathsf{fma}\left(\left(\left(\left(\color{blue}{\left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right)} \cdot a\right) \cdot \frac{1}{32400}\right) \cdot angle\right) \cdot angle, a, {\left(1 \cdot b\right)}^{2}\right) \]
                14. lower-PI.f64N/A

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

                  \[\leadsto \mathsf{fma}\left(\left(\left(\left(\left(\mathsf{PI}\left(\right) \cdot \color{blue}{\mathsf{PI}\left(\right)}\right) \cdot a\right) \cdot 3.08641975308642 \cdot 10^{-5}\right) \cdot angle\right) \cdot angle, a, {\left(1 \cdot b\right)}^{2}\right) \]
              6. Applied rewrites91.6%

                \[\leadsto \mathsf{fma}\left(\color{blue}{\left(\left(\left(\left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right) \cdot a\right) \cdot 3.08641975308642 \cdot 10^{-5}\right) \cdot angle\right) \cdot angle}, a, {\left(1 \cdot b\right)}^{2}\right) \]
            5. Recombined 2 regimes into one program.
            6. Final simplification69.2%

              \[\leadsto \begin{array}{l} \mathbf{if}\;a \leq 2.1 \cdot 10^{+30}:\\ \;\;\;\;b \cdot b\\ \mathbf{else}:\\ \;\;\;\;\mathsf{fma}\left(\left(\left(3.08641975308642 \cdot 10^{-5} \cdot \left(\left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right) \cdot a\right)\right) \cdot angle\right) \cdot angle, a, {\left(1 \cdot b\right)}^{2}\right)\\ \end{array} \]
            7. Add Preprocessing

            Alternative 6: 63.4% accurate, 9.1× speedup?

            \[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;a \leq 2.1 \cdot 10^{+30}:\\ \;\;\;\;b \cdot b\\ \mathbf{elif}\;a \leq 1.3 \cdot 10^{+144}:\\ \;\;\;\;\mathsf{fma}\left(\left(\left(\left(a \cdot a\right) \cdot 3.08641975308642 \cdot 10^{-5}\right) \cdot \mathsf{PI}\left(\right)\right) \cdot \mathsf{PI}\left(\right), angle \cdot angle, b \cdot b\right)\\ \mathbf{else}:\\ \;\;\;\;\left(\left(angle \cdot a\right) \cdot \left(\left(3.08641975308642 \cdot 10^{-5} \cdot a\right) \cdot angle\right)\right) \cdot \left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right)\\ \end{array} \end{array} \]
            (FPCore (a b angle)
             :precision binary64
             (if (<= a 2.1e+30)
               (* b b)
               (if (<= a 1.3e+144)
                 (fma
                  (* (* (* (* a a) 3.08641975308642e-5) (PI)) (PI))
                  (* angle angle)
                  (* b b))
                 (* (* (* angle a) (* (* 3.08641975308642e-5 a) angle)) (* (PI) (PI))))))
            \begin{array}{l}
            
            \\
            \begin{array}{l}
            \mathbf{if}\;a \leq 2.1 \cdot 10^{+30}:\\
            \;\;\;\;b \cdot b\\
            
            \mathbf{elif}\;a \leq 1.3 \cdot 10^{+144}:\\
            \;\;\;\;\mathsf{fma}\left(\left(\left(\left(a \cdot a\right) \cdot 3.08641975308642 \cdot 10^{-5}\right) \cdot \mathsf{PI}\left(\right)\right) \cdot \mathsf{PI}\left(\right), angle \cdot angle, b \cdot b\right)\\
            
            \mathbf{else}:\\
            \;\;\;\;\left(\left(angle \cdot a\right) \cdot \left(\left(3.08641975308642 \cdot 10^{-5} \cdot a\right) \cdot angle\right)\right) \cdot \left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right)\\
            
            
            \end{array}
            \end{array}
            
            Derivation
            1. Split input into 3 regimes
            2. if a < 2.1e30

              1. Initial program 75.5%

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

                \[\leadsto \color{blue}{{b}^{2}} \]
              4. Step-by-step derivation
                1. unpow2N/A

                  \[\leadsto \color{blue}{b \cdot b} \]
                2. lower-*.f6462.7

                  \[\leadsto \color{blue}{b \cdot b} \]
              5. Applied rewrites62.7%

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

              if 2.1e30 < a < 1.2999999999999999e144

              1. Initial program 86.4%

                \[{\left(a \cdot \sin \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} + {\left(b \cdot \cos \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\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({b}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right) + \frac{1}{32400} \cdot \left({a}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right)\right) + {b}^{2}} \]
              4. Step-by-step derivation
                1. *-commutativeN/A

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

                  \[\leadsto \color{blue}{\mathsf{fma}\left(\frac{-1}{32400} \cdot \left({b}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right) + \frac{1}{32400} \cdot \left({a}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right), {angle}^{2}, {b}^{2}\right)} \]
              5. Applied rewrites38.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}, a \cdot a, -3.08641975308642 \cdot 10^{-5} \cdot \left(b \cdot b\right)\right), angle \cdot angle, b \cdot b\right)} \]
              6. Taylor expanded in a around inf

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

                  \[\leadsto \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), \color{blue}{angle} \cdot angle, b \cdot b\right) \]

                if 1.2999999999999999e144 < a

                1. Initial program 99.7%

                  \[{\left(a \cdot \sin \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} + {\left(b \cdot \cos \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\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({b}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right) + \frac{1}{32400} \cdot \left({a}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right)\right) + {b}^{2}} \]
                4. Step-by-step derivation
                  1. *-commutativeN/A

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

                    \[\leadsto \color{blue}{\mathsf{fma}\left(\frac{-1}{32400} \cdot \left({b}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right) + \frac{1}{32400} \cdot \left({a}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right), {angle}^{2}, {b}^{2}\right)} \]
                5. Applied rewrites58.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}, a \cdot a, -3.08641975308642 \cdot 10^{-5} \cdot \left(b \cdot b\right)\right), angle \cdot angle, b \cdot b\right)} \]
                6. Taylor expanded in a around inf

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

                    \[\leadsto \left(3.08641975308642 \cdot 10^{-5} \cdot \left(a \cdot \left(\left(angle \cdot angle\right) \cdot a\right)\right)\right) \cdot \color{blue}{\left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right)} \]
                  2. Step-by-step derivation
                    1. Applied rewrites92.6%

                      \[\leadsto \left(\left(\left(3.08641975308642 \cdot 10^{-5} \cdot a\right) \cdot angle\right) \cdot \left(a \cdot angle\right)\right) \cdot \left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right) \]
                  3. Recombined 3 regimes into one program.
                  4. Final simplification68.5%

                    \[\leadsto \begin{array}{l} \mathbf{if}\;a \leq 2.1 \cdot 10^{+30}:\\ \;\;\;\;b \cdot b\\ \mathbf{elif}\;a \leq 1.3 \cdot 10^{+144}:\\ \;\;\;\;\mathsf{fma}\left(\left(\left(\left(a \cdot a\right) \cdot 3.08641975308642 \cdot 10^{-5}\right) \cdot \mathsf{PI}\left(\right)\right) \cdot \mathsf{PI}\left(\right), angle \cdot angle, b \cdot b\right)\\ \mathbf{else}:\\ \;\;\;\;\left(\left(angle \cdot a\right) \cdot \left(\left(3.08641975308642 \cdot 10^{-5} \cdot a\right) \cdot angle\right)\right) \cdot \left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right)\\ \end{array} \]
                  5. Add Preprocessing

                  Alternative 7: 62.4% accurate, 12.1× speedup?

                  \[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;a \leq 4.8 \cdot 10^{+139}:\\ \;\;\;\;b \cdot b\\ \mathbf{else}:\\ \;\;\;\;\left(\left(angle \cdot a\right) \cdot \left(\left(3.08641975308642 \cdot 10^{-5} \cdot a\right) \cdot angle\right)\right) \cdot \left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right)\\ \end{array} \end{array} \]
                  (FPCore (a b angle)
                   :precision binary64
                   (if (<= a 4.8e+139)
                     (* b b)
                     (* (* (* angle a) (* (* 3.08641975308642e-5 a) angle)) (* (PI) (PI)))))
                  \begin{array}{l}
                  
                  \\
                  \begin{array}{l}
                  \mathbf{if}\;a \leq 4.8 \cdot 10^{+139}:\\
                  \;\;\;\;b \cdot b\\
                  
                  \mathbf{else}:\\
                  \;\;\;\;\left(\left(angle \cdot a\right) \cdot \left(\left(3.08641975308642 \cdot 10^{-5} \cdot a\right) \cdot angle\right)\right) \cdot \left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right)\\
                  
                  
                  \end{array}
                  \end{array}
                  
                  Derivation
                  1. Split input into 2 regimes
                  2. if a < 4.80000000000000016e139

                    1. Initial program 76.4%

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

                      \[\leadsto \color{blue}{{b}^{2}} \]
                    4. Step-by-step derivation
                      1. unpow2N/A

                        \[\leadsto \color{blue}{b \cdot b} \]
                      2. lower-*.f6462.8

                        \[\leadsto \color{blue}{b \cdot b} \]
                    5. Applied rewrites62.8%

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

                    if 4.80000000000000016e139 < a

                    1. Initial program 99.7%

                      \[{\left(a \cdot \sin \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} + {\left(b \cdot \cos \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\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({b}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right) + \frac{1}{32400} \cdot \left({a}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right)\right) + {b}^{2}} \]
                    4. Step-by-step derivation
                      1. *-commutativeN/A

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

                        \[\leadsto \color{blue}{\mathsf{fma}\left(\frac{-1}{32400} \cdot \left({b}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right) + \frac{1}{32400} \cdot \left({a}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right), {angle}^{2}, {b}^{2}\right)} \]
                    5. Applied rewrites58.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}, a \cdot a, -3.08641975308642 \cdot 10^{-5} \cdot \left(b \cdot b\right)\right), angle \cdot angle, b \cdot b\right)} \]
                    6. Taylor expanded in a around inf

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

                        \[\leadsto \left(3.08641975308642 \cdot 10^{-5} \cdot \left(a \cdot \left(\left(angle \cdot angle\right) \cdot a\right)\right)\right) \cdot \color{blue}{\left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right)} \]
                      2. Step-by-step derivation
                        1. Applied rewrites92.6%

                          \[\leadsto \left(\left(\left(3.08641975308642 \cdot 10^{-5} \cdot a\right) \cdot angle\right) \cdot \left(a \cdot angle\right)\right) \cdot \left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right) \]
                      3. Recombined 2 regimes into one program.
                      4. Final simplification67.5%

                        \[\leadsto \begin{array}{l} \mathbf{if}\;a \leq 4.8 \cdot 10^{+139}:\\ \;\;\;\;b \cdot b\\ \mathbf{else}:\\ \;\;\;\;\left(\left(angle \cdot a\right) \cdot \left(\left(3.08641975308642 \cdot 10^{-5} \cdot a\right) \cdot angle\right)\right) \cdot \left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right)\\ \end{array} \]
                      5. Add Preprocessing

                      Alternative 8: 62.1% accurate, 12.1× speedup?

                      \[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;a \leq 4.8 \cdot 10^{+139}:\\ \;\;\;\;b \cdot b\\ \mathbf{else}:\\ \;\;\;\;\left(\left(\left(\left(angle \cdot a\right) \cdot angle\right) \cdot a\right) \cdot 3.08641975308642 \cdot 10^{-5}\right) \cdot \left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right)\\ \end{array} \end{array} \]
                      (FPCore (a b angle)
                       :precision binary64
                       (if (<= a 4.8e+139)
                         (* b b)
                         (* (* (* (* (* angle a) angle) a) 3.08641975308642e-5) (* (PI) (PI)))))
                      \begin{array}{l}
                      
                      \\
                      \begin{array}{l}
                      \mathbf{if}\;a \leq 4.8 \cdot 10^{+139}:\\
                      \;\;\;\;b \cdot b\\
                      
                      \mathbf{else}:\\
                      \;\;\;\;\left(\left(\left(\left(angle \cdot a\right) \cdot angle\right) \cdot a\right) \cdot 3.08641975308642 \cdot 10^{-5}\right) \cdot \left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right)\\
                      
                      
                      \end{array}
                      \end{array}
                      
                      Derivation
                      1. Split input into 2 regimes
                      2. if a < 4.80000000000000016e139

                        1. Initial program 76.4%

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

                          \[\leadsto \color{blue}{{b}^{2}} \]
                        4. Step-by-step derivation
                          1. unpow2N/A

                            \[\leadsto \color{blue}{b \cdot b} \]
                          2. lower-*.f6462.8

                            \[\leadsto \color{blue}{b \cdot b} \]
                        5. Applied rewrites62.8%

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

                        if 4.80000000000000016e139 < a

                        1. Initial program 99.7%

                          \[{\left(a \cdot \sin \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} + {\left(b \cdot \cos \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\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({b}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right) + \frac{1}{32400} \cdot \left({a}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right)\right) + {b}^{2}} \]
                        4. Step-by-step derivation
                          1. *-commutativeN/A

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

                            \[\leadsto \color{blue}{\mathsf{fma}\left(\frac{-1}{32400} \cdot \left({b}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right) + \frac{1}{32400} \cdot \left({a}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right), {angle}^{2}, {b}^{2}\right)} \]
                        5. Applied rewrites58.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}, a \cdot a, -3.08641975308642 \cdot 10^{-5} \cdot \left(b \cdot b\right)\right), angle \cdot angle, b \cdot b\right)} \]
                        6. Taylor expanded in a around inf

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

                            \[\leadsto \left(3.08641975308642 \cdot 10^{-5} \cdot \left(a \cdot \left(\left(angle \cdot angle\right) \cdot a\right)\right)\right) \cdot \color{blue}{\left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right)} \]
                          2. Step-by-step derivation
                            1. Applied rewrites92.5%

                              \[\leadsto \left(3.08641975308642 \cdot 10^{-5} \cdot \left(a \cdot \left(\left(a \cdot angle\right) \cdot angle\right)\right)\right) \cdot \left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right) \]
                          3. Recombined 2 regimes into one program.
                          4. Final simplification67.4%

                            \[\leadsto \begin{array}{l} \mathbf{if}\;a \leq 4.8 \cdot 10^{+139}:\\ \;\;\;\;b \cdot b\\ \mathbf{else}:\\ \;\;\;\;\left(\left(\left(\left(angle \cdot a\right) \cdot angle\right) \cdot a\right) \cdot 3.08641975308642 \cdot 10^{-5}\right) \cdot \left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right)\\ \end{array} \]
                          5. Add Preprocessing

                          Alternative 9: 61.3% accurate, 12.1× speedup?

                          \[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;a \leq 4.8 \cdot 10^{+139}:\\ \;\;\;\;b \cdot b\\ \mathbf{else}:\\ \;\;\;\;\left(\left(\left(\left(a \cdot a\right) \cdot 3.08641975308642 \cdot 10^{-5}\right) \cdot angle\right) \cdot angle\right) \cdot \left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right)\\ \end{array} \end{array} \]
                          (FPCore (a b angle)
                           :precision binary64
                           (if (<= a 4.8e+139)
                             (* b b)
                             (* (* (* (* (* a a) 3.08641975308642e-5) angle) angle) (* (PI) (PI)))))
                          \begin{array}{l}
                          
                          \\
                          \begin{array}{l}
                          \mathbf{if}\;a \leq 4.8 \cdot 10^{+139}:\\
                          \;\;\;\;b \cdot b\\
                          
                          \mathbf{else}:\\
                          \;\;\;\;\left(\left(\left(\left(a \cdot a\right) \cdot 3.08641975308642 \cdot 10^{-5}\right) \cdot angle\right) \cdot angle\right) \cdot \left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right)\\
                          
                          
                          \end{array}
                          \end{array}
                          
                          Derivation
                          1. Split input into 2 regimes
                          2. if a < 4.80000000000000016e139

                            1. Initial program 76.4%

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

                              \[\leadsto \color{blue}{{b}^{2}} \]
                            4. Step-by-step derivation
                              1. unpow2N/A

                                \[\leadsto \color{blue}{b \cdot b} \]
                              2. lower-*.f6462.8

                                \[\leadsto \color{blue}{b \cdot b} \]
                            5. Applied rewrites62.8%

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

                            if 4.80000000000000016e139 < a

                            1. Initial program 99.7%

                              \[{\left(a \cdot \sin \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} + {\left(b \cdot \cos \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\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({b}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right) + \frac{1}{32400} \cdot \left({a}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right)\right) + {b}^{2}} \]
                            4. Step-by-step derivation
                              1. *-commutativeN/A

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

                                \[\leadsto \color{blue}{\mathsf{fma}\left(\frac{-1}{32400} \cdot \left({b}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right) + \frac{1}{32400} \cdot \left({a}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right), {angle}^{2}, {b}^{2}\right)} \]
                            5. Applied rewrites58.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}, a \cdot a, -3.08641975308642 \cdot 10^{-5} \cdot \left(b \cdot b\right)\right), angle \cdot angle, b \cdot b\right)} \]
                            6. Taylor expanded in a around inf

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

                                \[\leadsto \left(3.08641975308642 \cdot 10^{-5} \cdot \left(a \cdot \left(\left(angle \cdot angle\right) \cdot a\right)\right)\right) \cdot \color{blue}{\left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right)} \]
                              2. Taylor expanded in a around 0

                                \[\leadsto \left(\frac{1}{32400} \cdot \left({a}^{2} \cdot {angle}^{2}\right)\right) \cdot \left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right) \]
                              3. Step-by-step derivation
                                1. Applied rewrites85.8%

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

                              Alternative 10: 60.2% accurate, 12.1× speedup?

                              \[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;a \leq 2.95 \cdot 10^{+140}:\\ \;\;\;\;b \cdot b\\ \mathbf{else}:\\ \;\;\;\;\left(\left(\left(a \cdot a\right) \cdot \mathsf{PI}\left(\right)\right) \cdot \mathsf{PI}\left(\right)\right) \cdot \left(\left(angle \cdot angle\right) \cdot 3.08641975308642 \cdot 10^{-5}\right)\\ \end{array} \end{array} \]
                              (FPCore (a b angle)
                               :precision binary64
                               (if (<= a 2.95e+140)
                                 (* b b)
                                 (* (* (* (* a a) (PI)) (PI)) (* (* angle angle) 3.08641975308642e-5))))
                              \begin{array}{l}
                              
                              \\
                              \begin{array}{l}
                              \mathbf{if}\;a \leq 2.95 \cdot 10^{+140}:\\
                              \;\;\;\;b \cdot b\\
                              
                              \mathbf{else}:\\
                              \;\;\;\;\left(\left(\left(a \cdot a\right) \cdot \mathsf{PI}\left(\right)\right) \cdot \mathsf{PI}\left(\right)\right) \cdot \left(\left(angle \cdot angle\right) \cdot 3.08641975308642 \cdot 10^{-5}\right)\\
                              
                              
                              \end{array}
                              \end{array}
                              
                              Derivation
                              1. Split input into 2 regimes
                              2. if a < 2.9500000000000001e140

                                1. Initial program 76.4%

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

                                  \[\leadsto \color{blue}{{b}^{2}} \]
                                4. Step-by-step derivation
                                  1. unpow2N/A

                                    \[\leadsto \color{blue}{b \cdot b} \]
                                  2. lower-*.f6462.8

                                    \[\leadsto \color{blue}{b \cdot b} \]
                                5. Applied rewrites62.8%

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

                                if 2.9500000000000001e140 < a

                                1. Initial program 99.7%

                                  \[{\left(a \cdot \sin \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} + {\left(b \cdot \cos \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\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({b}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right) + \frac{1}{32400} \cdot \left({a}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right)\right) + {b}^{2}} \]
                                4. Step-by-step derivation
                                  1. *-commutativeN/A

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

                                    \[\leadsto \color{blue}{\mathsf{fma}\left(\frac{-1}{32400} \cdot \left({b}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right) + \frac{1}{32400} \cdot \left({a}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right), {angle}^{2}, {b}^{2}\right)} \]
                                5. Applied rewrites58.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}, a \cdot a, -3.08641975308642 \cdot 10^{-5} \cdot \left(b \cdot b\right)\right), angle \cdot angle, b \cdot b\right)} \]
                                6. Taylor expanded in a around inf

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

                                    \[\leadsto \left(3.08641975308642 \cdot 10^{-5} \cdot \left(a \cdot \left(\left(angle \cdot angle\right) \cdot a\right)\right)\right) \cdot \color{blue}{\left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right)} \]
                                  2. Taylor expanded in a around 0

                                    \[\leadsto \frac{1}{32400} \cdot \left({a}^{2} \cdot \color{blue}{\left({angle}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right)}\right) \]
                                  3. Step-by-step derivation
                                    1. Applied rewrites80.5%

                                      \[\leadsto \left(\left(angle \cdot angle\right) \cdot 3.08641975308642 \cdot 10^{-5}\right) \cdot \left(\left(\left(a \cdot a\right) \cdot \mathsf{PI}\left(\right)\right) \cdot \color{blue}{\mathsf{PI}\left(\right)}\right) \]
                                  4. Recombined 2 regimes into one program.
                                  5. Final simplification65.6%

                                    \[\leadsto \begin{array}{l} \mathbf{if}\;a \leq 2.95 \cdot 10^{+140}:\\ \;\;\;\;b \cdot b\\ \mathbf{else}:\\ \;\;\;\;\left(\left(\left(a \cdot a\right) \cdot \mathsf{PI}\left(\right)\right) \cdot \mathsf{PI}\left(\right)\right) \cdot \left(\left(angle \cdot angle\right) \cdot 3.08641975308642 \cdot 10^{-5}\right)\\ \end{array} \]
                                  6. Add Preprocessing

                                  Alternative 11: 56.9% accurate, 74.7× speedup?

                                  \[\begin{array}{l} \\ b \cdot b \end{array} \]
                                  (FPCore (a b angle) :precision binary64 (* b b))
                                  double code(double a, double b, double angle) {
                                  	return b * b;
                                  }
                                  
                                  real(8) function code(a, b, angle)
                                      real(8), intent (in) :: a
                                      real(8), intent (in) :: b
                                      real(8), intent (in) :: angle
                                      code = b * b
                                  end function
                                  
                                  public static double code(double a, double b, double angle) {
                                  	return b * b;
                                  }
                                  
                                  def code(a, b, angle):
                                  	return b * b
                                  
                                  function code(a, b, angle)
                                  	return Float64(b * b)
                                  end
                                  
                                  function tmp = code(a, b, angle)
                                  	tmp = b * b;
                                  end
                                  
                                  code[a_, b_, angle_] := N[(b * b), $MachinePrecision]
                                  
                                  \begin{array}{l}
                                  
                                  \\
                                  b \cdot b
                                  \end{array}
                                  
                                  Derivation
                                  1. Initial program 80.0%

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

                                    \[\leadsto \color{blue}{{b}^{2}} \]
                                  4. Step-by-step derivation
                                    1. unpow2N/A

                                      \[\leadsto \color{blue}{b \cdot b} \]
                                    2. lower-*.f6457.6

                                      \[\leadsto \color{blue}{b \cdot b} \]
                                  5. Applied rewrites57.6%

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

                                  Reproduce

                                  ?
                                  herbie shell --seed 2024332 
                                  (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)))