ab-angle->ABCF C

Percentage Accurate: 78.8% → 78.8%
Time: 34.7s
Alternatives: 12
Speedup: 1.0×

Specification

?
\[\begin{array}{l} \\ \begin{array}{l} t_0 := \pi \cdot \frac{angle}{180}\\ {\left(a \cdot \cos t\_0\right)}^{2} + {\left(b \cdot \sin t\_0\right)}^{2} \end{array} \end{array} \]
(FPCore (a b angle)
 :precision binary64
 (let* ((t_0 (* PI (/ angle 180.0))))
   (+ (pow (* a (cos t_0)) 2.0) (pow (* b (sin t_0)) 2.0))))
double code(double a, double b, double angle) {
	double t_0 = ((double) M_PI) * (angle / 180.0);
	return pow((a * cos(t_0)), 2.0) + pow((b * sin(t_0)), 2.0);
}
public static double code(double a, double b, double angle) {
	double t_0 = Math.PI * (angle / 180.0);
	return Math.pow((a * Math.cos(t_0)), 2.0) + Math.pow((b * Math.sin(t_0)), 2.0);
}
def code(a, b, angle):
	t_0 = math.pi * (angle / 180.0)
	return math.pow((a * math.cos(t_0)), 2.0) + math.pow((b * math.sin(t_0)), 2.0)
function code(a, b, angle)
	t_0 = Float64(pi * Float64(angle / 180.0))
	return Float64((Float64(a * cos(t_0)) ^ 2.0) + (Float64(b * sin(t_0)) ^ 2.0))
end
function tmp = code(a, b, angle)
	t_0 = pi * (angle / 180.0);
	tmp = ((a * cos(t_0)) ^ 2.0) + ((b * sin(t_0)) ^ 2.0);
