b from scale-rotated-ellipse

Percentage Accurate: 0.1% → 37.4%
Time: 43.0s
Alternatives: 7
Speedup: 484.7×

Specification

?
\[\begin{array}{l} \\ \begin{array}{l} t_0 := \frac{angle}{180} \cdot \mathsf{PI}\left(\right)\\ t_1 := \sin t\_0\\ t_2 := \cos t\_0\\ t_3 := \frac{\frac{{\left(a \cdot t\_2\right)}^{2} + {\left(b \cdot t\_1\right)}^{2}}{y-scale}}{y-scale}\\ t_4 := \frac{\frac{{\left(a \cdot t\_1\right)}^{2} + {\left(b \cdot t\_2\right)}^{2}}{x-scale}}{x-scale}\\ t_5 := \left(b \cdot a\right) \cdot \left(b \cdot \left(-a\right)\right)\\ t_6 := \frac{4 \cdot t\_5}{{\left(x-scale \cdot y-scale\right)}^{2}}\\ \frac{-\sqrt{\left(\left(2 \cdot t\_6\right) \cdot t\_5\right) \cdot \left(\left(t\_4 + t\_3\right) - \sqrt{{\left(t\_4 - t\_3\right)}^{2} + {\left(\frac{\frac{\left(\left(2 \cdot \left({b}^{2} - {a}^{2}\right)\right) \cdot t\_1\right) \cdot t\_2}{x-scale}}{y-scale}\right)}^{2}}\right)}}{t\_6} \end{array} \end{array} \]
(FPCore (a b angle x-scale y-scale)
 :precision binary64
 (let* ((t_0 (* (/ angle 180.0) (PI)))
        (t_1 (sin t_0))
        (t_2 (cos t_0))
        (t_3
         (/ (/ (+ (pow (* a t_2) 2.0) (pow (* b t_1) 2.0)) y-scale) y-scale))
        (t_4
         (/ (/ (+ (pow (* a t_1) 2.0) (pow (* b t_2) 2.0)) x-scale) x-scale))
        (t_5 (* (* b a) (* b (- a))))
        (t_6 (/ (* 4.0 t_5) (pow (* x-scale y-scale) 2.0))))
   (/
    (-
     (sqrt
      (*
       (* (* 2.0 t_6) t_5)
       (-
        (+ t_4 t_3)
        (sqrt
         (+
          (pow (- t_4 t_3) 2.0)
          (pow
           (/
            (/ (* (* (* 2.0 (- (pow b 2.0) (pow a 2.0))) t_1) t_2) x-scale)
            y-scale)
           2.0)))))))
    t_6)))
\begin{array}{l}

\\
\begin{array}{l}
t_0 := \frac{angle}{180} \cdot \mathsf{PI}\left(\right)\\
t_1 := \sin t\_0\\
t_2 := \cos t\_0\\
t_3 := \frac{\frac{{\left(a \cdot t\_2\right)}^{2} + {\left(b \cdot t\_1\right)}^{2}}{y-scale}}{y-scale}\\
t_4 := \frac{\frac{{\left(a \cdot t\_1\right)}^{2} + {\left(b \cdot t\_2\right)}^{2}}{x-scale}}{x-scale}\\
t_5 := \left(b \cdot a\right) \cdot \left(b \cdot \left(-a\right)\right)\\
t_6 := \frac{4 \cdot t\_5}{{\left(x-scale \cdot y-scale\right)}^{2}}\\
\frac{-\sqrt{\left(\left(2 \cdot t\_6\right) \cdot t\_5\right) \cdot \left(\left(t\_4 + t\_3\right) - \sqrt{{\left(t\_4 - t\_3\right)}^{2} + {\left(\frac{\frac{\left(\left(2 \cdot \left({b}^{2} - {a}^{2}\right)\right) \cdot t\_1\right) \cdot t\_2}{x-scale}}{y-scale}\right)}^{2}}\right)}}{t\_6}
\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 7 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: 0.1% accurate, 1.0× speedup?

\[\begin{array}{l} \\ \begin{array}{l} t_0 := \frac{angle}{180} \cdot \mathsf{PI}\left(\right)\\ t_1 := \sin t\_0\\ t_2 := \cos t\_0\\ t_3 := \frac{\frac{{\left(a \cdot t\_2\right)}^{2} + {\left(b \cdot t\_1\right)}^{2}}{y-scale}}{y-scale}\\ t_4 := \frac{\frac{{\left(a \cdot t\_1\right)}^{2} + {\left(b \cdot t\_2\right)}^{2}}{x-scale}}{x-scale}\\ t_5 := \left(b \cdot a\right) \cdot \left(b \cdot \left(-a\right)\right)\\ t_6 := \frac{4 \cdot t\_5}{{\left(x-scale \cdot y-scale\right)}^{2}}\\ \frac{-\sqrt{\left(\left(2 \cdot t\_6\right) \cdot t\_5\right) \cdot \left(\left(t\_4 + t\_3\right) - \sqrt{{\left(t\_4 - t\_3\right)}^{2} + {\left(\frac{\frac{\left(\left(2 \cdot \left({b}^{2} - {a}^{2}\right)\right) \cdot t\_1\right) \cdot t\_2}{x-scale}}{y-scale}\right)}^{2}}\right)}}{t\_6} \end{array} \end{array} \]
(FPCore (a b angle x-scale y-scale)
 :precision binary64
 (let* ((t_0 (* (/ angle 180.0) (PI)))
        (t_1 (sin t_0))
        (t_2 (cos t_0))
        (t_3
         (/ (/ (+ (pow (* a t_2) 2.0) (pow (* b t_1) 2.0)) y-scale) y-scale))
        (t_4
         (/ (/ (+ (pow (* a t_1) 2.0) (pow (* b t_2) 2.0)) x-scale) x-scale))
        (t_5 (* (* b a) (* b (- a))))
        (t_6 (/ (* 4.0 t_5) (pow (* x-scale y-scale) 2.0))))
   (/
    (-
     (sqrt
      (*
       (* (* 2.0 t_6) t_5)
       (-
        (+ t_4 t_3)
        (sqrt
         (+
          (pow (- t_4 t_3) 2.0)
          (pow
           (/
            (/ (* (* (* 2.0 (- (pow b 2.0) (pow a 2.0))) t_1) t_2) x-scale)
            y-scale)
           2.0)))))))
    t_6)))
\begin{array}{l}

\\
\begin{array}{l}
t_0 := \frac{angle}{180} \cdot \mathsf{PI}\left(\right)\\
t_1 := \sin t\_0\\
t_2 := \cos t\_0\\
t_3 := \frac{\frac{{\left(a \cdot t\_2\right)}^{2} + {\left(b \cdot t\_1\right)}^{2}}{y-scale}}{y-scale}\\
t_4 := \frac{\frac{{\left(a \cdot t\_1\right)}^{2} + {\left(b \cdot t\_2\right)}^{2}}{x-scale}}{x-scale}\\
t_5 := \left(b \cdot a\right) \cdot \left(b \cdot \left(-a\right)\right)\\
t_6 := \frac{4 \cdot t\_5}{{\left(x-scale \cdot y-scale\right)}^{2}}\\
\frac{-\sqrt{\left(\left(2 \cdot t\_6\right) \cdot t\_5\right) \cdot \left(\left(t\_4 + t\_3\right) - \sqrt{{\left(t\_4 - t\_3\right)}^{2} + {\left(\frac{\frac{\left(\left(2 \cdot \left({b}^{2} - {a}^{2}\right)\right) \cdot t\_1\right) \cdot t\_2}{x-scale}}{y-scale}\right)}^{2}}\right)}}{t\_6}
\end{array}
\end{array}

Alternative 1: 37.4% accurate, 2.6× speedup?

\[\begin{array}{l} a_m = \left|a\right| \\ x-scale_m = \left|x-scale\right| \\ y-scale_m = \left|y-scale\right| \\ \begin{array}{l} t_0 := angle \cdot \mathsf{PI}\left(\right)\\ t_1 := 0.005555555555555556 \cdot t\_0\\ t_2 := a\_m \cdot \left(x-scale\_m \cdot \left(y-scale\_m \cdot \sqrt{8}\right)\right)\\ t_3 := \sin t\_1\\ t_4 := \cos t\_1\\ t_5 := \frac{{\left(t\_4 \cdot t\_3\right)}^{2}}{x-scale\_m \cdot x-scale\_m}\\ t_6 := \frac{{t\_4}^{2}}{y-scale\_m \cdot y-scale\_m}\\ \mathbf{if}\;y-scale\_m \leq 1.8 \cdot 10^{-151}:\\ \;\;\;\;0.25 \cdot \left(t\_2 \cdot \sqrt{\frac{{t\_3}^{2}}{x-scale\_m \cdot x-scale\_m} - 0.5 \cdot \frac{\mathsf{fma}\left(-2, t\_5, 4 \cdot t\_5\right)}{1 + -3.08641975308642 \cdot 10^{-5} \cdot {t\_0}^{2}}}\right)\\ \mathbf{elif}\;y-scale\_m \leq 1.55 \cdot 10^{-73}:\\ \;\;\;\;0.25 \cdot \left(t\_2 \cdot \sqrt{t\_6 - t\_6}\right)\\ \mathbf{else}:\\ \;\;\;\;a\_m \cdot x-scale\_m\\ \end{array} \end{array} \]
a_m = (fabs.f64 a)
x-scale_m = (fabs.f64 x-scale)
y-scale_m = (fabs.f64 y-scale)
(FPCore (a_m b angle x-scale_m y-scale_m)
 :precision binary64
 (let* ((t_0 (* angle (PI)))
        (t_1 (* 0.005555555555555556 t_0))
        (t_2 (* a_m (* x-scale_m (* y-scale_m (sqrt 8.0)))))
        (t_3 (sin t_1))
        (t_4 (cos t_1))
        (t_5 (/ (pow (* t_4 t_3) 2.0) (* x-scale_m x-scale_m)))
        (t_6 (/ (pow t_4 2.0) (* y-scale_m y-scale_m))))
   (if (<= y-scale_m 1.8e-151)
     (*
      0.25
      (*
       t_2
       (sqrt
        (-
         (/ (pow t_3 2.0) (* x-scale_m x-scale_m))
         (*
          0.5
          (/
           (fma -2.0 t_5 (* 4.0 t_5))
           (+ 1.0 (* -3.08641975308642e-5 (pow t_0 2.0)))))))))
     (if (<= y-scale_m 1.55e-73)
       (* 0.25 (* t_2 (sqrt (- t_6 t_6))))
       (* a_m x-scale_m)))))
\begin{array}{l}
a_m = \left|a\right|
\\
x-scale_m = \left|x-scale\right|
\\
y-scale_m = \left|y-scale\right|

\\
\begin{array}{l}
t_0 := angle \cdot \mathsf{PI}\left(\right)\\
t_1 := 0.005555555555555556 \cdot t\_0\\
t_2 := a\_m \cdot \left(x-scale\_m \cdot \left(y-scale\_m \cdot \sqrt{8}\right)\right)\\
t_3 := \sin t\_1\\
t_4 := \cos t\_1\\
t_5 := \frac{{\left(t\_4 \cdot t\_3\right)}^{2}}{x-scale\_m \cdot x-scale\_m}\\
t_6 := \frac{{t\_4}^{2}}{y-scale\_m \cdot y-scale\_m}\\
\mathbf{if}\;y-scale\_m \leq 1.8 \cdot 10^{-151}:\\
\;\;\;\;0.25 \cdot \left(t\_2 \cdot \sqrt{\frac{{t\_3}^{2}}{x-scale\_m \cdot x-scale\_m} - 0.5 \cdot \frac{\mathsf{fma}\left(-2, t\_5, 4 \cdot t\_5\right)}{1 + -3.08641975308642 \cdot 10^{-5} \cdot {t\_0}^{2}}}\right)\\

\mathbf{elif}\;y-scale\_m \leq 1.55 \cdot 10^{-73}:\\
\;\;\;\;0.25 \cdot \left(t\_2 \cdot \sqrt{t\_6 - t\_6}\right)\\

\mathbf{else}:\\
\;\;\;\;a\_m \cdot x-scale\_m\\


