
(FPCore (r a b) :precision binary64 (/ (* r (sin b)) (cos (+ a b))))
double code(double r, double a, double b) {
return (r * sin(b)) / cos((a + b));
}
real(8) function code(r, a, b)
real(8), intent (in) :: r
real(8), intent (in) :: a
real(8), intent (in) :: b
code = (r * sin(b)) / cos((a + b))
end function
public static double code(double r, double a, double b) {
return (r * Math.sin(b)) / Math.cos((a + b));
}
def code(r, a, b): return (r * math.sin(b)) / math.cos((a + b))
function code(r, a, b) return Float64(Float64(r * sin(b)) / cos(Float64(a + b))) end
function tmp = code(r, a, b) tmp = (r * sin(b)) / cos((a + b)); end
code[r_, a_, b_] := N[(N[(r * N[Sin[b], $MachinePrecision]), $MachinePrecision] / N[Cos[N[(a + b), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{r \cdot \sin b}{\cos \left(a + b\right)}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 23 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (r a b) :precision binary64 (/ (* r (sin b)) (cos (+ a b))))
double code(double r, double a, double b) {
return (r * sin(b)) / cos((a + b));
}
real(8) function code(r, a, b)
real(8), intent (in) :: r
real(8), intent (in) :: a
real(8), intent (in) :: b
code = (r * sin(b)) / cos((a + b))
end function
public static double code(double r, double a, double b) {
return (r * Math.sin(b)) / Math.cos((a + b));
}
def code(r, a, b): return (r * math.sin(b)) / math.cos((a + b))
function code(r, a, b) return Float64(Float64(r * sin(b)) / cos(Float64(a + b))) end
function tmp = code(r, a, b) tmp = (r * sin(b)) / cos((a + b)); end
code[r_, a_, b_] := N[(N[(r * N[Sin[b], $MachinePrecision]), $MachinePrecision] / N[Cos[N[(a + b), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{r \cdot \sin b}{\cos \left(a + b\right)}
\end{array}
(FPCore (r a b) :precision binary64 (/ (* r (sin b)) (fma (- 0.0 (sin b)) (sin a) (* (cos b) (cos a)))))
double code(double r, double a, double b) {
return (r * sin(b)) / fma((0.0 - sin(b)), sin(a), (cos(b) * cos(a)));
}
function code(r, a, b) return Float64(Float64(r * sin(b)) / fma(Float64(0.0 - sin(b)), sin(a), Float64(cos(b) * cos(a)))) end
code[r_, a_, b_] := N[(N[(r * N[Sin[b], $MachinePrecision]), $MachinePrecision] / N[(N[(0.0 - N[Sin[b], $MachinePrecision]), $MachinePrecision] * N[Sin[a], $MachinePrecision] + N[(N[Cos[b], $MachinePrecision] * N[Cos[a], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{r \cdot \sin b}{\mathsf{fma}\left(0 - \sin b, \sin a, \cos b \cdot \cos a\right)}
\end{array}
Initial program 78.5%
/-lowering-/.f64N/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
cos-lowering-cos.f64N/A
+-commutativeN/A
+-lowering-+.f6478.5%
Simplified78.5%
cos-sumN/A
sub-negN/A
+-commutativeN/A
distribute-lft-neg-inN/A
fma-defineN/A
fma-lowering-fma.f64N/A
neg-sub0N/A
--lowering--.f64N/A
sin-lowering-sin.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
cos-lowering-cos.f6499.6%
Applied egg-rr99.6%
remove-double-divN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
sin-lowering-sin.f6499.6%
Applied egg-rr99.6%
sub0-negN/A
remove-double-divN/A
neg-lowering-neg.f64N/A
sin-lowering-sin.f6499.6%
Applied egg-rr99.6%
Final simplification99.6%
(FPCore (r a b) :precision binary64 (/ (* r (sin b)) (- (* (cos b) (cos a)) (* (sin b) (sin a)))))
double code(double r, double a, double b) {
return (r * sin(b)) / ((cos(b) * cos(a)) - (sin(b) * sin(a)));
}
real(8) function code(r, a, b)
real(8), intent (in) :: r
real(8), intent (in) :: a
real(8), intent (in) :: b
code = (r * sin(b)) / ((cos(b) * cos(a)) - (sin(b) * sin(a)))
end function
public static double code(double r, double a, double b) {
return (r * Math.sin(b)) / ((Math.cos(b) * Math.cos(a)) - (Math.sin(b) * Math.sin(a)));
}
def code(r, a, b): return (r * math.sin(b)) / ((math.cos(b) * math.cos(a)) - (math.sin(b) * math.sin(a)))
function code(r, a, b) return Float64(Float64(r * sin(b)) / Float64(Float64(cos(b) * cos(a)) - Float64(sin(b) * sin(a)))) end
function tmp = code(r, a, b) tmp = (r * sin(b)) / ((cos(b) * cos(a)) - (sin(b) * sin(a))); end
code[r_, a_, b_] := N[(N[(r * N[Sin[b], $MachinePrecision]), $MachinePrecision] / N[(N[(N[Cos[b], $MachinePrecision] * N[Cos[a], $MachinePrecision]), $MachinePrecision] - N[(N[Sin[b], $MachinePrecision] * N[Sin[a], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{r \cdot \sin b}{\cos b \cdot \cos a - \sin b \cdot \sin a}
\end{array}
Initial program 78.5%
/-lowering-/.f64N/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
cos-lowering-cos.f64N/A
+-commutativeN/A
+-lowering-+.f6478.5%
Simplified78.5%
cos-sumN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
cos-lowering-cos.f64N/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
sin-lowering-sin.f6499.6%
Applied egg-rr99.6%
(FPCore (r a b) :precision binary64 (* (sin b) (/ r (- (* (cos b) (cos a)) (* (sin b) (sin a))))))
double code(double r, double a, double b) {
return sin(b) * (r / ((cos(b) * cos(a)) - (sin(b) * sin(a))));
}
real(8) function code(r, a, b)
real(8), intent (in) :: r
real(8), intent (in) :: a
real(8), intent (in) :: b
code = sin(b) * (r / ((cos(b) * cos(a)) - (sin(b) * sin(a))))
end function
public static double code(double r, double a, double b) {
return Math.sin(b) * (r / ((Math.cos(b) * Math.cos(a)) - (Math.sin(b) * Math.sin(a))));
}
def code(r, a, b): return math.sin(b) * (r / ((math.cos(b) * math.cos(a)) - (math.sin(b) * math.sin(a))))
function code(r, a, b) return Float64(sin(b) * Float64(r / Float64(Float64(cos(b) * cos(a)) - Float64(sin(b) * sin(a))))) end
function tmp = code(r, a, b) tmp = sin(b) * (r / ((cos(b) * cos(a)) - (sin(b) * sin(a)))); end
code[r_, a_, b_] := N[(N[Sin[b], $MachinePrecision] * N[(r / N[(N[(N[Cos[b], $MachinePrecision] * N[Cos[a], $MachinePrecision]), $MachinePrecision] - N[(N[Sin[b], $MachinePrecision] * N[Sin[a], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\sin b \cdot \frac{r}{\cos b \cdot \cos a - \sin b \cdot \sin a}
\end{array}
Initial program 78.5%
/-lowering-/.f64N/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
cos-lowering-cos.f64N/A
+-commutativeN/A
+-lowering-+.f6478.5%
Simplified78.5%
associate-*l/N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
cos-lowering-cos.f64N/A
+-lowering-+.f64N/A
sin-lowering-sin.f6478.5%
Applied egg-rr78.5%
cos-sumN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
cos-lowering-cos.f64N/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
sin-lowering-sin.f6499.5%
Applied egg-rr99.5%
Final simplification99.5%
(FPCore (r a b)
:precision binary64
(let* ((t_0 (* r (sin b))))
(if (<= a -6.8e-5)
(* r (/ (sin b) (cos a)))
(if (<= a 2.3e+15) (/ t_0 (cos b)) (/ t_0 (cos a))))))
double code(double r, double a, double b) {
double t_0 = r * sin(b);
double tmp;
if (a <= -6.8e-5) {
tmp = r * (sin(b) / cos(a));
} else if (a <= 2.3e+15) {
tmp = t_0 / cos(b);
} else {
tmp = t_0 / cos(a);
}
return tmp;
}
real(8) function code(r, a, b)
real(8), intent (in) :: r
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: t_0
real(8) :: tmp
t_0 = r * sin(b)
if (a <= (-6.8d-5)) then
tmp = r * (sin(b) / cos(a))
else if (a <= 2.3d+15) then
tmp = t_0 / cos(b)
else
tmp = t_0 / cos(a)
end if
code = tmp
end function
public static double code(double r, double a, double b) {
double t_0 = r * Math.sin(b);
double tmp;
if (a <= -6.8e-5) {
tmp = r * (Math.sin(b) / Math.cos(a));
} else if (a <= 2.3e+15) {
tmp = t_0 / Math.cos(b);
} else {
tmp = t_0 / Math.cos(a);
}
return tmp;
}
def code(r, a, b): t_0 = r * math.sin(b) tmp = 0 if a <= -6.8e-5: tmp = r * (math.sin(b) / math.cos(a)) elif a <= 2.3e+15: tmp = t_0 / math.cos(b) else: tmp = t_0 / math.cos(a) return tmp
function code(r, a, b) t_0 = Float64(r * sin(b)) tmp = 0.0 if (a <= -6.8e-5) tmp = Float64(r * Float64(sin(b) / cos(a))); elseif (a <= 2.3e+15) tmp = Float64(t_0 / cos(b)); else tmp = Float64(t_0 / cos(a)); end return tmp end
function tmp_2 = code(r, a, b) t_0 = r * sin(b); tmp = 0.0; if (a <= -6.8e-5) tmp = r * (sin(b) / cos(a)); elseif (a <= 2.3e+15) tmp = t_0 / cos(b); else tmp = t_0 / cos(a); end tmp_2 = tmp; end
code[r_, a_, b_] := Block[{t$95$0 = N[(r * N[Sin[b], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[a, -6.8e-5], N[(r * N[(N[Sin[b], $MachinePrecision] / N[Cos[a], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[a, 2.3e+15], N[(t$95$0 / N[Cos[b], $MachinePrecision]), $MachinePrecision], N[(t$95$0 / N[Cos[a], $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := r \cdot \sin b\\
\mathbf{if}\;a \leq -6.8 \cdot 10^{-5}:\\
\;\;\;\;r \cdot \frac{\sin b}{\cos a}\\
\mathbf{elif}\;a \leq 2.3 \cdot 10^{+15}:\\
\;\;\;\;\frac{t\_0}{\cos b}\\
\mathbf{else}:\\
\;\;\;\;\frac{t\_0}{\cos a}\\
\end{array}
\end{array}
if a < -6.7999999999999999e-5Initial program 59.8%
/-lowering-/.f64N/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
cos-lowering-cos.f64N/A
+-commutativeN/A
+-lowering-+.f6459.8%
Simplified59.8%
associate-*l/N/A
associate-/r/N/A
*-lft-identityN/A
div-invN/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
cos-lowering-cos.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
sin-lowering-sin.f6459.7%
Applied egg-rr59.7%
Taylor expanded in b around 0
cos-lowering-cos.f6459.0%
Simplified59.0%
clear-numN/A
un-div-invN/A
associate-/r/N/A
div-invN/A
associate-*l/N/A
remove-double-divN/A
*-lowering-*.f64N/A
un-div-invN/A
/-lowering-/.f64N/A
sin-lowering-sin.f64N/A
cos-lowering-cos.f6459.1%
Applied egg-rr59.1%
if -6.7999999999999999e-5 < a < 2.3e15Initial program 98.0%
/-lowering-/.f64N/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
cos-lowering-cos.f64N/A
+-commutativeN/A
+-lowering-+.f6498.0%
Simplified98.0%
Taylor expanded in a around 0
/-lowering-/.f64N/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
cos-lowering-cos.f6498.0%
Simplified98.0%
if 2.3e15 < a Initial program 57.0%
/-lowering-/.f64N/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
cos-lowering-cos.f64N/A
+-commutativeN/A
+-lowering-+.f6457.0%
Simplified57.0%
Taylor expanded in b around 0
cos-lowering-cos.f6456.6%
Simplified56.6%
Final simplification78.3%
(FPCore (r a b) :precision binary64 (if (<= a -0.000235) (* r (/ (sin b) (cos a))) (if (<= a 2.3e+15) (* (sin b) (/ r (cos b))) (/ (* r (sin b)) (cos a)))))
double code(double r, double a, double b) {
double tmp;
if (a <= -0.000235) {
tmp = r * (sin(b) / cos(a));
} else if (a <= 2.3e+15) {
tmp = sin(b) * (r / cos(b));
} else {
tmp = (r * sin(b)) / cos(a);
}
return tmp;
}
real(8) function code(r, a, b)
real(8), intent (in) :: r
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: tmp
if (a <= (-0.000235d0)) then
tmp = r * (sin(b) / cos(a))
else if (a <= 2.3d+15) then
tmp = sin(b) * (r / cos(b))
else
tmp = (r * sin(b)) / cos(a)
end if
code = tmp
end function
public static double code(double r, double a, double b) {
double tmp;
if (a <= -0.000235) {
tmp = r * (Math.sin(b) / Math.cos(a));
} else if (a <= 2.3e+15) {
tmp = Math.sin(b) * (r / Math.cos(b));
} else {
tmp = (r * Math.sin(b)) / Math.cos(a);
}
return tmp;
}
def code(r, a, b): tmp = 0 if a <= -0.000235: tmp = r * (math.sin(b) / math.cos(a)) elif a <= 2.3e+15: tmp = math.sin(b) * (r / math.cos(b)) else: tmp = (r * math.sin(b)) / math.cos(a) return tmp
function code(r, a, b) tmp = 0.0 if (a <= -0.000235) tmp = Float64(r * Float64(sin(b) / cos(a))); elseif (a <= 2.3e+15) tmp = Float64(sin(b) * Float64(r / cos(b))); else tmp = Float64(Float64(r * sin(b)) / cos(a)); end return tmp end
function tmp_2 = code(r, a, b) tmp = 0.0; if (a <= -0.000235) tmp = r * (sin(b) / cos(a)); elseif (a <= 2.3e+15) tmp = sin(b) * (r / cos(b)); else tmp = (r * sin(b)) / cos(a); end tmp_2 = tmp; end
code[r_, a_, b_] := If[LessEqual[a, -0.000235], N[(r * N[(N[Sin[b], $MachinePrecision] / N[Cos[a], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[a, 2.3e+15], N[(N[Sin[b], $MachinePrecision] * N[(r / N[Cos[b], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(r * N[Sin[b], $MachinePrecision]), $MachinePrecision] / N[Cos[a], $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq -0.000235:\\
\;\;\;\;r \cdot \frac{\sin b}{\cos a}\\
\mathbf{elif}\;a \leq 2.3 \cdot 10^{+15}:\\
\;\;\;\;\sin b \cdot \frac{r}{\cos b}\\
\mathbf{else}:\\
\;\;\;\;\frac{r \cdot \sin b}{\cos a}\\
\end{array}
\end{array}
if a < -2.34999999999999993e-4Initial program 59.8%
/-lowering-/.f64N/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
cos-lowering-cos.f64N/A
+-commutativeN/A
+-lowering-+.f6459.8%
Simplified59.8%
associate-*l/N/A
associate-/r/N/A
*-lft-identityN/A
div-invN/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
cos-lowering-cos.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
sin-lowering-sin.f6459.7%
Applied egg-rr59.7%
Taylor expanded in b around 0
cos-lowering-cos.f6459.0%
Simplified59.0%
clear-numN/A
un-div-invN/A
associate-/r/N/A
div-invN/A
associate-*l/N/A
remove-double-divN/A
*-lowering-*.f64N/A
un-div-invN/A
/-lowering-/.f64N/A
sin-lowering-sin.f64N/A
cos-lowering-cos.f6459.1%
Applied egg-rr59.1%
if -2.34999999999999993e-4 < a < 2.3e15Initial program 98.0%
/-lowering-/.f64N/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
cos-lowering-cos.f64N/A
+-commutativeN/A
+-lowering-+.f6498.0%
Simplified98.0%
associate-*l/N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
cos-lowering-cos.f64N/A
+-lowering-+.f64N/A
sin-lowering-sin.f6498.0%
Applied egg-rr98.0%
Taylor expanded in a around 0
/-lowering-/.f64N/A
cos-lowering-cos.f6498.0%
Simplified98.0%
if 2.3e15 < a Initial program 57.0%
/-lowering-/.f64N/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
cos-lowering-cos.f64N/A
+-commutativeN/A
+-lowering-+.f6457.0%
Simplified57.0%
Taylor expanded in b around 0
cos-lowering-cos.f6456.6%
Simplified56.6%
Final simplification78.2%
(FPCore (r a b) :precision binary64 (if (<= a -0.000155) (* r (/ (sin b) (cos a))) (if (<= a 2.3e+15) (* (sin b) (/ r (cos b))) (/ r (/ (cos a) (sin b))))))
double code(double r, double a, double b) {
double tmp;
if (a <= -0.000155) {
tmp = r * (sin(b) / cos(a));
} else if (a <= 2.3e+15) {
tmp = sin(b) * (r / cos(b));
} else {
tmp = r / (cos(a) / sin(b));
}
return tmp;
}
real(8) function code(r, a, b)
real(8), intent (in) :: r
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: tmp
if (a <= (-0.000155d0)) then
tmp = r * (sin(b) / cos(a))
else if (a <= 2.3d+15) then
tmp = sin(b) * (r / cos(b))
else
tmp = r / (cos(a) / sin(b))
end if
code = tmp
end function
public static double code(double r, double a, double b) {
double tmp;
if (a <= -0.000155) {
tmp = r * (Math.sin(b) / Math.cos(a));
} else if (a <= 2.3e+15) {
tmp = Math.sin(b) * (r / Math.cos(b));
} else {
tmp = r / (Math.cos(a) / Math.sin(b));
}
return tmp;
}
def code(r, a, b): tmp = 0 if a <= -0.000155: tmp = r * (math.sin(b) / math.cos(a)) elif a <= 2.3e+15: tmp = math.sin(b) * (r / math.cos(b)) else: tmp = r / (math.cos(a) / math.sin(b)) return tmp
function code(r, a, b) tmp = 0.0 if (a <= -0.000155) tmp = Float64(r * Float64(sin(b) / cos(a))); elseif (a <= 2.3e+15) tmp = Float64(sin(b) * Float64(r / cos(b))); else tmp = Float64(r / Float64(cos(a) / sin(b))); end return tmp end
function tmp_2 = code(r, a, b) tmp = 0.0; if (a <= -0.000155) tmp = r * (sin(b) / cos(a)); elseif (a <= 2.3e+15) tmp = sin(b) * (r / cos(b)); else tmp = r / (cos(a) / sin(b)); end tmp_2 = tmp; end
code[r_, a_, b_] := If[LessEqual[a, -0.000155], N[(r * N[(N[Sin[b], $MachinePrecision] / N[Cos[a], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[a, 2.3e+15], N[(N[Sin[b], $MachinePrecision] * N[(r / N[Cos[b], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(r / N[(N[Cos[a], $MachinePrecision] / N[Sin[b], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq -0.000155:\\
\;\;\;\;r \cdot \frac{\sin b}{\cos a}\\
\mathbf{elif}\;a \leq 2.3 \cdot 10^{+15}:\\
\;\;\;\;\sin b \cdot \frac{r}{\cos b}\\
\mathbf{else}:\\
\;\;\;\;\frac{r}{\frac{\cos a}{\sin b}}\\
\end{array}
\end{array}
if a < -1.55e-4Initial program 59.8%
/-lowering-/.f64N/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
cos-lowering-cos.f64N/A
+-commutativeN/A
+-lowering-+.f6459.8%
Simplified59.8%
associate-*l/N/A
associate-/r/N/A
*-lft-identityN/A
div-invN/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
cos-lowering-cos.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
sin-lowering-sin.f6459.7%
Applied egg-rr59.7%
Taylor expanded in b around 0
cos-lowering-cos.f6459.0%
Simplified59.0%
clear-numN/A
un-div-invN/A
associate-/r/N/A
div-invN/A
associate-*l/N/A
remove-double-divN/A
*-lowering-*.f64N/A
un-div-invN/A
/-lowering-/.f64N/A
sin-lowering-sin.f64N/A
cos-lowering-cos.f6459.1%
Applied egg-rr59.1%
if -1.55e-4 < a < 2.3e15Initial program 98.0%
/-lowering-/.f64N/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
cos-lowering-cos.f64N/A
+-commutativeN/A
+-lowering-+.f6498.0%
Simplified98.0%
associate-*l/N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
cos-lowering-cos.f64N/A
+-lowering-+.f64N/A
sin-lowering-sin.f6498.0%
Applied egg-rr98.0%
Taylor expanded in a around 0
/-lowering-/.f64N/A
cos-lowering-cos.f6498.0%
Simplified98.0%
if 2.3e15 < a Initial program 57.0%
/-lowering-/.f64N/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
cos-lowering-cos.f64N/A
+-commutativeN/A
+-lowering-+.f6457.0%
Simplified57.0%
associate-*l/N/A
associate-/r/N/A
*-lft-identityN/A
div-invN/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
cos-lowering-cos.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
sin-lowering-sin.f6456.8%
Applied egg-rr56.8%
Taylor expanded in b around 0
cos-lowering-cos.f6456.5%
Simplified56.5%
frac-timesN/A
*-lft-identityN/A
/-lowering-/.f64N/A
un-div-invN/A
/-lowering-/.f64N/A
cos-lowering-cos.f64N/A
sin-lowering-sin.f6456.6%
Applied egg-rr56.6%
Final simplification78.2%
(FPCore (r a b) :precision binary64 (let* ((t_0 (* r (/ (sin b) (cos a))))) (if (<= a -7.5e-6) t_0 (if (<= a 2.3e+15) (* (sin b) (/ r (cos b))) t_0))))
double code(double r, double a, double b) {
double t_0 = r * (sin(b) / cos(a));
double tmp;
if (a <= -7.5e-6) {
tmp = t_0;
} else if (a <= 2.3e+15) {
tmp = sin(b) * (r / cos(b));
} else {
tmp = t_0;
}
return tmp;
}
real(8) function code(r, a, b)
real(8), intent (in) :: r
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: t_0
real(8) :: tmp
t_0 = r * (sin(b) / cos(a))
if (a <= (-7.5d-6)) then
tmp = t_0
else if (a <= 2.3d+15) then
tmp = sin(b) * (r / cos(b))
else
tmp = t_0
end if
code = tmp
end function
public static double code(double r, double a, double b) {
double t_0 = r * (Math.sin(b) / Math.cos(a));
double tmp;
if (a <= -7.5e-6) {
tmp = t_0;
} else if (a <= 2.3e+15) {
tmp = Math.sin(b) * (r / Math.cos(b));
} else {
tmp = t_0;
}
return tmp;
}
def code(r, a, b): t_0 = r * (math.sin(b) / math.cos(a)) tmp = 0 if a <= -7.5e-6: tmp = t_0 elif a <= 2.3e+15: tmp = math.sin(b) * (r / math.cos(b)) else: tmp = t_0 return tmp
function code(r, a, b) t_0 = Float64(r * Float64(sin(b) / cos(a))) tmp = 0.0 if (a <= -7.5e-6) tmp = t_0; elseif (a <= 2.3e+15) tmp = Float64(sin(b) * Float64(r / cos(b))); else tmp = t_0; end return tmp end
function tmp_2 = code(r, a, b) t_0 = r * (sin(b) / cos(a)); tmp = 0.0; if (a <= -7.5e-6) tmp = t_0; elseif (a <= 2.3e+15) tmp = sin(b) * (r / cos(b)); else tmp = t_0; end tmp_2 = tmp; end
code[r_, a_, b_] := Block[{t$95$0 = N[(r * N[(N[Sin[b], $MachinePrecision] / N[Cos[a], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[a, -7.5e-6], t$95$0, If[LessEqual[a, 2.3e+15], N[(N[Sin[b], $MachinePrecision] * N[(r / N[Cos[b], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$0]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := r \cdot \frac{\sin b}{\cos a}\\
\mathbf{if}\;a \leq -7.5 \cdot 10^{-6}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;a \leq 2.3 \cdot 10^{+15}:\\
\;\;\;\;\sin b \cdot \frac{r}{\cos b}\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if a < -7.50000000000000019e-6 or 2.3e15 < a Initial program 58.4%
/-lowering-/.f64N/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
cos-lowering-cos.f64N/A
+-commutativeN/A
+-lowering-+.f6458.4%
Simplified58.4%
associate-*l/N/A
associate-/r/N/A
*-lft-identityN/A
div-invN/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
cos-lowering-cos.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
sin-lowering-sin.f6458.3%
Applied egg-rr58.3%
Taylor expanded in b around 0
cos-lowering-cos.f6457.8%
Simplified57.8%
clear-numN/A
un-div-invN/A
associate-/r/N/A
div-invN/A
associate-*l/N/A
remove-double-divN/A
*-lowering-*.f64N/A
un-div-invN/A
/-lowering-/.f64N/A
sin-lowering-sin.f64N/A
cos-lowering-cos.f6457.8%
Applied egg-rr57.8%
if -7.50000000000000019e-6 < a < 2.3e15Initial program 98.0%
/-lowering-/.f64N/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
cos-lowering-cos.f64N/A
+-commutativeN/A
+-lowering-+.f6498.0%
Simplified98.0%
associate-*l/N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
cos-lowering-cos.f64N/A
+-lowering-+.f64N/A
sin-lowering-sin.f6498.0%
Applied egg-rr98.0%
Taylor expanded in a around 0
/-lowering-/.f64N/A
cos-lowering-cos.f6498.0%
Simplified98.0%
Final simplification78.2%
(FPCore (r a b) :precision binary64 (let* ((t_0 (* (sin b) (/ r (cos a))))) (if (<= a -3.8e-5) t_0 (if (<= a 2.3e+15) (* (sin b) (/ r (cos b))) t_0))))
double code(double r, double a, double b) {
double t_0 = sin(b) * (r / cos(a));
double tmp;
if (a <= -3.8e-5) {
tmp = t_0;
} else if (a <= 2.3e+15) {
tmp = sin(b) * (r / cos(b));
} else {
tmp = t_0;
}
return tmp;
}
real(8) function code(r, a, b)
real(8), intent (in) :: r
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: t_0
real(8) :: tmp
t_0 = sin(b) * (r / cos(a))
if (a <= (-3.8d-5)) then
tmp = t_0
else if (a <= 2.3d+15) then
tmp = sin(b) * (r / cos(b))
else
tmp = t_0
end if
code = tmp
end function
public static double code(double r, double a, double b) {
double t_0 = Math.sin(b) * (r / Math.cos(a));
double tmp;
if (a <= -3.8e-5) {
tmp = t_0;
} else if (a <= 2.3e+15) {
tmp = Math.sin(b) * (r / Math.cos(b));
} else {
tmp = t_0;
}
return tmp;
}
def code(r, a, b): t_0 = math.sin(b) * (r / math.cos(a)) tmp = 0 if a <= -3.8e-5: tmp = t_0 elif a <= 2.3e+15: tmp = math.sin(b) * (r / math.cos(b)) else: tmp = t_0 return tmp
function code(r, a, b) t_0 = Float64(sin(b) * Float64(r / cos(a))) tmp = 0.0 if (a <= -3.8e-5) tmp = t_0; elseif (a <= 2.3e+15) tmp = Float64(sin(b) * Float64(r / cos(b))); else tmp = t_0; end return tmp end
function tmp_2 = code(r, a, b) t_0 = sin(b) * (r / cos(a)); tmp = 0.0; if (a <= -3.8e-5) tmp = t_0; elseif (a <= 2.3e+15) tmp = sin(b) * (r / cos(b)); else tmp = t_0; end tmp_2 = tmp; end
code[r_, a_, b_] := Block[{t$95$0 = N[(N[Sin[b], $MachinePrecision] * N[(r / N[Cos[a], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[a, -3.8e-5], t$95$0, If[LessEqual[a, 2.3e+15], N[(N[Sin[b], $MachinePrecision] * N[(r / N[Cos[b], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$0]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sin b \cdot \frac{r}{\cos a}\\
\mathbf{if}\;a \leq -3.8 \cdot 10^{-5}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;a \leq 2.3 \cdot 10^{+15}:\\
\;\;\;\;\sin b \cdot \frac{r}{\cos b}\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if a < -3.8000000000000002e-5 or 2.3e15 < a Initial program 58.4%
/-lowering-/.f64N/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
cos-lowering-cos.f64N/A
+-commutativeN/A
+-lowering-+.f6458.4%
Simplified58.4%
associate-*l/N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
cos-lowering-cos.f64N/A
+-lowering-+.f64N/A
sin-lowering-sin.f6458.3%
Applied egg-rr58.3%
Taylor expanded in b around 0
/-lowering-/.f64N/A
cos-lowering-cos.f6457.8%
Simplified57.8%
if -3.8000000000000002e-5 < a < 2.3e15Initial program 98.0%
/-lowering-/.f64N/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
cos-lowering-cos.f64N/A
+-commutativeN/A
+-lowering-+.f6498.0%
Simplified98.0%
associate-*l/N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
cos-lowering-cos.f64N/A
+-lowering-+.f64N/A
sin-lowering-sin.f6498.0%
Applied egg-rr98.0%
Taylor expanded in a around 0
/-lowering-/.f64N/A
cos-lowering-cos.f6498.0%
Simplified98.0%
Final simplification78.2%
(FPCore (r a b) :precision binary64 (/ (* r (sin b)) (cos (+ b a))))
double code(double r, double a, double b) {
return (r * sin(b)) / cos((b + a));
}
real(8) function code(r, a, b)
real(8), intent (in) :: r
real(8), intent (in) :: a
real(8), intent (in) :: b
code = (r * sin(b)) / cos((b + a))
end function
public static double code(double r, double a, double b) {
return (r * Math.sin(b)) / Math.cos((b + a));
}
def code(r, a, b): return (r * math.sin(b)) / math.cos((b + a))
function code(r, a, b) return Float64(Float64(r * sin(b)) / cos(Float64(b + a))) end
function tmp = code(r, a, b) tmp = (r * sin(b)) / cos((b + a)); end
code[r_, a_, b_] := N[(N[(r * N[Sin[b], $MachinePrecision]), $MachinePrecision] / N[Cos[N[(b + a), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{r \cdot \sin b}{\cos \left(b + a\right)}
\end{array}
Initial program 78.5%
Final simplification78.5%
(FPCore (r a b) :precision binary64 (* r (/ (sin b) (cos (+ b a)))))
double code(double r, double a, double b) {
return r * (sin(b) / cos((b + a)));
}
real(8) function code(r, a, b)
real(8), intent (in) :: r
real(8), intent (in) :: a
real(8), intent (in) :: b
code = r * (sin(b) / cos((b + a)))
end function
public static double code(double r, double a, double b) {
return r * (Math.sin(b) / Math.cos((b + a)));
}
def code(r, a, b): return r * (math.sin(b) / math.cos((b + a)))
function code(r, a, b) return Float64(r * Float64(sin(b) / cos(Float64(b + a)))) end
function tmp = code(r, a, b) tmp = r * (sin(b) / cos((b + a))); end
code[r_, a_, b_] := N[(r * N[(N[Sin[b], $MachinePrecision] / N[Cos[N[(b + a), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
r \cdot \frac{\sin b}{\cos \left(b + a\right)}
\end{array}
Initial program 78.5%
/-lowering-/.f64N/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
cos-lowering-cos.f64N/A
+-commutativeN/A
+-lowering-+.f6478.5%
Simplified78.5%
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
sin-lowering-sin.f64N/A
cos-lowering-cos.f64N/A
+-lowering-+.f6478.5%
Applied egg-rr78.5%
Final simplification78.5%
(FPCore (r a b) :precision binary64 (* (sin b) (/ r (cos (+ b a)))))
double code(double r, double a, double b) {
return sin(b) * (r / cos((b + a)));
}
real(8) function code(r, a, b)
real(8), intent (in) :: r
real(8), intent (in) :: a
real(8), intent (in) :: b
code = sin(b) * (r / cos((b + a)))
end function
public static double code(double r, double a, double b) {
return Math.sin(b) * (r / Math.cos((b + a)));
}
def code(r, a, b): return math.sin(b) * (r / math.cos((b + a)))
function code(r, a, b) return Float64(sin(b) * Float64(r / cos(Float64(b + a)))) end
function tmp = code(r, a, b) tmp = sin(b) * (r / cos((b + a))); end
code[r_, a_, b_] := N[(N[Sin[b], $MachinePrecision] * N[(r / N[Cos[N[(b + a), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\sin b \cdot \frac{r}{\cos \left(b + a\right)}
\end{array}
Initial program 78.5%
/-lowering-/.f64N/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
cos-lowering-cos.f64N/A
+-commutativeN/A
+-lowering-+.f6478.5%
Simplified78.5%
associate-*l/N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
cos-lowering-cos.f64N/A
+-lowering-+.f64N/A
sin-lowering-sin.f6478.5%
Applied egg-rr78.5%
Final simplification78.5%
(FPCore (r a b) :precision binary64 (* (sin b) (/ r (cos a))))
double code(double r, double a, double b) {
return sin(b) * (r / cos(a));
}
real(8) function code(r, a, b)
real(8), intent (in) :: r
real(8), intent (in) :: a
real(8), intent (in) :: b
code = sin(b) * (r / cos(a))
end function
public static double code(double r, double a, double b) {
return Math.sin(b) * (r / Math.cos(a));
}
def code(r, a, b): return math.sin(b) * (r / math.cos(a))
function code(r, a, b) return Float64(sin(b) * Float64(r / cos(a))) end
function tmp = code(r, a, b) tmp = sin(b) * (r / cos(a)); end
code[r_, a_, b_] := N[(N[Sin[b], $MachinePrecision] * N[(r / N[Cos[a], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\sin b \cdot \frac{r}{\cos a}
\end{array}
Initial program 78.5%
/-lowering-/.f64N/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
cos-lowering-cos.f64N/A
+-commutativeN/A
+-lowering-+.f6478.5%
Simplified78.5%
associate-*l/N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
cos-lowering-cos.f64N/A
+-lowering-+.f64N/A
sin-lowering-sin.f6478.5%
Applied egg-rr78.5%
Taylor expanded in b around 0
/-lowering-/.f64N/A
cos-lowering-cos.f6454.0%
Simplified54.0%
Final simplification54.0%
(FPCore (r a b)
:precision binary64
(if (<= b -9e+27)
(* r (sin b))
(if (<= b 3.8e+25)
(/
(*
r
(*
b
(+
1.0
(*
(* b b)
(+
-0.16666666666666666
(*
(* b b)
(+ 0.008333333333333333 (* (* b b) -0.0001984126984126984))))))))
(cos (+ b a)))
(/ r (/ 1.0 (sin b))))))
double code(double r, double a, double b) {
double tmp;
if (b <= -9e+27) {
tmp = r * sin(b);
} else if (b <= 3.8e+25) {
tmp = (r * (b * (1.0 + ((b * b) * (-0.16666666666666666 + ((b * b) * (0.008333333333333333 + ((b * b) * -0.0001984126984126984)))))))) / cos((b + a));
} else {
tmp = r / (1.0 / sin(b));
}
return tmp;
}
real(8) function code(r, a, b)
real(8), intent (in) :: r
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: tmp
if (b <= (-9d+27)) then
tmp = r * sin(b)
else if (b <= 3.8d+25) then
tmp = (r * (b * (1.0d0 + ((b * b) * ((-0.16666666666666666d0) + ((b * b) * (0.008333333333333333d0 + ((b * b) * (-0.0001984126984126984d0))))))))) / cos((b + a))
else
tmp = r / (1.0d0 / sin(b))
end if
code = tmp
end function
public static double code(double r, double a, double b) {
double tmp;
if (b <= -9e+27) {
tmp = r * Math.sin(b);
} else if (b <= 3.8e+25) {
tmp = (r * (b * (1.0 + ((b * b) * (-0.16666666666666666 + ((b * b) * (0.008333333333333333 + ((b * b) * -0.0001984126984126984)))))))) / Math.cos((b + a));
} else {
tmp = r / (1.0 / Math.sin(b));
}
return tmp;
}
def code(r, a, b): tmp = 0 if b <= -9e+27: tmp = r * math.sin(b) elif b <= 3.8e+25: tmp = (r * (b * (1.0 + ((b * b) * (-0.16666666666666666 + ((b * b) * (0.008333333333333333 + ((b * b) * -0.0001984126984126984)))))))) / math.cos((b + a)) else: tmp = r / (1.0 / math.sin(b)) return tmp
function code(r, a, b) tmp = 0.0 if (b <= -9e+27) tmp = Float64(r * sin(b)); elseif (b <= 3.8e+25) tmp = Float64(Float64(r * Float64(b * Float64(1.0 + Float64(Float64(b * b) * Float64(-0.16666666666666666 + Float64(Float64(b * b) * Float64(0.008333333333333333 + Float64(Float64(b * b) * -0.0001984126984126984)))))))) / cos(Float64(b + a))); else tmp = Float64(r / Float64(1.0 / sin(b))); end return tmp end
function tmp_2 = code(r, a, b) tmp = 0.0; if (b <= -9e+27) tmp = r * sin(b); elseif (b <= 3.8e+25) tmp = (r * (b * (1.0 + ((b * b) * (-0.16666666666666666 + ((b * b) * (0.008333333333333333 + ((b * b) * -0.0001984126984126984)))))))) / cos((b + a)); else tmp = r / (1.0 / sin(b)); end tmp_2 = tmp; end
code[r_, a_, b_] := If[LessEqual[b, -9e+27], N[(r * N[Sin[b], $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 3.8e+25], N[(N[(r * N[(b * N[(1.0 + N[(N[(b * b), $MachinePrecision] * N[(-0.16666666666666666 + N[(N[(b * b), $MachinePrecision] * N[(0.008333333333333333 + N[(N[(b * b), $MachinePrecision] * -0.0001984126984126984), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[Cos[N[(b + a), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], N[(r / N[(1.0 / N[Sin[b], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq -9 \cdot 10^{+27}:\\
\;\;\;\;r \cdot \sin b\\
\mathbf{elif}\;b \leq 3.8 \cdot 10^{+25}:\\
\;\;\;\;\frac{r \cdot \left(b \cdot \left(1 + \left(b \cdot b\right) \cdot \left(-0.16666666666666666 + \left(b \cdot b\right) \cdot \left(0.008333333333333333 + \left(b \cdot b\right) \cdot -0.0001984126984126984\right)\right)\right)\right)}{\cos \left(b + a\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{r}{\frac{1}{\sin b}}\\
\end{array}
\end{array}
if b < -8.9999999999999998e27Initial program 60.1%
/-lowering-/.f64N/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
cos-lowering-cos.f64N/A
+-commutativeN/A
+-lowering-+.f6460.1%
Simplified60.1%
associate-*l/N/A
associate-/r/N/A
*-lft-identityN/A
div-invN/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
cos-lowering-cos.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
sin-lowering-sin.f6460.0%
Applied egg-rr60.0%
Taylor expanded in b around 0
cos-lowering-cos.f6410.2%
Simplified10.2%
Taylor expanded in a around 0
*-lowering-*.f64N/A
sin-lowering-sin.f6410.5%
Simplified10.5%
if -8.9999999999999998e27 < b < 3.8e25Initial program 97.9%
/-lowering-/.f64N/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
cos-lowering-cos.f64N/A
+-commutativeN/A
+-lowering-+.f6497.9%
Simplified97.9%
Taylor expanded in b around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6493.8%
Simplified93.8%
if 3.8e25 < b Initial program 54.0%
/-lowering-/.f64N/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
cos-lowering-cos.f64N/A
+-commutativeN/A
+-lowering-+.f6454.0%
Simplified54.0%
associate-*l/N/A
associate-/r/N/A
*-lft-identityN/A
div-invN/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
cos-lowering-cos.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
sin-lowering-sin.f6453.8%
Applied egg-rr53.8%
Taylor expanded in b around 0
cos-lowering-cos.f6411.4%
Simplified11.4%
Taylor expanded in a around 0
Simplified10.9%
Final simplification54.5%
(FPCore (r a b)
:precision binary64
(if (<= b -2e+17)
(* r (sin b))
(if (<= b 15.5)
(/
(*
b
(+
r
(*
(* b b)
(* r (+ -0.16666666666666666 (* (* b b) 0.008333333333333333))))))
(cos (+ b a)))
(/ r (/ 1.0 (sin b))))))
double code(double r, double a, double b) {
double tmp;
if (b <= -2e+17) {
tmp = r * sin(b);
} else if (b <= 15.5) {
tmp = (b * (r + ((b * b) * (r * (-0.16666666666666666 + ((b * b) * 0.008333333333333333)))))) / cos((b + a));
} else {
tmp = r / (1.0 / sin(b));
}
return tmp;
}
real(8) function code(r, a, b)
real(8), intent (in) :: r
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: tmp
if (b <= (-2d+17)) then
tmp = r * sin(b)
else if (b <= 15.5d0) then
tmp = (b * (r + ((b * b) * (r * ((-0.16666666666666666d0) + ((b * b) * 0.008333333333333333d0)))))) / cos((b + a))
else
tmp = r / (1.0d0 / sin(b))
end if
code = tmp
end function
public static double code(double r, double a, double b) {
double tmp;
if (b <= -2e+17) {
tmp = r * Math.sin(b);
} else if (b <= 15.5) {
tmp = (b * (r + ((b * b) * (r * (-0.16666666666666666 + ((b * b) * 0.008333333333333333)))))) / Math.cos((b + a));
} else {
tmp = r / (1.0 / Math.sin(b));
}
return tmp;
}
def code(r, a, b): tmp = 0 if b <= -2e+17: tmp = r * math.sin(b) elif b <= 15.5: tmp = (b * (r + ((b * b) * (r * (-0.16666666666666666 + ((b * b) * 0.008333333333333333)))))) / math.cos((b + a)) else: tmp = r / (1.0 / math.sin(b)) return tmp
function code(r, a, b) tmp = 0.0 if (b <= -2e+17) tmp = Float64(r * sin(b)); elseif (b <= 15.5) tmp = Float64(Float64(b * Float64(r + Float64(Float64(b * b) * Float64(r * Float64(-0.16666666666666666 + Float64(Float64(b * b) * 0.008333333333333333)))))) / cos(Float64(b + a))); else tmp = Float64(r / Float64(1.0 / sin(b))); end return tmp end
function tmp_2 = code(r, a, b) tmp = 0.0; if (b <= -2e+17) tmp = r * sin(b); elseif (b <= 15.5) tmp = (b * (r + ((b * b) * (r * (-0.16666666666666666 + ((b * b) * 0.008333333333333333)))))) / cos((b + a)); else tmp = r / (1.0 / sin(b)); end tmp_2 = tmp; end
code[r_, a_, b_] := If[LessEqual[b, -2e+17], N[(r * N[Sin[b], $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 15.5], N[(N[(b * N[(r + N[(N[(b * b), $MachinePrecision] * N[(r * N[(-0.16666666666666666 + N[(N[(b * b), $MachinePrecision] * 0.008333333333333333), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[Cos[N[(b + a), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], N[(r / N[(1.0 / N[Sin[b], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq -2 \cdot 10^{+17}:\\
\;\;\;\;r \cdot \sin b\\
\mathbf{elif}\;b \leq 15.5:\\
\;\;\;\;\frac{b \cdot \left(r + \left(b \cdot b\right) \cdot \left(r \cdot \left(-0.16666666666666666 + \left(b \cdot b\right) \cdot 0.008333333333333333\right)\right)\right)}{\cos \left(b + a\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{r}{\frac{1}{\sin b}}\\
\end{array}
\end{array}
if b < -2e17Initial program 60.8%
/-lowering-/.f64N/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
cos-lowering-cos.f64N/A
+-commutativeN/A
+-lowering-+.f6460.8%
Simplified60.8%
associate-*l/N/A
associate-/r/N/A
*-lft-identityN/A
div-invN/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
cos-lowering-cos.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
sin-lowering-sin.f6460.6%
Applied egg-rr60.6%
Taylor expanded in b around 0
cos-lowering-cos.f6410.0%
Simplified10.0%
Taylor expanded in a around 0
*-lowering-*.f64N/A
sin-lowering-sin.f6410.4%
Simplified10.4%
if -2e17 < b < 15.5Initial program 98.4%
/-lowering-/.f64N/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
cos-lowering-cos.f64N/A
+-commutativeN/A
+-lowering-+.f6498.4%
Simplified98.4%
Taylor expanded in b around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
associate-*r*N/A
distribute-rgt-outN/A
+-commutativeN/A
metadata-evalN/A
sub-negN/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6495.7%
Simplified95.7%
if 15.5 < b Initial program 54.2%
/-lowering-/.f64N/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
cos-lowering-cos.f64N/A
+-commutativeN/A
+-lowering-+.f6454.2%
Simplified54.2%
associate-*l/N/A
associate-/r/N/A
*-lft-identityN/A
div-invN/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
cos-lowering-cos.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
sin-lowering-sin.f6454.0%
Applied egg-rr54.0%
Taylor expanded in b around 0
cos-lowering-cos.f6411.5%
Simplified11.5%
Taylor expanded in a around 0
Simplified10.6%
Final simplification54.4%
(FPCore (r a b)
:precision binary64
(if (<= b -8.5e+16)
(* r (sin b))
(if (<= b 11.0)
(*
(/ r (cos (+ b a)))
(*
b
(+
1.0
(*
b
(* b (+ -0.16666666666666666 (* b (* b 0.008333333333333333))))))))
(/ r (/ 1.0 (sin b))))))
double code(double r, double a, double b) {
double tmp;
if (b <= -8.5e+16) {
tmp = r * sin(b);
} else if (b <= 11.0) {
tmp = (r / cos((b + a))) * (b * (1.0 + (b * (b * (-0.16666666666666666 + (b * (b * 0.008333333333333333)))))));
} else {
tmp = r / (1.0 / sin(b));
}
return tmp;
}
real(8) function code(r, a, b)
real(8), intent (in) :: r
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: tmp
if (b <= (-8.5d+16)) then
tmp = r * sin(b)
else if (b <= 11.0d0) then
tmp = (r / cos((b + a))) * (b * (1.0d0 + (b * (b * ((-0.16666666666666666d0) + (b * (b * 0.008333333333333333d0)))))))
else
tmp = r / (1.0d0 / sin(b))
end if
code = tmp
end function
public static double code(double r, double a, double b) {
double tmp;
if (b <= -8.5e+16) {
tmp = r * Math.sin(b);
} else if (b <= 11.0) {
tmp = (r / Math.cos((b + a))) * (b * (1.0 + (b * (b * (-0.16666666666666666 + (b * (b * 0.008333333333333333)))))));
} else {
tmp = r / (1.0 / Math.sin(b));
}
return tmp;
}
def code(r, a, b): tmp = 0 if b <= -8.5e+16: tmp = r * math.sin(b) elif b <= 11.0: tmp = (r / math.cos((b + a))) * (b * (1.0 + (b * (b * (-0.16666666666666666 + (b * (b * 0.008333333333333333))))))) else: tmp = r / (1.0 / math.sin(b)) return tmp
function code(r, a, b) tmp = 0.0 if (b <= -8.5e+16) tmp = Float64(r * sin(b)); elseif (b <= 11.0) tmp = Float64(Float64(r / cos(Float64(b + a))) * Float64(b * Float64(1.0 + Float64(b * Float64(b * Float64(-0.16666666666666666 + Float64(b * Float64(b * 0.008333333333333333)))))))); else tmp = Float64(r / Float64(1.0 / sin(b))); end return tmp end
function tmp_2 = code(r, a, b) tmp = 0.0; if (b <= -8.5e+16) tmp = r * sin(b); elseif (b <= 11.0) tmp = (r / cos((b + a))) * (b * (1.0 + (b * (b * (-0.16666666666666666 + (b * (b * 0.008333333333333333))))))); else tmp = r / (1.0 / sin(b)); end tmp_2 = tmp; end
code[r_, a_, b_] := If[LessEqual[b, -8.5e+16], N[(r * N[Sin[b], $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 11.0], N[(N[(r / N[Cos[N[(b + a), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] * N[(b * N[(1.0 + N[(b * N[(b * N[(-0.16666666666666666 + N[(b * N[(b * 0.008333333333333333), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(r / N[(1.0 / N[Sin[b], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq -8.5 \cdot 10^{+16}:\\
\;\;\;\;r \cdot \sin b\\
\mathbf{elif}\;b \leq 11:\\
\;\;\;\;\frac{r}{\cos \left(b + a\right)} \cdot \left(b \cdot \left(1 + b \cdot \left(b \cdot \left(-0.16666666666666666 + b \cdot \left(b \cdot 0.008333333333333333\right)\right)\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{r}{\frac{1}{\sin b}}\\
\end{array}
\end{array}
if b < -8.5e16Initial program 60.8%
/-lowering-/.f64N/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
cos-lowering-cos.f64N/A
+-commutativeN/A
+-lowering-+.f6460.8%
Simplified60.8%
associate-*l/N/A
associate-/r/N/A
*-lft-identityN/A
div-invN/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
cos-lowering-cos.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
sin-lowering-sin.f6460.6%
Applied egg-rr60.6%
Taylor expanded in b around 0
cos-lowering-cos.f6410.0%
Simplified10.0%
Taylor expanded in a around 0
*-lowering-*.f64N/A
sin-lowering-sin.f6410.4%
Simplified10.4%
if -8.5e16 < b < 11Initial program 98.4%
/-lowering-/.f64N/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
cos-lowering-cos.f64N/A
+-commutativeN/A
+-lowering-+.f6498.4%
Simplified98.4%
associate-*l/N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
cos-lowering-cos.f64N/A
+-lowering-+.f64N/A
sin-lowering-sin.f6498.3%
Applied egg-rr98.3%
Taylor expanded in b around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6495.6%
Simplified95.6%
if 11 < b Initial program 54.2%
/-lowering-/.f64N/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
cos-lowering-cos.f64N/A
+-commutativeN/A
+-lowering-+.f6454.2%
Simplified54.2%
associate-*l/N/A
associate-/r/N/A
*-lft-identityN/A
div-invN/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
cos-lowering-cos.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
sin-lowering-sin.f6454.0%
Applied egg-rr54.0%
Taylor expanded in b around 0
cos-lowering-cos.f6411.5%
Simplified11.5%
Taylor expanded in a around 0
Simplified10.6%
Final simplification54.4%
(FPCore (r a b)
:precision binary64
(if (<= b -9e+27)
(* r (sin b))
(if (<= b 3.8e+25)
(/ (* (+ 1.0 (* b (* b -0.16666666666666666))) (* r b)) (cos (+ b a)))
(/ r (/ 1.0 (sin b))))))
double code(double r, double a, double b) {
double tmp;
if (b <= -9e+27) {
tmp = r * sin(b);
} else if (b <= 3.8e+25) {
tmp = ((1.0 + (b * (b * -0.16666666666666666))) * (r * b)) / cos((b + a));
} else {
tmp = r / (1.0 / sin(b));
}
return tmp;
}
real(8) function code(r, a, b)
real(8), intent (in) :: r
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: tmp
if (b <= (-9d+27)) then
tmp = r * sin(b)
else if (b <= 3.8d+25) then
tmp = ((1.0d0 + (b * (b * (-0.16666666666666666d0)))) * (r * b)) / cos((b + a))
else
tmp = r / (1.0d0 / sin(b))
end if
code = tmp
end function
public static double code(double r, double a, double b) {
double tmp;
if (b <= -9e+27) {
tmp = r * Math.sin(b);
} else if (b <= 3.8e+25) {
tmp = ((1.0 + (b * (b * -0.16666666666666666))) * (r * b)) / Math.cos((b + a));
} else {
tmp = r / (1.0 / Math.sin(b));
}
return tmp;
}
def code(r, a, b): tmp = 0 if b <= -9e+27: tmp = r * math.sin(b) elif b <= 3.8e+25: tmp = ((1.0 + (b * (b * -0.16666666666666666))) * (r * b)) / math.cos((b + a)) else: tmp = r / (1.0 / math.sin(b)) return tmp
function code(r, a, b) tmp = 0.0 if (b <= -9e+27) tmp = Float64(r * sin(b)); elseif (b <= 3.8e+25) tmp = Float64(Float64(Float64(1.0 + Float64(b * Float64(b * -0.16666666666666666))) * Float64(r * b)) / cos(Float64(b + a))); else tmp = Float64(r / Float64(1.0 / sin(b))); end return tmp end
function tmp_2 = code(r, a, b) tmp = 0.0; if (b <= -9e+27) tmp = r * sin(b); elseif (b <= 3.8e+25) tmp = ((1.0 + (b * (b * -0.16666666666666666))) * (r * b)) / cos((b + a)); else tmp = r / (1.0 / sin(b)); end tmp_2 = tmp; end
code[r_, a_, b_] := If[LessEqual[b, -9e+27], N[(r * N[Sin[b], $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 3.8e+25], N[(N[(N[(1.0 + N[(b * N[(b * -0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(r * b), $MachinePrecision]), $MachinePrecision] / N[Cos[N[(b + a), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], N[(r / N[(1.0 / N[Sin[b], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq -9 \cdot 10^{+27}:\\
\;\;\;\;r \cdot \sin b\\
\mathbf{elif}\;b \leq 3.8 \cdot 10^{+25}:\\
\;\;\;\;\frac{\left(1 + b \cdot \left(b \cdot -0.16666666666666666\right)\right) \cdot \left(r \cdot b\right)}{\cos \left(b + a\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{r}{\frac{1}{\sin b}}\\
\end{array}
\end{array}
if b < -8.9999999999999998e27Initial program 60.1%
/-lowering-/.f64N/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
cos-lowering-cos.f64N/A
+-commutativeN/A
+-lowering-+.f6460.1%
Simplified60.1%
associate-*l/N/A
associate-/r/N/A
*-lft-identityN/A
div-invN/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
cos-lowering-cos.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
sin-lowering-sin.f6460.0%
Applied egg-rr60.0%
Taylor expanded in b around 0
cos-lowering-cos.f6410.2%
Simplified10.2%
Taylor expanded in a around 0
*-lowering-*.f64N/A
sin-lowering-sin.f6410.5%
Simplified10.5%
if -8.9999999999999998e27 < b < 3.8e25Initial program 97.9%
/-lowering-/.f64N/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
cos-lowering-cos.f64N/A
+-commutativeN/A
+-lowering-+.f6497.9%
Simplified97.9%
Taylor expanded in b around 0
+-commutativeN/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
distribute-rgt-inN/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
distribute-lft1-inN/A
+-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6493.3%
Simplified93.3%
if 3.8e25 < b Initial program 54.0%
/-lowering-/.f64N/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
cos-lowering-cos.f64N/A
+-commutativeN/A
+-lowering-+.f6454.0%
Simplified54.0%
associate-*l/N/A
associate-/r/N/A
*-lft-identityN/A
div-invN/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
cos-lowering-cos.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
sin-lowering-sin.f6453.8%
Applied egg-rr53.8%
Taylor expanded in b around 0
cos-lowering-cos.f6411.4%
Simplified11.4%
Taylor expanded in a around 0
Simplified10.9%
Final simplification54.3%
(FPCore (r a b)
:precision binary64
(if (<= b -9e+27)
(* r (sin b))
(if (<= b 3.8e+25)
(/ (* r (* b (+ 1.0 (* b (* b -0.16666666666666666))))) (cos (+ b a)))
(/ r (/ 1.0 (sin b))))))
double code(double r, double a, double b) {
double tmp;
if (b <= -9e+27) {
tmp = r * sin(b);
} else if (b <= 3.8e+25) {
tmp = (r * (b * (1.0 + (b * (b * -0.16666666666666666))))) / cos((b + a));
} else {
tmp = r / (1.0 / sin(b));
}
return tmp;
}
real(8) function code(r, a, b)
real(8), intent (in) :: r
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: tmp
if (b <= (-9d+27)) then
tmp = r * sin(b)
else if (b <= 3.8d+25) then
tmp = (r * (b * (1.0d0 + (b * (b * (-0.16666666666666666d0)))))) / cos((b + a))
else
tmp = r / (1.0d0 / sin(b))
end if
code = tmp
end function
public static double code(double r, double a, double b) {
double tmp;
if (b <= -9e+27) {
tmp = r * Math.sin(b);
} else if (b <= 3.8e+25) {
tmp = (r * (b * (1.0 + (b * (b * -0.16666666666666666))))) / Math.cos((b + a));
} else {
tmp = r / (1.0 / Math.sin(b));
}
return tmp;
}
def code(r, a, b): tmp = 0 if b <= -9e+27: tmp = r * math.sin(b) elif b <= 3.8e+25: tmp = (r * (b * (1.0 + (b * (b * -0.16666666666666666))))) / math.cos((b + a)) else: tmp = r / (1.0 / math.sin(b)) return tmp
function code(r, a, b) tmp = 0.0 if (b <= -9e+27) tmp = Float64(r * sin(b)); elseif (b <= 3.8e+25) tmp = Float64(Float64(r * Float64(b * Float64(1.0 + Float64(b * Float64(b * -0.16666666666666666))))) / cos(Float64(b + a))); else tmp = Float64(r / Float64(1.0 / sin(b))); end return tmp end
function tmp_2 = code(r, a, b) tmp = 0.0; if (b <= -9e+27) tmp = r * sin(b); elseif (b <= 3.8e+25) tmp = (r * (b * (1.0 + (b * (b * -0.16666666666666666))))) / cos((b + a)); else tmp = r / (1.0 / sin(b)); end tmp_2 = tmp; end
code[r_, a_, b_] := If[LessEqual[b, -9e+27], N[(r * N[Sin[b], $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 3.8e+25], N[(N[(r * N[(b * N[(1.0 + N[(b * N[(b * -0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[Cos[N[(b + a), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], N[(r / N[(1.0 / N[Sin[b], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq -9 \cdot 10^{+27}:\\
\;\;\;\;r \cdot \sin b\\
\mathbf{elif}\;b \leq 3.8 \cdot 10^{+25}:\\
\;\;\;\;\frac{r \cdot \left(b \cdot \left(1 + b \cdot \left(b \cdot -0.16666666666666666\right)\right)\right)}{\cos \left(b + a\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{r}{\frac{1}{\sin b}}\\
\end{array}
\end{array}
if b < -8.9999999999999998e27Initial program 60.1%
/-lowering-/.f64N/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
cos-lowering-cos.f64N/A
+-commutativeN/A
+-lowering-+.f6460.1%
Simplified60.1%
associate-*l/N/A
associate-/r/N/A
*-lft-identityN/A
div-invN/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
cos-lowering-cos.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
sin-lowering-sin.f6460.0%
Applied egg-rr60.0%
Taylor expanded in b around 0
cos-lowering-cos.f6410.2%
Simplified10.2%
Taylor expanded in a around 0
*-lowering-*.f64N/A
sin-lowering-sin.f6410.5%
Simplified10.5%
if -8.9999999999999998e27 < b < 3.8e25Initial program 97.9%
/-lowering-/.f64N/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
cos-lowering-cos.f64N/A
+-commutativeN/A
+-lowering-+.f6497.9%
Simplified97.9%
Taylor expanded in b around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6493.3%
Simplified93.3%
if 3.8e25 < b Initial program 54.0%
/-lowering-/.f64N/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
cos-lowering-cos.f64N/A
+-commutativeN/A
+-lowering-+.f6454.0%
Simplified54.0%
associate-*l/N/A
associate-/r/N/A
*-lft-identityN/A
div-invN/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
cos-lowering-cos.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
sin-lowering-sin.f6453.8%
Applied egg-rr53.8%
Taylor expanded in b around 0
cos-lowering-cos.f6411.4%
Simplified11.4%
Taylor expanded in a around 0
Simplified10.9%
Final simplification54.3%
(FPCore (r a b)
:precision binary64
(if (<= b -1e+28)
(* r (sin b))
(if (<= b 4e+25)
(* (/ r (cos (+ b a))) (* b (+ 1.0 (* b (* b -0.16666666666666666)))))
(/ r (/ 1.0 (sin b))))))
double code(double r, double a, double b) {
double tmp;
if (b <= -1e+28) {
tmp = r * sin(b);
} else if (b <= 4e+25) {
tmp = (r / cos((b + a))) * (b * (1.0 + (b * (b * -0.16666666666666666))));
} else {
tmp = r / (1.0 / sin(b));
}
return tmp;
}
real(8) function code(r, a, b)
real(8), intent (in) :: r
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: tmp
if (b <= (-1d+28)) then
tmp = r * sin(b)
else if (b <= 4d+25) then
tmp = (r / cos((b + a))) * (b * (1.0d0 + (b * (b * (-0.16666666666666666d0)))))
else
tmp = r / (1.0d0 / sin(b))
end if
code = tmp
end function
public static double code(double r, double a, double b) {
double tmp;
if (b <= -1e+28) {
tmp = r * Math.sin(b);
} else if (b <= 4e+25) {
tmp = (r / Math.cos((b + a))) * (b * (1.0 + (b * (b * -0.16666666666666666))));
} else {
tmp = r / (1.0 / Math.sin(b));
}
return tmp;
}
def code(r, a, b): tmp = 0 if b <= -1e+28: tmp = r * math.sin(b) elif b <= 4e+25: tmp = (r / math.cos((b + a))) * (b * (1.0 + (b * (b * -0.16666666666666666)))) else: tmp = r / (1.0 / math.sin(b)) return tmp
function code(r, a, b) tmp = 0.0 if (b <= -1e+28) tmp = Float64(r * sin(b)); elseif (b <= 4e+25) tmp = Float64(Float64(r / cos(Float64(b + a))) * Float64(b * Float64(1.0 + Float64(b * Float64(b * -0.16666666666666666))))); else tmp = Float64(r / Float64(1.0 / sin(b))); end return tmp end
function tmp_2 = code(r, a, b) tmp = 0.0; if (b <= -1e+28) tmp = r * sin(b); elseif (b <= 4e+25) tmp = (r / cos((b + a))) * (b * (1.0 + (b * (b * -0.16666666666666666)))); else tmp = r / (1.0 / sin(b)); end tmp_2 = tmp; end
code[r_, a_, b_] := If[LessEqual[b, -1e+28], N[(r * N[Sin[b], $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 4e+25], N[(N[(r / N[Cos[N[(b + a), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] * N[(b * N[(1.0 + N[(b * N[(b * -0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(r / N[(1.0 / N[Sin[b], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq -1 \cdot 10^{+28}:\\
\;\;\;\;r \cdot \sin b\\
\mathbf{elif}\;b \leq 4 \cdot 10^{+25}:\\
\;\;\;\;\frac{r}{\cos \left(b + a\right)} \cdot \left(b \cdot \left(1 + b \cdot \left(b \cdot -0.16666666666666666\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{r}{\frac{1}{\sin b}}\\
\end{array}
\end{array}
if b < -9.99999999999999958e27Initial program 60.1%
/-lowering-/.f64N/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
cos-lowering-cos.f64N/A
+-commutativeN/A
+-lowering-+.f6460.1%
Simplified60.1%
associate-*l/N/A
associate-/r/N/A
*-lft-identityN/A
div-invN/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
cos-lowering-cos.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
sin-lowering-sin.f6460.0%
Applied egg-rr60.0%
Taylor expanded in b around 0
cos-lowering-cos.f6410.2%
Simplified10.2%
Taylor expanded in a around 0
*-lowering-*.f64N/A
sin-lowering-sin.f6410.5%
Simplified10.5%
if -9.99999999999999958e27 < b < 4.00000000000000036e25Initial program 97.9%
/-lowering-/.f64N/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
cos-lowering-cos.f64N/A
+-commutativeN/A
+-lowering-+.f6497.9%
Simplified97.9%
associate-*l/N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
cos-lowering-cos.f64N/A
+-lowering-+.f64N/A
sin-lowering-sin.f6497.8%
Applied egg-rr97.8%
Taylor expanded in b around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6493.2%
Simplified93.2%
if 4.00000000000000036e25 < b Initial program 54.0%
/-lowering-/.f64N/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
cos-lowering-cos.f64N/A
+-commutativeN/A
+-lowering-+.f6454.0%
Simplified54.0%
associate-*l/N/A
associate-/r/N/A
*-lft-identityN/A
div-invN/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
cos-lowering-cos.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
sin-lowering-sin.f6453.8%
Applied egg-rr53.8%
Taylor expanded in b around 0
cos-lowering-cos.f6411.4%
Simplified11.4%
Taylor expanded in a around 0
Simplified10.9%
Final simplification54.2%
(FPCore (r a b) :precision binary64 (if (<= b -9e+27) (* r (sin b)) (if (<= b 4e+25) (/ (* r b) (cos a)) (/ r (/ 1.0 (sin b))))))
double code(double r, double a, double b) {
double tmp;
if (b <= -9e+27) {
tmp = r * sin(b);
} else if (b <= 4e+25) {
tmp = (r * b) / cos(a);
} else {
tmp = r / (1.0 / sin(b));
}
return tmp;
}
real(8) function code(r, a, b)
real(8), intent (in) :: r
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: tmp
if (b <= (-9d+27)) then
tmp = r * sin(b)
else if (b <= 4d+25) then
tmp = (r * b) / cos(a)
else
tmp = r / (1.0d0 / sin(b))
end if
code = tmp
end function
public static double code(double r, double a, double b) {
double tmp;
if (b <= -9e+27) {
tmp = r * Math.sin(b);
} else if (b <= 4e+25) {
tmp = (r * b) / Math.cos(a);
} else {
tmp = r / (1.0 / Math.sin(b));
}
return tmp;
}
def code(r, a, b): tmp = 0 if b <= -9e+27: tmp = r * math.sin(b) elif b <= 4e+25: tmp = (r * b) / math.cos(a) else: tmp = r / (1.0 / math.sin(b)) return tmp
function code(r, a, b) tmp = 0.0 if (b <= -9e+27) tmp = Float64(r * sin(b)); elseif (b <= 4e+25) tmp = Float64(Float64(r * b) / cos(a)); else tmp = Float64(r / Float64(1.0 / sin(b))); end return tmp end
function tmp_2 = code(r, a, b) tmp = 0.0; if (b <= -9e+27) tmp = r * sin(b); elseif (b <= 4e+25) tmp = (r * b) / cos(a); else tmp = r / (1.0 / sin(b)); end tmp_2 = tmp; end
code[r_, a_, b_] := If[LessEqual[b, -9e+27], N[(r * N[Sin[b], $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 4e+25], N[(N[(r * b), $MachinePrecision] / N[Cos[a], $MachinePrecision]), $MachinePrecision], N[(r / N[(1.0 / N[Sin[b], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq -9 \cdot 10^{+27}:\\
\;\;\;\;r \cdot \sin b\\
\mathbf{elif}\;b \leq 4 \cdot 10^{+25}:\\
\;\;\;\;\frac{r \cdot b}{\cos a}\\
\mathbf{else}:\\
\;\;\;\;\frac{r}{\frac{1}{\sin b}}\\
\end{array}
\end{array}
if b < -8.9999999999999998e27Initial program 60.1%
/-lowering-/.f64N/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
cos-lowering-cos.f64N/A
+-commutativeN/A
+-lowering-+.f6460.1%
Simplified60.1%
associate-*l/N/A
associate-/r/N/A
*-lft-identityN/A
div-invN/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
cos-lowering-cos.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
sin-lowering-sin.f6460.0%
Applied egg-rr60.0%
Taylor expanded in b around 0
cos-lowering-cos.f6410.2%
Simplified10.2%
Taylor expanded in a around 0
*-lowering-*.f64N/A
sin-lowering-sin.f6410.5%
Simplified10.5%
if -8.9999999999999998e27 < b < 4.00000000000000036e25Initial program 97.9%
/-lowering-/.f64N/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
cos-lowering-cos.f64N/A
+-commutativeN/A
+-lowering-+.f6497.9%
Simplified97.9%
Taylor expanded in b around 0
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
cos-lowering-cos.f6492.9%
Simplified92.9%
if 4.00000000000000036e25 < b Initial program 54.0%
/-lowering-/.f64N/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
cos-lowering-cos.f64N/A
+-commutativeN/A
+-lowering-+.f6454.0%
Simplified54.0%
associate-*l/N/A
associate-/r/N/A
*-lft-identityN/A
div-invN/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
cos-lowering-cos.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
sin-lowering-sin.f6453.8%
Applied egg-rr53.8%
Taylor expanded in b around 0
cos-lowering-cos.f6411.4%
Simplified11.4%
Taylor expanded in a around 0
Simplified10.9%
Final simplification54.0%
(FPCore (r a b) :precision binary64 (let* ((t_0 (* r (sin b)))) (if (<= b -4e+28) t_0 (if (<= b 4e+25) (/ (* r b) (cos a)) t_0))))
double code(double r, double a, double b) {
double t_0 = r * sin(b);
double tmp;
if (b <= -4e+28) {
tmp = t_0;
} else if (b <= 4e+25) {
tmp = (r * b) / cos(a);
} else {
tmp = t_0;
}
return tmp;
}
real(8) function code(r, a, b)
real(8), intent (in) :: r
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: t_0
real(8) :: tmp
t_0 = r * sin(b)
if (b <= (-4d+28)) then
tmp = t_0
else if (b <= 4d+25) then
tmp = (r * b) / cos(a)
else
tmp = t_0
end if
code = tmp
end function
public static double code(double r, double a, double b) {
double t_0 = r * Math.sin(b);
double tmp;
if (b <= -4e+28) {
tmp = t_0;
} else if (b <= 4e+25) {
tmp = (r * b) / Math.cos(a);
} else {
tmp = t_0;
}
return tmp;
}
def code(r, a, b): t_0 = r * math.sin(b) tmp = 0 if b <= -4e+28: tmp = t_0 elif b <= 4e+25: tmp = (r * b) / math.cos(a) else: tmp = t_0 return tmp
function code(r, a, b) t_0 = Float64(r * sin(b)) tmp = 0.0 if (b <= -4e+28) tmp = t_0; elseif (b <= 4e+25) tmp = Float64(Float64(r * b) / cos(a)); else tmp = t_0; end return tmp end
function tmp_2 = code(r, a, b) t_0 = r * sin(b); tmp = 0.0; if (b <= -4e+28) tmp = t_0; elseif (b <= 4e+25) tmp = (r * b) / cos(a); else tmp = t_0; end tmp_2 = tmp; end
code[r_, a_, b_] := Block[{t$95$0 = N[(r * N[Sin[b], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[b, -4e+28], t$95$0, If[LessEqual[b, 4e+25], N[(N[(r * b), $MachinePrecision] / N[Cos[a], $MachinePrecision]), $MachinePrecision], t$95$0]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := r \cdot \sin b\\
\mathbf{if}\;b \leq -4 \cdot 10^{+28}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;b \leq 4 \cdot 10^{+25}:\\
\;\;\;\;\frac{r \cdot b}{\cos a}\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if b < -3.99999999999999983e28 or 4.00000000000000036e25 < b Initial program 57.0%
/-lowering-/.f64N/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
cos-lowering-cos.f64N/A
+-commutativeN/A
+-lowering-+.f6457.0%
Simplified57.0%
associate-*l/N/A
associate-/r/N/A
*-lft-identityN/A
div-invN/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
cos-lowering-cos.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
sin-lowering-sin.f6456.8%
Applied egg-rr56.8%
Taylor expanded in b around 0
cos-lowering-cos.f6410.8%
Simplified10.8%
Taylor expanded in a around 0
*-lowering-*.f64N/A
sin-lowering-sin.f6410.7%
Simplified10.7%
if -3.99999999999999983e28 < b < 4.00000000000000036e25Initial program 97.9%
/-lowering-/.f64N/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
cos-lowering-cos.f64N/A
+-commutativeN/A
+-lowering-+.f6497.9%
Simplified97.9%
Taylor expanded in b around 0
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
cos-lowering-cos.f6492.9%
Simplified92.9%
(FPCore (r a b) :precision binary64 (let* ((t_0 (* r (sin b)))) (if (<= b -4e+28) t_0 (if (<= b 4e+25) (* b (/ r (cos a))) t_0))))
double code(double r, double a, double b) {
double t_0 = r * sin(b);
double tmp;
if (b <= -4e+28) {
tmp = t_0;
} else if (b <= 4e+25) {
tmp = b * (r / cos(a));
} else {
tmp = t_0;
}
return tmp;
}
real(8) function code(r, a, b)
real(8), intent (in) :: r
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: t_0
real(8) :: tmp
t_0 = r * sin(b)
if (b <= (-4d+28)) then
tmp = t_0
else if (b <= 4d+25) then
tmp = b * (r / cos(a))
else
tmp = t_0
end if
code = tmp
end function
public static double code(double r, double a, double b) {
double t_0 = r * Math.sin(b);
double tmp;
if (b <= -4e+28) {
tmp = t_0;
} else if (b <= 4e+25) {
tmp = b * (r / Math.cos(a));
} else {
tmp = t_0;
}
return tmp;
}
def code(r, a, b): t_0 = r * math.sin(b) tmp = 0 if b <= -4e+28: tmp = t_0 elif b <= 4e+25: tmp = b * (r / math.cos(a)) else: tmp = t_0 return tmp
function code(r, a, b) t_0 = Float64(r * sin(b)) tmp = 0.0 if (b <= -4e+28) tmp = t_0; elseif (b <= 4e+25) tmp = Float64(b * Float64(r / cos(a))); else tmp = t_0; end return tmp end
function tmp_2 = code(r, a, b) t_0 = r * sin(b); tmp = 0.0; if (b <= -4e+28) tmp = t_0; elseif (b <= 4e+25) tmp = b * (r / cos(a)); else tmp = t_0; end tmp_2 = tmp; end
code[r_, a_, b_] := Block[{t$95$0 = N[(r * N[Sin[b], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[b, -4e+28], t$95$0, If[LessEqual[b, 4e+25], N[(b * N[(r / N[Cos[a], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$0]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := r \cdot \sin b\\
\mathbf{if}\;b \leq -4 \cdot 10^{+28}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;b \leq 4 \cdot 10^{+25}:\\
\;\;\;\;b \cdot \frac{r}{\cos a}\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if b < -3.99999999999999983e28 or 4.00000000000000036e25 < b Initial program 57.0%
/-lowering-/.f64N/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
cos-lowering-cos.f64N/A
+-commutativeN/A
+-lowering-+.f6457.0%
Simplified57.0%
associate-*l/N/A
associate-/r/N/A
*-lft-identityN/A
div-invN/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
cos-lowering-cos.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
sin-lowering-sin.f6456.8%
Applied egg-rr56.8%
Taylor expanded in b around 0
cos-lowering-cos.f6410.8%
Simplified10.8%
Taylor expanded in a around 0
*-lowering-*.f64N/A
sin-lowering-sin.f6410.7%
Simplified10.7%
if -3.99999999999999983e28 < b < 4.00000000000000036e25Initial program 97.9%
/-lowering-/.f64N/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
cos-lowering-cos.f64N/A
+-commutativeN/A
+-lowering-+.f6497.9%
Simplified97.9%
Taylor expanded in b around 0
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
cos-lowering-cos.f6492.9%
Simplified92.9%
associate-*l/N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
cos-lowering-cos.f6492.8%
Applied egg-rr92.8%
Final simplification54.0%
(FPCore (r a b) :precision binary64 (* r (sin b)))
double code(double r, double a, double b) {
return r * sin(b);
}
real(8) function code(r, a, b)
real(8), intent (in) :: r
real(8), intent (in) :: a
real(8), intent (in) :: b
code = r * sin(b)
end function
public static double code(double r, double a, double b) {
return r * Math.sin(b);
}
def code(r, a, b): return r * math.sin(b)
function code(r, a, b) return Float64(r * sin(b)) end
function tmp = code(r, a, b) tmp = r * sin(b); end
code[r_, a_, b_] := N[(r * N[Sin[b], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
r \cdot \sin b
\end{array}
Initial program 78.5%
/-lowering-/.f64N/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
cos-lowering-cos.f64N/A
+-commutativeN/A
+-lowering-+.f6478.5%
Simplified78.5%
associate-*l/N/A
associate-/r/N/A
*-lft-identityN/A
div-invN/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
cos-lowering-cos.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
sin-lowering-sin.f6478.3%
Applied egg-rr78.3%
Taylor expanded in b around 0
cos-lowering-cos.f6453.9%
Simplified53.9%
Taylor expanded in a around 0
*-lowering-*.f64N/A
sin-lowering-sin.f6437.2%
Simplified37.2%
(FPCore (r a b) :precision binary64 (* r b))
double code(double r, double a, double b) {
return r * b;
}
real(8) function code(r, a, b)
real(8), intent (in) :: r
real(8), intent (in) :: a
real(8), intent (in) :: b
code = r * b
end function
public static double code(double r, double a, double b) {
return r * b;
}
def code(r, a, b): return r * b
function code(r, a, b) return Float64(r * b) end
function tmp = code(r, a, b) tmp = r * b; end
code[r_, a_, b_] := N[(r * b), $MachinePrecision]
\begin{array}{l}
\\
r \cdot b
\end{array}
Initial program 78.5%
/-lowering-/.f64N/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
cos-lowering-cos.f64N/A
+-commutativeN/A
+-lowering-+.f6478.5%
Simplified78.5%
Taylor expanded in b around 0
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
cos-lowering-cos.f6450.7%
Simplified50.7%
Taylor expanded in a around 0
*-commutativeN/A
*-lowering-*.f6433.9%
Simplified33.9%
herbie shell --seed 2024159
(FPCore (r a b)
:name "rsin A (should all be same)"
:precision binary64
(/ (* r (sin b)) (cos (+ a b))))