
(FPCore (a b angle x-scale y-scale)
:precision binary64
(let* ((t_0 (* (/ angle 180.0) (PI)))
(t_1 (sin t_0))
(t_2 (cos t_0))
(t_3
(/ (/ (+ (pow (* a t_2) 2.0) (pow (* b t_1) 2.0)) y-scale) y-scale))
(t_4
(/ (/ (+ (pow (* a t_1) 2.0) (pow (* b t_2) 2.0)) x-scale) x-scale))
(t_5 (* (* b a) (* b (- a))))
(t_6 (/ (* 4.0 t_5) (pow (* x-scale y-scale) 2.0))))
(/
(-
(sqrt
(*
(* (* 2.0 t_6) t_5)
(+
(+ t_4 t_3)
(sqrt
(+
(pow (- t_4 t_3) 2.0)
(pow
(/
(/ (* (* (* 2.0 (- (pow b 2.0) (pow a 2.0))) t_1) t_2) x-scale)
y-scale)
2.0)))))))
t_6)))\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{angle}{180} \cdot \mathsf{PI}\left(\right)\\
t_1 := \sin t\_0\\
t_2 := \cos t\_0\\
t_3 := \frac{\frac{{\left(a \cdot t\_2\right)}^{2} + {\left(b \cdot t\_1\right)}^{2}}{y-scale}}{y-scale}\\
t_4 := \frac{\frac{{\left(a \cdot t\_1\right)}^{2} + {\left(b \cdot t\_2\right)}^{2}}{x-scale}}{x-scale}\\
t_5 := \left(b \cdot a\right) \cdot \left(b \cdot \left(-a\right)\right)\\
t_6 := \frac{4 \cdot t\_5}{{\left(x-scale \cdot y-scale\right)}^{2}}\\
\frac{-\sqrt{\left(\left(2 \cdot t\_6\right) \cdot t\_5\right) \cdot \left(\left(t\_4 + t\_3\right) + \sqrt{{\left(t\_4 - t\_3\right)}^{2} + {\left(\frac{\frac{\left(\left(2 \cdot \left({b}^{2} - {a}^{2}\right)\right) \cdot t\_1\right) \cdot t\_2}{x-scale}}{y-scale}\right)}^{2}}\right)}}{t\_6}
\end{array}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 7 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (a b angle x-scale y-scale)
:precision binary64
(let* ((t_0 (* (/ angle 180.0) (PI)))
(t_1 (sin t_0))
(t_2 (cos t_0))
(t_3
(/ (/ (+ (pow (* a t_2) 2.0) (pow (* b t_1) 2.0)) y-scale) y-scale))
(t_4
(/ (/ (+ (pow (* a t_1) 2.0) (pow (* b t_2) 2.0)) x-scale) x-scale))
(t_5 (* (* b a) (* b (- a))))
(t_6 (/ (* 4.0 t_5) (pow (* x-scale y-scale) 2.0))))
(/
(-
(sqrt
(*
(* (* 2.0 t_6) t_5)
(+
(+ t_4 t_3)
(sqrt
(+
(pow (- t_4 t_3) 2.0)
(pow
(/
(/ (* (* (* 2.0 (- (pow b 2.0) (pow a 2.0))) t_1) t_2) x-scale)
y-scale)
2.0)))))))
t_6)))\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{angle}{180} \cdot \mathsf{PI}\left(\right)\\
t_1 := \sin t\_0\\
t_2 := \cos t\_0\\
t_3 := \frac{\frac{{\left(a \cdot t\_2\right)}^{2} + {\left(b \cdot t\_1\right)}^{2}}{y-scale}}{y-scale}\\
t_4 := \frac{\frac{{\left(a \cdot t\_1\right)}^{2} + {\left(b \cdot t\_2\right)}^{2}}{x-scale}}{x-scale}\\
t_5 := \left(b \cdot a\right) \cdot \left(b \cdot \left(-a\right)\right)\\
t_6 := \frac{4 \cdot t\_5}{{\left(x-scale \cdot y-scale\right)}^{2}}\\
\frac{-\sqrt{\left(\left(2 \cdot t\_6\right) \cdot t\_5\right) \cdot \left(\left(t\_4 + t\_3\right) + \sqrt{{\left(t\_4 - t\_3\right)}^{2} + {\left(\frac{\frac{\left(\left(2 \cdot \left({b}^{2} - {a}^{2}\right)\right) \cdot t\_1\right) \cdot t\_2}{x-scale}}{y-scale}\right)}^{2}}\right)}}{t\_6}
\end{array}
\end{array}
y-scale_m = (fabs.f64 y-scale)
x-scale_m = (fabs.f64 x-scale)
(FPCore (a b angle x-scale_m y-scale_m)
:precision binary64
(let* ((t_0 (* (* (PI) angle) 0.005555555555555556))
(t_1 (* (PI) (* angle 0.005555555555555556))))
(if (<= y-scale_m 1.85e+84)
(*
(*
(* 0.25 (* (sqrt 8.0) x-scale_m))
(hypot (* b (sin t_1)) (* a (cos t_1))))
(sqrt 2.0))
(*
(hypot (* (cos t_0) b) (* (sin t_0) a))
(* (* (* (sqrt 2.0) y-scale_m) (sqrt 8.0)) 0.25)))))\begin{array}{l}
y-scale_m = \left|y-scale\right|
\\
x-scale_m = \left|x-scale\right|
\\
\begin{array}{l}
t_0 := \left(\mathsf{PI}\left(\right) \cdot angle\right) \cdot 0.005555555555555556\\
t_1 := \mathsf{PI}\left(\right) \cdot \left(angle \cdot 0.005555555555555556\right)\\
\mathbf{if}\;y-scale\_m \leq 1.85 \cdot 10^{+84}:\\
\;\;\;\;\left(\left(0.25 \cdot \left(\sqrt{8} \cdot x-scale\_m\right)\right) \cdot \mathsf{hypot}\left(b \cdot \sin t\_1, a \cdot \cos t\_1\right)\right) \cdot \sqrt{2}\\
\mathbf{else}:\\
\;\;\;\;\mathsf{hypot}\left(\cos t\_0 \cdot b, \sin t\_0 \cdot a\right) \cdot \left(\left(\left(\sqrt{2} \cdot y-scale\_m\right) \cdot \sqrt{8}\right) \cdot 0.25\right)\\
\end{array}
\end{array}
if y-scale < 1.85e84Initial program 2.2%
Taylor expanded in y-scale around 0
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-sqrt.f64N/A
lower-sqrt.f64N/A
distribute-lft-outN/A
lower-*.f64N/A
Applied rewrites23.5%
Applied rewrites26.0%
if 1.85e84 < y-scale Initial program 6.9%
Taylor expanded in y-scale around inf
Applied rewrites25.9%
Taylor expanded in b around -inf
Applied rewrites30.8%
Taylor expanded in x-scale around 0
Applied rewrites72.5%
Final simplification36.7%
y-scale_m = (fabs.f64 y-scale)
x-scale_m = (fabs.f64 x-scale)
(FPCore (a b angle x-scale_m y-scale_m)
:precision binary64
(let* ((t_0 (* (* (PI) angle) 0.005555555555555556))
(t_1 (cos t_0))
(t_2 (sin t_0)))
(if (<= y-scale_m 1.85e+84)
(*
(hypot (* t_2 b) (* t_1 a))
(* (* (sqrt 8.0) (sqrt 2.0)) (* 0.25 x-scale_m)))
(*
(hypot (* t_1 b) (* t_2 a))
(* (* (* (sqrt 2.0) y-scale_m) (sqrt 8.0)) 0.25)))))\begin{array}{l}
y-scale_m = \left|y-scale\right|
\\
x-scale_m = \left|x-scale\right|
\\
\begin{array}{l}
t_0 := \left(\mathsf{PI}\left(\right) \cdot angle\right) \cdot 0.005555555555555556\\
t_1 := \cos t\_0\\
t_2 := \sin t\_0\\
\mathbf{if}\;y-scale\_m \leq 1.85 \cdot 10^{+84}:\\
\;\;\;\;\mathsf{hypot}\left(t\_2 \cdot b, t\_1 \cdot a\right) \cdot \left(\left(\sqrt{8} \cdot \sqrt{2}\right) \cdot \left(0.25 \cdot x-scale\_m\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{hypot}\left(t\_1 \cdot b, t\_2 \cdot a\right) \cdot \left(\left(\left(\sqrt{2} \cdot y-scale\_m\right) \cdot \sqrt{8}\right) \cdot 0.25\right)\\
\end{array}
\end{array}
if y-scale < 1.85e84Initial program 2.2%
Taylor expanded in y-scale around 0
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-sqrt.f64N/A
lower-sqrt.f64N/A
distribute-lft-outN/A
lower-*.f64N/A
Applied rewrites23.5%
Applied rewrites26.0%
Taylor expanded in angle around inf
Applied rewrites25.9%
if 1.85e84 < y-scale Initial program 6.9%
Taylor expanded in y-scale around inf
Applied rewrites25.9%
Taylor expanded in b around -inf
Applied rewrites30.8%
Taylor expanded in x-scale around 0
Applied rewrites72.5%
Final simplification36.6%
y-scale_m = (fabs.f64 y-scale)
x-scale_m = (fabs.f64 x-scale)
(FPCore (a b angle x-scale_m y-scale_m)
:precision binary64
(let* ((t_0 (* (* (PI) angle) 0.005555555555555556)))
(if (<= y-scale_m 1.85e+84)
(*
(hypot (* b (sin (* (PI) (* angle 0.005555555555555556)))) (* 1.0 a))
(* (* 0.25 (* (sqrt 8.0) x-scale_m)) (sqrt 2.0)))
(*
(hypot (* (cos t_0) b) (* (sin t_0) a))
(* (* (* (sqrt 2.0) y-scale_m) (sqrt 8.0)) 0.25)))))\begin{array}{l}
y-scale_m = \left|y-scale\right|
\\
x-scale_m = \left|x-scale\right|
\\
\begin{array}{l}
t_0 := \left(\mathsf{PI}\left(\right) \cdot angle\right) \cdot 0.005555555555555556\\
\mathbf{if}\;y-scale\_m \leq 1.85 \cdot 10^{+84}:\\
\;\;\;\;\mathsf{hypot}\left(b \cdot \sin \left(\mathsf{PI}\left(\right) \cdot \left(angle \cdot 0.005555555555555556\right)\right), 1 \cdot a\right) \cdot \left(\left(0.25 \cdot \left(\sqrt{8} \cdot x-scale\_m\right)\right) \cdot \sqrt{2}\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{hypot}\left(\cos t\_0 \cdot b, \sin t\_0 \cdot a\right) \cdot \left(\left(\left(\sqrt{2} \cdot y-scale\_m\right) \cdot \sqrt{8}\right) \cdot 0.25\right)\\
\end{array}
\end{array}
if y-scale < 1.85e84Initial program 2.2%
Taylor expanded in y-scale around 0
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-sqrt.f64N/A
lower-sqrt.f64N/A
distribute-lft-outN/A
lower-*.f64N/A
Applied rewrites23.5%
Applied rewrites26.0%
Taylor expanded in angle around 0
Applied rewrites26.0%
if 1.85e84 < y-scale Initial program 6.9%
Taylor expanded in y-scale around inf
Applied rewrites25.9%
Taylor expanded in b around -inf
Applied rewrites30.8%
Taylor expanded in x-scale around 0
Applied rewrites72.5%
Final simplification36.7%
y-scale_m = (fabs.f64 y-scale)
x-scale_m = (fabs.f64 x-scale)
(FPCore (a b angle x-scale_m y-scale_m)
:precision binary64
(if (<= y-scale_m 3.1e+84)
(*
(hypot (* b (sin (* (PI) (* angle 0.005555555555555556)))) (* 1.0 a))
(* (* 0.25 (* (sqrt 8.0) x-scale_m)) (sqrt 2.0)))
(* b y-scale_m)))\begin{array}{l}
y-scale_m = \left|y-scale\right|
\\
x-scale_m = \left|x-scale\right|
\\
\begin{array}{l}
\mathbf{if}\;y-scale\_m \leq 3.1 \cdot 10^{+84}:\\
\;\;\;\;\mathsf{hypot}\left(b \cdot \sin \left(\mathsf{PI}\left(\right) \cdot \left(angle \cdot 0.005555555555555556\right)\right), 1 \cdot a\right) \cdot \left(\left(0.25 \cdot \left(\sqrt{8} \cdot x-scale\_m\right)\right) \cdot \sqrt{2}\right)\\
\mathbf{else}:\\
\;\;\;\;b \cdot y-scale\_m\\
\end{array}
\end{array}
if y-scale < 3.10000000000000003e84Initial program 2.2%
Taylor expanded in y-scale around 0
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-sqrt.f64N/A
lower-sqrt.f64N/A
distribute-lft-outN/A
lower-*.f64N/A
Applied rewrites23.5%
Applied rewrites26.0%
Taylor expanded in angle around 0
Applied rewrites26.0%
if 3.10000000000000003e84 < y-scale Initial program 6.9%
Taylor expanded in angle around 0
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
lower-sqrt.f64N/A
lower-sqrt.f6428.9
Applied rewrites28.9%
Applied rewrites29.0%
Taylor expanded in b around 0
Applied rewrites29.0%
Final simplification26.7%
y-scale_m = (fabs.f64 y-scale)
x-scale_m = (fabs.f64 x-scale)
(FPCore (a b angle x-scale_m y-scale_m)
:precision binary64
(if (<= y-scale_m 2.4e+85)
(*
(hypot (* (* (* (PI) angle) 0.005555555555555556) b) (* 1.0 a))
(* (* 0.25 (* (sqrt 8.0) x-scale_m)) (sqrt 2.0)))
(* b y-scale_m)))\begin{array}{l}
y-scale_m = \left|y-scale\right|
\\
x-scale_m = \left|x-scale\right|
\\
\begin{array}{l}
\mathbf{if}\;y-scale\_m \leq 2.4 \cdot 10^{+85}:\\
\;\;\;\;\mathsf{hypot}\left(\left(\left(\mathsf{PI}\left(\right) \cdot angle\right) \cdot 0.005555555555555556\right) \cdot b, 1 \cdot a\right) \cdot \left(\left(0.25 \cdot \left(\sqrt{8} \cdot x-scale\_m\right)\right) \cdot \sqrt{2}\right)\\
\mathbf{else}:\\
\;\;\;\;b \cdot y-scale\_m\\
\end{array}
\end{array}
if y-scale < 2.39999999999999997e85Initial program 2.2%
Taylor expanded in y-scale around 0
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-sqrt.f64N/A
lower-sqrt.f64N/A
distribute-lft-outN/A
lower-*.f64N/A
Applied rewrites23.5%
Applied rewrites26.0%
Taylor expanded in angle around 0
Applied rewrites25.4%
Taylor expanded in angle around 0
Applied rewrites25.4%
if 2.39999999999999997e85 < y-scale Initial program 6.9%
Taylor expanded in angle around 0
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
lower-sqrt.f64N/A
lower-sqrt.f6428.9
Applied rewrites28.9%
Applied rewrites29.0%
Taylor expanded in b around 0
Applied rewrites29.0%
Final simplification26.2%
y-scale_m = (fabs.f64 y-scale) x-scale_m = (fabs.f64 x-scale) (FPCore (a b angle x-scale_m y-scale_m) :precision binary64 (if (<= y-scale_m 1.6e+82) (* (* (* x-scale_m (sqrt 2.0)) (sqrt 8.0)) (* 0.25 a)) (* b y-scale_m)))
y-scale_m = fabs(y_45_scale);
x-scale_m = fabs(x_45_scale);
double code(double a, double b, double angle, double x_45_scale_m, double y_45_scale_m) {
double tmp;
if (y_45_scale_m <= 1.6e+82) {
tmp = ((x_45_scale_m * sqrt(2.0)) * sqrt(8.0)) * (0.25 * a);
} else {
tmp = b * y_45_scale_m;
}
return tmp;
}
y-scale_m = abs(y_45scale)
x-scale_m = abs(x_45scale)
real(8) function code(a, b, angle, x_45scale_m, y_45scale_m)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: angle
real(8), intent (in) :: x_45scale_m
real(8), intent (in) :: y_45scale_m
real(8) :: tmp
if (y_45scale_m <= 1.6d+82) then
tmp = ((x_45scale_m * sqrt(2.0d0)) * sqrt(8.0d0)) * (0.25d0 * a)
else
tmp = b * y_45scale_m
end if
code = tmp
end function
y-scale_m = Math.abs(y_45_scale);
x-scale_m = Math.abs(x_45_scale);
public static double code(double a, double b, double angle, double x_45_scale_m, double y_45_scale_m) {
double tmp;
if (y_45_scale_m <= 1.6e+82) {
tmp = ((x_45_scale_m * Math.sqrt(2.0)) * Math.sqrt(8.0)) * (0.25 * a);
} else {
tmp = b * y_45_scale_m;
}
return tmp;
}
y-scale_m = math.fabs(y_45_scale) x-scale_m = math.fabs(x_45_scale) def code(a, b, angle, x_45_scale_m, y_45_scale_m): tmp = 0 if y_45_scale_m <= 1.6e+82: tmp = ((x_45_scale_m * math.sqrt(2.0)) * math.sqrt(8.0)) * (0.25 * a) else: tmp = b * y_45_scale_m return tmp
y-scale_m = abs(y_45_scale) x-scale_m = abs(x_45_scale) function code(a, b, angle, x_45_scale_m, y_45_scale_m) tmp = 0.0 if (y_45_scale_m <= 1.6e+82) tmp = Float64(Float64(Float64(x_45_scale_m * sqrt(2.0)) * sqrt(8.0)) * Float64(0.25 * a)); else tmp = Float64(b * y_45_scale_m); end return tmp end
y-scale_m = abs(y_45_scale); x-scale_m = abs(x_45_scale); function tmp_2 = code(a, b, angle, x_45_scale_m, y_45_scale_m) tmp = 0.0; if (y_45_scale_m <= 1.6e+82) tmp = ((x_45_scale_m * sqrt(2.0)) * sqrt(8.0)) * (0.25 * a); else tmp = b * y_45_scale_m; end tmp_2 = tmp; end
y-scale_m = N[Abs[y$45$scale], $MachinePrecision] x-scale_m = N[Abs[x$45$scale], $MachinePrecision] code[a_, b_, angle_, x$45$scale$95$m_, y$45$scale$95$m_] := If[LessEqual[y$45$scale$95$m, 1.6e+82], N[(N[(N[(x$45$scale$95$m * N[Sqrt[2.0], $MachinePrecision]), $MachinePrecision] * N[Sqrt[8.0], $MachinePrecision]), $MachinePrecision] * N[(0.25 * a), $MachinePrecision]), $MachinePrecision], N[(b * y$45$scale$95$m), $MachinePrecision]]
\begin{array}{l}
y-scale_m = \left|y-scale\right|
\\
x-scale_m = \left|x-scale\right|
\\
\begin{array}{l}
\mathbf{if}\;y-scale\_m \leq 1.6 \cdot 10^{+82}:\\
\;\;\;\;\left(\left(x-scale\_m \cdot \sqrt{2}\right) \cdot \sqrt{8}\right) \cdot \left(0.25 \cdot a\right)\\
\mathbf{else}:\\
\;\;\;\;b \cdot y-scale\_m\\
\end{array}
\end{array}
if y-scale < 1.59999999999999987e82Initial program 2.2%
Taylor expanded in b around 0
Applied rewrites6.7%
Taylor expanded in angle around 0
Applied rewrites18.9%
if 1.59999999999999987e82 < y-scale Initial program 6.9%
Taylor expanded in angle around 0
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
lower-sqrt.f64N/A
lower-sqrt.f6428.9
Applied rewrites28.9%
Applied rewrites29.0%
Taylor expanded in b around 0
Applied rewrites29.0%
Final simplification21.2%
y-scale_m = (fabs.f64 y-scale) x-scale_m = (fabs.f64 x-scale) (FPCore (a b angle x-scale_m y-scale_m) :precision binary64 (* b y-scale_m))
y-scale_m = fabs(y_45_scale);
x-scale_m = fabs(x_45_scale);
double code(double a, double b, double angle, double x_45_scale_m, double y_45_scale_m) {
return b * y_45_scale_m;
}
y-scale_m = abs(y_45scale)
x-scale_m = abs(x_45scale)
real(8) function code(a, b, angle, x_45scale_m, y_45scale_m)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: angle
real(8), intent (in) :: x_45scale_m
real(8), intent (in) :: y_45scale_m
code = b * y_45scale_m
end function
y-scale_m = Math.abs(y_45_scale);
x-scale_m = Math.abs(x_45_scale);
public static double code(double a, double b, double angle, double x_45_scale_m, double y_45_scale_m) {
return b * y_45_scale_m;
}
y-scale_m = math.fabs(y_45_scale) x-scale_m = math.fabs(x_45_scale) def code(a, b, angle, x_45_scale_m, y_45_scale_m): return b * y_45_scale_m
y-scale_m = abs(y_45_scale) x-scale_m = abs(x_45_scale) function code(a, b, angle, x_45_scale_m, y_45_scale_m) return Float64(b * y_45_scale_m) end
y-scale_m = abs(y_45_scale); x-scale_m = abs(x_45_scale); function tmp = code(a, b, angle, x_45_scale_m, y_45_scale_m) tmp = b * y_45_scale_m; end
y-scale_m = N[Abs[y$45$scale], $MachinePrecision] x-scale_m = N[Abs[x$45$scale], $MachinePrecision] code[a_, b_, angle_, x$45$scale$95$m_, y$45$scale$95$m_] := N[(b * y$45$scale$95$m), $MachinePrecision]
\begin{array}{l}
y-scale_m = \left|y-scale\right|
\\
x-scale_m = \left|x-scale\right|
\\
b \cdot y-scale\_m
\end{array}
Initial program 3.3%
Taylor expanded in angle around 0
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
lower-sqrt.f64N/A
lower-sqrt.f6419.0
Applied rewrites19.0%
Applied rewrites19.1%
Taylor expanded in b around 0
Applied rewrites19.1%
herbie shell --seed 2024284
(FPCore (a b angle x-scale y-scale)
:name "a 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))))