b from scale-rotated-ellipse

Percentage Accurate: 0.1% → 43.3%

Time: 29.5s

Alternatives: 5

Speedup: 484.7×

• could not determine a ground truth

  1. a = -2.6529768731522044e+307
  2. b = -5.595901815053105e-304
  3. angle = 1.5492735553047592e+141
  4. x-scale = -3.4670173379226445e+290
  5. y-scale = 6.209880176400302e-115

Specification

Math

\[\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

(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)))

Initial Program: 0.1% accurate, 1.0× speedup

Math

\[\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

(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)))

Alternative 1: 43.3% accurate, 10.9× speedup

Math

\[\begin{array}{l} x-scale_m = \left|x-scale\right| \\ b_m = \left|b\right| \\ a_m = \left|a\right| \\ \begin{array}{l} \mathbf{if}\;b\_m \leq 5.2 \cdot 10^{-75}:\\ \;\;\;\;\sqrt{\left({\sin \left(\left(\mathsf{PI}\left(\right) \cdot angle\right) \cdot 0.005555555555555556\right)}^{2} \cdot \left(b\_m \cdot b\_m\right)\right) \cdot 0} \cdot \left(\left(\sqrt{8} \cdot x-scale\_m\right) \cdot 0.25\right)\\ \mathbf{else}:\\ \;\;\;\;x-scale\_m \cdot a\_m\\ \end{array} \end{array} \]

FPCore

x-scale_m = (fabs.f64 x-scale)
b_m = (fabs.f64 b)
a_m = (fabs.f64 a)
(FPCore (a_m b_m angle x-scale_m y-scale)
 :precision binary64
 (if (<= b_m 5.2e-75)
   (*
    (sqrt
     (*
      (* (pow (sin (* (* (PI) angle) 0.005555555555555556)) 2.0) (* b_m b_m))
      0.0))
    (* (* (sqrt 8.0) x-scale_m) 0.25))
   (* x-scale_m a_m)))

Derivation

1. Split input into 2 regimes