\end{array}
\end{array}
Derivation
  1. Split input into 3 regimes
  2. if y-scale < 1.80000000000000016e-151

    1. Initial program 0.0%

      \[\frac{-\sqrt{\left(\left(2 \cdot \frac{4 \cdot \left(\left(b \cdot a\right) \cdot \left(b \cdot \left(-a\right)\right)\right)}{{\left(x-scale \cdot y-scale\right)}^{2}}\right) \cdot \left(\left(b \cdot a\right) \cdot \left(b \cdot \left(-a\right)\right)\right)\right) \cdot \left(\left(\frac{\frac{{\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}}{x-scale}}{x-scale} + \frac{\frac{{\left(a \cdot \cos \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} + {\left(b \cdot \sin \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale}}{y-scale}\right) - \sqrt{{\left(\frac{\frac{{\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}}{x-scale}}{x-scale} - \frac{\frac{{\left(a \cdot \cos \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} + {\left(b \cdot \sin \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale}}{y-scale}\right)}^{2} + {\left(\frac{\frac{\left(\left(2 \cdot \left({b}^{2} - {a}^{2}\right)\right) \cdot \sin \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \cos \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)}{x-scale}}{y-scale}\right)}^{2}}\right)}}{\frac{4 \cdot \left(\left(b \cdot a\right) \cdot \left(b \cdot \left(-a\right)\right)\right)}{{\left(x-scale \cdot y-scale\right)}^{2}}} \]
    2. Add Preprocessing
    3. Taylor expanded in a around inf

      \[\leadsto \color{blue}{\frac{1}{4} \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\left(\frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{y-scale}^{2}} + \frac{{\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{x-scale}^{2}}\right) - \sqrt{4 \cdot \frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} \cdot {\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{x-scale}^{2} \cdot {y-scale}^{2}} + {\left(\frac{{\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{x-scale}^{2}} - \frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{y-scale}^{2}}\right)}^{2}}}\right)} \]
    4. Applied rewrites2.1%

      \[\leadsto \color{blue}{0.25 \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\left(\frac{{\cos \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale \cdot y-scale} + \frac{{\sin \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{x-scale \cdot x-scale}\right) - \sqrt{\mathsf{fma}\left(4, \frac{{\left(\cos \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \sin \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)\right)}^{2}}{{\left(x-scale \cdot y-scale\right)}^{2}}, {\left(\frac{{\sin \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{x-scale \cdot x-scale} - \frac{{\cos \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale \cdot y-scale}\right)}^{2}\right)}}\right)} \]
    5. Taylor expanded in y-scale around 0

      \[\leadsto \frac{1}{4} \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\frac{{\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{x-scale}^{2}} - \frac{1}{2} \cdot \frac{-2 \cdot \frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} \cdot {\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{x-scale}^{2}} + 4 \cdot \frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} \cdot {\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{x-scale}^{2}}}{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}}\right) \]
    6. Step-by-step derivation
      1. Applied rewrites21.1%

        \[\leadsto 0.25 \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\frac{{\sin \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{x-scale \cdot x-scale} - 0.5 \cdot \frac{\mathsf{fma}\left(-2, \frac{{\left(\cos \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \sin \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)\right)}^{2}}{x-scale \cdot x-scale}, 4 \cdot \frac{{\left(\cos \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \sin \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)\right)}^{2}}{x-scale \cdot x-scale}\right)}{{\cos \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}}\right) \]
      2. Taylor expanded in angle around 0

        \[\leadsto \frac{1}{4} \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\frac{{\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{x-scale \cdot x-scale} - \frac{1}{2} \cdot \frac{\mathsf{fma}\left(-2, \frac{{\left(\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)\right)}^{2}}{x-scale \cdot x-scale}, 4 \cdot \frac{{\left(\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)\right)}^{2}}{x-scale \cdot x-scale}\right)}{1 + \frac{-1}{32400} \cdot \left({angle}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right)}}\right) \]
      3. Step-by-step derivation
        1. lower-+.f64N/A

          \[\leadsto \frac{1}{4} \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\frac{{\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{x-scale \cdot x-scale} - \frac{1}{2} \cdot \frac{\mathsf{fma}\left(-2, \frac{{\left(\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)\right)}^{2}}{x-scale \cdot x-scale}, 4 \cdot \frac{{\left(\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)\right)}^{2}}{x-scale \cdot x-scale}\right)}{1 + \frac{-1}{32400} \cdot \left({angle}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right)}}\right) \]
        2. lower-*.f64N/A

          \[\leadsto \frac{1}{4} \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\frac{{\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{x-scale \cdot x-scale} - \frac{1}{2} \cdot \frac{\mathsf{fma}\left(-2, \frac{{\left(\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)\right)}^{2}}{x-scale \cdot x-scale}, 4 \cdot \frac{{\left(\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)\right)}^{2}}{x-scale \cdot x-scale}\right)}{1 + \frac{-1}{32400} \cdot \left({angle}^{2} \cdot {\mathsf{PI}\left(\right)}^{2}\right)}}\right) \]
        3. pow-prod-downN/A

          \[\leadsto \frac{1}{4} \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\frac{{\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{x-scale \cdot x-scale} - \frac{1}{2} \cdot \frac{\mathsf{fma}\left(-2, \frac{{\left(\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)\right)}^{2}}{x-scale \cdot x-scale}, 4 \cdot \frac{{\left(\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)\right)}^{2}}{x-scale \cdot x-scale}\right)}{1 + \frac{-1}{32400} \cdot {\left(angle \cdot \mathsf{PI}\left(\right)\right)}^{2}}}\right) \]
        4. lower-pow.f64N/A

          \[\leadsto \frac{1}{4} \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\frac{{\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{x-scale \cdot x-scale} - \frac{1}{2} \cdot \frac{\mathsf{fma}\left(-2, \frac{{\left(\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)\right)}^{2}}{x-scale \cdot x-scale}, 4 \cdot \frac{{\left(\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)\right)}^{2}}{x-scale \cdot x-scale}\right)}{1 + \frac{-1}{32400} \cdot {\left(angle \cdot \mathsf{PI}\left(\right)\right)}^{2}}}\right) \]
        5. lift-*.f64N/A

          \[\leadsto \frac{1}{4} \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\frac{{\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{x-scale \cdot x-scale} - \frac{1}{2} \cdot \frac{\mathsf{fma}\left(-2, \frac{{\left(\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)\right)}^{2}}{x-scale \cdot x-scale}, 4 \cdot \frac{{\left(\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)\right)}^{2}}{x-scale \cdot x-scale}\right)}{1 + \frac{-1}{32400} \cdot {\left(angle \cdot \mathsf{PI}\left(\right)\right)}^{2}}}\right) \]
        6. lift-PI.f6420.7

          \[\leadsto 0.25 \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\frac{{\sin \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{x-scale \cdot x-scale} - 0.5 \cdot \frac{\mathsf{fma}\left(-2, \frac{{\left(\cos \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \sin \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)\right)}^{2}}{x-scale \cdot x-scale}, 4 \cdot \frac{{\left(\cos \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \sin \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)\right)}^{2}}{x-scale \cdot x-scale}\right)}{1 + -3.08641975308642 \cdot 10^{-5} \cdot {\left(angle \cdot \mathsf{PI}\left(\right)\right)}^{2}}}\right) \]
      4. Applied rewrites20.7%

        \[\leadsto 0.25 \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\frac{{\sin \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{x-scale \cdot x-scale} - 0.5 \cdot \frac{\mathsf{fma}\left(-2, \frac{{\left(\cos \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \sin \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)\right)}^{2}}{x-scale \cdot x-scale}, 4 \cdot \frac{{\left(\cos \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \sin \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)\right)}^{2}}{x-scale \cdot x-scale}\right)}{1 + -3.08641975308642 \cdot 10^{-5} \cdot {\left(angle \cdot \mathsf{PI}\left(\right)\right)}^{2}}}\right) \]

      if 1.80000000000000016e-151 < y-scale < 1.54999999999999985e-73

      1. Initial program 0.0%

        \[\frac{-\sqrt{\left(\left(2 \cdot \frac{4 \cdot \left(\left(b \cdot a\right) \cdot \left(b \cdot \left(-a\right)\right)\right)}{{\left(x-scale \cdot y-scale\right)}^{2}}\right) \cdot \left(\left(b \cdot a\right) \cdot \left(b \cdot \left(-a\right)\right)\right)\right) \cdot \left(\left(\frac{\frac{{\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}}{x-scale}}{x-scale} + \frac{\frac{{\left(a \cdot \cos \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} + {\left(b \cdot \sin \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale}}{y-scale}\right) - \sqrt{{\left(\frac{\frac{{\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}}{x-scale}}{x-scale} - \frac{\frac{{\left(a \cdot \cos \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} + {\left(b \cdot \sin \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale}}{y-scale}\right)}^{2} + {\left(\frac{\frac{\left(\left(2 \cdot \left({b}^{2} - {a}^{2}\right)\right) \cdot \sin \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \cos \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)}{x-scale}}{y-scale}\right)}^{2}}\right)}}{\frac{4 \cdot \left(\left(b \cdot a\right) \cdot \left(b \cdot \left(-a\right)\right)\right)}{{\left(x-scale \cdot y-scale\right)}^{2}}} \]
      2. Add Preprocessing
      3. Taylor expanded in a around inf

        \[\leadsto \color{blue}{\frac{1}{4} \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\left(\frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{y-scale}^{2}} + \frac{{\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{x-scale}^{2}}\right) - \sqrt{4 \cdot \frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} \cdot {\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{x-scale}^{2} \cdot {y-scale}^{2}} + {\left(\frac{{\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{x-scale}^{2}} - \frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{y-scale}^{2}}\right)}^{2}}}\right)} \]
      4. Applied rewrites0.0%

        \[\leadsto \color{blue}{0.25 \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\left(\frac{{\cos \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale \cdot y-scale} + \frac{{\sin \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{x-scale \cdot x-scale}\right) - \sqrt{\mathsf{fma}\left(4, \frac{{\left(\cos \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \sin \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)\right)}^{2}}{{\left(x-scale \cdot y-scale\right)}^{2}}, {\left(\frac{{\sin \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{x-scale \cdot x-scale} - \frac{{\cos \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale \cdot y-scale}\right)}^{2}\right)}}\right)} \]
      5. Taylor expanded in x-scale around inf

        \[\leadsto \frac{1}{4} \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\left(\frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale \cdot y-scale} + \frac{{\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{x-scale \cdot x-scale}\right) - \frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{y-scale}^{2}}}\right) \]
      6. Step-by-step derivation
        1. lift-*.f64N/A

          \[\leadsto \frac{1}{4} \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\left(\frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale \cdot y-scale} + \frac{{\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{x-scale \cdot x-scale}\right) - \frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{y-scale}^{2}}}\right) \]
        2. lift-PI.f64N/A

          \[\leadsto \frac{1}{4} \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\left(\frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale \cdot y-scale} + \frac{{\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{x-scale \cdot x-scale}\right) - \frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{y-scale}^{2}}}\right) \]
        3. lift-*.f64N/A

          \[\leadsto \frac{1}{4} \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\left(\frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale \cdot y-scale} + \frac{{\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{x-scale \cdot x-scale}\right) - \frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{y-scale}^{2}}}\right) \]
        4. lift-cos.f64N/A

          \[\leadsto \frac{1}{4} \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\left(\frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale \cdot y-scale} + \frac{{\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{x-scale \cdot x-scale}\right) - \frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{y-scale}^{2}}}\right) \]
        5. lift-pow.f64N/A

          \[\leadsto \frac{1}{4} \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\left(\frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale \cdot y-scale} + \frac{{\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{x-scale \cdot x-scale}\right) - \frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{y-scale}^{2}}}\right) \]
        6. pow2N/A

          \[\leadsto \frac{1}{4} \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\left(\frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale \cdot y-scale} + \frac{{\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{x-scale \cdot x-scale}\right) - \frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale \cdot y-scale}}\right) \]
        7. lift-/.f64N/A

          \[\leadsto \frac{1}{4} \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\left(\frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale \cdot y-scale} + \frac{{\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{x-scale \cdot x-scale}\right) - \frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale \cdot y-scale}}\right) \]
        8. lift-*.f6418.8

          \[\leadsto 0.25 \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\left(\frac{{\cos \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale \cdot y-scale} + \frac{{\sin \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{x-scale \cdot x-scale}\right) - \frac{{\cos \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale \cdot y-scale}}\right) \]
      7. Applied rewrites18.8%

        \[\leadsto 0.25 \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\left(\frac{{\cos \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale \cdot y-scale} + \frac{{\sin \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{x-scale \cdot x-scale}\right) - \frac{{\cos \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale \cdot y-scale}}\right) \]
      8. Taylor expanded in x-scale around inf

        \[\leadsto \frac{1}{4} \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{y-scale}^{2}} - \frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale \cdot y-scale}}\right) \]
      9. Step-by-step derivation
        1. lift-*.f64N/A

          \[\leadsto \frac{1}{4} \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{y-scale}^{2}} - \frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale \cdot y-scale}}\right) \]
        2. lift-PI.f64N/A

          \[\leadsto \frac{1}{4} \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{y-scale}^{2}} - \frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale \cdot y-scale}}\right) \]
        3. lift-*.f64N/A

          \[\leadsto \frac{1}{4} \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{y-scale}^{2}} - \frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale \cdot y-scale}}\right) \]
        4. lift-cos.f64N/A

          \[\leadsto \frac{1}{4} \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{y-scale}^{2}} - \frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale \cdot y-scale}}\right) \]
        5. lift-pow.f64N/A

          \[\leadsto \frac{1}{4} \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{y-scale}^{2}} - \frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale \cdot y-scale}}\right) \]
        6. pow2N/A

          \[\leadsto \frac{1}{4} \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale \cdot y-scale} - \frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale \cdot y-scale}}\right) \]
        7. lift-/.f64N/A

          \[\leadsto \frac{1}{4} \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale \cdot y-scale} - \frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale \cdot y-scale}}\right) \]
        8. lift-*.f6444.8

          \[\leadsto 0.25 \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\frac{{\cos \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale \cdot y-scale} - \frac{{\cos \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale \cdot y-scale}}\right) \]
      10. Applied rewrites44.8%

        \[\leadsto 0.25 \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\frac{{\cos \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale \cdot y-scale} - \frac{{\cos \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale \cdot y-scale}}\right) \]

      if 1.54999999999999985e-73 < y-scale

      1. Initial program 0.1%

        \[\frac{-\sqrt{\left(\left(2 \cdot \frac{4 \cdot \left(\left(b \cdot a\right) \cdot \left(b \cdot \left(-a\right)\right)\right)}{{\left(x-scale \cdot y-scale\right)}^{2}}\right) \cdot \left(\left(b \cdot a\right) \cdot \left(b \cdot \left(-a\right)\right)\right)\right) \cdot \left(\left(\frac{\frac{{\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}}{x-scale}}{x-scale} + \frac{\frac{{\left(a \cdot \cos \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} + {\left(b \cdot \sin \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale}}{y-scale}\right) - \sqrt{{\left(\frac{\frac{{\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}}{x-scale}}{x-scale} - \frac{\frac{{\left(a \cdot \cos \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} + {\left(b \cdot \sin \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale}}{y-scale}\right)}^{2} + {\left(\frac{\frac{\left(\left(2 \cdot \left({b}^{2} - {a}^{2}\right)\right) \cdot \sin \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \cos \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)}{x-scale}}{y-scale}\right)}^{2}}\right)}}{\frac{4 \cdot \left(\left(b \cdot a\right) \cdot \left(b \cdot \left(-a\right)\right)\right)}{{\left(x-scale \cdot y-scale\right)}^{2}}} \]
      2. Add Preprocessing
      3. Taylor expanded in angle around 0

        \[\leadsto \color{blue}{\frac{1}{4} \cdot \left(a \cdot \left(x-scale \cdot \left(\sqrt{2} \cdot \sqrt{8}\right)\right)\right)} \]
      4. Step-by-step derivation
        1. Applied rewrites23.9%

          \[\leadsto \color{blue}{0.25 \cdot \left(a \cdot \left(x-scale \cdot 4\right)\right)} \]
        2. Taylor expanded in a around 0

          \[\leadsto a \cdot \color{blue}{x-scale} \]
        3. Step-by-step derivation
          1. lower-*.f6423.9

            \[\leadsto a \cdot x-scale \]
        4. Applied rewrites23.9%

          \[\leadsto a \cdot \color{blue}{x-scale} \]
      5. Recombined 3 regimes into one program.
      6. Add Preprocessing

      Alternative 2: 36.0% accurate, 2.2× speedup?

      \[\begin{array}{l} a_m = \left|a\right| \\ x-scale_m = \left|x-scale\right| \\ y-scale_m = \left|y-scale\right| \\ \begin{array}{l} t_0 := 0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\\ t_1 := \sin t\_0\\ t_2 := {t\_1}^{2}\\ t_3 := \cos t\_0\\ t_4 := {t\_3}^{2}\\ \mathbf{if}\;y-scale\_m \leq 1.08 \cdot 10^{-120}:\\ \;\;\;\;0.25 \cdot \left(\left(a\_m \cdot \left(x-scale\_m \cdot \left(y-scale\_m \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\frac{t\_2}{x-scale\_m \cdot x-scale\_m} - 0.5 \cdot \frac{\mathsf{fma}\left(-2, \frac{{\left(t\_3 \cdot t\_1\right)}^{2}}{x-scale\_m \cdot x-scale\_m}, 4 \cdot \frac{t\_4 \cdot t\_2}{x-scale\_m \cdot x-scale\_m}\right)}{t\_4}}\right)\\ \mathbf{else}:\\ \;\;\;\;a\_m \cdot x-scale\_m\\ \end{array} \end{array} \]
      a_m = (fabs.f64 a)
      x-scale_m = (fabs.f64 x-scale)
      y-scale_m = (fabs.f64 y-scale)
      (FPCore (a_m b angle x-scale_m y-scale_m)
       :precision binary64
       (let* ((t_0 (* 0.005555555555555556 (* angle (PI))))
              (t_1 (sin t_0))
              (t_2 (pow t_1 2.0))
              (t_3 (cos t_0))
              (t_4 (pow t_3 2.0)))
         (if (<= y-scale_m 1.08e-120)
           (*
            0.25
            (*
             (* a_m (* x-scale_m (* y-scale_m (sqrt 8.0))))
             (sqrt
              (-
               (/ t_2 (* x-scale_m x-scale_m))
               (*
                0.5
                (/
                 (fma
                  -2.0
                  (/ (pow (* t_3 t_1) 2.0) (* x-scale_m x-scale_m))
                  (* 4.0 (/ (* t_4 t_2) (* x-scale_m x-scale_m))))
                 t_4))))))
           (* a_m x-scale_m))))
      \begin{array}{l}
      a_m = \left|a\right|
      \\
      x-scale_m = \left|x-scale\right|
      \\
      y-scale_m = \left|y-scale\right|
      
      \\
      \begin{array}{l}
      t_0 := 0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\\
      t_1 := \sin t\_0\\
      t_2 := {t\_1}^{2}\\
      t_3 := \cos t\_0\\
      t_4 := {t\_3}^{2}\\
      \mathbf{if}\;y-scale\_m \leq 1.08 \cdot 10^{-120}:\\
      \;\;\;\;0.25 \cdot \left(\left(a\_m \cdot \left(x-scale\_m \cdot \left(y-scale\_m \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\frac{t\_2}{x-scale\_m \cdot x-scale\_m} - 0.5 \cdot \frac{\mathsf{fma}\left(-2, \frac{{\left(t\_3 \cdot t\_1\right)}^{2}}{x-scale\_m \cdot x-scale\_m}, 4 \cdot \frac{t\_4 \cdot t\_2}{x-scale\_m \cdot x-scale\_m}\right)}{t\_4}}\right)\\
      
      \mathbf{else}:\\
      \;\;\;\;a\_m \cdot x-scale\_m\\
      
      
      \end{array}
      \end{array}
      
      Derivation
      1. Split input into 2 regimes
      2. if y-scale < 1.0800000000000001e-120

        1. Initial program 0.0%

          \[\frac{-\sqrt{\left(\left(2 \cdot \frac{4 \cdot \left(\left(b \cdot a\right) \cdot \left(b \cdot \left(-a\right)\right)\right)}{{\left(x-scale \cdot y-scale\right)}^{2}}\right) \cdot \left(\left(b \cdot a\right) \cdot \left(b \cdot \left(-a\right)\right)\right)\right) \cdot \left(\left(\frac{\frac{{\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}}{x-scale}}{x-scale} + \frac{\frac{{\left(a \cdot \cos \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} + {\left(b \cdot \sin \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale}}{y-scale}\right) - \sqrt{{\left(\frac{\frac{{\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}}{x-scale}}{x-scale} - \frac{\frac{{\left(a \cdot \cos \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} + {\left(b \cdot \sin \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale}}{y-scale}\right)}^{2} + {\left(\frac{\frac{\left(\left(2 \cdot \left({b}^{2} - {a}^{2}\right)\right) \cdot \sin \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \cos \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)}{x-scale}}{y-scale}\right)}^{2}}\right)}}{\frac{4 \cdot \left(\left(b \cdot a\right) \cdot \left(b \cdot \left(-a\right)\right)\right)}{{\left(x-scale \cdot y-scale\right)}^{2}}} \]
        2. Add Preprocessing
        3. Taylor expanded in a around inf

          \[\leadsto \color{blue}{\frac{1}{4} \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\left(\frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{y-scale}^{2}} + \frac{{\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{x-scale}^{2}}\right) - \sqrt{4 \cdot \frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} \cdot {\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{x-scale}^{2} \cdot {y-scale}^{2}} + {\left(\frac{{\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{x-scale}^{2}} - \frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{y-scale}^{2}}\right)}^{2}}}\right)} \]
        4. Applied rewrites2.1%

          \[\leadsto \color{blue}{0.25 \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\left(\frac{{\cos \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale \cdot y-scale} + \frac{{\sin \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{x-scale \cdot x-scale}\right) - \sqrt{\mathsf{fma}\left(4, \frac{{\left(\cos \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \sin \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)\right)}^{2}}{{\left(x-scale \cdot y-scale\right)}^{2}}, {\left(\frac{{\sin \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{x-scale \cdot x-scale} - \frac{{\cos \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale \cdot y-scale}\right)}^{2}\right)}}\right)} \]
        5. Taylor expanded in y-scale around 0

          \[\leadsto \frac{1}{4} \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\frac{{\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{x-scale}^{2}} - \frac{1}{2} \cdot \frac{-2 \cdot \frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} \cdot {\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{x-scale}^{2}} + 4 \cdot \frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} \cdot {\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{x-scale}^{2}}}{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}}\right) \]
        6. Step-by-step derivation
          1. Applied rewrites20.5%

            \[\leadsto 0.25 \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\frac{{\sin \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{x-scale \cdot x-scale} - 0.5 \cdot \frac{\mathsf{fma}\left(-2, \frac{{\left(\cos \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \sin \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)\right)}^{2}}{x-scale \cdot x-scale}, 4 \cdot \frac{{\left(\cos \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \sin \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)\right)}^{2}}{x-scale \cdot x-scale}\right)}{{\cos \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}}\right) \]
          2. Step-by-step derivation
            1. lift-pow.f64N/A

              \[\leadsto \frac{1}{4} \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\frac{{\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{x-scale \cdot x-scale} - \frac{1}{2} \cdot \frac{\mathsf{fma}\left(-2, \frac{{\left(\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)\right)}^{2}}{x-scale \cdot x-scale}, 4 \cdot \frac{{\left(\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)\right)}^{2}}{x-scale \cdot x-scale}\right)}{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}}\right) \]
            2. lift-*.f64N/A

              \[\leadsto \frac{1}{4} \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\frac{{\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{x-scale \cdot x-scale} - \frac{1}{2} \cdot \frac{\mathsf{fma}\left(-2, \frac{{\left(\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)\right)}^{2}}{x-scale \cdot x-scale}, 4 \cdot \frac{{\left(\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)\right)}^{2}}{x-scale \cdot x-scale}\right)}{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}}\right) \]
            3. lift-cos.f64N/A

              \[\leadsto \frac{1}{4} \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\frac{{\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{x-scale \cdot x-scale} - \frac{1}{2} \cdot \frac{\mathsf{fma}\left(-2, \frac{{\left(\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)\right)}^{2}}{x-scale \cdot x-scale}, 4 \cdot \frac{{\left(\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)\right)}^{2}}{x-scale \cdot x-scale}\right)}{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}}\right) \]
            4. lift-*.f64N/A

              \[\leadsto \frac{1}{4} \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\frac{{\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{x-scale \cdot x-scale} - \frac{1}{2} \cdot \frac{\mathsf{fma}\left(-2, \frac{{\left(\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)\right)}^{2}}{x-scale \cdot x-scale}, 4 \cdot \frac{{\left(\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)\right)}^{2}}{x-scale \cdot x-scale}\right)}{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}}\right) \]
            5. lift-PI.f64N/A

              \[\leadsto \frac{1}{4} \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\frac{{\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{x-scale \cdot x-scale} - \frac{1}{2} \cdot \frac{\mathsf{fma}\left(-2, \frac{{\left(\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)\right)}^{2}}{x-scale \cdot x-scale}, 4 \cdot \frac{{\left(\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)\right)}^{2}}{x-scale \cdot x-scale}\right)}{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}}\right) \]
            6. lift-*.f64N/A

              \[\leadsto \frac{1}{4} \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\frac{{\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{x-scale \cdot x-scale} - \frac{1}{2} \cdot \frac{\mathsf{fma}\left(-2, \frac{{\left(\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)\right)}^{2}}{x-scale \cdot x-scale}, 4 \cdot \frac{{\left(\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)\right)}^{2}}{x-scale \cdot x-scale}\right)}{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}}\right) \]
            7. lift-sin.f64N/A

              \[\leadsto \frac{1}{4} \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\frac{{\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{x-scale \cdot x-scale} - \frac{1}{2} \cdot \frac{\mathsf{fma}\left(-2, \frac{{\left(\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)\right)}^{2}}{x-scale \cdot x-scale}, 4 \cdot \frac{{\left(\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)\right)}^{2}}{x-scale \cdot x-scale}\right)}{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}}\right) \]
            8. lift-*.f64N/A

              \[\leadsto \frac{1}{4} \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\frac{{\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{x-scale \cdot x-scale} - \frac{1}{2} \cdot \frac{\mathsf{fma}\left(-2, \frac{{\left(\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)\right)}^{2}}{x-scale \cdot x-scale}, 4 \cdot \frac{{\left(\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)\right)}^{2}}{x-scale \cdot x-scale}\right)}{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}}\right) \]
            9. lift-PI.f64N/A

              \[\leadsto \frac{1}{4} \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\frac{{\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{x-scale \cdot x-scale} - \frac{1}{2} \cdot \frac{\mathsf{fma}\left(-2, \frac{{\left(\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)\right)}^{2}}{x-scale \cdot x-scale}, 4 \cdot \frac{{\left(\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)\right)}^{2}}{x-scale \cdot x-scale}\right)}{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}}\right) \]
            10. lift-*.f64N/A

              \[\leadsto \frac{1}{4} \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\frac{{\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{x-scale \cdot x-scale} - \frac{1}{2} \cdot \frac{\mathsf{fma}\left(-2, \frac{{\left(\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)\right)}^{2}}{x-scale \cdot x-scale}, 4 \cdot \frac{{\left(\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)\right)}^{2}}{x-scale \cdot x-scale}\right)}{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}}\right) \]
            11. unpow-prod-downN/A

              \[\leadsto \frac{1}{4} \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\frac{{\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{x-scale \cdot x-scale} - \frac{1}{2} \cdot \frac{\mathsf{fma}\left(-2, \frac{{\left(\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)\right)}^{2}}{x-scale \cdot x-scale}, 4 \cdot \frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} \cdot {\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{x-scale \cdot x-scale}\right)}{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}}\right) \]
            12. lower-*.f64N/A

              \[\leadsto \frac{1}{4} \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\frac{{\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{x-scale \cdot x-scale} - \frac{1}{2} \cdot \frac{\mathsf{fma}\left(-2, \frac{{\left(\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)\right)}^{2}}{x-scale \cdot x-scale}, 4 \cdot \frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} \cdot {\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{x-scale \cdot x-scale}\right)}{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}}\right) \]
          3. Applied rewrites21.0%

            \[\leadsto 0.25 \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\frac{{\sin \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{x-scale \cdot x-scale} - 0.5 \cdot \frac{\mathsf{fma}\left(-2, \frac{{\left(\cos \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \sin \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)\right)}^{2}}{x-scale \cdot x-scale}, 4 \cdot \frac{{\cos \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} \cdot {\sin \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{x-scale \cdot x-scale}\right)}{{\cos \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}}\right) \]

          if 1.0800000000000001e-120 < y-scale

          1. Initial program 0.1%

            \[\frac{-\sqrt{\left(\left(2 \cdot \frac{4 \cdot \left(\left(b \cdot a\right) \cdot \left(b \cdot \left(-a\right)\right)\right)}{{\left(x-scale \cdot y-scale\right)}^{2}}\right) \cdot \left(\left(b \cdot a\right) \cdot \left(b \cdot \left(-a\right)\right)\right)\right) \cdot \left(\left(\frac{\frac{{\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}}{x-scale}}{x-scale} + \frac{\frac{{\left(a \cdot \cos \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} + {\left(b \cdot \sin \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale}}{y-scale}\right) - \sqrt{{\left(\frac{\frac{{\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}}{x-scale}}{x-scale} - \frac{\frac{{\left(a \cdot \cos \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} + {\left(b \cdot \sin \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale}}{y-scale}\right)}^{2} + {\left(\frac{\frac{\left(\left(2 \cdot \left({b}^{2} - {a}^{2}\right)\right) \cdot \sin \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \cos \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)}{x-scale}}{y-scale}\right)}^{2}}\right)}}{\frac{4 \cdot \left(\left(b \cdot a\right) \cdot \left(b \cdot \left(-a\right)\right)\right)}{{\left(x-scale \cdot y-scale\right)}^{2}}} \]
          2. Add Preprocessing
          3. Taylor expanded in angle around 0

            \[\leadsto \color{blue}{\frac{1}{4} \cdot \left(a \cdot \left(x-scale \cdot \left(\sqrt{2} \cdot \sqrt{8}\right)\right)\right)} \]
          4. Step-by-step derivation
            1. Applied rewrites25.5%

              \[\leadsto \color{blue}{0.25 \cdot \left(a \cdot \left(x-scale \cdot 4\right)\right)} \]
            2. Taylor expanded in a around 0

              \[\leadsto a \cdot \color{blue}{x-scale} \]
            3. Step-by-step derivation
              1. lower-*.f6425.5

                \[\leadsto a \cdot x-scale \]
            4. Applied rewrites25.5%

              \[\leadsto a \cdot \color{blue}{x-scale} \]
          5. Recombined 2 regimes into one program.
          6. Add Preprocessing

          Alternative 3: 36.2% accurate, 2.5× speedup?

          \[\begin{array}{l} a_m = \left|a\right| \\ x-scale_m = \left|x-scale\right| \\ y-scale_m = \left|y-scale\right| \\ \begin{array}{l} t_0 := 0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\\ t_1 := \cos t\_0\\ t_2 := \sin t\_0\\ t_3 := {\left(t\_1 \cdot t\_2\right)}^{2}\\ \mathbf{if}\;y-scale\_m \leq 1.08 \cdot 10^{-120}:\\ \;\;\;\;0.25 \cdot \left(\left(a\_m \cdot \left(x-scale\_m \cdot \left(y-scale\_m \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\frac{{t\_2}^{2} - 0.5 \cdot \frac{\mathsf{fma}\left(-2, t\_3, 4 \cdot t\_3\right)}{{t\_1}^{2}}}{x-scale\_m \cdot x-scale\_m}}\right)\\ \mathbf{else}:\\ \;\;\;\;a\_m \cdot x-scale\_m\\ \end{array} \end{array} \]
          a_m = (fabs.f64 a)
          x-scale_m = (fabs.f64 x-scale)
          y-scale_m = (fabs.f64 y-scale)
          (FPCore (a_m b angle x-scale_m y-scale_m)
           :precision binary64
           (let* ((t_0 (* 0.005555555555555556 (* angle (PI))))
                  (t_1 (cos t_0))
                  (t_2 (sin t_0))
                  (t_3 (pow (* t_1 t_2) 2.0)))
             (if (<= y-scale_m 1.08e-120)
               (*
                0.25
                (*
                 (* a_m (* x-scale_m (* y-scale_m (sqrt 8.0))))
                 (sqrt
                  (/
                   (- (pow t_2 2.0) (* 0.5 (/ (fma -2.0 t_3 (* 4.0 t_3)) (pow t_1 2.0))))
                   (* x-scale_m x-scale_m)))))
               (* a_m x-scale_m))))
          \begin{array}{l}
          a_m = \left|a\right|
          \\
          x-scale_m = \left|x-scale\right|
          \\
          y-scale_m = \left|y-scale\right|
          
          \\
          \begin{array}{l}
          t_0 := 0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\\
          t_1 := \cos t\_0\\
          t_2 := \sin t\_0\\
          t_3 := {\left(t\_1 \cdot t\_2\right)}^{2}\\
          \mathbf{if}\;y-scale\_m \leq 1.08 \cdot 10^{-120}:\\
          \;\;\;\;0.25 \cdot \left(\left(a\_m \cdot \left(x-scale\_m \cdot \left(y-scale\_m \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\frac{{t\_2}^{2} - 0.5 \cdot \frac{\mathsf{fma}\left(-2, t\_3, 4 \cdot t\_3\right)}{{t\_1}^{2}}}{x-scale\_m \cdot x-scale\_m}}\right)\\
          
          \mathbf{else}:\\
          \;\;\;\;a\_m \cdot x-scale\_m\\
          
          
          \end{array}
          \end{array}
          
          Derivation
          1. Split input into 2 regimes
          2. if y-scale < 1.0800000000000001e-120

            1. Initial program 0.0%

              \[\frac{-\sqrt{\left(\left(2 \cdot \frac{4 \cdot \left(\left(b \cdot a\right) \cdot \left(b \cdot \left(-a\right)\right)\right)}{{\left(x-scale \cdot y-scale\right)}^{2}}\right) \cdot \left(\left(b \cdot a\right) \cdot \left(b \cdot \left(-a\right)\right)\right)\right) \cdot \left(\left(\frac{\frac{{\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}}{x-scale}}{x-scale} + \frac{\frac{{\left(a \cdot \cos \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} + {\left(b \cdot \sin \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale}}{y-scale}\right) - \sqrt{{\left(\frac{\frac{{\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}}{x-scale}}{x-scale} - \frac{\frac{{\left(a \cdot \cos \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} + {\left(b \cdot \sin \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale}}{y-scale}\right)}^{2} + {\left(\frac{\frac{\left(\left(2 \cdot \left({b}^{2} - {a}^{2}\right)\right) \cdot \sin \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \cos \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)}{x-scale}}{y-scale}\right)}^{2}}\right)}}{\frac{4 \cdot \left(\left(b \cdot a\right) \cdot \left(b \cdot \left(-a\right)\right)\right)}{{\left(x-scale \cdot y-scale\right)}^{2}}} \]
            2. Add Preprocessing
            3. Taylor expanded in a around inf

              \[\leadsto \color{blue}{\frac{1}{4} \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\left(\frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{y-scale}^{2}} + \frac{{\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{x-scale}^{2}}\right) - \sqrt{4 \cdot \frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} \cdot {\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{x-scale}^{2} \cdot {y-scale}^{2}} + {\left(\frac{{\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{x-scale}^{2}} - \frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{y-scale}^{2}}\right)}^{2}}}\right)} \]
            4. Applied rewrites2.1%

              \[\leadsto \color{blue}{0.25 \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\left(\frac{{\cos \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale \cdot y-scale} + \frac{{\sin \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{x-scale \cdot x-scale}\right) - \sqrt{\mathsf{fma}\left(4, \frac{{\left(\cos \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \sin \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)\right)}^{2}}{{\left(x-scale \cdot y-scale\right)}^{2}}, {\left(\frac{{\sin \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{x-scale \cdot x-scale} - \frac{{\cos \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale \cdot y-scale}\right)}^{2}\right)}}\right)} \]
            5. Taylor expanded in y-scale around 0

              \[\leadsto \frac{1}{4} \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\frac{{\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{x-scale}^{2}} - \frac{1}{2} \cdot \frac{-2 \cdot \frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} \cdot {\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{x-scale}^{2}} + 4 \cdot \frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} \cdot {\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{x-scale}^{2}}}{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}}\right) \]
            6. Step-by-step derivation
              1. Applied rewrites20.5%

                \[\leadsto 0.25 \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\frac{{\sin \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{x-scale \cdot x-scale} - 0.5 \cdot \frac{\mathsf{fma}\left(-2, \frac{{\left(\cos \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \sin \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)\right)}^{2}}{x-scale \cdot x-scale}, 4 \cdot \frac{{\left(\cos \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \sin \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)\right)}^{2}}{x-scale \cdot x-scale}\right)}{{\cos \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}}\right) \]
              2. Taylor expanded in x-scale around 0

                \[\leadsto \frac{1}{4} \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\frac{{\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} - \frac{1}{2} \cdot \frac{-2 \cdot \left({\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} \cdot {\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}\right) + 4 \cdot \left({\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} \cdot {\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}\right)}{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}}{{x-scale}^{2}}}\right) \]
              3. Step-by-step derivation
                1. Applied rewrites21.7%

                  \[\leadsto 0.25 \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\frac{{\sin \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} - 0.5 \cdot \frac{\mathsf{fma}\left(-2, {\left(\cos \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \sin \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)\right)}^{2}, 4 \cdot {\left(\cos \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \sin \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)\right)}^{2}\right)}{{\cos \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}}{x-scale \cdot x-scale}}\right) \]

                if 1.0800000000000001e-120 < y-scale

                1. Initial program 0.1%

                  \[\frac{-\sqrt{\left(\left(2 \cdot \frac{4 \cdot \left(\left(b \cdot a\right) \cdot \left(b \cdot \left(-a\right)\right)\right)}{{\left(x-scale \cdot y-scale\right)}^{2}}\right) \cdot \left(\left(b \cdot a\right) \cdot \left(b \cdot \left(-a\right)\right)\right)\right) \cdot \left(\left(\frac{\frac{{\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}}{x-scale}}{x-scale} + \frac{\frac{{\left(a \cdot \cos \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} + {\left(b \cdot \sin \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale}}{y-scale}\right) - \sqrt{{\left(\frac{\frac{{\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}}{x-scale}}{x-scale} - \frac{\frac{{\left(a \cdot \cos \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} + {\left(b \cdot \sin \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale}}{y-scale}\right)}^{2} + {\left(\frac{\frac{\left(\left(2 \cdot \left({b}^{2} - {a}^{2}\right)\right) \cdot \sin \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \cos \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)}{x-scale}}{y-scale}\right)}^{2}}\right)}}{\frac{4 \cdot \left(\left(b \cdot a\right) \cdot \left(b \cdot \left(-a\right)\right)\right)}{{\left(x-scale \cdot y-scale\right)}^{2}}} \]
                2. Add Preprocessing
                3. Taylor expanded in angle around 0

                  \[\leadsto \color{blue}{\frac{1}{4} \cdot \left(a \cdot \left(x-scale \cdot \left(\sqrt{2} \cdot \sqrt{8}\right)\right)\right)} \]
                4. Step-by-step derivation
                  1. Applied rewrites25.5%

                    \[\leadsto \color{blue}{0.25 \cdot \left(a \cdot \left(x-scale \cdot 4\right)\right)} \]
                  2. Taylor expanded in a around 0

                    \[\leadsto a \cdot \color{blue}{x-scale} \]
                  3. Step-by-step derivation
                    1. lower-*.f6425.5

                      \[\leadsto a \cdot x-scale \]
                  4. Applied rewrites25.5%

                    \[\leadsto a \cdot \color{blue}{x-scale} \]
                5. Recombined 2 regimes into one program.
                6. Add Preprocessing

                Alternative 4: 37.5% accurate, 2.9× speedup?

                \[\begin{array}{l} a_m = \left|a\right| \\ x-scale_m = \left|x-scale\right| \\ y-scale_m = \left|y-scale\right| \\ \begin{array}{l} t_0 := 0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\\ t_1 := a\_m \cdot \left(x-scale\_m \cdot \left(y-scale\_m \cdot \sqrt{8}\right)\right)\\ t_2 := \sin t\_0\\ t_3 := \cos t\_0\\ t_4 := \frac{{t\_3}^{2}}{y-scale\_m \cdot y-scale\_m}\\ t_5 := \frac{{\left(t\_3 \cdot t\_2\right)}^{2}}{x-scale\_m \cdot x-scale\_m}\\ \mathbf{if}\;y-scale\_m \leq 1.8 \cdot 10^{-151}:\\ \;\;\;\;0.25 \cdot \left(t\_1 \cdot \sqrt{\frac{{t\_2}^{2}}{x-scale\_m \cdot x-scale\_m} - 0.5 \cdot \frac{\mathsf{fma}\left(-2, t\_5, 4 \cdot t\_5\right)}{1}}\right)\\ \mathbf{elif}\;y-scale\_m \leq 1.55 \cdot 10^{-73}:\\ \;\;\;\;0.25 \cdot \left(t\_1 \cdot \sqrt{t\_4 - t\_4}\right)\\ \mathbf{else}:\\ \;\;\;\;a\_m \cdot x-scale\_m\\ \end{array} \end{array} \]
                a_m = (fabs.f64 a)
                x-scale_m = (fabs.f64 x-scale)
                y-scale_m = (fabs.f64 y-scale)
                (FPCore (a_m b angle x-scale_m y-scale_m)
                 :precision binary64
                 (let* ((t_0 (* 0.005555555555555556 (* angle (PI))))
                        (t_1 (* a_m (* x-scale_m (* y-scale_m (sqrt 8.0)))))
                        (t_2 (sin t_0))
                        (t_3 (cos t_0))
                        (t_4 (/ (pow t_3 2.0) (* y-scale_m y-scale_m)))
                        (t_5 (/ (pow (* t_3 t_2) 2.0) (* x-scale_m x-scale_m))))
                   (if (<= y-scale_m 1.8e-151)
                     (*
                      0.25
                      (*
                       t_1
                       (sqrt
                        (-
                         (/ (pow t_2 2.0) (* x-scale_m x-scale_m))
                         (* 0.5 (/ (fma -2.0 t_5 (* 4.0 t_5)) 1.0))))))
                     (if (<= y-scale_m 1.55e-73)
                       (* 0.25 (* t_1 (sqrt (- t_4 t_4))))
                       (* a_m x-scale_m)))))
                \begin{array}{l}
                a_m = \left|a\right|
                \\
                x-scale_m = \left|x-scale\right|
                \\
                y-scale_m = \left|y-scale\right|
                
                \\
                \begin{array}{l}
                t_0 := 0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\\
                t_1 := a\_m \cdot \left(x-scale\_m \cdot \left(y-scale\_m \cdot \sqrt{8}\right)\right)\\
                t_2 := \sin t\_0\\
                t_3 := \cos t\_0\\
                t_4 := \frac{{t\_3}^{2}}{y-scale\_m \cdot y-scale\_m}\\
                t_5 := \frac{{\left(t\_3 \cdot t\_2\right)}^{2}}{x-scale\_m \cdot x-scale\_m}\\
                \mathbf{if}\;y-scale\_m \leq 1.8 \cdot 10^{-151}:\\
                \;\;\;\;0.25 \cdot \left(t\_1 \cdot \sqrt{\frac{{t\_2}^{2}}{x-scale\_m \cdot x-scale\_m} - 0.5 \cdot \frac{\mathsf{fma}\left(-2, t\_5, 4 \cdot t\_5\right)}{1}}\right)\\
                
                \mathbf{elif}\;y-scale\_m \leq 1.55 \cdot 10^{-73}:\\
                \;\;\;\;0.25 \cdot \left(t\_1 \cdot \sqrt{t\_4 - t\_4}\right)\\
                
                \mathbf{else}:\\
                \;\;\;\;a\_m \cdot x-scale\_m\\
                
                
                \end{array}
                \end{array}
                
                Derivation
                1. Split input into 3 regimes
                2. if y-scale < 1.80000000000000016e-151

                  1. Initial program 0.0%

                    \[\frac{-\sqrt{\left(\left(2 \cdot \frac{4 \cdot \left(\left(b \cdot a\right) \cdot \left(b \cdot \left(-a\right)\right)\right)}{{\left(x-scale \cdot y-scale\right)}^{2}}\right) \cdot \left(\left(b \cdot a\right) \cdot \left(b \cdot \left(-a\right)\right)\right)\right) \cdot \left(\left(\frac{\frac{{\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}}{x-scale}}{x-scale} + \frac{\frac{{\left(a \cdot \cos \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} + {\left(b \cdot \sin \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale}}{y-scale}\right) - \sqrt{{\left(\frac{\frac{{\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}}{x-scale}}{x-scale} - \frac{\frac{{\left(a \cdot \cos \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} + {\left(b \cdot \sin \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale}}{y-scale}\right)}^{2} + {\left(\frac{\frac{\left(\left(2 \cdot \left({b}^{2} - {a}^{2}\right)\right) \cdot \sin \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \cos \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)}{x-scale}}{y-scale}\right)}^{2}}\right)}}{\frac{4 \cdot \left(\left(b \cdot a\right) \cdot \left(b \cdot \left(-a\right)\right)\right)}{{\left(x-scale \cdot y-scale\right)}^{2}}} \]
                  2. Add Preprocessing
                  3. Taylor expanded in a around inf

                    \[\leadsto \color{blue}{\frac{1}{4} \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\left(\frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{y-scale}^{2}} + \frac{{\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{x-scale}^{2}}\right) - \sqrt{4 \cdot \frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} \cdot {\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{x-scale}^{2} \cdot {y-scale}^{2}} + {\left(\frac{{\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{x-scale}^{2}} - \frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{y-scale}^{2}}\right)}^{2}}}\right)} \]
                  4. Applied rewrites2.1%

                    \[\leadsto \color{blue}{0.25 \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\left(\frac{{\cos \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale \cdot y-scale} + \frac{{\sin \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{x-scale \cdot x-scale}\right) - \sqrt{\mathsf{fma}\left(4, \frac{{\left(\cos \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \sin \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)\right)}^{2}}{{\left(x-scale \cdot y-scale\right)}^{2}}, {\left(\frac{{\sin \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{x-scale \cdot x-scale} - \frac{{\cos \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale \cdot y-scale}\right)}^{2}\right)}}\right)} \]
                  5. Taylor expanded in y-scale around 0

                    \[\leadsto \frac{1}{4} \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\frac{{\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{x-scale}^{2}} - \frac{1}{2} \cdot \frac{-2 \cdot \frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} \cdot {\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{x-scale}^{2}} + 4 \cdot \frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} \cdot {\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{x-scale}^{2}}}{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}}\right) \]
                  6. Step-by-step derivation
                    1. Applied rewrites21.1%

                      \[\leadsto 0.25 \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\frac{{\sin \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{x-scale \cdot x-scale} - 0.5 \cdot \frac{\mathsf{fma}\left(-2, \frac{{\left(\cos \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \sin \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)\right)}^{2}}{x-scale \cdot x-scale}, 4 \cdot \frac{{\left(\cos \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \sin \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)\right)}^{2}}{x-scale \cdot x-scale}\right)}{{\cos \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}}\right) \]
                    2. Taylor expanded in angle around 0

                      \[\leadsto \frac{1}{4} \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\frac{{\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{x-scale \cdot x-scale} - \frac{1}{2} \cdot \frac{\mathsf{fma}\left(-2, \frac{{\left(\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)\right)}^{2}}{x-scale \cdot x-scale}, 4 \cdot \frac{{\left(\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)\right)}^{2}}{x-scale \cdot x-scale}\right)}{1}}\right) \]
                    3. Step-by-step derivation
                      1. Applied rewrites20.7%

                        \[\leadsto 0.25 \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\frac{{\sin \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{x-scale \cdot x-scale} - 0.5 \cdot \frac{\mathsf{fma}\left(-2, \frac{{\left(\cos \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \sin \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)\right)}^{2}}{x-scale \cdot x-scale}, 4 \cdot \frac{{\left(\cos \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \sin \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)\right)}^{2}}{x-scale \cdot x-scale}\right)}{1}}\right) \]

                      if 1.80000000000000016e-151 < y-scale < 1.54999999999999985e-73

                      1. Initial program 0.0%

                        \[\frac{-\sqrt{\left(\left(2 \cdot \frac{4 \cdot \left(\left(b \cdot a\right) \cdot \left(b \cdot \left(-a\right)\right)\right)}{{\left(x-scale \cdot y-scale\right)}^{2}}\right) \cdot \left(\left(b \cdot a\right) \cdot \left(b \cdot \left(-a\right)\right)\right)\right) \cdot \left(\left(\frac{\frac{{\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}}{x-scale}}{x-scale} + \frac{\frac{{\left(a \cdot \cos \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} + {\left(b \cdot \sin \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale}}{y-scale}\right) - \sqrt{{\left(\frac{\frac{{\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}}{x-scale}}{x-scale} - \frac{\frac{{\left(a \cdot \cos \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} + {\left(b \cdot \sin \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale}}{y-scale}\right)}^{2} + {\left(\frac{\frac{\left(\left(2 \cdot \left({b}^{2} - {a}^{2}\right)\right) \cdot \sin \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \cos \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)}{x-scale}}{y-scale}\right)}^{2}}\right)}}{\frac{4 \cdot \left(\left(b \cdot a\right) \cdot \left(b \cdot \left(-a\right)\right)\right)}{{\left(x-scale \cdot y-scale\right)}^{2}}} \]
                      2. Add Preprocessing
                      3. Taylor expanded in a around inf

                        \[\leadsto \color{blue}{\frac{1}{4} \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\left(\frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{y-scale}^{2}} + \frac{{\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{x-scale}^{2}}\right) - \sqrt{4 \cdot \frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} \cdot {\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{x-scale}^{2} \cdot {y-scale}^{2}} + {\left(\frac{{\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{x-scale}^{2}} - \frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{y-scale}^{2}}\right)}^{2}}}\right)} \]
                      4. Applied rewrites0.0%

                        \[\leadsto \color{blue}{0.25 \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\left(\frac{{\cos \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale \cdot y-scale} + \frac{{\sin \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{x-scale \cdot x-scale}\right) - \sqrt{\mathsf{fma}\left(4, \frac{{\left(\cos \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \sin \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)\right)}^{2}}{{\left(x-scale \cdot y-scale\right)}^{2}}, {\left(\frac{{\sin \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{x-scale \cdot x-scale} - \frac{{\cos \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale \cdot y-scale}\right)}^{2}\right)}}\right)} \]
                      5. Taylor expanded in x-scale around inf

                        \[\leadsto \frac{1}{4} \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\left(\frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale \cdot y-scale} + \frac{{\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{x-scale \cdot x-scale}\right) - \frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{y-scale}^{2}}}\right) \]
                      6. Step-by-step derivation
                        1. lift-*.f64N/A

                          \[\leadsto \frac{1}{4} \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\left(\frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale \cdot y-scale} + \frac{{\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{x-scale \cdot x-scale}\right) - \frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{y-scale}^{2}}}\right) \]
                        2. lift-PI.f64N/A

                          \[\leadsto \frac{1}{4} \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\left(\frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale \cdot y-scale} + \frac{{\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{x-scale \cdot x-scale}\right) - \frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{y-scale}^{2}}}\right) \]
                        3. lift-*.f64N/A

                          \[\leadsto \frac{1}{4} \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\left(\frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale \cdot y-scale} + \frac{{\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{x-scale \cdot x-scale}\right) - \frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{y-scale}^{2}}}\right) \]
                        4. lift-cos.f64N/A

                          \[\leadsto \frac{1}{4} \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\left(\frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale \cdot y-scale} + \frac{{\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{x-scale \cdot x-scale}\right) - \frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{y-scale}^{2}}}\right) \]
                        5. lift-pow.f64N/A

                          \[\leadsto \frac{1}{4} \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\left(\frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale \cdot y-scale} + \frac{{\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{x-scale \cdot x-scale}\right) - \frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{y-scale}^{2}}}\right) \]
                        6. pow2N/A

                          \[\leadsto \frac{1}{4} \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\left(\frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale \cdot y-scale} + \frac{{\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{x-scale \cdot x-scale}\right) - \frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale \cdot y-scale}}\right) \]
                        7. lift-/.f64N/A

                          \[\leadsto \frac{1}{4} \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\left(\frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale \cdot y-scale} + \frac{{\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{x-scale \cdot x-scale}\right) - \frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale \cdot y-scale}}\right) \]
                        8. lift-*.f6418.8

                          \[\leadsto 0.25 \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\left(\frac{{\cos \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale \cdot y-scale} + \frac{{\sin \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{x-scale \cdot x-scale}\right) - \frac{{\cos \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale \cdot y-scale}}\right) \]
                      7. Applied rewrites18.8%

                        \[\leadsto 0.25 \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\left(\frac{{\cos \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale \cdot y-scale} + \frac{{\sin \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{x-scale \cdot x-scale}\right) - \frac{{\cos \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale \cdot y-scale}}\right) \]
                      8. Taylor expanded in x-scale around inf

                        \[\leadsto \frac{1}{4} \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{y-scale}^{2}} - \frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale \cdot y-scale}}\right) \]
                      9. Step-by-step derivation
                        1. lift-*.f64N/A

                          \[\leadsto \frac{1}{4} \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{y-scale}^{2}} - \frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale \cdot y-scale}}\right) \]
                        2. lift-PI.f64N/A

                          \[\leadsto \frac{1}{4} \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{y-scale}^{2}} - \frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale \cdot y-scale}}\right) \]
                        3. lift-*.f64N/A

                          \[\leadsto \frac{1}{4} \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{y-scale}^{2}} - \frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale \cdot y-scale}}\right) \]
                        4. lift-cos.f64N/A

                          \[\leadsto \frac{1}{4} \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{y-scale}^{2}} - \frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale \cdot y-scale}}\right) \]
                        5. lift-pow.f64N/A

                          \[\leadsto \frac{1}{4} \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{y-scale}^{2}} - \frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale \cdot y-scale}}\right) \]
                        6. pow2N/A

                          \[\leadsto \frac{1}{4} \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale \cdot y-scale} - \frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale \cdot y-scale}}\right) \]
                        7. lift-/.f64N/A

                          \[\leadsto \frac{1}{4} \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale \cdot y-scale} - \frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale \cdot y-scale}}\right) \]
                        8. lift-*.f6444.8

                          \[\leadsto 0.25 \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\frac{{\cos \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale \cdot y-scale} - \frac{{\cos \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale \cdot y-scale}}\right) \]
                      10. Applied rewrites44.8%

                        \[\leadsto 0.25 \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\frac{{\cos \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale \cdot y-scale} - \frac{{\cos \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale \cdot y-scale}}\right) \]

                      if 1.54999999999999985e-73 < y-scale

                      1. Initial program 0.1%

                        \[\frac{-\sqrt{\left(\left(2 \cdot \frac{4 \cdot \left(\left(b \cdot a\right) \cdot \left(b \cdot \left(-a\right)\right)\right)}{{\left(x-scale \cdot y-scale\right)}^{2}}\right) \cdot \left(\left(b \cdot a\right) \cdot \left(b \cdot \left(-a\right)\right)\right)\right) \cdot \left(\left(\frac{\frac{{\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}}{x-scale}}{x-scale} + \frac{\frac{{\left(a \cdot \cos \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} + {\left(b \cdot \sin \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale}}{y-scale}\right) - \sqrt{{\left(\frac{\frac{{\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}}{x-scale}}{x-scale} - \frac{\frac{{\left(a \cdot \cos \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} + {\left(b \cdot \sin \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale}}{y-scale}\right)}^{2} + {\left(\frac{\frac{\left(\left(2 \cdot \left({b}^{2} - {a}^{2}\right)\right) \cdot \sin \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \cos \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)}{x-scale}}{y-scale}\right)}^{2}}\right)}}{\frac{4 \cdot \left(\left(b \cdot a\right) \cdot \left(b \cdot \left(-a\right)\right)\right)}{{\left(x-scale \cdot y-scale\right)}^{2}}} \]
                      2. Add Preprocessing
                      3. Taylor expanded in angle around 0

                        \[\leadsto \color{blue}{\frac{1}{4} \cdot \left(a \cdot \left(x-scale \cdot \left(\sqrt{2} \cdot \sqrt{8}\right)\right)\right)} \]
                      4. Step-by-step derivation
                        1. Applied rewrites23.9%

                          \[\leadsto \color{blue}{0.25 \cdot \left(a \cdot \left(x-scale \cdot 4\right)\right)} \]
                        2. Taylor expanded in a around 0

                          \[\leadsto a \cdot \color{blue}{x-scale} \]
                        3. Step-by-step derivation
                          1. lower-*.f6423.9

                            \[\leadsto a \cdot x-scale \]
                        4. Applied rewrites23.9%

                          \[\leadsto a \cdot \color{blue}{x-scale} \]
                      5. Recombined 3 regimes into one program.
                      6. Add Preprocessing

                      Alternative 5: 34.5% accurate, 24.9× speedup?

                      \[\begin{array}{l} a_m = \left|a\right| \\ x-scale_m = \left|x-scale\right| \\ y-scale_m = \left|y-scale\right| \\ \begin{array}{l} t_0 := \mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\\ \mathbf{if}\;y-scale\_m \leq 6.8 \cdot 10^{-121}:\\ \;\;\;\;0.25 \cdot \left(\left(a\_m \cdot \left(x-scale\_m \cdot \left(y-scale\_m \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\left(angle \cdot angle\right) \cdot \frac{3.08641975308642 \cdot 10^{-5} \cdot t\_0 - 0.5 \cdot \mathsf{fma}\left(-6.17283950617284 \cdot 10^{-5}, t\_0, 0.0001234567901234568 \cdot t\_0\right)}{x-scale\_m \cdot x-scale\_m}}\right)\\ \mathbf{else}:\\ \;\;\;\;a\_m \cdot x-scale\_m\\ \end{array} \end{array} \]
                      a_m = (fabs.f64 a)
                      x-scale_m = (fabs.f64 x-scale)
                      y-scale_m = (fabs.f64 y-scale)
                      (FPCore (a_m b angle x-scale_m y-scale_m)
                       :precision binary64
                       (let* ((t_0 (* (PI) (PI))))
                         (if (<= y-scale_m 6.8e-121)
                           (*
                            0.25
                            (*
                             (* a_m (* x-scale_m (* y-scale_m (sqrt 8.0))))
                             (sqrt
                              (*
                               (* angle angle)
                               (/
                                (-
                                 (* 3.08641975308642e-5 t_0)
                                 (*
                                  0.5
                                  (fma -6.17283950617284e-5 t_0 (* 0.0001234567901234568 t_0))))
                                (* x-scale_m x-scale_m))))))
                           (* a_m x-scale_m))))
                      \begin{array}{l}
                      a_m = \left|a\right|
                      \\
                      x-scale_m = \left|x-scale\right|
                      \\
                      y-scale_m = \left|y-scale\right|
                      
                      \\
                      \begin{array}{l}
                      t_0 := \mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\\
                      \mathbf{if}\;y-scale\_m \leq 6.8 \cdot 10^{-121}:\\
                      \;\;\;\;0.25 \cdot \left(\left(a\_m \cdot \left(x-scale\_m \cdot \left(y-scale\_m \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\left(angle \cdot angle\right) \cdot \frac{3.08641975308642 \cdot 10^{-5} \cdot t\_0 - 0.5 \cdot \mathsf{fma}\left(-6.17283950617284 \cdot 10^{-5}, t\_0, 0.0001234567901234568 \cdot t\_0\right)}{x-scale\_m \cdot x-scale\_m}}\right)\\
                      
                      \mathbf{else}:\\
                      \;\;\;\;a\_m \cdot x-scale\_m\\
                      
                      
                      \end{array}
                      \end{array}
                      
                      Derivation
                      1. Split input into 2 regimes
                      2. if y-scale < 6.80000000000000003e-121

                        1. Initial program 0.0%

                          \[\frac{-\sqrt{\left(\left(2 \cdot \frac{4 \cdot \left(\left(b \cdot a\right) \cdot \left(b \cdot \left(-a\right)\right)\right)}{{\left(x-scale \cdot y-scale\right)}^{2}}\right) \cdot \left(\left(b \cdot a\right) \cdot \left(b \cdot \left(-a\right)\right)\right)\right) \cdot \left(\left(\frac{\frac{{\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}}{x-scale}}{x-scale} + \frac{\frac{{\left(a \cdot \cos \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} + {\left(b \cdot \sin \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale}}{y-scale}\right) - \sqrt{{\left(\frac{\frac{{\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}}{x-scale}}{x-scale} - \frac{\frac{{\left(a \cdot \cos \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} + {\left(b \cdot \sin \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale}}{y-scale}\right)}^{2} + {\left(\frac{\frac{\left(\left(2 \cdot \left({b}^{2} - {a}^{2}\right)\right) \cdot \sin \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \cos \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)}{x-scale}}{y-scale}\right)}^{2}}\right)}}{\frac{4 \cdot \left(\left(b \cdot a\right) \cdot \left(b \cdot \left(-a\right)\right)\right)}{{\left(x-scale \cdot y-scale\right)}^{2}}} \]
                        2. Add Preprocessing
                        3. Taylor expanded in a around inf

                          \[\leadsto \color{blue}{\frac{1}{4} \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\left(\frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{y-scale}^{2}} + \frac{{\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{x-scale}^{2}}\right) - \sqrt{4 \cdot \frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} \cdot {\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{x-scale}^{2} \cdot {y-scale}^{2}} + {\left(\frac{{\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{x-scale}^{2}} - \frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{y-scale}^{2}}\right)}^{2}}}\right)} \]
                        4. Applied rewrites2.1%

                          \[\leadsto \color{blue}{0.25 \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\left(\frac{{\cos \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale \cdot y-scale} + \frac{{\sin \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{x-scale \cdot x-scale}\right) - \sqrt{\mathsf{fma}\left(4, \frac{{\left(\cos \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \sin \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)\right)}^{2}}{{\left(x-scale \cdot y-scale\right)}^{2}}, {\left(\frac{{\sin \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{x-scale \cdot x-scale} - \frac{{\cos \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale \cdot y-scale}\right)}^{2}\right)}}\right)} \]
                        5. Taylor expanded in y-scale around 0

                          \[\leadsto \frac{1}{4} \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\frac{{\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{x-scale}^{2}} - \frac{1}{2} \cdot \frac{-2 \cdot \frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} \cdot {\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{x-scale}^{2}} + 4 \cdot \frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} \cdot {\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{x-scale}^{2}}}{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}}\right) \]
                        6. Step-by-step derivation
                          1. Applied rewrites20.5%

                            \[\leadsto 0.25 \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\frac{{\sin \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{x-scale \cdot x-scale} - 0.5 \cdot \frac{\mathsf{fma}\left(-2, \frac{{\left(\cos \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \sin \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)\right)}^{2}}{x-scale \cdot x-scale}, 4 \cdot \frac{{\left(\cos \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \sin \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)\right)}^{2}}{x-scale \cdot x-scale}\right)}{{\cos \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}}\right) \]
                          2. Taylor expanded in angle around 0

                            \[\leadsto \frac{1}{4} \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{{angle}^{2} \cdot \left(\frac{1}{32400} \cdot \frac{{\mathsf{PI}\left(\right)}^{2}}{{x-scale}^{2}} - \frac{1}{2} \cdot \left(\frac{-1}{16200} \cdot \frac{{\mathsf{PI}\left(\right)}^{2}}{{x-scale}^{2}} + \frac{1}{8100} \cdot \frac{{\mathsf{PI}\left(\right)}^{2}}{{x-scale}^{2}}\right)\right)}\right) \]
                          3. Step-by-step derivation
                            1. lower-*.f64N/A

                              \[\leadsto \frac{1}{4} \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{{angle}^{2} \cdot \left(\frac{1}{32400} \cdot \frac{{\mathsf{PI}\left(\right)}^{2}}{{x-scale}^{2}} - \frac{1}{2} \cdot \left(\frac{-1}{16200} \cdot \frac{{\mathsf{PI}\left(\right)}^{2}}{{x-scale}^{2}} + \frac{1}{8100} \cdot \frac{{\mathsf{PI}\left(\right)}^{2}}{{x-scale}^{2}}\right)\right)}\right) \]
                            2. unpow2N/A

                              \[\leadsto \frac{1}{4} \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\left(angle \cdot angle\right) \cdot \left(\frac{1}{32400} \cdot \frac{{\mathsf{PI}\left(\right)}^{2}}{{x-scale}^{2}} - \frac{1}{2} \cdot \left(\frac{-1}{16200} \cdot \frac{{\mathsf{PI}\left(\right)}^{2}}{{x-scale}^{2}} + \frac{1}{8100} \cdot \frac{{\mathsf{PI}\left(\right)}^{2}}{{x-scale}^{2}}\right)\right)}\right) \]
                            3. lower-*.f64N/A

                              \[\leadsto \frac{1}{4} \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\left(angle \cdot angle\right) \cdot \left(\frac{1}{32400} \cdot \frac{{\mathsf{PI}\left(\right)}^{2}}{{x-scale}^{2}} - \frac{1}{2} \cdot \left(\frac{-1}{16200} \cdot \frac{{\mathsf{PI}\left(\right)}^{2}}{{x-scale}^{2}} + \frac{1}{8100} \cdot \frac{{\mathsf{PI}\left(\right)}^{2}}{{x-scale}^{2}}\right)\right)}\right) \]
                            4. lower--.f64N/A

                              \[\leadsto \frac{1}{4} \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\left(angle \cdot angle\right) \cdot \left(\frac{1}{32400} \cdot \frac{{\mathsf{PI}\left(\right)}^{2}}{{x-scale}^{2}} - \frac{1}{2} \cdot \left(\frac{-1}{16200} \cdot \frac{{\mathsf{PI}\left(\right)}^{2}}{{x-scale}^{2}} + \frac{1}{8100} \cdot \frac{{\mathsf{PI}\left(\right)}^{2}}{{x-scale}^{2}}\right)\right)}\right) \]
                          4. Applied rewrites18.6%

                            \[\leadsto 0.25 \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\left(angle \cdot angle\right) \cdot \left(3.08641975308642 \cdot 10^{-5} \cdot \frac{\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)}{x-scale \cdot x-scale} - 0.5 \cdot \mathsf{fma}\left(-6.17283950617284 \cdot 10^{-5}, \frac{\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)}{x-scale \cdot x-scale}, 0.0001234567901234568 \cdot \frac{\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)}{x-scale \cdot x-scale}\right)\right)}\right) \]
                          5. Taylor expanded in x-scale around 0

                            \[\leadsto \frac{1}{4} \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\left(angle \cdot angle\right) \cdot \frac{\frac{1}{32400} \cdot {\mathsf{PI}\left(\right)}^{2} - \frac{1}{2} \cdot \left(\frac{-1}{16200} \cdot {\mathsf{PI}\left(\right)}^{2} + \frac{1}{8100} \cdot {\mathsf{PI}\left(\right)}^{2}\right)}{{x-scale}^{2}}}\right) \]
                          6. Step-by-step derivation
                            1. lower-/.f64N/A

                              \[\leadsto \frac{1}{4} \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\left(angle \cdot angle\right) \cdot \frac{\frac{1}{32400} \cdot {\mathsf{PI}\left(\right)}^{2} - \frac{1}{2} \cdot \left(\frac{-1}{16200} \cdot {\mathsf{PI}\left(\right)}^{2} + \frac{1}{8100} \cdot {\mathsf{PI}\left(\right)}^{2}\right)}{{x-scale}^{2}}}\right) \]
                          7. Applied rewrites19.3%

                            \[\leadsto 0.25 \cdot \left(\left(a \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\left(angle \cdot angle\right) \cdot \frac{3.08641975308642 \cdot 10^{-5} \cdot \left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right) - 0.5 \cdot \mathsf{fma}\left(-6.17283950617284 \cdot 10^{-5}, \mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right), 0.0001234567901234568 \cdot \left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right)\right)}{x-scale \cdot x-scale}}\right) \]

                          if 6.80000000000000003e-121 < y-scale

                          1. Initial program 0.1%

                            \[\frac{-\sqrt{\left(\left(2 \cdot \frac{4 \cdot \left(\left(b \cdot a\right) \cdot \left(b \cdot \left(-a\right)\right)\right)}{{\left(x-scale \cdot y-scale\right)}^{2}}\right) \cdot \left(\left(b \cdot a\right) \cdot \left(b \cdot \left(-a\right)\right)\right)\right) \cdot \left(\left(\frac{\frac{{\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}}{x-scale}}{x-scale} + \frac{\frac{{\left(a \cdot \cos \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} + {\left(b \cdot \sin \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale}}{y-scale}\right) - \sqrt{{\left(\frac{\frac{{\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}}{x-scale}}{x-scale} - \frac{\frac{{\left(a \cdot \cos \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} + {\left(b \cdot \sin \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale}}{y-scale}\right)}^{2} + {\left(\frac{\frac{\left(\left(2 \cdot \left({b}^{2} - {a}^{2}\right)\right) \cdot \sin \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \cos \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)}{x-scale}}{y-scale}\right)}^{2}}\right)}}{\frac{4 \cdot \left(\left(b \cdot a\right) \cdot \left(b \cdot \left(-a\right)\right)\right)}{{\left(x-scale \cdot y-scale\right)}^{2}}} \]
                          2. Add Preprocessing
                          3. Taylor expanded in angle around 0

                            \[\leadsto \color{blue}{\frac{1}{4} \cdot \left(a \cdot \left(x-scale \cdot \left(\sqrt{2} \cdot \sqrt{8}\right)\right)\right)} \]
                          4. Step-by-step derivation
                            1. Applied rewrites25.5%

                              \[\leadsto \color{blue}{0.25 \cdot \left(a \cdot \left(x-scale \cdot 4\right)\right)} \]
                            2. Taylor expanded in a around 0

                              \[\leadsto a \cdot \color{blue}{x-scale} \]
                            3. Step-by-step derivation
                              1. lower-*.f6425.5

                                \[\leadsto a \cdot x-scale \]
                            4. Applied rewrites25.5%

                              \[\leadsto a \cdot \color{blue}{x-scale} \]
                          5. Recombined 2 regimes into one program.
                          6. Add Preprocessing

                          Alternative 6: 34.9% accurate, 55.9× speedup?

                          \[\begin{array}{l} a_m = \left|a\right| \\ x-scale_m = \left|x-scale\right| \\ y-scale_m = \left|y-scale\right| \\ \begin{array}{l} \mathbf{if}\;y-scale\_m \leq 7.5 \cdot 10^{-28}:\\ \;\;\;\;0.25 \cdot \left(\left(x-scale\_m \cdot \sqrt{8}\right) \cdot \sqrt{2 \cdot \left(a\_m \cdot a\_m\right)}\right)\\ \mathbf{else}:\\ \;\;\;\;a\_m \cdot x-scale\_m\\ \end{array} \end{array} \]
                          a_m = (fabs.f64 a)
                          x-scale_m = (fabs.f64 x-scale)
                          y-scale_m = (fabs.f64 y-scale)
                          (FPCore (a_m b angle x-scale_m y-scale_m)
                           :precision binary64
                           (if (<= y-scale_m 7.5e-28)
                             (* 0.25 (* (* x-scale_m (sqrt 8.0)) (sqrt (* 2.0 (* a_m a_m)))))
                             (* a_m x-scale_m)))
                          a_m = fabs(a);
                          x-scale_m = fabs(x_45_scale);
                          y-scale_m = fabs(y_45_scale);
                          double code(double a_m, double b, double angle, double x_45_scale_m, double y_45_scale_m) {
                          	double tmp;
                          	if (y_45_scale_m <= 7.5e-28) {
                          		tmp = 0.25 * ((x_45_scale_m * sqrt(8.0)) * sqrt((2.0 * (a_m * a_m))));
                          	} else {
                          		tmp = a_m * x_45_scale_m;
                          	}
                          	return tmp;
                          }
                          
                          a_m =     private
                          x-scale_m =     private
                          y-scale_m =     private
                          module fmin_fmax_functions
                              implicit none
                              private
                              public fmax
                              public fmin
                          
                              interface fmax
                                  module procedure fmax88
                                  module procedure fmax44
                                  module procedure fmax84
                                  module procedure fmax48
                              end interface
                              interface fmin
                                  module procedure fmin88
                                  module procedure fmin44
                                  module procedure fmin84
                                  module procedure fmin48
                              end interface
                          contains
                              real(8) function fmax88(x, y) result (res)
                                  real(8), intent (in) :: x
                                  real(8), intent (in) :: y
                                  res = merge(y, merge(x, max(x, y), y /= y), x /= x)
                              end function
                              real(4) function fmax44(x, y) result (res)
                                  real(4), intent (in) :: x
                                  real(4), intent (in) :: y
                                  res = merge(y, merge(x, max(x, y), y /= y), x /= x)
                              end function
                              real(8) function fmax84(x, y) result(res)
                                  real(8), intent (in) :: x
                                  real(4), intent (in) :: y
                                  res = merge(dble(y), merge(x, max(x, dble(y)), y /= y), x /= x)
                              end function
                              real(8) function fmax48(x, y) result(res)
                                  real(4), intent (in) :: x
                                  real(8), intent (in) :: y
                                  res = merge(y, merge(dble(x), max(dble(x), y), y /= y), x /= x)
                              end function
                              real(8) function fmin88(x, y) result (res)
                                  real(8), intent (in) :: x
                                  real(8), intent (in) :: y
                                  res = merge(y, merge(x, min(x, y), y /= y), x /= x)
                              end function
                              real(4) function fmin44(x, y) result (res)
                                  real(4), intent (in) :: x
                                  real(4), intent (in) :: y
                                  res = merge(y, merge(x, min(x, y), y /= y), x /= x)
                              end function
                              real(8) function fmin84(x, y) result(res)
                                  real(8), intent (in) :: x
                                  real(4), intent (in) :: y
                                  res = merge(dble(y), merge(x, min(x, dble(y)), y /= y), x /= x)
                              end function
                              real(8) function fmin48(x, y) result(res)
                                  real(4), intent (in) :: x
                                  real(8), intent (in) :: y
                                  res = merge(y, merge(dble(x), min(dble(x), y), y /= y), x /= x)
                              end function
                          end module
                          
                          real(8) function code(a_m, b, angle, x_45scale_m, y_45scale_m)
                          use fmin_fmax_functions
                              real(8), intent (in) :: a_m
                              real(8), intent (in) :: b
                              real(8), intent (in) :: angle
                              real(8), intent (in) :: x_45scale_m
                              real(8), intent (in) :: y_45scale_m
                              real(8) :: tmp
                              if (y_45scale_m <= 7.5d-28) then
                                  tmp = 0.25d0 * ((x_45scale_m * sqrt(8.0d0)) * sqrt((2.0d0 * (a_m * a_m))))
                              else
                                  tmp = a_m * x_45scale_m
                              end if
                              code = tmp
                          end function
                          
                          a_m = Math.abs(a);
                          x-scale_m = Math.abs(x_45_scale);
                          y-scale_m = Math.abs(y_45_scale);
                          public static double code(double a_m, double b, double angle, double x_45_scale_m, double y_45_scale_m) {
                          	double tmp;
                          	if (y_45_scale_m <= 7.5e-28) {
                          		tmp = 0.25 * ((x_45_scale_m * Math.sqrt(8.0)) * Math.sqrt((2.0 * (a_m * a_m))));
                          	} else {
                          		tmp = a_m * x_45_scale_m;
                          	}
                          	return tmp;
                          }
                          
                          a_m = math.fabs(a)
                          x-scale_m = math.fabs(x_45_scale)
                          y-scale_m = math.fabs(y_45_scale)
                          def code(a_m, b, angle, x_45_scale_m, y_45_scale_m):
                          	tmp = 0
                          	if y_45_scale_m <= 7.5e-28:
                          		tmp = 0.25 * ((x_45_scale_m * math.sqrt(8.0)) * math.sqrt((2.0 * (a_m * a_m))))
                          	else:
                          		tmp = a_m * x_45_scale_m
                          	return tmp
                          
                          a_m = abs(a)
                          x-scale_m = abs(x_45_scale)
                          y-scale_m = abs(y_45_scale)
                          function code(a_m, b, angle, x_45_scale_m, y_45_scale_m)
                          	tmp = 0.0
                          	if (y_45_scale_m <= 7.5e-28)
                          		tmp = Float64(0.25 * Float64(Float64(x_45_scale_m * sqrt(8.0)) * sqrt(Float64(2.0 * Float64(a_m * a_m)))));
                          	else
                          		tmp = Float64(a_m * x_45_scale_m);
                          	end
                          	return tmp
                          end
                          
                          a_m = abs(a);
                          x-scale_m = abs(x_45_scale);
                          y-scale_m = abs(y_45_scale);
                          function tmp_2 = code(a_m, b, angle, x_45_scale_m, y_45_scale_m)
                          	tmp = 0.0;
                          	if (y_45_scale_m <= 7.5e-28)
                          		tmp = 0.25 * ((x_45_scale_m * sqrt(8.0)) * sqrt((2.0 * (a_m * a_m))));
                          	else
                          		tmp = a_m * x_45_scale_m;
                          	end
                          	tmp_2 = tmp;
                          end
                          
                          a_m = N[Abs[a], $MachinePrecision]
                          x-scale_m = N[Abs[x$45$scale], $MachinePrecision]
                          y-scale_m = N[Abs[y$45$scale], $MachinePrecision]
                          code[a$95$m_, b_, angle_, x$45$scale$95$m_, y$45$scale$95$m_] := If[LessEqual[y$45$scale$95$m, 7.5e-28], N[(0.25 * N[(N[(x$45$scale$95$m * N[Sqrt[8.0], $MachinePrecision]), $MachinePrecision] * N[Sqrt[N[(2.0 * N[(a$95$m * a$95$m), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(a$95$m * x$45$scale$95$m), $MachinePrecision]]
                          
                          \begin{array}{l}
                          a_m = \left|a\right|
                          \\
                          x-scale_m = \left|x-scale\right|
                          \\
                          y-scale_m = \left|y-scale\right|
                          
                          \\
                          \begin{array}{l}
                          \mathbf{if}\;y-scale\_m \leq 7.5 \cdot 10^{-28}:\\
                          \;\;\;\;0.25 \cdot \left(\left(x-scale\_m \cdot \sqrt{8}\right) \cdot \sqrt{2 \cdot \left(a\_m \cdot a\_m\right)}\right)\\
                          
                          \mathbf{else}:\\
                          \;\;\;\;a\_m \cdot x-scale\_m\\
                          
                          
                          \end{array}
                          \end{array}
                          
                          Derivation
                          1. Split input into 2 regimes
                          2. if y-scale < 7.5000000000000003e-28

                            1. Initial program 0.0%

                              \[\frac{-\sqrt{\left(\left(2 \cdot \frac{4 \cdot \left(\left(b \cdot a\right) \cdot \left(b \cdot \left(-a\right)\right)\right)}{{\left(x-scale \cdot y-scale\right)}^{2}}\right) \cdot \left(\left(b \cdot a\right) \cdot \left(b \cdot \left(-a\right)\right)\right)\right) \cdot \left(\left(\frac{\frac{{\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}}{x-scale}}{x-scale} + \frac{\frac{{\left(a \cdot \cos \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} + {\left(b \cdot \sin \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale}}{y-scale}\right) - \sqrt{{\left(\frac{\frac{{\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}}{x-scale}}{x-scale} - \frac{\frac{{\left(a \cdot \cos \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} + {\left(b \cdot \sin \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale}}{y-scale}\right)}^{2} + {\left(\frac{\frac{\left(\left(2 \cdot \left({b}^{2} - {a}^{2}\right)\right) \cdot \sin \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \cos \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)}{x-scale}}{y-scale}\right)}^{2}}\right)}}{\frac{4 \cdot \left(\left(b \cdot a\right) \cdot \left(b \cdot \left(-a\right)\right)\right)}{{\left(x-scale \cdot y-scale\right)}^{2}}} \]
                            2. Add Preprocessing
                            3. Taylor expanded in y-scale around inf

                              \[\leadsto \color{blue}{\frac{1}{4} \cdot \left(\left(x-scale \cdot \sqrt{8}\right) \cdot \sqrt{\left({a}^{2} \cdot {\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} + {b}^{2} \cdot {\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}\right) - \frac{1}{2} \cdot \frac{-2 \cdot \left(\left({a}^{2} \cdot {\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} + {b}^{2} \cdot {\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}\right) \cdot \left(\frac{{a}^{2} \cdot {\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{x-scale}^{2}} + \frac{{b}^{2} \cdot {\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{x-scale}^{2}}\right)\right) + 4 \cdot \frac{{\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} \cdot \left({\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} \cdot {\left({b}^{2} - {a}^{2}\right)}^{2}\right)}{{x-scale}^{2}}}{\frac{{a}^{2} \cdot {\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{x-scale}^{2}} + \frac{{b}^{2} \cdot {\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{x-scale}^{2}}}}\right)} \]
                            4. Applied rewrites1.5%

                              \[\leadsto \color{blue}{0.25 \cdot \left(\left(x-scale \cdot \sqrt{8}\right) \cdot \sqrt{\left({\left(a \cdot \cos \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)\right)}^{2} + {\left(b \cdot \sin \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)\right)}^{2}\right) - 0.5 \cdot \frac{\mathsf{fma}\left(-2, \left({\left(a \cdot \cos \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)\right)}^{2} + {\left(b \cdot \sin \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)\right)}^{2}\right) \cdot \frac{{\left(a \cdot \sin \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)\right)}^{2} + {\left(b \cdot \cos \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)\right)}^{2}}{x-scale \cdot x-scale}, 4 \cdot \frac{{\cos \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} \cdot {\left(\sin \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \left(b \cdot b - a \cdot a\right)\right)}^{2}}{x-scale \cdot x-scale}\right)}{\frac{{\left(a \cdot \sin \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)\right)}^{2} + {\left(b \cdot \cos \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)\right)}^{2}}{x-scale \cdot x-scale}}}\right)} \]
                            5. Taylor expanded in angle around 0

                              \[\leadsto \frac{1}{4} \cdot \left(\left(x-scale \cdot \sqrt{8}\right) \cdot \sqrt{{a}^{2} - -1 \cdot {a}^{2}}\right) \]
                            6. Step-by-step derivation
                              1. lower--.f64N/A

                                \[\leadsto \frac{1}{4} \cdot \left(\left(x-scale \cdot \sqrt{8}\right) \cdot \sqrt{{a}^{2} - -1 \cdot {a}^{2}}\right) \]
                              2. pow2N/A

                                \[\leadsto \frac{1}{4} \cdot \left(\left(x-scale \cdot \sqrt{8}\right) \cdot \sqrt{a \cdot a - -1 \cdot {a}^{2}}\right) \]
                              3. lift-*.f64N/A

                                \[\leadsto \frac{1}{4} \cdot \left(\left(x-scale \cdot \sqrt{8}\right) \cdot \sqrt{a \cdot a - -1 \cdot {a}^{2}}\right) \]
                              4. lower-*.f64N/A

                                \[\leadsto \frac{1}{4} \cdot \left(\left(x-scale \cdot \sqrt{8}\right) \cdot \sqrt{a \cdot a - -1 \cdot {a}^{2}}\right) \]
                              5. pow2N/A

                                \[\leadsto \frac{1}{4} \cdot \left(\left(x-scale \cdot \sqrt{8}\right) \cdot \sqrt{a \cdot a - -1 \cdot \left(a \cdot a\right)}\right) \]
                              6. lift-*.f6425.0

                                \[\leadsto 0.25 \cdot \left(\left(x-scale \cdot \sqrt{8}\right) \cdot \sqrt{a \cdot a - -1 \cdot \left(a \cdot a\right)}\right) \]
                            7. Applied rewrites25.0%

                              \[\leadsto 0.25 \cdot \left(\left(x-scale \cdot \sqrt{8}\right) \cdot \sqrt{a \cdot a - -1 \cdot \left(a \cdot a\right)}\right) \]
                            8. Taylor expanded in a around 0

                              \[\leadsto \frac{1}{4} \cdot \left(\left(x-scale \cdot \sqrt{8}\right) \cdot \sqrt{2 \cdot {a}^{2}}\right) \]
                            9. Step-by-step derivation
                              1. lower-*.f64N/A

                                \[\leadsto \frac{1}{4} \cdot \left(\left(x-scale \cdot \sqrt{8}\right) \cdot \sqrt{2 \cdot {a}^{2}}\right) \]
                              2. pow2N/A

                                \[\leadsto \frac{1}{4} \cdot \left(\left(x-scale \cdot \sqrt{8}\right) \cdot \sqrt{2 \cdot \left(a \cdot a\right)}\right) \]
                              3. lift-*.f6425.0

                                \[\leadsto 0.25 \cdot \left(\left(x-scale \cdot \sqrt{8}\right) \cdot \sqrt{2 \cdot \left(a \cdot a\right)}\right) \]
                            10. Applied rewrites25.0%

                              \[\leadsto 0.25 \cdot \left(\left(x-scale \cdot \sqrt{8}\right) \cdot \sqrt{2 \cdot \left(a \cdot a\right)}\right) \]

                            if 7.5000000000000003e-28 < y-scale

                            1. Initial program 0.1%

                              \[\frac{-\sqrt{\left(\left(2 \cdot \frac{4 \cdot \left(\left(b \cdot a\right) \cdot \left(b \cdot \left(-a\right)\right)\right)}{{\left(x-scale \cdot y-scale\right)}^{2}}\right) \cdot \left(\left(b \cdot a\right) \cdot \left(b \cdot \left(-a\right)\right)\right)\right) \cdot \left(\left(\frac{\frac{{\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}}{x-scale}}{x-scale} + \frac{\frac{{\left(a \cdot \cos \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} + {\left(b \cdot \sin \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale}}{y-scale}\right) - \sqrt{{\left(\frac{\frac{{\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}}{x-scale}}{x-scale} - \frac{\frac{{\left(a \cdot \cos \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} + {\left(b \cdot \sin \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale}}{y-scale}\right)}^{2} + {\left(\frac{\frac{\left(\left(2 \cdot \left({b}^{2} - {a}^{2}\right)\right) \cdot \sin \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \cos \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)}{x-scale}}{y-scale}\right)}^{2}}\right)}}{\frac{4 \cdot \left(\left(b \cdot a\right) \cdot \left(b \cdot \left(-a\right)\right)\right)}{{\left(x-scale \cdot y-scale\right)}^{2}}} \]
                            2. Add Preprocessing
                            3. Taylor expanded in angle around 0

                              \[\leadsto \color{blue}{\frac{1}{4} \cdot \left(a \cdot \left(x-scale \cdot \left(\sqrt{2} \cdot \sqrt{8}\right)\right)\right)} \]
                            4. Step-by-step derivation
                              1. Applied rewrites25.7%

                                \[\leadsto \color{blue}{0.25 \cdot \left(a \cdot \left(x-scale \cdot 4\right)\right)} \]
                              2. Taylor expanded in a around 0

                                \[\leadsto a \cdot \color{blue}{x-scale} \]
                              3. Step-by-step derivation
                                1. lower-*.f6425.7

                                  \[\leadsto a \cdot x-scale \]
                              4. Applied rewrites25.7%

                                \[\leadsto a \cdot \color{blue}{x-scale} \]
                            5. Recombined 2 regimes into one program.
                            6. Add Preprocessing

                            Alternative 7: 32.1% accurate, 484.7× speedup?

                            \[\begin{array}{l} a_m = \left|a\right| \\ x-scale_m = \left|x-scale\right| \\ y-scale_m = \left|y-scale\right| \\ a\_m \cdot x-scale\_m \end{array} \]
                            a_m = (fabs.f64 a)
                            x-scale_m = (fabs.f64 x-scale)
                            y-scale_m = (fabs.f64 y-scale)
                            (FPCore (a_m b angle x-scale_m y-scale_m)
                             :precision binary64
                             (* a_m x-scale_m))
                            a_m = fabs(a);
                            x-scale_m = fabs(x_45_scale);
                            y-scale_m = fabs(y_45_scale);
                            double code(double a_m, double b, double angle, double x_45_scale_m, double y_45_scale_m) {
                            	return a_m * x_45_scale_m;
                            }
                            
                            a_m =     private
                            x-scale_m =     private
                            y-scale_m =     private
                            module fmin_fmax_functions
                                implicit none
                                private
                                public fmax
                                public fmin
                            
                                interface fmax
                                    module procedure fmax88
                                    module procedure fmax44
                                    module procedure fmax84
                                    module procedure fmax48
                                end interface
                                interface fmin
                                    module procedure fmin88
                                    module procedure fmin44
                                    module procedure fmin84
                                    module procedure fmin48
                                end interface
                            contains
                                real(8) function fmax88(x, y) result (res)
                                    real(8), intent (in) :: x
                                    real(8), intent (in) :: y
                                    res = merge(y, merge(x, max(x, y), y /= y), x /= x)
                                end function
                                real(4) function fmax44(x, y) result (res)
                                    real(4), intent (in) :: x
                                    real(4), intent (in) :: y
                                    res = merge(y, merge(x, max(x, y), y /= y), x /= x)
                                end function
                                real(8) function fmax84(x, y) result(res)
                                    real(8), intent (in) :: x
                                    real(4), intent (in) :: y
                                    res = merge(dble(y), merge(x, max(x, dble(y)), y /= y), x /= x)
                                end function
                                real(8) function fmax48(x, y) result(res)
                                    real(4), intent (in) :: x
                                    real(8), intent (in) :: y
                                    res = merge(y, merge(dble(x), max(dble(x), y), y /= y), x /= x)
                                end function
                                real(8) function fmin88(x, y) result (res)
                                    real(8), intent (in) :: x
                                    real(8), intent (in) :: y
                                    res = merge(y, merge(x, min(x, y), y /= y), x /= x)
                                end function
                                real(4) function fmin44(x, y) result (res)
                                    real(4), intent (in) :: x
                                    real(4), intent (in) :: y
                                    res = merge(y, merge(x, min(x, y), y /= y), x /= x)
                                end function
                                real(8) function fmin84(x, y) result(res)
                                    real(8), intent (in) :: x
                                    real(4), intent (in) :: y
                                    res = merge(dble(y), merge(x, min(x, dble(y)), y /= y), x /= x)
                                end function
                                real(8) function fmin48(x, y) result(res)
                                    real(4), intent (in) :: x
                                    real(8), intent (in) :: y
                                    res = merge(y, merge(dble(x), min(dble(x), y), y /= y), x /= x)
                                end function
                            end module
                            
                            real(8) function code(a_m, b, angle, x_45scale_m, y_45scale_m)
                            use fmin_fmax_functions
                                real(8), intent (in) :: a_m
                                real(8), intent (in) :: b
                                real(8), intent (in) :: angle
                                real(8), intent (in) :: x_45scale_m
                                real(8), intent (in) :: y_45scale_m
                                code = a_m * x_45scale_m
                            end function
                            
                            a_m = Math.abs(a);
                            x-scale_m = Math.abs(x_45_scale);
                            y-scale_m = Math.abs(y_45_scale);
                            public static double code(double a_m, double b, double angle, double x_45_scale_m, double y_45_scale_m) {
                            	return a_m * x_45_scale_m;
                            }
                            
                            a_m = math.fabs(a)
                            x-scale_m = math.fabs(x_45_scale)
                            y-scale_m = math.fabs(y_45_scale)
                            def code(a_m, b, angle, x_45_scale_m, y_45_scale_m):
                            	return a_m * x_45_scale_m
                            
                            a_m = abs(a)
                            x-scale_m = abs(x_45_scale)
                            y-scale_m = abs(y_45_scale)
                            function code(a_m, b, angle, x_45_scale_m, y_45_scale_m)
                            	return Float64(a_m * x_45_scale_m)
                            end
                            
                            a_m = abs(a);
                            x-scale_m = abs(x_45_scale);
                            y-scale_m = abs(y_45_scale);
                            function tmp = code(a_m, b, angle, x_45_scale_m, y_45_scale_m)
                            	tmp = a_m * x_45_scale_m;
                            end
                            
                            a_m = N[Abs[a], $MachinePrecision]
                            x-scale_m = N[Abs[x$45$scale], $MachinePrecision]
                            y-scale_m = N[Abs[y$45$scale], $MachinePrecision]
                            code[a$95$m_, b_, angle_, x$45$scale$95$m_, y$45$scale$95$m_] := N[(a$95$m * x$45$scale$95$m), $MachinePrecision]
                            
                            \begin{array}{l}
                            a_m = \left|a\right|
                            \\
                            x-scale_m = \left|x-scale\right|
                            \\
                            y-scale_m = \left|y-scale\right|
                            
                            \\
                            a\_m \cdot x-scale\_m
                            \end{array}
                            
                            Derivation
                            1. Initial program 0.1%

                              \[\frac{-\sqrt{\left(\left(2 \cdot \frac{4 \cdot \left(\left(b \cdot a\right) \cdot \left(b \cdot \left(-a\right)\right)\right)}{{\left(x-scale \cdot y-scale\right)}^{2}}\right) \cdot \left(\left(b \cdot a\right) \cdot \left(b \cdot \left(-a\right)\right)\right)\right) \cdot \left(\left(\frac{\frac{{\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}}{x-scale}}{x-scale} + \frac{\frac{{\left(a \cdot \cos \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} + {\left(b \cdot \sin \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale}}{y-scale}\right) - \sqrt{{\left(\frac{\frac{{\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}}{x-scale}}{x-scale} - \frac{\frac{{\left(a \cdot \cos \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} + {\left(b \cdot \sin \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{y-scale}}{y-scale}\right)}^{2} + {\left(\frac{\frac{\left(\left(2 \cdot \left({b}^{2} - {a}^{2}\right)\right) \cdot \sin \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \cos \left(\frac{angle}{180} \cdot \mathsf{PI}\left(\right)\right)}{x-scale}}{y-scale}\right)}^{2}}\right)}}{\frac{4 \cdot \left(\left(b \cdot a\right) \cdot \left(b \cdot \left(-a\right)\right)\right)}{{\left(x-scale \cdot y-scale\right)}^{2}}} \]
                            2. Add Preprocessing
                            3. Taylor expanded in angle around 0

                              \[\leadsto \color{blue}{\frac{1}{4} \cdot \left(a \cdot \left(x-scale \cdot \left(\sqrt{2} \cdot \sqrt{8}\right)\right)\right)} \]
                            4. Step-by-step derivation
                              1. Applied rewrites23.4%

                                \[\leadsto \color{blue}{0.25 \cdot \left(a \cdot \left(x-scale \cdot 4\right)\right)} \]
                              2. Taylor expanded in a around 0

                                \[\leadsto a \cdot \color{blue}{x-scale} \]
                              3. Step-by-step derivation
                                1. lower-*.f6423.4

                                  \[\leadsto a \cdot x-scale \]
                              4. Applied rewrites23.4%

                                \[\leadsto a \cdot \color{blue}{x-scale} \]
                              5. Add Preprocessing

                              Reproduce

                              ?
                              herbie shell --seed 2025051 
                              (FPCore (a b angle x-scale y-scale)
                                :name "b from scale-rotated-ellipse"
                                :precision binary64
                                (/ (- (sqrt (* (* (* 2.0 (/ (* 4.0 (* (* b a) (* b (- a)))) (pow (* x-scale y-scale) 2.0))) (* (* b a) (* b (- a)))) (- (+ (/ (/ (+ (pow (* a (sin (* (/ angle 180.0) (PI)))) 2.0) (pow (* b (cos (* (/ angle 180.0) (PI)))) 2.0)) x-scale) x-scale) (/ (/ (+ (pow (* a (cos (* (/ angle 180.0) (PI)))) 2.0) (pow (* b (sin (* (/ angle 180.0) (PI)))) 2.0)) y-scale) y-scale)) (sqrt (+ (pow (- (/ (/ (+ (pow (* a (sin (* (/ angle 180.0) (PI)))) 2.0) (pow (* b (cos (* (/ angle 180.0) (PI)))) 2.0)) x-scale) x-scale) (/ (/ (+ (pow (* a (cos (* (/ angle 180.0) (PI)))) 2.0) (pow (* b (sin (* (/ angle 180.0) (PI)))) 2.0)) y-scale) y-scale)) 2.0) (pow (/ (/ (* (* (* 2.0 (- (pow b 2.0) (pow a 2.0))) (sin (* (/ angle 180.0) (PI)))) (cos (* (/ angle 180.0) (PI)))) x-scale) y-scale) 2.0))))))) (/ (* 4.0 (* (* b a) (* b (- a)))) (pow (* x-scale y-scale) 2.0))))