Example from Robby

Percentage Accurate: 99.8% → 99.8%
Time: 12.1s
Alternatives: 12
Speedup: N/A×

Specification

?
\[\begin{array}{l} \\ \begin{array}{l} t_1 := \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\\ \left|\left(ew \cdot \sin t\right) \cdot \cos t\_1 + \left(eh \cdot \cos t\right) \cdot \sin t\_1\right| \end{array} \end{array} \]
(FPCore (eh ew t)
 :precision binary64
 (let* ((t_1 (atan (/ (/ eh ew) (tan t)))))
   (fabs (+ (* (* ew (sin t)) (cos t_1)) (* (* eh (cos t)) (sin t_1))))))
double code(double eh, double ew, double t) {
	double t_1 = atan(((eh / ew) / tan(t)));
	return fabs((((ew * sin(t)) * cos(t_1)) + ((eh * cos(t)) * sin(t_1))));
}
module fmin_fmax_functions
    implicit none
    private
    public fmax
    public fmin

    interface fmax
        module procedure fmax88
        module procedure fmax44
        module procedure fmax84
        module procedure fmax48
    end interface
    interface fmin
        module procedure fmin88
        module procedure fmin44
        module procedure fmin84
        module procedure fmin48
    end interface
contains
    real(8) function fmax88(x, y) result (res)
        real(8), intent (in) :: x
        real(8), intent (in) :: y
        res = merge(y, merge(x, max(x, y), y /= y), x /= x)
    end function
    real(4) function fmax44(x, y) result (res)
        real(4), intent (in) :: x
        real(4), intent (in) :: y
        res = merge(y, merge(x, max(x, y), y /= y), x /= x)
    end function
    real(8) function fmax84(x, y) result(res)
        real(8), intent (in) :: x
        real(4), intent (in) :: y
        res = merge(dble(y), merge(x, max(x, dble(y)), y /= y), x /= x)
    end function
    real(8) function fmax48(x, y) result(res)
        real(4), intent (in) :: x
        real(8), intent (in) :: y
        res = merge(y, merge(dble(x), max(dble(x), y), y /= y), x /= x)
    end function
    real(8) function fmin88(x, y) result (res)
        real(8), intent (in) :: x
        real(8), intent (in) :: y
        res = merge(y, merge(x, min(x, y), y /= y), x /= x)
    end function
    real(4) function fmin44(x, y) result (res)
        real(4), intent (in) :: x
        real(4), intent (in) :: y
        res = merge(y, merge(x, min(x, y), y /= y), x /= x)
    end function
    real(8) function fmin84(x, y) result(res)
        real(8), intent (in) :: x
        real(4), intent (in) :: y
        res = merge(dble(y), merge(x, min(x, dble(y)), y /= y), x /= x)
    end function
    real(8) function fmin48(x, y) result(res)
        real(4), intent (in) :: x
        real(8), intent (in) :: y
        res = merge(y, merge(dble(x), min(dble(x), y), y /= y), x /= x)
    end function
end module

real(8) function code(eh, ew, t)
use fmin_fmax_functions
    real(8), intent (in) :: eh
    real(8), intent (in) :: ew
    real(8), intent (in) :: t
    real(8) :: t_1
    t_1 = atan(((eh / ew) / tan(t)))
    code = abs((((ew * sin(t)) * cos(t_1)) + ((eh * cos(t)) * sin(t_1))))
end function
public static double code(double eh, double ew, double t) {
	double t_1 = Math.atan(((eh / ew) / Math.tan(t)));
	return Math.abs((((ew * Math.sin(t)) * Math.cos(t_1)) + ((eh * Math.cos(t)) * Math.sin(t_1))));
}
def code(eh, ew, t):
	t_1 = math.atan(((eh / ew) / math.tan(t)))
	return math.fabs((((ew * math.sin(t)) * math.cos(t_1)) + ((eh * math.cos(t)) * math.sin(t_1))))
function code(eh, ew, t)
	t_1 = atan(Float64(Float64(eh / ew) / tan(t)))
	return abs(Float64(Float64(Float64(ew * sin(t)) * cos(t_1)) + Float64(Float64(eh * cos(t)) * sin(t_1))))
end
function tmp = code(eh, ew, t)
	t_1 = atan(((eh / ew) / tan(t)));
	tmp = abs((((ew * sin(t)) * cos(t_1)) + ((eh * cos(t)) * sin(t_1))));