2. if b < 5.2e-75

   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 x-scale around 0

      \[\leadsto \color{blue}{\frac{1}{4} \cdot \left(\left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right) \cdot \sqrt{\left(\frac{{a}^{2} \cdot {\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{y-scale}^{2}} + \frac{{b}^{2} \cdot {\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{y-scale}^{2}}\right) - \frac{1}{2} \cdot \frac{-2 \cdot \left(\left({a}^{2} \cdot {\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} + {b}^{2} \cdot {\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}\right) \cdot \left(\frac{{a}^{2} \cdot {\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{y-scale}^{2}} + \frac{{b}^{2} \cdot {\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{y-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)}{{y-scale}^{2}}}{{a}^{2} \cdot {\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} + {b}^{2} \cdot {\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}}\right)} \]

   5. Applied rewrites 4.9%

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

   6. Taylor expanded in a around 0

      \[\leadsto \frac{1}{4} \cdot \color{blue}{\left(\left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right) \cdot \sqrt{\frac{-1}{2} \cdot \frac{-2 \cdot \frac{{b}^{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)}{{y-scale}^{2}} + 4 \cdot \frac{{b}^{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)}{{y-scale}^{2}}}{{b}^{2} \cdot {\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}} + \frac{{b}^{2} \cdot {\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{y-scale}^{2}}}\right)} \]

   8. Applied rewrites 3.9%

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

   9. Taylor expanded in y-scale around 0

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

   11. Applied rewrites 37.5%

       \[\leadsto \left(0.25 \cdot \left(x-scale \cdot \sqrt{8}\right)\right) \cdot \sqrt{0 \cdot \left(\left(b \cdot b\right) \cdot {\sin \left(0.005555555555555556 \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}\right)} \]

   if 5.2e-75 < b

   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. associate-*r* N/A

         \[\leadsto \color{blue}{\left(\frac{1}{4} \cdot a\right) \cdot \left(x-scale \cdot \left(\sqrt{2} \cdot \sqrt{8}\right)\right)} \]

      2. lower-*.f64 N/A

         \[\leadsto \color{blue}{\left(\frac{1}{4} \cdot a\right) \cdot \left(x-scale \cdot \left(\sqrt{2} \cdot \sqrt{8}\right)\right)} \]

      3. lower-*.f64 N/A

         \[\leadsto \color{blue}{\left(\frac{1}{4} \cdot a\right)} \cdot \left(x-scale \cdot \left(\sqrt{2} \cdot \sqrt{8}\right)\right) \]

      4. associate-*r* N/A

         \[\leadsto \left(\frac{1}{4} \cdot a\right) \cdot \color{blue}{\left(\left(x-scale \cdot \sqrt{2}\right) \cdot \sqrt{8}\right)} \]

      5. lower-*.f64 N/A

         \[\leadsto \left(\frac{1}{4} \cdot a\right) \cdot \color{blue}{\left(\left(x-scale \cdot \sqrt{2}\right) \cdot \sqrt{8}\right)} \]

      6. lower-*.f64 N/A

         \[\leadsto \left(\frac{1}{4} \cdot a\right) \cdot \left(\color{blue}{\left(x-scale \cdot \sqrt{2}\right)} \cdot \sqrt{8}\right) \]

      7. lower-sqrt.f64 N/A

         \[\leadsto \left(\frac{1}{4} \cdot a\right) \cdot \left(\left(x-scale \cdot \color{blue}{\sqrt{2}}\right) \cdot \sqrt{8}\right) \]

      8. lower-sqrt.f64 25.8

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

   5. Applied rewrites 25.8%

      \[\leadsto \color{blue}{\left(0.25 \cdot a\right) \cdot \left(\left(x-scale \cdot \sqrt{2}\right) \cdot \sqrt{8}\right)} \]
   6. Step-by-step derivation

   7. Applied rewrites 25.9%

      \[\leadsto x-scale \cdot \color{blue}{\left(4 \cdot \left(0.25 \cdot a\right)\right)} \]
   8. Step-by-step derivation

   9. Applied rewrites 25.9%

      \[\leadsto x-scale \cdot a \]
3. Recombined 2 regimes into one program.

4. Final simplification 34.3%

   \[\leadsto \begin{array}{l} \mathbf{if}\;b \leq 5.2 \cdot 10^{-75}:\\ \;\;\;\;\sqrt{\left({\sin \left(\left(\mathsf{PI}\left(\right) \cdot angle\right) \cdot 0.005555555555555556\right)}^{2} \cdot \left(b \cdot b\right)\right) \cdot 0} \cdot \left(\left(\sqrt{8} \cdot x-scale\right) \cdot 0.25\right)\\ \mathbf{else}:\\ \;\;\;\;x-scale \cdot a\\ \end{array} \]
5. Add Preprocessing

Alternative 2: 37.8% accurate, 11.1× speedup

Math

\[\begin{array}{l} x-scale_m = \left|x-scale\right| \\ b_m = \left|b\right| \\ a_m = \left|a\right| \\ \begin{array}{l} \mathbf{if}\;y-scale \leq 6.2 \cdot 10^{-42}:\\ \;\;\;\;\sqrt{{\sin \left(\left(\mathsf{PI}\left(\right) \cdot angle\right) \cdot 0.005555555555555556\right)}^{2} \cdot 0} \cdot \left(\left(\left(\sqrt{8} \cdot y-scale\right) \cdot a\_m\right) \cdot 0.25\right)\\ \mathbf{else}:\\ \;\;\;\;x-scale\_m \cdot a\_m\\ \end{array} \end{array} \]

FPCore

x-scale_m = (fabs.f64 x-scale)
b_m = (fabs.f64 b)
a_m = (fabs.f64 a)
(FPCore (a_m b_m angle x-scale_m y-scale)
 :precision binary64
 (if (<= y-scale 6.2e-42)
   (*
    (sqrt (* (pow (sin (* (* (PI) angle) 0.005555555555555556)) 2.0) 0.0))
    (* (* (* (sqrt 8.0) y-scale) a_m) 0.25))
   (* x-scale_m a_m)))

Derivation

1. Split input into 2 regimes

2. if y-scale < 6.2000000000000005e-42

   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 b around 0

      \[\leadsto \color{blue}{\frac{1}{4} \cdot \left(\left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right) \cdot \sqrt{\left(\frac{{a}^{2} \cdot {\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{y-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}}\right) - \sqrt{4 \cdot \frac{{a}^{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)}{{x-scale}^{2} \cdot {y-scale}^{2}} + {\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{{a}^{2} \cdot {\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{y-scale}^{2}}\right)}^{2}}}\right)} \]

   5. Applied rewrites 4.9%

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

   6. Taylor expanded in x-scale around inf

      \[\leadsto \frac{1}{4} \cdot \color{blue}{\left(\left(y-scale \cdot \sqrt{8}\right) \cdot \sqrt{{a}^{2} \cdot {\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} - \frac{1}{2} \cdot \frac{{y-scale}^{2} \cdot \left(-2 \cdot \frac{{a}^{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)}{{y-scale}^{2}} + 4 \cdot \frac{{a}^{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)}{{y-scale}^{2}}\right)}{{a}^{2} \cdot {\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}}\right)} \]

   8. Applied rewrites 2.1%

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

   9. Taylor expanded in a around 0

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

   11. Applied rewrites 38.2%

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

   if 6.2000000000000005e-42 < 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. associate-*r* N/A

         \[\leadsto \color{blue}{\left(\frac{1}{4} \cdot a\right) \cdot \left(x-scale \cdot \left(\sqrt{2} \cdot \sqrt{8}\right)\right)} \]

      2. lower-*.f64 N/A

         \[\leadsto \color{blue}{\left(\frac{1}{4} \cdot a\right) \cdot \left(x-scale \cdot \left(\sqrt{2} \cdot \sqrt{8}\right)\right)} \]

      3. lower-*.f64 N/A

         \[\leadsto \color{blue}{\left(\frac{1}{4} \cdot a\right)} \cdot \left(x-scale \cdot \left(\sqrt{2} \cdot \sqrt{8}\right)\right) \]

      4. associate-*r* N/A

         \[\leadsto \left(\frac{1}{4} \cdot a\right) \cdot \color{blue}{\left(\left(x-scale \cdot \sqrt{2}\right) \cdot \sqrt{8}\right)} \]

      5. lower-*.f64 N/A

         \[\leadsto \left(\frac{1}{4} \cdot a\right) \cdot \color{blue}{\left(\left(x-scale \cdot \sqrt{2}\right) \cdot \sqrt{8}\right)} \]

      6. lower-*.f64 N/A

         \[\leadsto \left(\frac{1}{4} \cdot a\right) \cdot \left(\color{blue}{\left(x-scale \cdot \sqrt{2}\right)} \cdot \sqrt{8}\right) \]

      7. lower-sqrt.f64 N/A

         \[\leadsto \left(\frac{1}{4} \cdot a\right) \cdot \left(\left(x-scale \cdot \color{blue}{\sqrt{2}}\right) \cdot \sqrt{8}\right) \]

      8. lower-sqrt.f64 25.7

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

   5. Applied rewrites 25.7%

      \[\leadsto \color{blue}{\left(0.25 \cdot a\right) \cdot \left(\left(x-scale \cdot \sqrt{2}\right) \cdot \sqrt{8}\right)} \]
   6. Step-by-step derivation

   7. Applied rewrites 25.9%

      \[\leadsto x-scale \cdot \color{blue}{\left(4 \cdot \left(0.25 \cdot a\right)\right)} \]
   8. Step-by-step derivation

   9. Applied rewrites 25.9%

      \[\leadsto x-scale \cdot a \]
3. Recombined 2 regimes into one program.

4. Final simplification 34.8%

   \[\leadsto \begin{array}{l} \mathbf{if}\;y-scale \leq 6.2 \cdot 10^{-42}:\\ \;\;\;\;\sqrt{{\sin \left(\left(\mathsf{PI}\left(\right) \cdot angle\right) \cdot 0.005555555555555556\right)}^{2} \cdot 0} \cdot \left(\left(\left(\sqrt{8} \cdot y-scale\right) \cdot a\right) \cdot 0.25\right)\\ \mathbf{else}:\\ \;\;\;\;x-scale \cdot a\\ \end{array} \]
5. Add Preprocessing

Alternative 3: 32.9% accurate, 30.9× speedup

Math

\[\begin{array}{l} x-scale_m = \left|x-scale\right| \\ b_m = \left|b\right| \\ a_m = \left|a\right| \\ \begin{array}{l} t_0 := \sqrt{8} \cdot y-scale\\ t_1 := \mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\\ \mathbf{if}\;y-scale \leq 2 \cdot 10^{-147}:\\ \;\;\;\;\sqrt{\left(t\_1 \cdot \left(a\_m \cdot a\_m\right)\right) \cdot 0} \cdot \left(\left(t\_0 \cdot angle\right) \cdot 0.25\right)\\ \mathbf{elif}\;y-scale \leq 5.3 \cdot 10^{-92}:\\ \;\;\;\;\sqrt{\frac{t\_1 \cdot \left(b\_m \cdot b\_m\right)}{y-scale \cdot y-scale} \cdot 0} \cdot \left(\left(\left(t\_0 \cdot x-scale\_m\right) \cdot angle\right) \cdot 0.25\right)\\ \mathbf{else}:\\ \;\;\;\;x-scale\_m \cdot a\_m\\ \end{array} \end{array} \]

FPCore

x-scale_m = (fabs.f64 x-scale)
b_m = (fabs.f64 b)
a_m = (fabs.f64 a)
(FPCore (a_m b_m angle x-scale_m y-scale)
 :precision binary64
 (let* ((t_0 (* (sqrt 8.0) y-scale)) (t_1 (* (PI) (PI))))
   (if (<= y-scale 2e-147)
     (* (sqrt (* (* t_1 (* a_m a_m)) 0.0)) (* (* t_0 angle) 0.25))
     (if (<= y-scale 5.3e-92)
       (*
        (sqrt (* (/ (* t_1 (* b_m b_m)) (* y-scale y-scale)) 0.0))
        (* (* (* t_0 x-scale_m) angle) 0.25))
       (* x-scale_m a_m)))))

Derivation

1. Split input into 3 regimes

2. if y-scale < 1.9999999999999999e-147

   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 b around 0

      \[\leadsto \color{blue}{\frac{1}{4} \cdot \left(\left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right) \cdot \sqrt{\left(\frac{{a}^{2} \cdot {\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{y-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}}\right) - \sqrt{4 \cdot \frac{{a}^{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)}{{x-scale}^{2} \cdot {y-scale}^{2}} + {\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{{a}^{2} \cdot {\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{y-scale}^{2}}\right)}^{2}}}\right)} \]

   5. Applied rewrites 5.5%

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

   6. Taylor expanded in x-scale around inf

      \[\leadsto \frac{1}{4} \cdot \color{blue}{\left(\left(y-scale \cdot \sqrt{8}\right) \cdot \sqrt{{a}^{2} \cdot {\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} - \frac{1}{2} \cdot \frac{{y-scale}^{2} \cdot \left(-2 \cdot \frac{{a}^{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)}{{y-scale}^{2}} + 4 \cdot \frac{{a}^{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)}{{y-scale}^{2}}\right)}{{a}^{2} \cdot {\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}}\right)} \]

   8. Applied rewrites 1.0%

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

   9. Taylor expanded in angle around 0

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

   11. Applied rewrites 30.6%

       \[\leadsto \left(0.25 \cdot \left(angle \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\left(\left(a \cdot a\right) \cdot \left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right)\right) \cdot 0} \]

   if 1.9999999999999999e-147 < y-scale < 5.30000000000000029e-92

   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 x-scale around 0

      \[\leadsto \color{blue}{\frac{1}{4} \cdot \left(\left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right) \cdot \sqrt{\left(\frac{{a}^{2} \cdot {\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{y-scale}^{2}} + \frac{{b}^{2} \cdot {\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{y-scale}^{2}}\right) - \frac{1}{2} \cdot \frac{-2 \cdot \left(\left({a}^{2} \cdot {\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} + {b}^{2} \cdot {\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}\right) \cdot \left(\frac{{a}^{2} \cdot {\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{y-scale}^{2}} + \frac{{b}^{2} \cdot {\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{y-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)}{{y-scale}^{2}}}{{a}^{2} \cdot {\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} + {b}^{2} \cdot {\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}}\right)} \]

   5. Applied rewrites 1.2%

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

   6. Taylor expanded in a around 0

      \[\leadsto \frac{1}{4} \cdot \color{blue}{\left(\left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right) \cdot \sqrt{\frac{-1}{2} \cdot \frac{-2 \cdot \frac{{b}^{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)}{{y-scale}^{2}} + 4 \cdot \frac{{b}^{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)}{{y-scale}^{2}}}{{b}^{2} \cdot {\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}} + \frac{{b}^{2} \cdot {\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{y-scale}^{2}}}\right)} \]

   8. Applied rewrites 1.5%

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

   9. Taylor expanded in angle around 0

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

   11. Applied rewrites 34.3%

       \[\leadsto \left(0.25 \cdot \left(angle \cdot \left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right)\right) \cdot \sqrt{\frac{\left(b \cdot b\right) \cdot \left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right)}{y-scale \cdot y-scale} \cdot 0} \]

   if 5.30000000000000029e-92 < 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. associate-*r* N/A

         \[\leadsto \color{blue}{\left(\frac{1}{4} \cdot a\right) \cdot \left(x-scale \cdot \left(\sqrt{2} \cdot \sqrt{8}\right)\right)} \]

      2. lower-*.f64 N/A

         \[\leadsto \color{blue}{\left(\frac{1}{4} \cdot a\right) \cdot \left(x-scale \cdot \left(\sqrt{2} \cdot \sqrt{8}\right)\right)} \]

      3. lower-*.f64 N/A

         \[\leadsto \color{blue}{\left(\frac{1}{4} \cdot a\right)} \cdot \left(x-scale \cdot \left(\sqrt{2} \cdot \sqrt{8}\right)\right) \]

      4. associate-*r* N/A

         \[\leadsto \left(\frac{1}{4} \cdot a\right) \cdot \color{blue}{\left(\left(x-scale \cdot \sqrt{2}\right) \cdot \sqrt{8}\right)} \]

      5. lower-*.f64 N/A

         \[\leadsto \left(\frac{1}{4} \cdot a\right) \cdot \color{blue}{\left(\left(x-scale \cdot \sqrt{2}\right) \cdot \sqrt{8}\right)} \]

      6. lower-*.f64 N/A

         \[\leadsto \left(\frac{1}{4} \cdot a\right) \cdot \left(\color{blue}{\left(x-scale \cdot \sqrt{2}\right)} \cdot \sqrt{8}\right) \]

      7. lower-sqrt.f64 N/A

         \[\leadsto \left(\frac{1}{4} \cdot a\right) \cdot \left(\left(x-scale \cdot \color{blue}{\sqrt{2}}\right) \cdot \sqrt{8}\right) \]

      8. lower-sqrt.f64 27.1

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

   5. Applied rewrites 27.1%

      \[\leadsto \color{blue}{\left(0.25 \cdot a\right) \cdot \left(\left(x-scale \cdot \sqrt{2}\right) \cdot \sqrt{8}\right)} \]
   6. Step-by-step derivation

   7. Applied rewrites 27.2%

      \[\leadsto x-scale \cdot \color{blue}{\left(4 \cdot \left(0.25 \cdot a\right)\right)} \]
   8. Step-by-step derivation

   9. Applied rewrites 27.2%

      \[\leadsto x-scale \cdot a \]
3. Recombined 3 regimes into one program.

4. Final simplification 29.8%

   \[\leadsto \begin{array}{l} \mathbf{if}\;y-scale \leq 2 \cdot 10^{-147}:\\ \;\;\;\;\sqrt{\left(\left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right) \cdot \left(a \cdot a\right)\right) \cdot 0} \cdot \left(\left(\left(\sqrt{8} \cdot y-scale\right) \cdot angle\right) \cdot 0.25\right)\\ \mathbf{elif}\;y-scale \leq 5.3 \cdot 10^{-92}:\\ \;\;\;\;\sqrt{\frac{\left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right) \cdot \left(b \cdot b\right)}{y-scale \cdot y-scale} \cdot 0} \cdot \left(\left(\left(\left(\sqrt{8} \cdot y-scale\right) \cdot x-scale\right) \cdot angle\right) \cdot 0.25\right)\\ \mathbf{else}:\\ \;\;\;\;x-scale \cdot a\\ \end{array} \]
5. Add Preprocessing

Alternative 4: 32.7% accurate, 43.4× speedup

Math

\[\begin{array}{l} x-scale_m = \left|x-scale\right| \\ b_m = \left|b\right| \\ a_m = \left|a\right| \\ \begin{array}{l} \mathbf{if}\;y-scale \leq 2.2 \cdot 10^{-147}:\\ \;\;\;\;\sqrt{\left(\left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right) \cdot \left(a\_m \cdot a\_m\right)\right) \cdot 0} \cdot \left(\left(\left(\sqrt{8} \cdot y-scale\right) \cdot angle\right) \cdot 0.25\right)\\ \mathbf{else}:\\ \;\;\;\;x-scale\_m \cdot a\_m\\ \end{array} \end{array} \]

FPCore

x-scale_m = (fabs.f64 x-scale)
b_m = (fabs.f64 b)
a_m = (fabs.f64 a)
(FPCore (a_m b_m angle x-scale_m y-scale)
 :precision binary64
 (if (<= y-scale 2.2e-147)
   (*
    (sqrt (* (* (* (PI) (PI)) (* a_m a_m)) 0.0))
    (* (* (* (sqrt 8.0) y-scale) angle) 0.25))
   (* x-scale_m a_m)))

Derivation

1. Split input into 2 regimes

2. if y-scale < 2.2000000000000001e-147

   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 b around 0

      \[\leadsto \color{blue}{\frac{1}{4} \cdot \left(\left(x-scale \cdot \left(y-scale \cdot \sqrt{8}\right)\right) \cdot \sqrt{\left(\frac{{a}^{2} \cdot {\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{y-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}}\right) - \sqrt{4 \cdot \frac{{a}^{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)}{{x-scale}^{2} \cdot {y-scale}^{2}} + {\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{{a}^{2} \cdot {\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}{{y-scale}^{2}}\right)}^{2}}}\right)} \]

   5. Applied rewrites 5.5%

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

   6. Taylor expanded in x-scale around inf

      \[\leadsto \frac{1}{4} \cdot \color{blue}{\left(\left(y-scale \cdot \sqrt{8}\right) \cdot \sqrt{{a}^{2} \cdot {\sin \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2} - \frac{1}{2} \cdot \frac{{y-scale}^{2} \cdot \left(-2 \cdot \frac{{a}^{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)}{{y-scale}^{2}} + 4 \cdot \frac{{a}^{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)}{{y-scale}^{2}}\right)}{{a}^{2} \cdot {\cos \left(\frac{1}{180} \cdot \left(angle \cdot \mathsf{PI}\left(\right)\right)\right)}^{2}}}\right)} \]

   8. Applied rewrites 1.0%

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

   9. Taylor expanded in angle around 0

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

   11. Applied rewrites 30.6%

       \[\leadsto \left(0.25 \cdot \left(angle \cdot \left(y-scale \cdot \sqrt{8}\right)\right)\right) \cdot \sqrt{\left(\left(a \cdot a\right) \cdot \left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right)\right) \cdot 0} \]

   if 2.2000000000000001e-147 < 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. associate-*r* N/A

         \[\leadsto \color{blue}{\left(\frac{1}{4} \cdot a\right) \cdot \left(x-scale \cdot \left(\sqrt{2} \cdot \sqrt{8}\right)\right)} \]

      2. lower-*.f64 N/A

         \[\leadsto \color{blue}{\left(\frac{1}{4} \cdot a\right) \cdot \left(x-scale \cdot \left(\sqrt{2} \cdot \sqrt{8}\right)\right)} \]

      3. lower-*.f64 N/A

         \[\leadsto \color{blue}{\left(\frac{1}{4} \cdot a\right)} \cdot \left(x-scale \cdot \left(\sqrt{2} \cdot \sqrt{8}\right)\right) \]

      4. associate-*r* N/A

         \[\leadsto \left(\frac{1}{4} \cdot a\right) \cdot \color{blue}{\left(\left(x-scale \cdot \sqrt{2}\right) \cdot \sqrt{8}\right)} \]

      5. lower-*.f64 N/A

         \[\leadsto \left(\frac{1}{4} \cdot a\right) \cdot \color{blue}{\left(\left(x-scale \cdot \sqrt{2}\right) \cdot \sqrt{8}\right)} \]

      6. lower-*.f64 N/A

         \[\leadsto \left(\frac{1}{4} \cdot a\right) \cdot \left(\color{blue}{\left(x-scale \cdot \sqrt{2}\right)} \cdot \sqrt{8}\right) \]

      7. lower-sqrt.f64 N/A

         \[\leadsto \left(\frac{1}{4} \cdot a\right) \cdot \left(\left(x-scale \cdot \color{blue}{\sqrt{2}}\right) \cdot \sqrt{8}\right) \]

      8. lower-sqrt.f64 25.1

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

   5. Applied rewrites 25.1%

      \[\leadsto \color{blue}{\left(0.25 \cdot a\right) \cdot \left(\left(x-scale \cdot \sqrt{2}\right) \cdot \sqrt{8}\right)} \]
   6. Step-by-step derivation

   7. Applied rewrites 25.2%

      \[\leadsto x-scale \cdot \color{blue}{\left(4 \cdot \left(0.25 \cdot a\right)\right)} \]
   8. Step-by-step derivation

   9. Applied rewrites 25.2%

      \[\leadsto x-scale \cdot a \]
3. Recombined 2 regimes into one program.

4. Final simplification 28.7%

   \[\leadsto \begin{array}{l} \mathbf{if}\;y-scale \leq 2.2 \cdot 10^{-147}:\\ \;\;\;\;\sqrt{\left(\left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right) \cdot \left(a \cdot a\right)\right) \cdot 0} \cdot \left(\left(\left(\sqrt{8} \cdot y-scale\right) \cdot angle\right) \cdot 0.25\right)\\ \mathbf{else}:\\ \;\;\;\;x-scale \cdot a\\ \end{array} \]
5. Add Preprocessing

Alternative 5: 32.8% accurate, 484.7× speedup

Math

\[\begin{array}{l} x-scale_m = \left|x-scale\right| \\ b_m = \left|b\right| \\ a_m = \left|a\right| \\ x-scale\_m \cdot a\_m \end{array} \]

FPCore

x-scale_m = (fabs.f64 x-scale)
b_m = (fabs.f64 b)
a_m = (fabs.f64 a)
(FPCore (a_m b_m angle x-scale_m y-scale)
 :precision binary64
 (* x-scale_m a_m))

C

x-scale_m = fabs(x_45_scale);
b_m = fabs(b);
a_m = fabs(a);
double code(double a_m, double b_m, double angle, double x_45_scale_m, double y_45_scale) {
	return x_45_scale_m * a_m;
}

Fortran

x-scale_m = abs(x_45scale)
b_m = abs(b)
a_m = abs(a)
real(8) function code(a_m, b_m, angle, x_45scale_m, y_45scale)
    real(8), intent (in) :: a_m
    real(8), intent (in) :: b_m
    real(8), intent (in) :: angle
    real(8), intent (in) :: x_45scale_m
    real(8), intent (in) :: y_45scale
    code = x_45scale_m * a_m
end function

Java

x-scale_m = Math.abs(x_45_scale);
b_m = Math.abs(b);
a_m = Math.abs(a);
public static double code(double a_m, double b_m, double angle, double x_45_scale_m, double y_45_scale) {
	return x_45_scale_m * a_m;
}

Python

x-scale_m = math.fabs(x_45_scale)
b_m = math.fabs(b)
a_m = math.fabs(a)
def code(a_m, b_m, angle, x_45_scale_m, y_45_scale):
	return x_45_scale_m * a_m

Julia

x-scale_m = abs(x_45_scale)
b_m = abs(b)
a_m = abs(a)
function code(a_m, b_m, angle, x_45_scale_m, y_45_scale)
	return Float64(x_45_scale_m * a_m)
end

MATLAB

x-scale_m = abs(x_45_scale);
b_m = abs(b);
a_m = abs(a);
function tmp = code(a_m, b_m, angle, x_45_scale_m, y_45_scale)
	tmp = x_45_scale_m * a_m;
end

Wolfram

x-scale_m = N[Abs[x$45$scale], $MachinePrecision]
b_m = N[Abs[b], $MachinePrecision]
a_m = N[Abs[a], $MachinePrecision]
code[a$95$m_, b$95$m_, angle_, x$45$scale$95$m_, y$45$scale_] := N[(x$45$scale$95$m * a$95$m), $MachinePrecision]

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. associate-*r* N/A

      \[\leadsto \color{blue}{\left(\frac{1}{4} \cdot a\right) \cdot \left(x-scale \cdot \left(\sqrt{2} \cdot \sqrt{8}\right)\right)} \]

   2. lower-*.f64 N/A

      \[\leadsto \color{blue}{\left(\frac{1}{4} \cdot a\right) \cdot \left(x-scale \cdot \left(\sqrt{2} \cdot \sqrt{8}\right)\right)} \]

   3. lower-*.f64 N/A

      \[\leadsto \color{blue}{\left(\frac{1}{4} \cdot a\right)} \cdot \left(x-scale \cdot \left(\sqrt{2} \cdot \sqrt{8}\right)\right) \]

   4. associate-*r* N/A

      \[\leadsto \left(\frac{1}{4} \cdot a\right) \cdot \color{blue}{\left(\left(x-scale \cdot \sqrt{2}\right) \cdot \sqrt{8}\right)} \]

   5. lower-*.f64 N/A

      \[\leadsto \left(\frac{1}{4} \cdot a\right) \cdot \color{blue}{\left(\left(x-scale \cdot \sqrt{2}\right) \cdot \sqrt{8}\right)} \]

   6. lower-*.f64 N/A

      \[\leadsto \left(\frac{1}{4} \cdot a\right) \cdot \left(\color{blue}{\left(x-scale \cdot \sqrt{2}\right)} \cdot \sqrt{8}\right) \]

   7. lower-sqrt.f64 N/A

      \[\leadsto \left(\frac{1}{4} \cdot a\right) \cdot \left(\left(x-scale \cdot \color{blue}{\sqrt{2}}\right) \cdot \sqrt{8}\right) \]

   8. lower-sqrt.f64 21.5

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

5. Applied rewrites 21.5%

   \[\leadsto \color{blue}{\left(0.25 \cdot a\right) \cdot \left(\left(x-scale \cdot \sqrt{2}\right) \cdot \sqrt{8}\right)} \]
6. Step-by-step derivation

7. Applied rewrites 21.6%

   \[\leadsto x-scale \cdot \color{blue}{\left(4 \cdot \left(0.25 \cdot a\right)\right)} \]
8. Step-by-step derivation

9. Applied rewrites 21.6%

   \[\leadsto x-scale \cdot a \]
10. Add Preprocessing

Reproduce

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