end
code[a_, b_, angle_] := Block[{t$95$0 = N[(Pi * N[(angle / 180.0), $MachinePrecision]), $MachinePrecision]}, N[(N[Power[N[(a * N[Cos[t$95$0], $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision] + N[Power[N[(b * N[Sin[t$95$0], $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}

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

Sampling outcomes in binary64 precision:

Local Percentage Accuracy vs ?

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

Accuracy vs Speed?

Herbie found 12 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: 78.8% accurate, 1.0× speedup?

\[\begin{array}{l} \\ \begin{array}{l} t_0 := \pi \cdot \frac{angle}{180}\\ {\left(a \cdot \cos t\_0\right)}^{2} + {\left(b \cdot \sin t\_0\right)}^{2} \end{array} \end{array} \]
(FPCore (a b angle)
 :precision binary64
 (let* ((t_0 (* PI (/ angle 180.0))))
   (+ (pow (* a (cos t_0)) 2.0) (pow (* b (sin t_0)) 2.0))))
double code(double a, double b, double angle) {
	double t_0 = ((double) M_PI) * (angle / 180.0);
	return pow((a * cos(t_0)), 2.0) + pow((b * sin(t_0)), 2.0);
}
public static double code(double a, double b, double angle) {
	double t_0 = Math.PI * (angle / 180.0);
	return Math.pow((a * Math.cos(t_0)), 2.0) + Math.pow((b * Math.sin(t_0)), 2.0);
}
def code(a, b, angle):
	t_0 = math.pi * (angle / 180.0)
	return math.pow((a * math.cos(t_0)), 2.0) + math.pow((b * math.sin(t_0)), 2.0)
function code(a, b, angle)
	t_0 = Float64(pi * Float64(angle / 180.0))
	return Float64((Float64(a * cos(t_0)) ^ 2.0) + (Float64(b * sin(t_0)) ^ 2.0))
end
function tmp = code(a, b, angle)
	t_0 = pi * (angle / 180.0);
	tmp = ((a * cos(t_0)) ^ 2.0) + ((b * sin(t_0)) ^ 2.0);
end
code[a_, b_, angle_] := Block[{t$95$0 = N[(Pi * N[(angle / 180.0), $MachinePrecision]), $MachinePrecision]}, N[(N[Power[N[(a * N[Cos[t$95$0], $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision] + N[Power[N[(b * N[Sin[t$95$0], $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}

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

Alternative 1: 78.8% accurate, 0.6× speedup?

\[\begin{array}{l} angle_m = \left|angle\right| \\ {\left(a \cdot \cos \left(\pi \cdot {\left(e^{\log \left(angle\_m \cdot 0.005555555555555556\right) \cdot 0.3333333333333333}\right)}^{3}\right)\right)}^{2} + {\left(b \cdot \sin \left(\pi \cdot \left(angle\_m \cdot 0.005555555555555556\right)\right)\right)}^{2} \end{array} \]
angle_m = (fabs.f64 angle)
(FPCore (a b angle_m)
 :precision binary64
 (+
  (pow
   (*
    a
    (cos
     (*
      PI
      (pow
       (exp (* (log (* angle_m 0.005555555555555556)) 0.3333333333333333))
       3.0))))
   2.0)
  (pow (* b (sin (* PI (* angle_m 0.005555555555555556)))) 2.0)))
angle_m = fabs(angle);
double code(double a, double b, double angle_m) {
	return pow((a * cos((((double) M_PI) * pow(exp((log((angle_m * 0.005555555555555556)) * 0.3333333333333333)), 3.0)))), 2.0) + pow((b * sin((((double) M_PI) * (angle_m * 0.005555555555555556)))), 2.0);
}
angle_m = Math.abs(angle);
public static double code(double a, double b, double angle_m) {
	return Math.pow((a * Math.cos((Math.PI * Math.pow(Math.exp((Math.log((angle_m * 0.005555555555555556)) * 0.3333333333333333)), 3.0)))), 2.0) + Math.pow((b * Math.sin((Math.PI * (angle_m * 0.005555555555555556)))), 2.0);
}
angle_m = math.fabs(angle)
def code(a, b, angle_m):
	return math.pow((a * math.cos((math.pi * math.pow(math.exp((math.log((angle_m * 0.005555555555555556)) * 0.3333333333333333)), 3.0)))), 2.0) + math.pow((b * math.sin((math.pi * (angle_m * 0.005555555555555556)))), 2.0)
angle_m = abs(angle)
function code(a, b, angle_m)
	return Float64((Float64(a * cos(Float64(pi * (exp(Float64(log(Float64(angle_m * 0.005555555555555556)) * 0.3333333333333333)) ^ 3.0)))) ^ 2.0) + (Float64(b * sin(Float64(pi * Float64(angle_m * 0.005555555555555556)))) ^ 2.0))
end
angle_m = abs(angle);
function tmp = code(a, b, angle_m)
	tmp = ((a * cos((pi * (exp((log((angle_m * 0.005555555555555556)) * 0.3333333333333333)) ^ 3.0)))) ^ 2.0) + ((b * sin((pi * (angle_m * 0.005555555555555556)))) ^ 2.0);
end
angle_m = N[Abs[angle], $MachinePrecision]
code[a_, b_, angle$95$m_] := N[(N[Power[N[(a * N[Cos[N[(Pi * N[Power[N[Exp[N[(N[Log[N[(angle$95$m * 0.005555555555555556), $MachinePrecision]], $MachinePrecision] * 0.3333333333333333), $MachinePrecision]], $MachinePrecision], 3.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision] + N[Power[N[(b * N[Sin[N[(Pi * N[(angle$95$m * 0.005555555555555556), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
angle_m = \left|angle\right|

\\
{\left(a \cdot \cos \left(\pi \cdot {\left(e^{\log \left(angle\_m \cdot 0.005555555555555556\right) \cdot 0.3333333333333333}\right)}^{3}\right)\right)}^{2} + {\left(b \cdot \sin \left(\pi \cdot \left(angle\_m \cdot 0.005555555555555556\right)\right)\right)}^{2}
\end{array}
Derivation
  1. Initial program 83.5%

    \[{\left(a \cdot \cos \left(\pi \cdot \frac{angle}{180}\right)\right)}^{2} + {\left(b \cdot \sin \left(\pi \cdot \frac{angle}{180}\right)\right)}^{2} \]
  2. Step-by-step derivation
    1. Simplified83.5%

      \[\leadsto \color{blue}{{\left(a \cdot \cos \left(\pi \cdot \left(angle \cdot 0.005555555555555556\right)\right)\right)}^{2} + {\left(b \cdot \sin \left(\pi \cdot \left(angle \cdot 0.005555555555555556\right)\right)\right)}^{2}} \]
    2. Add Preprocessing
    3. Step-by-step derivation
      1. metadata-eval83.5%

        \[\leadsto {\left(a \cdot \cos \left(\pi \cdot \left(angle \cdot \color{blue}{\frac{1}{180}}\right)\right)\right)}^{2} + {\left(b \cdot \sin \left(\pi \cdot \left(angle \cdot 0.005555555555555556\right)\right)\right)}^{2} \]
      2. div-inv83.5%

        \[\leadsto {\left(a \cdot \cos \left(\pi \cdot \color{blue}{\frac{angle}{180}}\right)\right)}^{2} + {\left(b \cdot \sin \left(\pi \cdot \left(angle \cdot 0.005555555555555556\right)\right)\right)}^{2} \]
      3. add-cube-cbrt83.4%

        \[\leadsto {\left(a \cdot \cos \left(\pi \cdot \color{blue}{\left(\left(\sqrt[3]{\frac{angle}{180}} \cdot \sqrt[3]{\frac{angle}{180}}\right) \cdot \sqrt[3]{\frac{angle}{180}}\right)}\right)\right)}^{2} + {\left(b \cdot \sin \left(\pi \cdot \left(angle \cdot 0.005555555555555556\right)\right)\right)}^{2} \]
      4. pow383.5%

        \[\leadsto {\left(a \cdot \cos \left(\pi \cdot \color{blue}{{\left(\sqrt[3]{\frac{angle}{180}}\right)}^{3}}\right)\right)}^{2} + {\left(b \cdot \sin \left(\pi \cdot \left(angle \cdot 0.005555555555555556\right)\right)\right)}^{2} \]
      5. div-inv83.5%

        \[\leadsto {\left(a \cdot \cos \left(\pi \cdot {\left(\sqrt[3]{\color{blue}{angle \cdot \frac{1}{180}}}\right)}^{3}\right)\right)}^{2} + {\left(b \cdot \sin \left(\pi \cdot \left(angle \cdot 0.005555555555555556\right)\right)\right)}^{2} \]
      6. metadata-eval83.5%

        \[\leadsto {\left(a \cdot \cos \left(\pi \cdot {\left(\sqrt[3]{angle \cdot \color{blue}{0.005555555555555556}}\right)}^{3}\right)\right)}^{2} + {\left(b \cdot \sin \left(\pi \cdot \left(angle \cdot 0.005555555555555556\right)\right)\right)}^{2} \]
    4. Applied egg-rr83.5%

      \[\leadsto {\left(a \cdot \cos \left(\pi \cdot \color{blue}{{\left(\sqrt[3]{angle \cdot 0.005555555555555556}\right)}^{3}}\right)\right)}^{2} + {\left(b \cdot \sin \left(\pi \cdot \left(angle \cdot 0.005555555555555556\right)\right)\right)}^{2} \]
    5. Step-by-step derivation
      1. pow1/338.4%

        \[\leadsto {\left(a \cdot \cos \left(\pi \cdot {\color{blue}{\left({\left(angle \cdot 0.005555555555555556\right)}^{0.3333333333333333}\right)}}^{3}\right)\right)}^{2} + {\left(b \cdot \sin \left(\pi \cdot \left(angle \cdot 0.005555555555555556\right)\right)\right)}^{2} \]
      2. pow-to-exp38.5%

        \[\leadsto {\left(a \cdot \cos \left(\pi \cdot {\color{blue}{\left(e^{\log \left(angle \cdot 0.005555555555555556\right) \cdot 0.3333333333333333}\right)}}^{3}\right)\right)}^{2} + {\left(b \cdot \sin \left(\pi \cdot \left(angle \cdot 0.005555555555555556\right)\right)\right)}^{2} \]
    6. Applied egg-rr38.5%

      \[\leadsto {\left(a \cdot \cos \left(\pi \cdot {\color{blue}{\left(e^{\log \left(angle \cdot 0.005555555555555556\right) \cdot 0.3333333333333333}\right)}}^{3}\right)\right)}^{2} + {\left(b \cdot \sin \left(\pi \cdot \left(angle \cdot 0.005555555555555556\right)\right)\right)}^{2} \]
    7. Final simplification38.5%

      \[\leadsto {\left(a \cdot \cos \left(\pi \cdot {\left(e^{\log \left(angle \cdot 0.005555555555555556\right) \cdot 0.3333333333333333}\right)}^{3}\right)\right)}^{2} + {\left(b \cdot \sin \left(\pi \cdot \left(angle \cdot 0.005555555555555556\right)\right)\right)}^{2} \]
    8. Add Preprocessing

    Alternative 2: 78.9% accurate, 0.7× speedup?

    \[\begin{array}{l} angle_m = \left|angle\right| \\ {\left(a \cdot \cos \left(\pi \cdot \left(angle\_m \cdot 0.005555555555555556\right)\right)\right)}^{2} + {\left(b \cdot \mathsf{log1p}\left(\mathsf{expm1}\left(\sin \left(angle\_m \cdot \left(\pi \cdot 0.005555555555555556\right)\right)\right)\right)\right)}^{2} \end{array} \]
    angle_m = (fabs.f64 angle)
    (FPCore (a b angle_m)
     :precision binary64
     (+
      (pow (* a (cos (* PI (* angle_m 0.005555555555555556)))) 2.0)
      (pow
       (* b (log1p (expm1 (sin (* angle_m (* PI 0.005555555555555556))))))
       2.0)))
    angle_m = fabs(angle);
    double code(double a, double b, double angle_m) {
    	return pow((a * cos((((double) M_PI) * (angle_m * 0.005555555555555556)))), 2.0) + pow((b * log1p(expm1(sin((angle_m * (((double) M_PI) * 0.005555555555555556)))))), 2.0);
    }
    
    angle_m = Math.abs(angle);
    public static double code(double a, double b, double angle_m) {
    	return Math.pow((a * Math.cos((Math.PI * (angle_m * 0.005555555555555556)))), 2.0) + Math.pow((b * Math.log1p(Math.expm1(Math.sin((angle_m * (Math.PI * 0.005555555555555556)))))), 2.0);
    }
    
    angle_m = math.fabs(angle)
    def code(a, b, angle_m):
    	return math.pow((a * math.cos((math.pi * (angle_m * 0.005555555555555556)))), 2.0) + math.pow((b * math.log1p(math.expm1(math.sin((angle_m * (math.pi * 0.005555555555555556)))))), 2.0)
    
    angle_m = abs(angle)
    function code(a, b, angle_m)
    	return Float64((Float64(a * cos(Float64(pi * Float64(angle_m * 0.005555555555555556)))) ^ 2.0) + (Float64(b * log1p(expm1(sin(Float64(angle_m * Float64(pi * 0.005555555555555556)))))) ^ 2.0))
    end
    
    angle_m = N[Abs[angle], $MachinePrecision]
    code[a_, b_, angle$95$m_] := N[(N[Power[N[(a * N[Cos[N[(Pi * N[(angle$95$m * 0.005555555555555556), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision] + N[Power[N[(b * N[Log[1 + N[(Exp[N[Sin[N[(angle$95$m * N[(Pi * 0.005555555555555556), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]] - 1), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]
    
    \begin{array}{l}
    angle_m = \left|angle\right|
    
    \\
    {\left(a \cdot \cos \left(\pi \cdot \left(angle\_m \cdot 0.005555555555555556\right)\right)\right)}^{2} + {\left(b \cdot \mathsf{log1p}\left(\mathsf{expm1}\left(\sin \left(angle\_m \cdot \left(\pi \cdot 0.005555555555555556\right)\right)\right)\right)\right)}^{2}
    \end{array}
    
    Derivation
    1. Initial program 83.5%

      \[{\left(a \cdot \cos \left(\pi \cdot \frac{angle}{180}\right)\right)}^{2} + {\left(b \cdot \sin \left(\pi \cdot \frac{angle}{180}\right)\right)}^{2} \]
    2. Step-by-step derivation
      1. Simplified83.5%

        \[\leadsto \color{blue}{{\left(a \cdot \cos \left(\pi \cdot \left(angle \cdot 0.005555555555555556\right)\right)\right)}^{2} + {\left(b \cdot \sin \left(\pi \cdot \left(angle \cdot 0.005555555555555556\right)\right)\right)}^{2}} \]
      2. Add Preprocessing
      3. Step-by-step derivation
        1. log1p-expm1-u83.5%

          \[\leadsto {\left(a \cdot \cos \left(\pi \cdot \left(angle \cdot 0.005555555555555556\right)\right)\right)}^{2} + {\left(b \cdot \color{blue}{\mathsf{log1p}\left(\mathsf{expm1}\left(\sin \left(\pi \cdot \left(angle \cdot 0.005555555555555556\right)\right)\right)\right)}\right)}^{2} \]
        2. metadata-eval83.5%

          \[\leadsto {\left(a \cdot \cos \left(\pi \cdot \left(angle \cdot 0.005555555555555556\right)\right)\right)}^{2} + {\left(b \cdot \mathsf{log1p}\left(\mathsf{expm1}\left(\sin \left(\pi \cdot \left(angle \cdot \color{blue}{\frac{1}{180}}\right)\right)\right)\right)\right)}^{2} \]
        3. div-inv83.5%

          \[\leadsto {\left(a \cdot \cos \left(\pi \cdot \left(angle \cdot 0.005555555555555556\right)\right)\right)}^{2} + {\left(b \cdot \mathsf{log1p}\left(\mathsf{expm1}\left(\sin \left(\pi \cdot \color{blue}{\frac{angle}{180}}\right)\right)\right)\right)}^{2} \]
        4. *-commutative83.5%

          \[\leadsto {\left(a \cdot \cos \left(\pi \cdot \left(angle \cdot 0.005555555555555556\right)\right)\right)}^{2} + {\left(b \cdot \mathsf{log1p}\left(\mathsf{expm1}\left(\sin \color{blue}{\left(\frac{angle}{180} \cdot \pi\right)}\right)\right)\right)}^{2} \]
        5. div-inv83.5%

          \[\leadsto {\left(a \cdot \cos \left(\pi \cdot \left(angle \cdot 0.005555555555555556\right)\right)\right)}^{2} + {\left(b \cdot \mathsf{log1p}\left(\mathsf{expm1}\left(\sin \left(\color{blue}{\left(angle \cdot \frac{1}{180}\right)} \cdot \pi\right)\right)\right)\right)}^{2} \]
        6. metadata-eval83.5%

          \[\leadsto {\left(a \cdot \cos \left(\pi \cdot \left(angle \cdot 0.005555555555555556\right)\right)\right)}^{2} + {\left(b \cdot \mathsf{log1p}\left(\mathsf{expm1}\left(\sin \left(\left(angle \cdot \color{blue}{0.005555555555555556}\right) \cdot \pi\right)\right)\right)\right)}^{2} \]
        7. associate-*l*83.6%

          \[\leadsto {\left(a \cdot \cos \left(\pi \cdot \left(angle \cdot 0.005555555555555556\right)\right)\right)}^{2} + {\left(b \cdot \mathsf{log1p}\left(\mathsf{expm1}\left(\sin \color{blue}{\left(angle \cdot \left(0.005555555555555556 \cdot \pi\right)\right)}\right)\right)\right)}^{2} \]
      4. Applied egg-rr83.6%

        \[\leadsto {\left(a \cdot \cos \left(\pi \cdot \left(angle \cdot 0.005555555555555556\right)\right)\right)}^{2} + {\left(b \cdot \color{blue}{\mathsf{log1p}\left(\mathsf{expm1}\left(\sin \left(angle \cdot \left(0.005555555555555556 \cdot \pi\right)\right)\right)\right)}\right)}^{2} \]
      5. Final simplification83.6%

        \[\leadsto {\left(a \cdot \cos \left(\pi \cdot \left(angle \cdot 0.005555555555555556\right)\right)\right)}^{2} + {\left(b \cdot \mathsf{log1p}\left(\mathsf{expm1}\left(\sin \left(angle \cdot \left(\pi \cdot 0.005555555555555556\right)\right)\right)\right)\right)}^{2} \]
      6. Add Preprocessing

      Alternative 3: 78.8% accurate, 1.0× speedup?

      \[\begin{array}{l} angle_m = \left|angle\right| \\ \begin{array}{l} t_0 := \pi \cdot \left(angle\_m \cdot 0.005555555555555556\right)\\ {\left(b \cdot \sin t\_0\right)}^{2} + {\left(a \cdot \cos t\_0\right)}^{2} \end{array} \end{array} \]
      angle_m = (fabs.f64 angle)
      (FPCore (a b angle_m)
       :precision binary64
       (let* ((t_0 (* PI (* angle_m 0.005555555555555556))))
         (+ (pow (* b (sin t_0)) 2.0) (pow (* a (cos t_0)) 2.0))))
      angle_m = fabs(angle);
      double code(double a, double b, double angle_m) {
      	double t_0 = ((double) M_PI) * (angle_m * 0.005555555555555556);
      	return pow((b * sin(t_0)), 2.0) + pow((a * cos(t_0)), 2.0);
      }
      
      angle_m = Math.abs(angle);
      public static double code(double a, double b, double angle_m) {
      	double t_0 = Math.PI * (angle_m * 0.005555555555555556);
      	return Math.pow((b * Math.sin(t_0)), 2.0) + Math.pow((a * Math.cos(t_0)), 2.0);
      }
      
      angle_m = math.fabs(angle)
      def code(a, b, angle_m):
      	t_0 = math.pi * (angle_m * 0.005555555555555556)
      	return math.pow((b * math.sin(t_0)), 2.0) + math.pow((a * math.cos(t_0)), 2.0)
      
      angle_m = abs(angle)
      function code(a, b, angle_m)
      	t_0 = Float64(pi * Float64(angle_m * 0.005555555555555556))
      	return Float64((Float64(b * sin(t_0)) ^ 2.0) + (Float64(a * cos(t_0)) ^ 2.0))
      end
      
      angle_m = abs(angle);
      function tmp = code(a, b, angle_m)
      	t_0 = pi * (angle_m * 0.005555555555555556);
      	tmp = ((b * sin(t_0)) ^ 2.0) + ((a * cos(t_0)) ^ 2.0);
      end
      
      angle_m = N[Abs[angle], $MachinePrecision]
      code[a_, b_, angle$95$m_] := Block[{t$95$0 = N[(Pi * N[(angle$95$m * 0.005555555555555556), $MachinePrecision]), $MachinePrecision]}, N[(N[Power[N[(b * N[Sin[t$95$0], $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision] + N[Power[N[(a * N[Cos[t$95$0], $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]]
      
      \begin{array}{l}
      angle_m = \left|angle\right|
      
      \\
      \begin{array}{l}
      t_0 := \pi \cdot \left(angle\_m \cdot 0.005555555555555556\right)\\
      {\left(b \cdot \sin t\_0\right)}^{2} + {\left(a \cdot \cos t\_0\right)}^{2}
      \end{array}
      \end{array}
      
      Derivation
      1. Initial program 83.5%

        \[{\left(a \cdot \cos \left(\pi \cdot \frac{angle}{180}\right)\right)}^{2} + {\left(b \cdot \sin \left(\pi \cdot \frac{angle}{180}\right)\right)}^{2} \]
      2. Step-by-step derivation
        1. Simplified83.5%

          \[\leadsto \color{blue}{{\left(a \cdot \cos \left(\pi \cdot \left(angle \cdot 0.005555555555555556\right)\right)\right)}^{2} + {\left(b \cdot \sin \left(\pi \cdot \left(angle \cdot 0.005555555555555556\right)\right)\right)}^{2}} \]
        2. Add Preprocessing
        3. Final simplification83.5%

          \[\leadsto {\left(b \cdot \sin \left(\pi \cdot \left(angle \cdot 0.005555555555555556\right)\right)\right)}^{2} + {\left(a \cdot \cos \left(\pi \cdot \left(angle \cdot 0.005555555555555556\right)\right)\right)}^{2} \]
        4. Add Preprocessing

        Alternative 4: 78.8% accurate, 1.0× speedup?

        \[\begin{array}{l} angle_m = \left|angle\right| \\ {\left(a \cdot \cos \left(\pi \cdot \left(angle\_m \cdot 0.005555555555555556\right)\right)\right)}^{2} + {\left(b \cdot \sin \left(\frac{\pi}{\frac{180}{angle\_m}}\right)\right)}^{2} \end{array} \]
        angle_m = (fabs.f64 angle)
        (FPCore (a b angle_m)
         :precision binary64
         (+
          (pow (* a (cos (* PI (* angle_m 0.005555555555555556)))) 2.0)
          (pow (* b (sin (/ PI (/ 180.0 angle_m)))) 2.0)))
        angle_m = fabs(angle);
        double code(double a, double b, double angle_m) {
        	return pow((a * cos((((double) M_PI) * (angle_m * 0.005555555555555556)))), 2.0) + pow((b * sin((((double) M_PI) / (180.0 / angle_m)))), 2.0);
        }
        
        angle_m = Math.abs(angle);
        public static double code(double a, double b, double angle_m) {
        	return Math.pow((a * Math.cos((Math.PI * (angle_m * 0.005555555555555556)))), 2.0) + Math.pow((b * Math.sin((Math.PI / (180.0 / angle_m)))), 2.0);
        }
        
        angle_m = math.fabs(angle)
        def code(a, b, angle_m):
        	return math.pow((a * math.cos((math.pi * (angle_m * 0.005555555555555556)))), 2.0) + math.pow((b * math.sin((math.pi / (180.0 / angle_m)))), 2.0)
        
        angle_m = abs(angle)
        function code(a, b, angle_m)
        	return Float64((Float64(a * cos(Float64(pi * Float64(angle_m * 0.005555555555555556)))) ^ 2.0) + (Float64(b * sin(Float64(pi / Float64(180.0 / angle_m)))) ^ 2.0))
        end
        
        angle_m = abs(angle);
        function tmp = code(a, b, angle_m)
        	tmp = ((a * cos((pi * (angle_m * 0.005555555555555556)))) ^ 2.0) + ((b * sin((pi / (180.0 / angle_m)))) ^ 2.0);
        end
        
        angle_m = N[Abs[angle], $MachinePrecision]
        code[a_, b_, angle$95$m_] := N[(N[Power[N[(a * N[Cos[N[(Pi * N[(angle$95$m * 0.005555555555555556), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision] + N[Power[N[(b * N[Sin[N[(Pi / N[(180.0 / angle$95$m), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]
        
        \begin{array}{l}
        angle_m = \left|angle\right|
        
        \\
        {\left(a \cdot \cos \left(\pi \cdot \left(angle\_m \cdot 0.005555555555555556\right)\right)\right)}^{2} + {\left(b \cdot \sin \left(\frac{\pi}{\frac{180}{angle\_m}}\right)\right)}^{2}
        \end{array}
        
        Derivation
        1. Initial program 83.5%

          \[{\left(a \cdot \cos \left(\pi \cdot \frac{angle}{180}\right)\right)}^{2} + {\left(b \cdot \sin \left(\pi \cdot \frac{angle}{180}\right)\right)}^{2} \]
        2. Step-by-step derivation
          1. Simplified83.5%

            \[\leadsto \color{blue}{{\left(a \cdot \cos \left(\pi \cdot \left(angle \cdot 0.005555555555555556\right)\right)\right)}^{2} + {\left(b \cdot \sin \left(\pi \cdot \left(angle \cdot 0.005555555555555556\right)\right)\right)}^{2}} \]
          2. Add Preprocessing
          3. Step-by-step derivation
            1. metadata-eval82.8%

              \[\leadsto {\left(a \cdot 1\right)}^{2} + {\left(b \cdot \sin \left(\pi \cdot \left(angle \cdot \color{blue}{\frac{1}{180}}\right)\right)\right)}^{2} \]
            2. div-inv82.8%

              \[\leadsto {\left(a \cdot 1\right)}^{2} + {\left(b \cdot \sin \left(\pi \cdot \color{blue}{\frac{angle}{180}}\right)\right)}^{2} \]
            3. clear-num82.8%

              \[\leadsto {\left(a \cdot 1\right)}^{2} + {\left(b \cdot \sin \left(\pi \cdot \color{blue}{\frac{1}{\frac{180}{angle}}}\right)\right)}^{2} \]
            4. un-div-inv82.8%

              \[\leadsto {\left(a \cdot 1\right)}^{2} + {\left(b \cdot \sin \color{blue}{\left(\frac{\pi}{\frac{180}{angle}}\right)}\right)}^{2} \]
          4. Applied egg-rr83.5%

            \[\leadsto {\left(a \cdot \cos \left(\pi \cdot \left(angle \cdot 0.005555555555555556\right)\right)\right)}^{2} + {\left(b \cdot \sin \color{blue}{\left(\frac{\pi}{\frac{180}{angle}}\right)}\right)}^{2} \]
          5. Final simplification83.5%

            \[\leadsto {\left(a \cdot \cos \left(\pi \cdot \left(angle \cdot 0.005555555555555556\right)\right)\right)}^{2} + {\left(b \cdot \sin \left(\frac{\pi}{\frac{180}{angle}}\right)\right)}^{2} \]
          6. Add Preprocessing

          Alternative 5: 78.8% accurate, 1.0× speedup?

          \[\begin{array}{l} angle_m = \left|angle\right| \\ \begin{array}{l} t_0 := \pi \cdot \frac{angle\_m}{180}\\ {\left(a \cdot \cos t\_0\right)}^{2} + {\left(b \cdot \sin t\_0\right)}^{2} \end{array} \end{array} \]
          angle_m = (fabs.f64 angle)
          (FPCore (a b angle_m)
           :precision binary64
           (let* ((t_0 (* PI (/ angle_m 180.0))))
             (+ (pow (* a (cos t_0)) 2.0) (pow (* b (sin t_0)) 2.0))))
          angle_m = fabs(angle);
          double code(double a, double b, double angle_m) {
          	double t_0 = ((double) M_PI) * (angle_m / 180.0);
          	return pow((a * cos(t_0)), 2.0) + pow((b * sin(t_0)), 2.0);
          }
          
          angle_m = Math.abs(angle);
          public static double code(double a, double b, double angle_m) {
          	double t_0 = Math.PI * (angle_m / 180.0);
          	return Math.pow((a * Math.cos(t_0)), 2.0) + Math.pow((b * Math.sin(t_0)), 2.0);
          }
          
          angle_m = math.fabs(angle)
          def code(a, b, angle_m):
          	t_0 = math.pi * (angle_m / 180.0)
          	return math.pow((a * math.cos(t_0)), 2.0) + math.pow((b * math.sin(t_0)), 2.0)
          
          angle_m = abs(angle)
          function code(a, b, angle_m)
          	t_0 = Float64(pi * Float64(angle_m / 180.0))
          	return Float64((Float64(a * cos(t_0)) ^ 2.0) + (Float64(b * sin(t_0)) ^ 2.0))
          end
          
          angle_m = abs(angle);
          function tmp = code(a, b, angle_m)
          	t_0 = pi * (angle_m / 180.0);
          	tmp = ((a * cos(t_0)) ^ 2.0) + ((b * sin(t_0)) ^ 2.0);
          end
          
          angle_m = N[Abs[angle], $MachinePrecision]
          code[a_, b_, angle$95$m_] := Block[{t$95$0 = N[(Pi * N[(angle$95$m / 180.0), $MachinePrecision]), $MachinePrecision]}, N[(N[Power[N[(a * N[Cos[t$95$0], $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision] + N[Power[N[(b * N[Sin[t$95$0], $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]]
          
          \begin{array}{l}
          angle_m = \left|angle\right|
          
          \\
          \begin{array}{l}
          t_0 := \pi \cdot \frac{angle\_m}{180}\\
          {\left(a \cdot \cos t\_0\right)}^{2} + {\left(b \cdot \sin t\_0\right)}^{2}
          \end{array}
          \end{array}
          
          Derivation
          1. Initial program 83.5%

            \[{\left(a \cdot \cos \left(\pi \cdot \frac{angle}{180}\right)\right)}^{2} + {\left(b \cdot \sin \left(\pi \cdot \frac{angle}{180}\right)\right)}^{2} \]
          2. Add Preprocessing
          3. Final simplification83.5%

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

          Alternative 6: 78.8% accurate, 1.3× speedup?

          \[\begin{array}{l} angle_m = \left|angle\right| \\ {a}^{2} + {\left(b \cdot \sin \left(\pi \cdot \frac{1}{\frac{180}{angle\_m}}\right)\right)}^{2} \end{array} \]
          angle_m = (fabs.f64 angle)
          (FPCore (a b angle_m)
           :precision binary64
           (+ (pow a 2.0) (pow (* b (sin (* PI (/ 1.0 (/ 180.0 angle_m))))) 2.0)))
          angle_m = fabs(angle);
          double code(double a, double b, double angle_m) {
          	return pow(a, 2.0) + pow((b * sin((((double) M_PI) * (1.0 / (180.0 / angle_m))))), 2.0);
          }
          
          angle_m = Math.abs(angle);
          public static double code(double a, double b, double angle_m) {
          	return Math.pow(a, 2.0) + Math.pow((b * Math.sin((Math.PI * (1.0 / (180.0 / angle_m))))), 2.0);
          }
          
          angle_m = math.fabs(angle)
          def code(a, b, angle_m):
          	return math.pow(a, 2.0) + math.pow((b * math.sin((math.pi * (1.0 / (180.0 / angle_m))))), 2.0)
          
          angle_m = abs(angle)
          function code(a, b, angle_m)
          	return Float64((a ^ 2.0) + (Float64(b * sin(Float64(pi * Float64(1.0 / Float64(180.0 / angle_m))))) ^ 2.0))
          end
          
          angle_m = abs(angle);
          function tmp = code(a, b, angle_m)
          	tmp = (a ^ 2.0) + ((b * sin((pi * (1.0 / (180.0 / angle_m))))) ^ 2.0);
          end
          
          angle_m = N[Abs[angle], $MachinePrecision]
          code[a_, b_, angle$95$m_] := N[(N[Power[a, 2.0], $MachinePrecision] + N[Power[N[(b * N[Sin[N[(Pi * N[(1.0 / N[(180.0 / angle$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]
          
          \begin{array}{l}
          angle_m = \left|angle\right|
          
          \\
          {a}^{2} + {\left(b \cdot \sin \left(\pi \cdot \frac{1}{\frac{180}{angle\_m}}\right)\right)}^{2}
          \end{array}
          
          Derivation
          1. Initial program 83.5%

            \[{\left(a \cdot \cos \left(\pi \cdot \frac{angle}{180}\right)\right)}^{2} + {\left(b \cdot \sin \left(\pi \cdot \frac{angle}{180}\right)\right)}^{2} \]
          2. Step-by-step derivation
            1. Simplified83.5%

              \[\leadsto \color{blue}{{\left(a \cdot \cos \left(\pi \cdot \left(angle \cdot 0.005555555555555556\right)\right)\right)}^{2} + {\left(b \cdot \sin \left(\pi \cdot \left(angle \cdot 0.005555555555555556\right)\right)\right)}^{2}} \]
            2. Add Preprocessing
            3. Taylor expanded in angle around 0 82.8%

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

                \[\leadsto {\left(a \cdot 1\right)}^{2} + {\left(b \cdot \sin \left(\pi \cdot \left(angle \cdot \color{blue}{\frac{1}{180}}\right)\right)\right)}^{2} \]
              2. div-inv82.8%

                \[\leadsto {\left(a \cdot 1\right)}^{2} + {\left(b \cdot \sin \left(\pi \cdot \color{blue}{\frac{angle}{180}}\right)\right)}^{2} \]
              3. clear-num82.8%

                \[\leadsto {\left(a \cdot 1\right)}^{2} + {\left(b \cdot \sin \left(\pi \cdot \color{blue}{\frac{1}{\frac{180}{angle}}}\right)\right)}^{2} \]
            5. Applied egg-rr82.8%

              \[\leadsto {\left(a \cdot 1\right)}^{2} + {\left(b \cdot \sin \left(\pi \cdot \color{blue}{\frac{1}{\frac{180}{angle}}}\right)\right)}^{2} \]
            6. Final simplification82.8%

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

            Alternative 7: 78.8% accurate, 1.3× speedup?

            \[\begin{array}{l} angle_m = \left|angle\right| \\ {\left(b \cdot \sin \left(\pi \cdot \left(angle\_m \cdot 0.005555555555555556\right)\right)\right)}^{2} + {a}^{2} \end{array} \]
            angle_m = (fabs.f64 angle)
            (FPCore (a b angle_m)
             :precision binary64
             (+ (pow (* b (sin (* PI (* angle_m 0.005555555555555556)))) 2.0) (pow a 2.0)))
            angle_m = fabs(angle);
            double code(double a, double b, double angle_m) {
            	return pow((b * sin((((double) M_PI) * (angle_m * 0.005555555555555556)))), 2.0) + pow(a, 2.0);
            }
            
            angle_m = Math.abs(angle);
            public static double code(double a, double b, double angle_m) {
            	return Math.pow((b * Math.sin((Math.PI * (angle_m * 0.005555555555555556)))), 2.0) + Math.pow(a, 2.0);
            }
            
            angle_m = math.fabs(angle)
            def code(a, b, angle_m):
            	return math.pow((b * math.sin((math.pi * (angle_m * 0.005555555555555556)))), 2.0) + math.pow(a, 2.0)
            
            angle_m = abs(angle)
            function code(a, b, angle_m)
            	return Float64((Float64(b * sin(Float64(pi * Float64(angle_m * 0.005555555555555556)))) ^ 2.0) + (a ^ 2.0))
            end
            
            angle_m = abs(angle);
            function tmp = code(a, b, angle_m)
            	tmp = ((b * sin((pi * (angle_m * 0.005555555555555556)))) ^ 2.0) + (a ^ 2.0);
            end
            
            angle_m = N[Abs[angle], $MachinePrecision]
            code[a_, b_, angle$95$m_] := N[(N[Power[N[(b * N[Sin[N[(Pi * N[(angle$95$m * 0.005555555555555556), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision] + N[Power[a, 2.0], $MachinePrecision]), $MachinePrecision]
            
            \begin{array}{l}
            angle_m = \left|angle\right|
            
            \\
            {\left(b \cdot \sin \left(\pi \cdot \left(angle\_m \cdot 0.005555555555555556\right)\right)\right)}^{2} + {a}^{2}
            \end{array}
            
            Derivation
            1. Initial program 83.5%

              \[{\left(a \cdot \cos \left(\pi \cdot \frac{angle}{180}\right)\right)}^{2} + {\left(b \cdot \sin \left(\pi \cdot \frac{angle}{180}\right)\right)}^{2} \]
            2. Step-by-step derivation
              1. Simplified83.5%

                \[\leadsto \color{blue}{{\left(a \cdot \cos \left(\pi \cdot \left(angle \cdot 0.005555555555555556\right)\right)\right)}^{2} + {\left(b \cdot \sin \left(\pi \cdot \left(angle \cdot 0.005555555555555556\right)\right)\right)}^{2}} \]
              2. Add Preprocessing
              3. Taylor expanded in angle around 0 82.8%

                \[\leadsto {\left(a \cdot \color{blue}{1}\right)}^{2} + {\left(b \cdot \sin \left(\pi \cdot \left(angle \cdot 0.005555555555555556\right)\right)\right)}^{2} \]
              4. Final simplification82.8%

                \[\leadsto {\left(b \cdot \sin \left(\pi \cdot \left(angle \cdot 0.005555555555555556\right)\right)\right)}^{2} + {a}^{2} \]
              5. Add Preprocessing

              Alternative 8: 78.8% accurate, 1.3× speedup?

              \[\begin{array}{l} angle_m = \left|angle\right| \\ {\left(b \cdot \sin \left(\frac{\pi}{\frac{180}{angle\_m}}\right)\right)}^{2} + {a}^{2} \end{array} \]
              angle_m = (fabs.f64 angle)
              (FPCore (a b angle_m)
               :precision binary64
               (+ (pow (* b (sin (/ PI (/ 180.0 angle_m)))) 2.0) (pow a 2.0)))
              angle_m = fabs(angle);
              double code(double a, double b, double angle_m) {
              	return pow((b * sin((((double) M_PI) / (180.0 / angle_m)))), 2.0) + pow(a, 2.0);
              }
              
              angle_m = Math.abs(angle);
              public static double code(double a, double b, double angle_m) {
              	return Math.pow((b * Math.sin((Math.PI / (180.0 / angle_m)))), 2.0) + Math.pow(a, 2.0);
              }
              
              angle_m = math.fabs(angle)
              def code(a, b, angle_m):
              	return math.pow((b * math.sin((math.pi / (180.0 / angle_m)))), 2.0) + math.pow(a, 2.0)
              
              angle_m = abs(angle)
              function code(a, b, angle_m)
              	return Float64((Float64(b * sin(Float64(pi / Float64(180.0 / angle_m)))) ^ 2.0) + (a ^ 2.0))
              end
              
              angle_m = abs(angle);
              function tmp = code(a, b, angle_m)
              	tmp = ((b * sin((pi / (180.0 / angle_m)))) ^ 2.0) + (a ^ 2.0);
              end
              
              angle_m = N[Abs[angle], $MachinePrecision]
              code[a_, b_, angle$95$m_] := N[(N[Power[N[(b * N[Sin[N[(Pi / N[(180.0 / angle$95$m), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision] + N[Power[a, 2.0], $MachinePrecision]), $MachinePrecision]
              
              \begin{array}{l}
              angle_m = \left|angle\right|
              
              \\
              {\left(b \cdot \sin \left(\frac{\pi}{\frac{180}{angle\_m}}\right)\right)}^{2} + {a}^{2}
              \end{array}
              
              Derivation
              1. Initial program 83.5%

                \[{\left(a \cdot \cos \left(\pi \cdot \frac{angle}{180}\right)\right)}^{2} + {\left(b \cdot \sin \left(\pi \cdot \frac{angle}{180}\right)\right)}^{2} \]
              2. Step-by-step derivation
                1. Simplified83.5%

                  \[\leadsto \color{blue}{{\left(a \cdot \cos \left(\pi \cdot \left(angle \cdot 0.005555555555555556\right)\right)\right)}^{2} + {\left(b \cdot \sin \left(\pi \cdot \left(angle \cdot 0.005555555555555556\right)\right)\right)}^{2}} \]
                2. Add Preprocessing
                3. Taylor expanded in angle around 0 82.8%

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

                    \[\leadsto {\left(a \cdot 1\right)}^{2} + {\left(b \cdot \sin \left(\pi \cdot \left(angle \cdot \color{blue}{\frac{1}{180}}\right)\right)\right)}^{2} \]
                  2. div-inv82.8%

                    \[\leadsto {\left(a \cdot 1\right)}^{2} + {\left(b \cdot \sin \left(\pi \cdot \color{blue}{\frac{angle}{180}}\right)\right)}^{2} \]
                  3. clear-num82.8%

                    \[\leadsto {\left(a \cdot 1\right)}^{2} + {\left(b \cdot \sin \left(\pi \cdot \color{blue}{\frac{1}{\frac{180}{angle}}}\right)\right)}^{2} \]
                  4. un-div-inv82.8%

                    \[\leadsto {\left(a \cdot 1\right)}^{2} + {\left(b \cdot \sin \color{blue}{\left(\frac{\pi}{\frac{180}{angle}}\right)}\right)}^{2} \]
                5. Applied egg-rr82.8%

                  \[\leadsto {\left(a \cdot 1\right)}^{2} + {\left(b \cdot \sin \color{blue}{\left(\frac{\pi}{\frac{180}{angle}}\right)}\right)}^{2} \]
                6. Final simplification82.8%

                  \[\leadsto {\left(b \cdot \sin \left(\frac{\pi}{\frac{180}{angle}}\right)\right)}^{2} + {a}^{2} \]
                7. Add Preprocessing

                Alternative 9: 72.2% accurate, 3.5× speedup?

                \[\begin{array}{l} angle_m = \left|angle\right| \\ {a}^{2} + \left(angle\_m \cdot 0.005555555555555556\right) \cdot \left(\left(\pi \cdot b\right) \cdot \left(0.005555555555555556 \cdot \left(angle\_m \cdot \left(\pi \cdot b\right)\right)\right)\right) \end{array} \]
                angle_m = (fabs.f64 angle)
                (FPCore (a b angle_m)
                 :precision binary64
                 (+
                  (pow a 2.0)
                  (*
                   (* angle_m 0.005555555555555556)
                   (* (* PI b) (* 0.005555555555555556 (* angle_m (* PI b)))))))
                angle_m = fabs(angle);
                double code(double a, double b, double angle_m) {
                	return pow(a, 2.0) + ((angle_m * 0.005555555555555556) * ((((double) M_PI) * b) * (0.005555555555555556 * (angle_m * (((double) M_PI) * b)))));
                }
                
                angle_m = Math.abs(angle);
                public static double code(double a, double b, double angle_m) {
                	return Math.pow(a, 2.0) + ((angle_m * 0.005555555555555556) * ((Math.PI * b) * (0.005555555555555556 * (angle_m * (Math.PI * b)))));
                }
                
                angle_m = math.fabs(angle)
                def code(a, b, angle_m):
                	return math.pow(a, 2.0) + ((angle_m * 0.005555555555555556) * ((math.pi * b) * (0.005555555555555556 * (angle_m * (math.pi * b)))))
                
                angle_m = abs(angle)
                function code(a, b, angle_m)
                	return Float64((a ^ 2.0) + Float64(Float64(angle_m * 0.005555555555555556) * Float64(Float64(pi * b) * Float64(0.005555555555555556 * Float64(angle_m * Float64(pi * b))))))
                end
                
                angle_m = abs(angle);
                function tmp = code(a, b, angle_m)
                	tmp = (a ^ 2.0) + ((angle_m * 0.005555555555555556) * ((pi * b) * (0.005555555555555556 * (angle_m * (pi * b)))));
                end
                
                angle_m = N[Abs[angle], $MachinePrecision]
                code[a_, b_, angle$95$m_] := N[(N[Power[a, 2.0], $MachinePrecision] + N[(N[(angle$95$m * 0.005555555555555556), $MachinePrecision] * N[(N[(Pi * b), $MachinePrecision] * N[(0.005555555555555556 * N[(angle$95$m * N[(Pi * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
                
                \begin{array}{l}
                angle_m = \left|angle\right|
                
                \\
                {a}^{2} + \left(angle\_m \cdot 0.005555555555555556\right) \cdot \left(\left(\pi \cdot b\right) \cdot \left(0.005555555555555556 \cdot \left(angle\_m \cdot \left(\pi \cdot b\right)\right)\right)\right)
                \end{array}
                
                Derivation
                1. Initial program 83.5%

                  \[{\left(a \cdot \cos \left(\pi \cdot \frac{angle}{180}\right)\right)}^{2} + {\left(b \cdot \sin \left(\pi \cdot \frac{angle}{180}\right)\right)}^{2} \]
                2. Step-by-step derivation
                  1. Simplified83.5%

                    \[\leadsto \color{blue}{{\left(a \cdot \cos \left(\pi \cdot \left(angle \cdot 0.005555555555555556\right)\right)\right)}^{2} + {\left(b \cdot \sin \left(\pi \cdot \left(angle \cdot 0.005555555555555556\right)\right)\right)}^{2}} \]
                  2. Add Preprocessing
                  3. Taylor expanded in angle around 0 82.8%

                    \[\leadsto {\left(a \cdot \color{blue}{1}\right)}^{2} + {\left(b \cdot \sin \left(\pi \cdot \left(angle \cdot 0.005555555555555556\right)\right)\right)}^{2} \]
                  4. Taylor expanded in angle around 0 78.5%

                    \[\leadsto {\left(a \cdot 1\right)}^{2} + {\color{blue}{\left(0.005555555555555556 \cdot \left(angle \cdot \left(b \cdot \pi\right)\right)\right)}}^{2} \]
                  5. Step-by-step derivation
                    1. unpow278.5%

                      \[\leadsto {\left(a \cdot 1\right)}^{2} + \color{blue}{\left(0.005555555555555556 \cdot \left(angle \cdot \left(b \cdot \pi\right)\right)\right) \cdot \left(0.005555555555555556 \cdot \left(angle \cdot \left(b \cdot \pi\right)\right)\right)} \]
                    2. associate-*r*78.5%

                      \[\leadsto {\left(a \cdot 1\right)}^{2} + \color{blue}{\left(\left(0.005555555555555556 \cdot angle\right) \cdot \left(b \cdot \pi\right)\right)} \cdot \left(0.005555555555555556 \cdot \left(angle \cdot \left(b \cdot \pi\right)\right)\right) \]
                    3. *-commutative78.5%

                      \[\leadsto {\left(a \cdot 1\right)}^{2} + \left(\color{blue}{\left(angle \cdot 0.005555555555555556\right)} \cdot \left(b \cdot \pi\right)\right) \cdot \left(0.005555555555555556 \cdot \left(angle \cdot \left(b \cdot \pi\right)\right)\right) \]
                    4. associate-*l*77.5%

                      \[\leadsto {\left(a \cdot 1\right)}^{2} + \color{blue}{\left(angle \cdot 0.005555555555555556\right) \cdot \left(\left(b \cdot \pi\right) \cdot \left(0.005555555555555556 \cdot \left(angle \cdot \left(b \cdot \pi\right)\right)\right)\right)} \]
                    5. *-commutative77.5%

                      \[\leadsto {\left(a \cdot 1\right)}^{2} + \left(angle \cdot 0.005555555555555556\right) \cdot \left(\color{blue}{\left(\pi \cdot b\right)} \cdot \left(0.005555555555555556 \cdot \left(angle \cdot \left(b \cdot \pi\right)\right)\right)\right) \]
                    6. associate-*r*77.4%

                      \[\leadsto {\left(a \cdot 1\right)}^{2} + \left(angle \cdot 0.005555555555555556\right) \cdot \left(\left(\pi \cdot b\right) \cdot \color{blue}{\left(\left(0.005555555555555556 \cdot angle\right) \cdot \left(b \cdot \pi\right)\right)}\right) \]
                    7. *-commutative77.4%

                      \[\leadsto {\left(a \cdot 1\right)}^{2} + \left(angle \cdot 0.005555555555555556\right) \cdot \left(\left(\pi \cdot b\right) \cdot \left(\color{blue}{\left(angle \cdot 0.005555555555555556\right)} \cdot \left(b \cdot \pi\right)\right)\right) \]
                    8. *-commutative77.4%

                      \[\leadsto {\left(a \cdot 1\right)}^{2} + \left(angle \cdot 0.005555555555555556\right) \cdot \left(\left(\pi \cdot b\right) \cdot \left(\left(angle \cdot 0.005555555555555556\right) \cdot \color{blue}{\left(\pi \cdot b\right)}\right)\right) \]
                  6. Applied egg-rr77.4%

                    \[\leadsto {\left(a \cdot 1\right)}^{2} + \color{blue}{\left(angle \cdot 0.005555555555555556\right) \cdot \left(\left(\pi \cdot b\right) \cdot \left(\left(angle \cdot 0.005555555555555556\right) \cdot \left(\pi \cdot b\right)\right)\right)} \]
                  7. Taylor expanded in angle around 0 77.5%

                    \[\leadsto {\left(a \cdot 1\right)}^{2} + \left(angle \cdot 0.005555555555555556\right) \cdot \left(\left(\pi \cdot b\right) \cdot \color{blue}{\left(0.005555555555555556 \cdot \left(angle \cdot \left(b \cdot \pi\right)\right)\right)}\right) \]
                  8. Step-by-step derivation
                    1. *-commutative77.5%

                      \[\leadsto {\left(a \cdot 1\right)}^{2} + \left(angle \cdot 0.005555555555555556\right) \cdot \left(\left(\pi \cdot b\right) \cdot \left(0.005555555555555556 \cdot \left(angle \cdot \color{blue}{\left(\pi \cdot b\right)}\right)\right)\right) \]
                  9. Simplified77.5%

                    \[\leadsto {\left(a \cdot 1\right)}^{2} + \left(angle \cdot 0.005555555555555556\right) \cdot \left(\left(\pi \cdot b\right) \cdot \color{blue}{\left(0.005555555555555556 \cdot \left(angle \cdot \left(\pi \cdot b\right)\right)\right)}\right) \]
                  10. Final simplification77.5%

                    \[\leadsto {a}^{2} + \left(angle \cdot 0.005555555555555556\right) \cdot \left(\left(\pi \cdot b\right) \cdot \left(0.005555555555555556 \cdot \left(angle \cdot \left(\pi \cdot b\right)\right)\right)\right) \]
                  11. Add Preprocessing

                  Alternative 10: 72.2% accurate, 3.5× speedup?

                  \[\begin{array}{l} angle_m = \left|angle\right| \\ {a}^{2} + \left(angle\_m \cdot 0.005555555555555556\right) \cdot \left(\left(\pi \cdot b\right) \cdot \left(b \cdot \left(angle\_m \cdot \left(\pi \cdot 0.005555555555555556\right)\right)\right)\right) \end{array} \]
                  angle_m = (fabs.f64 angle)
                  (FPCore (a b angle_m)
                   :precision binary64
                   (+
                    (pow a 2.0)
                    (*
                     (* angle_m 0.005555555555555556)
                     (* (* PI b) (* b (* angle_m (* PI 0.005555555555555556)))))))
                  angle_m = fabs(angle);
                  double code(double a, double b, double angle_m) {
                  	return pow(a, 2.0) + ((angle_m * 0.005555555555555556) * ((((double) M_PI) * b) * (b * (angle_m * (((double) M_PI) * 0.005555555555555556)))));
                  }
                  
                  angle_m = Math.abs(angle);
                  public static double code(double a, double b, double angle_m) {
                  	return Math.pow(a, 2.0) + ((angle_m * 0.005555555555555556) * ((Math.PI * b) * (b * (angle_m * (Math.PI * 0.005555555555555556)))));
                  }
                  
                  angle_m = math.fabs(angle)
                  def code(a, b, angle_m):
                  	return math.pow(a, 2.0) + ((angle_m * 0.005555555555555556) * ((math.pi * b) * (b * (angle_m * (math.pi * 0.005555555555555556)))))
                  
                  angle_m = abs(angle)
                  function code(a, b, angle_m)
                  	return Float64((a ^ 2.0) + Float64(Float64(angle_m * 0.005555555555555556) * Float64(Float64(pi * b) * Float64(b * Float64(angle_m * Float64(pi * 0.005555555555555556))))))
                  end
                  
                  angle_m = abs(angle);
                  function tmp = code(a, b, angle_m)
                  	tmp = (a ^ 2.0) + ((angle_m * 0.005555555555555556) * ((pi * b) * (b * (angle_m * (pi * 0.005555555555555556)))));
                  end
                  
                  angle_m = N[Abs[angle], $MachinePrecision]
                  code[a_, b_, angle$95$m_] := N[(N[Power[a, 2.0], $MachinePrecision] + N[(N[(angle$95$m * 0.005555555555555556), $MachinePrecision] * N[(N[(Pi * b), $MachinePrecision] * N[(b * N[(angle$95$m * N[(Pi * 0.005555555555555556), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
                  
                  \begin{array}{l}
                  angle_m = \left|angle\right|
                  
                  \\
                  {a}^{2} + \left(angle\_m \cdot 0.005555555555555556\right) \cdot \left(\left(\pi \cdot b\right) \cdot \left(b \cdot \left(angle\_m \cdot \left(\pi \cdot 0.005555555555555556\right)\right)\right)\right)
                  \end{array}
                  
                  Derivation
                  1. Initial program 83.5%

                    \[{\left(a \cdot \cos \left(\pi \cdot \frac{angle}{180}\right)\right)}^{2} + {\left(b \cdot \sin \left(\pi \cdot \frac{angle}{180}\right)\right)}^{2} \]
                  2. Step-by-step derivation
                    1. Simplified83.5%

                      \[\leadsto \color{blue}{{\left(a \cdot \cos \left(\pi \cdot \left(angle \cdot 0.005555555555555556\right)\right)\right)}^{2} + {\left(b \cdot \sin \left(\pi \cdot \left(angle \cdot 0.005555555555555556\right)\right)\right)}^{2}} \]
                    2. Add Preprocessing
                    3. Taylor expanded in angle around 0 82.8%

                      \[\leadsto {\left(a \cdot \color{blue}{1}\right)}^{2} + {\left(b \cdot \sin \left(\pi \cdot \left(angle \cdot 0.005555555555555556\right)\right)\right)}^{2} \]
                    4. Taylor expanded in angle around 0 78.5%

                      \[\leadsto {\left(a \cdot 1\right)}^{2} + {\color{blue}{\left(0.005555555555555556 \cdot \left(angle \cdot \left(b \cdot \pi\right)\right)\right)}}^{2} \]
                    5. Step-by-step derivation
                      1. unpow278.5%

                        \[\leadsto {\left(a \cdot 1\right)}^{2} + \color{blue}{\left(0.005555555555555556 \cdot \left(angle \cdot \left(b \cdot \pi\right)\right)\right) \cdot \left(0.005555555555555556 \cdot \left(angle \cdot \left(b \cdot \pi\right)\right)\right)} \]
                      2. associate-*r*78.5%

                        \[\leadsto {\left(a \cdot 1\right)}^{2} + \color{blue}{\left(\left(0.005555555555555556 \cdot angle\right) \cdot \left(b \cdot \pi\right)\right)} \cdot \left(0.005555555555555556 \cdot \left(angle \cdot \left(b \cdot \pi\right)\right)\right) \]
                      3. *-commutative78.5%

                        \[\leadsto {\left(a \cdot 1\right)}^{2} + \left(\color{blue}{\left(angle \cdot 0.005555555555555556\right)} \cdot \left(b \cdot \pi\right)\right) \cdot \left(0.005555555555555556 \cdot \left(angle \cdot \left(b \cdot \pi\right)\right)\right) \]
                      4. associate-*l*77.5%

                        \[\leadsto {\left(a \cdot 1\right)}^{2} + \color{blue}{\left(angle \cdot 0.005555555555555556\right) \cdot \left(\left(b \cdot \pi\right) \cdot \left(0.005555555555555556 \cdot \left(angle \cdot \left(b \cdot \pi\right)\right)\right)\right)} \]
                      5. *-commutative77.5%

                        \[\leadsto {\left(a \cdot 1\right)}^{2} + \left(angle \cdot 0.005555555555555556\right) \cdot \left(\color{blue}{\left(\pi \cdot b\right)} \cdot \left(0.005555555555555556 \cdot \left(angle \cdot \left(b \cdot \pi\right)\right)\right)\right) \]
                      6. associate-*r*77.4%

                        \[\leadsto {\left(a \cdot 1\right)}^{2} + \left(angle \cdot 0.005555555555555556\right) \cdot \left(\left(\pi \cdot b\right) \cdot \color{blue}{\left(\left(0.005555555555555556 \cdot angle\right) \cdot \left(b \cdot \pi\right)\right)}\right) \]
                      7. *-commutative77.4%

                        \[\leadsto {\left(a \cdot 1\right)}^{2} + \left(angle \cdot 0.005555555555555556\right) \cdot \left(\left(\pi \cdot b\right) \cdot \left(\color{blue}{\left(angle \cdot 0.005555555555555556\right)} \cdot \left(b \cdot \pi\right)\right)\right) \]
                      8. *-commutative77.4%

                        \[\leadsto {\left(a \cdot 1\right)}^{2} + \left(angle \cdot 0.005555555555555556\right) \cdot \left(\left(\pi \cdot b\right) \cdot \left(\left(angle \cdot 0.005555555555555556\right) \cdot \color{blue}{\left(\pi \cdot b\right)}\right)\right) \]
                    6. Applied egg-rr77.4%

                      \[\leadsto {\left(a \cdot 1\right)}^{2} + \color{blue}{\left(angle \cdot 0.005555555555555556\right) \cdot \left(\left(\pi \cdot b\right) \cdot \left(\left(angle \cdot 0.005555555555555556\right) \cdot \left(\pi \cdot b\right)\right)\right)} \]
                    7. Taylor expanded in angle around 0 77.5%

                      \[\leadsto {\left(a \cdot 1\right)}^{2} + \left(angle \cdot 0.005555555555555556\right) \cdot \left(\left(\pi \cdot b\right) \cdot \color{blue}{\left(0.005555555555555556 \cdot \left(angle \cdot \left(b \cdot \pi\right)\right)\right)}\right) \]
                    8. Step-by-step derivation
                      1. associate-*r*77.4%

                        \[\leadsto {\left(a \cdot 1\right)}^{2} + \left(angle \cdot 0.005555555555555556\right) \cdot \left(\left(\pi \cdot b\right) \cdot \color{blue}{\left(\left(0.005555555555555556 \cdot angle\right) \cdot \left(b \cdot \pi\right)\right)}\right) \]
                      2. *-commutative77.4%

                        \[\leadsto {\left(a \cdot 1\right)}^{2} + \left(angle \cdot 0.005555555555555556\right) \cdot \left(\left(\pi \cdot b\right) \cdot \left(\color{blue}{\left(angle \cdot 0.005555555555555556\right)} \cdot \left(b \cdot \pi\right)\right)\right) \]
                      3. associate-*r*77.4%

                        \[\leadsto {\left(a \cdot 1\right)}^{2} + \left(angle \cdot 0.005555555555555556\right) \cdot \left(\left(\pi \cdot b\right) \cdot \color{blue}{\left(\left(\left(angle \cdot 0.005555555555555556\right) \cdot b\right) \cdot \pi\right)}\right) \]
                      4. *-commutative77.4%

                        \[\leadsto {\left(a \cdot 1\right)}^{2} + \left(angle \cdot 0.005555555555555556\right) \cdot \left(\left(\pi \cdot b\right) \cdot \left(\color{blue}{\left(b \cdot \left(angle \cdot 0.005555555555555556\right)\right)} \cdot \pi\right)\right) \]
                      5. associate-*l*77.4%

                        \[\leadsto {\left(a \cdot 1\right)}^{2} + \left(angle \cdot 0.005555555555555556\right) \cdot \left(\left(\pi \cdot b\right) \cdot \color{blue}{\left(b \cdot \left(\left(angle \cdot 0.005555555555555556\right) \cdot \pi\right)\right)}\right) \]
                      6. associate-*r*77.5%

                        \[\leadsto {\left(a \cdot 1\right)}^{2} + \left(angle \cdot 0.005555555555555556\right) \cdot \left(\left(\pi \cdot b\right) \cdot \left(b \cdot \color{blue}{\left(angle \cdot \left(0.005555555555555556 \cdot \pi\right)\right)}\right)\right) \]
                      7. *-commutative77.5%

                        \[\leadsto {\left(a \cdot 1\right)}^{2} + \left(angle \cdot 0.005555555555555556\right) \cdot \left(\left(\pi \cdot b\right) \cdot \left(b \cdot \left(angle \cdot \color{blue}{\left(\pi \cdot 0.005555555555555556\right)}\right)\right)\right) \]
                    9. Simplified77.5%

                      \[\leadsto {\left(a \cdot 1\right)}^{2} + \left(angle \cdot 0.005555555555555556\right) \cdot \left(\left(\pi \cdot b\right) \cdot \color{blue}{\left(b \cdot \left(angle \cdot \left(\pi \cdot 0.005555555555555556\right)\right)\right)}\right) \]
                    10. Final simplification77.5%

                      \[\leadsto {a}^{2} + \left(angle \cdot 0.005555555555555556\right) \cdot \left(\left(\pi \cdot b\right) \cdot \left(b \cdot \left(angle \cdot \left(\pi \cdot 0.005555555555555556\right)\right)\right)\right) \]
                    11. Add Preprocessing

                    Alternative 11: 73.6% accurate, 3.5× speedup?

                    \[\begin{array}{l} angle_m = \left|angle\right| \\ \begin{array}{l} t_0 := \left(angle\_m \cdot 0.005555555555555556\right) \cdot \left(\pi \cdot b\right)\\ {a}^{2} + t\_0 \cdot t\_0 \end{array} \end{array} \]
                    angle_m = (fabs.f64 angle)
                    (FPCore (a b angle_m)
                     :precision binary64
                     (let* ((t_0 (* (* angle_m 0.005555555555555556) (* PI b))))
                       (+ (pow a 2.0) (* t_0 t_0))))
                    angle_m = fabs(angle);
                    double code(double a, double b, double angle_m) {
                    	double t_0 = (angle_m * 0.005555555555555556) * (((double) M_PI) * b);
                    	return pow(a, 2.0) + (t_0 * t_0);
                    }
                    
                    angle_m = Math.abs(angle);
                    public static double code(double a, double b, double angle_m) {
                    	double t_0 = (angle_m * 0.005555555555555556) * (Math.PI * b);
                    	return Math.pow(a, 2.0) + (t_0 * t_0);
                    }
                    
                    angle_m = math.fabs(angle)
                    def code(a, b, angle_m):
                    	t_0 = (angle_m * 0.005555555555555556) * (math.pi * b)
                    	return math.pow(a, 2.0) + (t_0 * t_0)
                    
                    angle_m = abs(angle)
                    function code(a, b, angle_m)
                    	t_0 = Float64(Float64(angle_m * 0.005555555555555556) * Float64(pi * b))
                    	return Float64((a ^ 2.0) + Float64(t_0 * t_0))
                    end
                    
                    angle_m = abs(angle);
                    function tmp = code(a, b, angle_m)
                    	t_0 = (angle_m * 0.005555555555555556) * (pi * b);
                    	tmp = (a ^ 2.0) + (t_0 * t_0);
                    end
                    
                    angle_m = N[Abs[angle], $MachinePrecision]
                    code[a_, b_, angle$95$m_] := Block[{t$95$0 = N[(N[(angle$95$m * 0.005555555555555556), $MachinePrecision] * N[(Pi * b), $MachinePrecision]), $MachinePrecision]}, N[(N[Power[a, 2.0], $MachinePrecision] + N[(t$95$0 * t$95$0), $MachinePrecision]), $MachinePrecision]]
                    
                    \begin{array}{l}
                    angle_m = \left|angle\right|
                    
                    \\
                    \begin{array}{l}
                    t_0 := \left(angle\_m \cdot 0.005555555555555556\right) \cdot \left(\pi \cdot b\right)\\
                    {a}^{2} + t\_0 \cdot t\_0
                    \end{array}
                    \end{array}
                    
                    Derivation
                    1. Initial program 83.5%

                      \[{\left(a \cdot \cos \left(\pi \cdot \frac{angle}{180}\right)\right)}^{2} + {\left(b \cdot \sin \left(\pi \cdot \frac{angle}{180}\right)\right)}^{2} \]
                    2. Step-by-step derivation
                      1. Simplified83.5%

                        \[\leadsto \color{blue}{{\left(a \cdot \cos \left(\pi \cdot \left(angle \cdot 0.005555555555555556\right)\right)\right)}^{2} + {\left(b \cdot \sin \left(\pi \cdot \left(angle \cdot 0.005555555555555556\right)\right)\right)}^{2}} \]
                      2. Add Preprocessing
                      3. Taylor expanded in angle around 0 82.8%

                        \[\leadsto {\left(a \cdot \color{blue}{1}\right)}^{2} + {\left(b \cdot \sin \left(\pi \cdot \left(angle \cdot 0.005555555555555556\right)\right)\right)}^{2} \]
                      4. Taylor expanded in angle around 0 78.5%

                        \[\leadsto {\left(a \cdot 1\right)}^{2} + {\color{blue}{\left(0.005555555555555556 \cdot \left(angle \cdot \left(b \cdot \pi\right)\right)\right)}}^{2} \]
                      5. Step-by-step derivation
                        1. unpow278.5%

                          \[\leadsto {\left(a \cdot 1\right)}^{2} + \color{blue}{\left(0.005555555555555556 \cdot \left(angle \cdot \left(b \cdot \pi\right)\right)\right) \cdot \left(0.005555555555555556 \cdot \left(angle \cdot \left(b \cdot \pi\right)\right)\right)} \]
                        2. associate-*r*78.5%

                          \[\leadsto {\left(a \cdot 1\right)}^{2} + \color{blue}{\left(\left(0.005555555555555556 \cdot angle\right) \cdot \left(b \cdot \pi\right)\right)} \cdot \left(0.005555555555555556 \cdot \left(angle \cdot \left(b \cdot \pi\right)\right)\right) \]
                        3. *-commutative78.5%

                          \[\leadsto {\left(a \cdot 1\right)}^{2} + \left(\color{blue}{\left(angle \cdot 0.005555555555555556\right)} \cdot \left(b \cdot \pi\right)\right) \cdot \left(0.005555555555555556 \cdot \left(angle \cdot \left(b \cdot \pi\right)\right)\right) \]
                        4. *-commutative78.5%

                          \[\leadsto {\left(a \cdot 1\right)}^{2} + \left(\left(angle \cdot 0.005555555555555556\right) \cdot \color{blue}{\left(\pi \cdot b\right)}\right) \cdot \left(0.005555555555555556 \cdot \left(angle \cdot \left(b \cdot \pi\right)\right)\right) \]
                        5. associate-*r*78.5%

                          \[\leadsto {\left(a \cdot 1\right)}^{2} + \left(\left(angle \cdot 0.005555555555555556\right) \cdot \left(\pi \cdot b\right)\right) \cdot \color{blue}{\left(\left(0.005555555555555556 \cdot angle\right) \cdot \left(b \cdot \pi\right)\right)} \]
                        6. *-commutative78.5%

                          \[\leadsto {\left(a \cdot 1\right)}^{2} + \left(\left(angle \cdot 0.005555555555555556\right) \cdot \left(\pi \cdot b\right)\right) \cdot \left(\color{blue}{\left(angle \cdot 0.005555555555555556\right)} \cdot \left(b \cdot \pi\right)\right) \]
                        7. *-commutative78.5%

                          \[\leadsto {\left(a \cdot 1\right)}^{2} + \left(\left(angle \cdot 0.005555555555555556\right) \cdot \left(\pi \cdot b\right)\right) \cdot \left(\left(angle \cdot 0.005555555555555556\right) \cdot \color{blue}{\left(\pi \cdot b\right)}\right) \]
                      6. Applied egg-rr78.5%

                        \[\leadsto {\left(a \cdot 1\right)}^{2} + \color{blue}{\left(\left(angle \cdot 0.005555555555555556\right) \cdot \left(\pi \cdot b\right)\right) \cdot \left(\left(angle \cdot 0.005555555555555556\right) \cdot \left(\pi \cdot b\right)\right)} \]
                      7. Final simplification78.5%

                        \[\leadsto {a}^{2} + \left(\left(angle \cdot 0.005555555555555556\right) \cdot \left(\pi \cdot b\right)\right) \cdot \left(\left(angle \cdot 0.005555555555555556\right) \cdot \left(\pi \cdot b\right)\right) \]
                      8. Add Preprocessing

                      Alternative 12: 73.7% accurate, 3.5× speedup?

                      \[\begin{array}{l} angle_m = \left|angle\right| \\ {a}^{2} + \left(0.005555555555555556 \cdot \left(\left(angle\_m \cdot 0.005555555555555556\right) \cdot \left(\pi \cdot b\right)\right)\right) \cdot \left(\pi \cdot \left(angle\_m \cdot b\right)\right) \end{array} \]
                      angle_m = (fabs.f64 angle)
                      (FPCore (a b angle_m)
                       :precision binary64
                       (+
                        (pow a 2.0)
                        (*
                         (* 0.005555555555555556 (* (* angle_m 0.005555555555555556) (* PI b)))
                         (* PI (* angle_m b)))))
                      angle_m = fabs(angle);
                      double code(double a, double b, double angle_m) {
                      	return pow(a, 2.0) + ((0.005555555555555556 * ((angle_m * 0.005555555555555556) * (((double) M_PI) * b))) * (((double) M_PI) * (angle_m * b)));
                      }
                      
                      angle_m = Math.abs(angle);
                      public static double code(double a, double b, double angle_m) {
                      	return Math.pow(a, 2.0) + ((0.005555555555555556 * ((angle_m * 0.005555555555555556) * (Math.PI * b))) * (Math.PI * (angle_m * b)));
                      }
                      
                      angle_m = math.fabs(angle)
                      def code(a, b, angle_m):
                      	return math.pow(a, 2.0) + ((0.005555555555555556 * ((angle_m * 0.005555555555555556) * (math.pi * b))) * (math.pi * (angle_m * b)))
                      
                      angle_m = abs(angle)
                      function code(a, b, angle_m)
                      	return Float64((a ^ 2.0) + Float64(Float64(0.005555555555555556 * Float64(Float64(angle_m * 0.005555555555555556) * Float64(pi * b))) * Float64(pi * Float64(angle_m * b))))
                      end
                      
                      angle_m = abs(angle);
                      function tmp = code(a, b, angle_m)
                      	tmp = (a ^ 2.0) + ((0.005555555555555556 * ((angle_m * 0.005555555555555556) * (pi * b))) * (pi * (angle_m * b)));
                      end
                      
                      angle_m = N[Abs[angle], $MachinePrecision]
                      code[a_, b_, angle$95$m_] := N[(N[Power[a, 2.0], $MachinePrecision] + N[(N[(0.005555555555555556 * N[(N[(angle$95$m * 0.005555555555555556), $MachinePrecision] * N[(Pi * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(Pi * N[(angle$95$m * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
                      
                      \begin{array}{l}
                      angle_m = \left|angle\right|
                      
                      \\
                      {a}^{2} + \left(0.005555555555555556 \cdot \left(\left(angle\_m \cdot 0.005555555555555556\right) \cdot \left(\pi \cdot b\right)\right)\right) \cdot \left(\pi \cdot \left(angle\_m \cdot b\right)\right)
                      \end{array}
                      
                      Derivation
                      1. Initial program 83.5%

                        \[{\left(a \cdot \cos \left(\pi \cdot \frac{angle}{180}\right)\right)}^{2} + {\left(b \cdot \sin \left(\pi \cdot \frac{angle}{180}\right)\right)}^{2} \]
                      2. Step-by-step derivation
                        1. Simplified83.5%

                          \[\leadsto \color{blue}{{\left(a \cdot \cos \left(\pi \cdot \left(angle \cdot 0.005555555555555556\right)\right)\right)}^{2} + {\left(b \cdot \sin \left(\pi \cdot \left(angle \cdot 0.005555555555555556\right)\right)\right)}^{2}} \]
                        2. Add Preprocessing
                        3. Taylor expanded in angle around 0 82.8%

                          \[\leadsto {\left(a \cdot \color{blue}{1}\right)}^{2} + {\left(b \cdot \sin \left(\pi \cdot \left(angle \cdot 0.005555555555555556\right)\right)\right)}^{2} \]
                        4. Taylor expanded in angle around 0 78.5%

                          \[\leadsto {\left(a \cdot 1\right)}^{2} + {\color{blue}{\left(0.005555555555555556 \cdot \left(angle \cdot \left(b \cdot \pi\right)\right)\right)}}^{2} \]
                        5. Step-by-step derivation
                          1. unpow278.5%

                            \[\leadsto {\left(a \cdot 1\right)}^{2} + \color{blue}{\left(0.005555555555555556 \cdot \left(angle \cdot \left(b \cdot \pi\right)\right)\right) \cdot \left(0.005555555555555556 \cdot \left(angle \cdot \left(b \cdot \pi\right)\right)\right)} \]
                          2. associate-*r*78.5%

                            \[\leadsto {\left(a \cdot 1\right)}^{2} + \color{blue}{\left(\left(0.005555555555555556 \cdot \left(angle \cdot \left(b \cdot \pi\right)\right)\right) \cdot 0.005555555555555556\right) \cdot \left(angle \cdot \left(b \cdot \pi\right)\right)} \]
                          3. associate-*r*78.5%

                            \[\leadsto {\left(a \cdot 1\right)}^{2} + \left(\color{blue}{\left(\left(0.005555555555555556 \cdot angle\right) \cdot \left(b \cdot \pi\right)\right)} \cdot 0.005555555555555556\right) \cdot \left(angle \cdot \left(b \cdot \pi\right)\right) \]
                          4. *-commutative78.5%

                            \[\leadsto {\left(a \cdot 1\right)}^{2} + \left(\left(\color{blue}{\left(angle \cdot 0.005555555555555556\right)} \cdot \left(b \cdot \pi\right)\right) \cdot 0.005555555555555556\right) \cdot \left(angle \cdot \left(b \cdot \pi\right)\right) \]
                          5. *-commutative78.5%

                            \[\leadsto {\left(a \cdot 1\right)}^{2} + \left(\left(\left(angle \cdot 0.005555555555555556\right) \cdot \color{blue}{\left(\pi \cdot b\right)}\right) \cdot 0.005555555555555556\right) \cdot \left(angle \cdot \left(b \cdot \pi\right)\right) \]
                          6. associate-*r*78.5%

                            \[\leadsto {\left(a \cdot 1\right)}^{2} + \left(\left(\left(angle \cdot 0.005555555555555556\right) \cdot \left(\pi \cdot b\right)\right) \cdot 0.005555555555555556\right) \cdot \color{blue}{\left(\left(angle \cdot b\right) \cdot \pi\right)} \]
                          7. *-commutative78.5%

                            \[\leadsto {\left(a \cdot 1\right)}^{2} + \left(\left(\left(angle \cdot 0.005555555555555556\right) \cdot \left(\pi \cdot b\right)\right) \cdot 0.005555555555555556\right) \cdot \color{blue}{\left(\pi \cdot \left(angle \cdot b\right)\right)} \]
                        6. Applied egg-rr78.5%

                          \[\leadsto {\left(a \cdot 1\right)}^{2} + \color{blue}{\left(\left(\left(angle \cdot 0.005555555555555556\right) \cdot \left(\pi \cdot b\right)\right) \cdot 0.005555555555555556\right) \cdot \left(\pi \cdot \left(angle \cdot b\right)\right)} \]
                        7. Final simplification78.5%

                          \[\leadsto {a}^{2} + \left(0.005555555555555556 \cdot \left(\left(angle \cdot 0.005555555555555556\right) \cdot \left(\pi \cdot b\right)\right)\right) \cdot \left(\pi \cdot \left(angle \cdot b\right)\right) \]
                        8. Add Preprocessing

                        Reproduce

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