end
code[eh_, ew_, t_] := Block[{t$95$1 = N[ArcTan[N[(N[(eh / ew), $MachinePrecision] / N[Tan[t], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]}, N[Abs[N[(N[(N[(ew * N[Sin[t], $MachinePrecision]), $MachinePrecision] * N[Cos[t$95$1], $MachinePrecision]), $MachinePrecision] + N[(N[(eh * N[Cos[t], $MachinePrecision]), $MachinePrecision] * N[Sin[t$95$1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]
\begin{array}{l}

\\
\begin{array}{l}
t_1 := \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\\
\left|\left(ew \cdot \sin t\right) \cdot \cos t\_1 + \left(eh \cdot \cos t\right) \cdot \sin t\_1\right|
\end{array}
\end{array}

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: 99.8% accurate, 1.0× speedup?

\[\begin{array}{l} \\ \begin{array}{l} t_1 := \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\\ \left|\left(ew \cdot \sin t\right) \cdot \cos t\_1 + \left(eh \cdot \cos t\right) \cdot \sin t\_1\right| \end{array} \end{array} \]
(FPCore (eh ew t)
 :precision binary64
 (let* ((t_1 (atan (/ (/ eh ew) (tan t)))))
   (fabs (+ (* (* ew (sin t)) (cos t_1)) (* (* eh (cos t)) (sin t_1))))))
double code(double eh, double ew, double t) {
	double t_1 = atan(((eh / ew) / tan(t)));
	return fabs((((ew * sin(t)) * cos(t_1)) + ((eh * cos(t)) * sin(t_1))));
}
module fmin_fmax_functions
    implicit none
    private
    public fmax
    public fmin

    interface fmax
        module procedure fmax88
        module procedure fmax44
        module procedure fmax84
        module procedure fmax48
    end interface
    interface fmin
        module procedure fmin88
        module procedure fmin44
        module procedure fmin84
        module procedure fmin48
    end interface
contains
    real(8) function fmax88(x, y) result (res)
        real(8), intent (in) :: x
        real(8), intent (in) :: y
        res = merge(y, merge(x, max(x, y), y /= y), x /= x)
    end function
    real(4) function fmax44(x, y) result (res)
        real(4), intent (in) :: x
        real(4), intent (in) :: y
        res = merge(y, merge(x, max(x, y), y /= y), x /= x)
    end function
    real(8) function fmax84(x, y) result(res)
        real(8), intent (in) :: x
        real(4), intent (in) :: y
        res = merge(dble(y), merge(x, max(x, dble(y)), y /= y), x /= x)
    end function
    real(8) function fmax48(x, y) result(res)
        real(4), intent (in) :: x
        real(8), intent (in) :: y
        res = merge(y, merge(dble(x), max(dble(x), y), y /= y), x /= x)
    end function
    real(8) function fmin88(x, y) result (res)
        real(8), intent (in) :: x
        real(8), intent (in) :: y
        res = merge(y, merge(x, min(x, y), y /= y), x /= x)
    end function
    real(4) function fmin44(x, y) result (res)
        real(4), intent (in) :: x
        real(4), intent (in) :: y
        res = merge(y, merge(x, min(x, y), y /= y), x /= x)
    end function
    real(8) function fmin84(x, y) result(res)
        real(8), intent (in) :: x
        real(4), intent (in) :: y
        res = merge(dble(y), merge(x, min(x, dble(y)), y /= y), x /= x)
    end function
    real(8) function fmin48(x, y) result(res)
        real(4), intent (in) :: x
        real(8), intent (in) :: y
        res = merge(y, merge(dble(x), min(dble(x), y), y /= y), x /= x)
    end function
end module

real(8) function code(eh, ew, t)
use fmin_fmax_functions
    real(8), intent (in) :: eh
    real(8), intent (in) :: ew
    real(8), intent (in) :: t
    real(8) :: t_1
    t_1 = atan(((eh / ew) / tan(t)))
    code = abs((((ew * sin(t)) * cos(t_1)) + ((eh * cos(t)) * sin(t_1))))
end function
public static double code(double eh, double ew, double t) {
	double t_1 = Math.atan(((eh / ew) / Math.tan(t)));
	return Math.abs((((ew * Math.sin(t)) * Math.cos(t_1)) + ((eh * Math.cos(t)) * Math.sin(t_1))));
}
def code(eh, ew, t):
	t_1 = math.atan(((eh / ew) / math.tan(t)))
	return math.fabs((((ew * math.sin(t)) * math.cos(t_1)) + ((eh * math.cos(t)) * math.sin(t_1))))
function code(eh, ew, t)
	t_1 = atan(Float64(Float64(eh / ew) / tan(t)))
	return abs(Float64(Float64(Float64(ew * sin(t)) * cos(t_1)) + Float64(Float64(eh * cos(t)) * sin(t_1))))
end
function tmp = code(eh, ew, t)
	t_1 = atan(((eh / ew) / tan(t)));
	tmp = abs((((ew * sin(t)) * cos(t_1)) + ((eh * cos(t)) * sin(t_1))));
end
code[eh_, ew_, t_] := Block[{t$95$1 = N[ArcTan[N[(N[(eh / ew), $MachinePrecision] / N[Tan[t], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]}, N[Abs[N[(N[(N[(ew * N[Sin[t], $MachinePrecision]), $MachinePrecision] * N[Cos[t$95$1], $MachinePrecision]), $MachinePrecision] + N[(N[(eh * N[Cos[t], $MachinePrecision]), $MachinePrecision] * N[Sin[t$95$1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]
\begin{array}{l}

\\
\begin{array}{l}
t_1 := \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\\
\left|\left(ew \cdot \sin t\right) \cdot \cos t\_1 + \left(eh \cdot \cos t\right) \cdot \sin t\_1\right|
\end{array}
\end{array}

Alternative 1: 99.8% accurate, 1.1× speedup?

\[\begin{array}{l} \\ \begin{array}{l} t_1 := \frac{\frac{eh}{ew}}{\tan t}\\ \left|\mathsf{fma}\left(ew, \sin t \cdot \frac{1}{\sqrt{1 + {t\_1}^{2}}}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} t\_1\right)\right| \end{array} \end{array} \]
(FPCore (eh ew t)
 :precision binary64
 (let* ((t_1 (/ (/ eh ew) (tan t))))
   (fabs
    (fma
     ew
     (* (sin t) (/ 1.0 (sqrt (+ 1.0 (pow t_1 2.0)))))
     (* (* (cos t) eh) (tanh (asinh t_1)))))))
double code(double eh, double ew, double t) {
	double t_1 = (eh / ew) / tan(t);
	return fabs(fma(ew, (sin(t) * (1.0 / sqrt((1.0 + pow(t_1, 2.0))))), ((cos(t) * eh) * tanh(asinh(t_1)))));
}
function code(eh, ew, t)
	t_1 = Float64(Float64(eh / ew) / tan(t))
	return abs(fma(ew, Float64(sin(t) * Float64(1.0 / sqrt(Float64(1.0 + (t_1 ^ 2.0))))), Float64(Float64(cos(t) * eh) * tanh(asinh(t_1)))))
end
code[eh_, ew_, t_] := Block[{t$95$1 = N[(N[(eh / ew), $MachinePrecision] / N[Tan[t], $MachinePrecision]), $MachinePrecision]}, N[Abs[N[(ew * N[(N[Sin[t], $MachinePrecision] * N[(1.0 / N[Sqrt[N[(1.0 + N[Power[t$95$1, 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(N[Cos[t], $MachinePrecision] * eh), $MachinePrecision] * N[Tanh[N[ArcSinh[t$95$1], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]
\begin{array}{l}

\\
\begin{array}{l}
t_1 := \frac{\frac{eh}{ew}}{\tan t}\\
\left|\mathsf{fma}\left(ew, \sin t \cdot \frac{1}{\sqrt{1 + {t\_1}^{2}}}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} t\_1\right)\right|
\end{array}
\end{array}
Derivation
  1. Initial program 99.8%

    \[\left|\left(ew \cdot \sin t\right) \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
  2. Add Preprocessing
  3. Step-by-step derivation
    1. lift-+.f64N/A

      \[\leadsto \left|\color{blue}{\left(ew \cdot \sin t\right) \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)}\right| \]
    2. lift-*.f64N/A

      \[\leadsto \left|\color{blue}{\left(ew \cdot \sin t\right) \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
    3. lift-*.f64N/A

      \[\leadsto \left|\color{blue}{\left(ew \cdot \sin t\right)} \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
    4. lift-sin.f64N/A

      \[\leadsto \left|\left(ew \cdot \color{blue}{\sin t}\right) \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
    5. associate-*l*N/A

      \[\leadsto \left|\color{blue}{ew \cdot \left(\sin t \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
    6. lift-*.f64N/A

      \[\leadsto \left|ew \cdot \left(\sin t \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right) + \color{blue}{\left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)}\right| \]
    7. lift-*.f64N/A

      \[\leadsto \left|ew \cdot \left(\sin t \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right) + \color{blue}{\left(eh \cdot \cos t\right)} \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
    8. lift-cos.f64N/A

      \[\leadsto \left|ew \cdot \left(\sin t \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right) + \left(eh \cdot \color{blue}{\cos t}\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
    9. lift-sin.f64N/A

      \[\leadsto \left|ew \cdot \left(\sin t \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right) + \left(eh \cdot \cos t\right) \cdot \color{blue}{\sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)}\right| \]
    10. lift-atan.f64N/A

      \[\leadsto \left|ew \cdot \left(\sin t \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right) + \left(eh \cdot \cos t\right) \cdot \sin \color{blue}{\tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)}\right| \]
    11. lift-/.f64N/A

      \[\leadsto \left|ew \cdot \left(\sin t \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \color{blue}{\left(\frac{\frac{eh}{ew}}{\tan t}\right)}\right| \]
    12. lift-/.f64N/A

      \[\leadsto \left|ew \cdot \left(\sin t \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\color{blue}{\frac{eh}{ew}}}{\tan t}\right)\right| \]
    13. lift-tan.f64N/A

      \[\leadsto \left|ew \cdot \left(\sin t \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\color{blue}{\tan t}}\right)\right| \]
  4. Applied rewrites99.8%

    \[\leadsto \left|\color{blue}{\mathsf{fma}\left(ew, \sin t \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right), \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)}\right| \]
  5. Step-by-step derivation
    1. lift-cos.f64N/A

      \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \color{blue}{\cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
    2. lift-atan.f64N/A

      \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \cos \color{blue}{\tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
    3. lift-/.f64N/A

      \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \cos \tan^{-1} \left(\frac{\color{blue}{\frac{eh}{ew}}}{\tan t}\right), \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
    4. lift-/.f64N/A

      \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \cos \tan^{-1} \color{blue}{\left(\frac{\frac{eh}{ew}}{\tan t}\right)}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
    5. lift-tan.f64N/A

      \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\color{blue}{\tan t}}\right), \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
    6. cos-atanN/A

      \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \color{blue}{\frac{1}{\sqrt{1 + \frac{\frac{eh}{ew}}{\tan t} \cdot \frac{\frac{eh}{ew}}{\tan t}}}}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
    7. unpow2N/A

      \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \frac{1}{\sqrt{1 + \color{blue}{{\left(\frac{\frac{eh}{ew}}{\tan t}\right)}^{2}}}}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
    8. lift-tan.f64N/A

      \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \frac{1}{\sqrt{1 + {\left(\frac{\frac{eh}{ew}}{\color{blue}{\tan t}}\right)}^{2}}}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
    9. lift-/.f64N/A

      \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \frac{1}{\sqrt{1 + {\color{blue}{\left(\frac{\frac{eh}{ew}}{\tan t}\right)}}^{2}}}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
    10. lift-/.f64N/A

      \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \frac{1}{\sqrt{1 + {\left(\frac{\color{blue}{\frac{eh}{ew}}}{\tan t}\right)}^{2}}}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
    11. lift-pow.f64N/A

      \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \frac{1}{\sqrt{1 + \color{blue}{{\left(\frac{\frac{eh}{ew}}{\tan t}\right)}^{2}}}}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
    12. lift-+.f64N/A

      \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \frac{1}{\sqrt{\color{blue}{1 + {\left(\frac{\frac{eh}{ew}}{\tan t}\right)}^{2}}}}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
    13. lift-sqrt.f64N/A

      \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \frac{1}{\color{blue}{\sqrt{1 + {\left(\frac{\frac{eh}{ew}}{\tan t}\right)}^{2}}}}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
    14. lift-/.f6499.8

      \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \color{blue}{\frac{1}{\sqrt{1 + {\left(\frac{\frac{eh}{ew}}{\tan t}\right)}^{2}}}}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
  6. Applied rewrites99.8%

    \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \color{blue}{\frac{1}{\sqrt{1 + {\left(\frac{\frac{eh}{ew}}{\tan t}\right)}^{2}}}}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
  7. Add Preprocessing

Alternative 2: 99.1% accurate, 1.1× speedup?

\[\begin{array}{l} \\ \left|\left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{\mathsf{fma}\left(-0.3333333333333333, \left(t \cdot t\right) \cdot eh, eh\right)}{ew}}{t}\right) + \left(ew \cdot \sin t\right) \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \end{array} \]
(FPCore (eh ew t)
 :precision binary64
 (fabs
  (+
   (*
    (* eh (cos t))
    (sin (atan (/ (/ (fma -0.3333333333333333 (* (* t t) eh) eh) ew) t))))
   (* (* ew (sin t)) (cos (atan (/ (/ eh ew) (tan t))))))))
double code(double eh, double ew, double t) {
	return fabs((((eh * cos(t)) * sin(atan(((fma(-0.3333333333333333, ((t * t) * eh), eh) / ew) / t)))) + ((ew * sin(t)) * cos(atan(((eh / ew) / tan(t)))))));
}
function code(eh, ew, t)
	return abs(Float64(Float64(Float64(eh * cos(t)) * sin(atan(Float64(Float64(fma(-0.3333333333333333, Float64(Float64(t * t) * eh), eh) / ew) / t)))) + Float64(Float64(ew * sin(t)) * cos(atan(Float64(Float64(eh / ew) / tan(t)))))))
end
code[eh_, ew_, t_] := N[Abs[N[(N[(N[(eh * N[Cos[t], $MachinePrecision]), $MachinePrecision] * N[Sin[N[ArcTan[N[(N[(N[(-0.3333333333333333 * N[(N[(t * t), $MachinePrecision] * eh), $MachinePrecision] + eh), $MachinePrecision] / ew), $MachinePrecision] / t), $MachinePrecision]], $MachinePrecision]], $MachinePrecision]), $MachinePrecision] + N[(N[(ew * N[Sin[t], $MachinePrecision]), $MachinePrecision] * N[Cos[N[ArcTan[N[(N[(eh / ew), $MachinePrecision] / N[Tan[t], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]
\begin{array}{l}

\\
\left|\left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{\mathsf{fma}\left(-0.3333333333333333, \left(t \cdot t\right) \cdot eh, eh\right)}{ew}}{t}\right) + \left(ew \cdot \sin t\right) \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right|
\end{array}
Derivation
  1. Initial program 99.8%

    \[\left|\left(ew \cdot \sin t\right) \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
  2. Add Preprocessing
  3. Taylor expanded in t around 0

    \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \color{blue}{\left(\frac{\frac{-1}{3} \cdot \frac{eh \cdot {t}^{2}}{ew} + \frac{eh}{ew}}{t}\right)}\right| \]
  4. Step-by-step derivation
    1. lower-/.f64N/A

      \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{-1}{3} \cdot \frac{eh \cdot {t}^{2}}{ew} + \frac{eh}{ew}}{\color{blue}{t}}\right)\right| \]
    2. associate-*r/N/A

      \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{\frac{-1}{3} \cdot \left(eh \cdot {t}^{2}\right)}{ew} + \frac{eh}{ew}}{t}\right)\right| \]
    3. div-add-revN/A

      \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{\frac{-1}{3} \cdot \left(eh \cdot {t}^{2}\right) + eh}{ew}}{t}\right)\right| \]
    4. lower-/.f64N/A

      \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{\frac{-1}{3} \cdot \left(eh \cdot {t}^{2}\right) + eh}{ew}}{t}\right)\right| \]
    5. lower-fma.f64N/A

      \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{\mathsf{fma}\left(\frac{-1}{3}, eh \cdot {t}^{2}, eh\right)}{ew}}{t}\right)\right| \]
    6. *-commutativeN/A

      \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{\mathsf{fma}\left(\frac{-1}{3}, {t}^{2} \cdot eh, eh\right)}{ew}}{t}\right)\right| \]
    7. lower-*.f64N/A

      \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{\mathsf{fma}\left(\frac{-1}{3}, {t}^{2} \cdot eh, eh\right)}{ew}}{t}\right)\right| \]
    8. unpow2N/A

      \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{\mathsf{fma}\left(\frac{-1}{3}, \left(t \cdot t\right) \cdot eh, eh\right)}{ew}}{t}\right)\right| \]
    9. lower-*.f6499.3

      \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{\mathsf{fma}\left(-0.3333333333333333, \left(t \cdot t\right) \cdot eh, eh\right)}{ew}}{t}\right)\right| \]
  5. Applied rewrites99.3%

    \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \color{blue}{\left(\frac{\frac{\mathsf{fma}\left(-0.3333333333333333, \left(t \cdot t\right) \cdot eh, eh\right)}{ew}}{t}\right)}\right| \]
  6. Final simplification99.3%

    \[\leadsto \left|\left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{\mathsf{fma}\left(-0.3333333333333333, \left(t \cdot t\right) \cdot eh, eh\right)}{ew}}{t}\right) + \left(ew \cdot \sin t\right) \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
  7. Add Preprocessing

Alternative 3: 87.5% accurate, 1.6× speedup?

\[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;eh \leq -6.1 \cdot 10^{+57}:\\ \;\;\;\;\left|eh \cdot \left(\cos t \cdot \sin \tan^{-1} \left(\frac{eh \cdot \cos t}{ew \cdot \sin t}\right)\right)\right|\\ \mathbf{elif}\;eh \leq 7.2 \cdot 10^{+72}:\\ \;\;\;\;\left|\mathsf{fma}\left(ew, \sin t, eh \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right|\\ \mathbf{else}:\\ \;\;\;\;\left|\left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\cos t}{ew} \cdot \frac{eh}{\sin t}\right)\right|\\ \end{array} \end{array} \]
(FPCore (eh ew t)
 :precision binary64
 (if (<= eh -6.1e+57)
   (fabs (* eh (* (cos t) (sin (atan (/ (* eh (cos t)) (* ew (sin t))))))))
   (if (<= eh 7.2e+72)
     (fabs (fma ew (sin t) (* eh (tanh (asinh (/ (/ eh ew) (tan t)))))))
     (fabs
      (* (* (cos t) eh) (tanh (asinh (* (/ (cos t) ew) (/ eh (sin t))))))))))
double code(double eh, double ew, double t) {
	double tmp;
	if (eh <= -6.1e+57) {
		tmp = fabs((eh * (cos(t) * sin(atan(((eh * cos(t)) / (ew * sin(t))))))));
	} else if (eh <= 7.2e+72) {
		tmp = fabs(fma(ew, sin(t), (eh * tanh(asinh(((eh / ew) / tan(t)))))));
	} else {
		tmp = fabs(((cos(t) * eh) * tanh(asinh(((cos(t) / ew) * (eh / sin(t)))))));
	}
	return tmp;
}
function code(eh, ew, t)
	tmp = 0.0
	if (eh <= -6.1e+57)
		tmp = abs(Float64(eh * Float64(cos(t) * sin(atan(Float64(Float64(eh * cos(t)) / Float64(ew * sin(t))))))));
	elseif (eh <= 7.2e+72)
		tmp = abs(fma(ew, sin(t), Float64(eh * tanh(asinh(Float64(Float64(eh / ew) / tan(t)))))));
	else
		tmp = abs(Float64(Float64(cos(t) * eh) * tanh(asinh(Float64(Float64(cos(t) / ew) * Float64(eh / sin(t)))))));
	end
	return tmp
end
code[eh_, ew_, t_] := If[LessEqual[eh, -6.1e+57], N[Abs[N[(eh * N[(N[Cos[t], $MachinePrecision] * N[Sin[N[ArcTan[N[(N[(eh * N[Cos[t], $MachinePrecision]), $MachinePrecision] / N[(ew * N[Sin[t], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], If[LessEqual[eh, 7.2e+72], N[Abs[N[(ew * N[Sin[t], $MachinePrecision] + N[(eh * N[Tanh[N[ArcSinh[N[(N[(eh / ew), $MachinePrecision] / N[Tan[t], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], N[Abs[N[(N[(N[Cos[t], $MachinePrecision] * eh), $MachinePrecision] * N[Tanh[N[ArcSinh[N[(N[(N[Cos[t], $MachinePrecision] / ew), $MachinePrecision] * N[(eh / N[Sin[t], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]]
\begin{array}{l}

\\
\begin{array}{l}
\mathbf{if}\;eh \leq -6.1 \cdot 10^{+57}:\\
\;\;\;\;\left|eh \cdot \left(\cos t \cdot \sin \tan^{-1} \left(\frac{eh \cdot \cos t}{ew \cdot \sin t}\right)\right)\right|\\

\mathbf{elif}\;eh \leq 7.2 \cdot 10^{+72}:\\
\;\;\;\;\left|\mathsf{fma}\left(ew, \sin t, eh \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right|\\

\mathbf{else}:\\
\;\;\;\;\left|\left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\cos t}{ew} \cdot \frac{eh}{\sin t}\right)\right|\\


\end{array}
\end{array}
Derivation
  1. Split input into 3 regimes
  2. if eh < -6.09999999999999975e57

    1. Initial program 99.9%

      \[\left|\left(ew \cdot \sin t\right) \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
    2. Add Preprocessing
    3. Taylor expanded in t around 0

      \[\leadsto \left|\color{blue}{eh \cdot \sin \tan^{-1} \left(\frac{eh \cdot \cos t}{ew \cdot \sin t}\right)}\right| \]
    4. Step-by-step derivation
      1. *-commutativeN/A

        \[\leadsto \left|\sin \tan^{-1} \left(\frac{eh \cdot \cos t}{ew \cdot \sin t}\right) \cdot \color{blue}{eh}\right| \]
      2. lower-*.f64N/A

        \[\leadsto \left|\sin \tan^{-1} \left(\frac{eh \cdot \cos t}{ew \cdot \sin t}\right) \cdot \color{blue}{eh}\right| \]
    5. Applied rewrites50.9%

      \[\leadsto \left|\color{blue}{\tanh \sinh^{-1} \left(\frac{\cos t}{ew} \cdot \frac{eh}{\sin t}\right) \cdot eh}\right| \]
    6. Taylor expanded in eh around inf

      \[\leadsto \left|\color{blue}{eh \cdot \left(\cos t \cdot \sin \tan^{-1} \left(\frac{eh \cdot \cos t}{ew \cdot \sin t}\right)\right)}\right| \]
    7. Step-by-step derivation
      1. lower-*.f64N/A

        \[\leadsto \left|eh \cdot \color{blue}{\left(\cos t \cdot \sin \tan^{-1} \left(\frac{eh \cdot \cos t}{ew \cdot \sin t}\right)\right)}\right| \]
      2. lower-*.f64N/A

        \[\leadsto \left|eh \cdot \left(\cos t \cdot \color{blue}{\sin \tan^{-1} \left(\frac{eh \cdot \cos t}{ew \cdot \sin t}\right)}\right)\right| \]
      3. lift-cos.f64N/A

        \[\leadsto \left|eh \cdot \left(\cos t \cdot \sin \color{blue}{\tan^{-1} \left(\frac{eh \cdot \cos t}{ew \cdot \sin t}\right)}\right)\right| \]
      4. lower-sin.f64N/A

        \[\leadsto \left|eh \cdot \left(\cos t \cdot \sin \tan^{-1} \left(\frac{eh \cdot \cos t}{ew \cdot \sin t}\right)\right)\right| \]
      5. lower-atan.f64N/A

        \[\leadsto \left|eh \cdot \left(\cos t \cdot \sin \tan^{-1} \left(\frac{eh \cdot \cos t}{ew \cdot \sin t}\right)\right)\right| \]
      6. lower-/.f64N/A

        \[\leadsto \left|eh \cdot \left(\cos t \cdot \sin \tan^{-1} \left(\frac{eh \cdot \cos t}{ew \cdot \sin t}\right)\right)\right| \]
      7. lift-cos.f64N/A

        \[\leadsto \left|eh \cdot \left(\cos t \cdot \sin \tan^{-1} \left(\frac{eh \cdot \cos t}{ew \cdot \sin t}\right)\right)\right| \]
      8. lift-*.f64N/A

        \[\leadsto \left|eh \cdot \left(\cos t \cdot \sin \tan^{-1} \left(\frac{eh \cdot \cos t}{ew \cdot \sin t}\right)\right)\right| \]
      9. lift-sin.f64N/A

        \[\leadsto \left|eh \cdot \left(\cos t \cdot \sin \tan^{-1} \left(\frac{eh \cdot \cos t}{ew \cdot \sin t}\right)\right)\right| \]
      10. lift-*.f6491.6

        \[\leadsto \left|eh \cdot \left(\cos t \cdot \sin \tan^{-1} \left(\frac{eh \cdot \cos t}{ew \cdot \sin t}\right)\right)\right| \]
    8. Applied rewrites91.6%

      \[\leadsto \left|\color{blue}{eh \cdot \left(\cos t \cdot \sin \tan^{-1} \left(\frac{eh \cdot \cos t}{ew \cdot \sin t}\right)\right)}\right| \]

    if -6.09999999999999975e57 < eh < 7.20000000000000069e72

    1. Initial program 99.8%

      \[\left|\left(ew \cdot \sin t\right) \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
    2. Add Preprocessing
    3. Step-by-step derivation
      1. lift-+.f64N/A

        \[\leadsto \left|\color{blue}{\left(ew \cdot \sin t\right) \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)}\right| \]
      2. lift-*.f64N/A

        \[\leadsto \left|\color{blue}{\left(ew \cdot \sin t\right) \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
      3. lift-*.f64N/A

        \[\leadsto \left|\color{blue}{\left(ew \cdot \sin t\right)} \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
      4. lift-sin.f64N/A

        \[\leadsto \left|\left(ew \cdot \color{blue}{\sin t}\right) \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
      5. associate-*l*N/A

        \[\leadsto \left|\color{blue}{ew \cdot \left(\sin t \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
      6. lift-*.f64N/A

        \[\leadsto \left|ew \cdot \left(\sin t \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right) + \color{blue}{\left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)}\right| \]
      7. lift-*.f64N/A

        \[\leadsto \left|ew \cdot \left(\sin t \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right) + \color{blue}{\left(eh \cdot \cos t\right)} \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
      8. lift-cos.f64N/A

        \[\leadsto \left|ew \cdot \left(\sin t \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right) + \left(eh \cdot \color{blue}{\cos t}\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
      9. lift-sin.f64N/A

        \[\leadsto \left|ew \cdot \left(\sin t \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right) + \left(eh \cdot \cos t\right) \cdot \color{blue}{\sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)}\right| \]
      10. lift-atan.f64N/A

        \[\leadsto \left|ew \cdot \left(\sin t \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right) + \left(eh \cdot \cos t\right) \cdot \sin \color{blue}{\tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)}\right| \]
      11. lift-/.f64N/A

        \[\leadsto \left|ew \cdot \left(\sin t \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \color{blue}{\left(\frac{\frac{eh}{ew}}{\tan t}\right)}\right| \]
      12. lift-/.f64N/A

        \[\leadsto \left|ew \cdot \left(\sin t \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\color{blue}{\frac{eh}{ew}}}{\tan t}\right)\right| \]
      13. lift-tan.f64N/A

        \[\leadsto \left|ew \cdot \left(\sin t \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\color{blue}{\tan t}}\right)\right| \]
    4. Applied rewrites99.8%

      \[\leadsto \left|\color{blue}{\mathsf{fma}\left(ew, \sin t \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right), \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)}\right| \]
    5. Step-by-step derivation
      1. lift-cos.f64N/A

        \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \color{blue}{\cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
      2. lift-atan.f64N/A

        \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \cos \color{blue}{\tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
      3. lift-/.f64N/A

        \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \cos \tan^{-1} \left(\frac{\color{blue}{\frac{eh}{ew}}}{\tan t}\right), \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
      4. lift-/.f64N/A

        \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \cos \tan^{-1} \color{blue}{\left(\frac{\frac{eh}{ew}}{\tan t}\right)}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
      5. lift-tan.f64N/A

        \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\color{blue}{\tan t}}\right), \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
      6. cos-atanN/A

        \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \color{blue}{\frac{1}{\sqrt{1 + \frac{\frac{eh}{ew}}{\tan t} \cdot \frac{\frac{eh}{ew}}{\tan t}}}}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
      7. unpow2N/A

        \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \frac{1}{\sqrt{1 + \color{blue}{{\left(\frac{\frac{eh}{ew}}{\tan t}\right)}^{2}}}}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
      8. lift-tan.f64N/A

        \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \frac{1}{\sqrt{1 + {\left(\frac{\frac{eh}{ew}}{\color{blue}{\tan t}}\right)}^{2}}}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
      9. lift-/.f64N/A

        \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \frac{1}{\sqrt{1 + {\color{blue}{\left(\frac{\frac{eh}{ew}}{\tan t}\right)}}^{2}}}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
      10. lift-/.f64N/A

        \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \frac{1}{\sqrt{1 + {\left(\frac{\color{blue}{\frac{eh}{ew}}}{\tan t}\right)}^{2}}}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
      11. lift-pow.f64N/A

        \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \frac{1}{\sqrt{1 + \color{blue}{{\left(\frac{\frac{eh}{ew}}{\tan t}\right)}^{2}}}}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
      12. lift-+.f64N/A

        \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \frac{1}{\sqrt{\color{blue}{1 + {\left(\frac{\frac{eh}{ew}}{\tan t}\right)}^{2}}}}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
      13. lift-sqrt.f64N/A

        \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \frac{1}{\color{blue}{\sqrt{1 + {\left(\frac{\frac{eh}{ew}}{\tan t}\right)}^{2}}}}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
      14. lift-/.f6499.8

        \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \color{blue}{\frac{1}{\sqrt{1 + {\left(\frac{\frac{eh}{ew}}{\tan t}\right)}^{2}}}}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
    6. Applied rewrites99.8%

      \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \color{blue}{\frac{1}{\sqrt{1 + {\left(\frac{\frac{eh}{ew}}{\tan t}\right)}^{2}}}}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
    7. Taylor expanded in t around 0

      \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \frac{1}{\sqrt{1 + {\left(\frac{\frac{eh}{ew}}{\tan t}\right)}^{2}}}, \color{blue}{eh} \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
    8. Step-by-step derivation
      1. Applied rewrites88.3%

        \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \frac{1}{\sqrt{1 + {\left(\frac{\frac{eh}{ew}}{\tan t}\right)}^{2}}}, \color{blue}{eh} \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
      2. Taylor expanded in eh around 0

        \[\leadsto \left|\mathsf{fma}\left(ew, \color{blue}{\sin t}, eh \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
      3. Step-by-step derivation
        1. unpow2N/A

          \[\leadsto \left|\mathsf{fma}\left(ew, \sin t, eh \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
        2. cos-atanN/A

          \[\leadsto \left|\mathsf{fma}\left(ew, \sin t, eh \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
        3. lift-sin.f6487.8

          \[\leadsto \left|\mathsf{fma}\left(ew, \sin t, eh \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
      4. Applied rewrites87.8%

        \[\leadsto \left|\mathsf{fma}\left(ew, \color{blue}{\sin t}, eh \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]

      if 7.20000000000000069e72 < eh

      1. Initial program 99.8%

        \[\left|\left(ew \cdot \sin t\right) \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
      2. Add Preprocessing
      3. Taylor expanded in eh around inf

        \[\leadsto \left|\color{blue}{eh \cdot \left(\cos t \cdot \sin \tan^{-1} \left(\frac{eh \cdot \cos t}{ew \cdot \sin t}\right)\right)}\right| \]
      4. Step-by-step derivation
        1. associate-*r*N/A

          \[\leadsto \left|\left(eh \cdot \cos t\right) \cdot \color{blue}{\sin \tan^{-1} \left(\frac{eh \cdot \cos t}{ew \cdot \sin t}\right)}\right| \]
        2. lower-*.f64N/A

          \[\leadsto \left|\left(eh \cdot \cos t\right) \cdot \color{blue}{\sin \tan^{-1} \left(\frac{eh \cdot \cos t}{ew \cdot \sin t}\right)}\right| \]
        3. *-commutativeN/A

          \[\leadsto \left|\left(\cos t \cdot eh\right) \cdot \sin \color{blue}{\tan^{-1} \left(\frac{eh \cdot \cos t}{ew \cdot \sin t}\right)}\right| \]
        4. lower-*.f64N/A

          \[\leadsto \left|\left(\cos t \cdot eh\right) \cdot \sin \color{blue}{\tan^{-1} \left(\frac{eh \cdot \cos t}{ew \cdot \sin t}\right)}\right| \]
        5. lift-cos.f64N/A

          \[\leadsto \left|\left(\cos t \cdot eh\right) \cdot \sin \tan^{-1} \color{blue}{\left(\frac{eh \cdot \cos t}{ew \cdot \sin t}\right)}\right| \]
        6. sin-atanN/A

          \[\leadsto \left|\left(\cos t \cdot eh\right) \cdot \frac{\frac{eh \cdot \cos t}{ew \cdot \sin t}}{\color{blue}{\sqrt{1 + \frac{eh \cdot \cos t}{ew \cdot \sin t} \cdot \frac{eh \cdot \cos t}{ew \cdot \sin t}}}}\right| \]
        7. tanh-asinh-revN/A

          \[\leadsto \left|\left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{eh \cdot \cos t}{ew \cdot \sin t}\right)\right| \]
        8. lower-tanh.f64N/A

          \[\leadsto \left|\left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{eh \cdot \cos t}{ew \cdot \sin t}\right)\right| \]
        9. lower-asinh.f64N/A

          \[\leadsto \left|\left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{eh \cdot \cos t}{ew \cdot \sin t}\right)\right| \]
        10. *-commutativeN/A

          \[\leadsto \left|\left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\cos t \cdot eh}{ew \cdot \sin t}\right)\right| \]
        11. times-fracN/A

          \[\leadsto \left|\left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\cos t}{ew} \cdot \frac{eh}{\sin t}\right)\right| \]
        12. lower-*.f64N/A

          \[\leadsto \left|\left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\cos t}{ew} \cdot \frac{eh}{\sin t}\right)\right| \]
        13. lower-/.f64N/A

          \[\leadsto \left|\left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\cos t}{ew} \cdot \frac{eh}{\sin t}\right)\right| \]
        14. lift-cos.f64N/A

          \[\leadsto \left|\left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\cos t}{ew} \cdot \frac{eh}{\sin t}\right)\right| \]
        15. lower-/.f64N/A

          \[\leadsto \left|\left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\cos t}{ew} \cdot \frac{eh}{\sin t}\right)\right| \]
      5. Applied rewrites93.2%

        \[\leadsto \left|\color{blue}{\left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\cos t}{ew} \cdot \frac{eh}{\sin t}\right)}\right| \]
    9. Recombined 3 regimes into one program.
    10. Add Preprocessing

    Alternative 4: 87.5% accurate, 1.6× speedup?

    \[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;eh \leq -6.1 \cdot 10^{+57} \lor \neg \left(eh \leq 7.2 \cdot 10^{+72}\right):\\ \;\;\;\;\left|eh \cdot \left(\cos t \cdot \sin \tan^{-1} \left(\frac{eh \cdot \cos t}{ew \cdot \sin t}\right)\right)\right|\\ \mathbf{else}:\\ \;\;\;\;\left|\mathsf{fma}\left(ew, \sin t, eh \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right|\\ \end{array} \end{array} \]
    (FPCore (eh ew t)
     :precision binary64
     (if (or (<= eh -6.1e+57) (not (<= eh 7.2e+72)))
       (fabs (* eh (* (cos t) (sin (atan (/ (* eh (cos t)) (* ew (sin t))))))))
       (fabs (fma ew (sin t) (* eh (tanh (asinh (/ (/ eh ew) (tan t)))))))))
    double code(double eh, double ew, double t) {
    	double tmp;
    	if ((eh <= -6.1e+57) || !(eh <= 7.2e+72)) {
    		tmp = fabs((eh * (cos(t) * sin(atan(((eh * cos(t)) / (ew * sin(t))))))));
    	} else {
    		tmp = fabs(fma(ew, sin(t), (eh * tanh(asinh(((eh / ew) / tan(t)))))));
    	}
    	return tmp;
    }
    
    function code(eh, ew, t)
    	tmp = 0.0
    	if ((eh <= -6.1e+57) || !(eh <= 7.2e+72))
    		tmp = abs(Float64(eh * Float64(cos(t) * sin(atan(Float64(Float64(eh * cos(t)) / Float64(ew * sin(t))))))));
    	else
    		tmp = abs(fma(ew, sin(t), Float64(eh * tanh(asinh(Float64(Float64(eh / ew) / tan(t)))))));
    	end
    	return tmp
    end
    
    code[eh_, ew_, t_] := If[Or[LessEqual[eh, -6.1e+57], N[Not[LessEqual[eh, 7.2e+72]], $MachinePrecision]], N[Abs[N[(eh * N[(N[Cos[t], $MachinePrecision] * N[Sin[N[ArcTan[N[(N[(eh * N[Cos[t], $MachinePrecision]), $MachinePrecision] / N[(ew * N[Sin[t], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], N[Abs[N[(ew * N[Sin[t], $MachinePrecision] + N[(eh * N[Tanh[N[ArcSinh[N[(N[(eh / ew), $MachinePrecision] / N[Tan[t], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]
    
    \begin{array}{l}
    
    \\
    \begin{array}{l}
    \mathbf{if}\;eh \leq -6.1 \cdot 10^{+57} \lor \neg \left(eh \leq 7.2 \cdot 10^{+72}\right):\\
    \;\;\;\;\left|eh \cdot \left(\cos t \cdot \sin \tan^{-1} \left(\frac{eh \cdot \cos t}{ew \cdot \sin t}\right)\right)\right|\\
    
    \mathbf{else}:\\
    \;\;\;\;\left|\mathsf{fma}\left(ew, \sin t, eh \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right|\\
    
    
    \end{array}
    \end{array}
    
    Derivation
    1. Split input into 2 regimes
    2. if eh < -6.09999999999999975e57 or 7.20000000000000069e72 < eh

      1. Initial program 99.8%

        \[\left|\left(ew \cdot \sin t\right) \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
      2. Add Preprocessing
      3. Taylor expanded in t around 0

        \[\leadsto \left|\color{blue}{eh \cdot \sin \tan^{-1} \left(\frac{eh \cdot \cos t}{ew \cdot \sin t}\right)}\right| \]
      4. Step-by-step derivation
        1. *-commutativeN/A

          \[\leadsto \left|\sin \tan^{-1} \left(\frac{eh \cdot \cos t}{ew \cdot \sin t}\right) \cdot \color{blue}{eh}\right| \]
        2. lower-*.f64N/A

          \[\leadsto \left|\sin \tan^{-1} \left(\frac{eh \cdot \cos t}{ew \cdot \sin t}\right) \cdot \color{blue}{eh}\right| \]
      5. Applied rewrites48.7%

        \[\leadsto \left|\color{blue}{\tanh \sinh^{-1} \left(\frac{\cos t}{ew} \cdot \frac{eh}{\sin t}\right) \cdot eh}\right| \]
      6. Taylor expanded in eh around inf

        \[\leadsto \left|\color{blue}{eh \cdot \left(\cos t \cdot \sin \tan^{-1} \left(\frac{eh \cdot \cos t}{ew \cdot \sin t}\right)\right)}\right| \]
      7. Step-by-step derivation
        1. lower-*.f64N/A

          \[\leadsto \left|eh \cdot \color{blue}{\left(\cos t \cdot \sin \tan^{-1} \left(\frac{eh \cdot \cos t}{ew \cdot \sin t}\right)\right)}\right| \]
        2. lower-*.f64N/A

          \[\leadsto \left|eh \cdot \left(\cos t \cdot \color{blue}{\sin \tan^{-1} \left(\frac{eh \cdot \cos t}{ew \cdot \sin t}\right)}\right)\right| \]
        3. lift-cos.f64N/A

          \[\leadsto \left|eh \cdot \left(\cos t \cdot \sin \color{blue}{\tan^{-1} \left(\frac{eh \cdot \cos t}{ew \cdot \sin t}\right)}\right)\right| \]
        4. lower-sin.f64N/A

          \[\leadsto \left|eh \cdot \left(\cos t \cdot \sin \tan^{-1} \left(\frac{eh \cdot \cos t}{ew \cdot \sin t}\right)\right)\right| \]
        5. lower-atan.f64N/A

          \[\leadsto \left|eh \cdot \left(\cos t \cdot \sin \tan^{-1} \left(\frac{eh \cdot \cos t}{ew \cdot \sin t}\right)\right)\right| \]
        6. lower-/.f64N/A

          \[\leadsto \left|eh \cdot \left(\cos t \cdot \sin \tan^{-1} \left(\frac{eh \cdot \cos t}{ew \cdot \sin t}\right)\right)\right| \]
        7. lift-cos.f64N/A

          \[\leadsto \left|eh \cdot \left(\cos t \cdot \sin \tan^{-1} \left(\frac{eh \cdot \cos t}{ew \cdot \sin t}\right)\right)\right| \]
        8. lift-*.f64N/A

          \[\leadsto \left|eh \cdot \left(\cos t \cdot \sin \tan^{-1} \left(\frac{eh \cdot \cos t}{ew \cdot \sin t}\right)\right)\right| \]
        9. lift-sin.f64N/A

          \[\leadsto \left|eh \cdot \left(\cos t \cdot \sin \tan^{-1} \left(\frac{eh \cdot \cos t}{ew \cdot \sin t}\right)\right)\right| \]
        10. lift-*.f6492.5

          \[\leadsto \left|eh \cdot \left(\cos t \cdot \sin \tan^{-1} \left(\frac{eh \cdot \cos t}{ew \cdot \sin t}\right)\right)\right| \]
      8. Applied rewrites92.5%

        \[\leadsto \left|\color{blue}{eh \cdot \left(\cos t \cdot \sin \tan^{-1} \left(\frac{eh \cdot \cos t}{ew \cdot \sin t}\right)\right)}\right| \]

      if -6.09999999999999975e57 < eh < 7.20000000000000069e72

      1. Initial program 99.8%

        \[\left|\left(ew \cdot \sin t\right) \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
      2. Add Preprocessing
      3. Step-by-step derivation
        1. lift-+.f64N/A

          \[\leadsto \left|\color{blue}{\left(ew \cdot \sin t\right) \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)}\right| \]
        2. lift-*.f64N/A

          \[\leadsto \left|\color{blue}{\left(ew \cdot \sin t\right) \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
        3. lift-*.f64N/A

          \[\leadsto \left|\color{blue}{\left(ew \cdot \sin t\right)} \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
        4. lift-sin.f64N/A

          \[\leadsto \left|\left(ew \cdot \color{blue}{\sin t}\right) \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
        5. associate-*l*N/A

          \[\leadsto \left|\color{blue}{ew \cdot \left(\sin t \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
        6. lift-*.f64N/A

          \[\leadsto \left|ew \cdot \left(\sin t \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right) + \color{blue}{\left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)}\right| \]
        7. lift-*.f64N/A

          \[\leadsto \left|ew \cdot \left(\sin t \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right) + \color{blue}{\left(eh \cdot \cos t\right)} \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
        8. lift-cos.f64N/A

          \[\leadsto \left|ew \cdot \left(\sin t \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right) + \left(eh \cdot \color{blue}{\cos t}\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
        9. lift-sin.f64N/A

          \[\leadsto \left|ew \cdot \left(\sin t \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right) + \left(eh \cdot \cos t\right) \cdot \color{blue}{\sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)}\right| \]
        10. lift-atan.f64N/A

          \[\leadsto \left|ew \cdot \left(\sin t \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right) + \left(eh \cdot \cos t\right) \cdot \sin \color{blue}{\tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)}\right| \]
        11. lift-/.f64N/A

          \[\leadsto \left|ew \cdot \left(\sin t \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \color{blue}{\left(\frac{\frac{eh}{ew}}{\tan t}\right)}\right| \]
        12. lift-/.f64N/A

          \[\leadsto \left|ew \cdot \left(\sin t \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\color{blue}{\frac{eh}{ew}}}{\tan t}\right)\right| \]
        13. lift-tan.f64N/A

          \[\leadsto \left|ew \cdot \left(\sin t \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\color{blue}{\tan t}}\right)\right| \]
      4. Applied rewrites99.8%

        \[\leadsto \left|\color{blue}{\mathsf{fma}\left(ew, \sin t \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right), \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)}\right| \]
      5. Step-by-step derivation
        1. lift-cos.f64N/A

          \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \color{blue}{\cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
        2. lift-atan.f64N/A

          \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \cos \color{blue}{\tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
        3. lift-/.f64N/A

          \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \cos \tan^{-1} \left(\frac{\color{blue}{\frac{eh}{ew}}}{\tan t}\right), \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
        4. lift-/.f64N/A

          \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \cos \tan^{-1} \color{blue}{\left(\frac{\frac{eh}{ew}}{\tan t}\right)}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
        5. lift-tan.f64N/A

          \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\color{blue}{\tan t}}\right), \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
        6. cos-atanN/A

          \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \color{blue}{\frac{1}{\sqrt{1 + \frac{\frac{eh}{ew}}{\tan t} \cdot \frac{\frac{eh}{ew}}{\tan t}}}}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
        7. unpow2N/A

          \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \frac{1}{\sqrt{1 + \color{blue}{{\left(\frac{\frac{eh}{ew}}{\tan t}\right)}^{2}}}}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
        8. lift-tan.f64N/A

          \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \frac{1}{\sqrt{1 + {\left(\frac{\frac{eh}{ew}}{\color{blue}{\tan t}}\right)}^{2}}}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
        9. lift-/.f64N/A

          \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \frac{1}{\sqrt{1 + {\color{blue}{\left(\frac{\frac{eh}{ew}}{\tan t}\right)}}^{2}}}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
        10. lift-/.f64N/A

          \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \frac{1}{\sqrt{1 + {\left(\frac{\color{blue}{\frac{eh}{ew}}}{\tan t}\right)}^{2}}}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
        11. lift-pow.f64N/A

          \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \frac{1}{\sqrt{1 + \color{blue}{{\left(\frac{\frac{eh}{ew}}{\tan t}\right)}^{2}}}}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
        12. lift-+.f64N/A

          \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \frac{1}{\sqrt{\color{blue}{1 + {\left(\frac{\frac{eh}{ew}}{\tan t}\right)}^{2}}}}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
        13. lift-sqrt.f64N/A

          \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \frac{1}{\color{blue}{\sqrt{1 + {\left(\frac{\frac{eh}{ew}}{\tan t}\right)}^{2}}}}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
        14. lift-/.f6499.8

          \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \color{blue}{\frac{1}{\sqrt{1 + {\left(\frac{\frac{eh}{ew}}{\tan t}\right)}^{2}}}}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
      6. Applied rewrites99.8%

        \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \color{blue}{\frac{1}{\sqrt{1 + {\left(\frac{\frac{eh}{ew}}{\tan t}\right)}^{2}}}}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
      7. Taylor expanded in t around 0

        \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \frac{1}{\sqrt{1 + {\left(\frac{\frac{eh}{ew}}{\tan t}\right)}^{2}}}, \color{blue}{eh} \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
      8. Step-by-step derivation
        1. Applied rewrites88.3%

          \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \frac{1}{\sqrt{1 + {\left(\frac{\frac{eh}{ew}}{\tan t}\right)}^{2}}}, \color{blue}{eh} \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
        2. Taylor expanded in eh around 0

          \[\leadsto \left|\mathsf{fma}\left(ew, \color{blue}{\sin t}, eh \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
        3. Step-by-step derivation
          1. unpow2N/A

            \[\leadsto \left|\mathsf{fma}\left(ew, \sin t, eh \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
          2. cos-atanN/A

            \[\leadsto \left|\mathsf{fma}\left(ew, \sin t, eh \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
          3. lift-sin.f6487.8

            \[\leadsto \left|\mathsf{fma}\left(ew, \sin t, eh \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
        4. Applied rewrites87.8%

          \[\leadsto \left|\mathsf{fma}\left(ew, \color{blue}{\sin t}, eh \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
      9. Recombined 2 regimes into one program.
      10. Final simplification89.7%

        \[\leadsto \begin{array}{l} \mathbf{if}\;eh \leq -6.1 \cdot 10^{+57} \lor \neg \left(eh \leq 7.2 \cdot 10^{+72}\right):\\ \;\;\;\;\left|eh \cdot \left(\cos t \cdot \sin \tan^{-1} \left(\frac{eh \cdot \cos t}{ew \cdot \sin t}\right)\right)\right|\\ \mathbf{else}:\\ \;\;\;\;\left|\mathsf{fma}\left(ew, \sin t, eh \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right|\\ \end{array} \]
      11. Add Preprocessing

      Alternative 5: 98.6% accurate, 1.6× speedup?

      \[\begin{array}{l} \\ \left|\mathsf{fma}\left(ew, \sin t, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \end{array} \]
      (FPCore (eh ew t)
       :precision binary64
       (fabs
        (fma ew (sin t) (* (* (cos t) eh) (tanh (asinh (/ (/ eh ew) (tan t))))))))
      double code(double eh, double ew, double t) {
      	return fabs(fma(ew, sin(t), ((cos(t) * eh) * tanh(asinh(((eh / ew) / tan(t)))))));
      }
      
      function code(eh, ew, t)
      	return abs(fma(ew, sin(t), Float64(Float64(cos(t) * eh) * tanh(asinh(Float64(Float64(eh / ew) / tan(t)))))))
      end
      
      code[eh_, ew_, t_] := N[Abs[N[(ew * N[Sin[t], $MachinePrecision] + N[(N[(N[Cos[t], $MachinePrecision] * eh), $MachinePrecision] * N[Tanh[N[ArcSinh[N[(N[(eh / ew), $MachinePrecision] / N[Tan[t], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]
      
      \begin{array}{l}
      
      \\
      \left|\mathsf{fma}\left(ew, \sin t, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right|
      \end{array}
      
      Derivation
      1. Initial program 99.8%

        \[\left|\left(ew \cdot \sin t\right) \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
      2. Add Preprocessing
      3. Step-by-step derivation
        1. lift-+.f64N/A

          \[\leadsto \left|\color{blue}{\left(ew \cdot \sin t\right) \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)}\right| \]
        2. lift-*.f64N/A

          \[\leadsto \left|\color{blue}{\left(ew \cdot \sin t\right) \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
        3. lift-*.f64N/A

          \[\leadsto \left|\color{blue}{\left(ew \cdot \sin t\right)} \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
        4. lift-sin.f64N/A

          \[\leadsto \left|\left(ew \cdot \color{blue}{\sin t}\right) \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
        5. associate-*l*N/A

          \[\leadsto \left|\color{blue}{ew \cdot \left(\sin t \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
        6. lift-*.f64N/A

          \[\leadsto \left|ew \cdot \left(\sin t \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right) + \color{blue}{\left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)}\right| \]
        7. lift-*.f64N/A

          \[\leadsto \left|ew \cdot \left(\sin t \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right) + \color{blue}{\left(eh \cdot \cos t\right)} \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
        8. lift-cos.f64N/A

          \[\leadsto \left|ew \cdot \left(\sin t \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right) + \left(eh \cdot \color{blue}{\cos t}\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
        9. lift-sin.f64N/A

          \[\leadsto \left|ew \cdot \left(\sin t \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right) + \left(eh \cdot \cos t\right) \cdot \color{blue}{\sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)}\right| \]
        10. lift-atan.f64N/A

          \[\leadsto \left|ew \cdot \left(\sin t \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right) + \left(eh \cdot \cos t\right) \cdot \sin \color{blue}{\tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)}\right| \]
        11. lift-/.f64N/A

          \[\leadsto \left|ew \cdot \left(\sin t \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \color{blue}{\left(\frac{\frac{eh}{ew}}{\tan t}\right)}\right| \]
        12. lift-/.f64N/A

          \[\leadsto \left|ew \cdot \left(\sin t \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\color{blue}{\frac{eh}{ew}}}{\tan t}\right)\right| \]
        13. lift-tan.f64N/A

          \[\leadsto \left|ew \cdot \left(\sin t \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\color{blue}{\tan t}}\right)\right| \]
      4. Applied rewrites99.8%

        \[\leadsto \left|\color{blue}{\mathsf{fma}\left(ew, \sin t \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right), \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)}\right| \]
      5. Step-by-step derivation
        1. lift-cos.f64N/A

          \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \color{blue}{\cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
        2. lift-atan.f64N/A

          \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \cos \color{blue}{\tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
        3. lift-/.f64N/A

          \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \cos \tan^{-1} \left(\frac{\color{blue}{\frac{eh}{ew}}}{\tan t}\right), \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
        4. lift-/.f64N/A

          \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \cos \tan^{-1} \color{blue}{\left(\frac{\frac{eh}{ew}}{\tan t}\right)}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
        5. lift-tan.f64N/A

          \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\color{blue}{\tan t}}\right), \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
        6. cos-atanN/A

          \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \color{blue}{\frac{1}{\sqrt{1 + \frac{\frac{eh}{ew}}{\tan t} \cdot \frac{\frac{eh}{ew}}{\tan t}}}}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
        7. unpow2N/A

          \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \frac{1}{\sqrt{1 + \color{blue}{{\left(\frac{\frac{eh}{ew}}{\tan t}\right)}^{2}}}}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
        8. lift-tan.f64N/A

          \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \frac{1}{\sqrt{1 + {\left(\frac{\frac{eh}{ew}}{\color{blue}{\tan t}}\right)}^{2}}}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
        9. lift-/.f64N/A

          \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \frac{1}{\sqrt{1 + {\color{blue}{\left(\frac{\frac{eh}{ew}}{\tan t}\right)}}^{2}}}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
        10. lift-/.f64N/A

          \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \frac{1}{\sqrt{1 + {\left(\frac{\color{blue}{\frac{eh}{ew}}}{\tan t}\right)}^{2}}}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
        11. lift-pow.f64N/A

          \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \frac{1}{\sqrt{1 + \color{blue}{{\left(\frac{\frac{eh}{ew}}{\tan t}\right)}^{2}}}}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
        12. lift-+.f64N/A

          \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \frac{1}{\sqrt{\color{blue}{1 + {\left(\frac{\frac{eh}{ew}}{\tan t}\right)}^{2}}}}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
        13. lift-sqrt.f64N/A

          \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \frac{1}{\color{blue}{\sqrt{1 + {\left(\frac{\frac{eh}{ew}}{\tan t}\right)}^{2}}}}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
        14. lift-/.f6499.8

          \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \color{blue}{\frac{1}{\sqrt{1 + {\left(\frac{\frac{eh}{ew}}{\tan t}\right)}^{2}}}}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
      6. Applied rewrites99.8%

        \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \color{blue}{\frac{1}{\sqrt{1 + {\left(\frac{\frac{eh}{ew}}{\tan t}\right)}^{2}}}}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
      7. Taylor expanded in eh around 0

        \[\leadsto \left|\mathsf{fma}\left(ew, \color{blue}{\sin t}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
      8. Step-by-step derivation
        1. unpow2N/A

          \[\leadsto \left|\mathsf{fma}\left(ew, \sin t, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
        2. cos-atanN/A

          \[\leadsto \left|\mathsf{fma}\left(ew, \sin t, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
        3. lift-sin.f6499.0

          \[\leadsto \left|\mathsf{fma}\left(ew, \sin t, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
      9. Applied rewrites99.0%

        \[\leadsto \left|\mathsf{fma}\left(ew, \color{blue}{\sin t}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
      10. Add Preprocessing

      Alternative 6: 81.0% accurate, 1.9× speedup?

      \[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;eh \leq -4.9 \cdot 10^{+163} \lor \neg \left(eh \leq 7.2 \cdot 10^{+196}\right):\\ \;\;\;\;\left|\mathsf{fma}\left(ew, \frac{ew \cdot \left(t \cdot t\right)}{eh}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{eh}{ew \cdot t}\right)\right)\right|\\ \mathbf{else}:\\ \;\;\;\;\left|\mathsf{fma}\left(ew, \sin t, eh \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right|\\ \end{array} \end{array} \]
      (FPCore (eh ew t)
       :precision binary64
       (if (or (<= eh -4.9e+163) (not (<= eh 7.2e+196)))
         (fabs
          (fma
           ew
           (/ (* ew (* t t)) eh)
           (* (* (cos t) eh) (tanh (asinh (/ eh (* ew t)))))))
         (fabs (fma ew (sin t) (* eh (tanh (asinh (/ (/ eh ew) (tan t)))))))))
      double code(double eh, double ew, double t) {
      	double tmp;
      	if ((eh <= -4.9e+163) || !(eh <= 7.2e+196)) {
      		tmp = fabs(fma(ew, ((ew * (t * t)) / eh), ((cos(t) * eh) * tanh(asinh((eh / (ew * t)))))));
      	} else {
      		tmp = fabs(fma(ew, sin(t), (eh * tanh(asinh(((eh / ew) / tan(t)))))));
      	}
      	return tmp;
      }
      
      function code(eh, ew, t)
      	tmp = 0.0
      	if ((eh <= -4.9e+163) || !(eh <= 7.2e+196))
      		tmp = abs(fma(ew, Float64(Float64(ew * Float64(t * t)) / eh), Float64(Float64(cos(t) * eh) * tanh(asinh(Float64(eh / Float64(ew * t)))))));
      	else
      		tmp = abs(fma(ew, sin(t), Float64(eh * tanh(asinh(Float64(Float64(eh / ew) / tan(t)))))));
      	end
      	return tmp
      end
      
      code[eh_, ew_, t_] := If[Or[LessEqual[eh, -4.9e+163], N[Not[LessEqual[eh, 7.2e+196]], $MachinePrecision]], N[Abs[N[(ew * N[(N[(ew * N[(t * t), $MachinePrecision]), $MachinePrecision] / eh), $MachinePrecision] + N[(N[(N[Cos[t], $MachinePrecision] * eh), $MachinePrecision] * N[Tanh[N[ArcSinh[N[(eh / N[(ew * t), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], N[Abs[N[(ew * N[Sin[t], $MachinePrecision] + N[(eh * N[Tanh[N[ArcSinh[N[(N[(eh / ew), $MachinePrecision] / N[Tan[t], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]
      
      \begin{array}{l}
      
      \\
      \begin{array}{l}
      \mathbf{if}\;eh \leq -4.9 \cdot 10^{+163} \lor \neg \left(eh \leq 7.2 \cdot 10^{+196}\right):\\
      \;\;\;\;\left|\mathsf{fma}\left(ew, \frac{ew \cdot \left(t \cdot t\right)}{eh}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{eh}{ew \cdot t}\right)\right)\right|\\
      
      \mathbf{else}:\\
      \;\;\;\;\left|\mathsf{fma}\left(ew, \sin t, eh \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right|\\
      
      
      \end{array}
      \end{array}
      
      Derivation
      1. Split input into 2 regimes
      2. if eh < -4.9e163 or 7.20000000000000015e196 < eh

        1. Initial program 99.8%

          \[\left|\left(ew \cdot \sin t\right) \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
        2. Add Preprocessing
        3. Step-by-step derivation
          1. lift-+.f64N/A

            \[\leadsto \left|\color{blue}{\left(ew \cdot \sin t\right) \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)}\right| \]
          2. lift-*.f64N/A

            \[\leadsto \left|\color{blue}{\left(ew \cdot \sin t\right) \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
          3. lift-*.f64N/A

            \[\leadsto \left|\color{blue}{\left(ew \cdot \sin t\right)} \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
          4. lift-sin.f64N/A

            \[\leadsto \left|\left(ew \cdot \color{blue}{\sin t}\right) \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
          5. associate-*l*N/A

            \[\leadsto \left|\color{blue}{ew \cdot \left(\sin t \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
          6. lift-*.f64N/A

            \[\leadsto \left|ew \cdot \left(\sin t \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right) + \color{blue}{\left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)}\right| \]
          7. lift-*.f64N/A

            \[\leadsto \left|ew \cdot \left(\sin t \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right) + \color{blue}{\left(eh \cdot \cos t\right)} \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
          8. lift-cos.f64N/A

            \[\leadsto \left|ew \cdot \left(\sin t \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right) + \left(eh \cdot \color{blue}{\cos t}\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
          9. lift-sin.f64N/A

            \[\leadsto \left|ew \cdot \left(\sin t \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right) + \left(eh \cdot \cos t\right) \cdot \color{blue}{\sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)}\right| \]
          10. lift-atan.f64N/A

            \[\leadsto \left|ew \cdot \left(\sin t \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right) + \left(eh \cdot \cos t\right) \cdot \sin \color{blue}{\tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)}\right| \]
          11. lift-/.f64N/A

            \[\leadsto \left|ew \cdot \left(\sin t \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \color{blue}{\left(\frac{\frac{eh}{ew}}{\tan t}\right)}\right| \]
          12. lift-/.f64N/A

            \[\leadsto \left|ew \cdot \left(\sin t \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\color{blue}{\frac{eh}{ew}}}{\tan t}\right)\right| \]
          13. lift-tan.f64N/A

            \[\leadsto \left|ew \cdot \left(\sin t \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\color{blue}{\tan t}}\right)\right| \]
        4. Applied rewrites99.8%

          \[\leadsto \left|\color{blue}{\mathsf{fma}\left(ew, \sin t \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right), \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)}\right| \]
        5. Step-by-step derivation
          1. lift-cos.f64N/A

            \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \color{blue}{\cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
          2. lift-atan.f64N/A

            \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \cos \color{blue}{\tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
          3. lift-/.f64N/A

            \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \cos \tan^{-1} \left(\frac{\color{blue}{\frac{eh}{ew}}}{\tan t}\right), \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
          4. lift-/.f64N/A

            \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \cos \tan^{-1} \color{blue}{\left(\frac{\frac{eh}{ew}}{\tan t}\right)}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
          5. lift-tan.f64N/A

            \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\color{blue}{\tan t}}\right), \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
          6. cos-atanN/A

            \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \color{blue}{\frac{1}{\sqrt{1 + \frac{\frac{eh}{ew}}{\tan t} \cdot \frac{\frac{eh}{ew}}{\tan t}}}}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
          7. unpow2N/A

            \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \frac{1}{\sqrt{1 + \color{blue}{{\left(\frac{\frac{eh}{ew}}{\tan t}\right)}^{2}}}}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
          8. lift-tan.f64N/A

            \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \frac{1}{\sqrt{1 + {\left(\frac{\frac{eh}{ew}}{\color{blue}{\tan t}}\right)}^{2}}}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
          9. lift-/.f64N/A

            \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \frac{1}{\sqrt{1 + {\color{blue}{\left(\frac{\frac{eh}{ew}}{\tan t}\right)}}^{2}}}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
          10. lift-/.f64N/A

            \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \frac{1}{\sqrt{1 + {\left(\frac{\color{blue}{\frac{eh}{ew}}}{\tan t}\right)}^{2}}}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
          11. lift-pow.f64N/A

            \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \frac{1}{\sqrt{1 + \color{blue}{{\left(\frac{\frac{eh}{ew}}{\tan t}\right)}^{2}}}}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
          12. lift-+.f64N/A

            \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \frac{1}{\sqrt{\color{blue}{1 + {\left(\frac{\frac{eh}{ew}}{\tan t}\right)}^{2}}}}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
          13. lift-sqrt.f64N/A

            \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \frac{1}{\color{blue}{\sqrt{1 + {\left(\frac{\frac{eh}{ew}}{\tan t}\right)}^{2}}}}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
          14. lift-/.f6499.8

            \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \color{blue}{\frac{1}{\sqrt{1 + {\left(\frac{\frac{eh}{ew}}{\tan t}\right)}^{2}}}}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
        6. Applied rewrites99.8%

          \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \color{blue}{\frac{1}{\sqrt{1 + {\left(\frac{\frac{eh}{ew}}{\tan t}\right)}^{2}}}}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
        7. Taylor expanded in t around 0

          \[\leadsto \left|\mathsf{fma}\left(ew, \color{blue}{\frac{ew \cdot {t}^{2}}{eh}}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
        8. Step-by-step derivation
          1. unpow2N/A

            \[\leadsto \left|\mathsf{fma}\left(ew, \frac{ew \cdot {t}^{2}}{eh}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
          2. cos-atanN/A

            \[\leadsto \left|\mathsf{fma}\left(ew, \frac{ew \cdot \color{blue}{{t}^{2}}}{eh}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
          3. lower-/.f64N/A

            \[\leadsto \left|\mathsf{fma}\left(ew, \frac{ew \cdot {t}^{2}}{\color{blue}{eh}}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
          4. lower-*.f64N/A

            \[\leadsto \left|\mathsf{fma}\left(ew, \frac{ew \cdot {t}^{2}}{eh}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
          5. pow2N/A

            \[\leadsto \left|\mathsf{fma}\left(ew, \frac{ew \cdot \left(t \cdot t\right)}{eh}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
          6. lift-*.f6469.1

            \[\leadsto \left|\mathsf{fma}\left(ew, \frac{ew \cdot \left(t \cdot t\right)}{eh}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
        9. Applied rewrites69.1%

          \[\leadsto \left|\mathsf{fma}\left(ew, \color{blue}{\frac{ew \cdot \left(t \cdot t\right)}{eh}}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
        10. Taylor expanded in t around 0

          \[\leadsto \left|\mathsf{fma}\left(ew, \frac{ew \cdot \left(t \cdot t\right)}{eh}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \color{blue}{\left(\frac{eh}{ew \cdot t}\right)}\right)\right| \]
        11. Step-by-step derivation
          1. lower-/.f64N/A

            \[\leadsto \left|\mathsf{fma}\left(ew, \frac{ew \cdot \left(t \cdot t\right)}{eh}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{eh}{\color{blue}{ew \cdot t}}\right)\right)\right| \]
          2. lower-*.f6469.1

            \[\leadsto \left|\mathsf{fma}\left(ew, \frac{ew \cdot \left(t \cdot t\right)}{eh}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{eh}{ew \cdot \color{blue}{t}}\right)\right)\right| \]
        12. Applied rewrites69.1%

          \[\leadsto \left|\mathsf{fma}\left(ew, \frac{ew \cdot \left(t \cdot t\right)}{eh}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \color{blue}{\left(\frac{eh}{ew \cdot t}\right)}\right)\right| \]

        if -4.9e163 < eh < 7.20000000000000015e196

        1. Initial program 99.8%

          \[\left|\left(ew \cdot \sin t\right) \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
        2. Add Preprocessing
        3. Step-by-step derivation
          1. lift-+.f64N/A

            \[\leadsto \left|\color{blue}{\left(ew \cdot \sin t\right) \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)}\right| \]
          2. lift-*.f64N/A

            \[\leadsto \left|\color{blue}{\left(ew \cdot \sin t\right) \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
          3. lift-*.f64N/A

            \[\leadsto \left|\color{blue}{\left(ew \cdot \sin t\right)} \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
          4. lift-sin.f64N/A

            \[\leadsto \left|\left(ew \cdot \color{blue}{\sin t}\right) \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
          5. associate-*l*N/A

            \[\leadsto \left|\color{blue}{ew \cdot \left(\sin t \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
          6. lift-*.f64N/A

            \[\leadsto \left|ew \cdot \left(\sin t \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right) + \color{blue}{\left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)}\right| \]
          7. lift-*.f64N/A

            \[\leadsto \left|ew \cdot \left(\sin t \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right) + \color{blue}{\left(eh \cdot \cos t\right)} \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
          8. lift-cos.f64N/A

            \[\leadsto \left|ew \cdot \left(\sin t \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right) + \left(eh \cdot \color{blue}{\cos t}\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
          9. lift-sin.f64N/A

            \[\leadsto \left|ew \cdot \left(\sin t \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right) + \left(eh \cdot \cos t\right) \cdot \color{blue}{\sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)}\right| \]
          10. lift-atan.f64N/A

            \[\leadsto \left|ew \cdot \left(\sin t \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right) + \left(eh \cdot \cos t\right) \cdot \sin \color{blue}{\tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)}\right| \]
          11. lift-/.f64N/A

            \[\leadsto \left|ew \cdot \left(\sin t \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \color{blue}{\left(\frac{\frac{eh}{ew}}{\tan t}\right)}\right| \]
          12. lift-/.f64N/A

            \[\leadsto \left|ew \cdot \left(\sin t \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\color{blue}{\frac{eh}{ew}}}{\tan t}\right)\right| \]
          13. lift-tan.f64N/A

            \[\leadsto \left|ew \cdot \left(\sin t \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\color{blue}{\tan t}}\right)\right| \]
        4. Applied rewrites99.8%

          \[\leadsto \left|\color{blue}{\mathsf{fma}\left(ew, \sin t \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right), \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)}\right| \]
        5. Step-by-step derivation
          1. lift-cos.f64N/A

            \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \color{blue}{\cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
          2. lift-atan.f64N/A

            \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \cos \color{blue}{\tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
          3. lift-/.f64N/A

            \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \cos \tan^{-1} \left(\frac{\color{blue}{\frac{eh}{ew}}}{\tan t}\right), \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
          4. lift-/.f64N/A

            \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \cos \tan^{-1} \color{blue}{\left(\frac{\frac{eh}{ew}}{\tan t}\right)}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
          5. lift-tan.f64N/A

            \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\color{blue}{\tan t}}\right), \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
          6. cos-atanN/A

            \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \color{blue}{\frac{1}{\sqrt{1 + \frac{\frac{eh}{ew}}{\tan t} \cdot \frac{\frac{eh}{ew}}{\tan t}}}}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
          7. unpow2N/A

            \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \frac{1}{\sqrt{1 + \color{blue}{{\left(\frac{\frac{eh}{ew}}{\tan t}\right)}^{2}}}}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
          8. lift-tan.f64N/A

            \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \frac{1}{\sqrt{1 + {\left(\frac{\frac{eh}{ew}}{\color{blue}{\tan t}}\right)}^{2}}}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
          9. lift-/.f64N/A

            \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \frac{1}{\sqrt{1 + {\color{blue}{\left(\frac{\frac{eh}{ew}}{\tan t}\right)}}^{2}}}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
          10. lift-/.f64N/A

            \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \frac{1}{\sqrt{1 + {\left(\frac{\color{blue}{\frac{eh}{ew}}}{\tan t}\right)}^{2}}}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
          11. lift-pow.f64N/A

            \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \frac{1}{\sqrt{1 + \color{blue}{{\left(\frac{\frac{eh}{ew}}{\tan t}\right)}^{2}}}}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
          12. lift-+.f64N/A

            \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \frac{1}{\sqrt{\color{blue}{1 + {\left(\frac{\frac{eh}{ew}}{\tan t}\right)}^{2}}}}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
          13. lift-sqrt.f64N/A

            \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \frac{1}{\color{blue}{\sqrt{1 + {\left(\frac{\frac{eh}{ew}}{\tan t}\right)}^{2}}}}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
          14. lift-/.f6499.8

            \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \color{blue}{\frac{1}{\sqrt{1 + {\left(\frac{\frac{eh}{ew}}{\tan t}\right)}^{2}}}}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
        6. Applied rewrites99.8%

          \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \color{blue}{\frac{1}{\sqrt{1 + {\left(\frac{\frac{eh}{ew}}{\tan t}\right)}^{2}}}}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
        7. Taylor expanded in t around 0

          \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \frac{1}{\sqrt{1 + {\left(\frac{\frac{eh}{ew}}{\tan t}\right)}^{2}}}, \color{blue}{eh} \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
        8. Step-by-step derivation
          1. Applied rewrites81.3%

            \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \frac{1}{\sqrt{1 + {\left(\frac{\frac{eh}{ew}}{\tan t}\right)}^{2}}}, \color{blue}{eh} \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
          2. Taylor expanded in eh around 0

            \[\leadsto \left|\mathsf{fma}\left(ew, \color{blue}{\sin t}, eh \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
          3. Step-by-step derivation
            1. unpow2N/A

              \[\leadsto \left|\mathsf{fma}\left(ew, \sin t, eh \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
            2. cos-atanN/A

              \[\leadsto \left|\mathsf{fma}\left(ew, \sin t, eh \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
            3. lift-sin.f6480.9

              \[\leadsto \left|\mathsf{fma}\left(ew, \sin t, eh \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
          4. Applied rewrites80.9%

            \[\leadsto \left|\mathsf{fma}\left(ew, \color{blue}{\sin t}, eh \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
        9. Recombined 2 regimes into one program.
        10. Final simplification78.1%

          \[\leadsto \begin{array}{l} \mathbf{if}\;eh \leq -4.9 \cdot 10^{+163} \lor \neg \left(eh \leq 7.2 \cdot 10^{+196}\right):\\ \;\;\;\;\left|\mathsf{fma}\left(ew, \frac{ew \cdot \left(t \cdot t\right)}{eh}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{eh}{ew \cdot t}\right)\right)\right|\\ \mathbf{else}:\\ \;\;\;\;\left|\mathsf{fma}\left(ew, \sin t, eh \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right|\\ \end{array} \]
        11. Add Preprocessing

        Alternative 7: 63.4% accurate, 2.4× speedup?

        \[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;eh \leq -1.15 \cdot 10^{+55} \lor \neg \left(eh \leq 2.35 \cdot 10^{-67}\right):\\ \;\;\;\;\left|\mathsf{fma}\left(ew, \frac{ew \cdot \left(t \cdot t\right)}{eh}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{eh}{ew \cdot t}\right)\right)\right|\\ \mathbf{else}:\\ \;\;\;\;\left|ew \cdot \sin t\right|\\ \end{array} \end{array} \]
        (FPCore (eh ew t)
         :precision binary64
         (if (or (<= eh -1.15e+55) (not (<= eh 2.35e-67)))
           (fabs
            (fma
             ew
             (/ (* ew (* t t)) eh)
             (* (* (cos t) eh) (tanh (asinh (/ eh (* ew t)))))))
           (fabs (* ew (sin t)))))
        double code(double eh, double ew, double t) {
        	double tmp;
        	if ((eh <= -1.15e+55) || !(eh <= 2.35e-67)) {
        		tmp = fabs(fma(ew, ((ew * (t * t)) / eh), ((cos(t) * eh) * tanh(asinh((eh / (ew * t)))))));
        	} else {
        		tmp = fabs((ew * sin(t)));
        	}
        	return tmp;
        }
        
        function code(eh, ew, t)
        	tmp = 0.0
        	if ((eh <= -1.15e+55) || !(eh <= 2.35e-67))
        		tmp = abs(fma(ew, Float64(Float64(ew * Float64(t * t)) / eh), Float64(Float64(cos(t) * eh) * tanh(asinh(Float64(eh / Float64(ew * t)))))));
        	else
        		tmp = abs(Float64(ew * sin(t)));
        	end
        	return tmp
        end
        
        code[eh_, ew_, t_] := If[Or[LessEqual[eh, -1.15e+55], N[Not[LessEqual[eh, 2.35e-67]], $MachinePrecision]], N[Abs[N[(ew * N[(N[(ew * N[(t * t), $MachinePrecision]), $MachinePrecision] / eh), $MachinePrecision] + N[(N[(N[Cos[t], $MachinePrecision] * eh), $MachinePrecision] * N[Tanh[N[ArcSinh[N[(eh / N[(ew * t), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], N[Abs[N[(ew * N[Sin[t], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]
        
        \begin{array}{l}
        
        \\
        \begin{array}{l}
        \mathbf{if}\;eh \leq -1.15 \cdot 10^{+55} \lor \neg \left(eh \leq 2.35 \cdot 10^{-67}\right):\\
        \;\;\;\;\left|\mathsf{fma}\left(ew, \frac{ew \cdot \left(t \cdot t\right)}{eh}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{eh}{ew \cdot t}\right)\right)\right|\\
        
        \mathbf{else}:\\
        \;\;\;\;\left|ew \cdot \sin t\right|\\
        
        
        \end{array}
        \end{array}
        
        Derivation
        1. Split input into 2 regimes
        2. if eh < -1.14999999999999994e55 or 2.35000000000000002e-67 < eh

          1. Initial program 99.8%

            \[\left|\left(ew \cdot \sin t\right) \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
          2. Add Preprocessing
          3. Step-by-step derivation
            1. lift-+.f64N/A

              \[\leadsto \left|\color{blue}{\left(ew \cdot \sin t\right) \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)}\right| \]
            2. lift-*.f64N/A

              \[\leadsto \left|\color{blue}{\left(ew \cdot \sin t\right) \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
            3. lift-*.f64N/A

              \[\leadsto \left|\color{blue}{\left(ew \cdot \sin t\right)} \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
            4. lift-sin.f64N/A

              \[\leadsto \left|\left(ew \cdot \color{blue}{\sin t}\right) \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
            5. associate-*l*N/A

              \[\leadsto \left|\color{blue}{ew \cdot \left(\sin t \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
            6. lift-*.f64N/A

              \[\leadsto \left|ew \cdot \left(\sin t \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right) + \color{blue}{\left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)}\right| \]
            7. lift-*.f64N/A

              \[\leadsto \left|ew \cdot \left(\sin t \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right) + \color{blue}{\left(eh \cdot \cos t\right)} \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
            8. lift-cos.f64N/A

              \[\leadsto \left|ew \cdot \left(\sin t \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right) + \left(eh \cdot \color{blue}{\cos t}\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
            9. lift-sin.f64N/A

              \[\leadsto \left|ew \cdot \left(\sin t \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right) + \left(eh \cdot \cos t\right) \cdot \color{blue}{\sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)}\right| \]
            10. lift-atan.f64N/A

              \[\leadsto \left|ew \cdot \left(\sin t \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right) + \left(eh \cdot \cos t\right) \cdot \sin \color{blue}{\tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)}\right| \]
            11. lift-/.f64N/A

              \[\leadsto \left|ew \cdot \left(\sin t \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \color{blue}{\left(\frac{\frac{eh}{ew}}{\tan t}\right)}\right| \]
            12. lift-/.f64N/A

              \[\leadsto \left|ew \cdot \left(\sin t \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\color{blue}{\frac{eh}{ew}}}{\tan t}\right)\right| \]
            13. lift-tan.f64N/A

              \[\leadsto \left|ew \cdot \left(\sin t \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\color{blue}{\tan t}}\right)\right| \]
          4. Applied rewrites99.8%

            \[\leadsto \left|\color{blue}{\mathsf{fma}\left(ew, \sin t \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right), \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)}\right| \]
          5. Step-by-step derivation
            1. lift-cos.f64N/A

              \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \color{blue}{\cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
            2. lift-atan.f64N/A

              \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \cos \color{blue}{\tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
            3. lift-/.f64N/A

              \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \cos \tan^{-1} \left(\frac{\color{blue}{\frac{eh}{ew}}}{\tan t}\right), \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
            4. lift-/.f64N/A

              \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \cos \tan^{-1} \color{blue}{\left(\frac{\frac{eh}{ew}}{\tan t}\right)}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
            5. lift-tan.f64N/A

              \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\color{blue}{\tan t}}\right), \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
            6. cos-atanN/A

              \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \color{blue}{\frac{1}{\sqrt{1 + \frac{\frac{eh}{ew}}{\tan t} \cdot \frac{\frac{eh}{ew}}{\tan t}}}}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
            7. unpow2N/A

              \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \frac{1}{\sqrt{1 + \color{blue}{{\left(\frac{\frac{eh}{ew}}{\tan t}\right)}^{2}}}}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
            8. lift-tan.f64N/A

              \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \frac{1}{\sqrt{1 + {\left(\frac{\frac{eh}{ew}}{\color{blue}{\tan t}}\right)}^{2}}}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
            9. lift-/.f64N/A

              \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \frac{1}{\sqrt{1 + {\color{blue}{\left(\frac{\frac{eh}{ew}}{\tan t}\right)}}^{2}}}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
            10. lift-/.f64N/A

              \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \frac{1}{\sqrt{1 + {\left(\frac{\color{blue}{\frac{eh}{ew}}}{\tan t}\right)}^{2}}}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
            11. lift-pow.f64N/A

              \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \frac{1}{\sqrt{1 + \color{blue}{{\left(\frac{\frac{eh}{ew}}{\tan t}\right)}^{2}}}}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
            12. lift-+.f64N/A

              \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \frac{1}{\sqrt{\color{blue}{1 + {\left(\frac{\frac{eh}{ew}}{\tan t}\right)}^{2}}}}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
            13. lift-sqrt.f64N/A

              \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \frac{1}{\color{blue}{\sqrt{1 + {\left(\frac{\frac{eh}{ew}}{\tan t}\right)}^{2}}}}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
            14. lift-/.f6499.8

              \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \color{blue}{\frac{1}{\sqrt{1 + {\left(\frac{\frac{eh}{ew}}{\tan t}\right)}^{2}}}}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
          6. Applied rewrites99.8%

            \[\leadsto \left|\mathsf{fma}\left(ew, \sin t \cdot \color{blue}{\frac{1}{\sqrt{1 + {\left(\frac{\frac{eh}{ew}}{\tan t}\right)}^{2}}}}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
          7. Taylor expanded in t around 0

            \[\leadsto \left|\mathsf{fma}\left(ew, \color{blue}{\frac{ew \cdot {t}^{2}}{eh}}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
          8. Step-by-step derivation
            1. unpow2N/A

              \[\leadsto \left|\mathsf{fma}\left(ew, \frac{ew \cdot {t}^{2}}{eh}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
            2. cos-atanN/A

              \[\leadsto \left|\mathsf{fma}\left(ew, \frac{ew \cdot \color{blue}{{t}^{2}}}{eh}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
            3. lower-/.f64N/A

              \[\leadsto \left|\mathsf{fma}\left(ew, \frac{ew \cdot {t}^{2}}{\color{blue}{eh}}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
            4. lower-*.f64N/A

              \[\leadsto \left|\mathsf{fma}\left(ew, \frac{ew \cdot {t}^{2}}{eh}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
            5. pow2N/A

              \[\leadsto \left|\mathsf{fma}\left(ew, \frac{ew \cdot \left(t \cdot t\right)}{eh}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
            6. lift-*.f6461.2

              \[\leadsto \left|\mathsf{fma}\left(ew, \frac{ew \cdot \left(t \cdot t\right)}{eh}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
          9. Applied rewrites61.2%

            \[\leadsto \left|\mathsf{fma}\left(ew, \color{blue}{\frac{ew \cdot \left(t \cdot t\right)}{eh}}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right)\right| \]
          10. Taylor expanded in t around 0

            \[\leadsto \left|\mathsf{fma}\left(ew, \frac{ew \cdot \left(t \cdot t\right)}{eh}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \color{blue}{\left(\frac{eh}{ew \cdot t}\right)}\right)\right| \]
          11. Step-by-step derivation
            1. lower-/.f64N/A

              \[\leadsto \left|\mathsf{fma}\left(ew, \frac{ew \cdot \left(t \cdot t\right)}{eh}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{eh}{\color{blue}{ew \cdot t}}\right)\right)\right| \]
            2. lower-*.f6461.2

              \[\leadsto \left|\mathsf{fma}\left(ew, \frac{ew \cdot \left(t \cdot t\right)}{eh}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{eh}{ew \cdot \color{blue}{t}}\right)\right)\right| \]
          12. Applied rewrites61.2%

            \[\leadsto \left|\mathsf{fma}\left(ew, \frac{ew \cdot \left(t \cdot t\right)}{eh}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \color{blue}{\left(\frac{eh}{ew \cdot t}\right)}\right)\right| \]

          if -1.14999999999999994e55 < eh < 2.35000000000000002e-67

          1. Initial program 99.7%

            \[\left|\left(ew \cdot \sin t\right) \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
          2. Add Preprocessing
          3. Step-by-step derivation
            1. lift-cos.f64N/A

              \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \color{blue}{\cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
            2. lift-atan.f64N/A

              \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \cos \color{blue}{\tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
            3. lift-/.f64N/A

              \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \cos \tan^{-1} \color{blue}{\left(\frac{\frac{eh}{ew}}{\tan t}\right)} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
            4. lift-/.f64N/A

              \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \cos \tan^{-1} \left(\frac{\color{blue}{\frac{eh}{ew}}}{\tan t}\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
            5. lift-tan.f64N/A

              \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\color{blue}{\tan t}}\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
            6. cos-atanN/A

              \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \color{blue}{\frac{1}{\sqrt{1 + \frac{\frac{eh}{ew}}{\tan t} \cdot \frac{\frac{eh}{ew}}{\tan t}}}} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
            7. lower-/.f64N/A

              \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \color{blue}{\frac{1}{\sqrt{1 + \frac{\frac{eh}{ew}}{\tan t} \cdot \frac{\frac{eh}{ew}}{\tan t}}}} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
            8. lower-sqrt.f64N/A

              \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \frac{1}{\color{blue}{\sqrt{1 + \frac{\frac{eh}{ew}}{\tan t} \cdot \frac{\frac{eh}{ew}}{\tan t}}}} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
            9. lower-+.f64N/A

              \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \frac{1}{\sqrt{\color{blue}{1 + \frac{\frac{eh}{ew}}{\tan t} \cdot \frac{\frac{eh}{ew}}{\tan t}}}} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
            10. pow2N/A

              \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \frac{1}{\sqrt{1 + \color{blue}{{\left(\frac{\frac{eh}{ew}}{\tan t}\right)}^{2}}}} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
            11. lower-pow.f64N/A

              \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \frac{1}{\sqrt{1 + \color{blue}{{\left(\frac{\frac{eh}{ew}}{\tan t}\right)}^{2}}}} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
            12. lift-/.f64N/A

              \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \frac{1}{\sqrt{1 + {\left(\frac{\color{blue}{\frac{eh}{ew}}}{\tan t}\right)}^{2}}} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
            13. lift-tan.f64N/A

              \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \frac{1}{\sqrt{1 + {\left(\frac{\frac{eh}{ew}}{\color{blue}{\tan t}}\right)}^{2}}} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
            14. lift-/.f6499.7

              \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \frac{1}{\sqrt{1 + {\color{blue}{\left(\frac{\frac{eh}{ew}}{\tan t}\right)}}^{2}}} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
          4. Applied rewrites99.7%

            \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \color{blue}{\frac{1}{\sqrt{1 + {\left(\frac{\frac{eh}{ew}}{\tan t}\right)}^{2}}}} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
          5. Taylor expanded in eh around 0

            \[\leadsto \left|\color{blue}{ew \cdot \sin t}\right| \]
          6. Step-by-step derivation
            1. unpow2N/A

              \[\leadsto \left|ew \cdot \sin t\right| \]
            2. cos-atanN/A

              \[\leadsto \left|ew \cdot \sin t\right| \]
            3. lift-sin.f64N/A

              \[\leadsto \left|ew \cdot \sin t\right| \]
            4. lift-*.f6468.8

              \[\leadsto \left|ew \cdot \color{blue}{\sin t}\right| \]
          7. Applied rewrites68.8%

            \[\leadsto \left|\color{blue}{ew \cdot \sin t}\right| \]
        3. Recombined 2 regimes into one program.
        4. Final simplification64.8%

          \[\leadsto \begin{array}{l} \mathbf{if}\;eh \leq -1.15 \cdot 10^{+55} \lor \neg \left(eh \leq 2.35 \cdot 10^{-67}\right):\\ \;\;\;\;\left|\mathsf{fma}\left(ew, \frac{ew \cdot \left(t \cdot t\right)}{eh}, \left(\cos t \cdot eh\right) \cdot \tanh \sinh^{-1} \left(\frac{eh}{ew \cdot t}\right)\right)\right|\\ \mathbf{else}:\\ \;\;\;\;\left|ew \cdot \sin t\right|\\ \end{array} \]
        5. Add Preprocessing

        Alternative 8: 58.8% accurate, 3.3× speedup?

        \[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;eh \leq -1.15 \cdot 10^{+55} \lor \neg \left(eh \leq 1.9 \cdot 10^{-79}\right):\\ \;\;\;\;\left|\tanh \sinh^{-1} \left(\frac{1}{ew} \cdot \frac{eh + 0.16666666666666666 \cdot \left(eh \cdot \left(t \cdot t\right)\right)}{t}\right) \cdot eh\right|\\ \mathbf{else}:\\ \;\;\;\;\left|ew \cdot \sin t\right|\\ \end{array} \end{array} \]
        (FPCore (eh ew t)
         :precision binary64
         (if (or (<= eh -1.15e+55) (not (<= eh 1.9e-79)))
           (fabs
            (*
             (tanh
              (asinh
               (* (/ 1.0 ew) (/ (+ eh (* 0.16666666666666666 (* eh (* t t)))) t))))
             eh))
           (fabs (* ew (sin t)))))
        double code(double eh, double ew, double t) {
        	double tmp;
        	if ((eh <= -1.15e+55) || !(eh <= 1.9e-79)) {
        		tmp = fabs((tanh(asinh(((1.0 / ew) * ((eh + (0.16666666666666666 * (eh * (t * t)))) / t)))) * eh));
        	} else {
        		tmp = fabs((ew * sin(t)));
        	}
        	return tmp;
        }
        
        def code(eh, ew, t):
        	tmp = 0
        	if (eh <= -1.15e+55) or not (eh <= 1.9e-79):
        		tmp = math.fabs((math.tanh(math.asinh(((1.0 / ew) * ((eh + (0.16666666666666666 * (eh * (t * t)))) / t)))) * eh))
        	else:
        		tmp = math.fabs((ew * math.sin(t)))
        	return tmp
        
        function code(eh, ew, t)
        	tmp = 0.0
        	if ((eh <= -1.15e+55) || !(eh <= 1.9e-79))
        		tmp = abs(Float64(tanh(asinh(Float64(Float64(1.0 / ew) * Float64(Float64(eh + Float64(0.16666666666666666 * Float64(eh * Float64(t * t)))) / t)))) * eh));
        	else
        		tmp = abs(Float64(ew * sin(t)));
        	end
        	return tmp
        end
        
        function tmp_2 = code(eh, ew, t)
        	tmp = 0.0;
        	if ((eh <= -1.15e+55) || ~((eh <= 1.9e-79)))
        		tmp = abs((tanh(asinh(((1.0 / ew) * ((eh + (0.16666666666666666 * (eh * (t * t)))) / t)))) * eh));
        	else
        		tmp = abs((ew * sin(t)));
        	end
        	tmp_2 = tmp;
        end
        
        code[eh_, ew_, t_] := If[Or[LessEqual[eh, -1.15e+55], N[Not[LessEqual[eh, 1.9e-79]], $MachinePrecision]], N[Abs[N[(N[Tanh[N[ArcSinh[N[(N[(1.0 / ew), $MachinePrecision] * N[(N[(eh + N[(0.16666666666666666 * N[(eh * N[(t * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], $MachinePrecision] * eh), $MachinePrecision]], $MachinePrecision], N[Abs[N[(ew * N[Sin[t], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]
        
        \begin{array}{l}
        
        \\
        \begin{array}{l}
        \mathbf{if}\;eh \leq -1.15 \cdot 10^{+55} \lor \neg \left(eh \leq 1.9 \cdot 10^{-79}\right):\\
        \;\;\;\;\left|\tanh \sinh^{-1} \left(\frac{1}{ew} \cdot \frac{eh + 0.16666666666666666 \cdot \left(eh \cdot \left(t \cdot t\right)\right)}{t}\right) \cdot eh\right|\\
        
        \mathbf{else}:\\
        \;\;\;\;\left|ew \cdot \sin t\right|\\
        
        
        \end{array}
        \end{array}
        
        Derivation
        1. Split input into 2 regimes
        2. if eh < -1.14999999999999994e55 or 1.9000000000000001e-79 < eh

          1. Initial program 99.8%

            \[\left|\left(ew \cdot \sin t\right) \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
          2. Add Preprocessing
          3. Taylor expanded in t around 0

            \[\leadsto \left|\color{blue}{eh \cdot \sin \tan^{-1} \left(\frac{eh \cdot \cos t}{ew \cdot \sin t}\right)}\right| \]
          4. Step-by-step derivation
            1. *-commutativeN/A

              \[\leadsto \left|\sin \tan^{-1} \left(\frac{eh \cdot \cos t}{ew \cdot \sin t}\right) \cdot \color{blue}{eh}\right| \]
            2. lower-*.f64N/A

              \[\leadsto \left|\sin \tan^{-1} \left(\frac{eh \cdot \cos t}{ew \cdot \sin t}\right) \cdot \color{blue}{eh}\right| \]
          5. Applied rewrites49.6%

            \[\leadsto \left|\color{blue}{\tanh \sinh^{-1} \left(\frac{\cos t}{ew} \cdot \frac{eh}{\sin t}\right) \cdot eh}\right| \]
          6. Taylor expanded in t around 0

            \[\leadsto \left|\tanh \sinh^{-1} \left(\frac{1}{ew} \cdot \frac{eh}{\sin t}\right) \cdot eh\right| \]
          7. Step-by-step derivation
            1. Applied rewrites49.6%

              \[\leadsto \left|\tanh \sinh^{-1} \left(\frac{1}{ew} \cdot \frac{eh}{\sin t}\right) \cdot eh\right| \]
            2. Taylor expanded in t around 0

              \[\leadsto \left|\tanh \sinh^{-1} \left(\frac{1}{ew} \cdot \frac{eh + \frac{1}{6} \cdot \left(eh \cdot {t}^{2}\right)}{t}\right) \cdot eh\right| \]
            3. Step-by-step derivation
              1. lower-/.f64N/A

                \[\leadsto \left|\tanh \sinh^{-1} \left(\frac{1}{ew} \cdot \frac{eh + \frac{1}{6} \cdot \left(eh \cdot {t}^{2}\right)}{t}\right) \cdot eh\right| \]
              2. lower-+.f64N/A

                \[\leadsto \left|\tanh \sinh^{-1} \left(\frac{1}{ew} \cdot \frac{eh + \frac{1}{6} \cdot \left(eh \cdot {t}^{2}\right)}{t}\right) \cdot eh\right| \]
              3. lower-*.f64N/A

                \[\leadsto \left|\tanh \sinh^{-1} \left(\frac{1}{ew} \cdot \frac{eh + \frac{1}{6} \cdot \left(eh \cdot {t}^{2}\right)}{t}\right) \cdot eh\right| \]
              4. lower-*.f64N/A

                \[\leadsto \left|\tanh \sinh^{-1} \left(\frac{1}{ew} \cdot \frac{eh + \frac{1}{6} \cdot \left(eh \cdot {t}^{2}\right)}{t}\right) \cdot eh\right| \]
              5. pow2N/A

                \[\leadsto \left|\tanh \sinh^{-1} \left(\frac{1}{ew} \cdot \frac{eh + \frac{1}{6} \cdot \left(eh \cdot \left(t \cdot t\right)\right)}{t}\right) \cdot eh\right| \]
              6. lift-*.f6449.8

                \[\leadsto \left|\tanh \sinh^{-1} \left(\frac{1}{ew} \cdot \frac{eh + 0.16666666666666666 \cdot \left(eh \cdot \left(t \cdot t\right)\right)}{t}\right) \cdot eh\right| \]
            4. Applied rewrites49.8%

              \[\leadsto \left|\tanh \sinh^{-1} \left(\frac{1}{ew} \cdot \frac{eh + 0.16666666666666666 \cdot \left(eh \cdot \left(t \cdot t\right)\right)}{t}\right) \cdot eh\right| \]

            if -1.14999999999999994e55 < eh < 1.9000000000000001e-79

            1. Initial program 99.8%

              \[\left|\left(ew \cdot \sin t\right) \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
            2. Add Preprocessing
            3. Step-by-step derivation
              1. lift-cos.f64N/A

                \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \color{blue}{\cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
              2. lift-atan.f64N/A

                \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \cos \color{blue}{\tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
              3. lift-/.f64N/A

                \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \cos \tan^{-1} \color{blue}{\left(\frac{\frac{eh}{ew}}{\tan t}\right)} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
              4. lift-/.f64N/A

                \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \cos \tan^{-1} \left(\frac{\color{blue}{\frac{eh}{ew}}}{\tan t}\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
              5. lift-tan.f64N/A

                \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\color{blue}{\tan t}}\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
              6. cos-atanN/A

                \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \color{blue}{\frac{1}{\sqrt{1 + \frac{\frac{eh}{ew}}{\tan t} \cdot \frac{\frac{eh}{ew}}{\tan t}}}} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
              7. lower-/.f64N/A

                \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \color{blue}{\frac{1}{\sqrt{1 + \frac{\frac{eh}{ew}}{\tan t} \cdot \frac{\frac{eh}{ew}}{\tan t}}}} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
              8. lower-sqrt.f64N/A

                \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \frac{1}{\color{blue}{\sqrt{1 + \frac{\frac{eh}{ew}}{\tan t} \cdot \frac{\frac{eh}{ew}}{\tan t}}}} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
              9. lower-+.f64N/A

                \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \frac{1}{\sqrt{\color{blue}{1 + \frac{\frac{eh}{ew}}{\tan t} \cdot \frac{\frac{eh}{ew}}{\tan t}}}} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
              10. pow2N/A

                \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \frac{1}{\sqrt{1 + \color{blue}{{\left(\frac{\frac{eh}{ew}}{\tan t}\right)}^{2}}}} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
              11. lower-pow.f64N/A

                \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \frac{1}{\sqrt{1 + \color{blue}{{\left(\frac{\frac{eh}{ew}}{\tan t}\right)}^{2}}}} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
              12. lift-/.f64N/A

                \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \frac{1}{\sqrt{1 + {\left(\frac{\color{blue}{\frac{eh}{ew}}}{\tan t}\right)}^{2}}} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
              13. lift-tan.f64N/A

                \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \frac{1}{\sqrt{1 + {\left(\frac{\frac{eh}{ew}}{\color{blue}{\tan t}}\right)}^{2}}} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
              14. lift-/.f6499.8

                \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \frac{1}{\sqrt{1 + {\color{blue}{\left(\frac{\frac{eh}{ew}}{\tan t}\right)}}^{2}}} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
            4. Applied rewrites99.8%

              \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \color{blue}{\frac{1}{\sqrt{1 + {\left(\frac{\frac{eh}{ew}}{\tan t}\right)}^{2}}}} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
            5. Taylor expanded in eh around 0

              \[\leadsto \left|\color{blue}{ew \cdot \sin t}\right| \]
            6. Step-by-step derivation
              1. unpow2N/A

                \[\leadsto \left|ew \cdot \sin t\right| \]
              2. cos-atanN/A

                \[\leadsto \left|ew \cdot \sin t\right| \]
              3. lift-sin.f64N/A

                \[\leadsto \left|ew \cdot \sin t\right| \]
              4. lift-*.f6470.0

                \[\leadsto \left|ew \cdot \color{blue}{\sin t}\right| \]
            7. Applied rewrites70.0%

              \[\leadsto \left|\color{blue}{ew \cdot \sin t}\right| \]
          8. Recombined 2 regimes into one program.
          9. Final simplification59.2%

            \[\leadsto \begin{array}{l} \mathbf{if}\;eh \leq -1.15 \cdot 10^{+55} \lor \neg \left(eh \leq 1.9 \cdot 10^{-79}\right):\\ \;\;\;\;\left|\tanh \sinh^{-1} \left(\frac{1}{ew} \cdot \frac{eh + 0.16666666666666666 \cdot \left(eh \cdot \left(t \cdot t\right)\right)}{t}\right) \cdot eh\right|\\ \mathbf{else}:\\ \;\;\;\;\left|ew \cdot \sin t\right|\\ \end{array} \]
          10. Add Preprocessing

          Alternative 9: 58.1% accurate, 3.5× speedup?

          \[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;eh \leq -1.15 \cdot 10^{+55} \lor \neg \left(eh \leq 2.7 \cdot 10^{-58}\right):\\ \;\;\;\;\left|\tanh \sinh^{-1} \left(\frac{1}{ew} \cdot \frac{eh}{t}\right) \cdot eh\right|\\ \mathbf{else}:\\ \;\;\;\;\left|ew \cdot \sin t\right|\\ \end{array} \end{array} \]
          (FPCore (eh ew t)
           :precision binary64
           (if (or (<= eh -1.15e+55) (not (<= eh 2.7e-58)))
             (fabs (* (tanh (asinh (* (/ 1.0 ew) (/ eh t)))) eh))
             (fabs (* ew (sin t)))))
          double code(double eh, double ew, double t) {
          	double tmp;
          	if ((eh <= -1.15e+55) || !(eh <= 2.7e-58)) {
          		tmp = fabs((tanh(asinh(((1.0 / ew) * (eh / t)))) * eh));
          	} else {
          		tmp = fabs((ew * sin(t)));
          	}
          	return tmp;
          }
          
          def code(eh, ew, t):
          	tmp = 0
          	if (eh <= -1.15e+55) or not (eh <= 2.7e-58):
          		tmp = math.fabs((math.tanh(math.asinh(((1.0 / ew) * (eh / t)))) * eh))
          	else:
          		tmp = math.fabs((ew * math.sin(t)))
          	return tmp
          
          function code(eh, ew, t)
          	tmp = 0.0
          	if ((eh <= -1.15e+55) || !(eh <= 2.7e-58))
          		tmp = abs(Float64(tanh(asinh(Float64(Float64(1.0 / ew) * Float64(eh / t)))) * eh));
          	else
          		tmp = abs(Float64(ew * sin(t)));
          	end
          	return tmp
          end
          
          function tmp_2 = code(eh, ew, t)
          	tmp = 0.0;
          	if ((eh <= -1.15e+55) || ~((eh <= 2.7e-58)))
          		tmp = abs((tanh(asinh(((1.0 / ew) * (eh / t)))) * eh));
          	else
          		tmp = abs((ew * sin(t)));
          	end
          	tmp_2 = tmp;
          end
          
          code[eh_, ew_, t_] := If[Or[LessEqual[eh, -1.15e+55], N[Not[LessEqual[eh, 2.7e-58]], $MachinePrecision]], N[Abs[N[(N[Tanh[N[ArcSinh[N[(N[(1.0 / ew), $MachinePrecision] * N[(eh / t), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], $MachinePrecision] * eh), $MachinePrecision]], $MachinePrecision], N[Abs[N[(ew * N[Sin[t], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]
          
          \begin{array}{l}
          
          \\
          \begin{array}{l}
          \mathbf{if}\;eh \leq -1.15 \cdot 10^{+55} \lor \neg \left(eh \leq 2.7 \cdot 10^{-58}\right):\\
          \;\;\;\;\left|\tanh \sinh^{-1} \left(\frac{1}{ew} \cdot \frac{eh}{t}\right) \cdot eh\right|\\
          
          \mathbf{else}:\\
          \;\;\;\;\left|ew \cdot \sin t\right|\\
          
          
          \end{array}
          \end{array}
          
          Derivation
          1. Split input into 2 regimes
          2. if eh < -1.14999999999999994e55 or 2.6999999999999999e-58 < eh

            1. Initial program 99.8%

              \[\left|\left(ew \cdot \sin t\right) \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
            2. Add Preprocessing
            3. Taylor expanded in t around 0

              \[\leadsto \left|\color{blue}{eh \cdot \sin \tan^{-1} \left(\frac{eh \cdot \cos t}{ew \cdot \sin t}\right)}\right| \]
            4. Step-by-step derivation
              1. *-commutativeN/A

                \[\leadsto \left|\sin \tan^{-1} \left(\frac{eh \cdot \cos t}{ew \cdot \sin t}\right) \cdot \color{blue}{eh}\right| \]
              2. lower-*.f64N/A

                \[\leadsto \left|\sin \tan^{-1} \left(\frac{eh \cdot \cos t}{ew \cdot \sin t}\right) \cdot \color{blue}{eh}\right| \]
            5. Applied rewrites50.1%

              \[\leadsto \left|\color{blue}{\tanh \sinh^{-1} \left(\frac{\cos t}{ew} \cdot \frac{eh}{\sin t}\right) \cdot eh}\right| \]
            6. Taylor expanded in t around 0

              \[\leadsto \left|\tanh \sinh^{-1} \left(\frac{1}{ew} \cdot \frac{eh}{\sin t}\right) \cdot eh\right| \]
            7. Step-by-step derivation
              1. Applied rewrites50.1%

                \[\leadsto \left|\tanh \sinh^{-1} \left(\frac{1}{ew} \cdot \frac{eh}{\sin t}\right) \cdot eh\right| \]
              2. Taylor expanded in t around 0

                \[\leadsto \left|\tanh \sinh^{-1} \left(\frac{1}{ew} \cdot \frac{eh}{t}\right) \cdot eh\right| \]
              3. Step-by-step derivation
                1. Applied rewrites47.8%

                  \[\leadsto \left|\tanh \sinh^{-1} \left(\frac{1}{ew} \cdot \frac{eh}{t}\right) \cdot eh\right| \]

                if -1.14999999999999994e55 < eh < 2.6999999999999999e-58

                1. Initial program 99.7%

                  \[\left|\left(ew \cdot \sin t\right) \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
                2. Add Preprocessing
                3. Step-by-step derivation
                  1. lift-cos.f64N/A

                    \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \color{blue}{\cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
                  2. lift-atan.f64N/A

                    \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \cos \color{blue}{\tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
                  3. lift-/.f64N/A

                    \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \cos \tan^{-1} \color{blue}{\left(\frac{\frac{eh}{ew}}{\tan t}\right)} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
                  4. lift-/.f64N/A

                    \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \cos \tan^{-1} \left(\frac{\color{blue}{\frac{eh}{ew}}}{\tan t}\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
                  5. lift-tan.f64N/A

                    \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\color{blue}{\tan t}}\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
                  6. cos-atanN/A

                    \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \color{blue}{\frac{1}{\sqrt{1 + \frac{\frac{eh}{ew}}{\tan t} \cdot \frac{\frac{eh}{ew}}{\tan t}}}} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
                  7. lower-/.f64N/A

                    \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \color{blue}{\frac{1}{\sqrt{1 + \frac{\frac{eh}{ew}}{\tan t} \cdot \frac{\frac{eh}{ew}}{\tan t}}}} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
                  8. lower-sqrt.f64N/A

                    \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \frac{1}{\color{blue}{\sqrt{1 + \frac{\frac{eh}{ew}}{\tan t} \cdot \frac{\frac{eh}{ew}}{\tan t}}}} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
                  9. lower-+.f64N/A

                    \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \frac{1}{\sqrt{\color{blue}{1 + \frac{\frac{eh}{ew}}{\tan t} \cdot \frac{\frac{eh}{ew}}{\tan t}}}} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
                  10. pow2N/A

                    \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \frac{1}{\sqrt{1 + \color{blue}{{\left(\frac{\frac{eh}{ew}}{\tan t}\right)}^{2}}}} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
                  11. lower-pow.f64N/A

                    \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \frac{1}{\sqrt{1 + \color{blue}{{\left(\frac{\frac{eh}{ew}}{\tan t}\right)}^{2}}}} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
                  12. lift-/.f64N/A

                    \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \frac{1}{\sqrt{1 + {\left(\frac{\color{blue}{\frac{eh}{ew}}}{\tan t}\right)}^{2}}} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
                  13. lift-tan.f64N/A

                    \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \frac{1}{\sqrt{1 + {\left(\frac{\frac{eh}{ew}}{\color{blue}{\tan t}}\right)}^{2}}} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
                  14. lift-/.f6499.7

                    \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \frac{1}{\sqrt{1 + {\color{blue}{\left(\frac{\frac{eh}{ew}}{\tan t}\right)}}^{2}}} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
                4. Applied rewrites99.7%

                  \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \color{blue}{\frac{1}{\sqrt{1 + {\left(\frac{\frac{eh}{ew}}{\tan t}\right)}^{2}}}} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
                5. Taylor expanded in eh around 0

                  \[\leadsto \left|\color{blue}{ew \cdot \sin t}\right| \]
                6. Step-by-step derivation
                  1. unpow2N/A

                    \[\leadsto \left|ew \cdot \sin t\right| \]
                  2. cos-atanN/A

                    \[\leadsto \left|ew \cdot \sin t\right| \]
                  3. lift-sin.f64N/A

                    \[\leadsto \left|ew \cdot \sin t\right| \]
                  4. lift-*.f6468.8

                    \[\leadsto \left|ew \cdot \color{blue}{\sin t}\right| \]
                7. Applied rewrites68.8%

                  \[\leadsto \left|\color{blue}{ew \cdot \sin t}\right| \]
              4. Recombined 2 regimes into one program.
              5. Final simplification57.9%

                \[\leadsto \begin{array}{l} \mathbf{if}\;eh \leq -1.15 \cdot 10^{+55} \lor \neg \left(eh \leq 2.7 \cdot 10^{-58}\right):\\ \;\;\;\;\left|\tanh \sinh^{-1} \left(\frac{1}{ew} \cdot \frac{eh}{t}\right) \cdot eh\right|\\ \mathbf{else}:\\ \;\;\;\;\left|ew \cdot \sin t\right|\\ \end{array} \]
              6. Add Preprocessing

              Alternative 10: 58.0% accurate, 3.7× speedup?

              \[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;eh \leq -1.15 \cdot 10^{+55} \lor \neg \left(eh \leq 2.7 \cdot 10^{-58}\right):\\ \;\;\;\;\left|\tanh \sinh^{-1} \left(\frac{eh}{ew \cdot t}\right) \cdot eh\right|\\ \mathbf{else}:\\ \;\;\;\;\left|ew \cdot \sin t\right|\\ \end{array} \end{array} \]
              (FPCore (eh ew t)
               :precision binary64
               (if (or (<= eh -1.15e+55) (not (<= eh 2.7e-58)))
                 (fabs (* (tanh (asinh (/ eh (* ew t)))) eh))
                 (fabs (* ew (sin t)))))
              double code(double eh, double ew, double t) {
              	double tmp;
              	if ((eh <= -1.15e+55) || !(eh <= 2.7e-58)) {
              		tmp = fabs((tanh(asinh((eh / (ew * t)))) * eh));
              	} else {
              		tmp = fabs((ew * sin(t)));
              	}
              	return tmp;
              }
              
              def code(eh, ew, t):
              	tmp = 0
              	if (eh <= -1.15e+55) or not (eh <= 2.7e-58):
              		tmp = math.fabs((math.tanh(math.asinh((eh / (ew * t)))) * eh))
              	else:
              		tmp = math.fabs((ew * math.sin(t)))
              	return tmp
              
              function code(eh, ew, t)
              	tmp = 0.0
              	if ((eh <= -1.15e+55) || !(eh <= 2.7e-58))
              		tmp = abs(Float64(tanh(asinh(Float64(eh / Float64(ew * t)))) * eh));
              	else
              		tmp = abs(Float64(ew * sin(t)));
              	end
              	return tmp
              end
              
              function tmp_2 = code(eh, ew, t)
              	tmp = 0.0;
              	if ((eh <= -1.15e+55) || ~((eh <= 2.7e-58)))
              		tmp = abs((tanh(asinh((eh / (ew * t)))) * eh));
              	else
              		tmp = abs((ew * sin(t)));
              	end
              	tmp_2 = tmp;
              end
              
              code[eh_, ew_, t_] := If[Or[LessEqual[eh, -1.15e+55], N[Not[LessEqual[eh, 2.7e-58]], $MachinePrecision]], N[Abs[N[(N[Tanh[N[ArcSinh[N[(eh / N[(ew * t), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], $MachinePrecision] * eh), $MachinePrecision]], $MachinePrecision], N[Abs[N[(ew * N[Sin[t], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]
              
              \begin{array}{l}
              
              \\
              \begin{array}{l}
              \mathbf{if}\;eh \leq -1.15 \cdot 10^{+55} \lor \neg \left(eh \leq 2.7 \cdot 10^{-58}\right):\\
              \;\;\;\;\left|\tanh \sinh^{-1} \left(\frac{eh}{ew \cdot t}\right) \cdot eh\right|\\
              
              \mathbf{else}:\\
              \;\;\;\;\left|ew \cdot \sin t\right|\\
              
              
              \end{array}
              \end{array}
              
              Derivation
              1. Split input into 2 regimes
              2. if eh < -1.14999999999999994e55 or 2.6999999999999999e-58 < eh

                1. Initial program 99.8%

                  \[\left|\left(ew \cdot \sin t\right) \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
                2. Add Preprocessing
                3. Taylor expanded in t around 0

                  \[\leadsto \left|\color{blue}{eh \cdot \sin \tan^{-1} \left(\frac{eh \cdot \cos t}{ew \cdot \sin t}\right)}\right| \]
                4. Step-by-step derivation
                  1. *-commutativeN/A

                    \[\leadsto \left|\sin \tan^{-1} \left(\frac{eh \cdot \cos t}{ew \cdot \sin t}\right) \cdot \color{blue}{eh}\right| \]
                  2. lower-*.f64N/A

                    \[\leadsto \left|\sin \tan^{-1} \left(\frac{eh \cdot \cos t}{ew \cdot \sin t}\right) \cdot \color{blue}{eh}\right| \]
                5. Applied rewrites50.1%

                  \[\leadsto \left|\color{blue}{\tanh \sinh^{-1} \left(\frac{\cos t}{ew} \cdot \frac{eh}{\sin t}\right) \cdot eh}\right| \]
                6. Taylor expanded in t around 0

                  \[\leadsto \left|\tanh \sinh^{-1} \left(\frac{eh}{ew \cdot t}\right) \cdot eh\right| \]
                7. Step-by-step derivation
                  1. lower-/.f64N/A

                    \[\leadsto \left|\tanh \sinh^{-1} \left(\frac{eh}{ew \cdot t}\right) \cdot eh\right| \]
                  2. lower-*.f6447.6

                    \[\leadsto \left|\tanh \sinh^{-1} \left(\frac{eh}{ew \cdot t}\right) \cdot eh\right| \]
                8. Applied rewrites47.6%

                  \[\leadsto \left|\tanh \sinh^{-1} \left(\frac{eh}{ew \cdot t}\right) \cdot eh\right| \]

                if -1.14999999999999994e55 < eh < 2.6999999999999999e-58

                1. Initial program 99.7%

                  \[\left|\left(ew \cdot \sin t\right) \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
                2. Add Preprocessing
                3. Step-by-step derivation
                  1. lift-cos.f64N/A

                    \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \color{blue}{\cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
                  2. lift-atan.f64N/A

                    \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \cos \color{blue}{\tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
                  3. lift-/.f64N/A

                    \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \cos \tan^{-1} \color{blue}{\left(\frac{\frac{eh}{ew}}{\tan t}\right)} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
                  4. lift-/.f64N/A

                    \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \cos \tan^{-1} \left(\frac{\color{blue}{\frac{eh}{ew}}}{\tan t}\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
                  5. lift-tan.f64N/A

                    \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\color{blue}{\tan t}}\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
                  6. cos-atanN/A

                    \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \color{blue}{\frac{1}{\sqrt{1 + \frac{\frac{eh}{ew}}{\tan t} \cdot \frac{\frac{eh}{ew}}{\tan t}}}} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
                  7. lower-/.f64N/A

                    \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \color{blue}{\frac{1}{\sqrt{1 + \frac{\frac{eh}{ew}}{\tan t} \cdot \frac{\frac{eh}{ew}}{\tan t}}}} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
                  8. lower-sqrt.f64N/A

                    \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \frac{1}{\color{blue}{\sqrt{1 + \frac{\frac{eh}{ew}}{\tan t} \cdot \frac{\frac{eh}{ew}}{\tan t}}}} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
                  9. lower-+.f64N/A

                    \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \frac{1}{\sqrt{\color{blue}{1 + \frac{\frac{eh}{ew}}{\tan t} \cdot \frac{\frac{eh}{ew}}{\tan t}}}} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
                  10. pow2N/A

                    \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \frac{1}{\sqrt{1 + \color{blue}{{\left(\frac{\frac{eh}{ew}}{\tan t}\right)}^{2}}}} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
                  11. lower-pow.f64N/A

                    \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \frac{1}{\sqrt{1 + \color{blue}{{\left(\frac{\frac{eh}{ew}}{\tan t}\right)}^{2}}}} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
                  12. lift-/.f64N/A

                    \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \frac{1}{\sqrt{1 + {\left(\frac{\color{blue}{\frac{eh}{ew}}}{\tan t}\right)}^{2}}} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
                  13. lift-tan.f64N/A

                    \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \frac{1}{\sqrt{1 + {\left(\frac{\frac{eh}{ew}}{\color{blue}{\tan t}}\right)}^{2}}} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
                  14. lift-/.f6499.7

                    \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \frac{1}{\sqrt{1 + {\color{blue}{\left(\frac{\frac{eh}{ew}}{\tan t}\right)}}^{2}}} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
                4. Applied rewrites99.7%

                  \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \color{blue}{\frac{1}{\sqrt{1 + {\left(\frac{\frac{eh}{ew}}{\tan t}\right)}^{2}}}} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
                5. Taylor expanded in eh around 0

                  \[\leadsto \left|\color{blue}{ew \cdot \sin t}\right| \]
                6. Step-by-step derivation
                  1. unpow2N/A

                    \[\leadsto \left|ew \cdot \sin t\right| \]
                  2. cos-atanN/A

                    \[\leadsto \left|ew \cdot \sin t\right| \]
                  3. lift-sin.f64N/A

                    \[\leadsto \left|ew \cdot \sin t\right| \]
                  4. lift-*.f6468.8

                    \[\leadsto \left|ew \cdot \color{blue}{\sin t}\right| \]
                7. Applied rewrites68.8%

                  \[\leadsto \left|\color{blue}{ew \cdot \sin t}\right| \]
              3. Recombined 2 regimes into one program.
              4. Final simplification57.8%

                \[\leadsto \begin{array}{l} \mathbf{if}\;eh \leq -1.15 \cdot 10^{+55} \lor \neg \left(eh \leq 2.7 \cdot 10^{-58}\right):\\ \;\;\;\;\left|\tanh \sinh^{-1} \left(\frac{eh}{ew \cdot t}\right) \cdot eh\right|\\ \mathbf{else}:\\ \;\;\;\;\left|ew \cdot \sin t\right|\\ \end{array} \]
              5. Add Preprocessing

              Alternative 11: 42.4% accurate, 8.1× speedup?

              \[\begin{array}{l} \\ \left|ew \cdot \sin t\right| \end{array} \]
              (FPCore (eh ew t) :precision binary64 (fabs (* ew (sin t))))
              double code(double eh, double ew, double t) {
              	return fabs((ew * sin(t)));
              }
              
              module fmin_fmax_functions
                  implicit none
                  private
                  public fmax
                  public fmin
              
                  interface fmax
                      module procedure fmax88
                      module procedure fmax44
                      module procedure fmax84
                      module procedure fmax48
                  end interface
                  interface fmin
                      module procedure fmin88
                      module procedure fmin44
                      module procedure fmin84
                      module procedure fmin48
                  end interface
              contains
                  real(8) function fmax88(x, y) result (res)
                      real(8), intent (in) :: x
                      real(8), intent (in) :: y
                      res = merge(y, merge(x, max(x, y), y /= y), x /= x)
                  end function
                  real(4) function fmax44(x, y) result (res)
                      real(4), intent (in) :: x
                      real(4), intent (in) :: y
                      res = merge(y, merge(x, max(x, y), y /= y), x /= x)
                  end function
                  real(8) function fmax84(x, y) result(res)
                      real(8), intent (in) :: x
                      real(4), intent (in) :: y
                      res = merge(dble(y), merge(x, max(x, dble(y)), y /= y), x /= x)
                  end function
                  real(8) function fmax48(x, y) result(res)
                      real(4), intent (in) :: x
                      real(8), intent (in) :: y
                      res = merge(y, merge(dble(x), max(dble(x), y), y /= y), x /= x)
                  end function
                  real(8) function fmin88(x, y) result (res)
                      real(8), intent (in) :: x
                      real(8), intent (in) :: y
                      res = merge(y, merge(x, min(x, y), y /= y), x /= x)
                  end function
                  real(4) function fmin44(x, y) result (res)
                      real(4), intent (in) :: x
                      real(4), intent (in) :: y
                      res = merge(y, merge(x, min(x, y), y /= y), x /= x)
                  end function
                  real(8) function fmin84(x, y) result(res)
                      real(8), intent (in) :: x
                      real(4), intent (in) :: y
                      res = merge(dble(y), merge(x, min(x, dble(y)), y /= y), x /= x)
                  end function
                  real(8) function fmin48(x, y) result(res)
                      real(4), intent (in) :: x
                      real(8), intent (in) :: y
                      res = merge(y, merge(dble(x), min(dble(x), y), y /= y), x /= x)
                  end function
              end module
              
              real(8) function code(eh, ew, t)
              use fmin_fmax_functions
                  real(8), intent (in) :: eh
                  real(8), intent (in) :: ew
                  real(8), intent (in) :: t
                  code = abs((ew * sin(t)))
              end function
              
              public static double code(double eh, double ew, double t) {
              	return Math.abs((ew * Math.sin(t)));
              }
              
              def code(eh, ew, t):
              	return math.fabs((ew * math.sin(t)))
              
              function code(eh, ew, t)
              	return abs(Float64(ew * sin(t)))
              end
              
              function tmp = code(eh, ew, t)
              	tmp = abs((ew * sin(t)));
              end
              
              code[eh_, ew_, t_] := N[Abs[N[(ew * N[Sin[t], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]
              
              \begin{array}{l}
              
              \\
              \left|ew \cdot \sin t\right|
              \end{array}
              
              Derivation
              1. Initial program 99.8%

                \[\left|\left(ew \cdot \sin t\right) \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
              2. Add Preprocessing
              3. Step-by-step derivation
                1. lift-cos.f64N/A

                  \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \color{blue}{\cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
                2. lift-atan.f64N/A

                  \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \cos \color{blue}{\tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
                3. lift-/.f64N/A

                  \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \cos \tan^{-1} \color{blue}{\left(\frac{\frac{eh}{ew}}{\tan t}\right)} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
                4. lift-/.f64N/A

                  \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \cos \tan^{-1} \left(\frac{\color{blue}{\frac{eh}{ew}}}{\tan t}\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
                5. lift-tan.f64N/A

                  \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\color{blue}{\tan t}}\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
                6. cos-atanN/A

                  \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \color{blue}{\frac{1}{\sqrt{1 + \frac{\frac{eh}{ew}}{\tan t} \cdot \frac{\frac{eh}{ew}}{\tan t}}}} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
                7. lower-/.f64N/A

                  \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \color{blue}{\frac{1}{\sqrt{1 + \frac{\frac{eh}{ew}}{\tan t} \cdot \frac{\frac{eh}{ew}}{\tan t}}}} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
                8. lower-sqrt.f64N/A

                  \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \frac{1}{\color{blue}{\sqrt{1 + \frac{\frac{eh}{ew}}{\tan t} \cdot \frac{\frac{eh}{ew}}{\tan t}}}} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
                9. lower-+.f64N/A

                  \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \frac{1}{\sqrt{\color{blue}{1 + \frac{\frac{eh}{ew}}{\tan t} \cdot \frac{\frac{eh}{ew}}{\tan t}}}} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
                10. pow2N/A

                  \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \frac{1}{\sqrt{1 + \color{blue}{{\left(\frac{\frac{eh}{ew}}{\tan t}\right)}^{2}}}} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
                11. lower-pow.f64N/A

                  \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \frac{1}{\sqrt{1 + \color{blue}{{\left(\frac{\frac{eh}{ew}}{\tan t}\right)}^{2}}}} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
                12. lift-/.f64N/A

                  \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \frac{1}{\sqrt{1 + {\left(\frac{\color{blue}{\frac{eh}{ew}}}{\tan t}\right)}^{2}}} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
                13. lift-tan.f64N/A

                  \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \frac{1}{\sqrt{1 + {\left(\frac{\frac{eh}{ew}}{\color{blue}{\tan t}}\right)}^{2}}} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
                14. lift-/.f6499.8

                  \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \frac{1}{\sqrt{1 + {\color{blue}{\left(\frac{\frac{eh}{ew}}{\tan t}\right)}}^{2}}} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
              4. Applied rewrites99.8%

                \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \color{blue}{\frac{1}{\sqrt{1 + {\left(\frac{\frac{eh}{ew}}{\tan t}\right)}^{2}}}} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
              5. Taylor expanded in eh around 0

                \[\leadsto \left|\color{blue}{ew \cdot \sin t}\right| \]
              6. Step-by-step derivation
                1. unpow2N/A

                  \[\leadsto \left|ew \cdot \sin t\right| \]
                2. cos-atanN/A

                  \[\leadsto \left|ew \cdot \sin t\right| \]
                3. lift-sin.f64N/A

                  \[\leadsto \left|ew \cdot \sin t\right| \]
                4. lift-*.f6441.0

                  \[\leadsto \left|ew \cdot \color{blue}{\sin t}\right| \]
              7. Applied rewrites41.0%

                \[\leadsto \left|\color{blue}{ew \cdot \sin t}\right| \]
              8. Add Preprocessing

              Alternative 12: 18.8% accurate, 108.8× speedup?

              \[\begin{array}{l} \\ \left|ew \cdot t\right| \end{array} \]
              (FPCore (eh ew t) :precision binary64 (fabs (* ew t)))
              double code(double eh, double ew, double t) {
              	return fabs((ew * t));
              }
              
              module fmin_fmax_functions
                  implicit none
                  private
                  public fmax
                  public fmin
              
                  interface fmax
                      module procedure fmax88
                      module procedure fmax44
                      module procedure fmax84
                      module procedure fmax48
                  end interface
                  interface fmin
                      module procedure fmin88
                      module procedure fmin44
                      module procedure fmin84
                      module procedure fmin48
                  end interface
              contains
                  real(8) function fmax88(x, y) result (res)
                      real(8), intent (in) :: x
                      real(8), intent (in) :: y
                      res = merge(y, merge(x, max(x, y), y /= y), x /= x)
                  end function
                  real(4) function fmax44(x, y) result (res)
                      real(4), intent (in) :: x
                      real(4), intent (in) :: y
                      res = merge(y, merge(x, max(x, y), y /= y), x /= x)
                  end function
                  real(8) function fmax84(x, y) result(res)
                      real(8), intent (in) :: x
                      real(4), intent (in) :: y
                      res = merge(dble(y), merge(x, max(x, dble(y)), y /= y), x /= x)
                  end function
                  real(8) function fmax48(x, y) result(res)
                      real(4), intent (in) :: x
                      real(8), intent (in) :: y
                      res = merge(y, merge(dble(x), max(dble(x), y), y /= y), x /= x)
                  end function
                  real(8) function fmin88(x, y) result (res)
                      real(8), intent (in) :: x
                      real(8), intent (in) :: y
                      res = merge(y, merge(x, min(x, y), y /= y), x /= x)
                  end function
                  real(4) function fmin44(x, y) result (res)
                      real(4), intent (in) :: x
                      real(4), intent (in) :: y
                      res = merge(y, merge(x, min(x, y), y /= y), x /= x)
                  end function
                  real(8) function fmin84(x, y) result(res)
                      real(8), intent (in) :: x
                      real(4), intent (in) :: y
                      res = merge(dble(y), merge(x, min(x, dble(y)), y /= y), x /= x)
                  end function
                  real(8) function fmin48(x, y) result(res)
                      real(4), intent (in) :: x
                      real(8), intent (in) :: y
                      res = merge(y, merge(dble(x), min(dble(x), y), y /= y), x /= x)
                  end function
              end module
              
              real(8) function code(eh, ew, t)
              use fmin_fmax_functions
                  real(8), intent (in) :: eh
                  real(8), intent (in) :: ew
                  real(8), intent (in) :: t
                  code = abs((ew * t))
              end function
              
              public static double code(double eh, double ew, double t) {
              	return Math.abs((ew * t));
              }
              
              def code(eh, ew, t):
              	return math.fabs((ew * t))
              
              function code(eh, ew, t)
              	return abs(Float64(ew * t))
              end
              
              function tmp = code(eh, ew, t)
              	tmp = abs((ew * t));
              end
              
              code[eh_, ew_, t_] := N[Abs[N[(ew * t), $MachinePrecision]], $MachinePrecision]
              
              \begin{array}{l}
              
              \\
              \left|ew \cdot t\right|
              \end{array}
              
              Derivation
              1. Initial program 99.8%

                \[\left|\left(ew \cdot \sin t\right) \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
              2. Add Preprocessing
              3. Step-by-step derivation
                1. lift-cos.f64N/A

                  \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \color{blue}{\cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
                2. lift-atan.f64N/A

                  \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \cos \color{blue}{\tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
                3. lift-/.f64N/A

                  \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \cos \tan^{-1} \color{blue}{\left(\frac{\frac{eh}{ew}}{\tan t}\right)} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
                4. lift-/.f64N/A

                  \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \cos \tan^{-1} \left(\frac{\color{blue}{\frac{eh}{ew}}}{\tan t}\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
                5. lift-tan.f64N/A

                  \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\color{blue}{\tan t}}\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
                6. cos-atanN/A

                  \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \color{blue}{\frac{1}{\sqrt{1 + \frac{\frac{eh}{ew}}{\tan t} \cdot \frac{\frac{eh}{ew}}{\tan t}}}} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
                7. lower-/.f64N/A

                  \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \color{blue}{\frac{1}{\sqrt{1 + \frac{\frac{eh}{ew}}{\tan t} \cdot \frac{\frac{eh}{ew}}{\tan t}}}} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
                8. lower-sqrt.f64N/A

                  \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \frac{1}{\color{blue}{\sqrt{1 + \frac{\frac{eh}{ew}}{\tan t} \cdot \frac{\frac{eh}{ew}}{\tan t}}}} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
                9. lower-+.f64N/A

                  \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \frac{1}{\sqrt{\color{blue}{1 + \frac{\frac{eh}{ew}}{\tan t} \cdot \frac{\frac{eh}{ew}}{\tan t}}}} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
                10. pow2N/A

                  \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \frac{1}{\sqrt{1 + \color{blue}{{\left(\frac{\frac{eh}{ew}}{\tan t}\right)}^{2}}}} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
                11. lower-pow.f64N/A

                  \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \frac{1}{\sqrt{1 + \color{blue}{{\left(\frac{\frac{eh}{ew}}{\tan t}\right)}^{2}}}} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
                12. lift-/.f64N/A

                  \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \frac{1}{\sqrt{1 + {\left(\frac{\color{blue}{\frac{eh}{ew}}}{\tan t}\right)}^{2}}} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
                13. lift-tan.f64N/A

                  \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \frac{1}{\sqrt{1 + {\left(\frac{\frac{eh}{ew}}{\color{blue}{\tan t}}\right)}^{2}}} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
                14. lift-/.f6499.8

                  \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \frac{1}{\sqrt{1 + {\color{blue}{\left(\frac{\frac{eh}{ew}}{\tan t}\right)}}^{2}}} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
              4. Applied rewrites99.8%

                \[\leadsto \left|\left(ew \cdot \sin t\right) \cdot \color{blue}{\frac{1}{\sqrt{1 + {\left(\frac{\frac{eh}{ew}}{\tan t}\right)}^{2}}}} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right| \]
              5. Taylor expanded in eh around 0

                \[\leadsto \left|\color{blue}{ew \cdot \sin t}\right| \]
              6. Step-by-step derivation
                1. unpow2N/A

                  \[\leadsto \left|ew \cdot \sin t\right| \]
                2. cos-atanN/A

                  \[\leadsto \left|ew \cdot \sin t\right| \]
                3. lift-sin.f64N/A

                  \[\leadsto \left|ew \cdot \sin t\right| \]
                4. lift-*.f6441.0

                  \[\leadsto \left|ew \cdot \color{blue}{\sin t}\right| \]
              7. Applied rewrites41.0%

                \[\leadsto \left|\color{blue}{ew \cdot \sin t}\right| \]
              8. Taylor expanded in t around 0

                \[\leadsto \left|ew \cdot t\right| \]
              9. Step-by-step derivation
                1. Applied rewrites18.3%

                  \[\leadsto \left|ew \cdot t\right| \]
                2. Add Preprocessing

                Reproduce

                ?
                herbie shell --seed 2025082 
                (FPCore (eh ew t)
                  :name "Example from Robby"
                  :precision binary64
                  (fabs (+ (* (* ew (sin t)) (cos (atan (/ (/ eh ew) (tan t))))) (* (* eh (cos t)) (sin (atan (/ (/ eh ew) (tan t))))))))