
(FPCore (re im) :precision binary64 (* (* 0.5 (cos re)) (- (exp (- 0.0 im)) (exp im))))
double code(double re, double im) {
return (0.5 * cos(re)) * (exp((0.0 - im)) - exp(im));
}
real(8) function code(re, im)
real(8), intent (in) :: re
real(8), intent (in) :: im
code = (0.5d0 * cos(re)) * (exp((0.0d0 - im)) - exp(im))
end function
public static double code(double re, double im) {
return (0.5 * Math.cos(re)) * (Math.exp((0.0 - im)) - Math.exp(im));
}
def code(re, im): return (0.5 * math.cos(re)) * (math.exp((0.0 - im)) - math.exp(im))
function code(re, im) return Float64(Float64(0.5 * cos(re)) * Float64(exp(Float64(0.0 - im)) - exp(im))) end
function tmp = code(re, im) tmp = (0.5 * cos(re)) * (exp((0.0 - im)) - exp(im)); end
code[re_, im_] := N[(N[(0.5 * N[Cos[re], $MachinePrecision]), $MachinePrecision] * N[(N[Exp[N[(0.0 - im), $MachinePrecision]], $MachinePrecision] - N[Exp[im], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(0.5 \cdot \cos re\right) \cdot \left(e^{0 - im} - e^{im}\right)
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 28 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (re im) :precision binary64 (* (* 0.5 (cos re)) (- (exp (- 0.0 im)) (exp im))))
double code(double re, double im) {
return (0.5 * cos(re)) * (exp((0.0 - im)) - exp(im));
}
real(8) function code(re, im)
real(8), intent (in) :: re
real(8), intent (in) :: im
code = (0.5d0 * cos(re)) * (exp((0.0d0 - im)) - exp(im))
end function
public static double code(double re, double im) {
return (0.5 * Math.cos(re)) * (Math.exp((0.0 - im)) - Math.exp(im));
}
def code(re, im): return (0.5 * math.cos(re)) * (math.exp((0.0 - im)) - math.exp(im))
function code(re, im) return Float64(Float64(0.5 * cos(re)) * Float64(exp(Float64(0.0 - im)) - exp(im))) end
function tmp = code(re, im) tmp = (0.5 * cos(re)) * (exp((0.0 - im)) - exp(im)); end
code[re_, im_] := N[(N[(0.5 * N[Cos[re], $MachinePrecision]), $MachinePrecision] * N[(N[Exp[N[(0.0 - im), $MachinePrecision]], $MachinePrecision] - N[Exp[im], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(0.5 \cdot \cos re\right) \cdot \left(e^{0 - im} - e^{im}\right)
\end{array}
im\_m = (fabs.f64 im)
im\_s = (copysign.f64 #s(literal 1 binary64) im)
(FPCore (im_s re im_m)
:precision binary64
(let* ((t_0 (- (exp (- 0.0 im_m)) (exp im_m))))
(*
im_s
(if (<= t_0 -0.2)
(* (* 0.5 (cos re)) (/ 1.0 (/ 1.0 t_0)))
(*
(* im_m (cos re))
(+
-1.0
(*
(* im_m im_m)
(+
-0.16666666666666666
(* (* im_m im_m) -0.008333333333333333)))))))))im\_m = fabs(im);
im\_s = copysign(1.0, im);
double code(double im_s, double re, double im_m) {
double t_0 = exp((0.0 - im_m)) - exp(im_m);
double tmp;
if (t_0 <= -0.2) {
tmp = (0.5 * cos(re)) * (1.0 / (1.0 / t_0));
} else {
tmp = (im_m * cos(re)) * (-1.0 + ((im_m * im_m) * (-0.16666666666666666 + ((im_m * im_m) * -0.008333333333333333))));
}
return im_s * tmp;
}
im\_m = abs(im)
im\_s = copysign(1.0d0, im)
real(8) function code(im_s, re, im_m)
real(8), intent (in) :: im_s
real(8), intent (in) :: re
real(8), intent (in) :: im_m
real(8) :: t_0
real(8) :: tmp
t_0 = exp((0.0d0 - im_m)) - exp(im_m)
if (t_0 <= (-0.2d0)) then
tmp = (0.5d0 * cos(re)) * (1.0d0 / (1.0d0 / t_0))
else
tmp = (im_m * cos(re)) * ((-1.0d0) + ((im_m * im_m) * ((-0.16666666666666666d0) + ((im_m * im_m) * (-0.008333333333333333d0)))))
end if
code = im_s * tmp
end function
im\_m = Math.abs(im);
im\_s = Math.copySign(1.0, im);
public static double code(double im_s, double re, double im_m) {
double t_0 = Math.exp((0.0 - im_m)) - Math.exp(im_m);
double tmp;
if (t_0 <= -0.2) {
tmp = (0.5 * Math.cos(re)) * (1.0 / (1.0 / t_0));
} else {
tmp = (im_m * Math.cos(re)) * (-1.0 + ((im_m * im_m) * (-0.16666666666666666 + ((im_m * im_m) * -0.008333333333333333))));
}
return im_s * tmp;
}
im\_m = math.fabs(im) im\_s = math.copysign(1.0, im) def code(im_s, re, im_m): t_0 = math.exp((0.0 - im_m)) - math.exp(im_m) tmp = 0 if t_0 <= -0.2: tmp = (0.5 * math.cos(re)) * (1.0 / (1.0 / t_0)) else: tmp = (im_m * math.cos(re)) * (-1.0 + ((im_m * im_m) * (-0.16666666666666666 + ((im_m * im_m) * -0.008333333333333333)))) return im_s * tmp
im\_m = abs(im) im\_s = copysign(1.0, im) function code(im_s, re, im_m) t_0 = Float64(exp(Float64(0.0 - im_m)) - exp(im_m)) tmp = 0.0 if (t_0 <= -0.2) tmp = Float64(Float64(0.5 * cos(re)) * Float64(1.0 / Float64(1.0 / t_0))); else tmp = Float64(Float64(im_m * cos(re)) * Float64(-1.0 + Float64(Float64(im_m * im_m) * Float64(-0.16666666666666666 + Float64(Float64(im_m * im_m) * -0.008333333333333333))))); end return Float64(im_s * tmp) end
im\_m = abs(im); im\_s = sign(im) * abs(1.0); function tmp_2 = code(im_s, re, im_m) t_0 = exp((0.0 - im_m)) - exp(im_m); tmp = 0.0; if (t_0 <= -0.2) tmp = (0.5 * cos(re)) * (1.0 / (1.0 / t_0)); else tmp = (im_m * cos(re)) * (-1.0 + ((im_m * im_m) * (-0.16666666666666666 + ((im_m * im_m) * -0.008333333333333333)))); end tmp_2 = im_s * tmp; end
im\_m = N[Abs[im], $MachinePrecision]
im\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[im]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[im$95$s_, re_, im$95$m_] := Block[{t$95$0 = N[(N[Exp[N[(0.0 - im$95$m), $MachinePrecision]], $MachinePrecision] - N[Exp[im$95$m], $MachinePrecision]), $MachinePrecision]}, N[(im$95$s * If[LessEqual[t$95$0, -0.2], N[(N[(0.5 * N[Cos[re], $MachinePrecision]), $MachinePrecision] * N[(1.0 / N[(1.0 / t$95$0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(im$95$m * N[Cos[re], $MachinePrecision]), $MachinePrecision] * N[(-1.0 + N[(N[(im$95$m * im$95$m), $MachinePrecision] * N[(-0.16666666666666666 + N[(N[(im$95$m * im$95$m), $MachinePrecision] * -0.008333333333333333), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]]
\begin{array}{l}
im\_m = \left|im\right|
\\
im\_s = \mathsf{copysign}\left(1, im\right)
\\
\begin{array}{l}
t_0 := e^{0 - im\_m} - e^{im\_m}\\
im\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_0 \leq -0.2:\\
\;\;\;\;\left(0.5 \cdot \cos re\right) \cdot \frac{1}{\frac{1}{t\_0}}\\
\mathbf{else}:\\
\;\;\;\;\left(im\_m \cdot \cos re\right) \cdot \left(-1 + \left(im\_m \cdot im\_m\right) \cdot \left(-0.16666666666666666 + \left(im\_m \cdot im\_m\right) \cdot -0.008333333333333333\right)\right)\\
\end{array}
\end{array}
\end{array}
if (-.f64 (exp.f64 (-.f64 #s(literal 0 binary64) im)) (exp.f64 im)) < -0.20000000000000001Initial program 99.9%
flip--N/A
clear-numN/A
exp-0N/A
/-lowering-/.f64N/A
exp-0N/A
clear-numN/A
exp-0N/A
flip--N/A
/-lowering-/.f64N/A
exp-0N/A
--lowering--.f64N/A
exp-lowering-exp.f64N/A
--lowering--.f64N/A
exp-lowering-exp.f6499.9%
Applied egg-rr99.9%
if -0.20000000000000001 < (-.f64 (exp.f64 (-.f64 #s(literal 0 binary64) im)) (exp.f64 im)) Initial program 45.6%
Taylor expanded in im around 0
*-lowering-*.f64N/A
*-commutativeN/A
+-commutativeN/A
distribute-lft-inN/A
associate-*r*N/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
distribute-rgt-outN/A
distribute-lft-outN/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
Simplified90.5%
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
cos-lowering-cos.f6490.5%
Applied egg-rr90.5%
Final simplification92.9%
im\_m = (fabs.f64 im)
im\_s = (copysign.f64 #s(literal 1 binary64) im)
(FPCore (im_s re im_m)
:precision binary64
(let* ((t_0 (- (exp (- 0.0 im_m)) (exp im_m))))
(*
im_s
(if (<= t_0 -0.2)
(* t_0 (* 0.5 (cos re)))
(*
(* im_m (cos re))
(+
-1.0
(*
(* im_m im_m)
(+
-0.16666666666666666
(* (* im_m im_m) -0.008333333333333333)))))))))im\_m = fabs(im);
im\_s = copysign(1.0, im);
double code(double im_s, double re, double im_m) {
double t_0 = exp((0.0 - im_m)) - exp(im_m);
double tmp;
if (t_0 <= -0.2) {
tmp = t_0 * (0.5 * cos(re));
} else {
tmp = (im_m * cos(re)) * (-1.0 + ((im_m * im_m) * (-0.16666666666666666 + ((im_m * im_m) * -0.008333333333333333))));
}
return im_s * tmp;
}
im\_m = abs(im)
im\_s = copysign(1.0d0, im)
real(8) function code(im_s, re, im_m)
real(8), intent (in) :: im_s
real(8), intent (in) :: re
real(8), intent (in) :: im_m
real(8) :: t_0
real(8) :: tmp
t_0 = exp((0.0d0 - im_m)) - exp(im_m)
if (t_0 <= (-0.2d0)) then
tmp = t_0 * (0.5d0 * cos(re))
else
tmp = (im_m * cos(re)) * ((-1.0d0) + ((im_m * im_m) * ((-0.16666666666666666d0) + ((im_m * im_m) * (-0.008333333333333333d0)))))
end if
code = im_s * tmp
end function
im\_m = Math.abs(im);
im\_s = Math.copySign(1.0, im);
public static double code(double im_s, double re, double im_m) {
double t_0 = Math.exp((0.0 - im_m)) - Math.exp(im_m);
double tmp;
if (t_0 <= -0.2) {
tmp = t_0 * (0.5 * Math.cos(re));
} else {
tmp = (im_m * Math.cos(re)) * (-1.0 + ((im_m * im_m) * (-0.16666666666666666 + ((im_m * im_m) * -0.008333333333333333))));
}
return im_s * tmp;
}
im\_m = math.fabs(im) im\_s = math.copysign(1.0, im) def code(im_s, re, im_m): t_0 = math.exp((0.0 - im_m)) - math.exp(im_m) tmp = 0 if t_0 <= -0.2: tmp = t_0 * (0.5 * math.cos(re)) else: tmp = (im_m * math.cos(re)) * (-1.0 + ((im_m * im_m) * (-0.16666666666666666 + ((im_m * im_m) * -0.008333333333333333)))) return im_s * tmp
im\_m = abs(im) im\_s = copysign(1.0, im) function code(im_s, re, im_m) t_0 = Float64(exp(Float64(0.0 - im_m)) - exp(im_m)) tmp = 0.0 if (t_0 <= -0.2) tmp = Float64(t_0 * Float64(0.5 * cos(re))); else tmp = Float64(Float64(im_m * cos(re)) * Float64(-1.0 + Float64(Float64(im_m * im_m) * Float64(-0.16666666666666666 + Float64(Float64(im_m * im_m) * -0.008333333333333333))))); end return Float64(im_s * tmp) end
im\_m = abs(im); im\_s = sign(im) * abs(1.0); function tmp_2 = code(im_s, re, im_m) t_0 = exp((0.0 - im_m)) - exp(im_m); tmp = 0.0; if (t_0 <= -0.2) tmp = t_0 * (0.5 * cos(re)); else tmp = (im_m * cos(re)) * (-1.0 + ((im_m * im_m) * (-0.16666666666666666 + ((im_m * im_m) * -0.008333333333333333)))); end tmp_2 = im_s * tmp; end
im\_m = N[Abs[im], $MachinePrecision]
im\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[im]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[im$95$s_, re_, im$95$m_] := Block[{t$95$0 = N[(N[Exp[N[(0.0 - im$95$m), $MachinePrecision]], $MachinePrecision] - N[Exp[im$95$m], $MachinePrecision]), $MachinePrecision]}, N[(im$95$s * If[LessEqual[t$95$0, -0.2], N[(t$95$0 * N[(0.5 * N[Cos[re], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(im$95$m * N[Cos[re], $MachinePrecision]), $MachinePrecision] * N[(-1.0 + N[(N[(im$95$m * im$95$m), $MachinePrecision] * N[(-0.16666666666666666 + N[(N[(im$95$m * im$95$m), $MachinePrecision] * -0.008333333333333333), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]]
\begin{array}{l}
im\_m = \left|im\right|
\\
im\_s = \mathsf{copysign}\left(1, im\right)
\\
\begin{array}{l}
t_0 := e^{0 - im\_m} - e^{im\_m}\\
im\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_0 \leq -0.2:\\
\;\;\;\;t\_0 \cdot \left(0.5 \cdot \cos re\right)\\
\mathbf{else}:\\
\;\;\;\;\left(im\_m \cdot \cos re\right) \cdot \left(-1 + \left(im\_m \cdot im\_m\right) \cdot \left(-0.16666666666666666 + \left(im\_m \cdot im\_m\right) \cdot -0.008333333333333333\right)\right)\\
\end{array}
\end{array}
\end{array}
if (-.f64 (exp.f64 (-.f64 #s(literal 0 binary64) im)) (exp.f64 im)) < -0.20000000000000001Initial program 99.9%
if -0.20000000000000001 < (-.f64 (exp.f64 (-.f64 #s(literal 0 binary64) im)) (exp.f64 im)) Initial program 45.6%
Taylor expanded in im around 0
*-lowering-*.f64N/A
*-commutativeN/A
+-commutativeN/A
distribute-lft-inN/A
associate-*r*N/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
distribute-rgt-outN/A
distribute-lft-outN/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
Simplified90.5%
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
cos-lowering-cos.f6490.5%
Applied egg-rr90.5%
Final simplification92.9%
im\_m = (fabs.f64 im)
im\_s = (copysign.f64 #s(literal 1 binary64) im)
(FPCore (im_s re im_m)
:precision binary64
(let* ((t_0
(+
-0.16666666666666666
(*
im_m
(*
im_m
(+
-0.008333333333333333
(* im_m (* im_m -0.0001984126984126984)))))))
(t_1 (* im_m (* im_m t_0)))
(t_2 (+ -0.008333333333333333 (* (* im_m im_m) 0.0001984126984126984)))
(t_3 (* (* im_m im_m) (* im_m im_m)))
(t_4
(+
1.0
(*
(* im_m im_m)
(+
-0.16666666666666666
(*
-0.16666666666666666
(* im_m (* im_m -0.16666666666666666))))))))
(*
im_s
(if (<= im_m 1.0)
(*
(* im_m (cos re))
(/ (+ -1.0 (* t_1 (* t_3 (* t_0 t_0)))) (+ 1.0 (* t_1 (- t_1 -1.0)))))
(if (<= im_m 2.5e+24)
(+ (* (exp im_m) -0.5) (/ 0.5 (exp im_m)))
(if (<= im_m 4.6e+44)
(/
(*
im_m
(*
(cos re)
(+
(* (+ -1.0 (* t_3 (* (* im_m im_m) -0.004629629629629629))) t_2)
(*
t_3
(*
t_4
(+ 6.944444444444444e-5 (* t_3 -3.936759889140842e-8)))))))
(* t_2 t_4))
(*
im_m
(*
(cos re)
(+
-1.0
(*
(* im_m im_m)
(+
-0.16666666666666666
(*
im_m
(*
im_m
(+
-0.008333333333333333
(* (* im_m im_m) -0.0001984126984126984)))))))))))))))im\_m = fabs(im);
im\_s = copysign(1.0, im);
double code(double im_s, double re, double im_m) {
double t_0 = -0.16666666666666666 + (im_m * (im_m * (-0.008333333333333333 + (im_m * (im_m * -0.0001984126984126984)))));
double t_1 = im_m * (im_m * t_0);
double t_2 = -0.008333333333333333 + ((im_m * im_m) * 0.0001984126984126984);
double t_3 = (im_m * im_m) * (im_m * im_m);
double t_4 = 1.0 + ((im_m * im_m) * (-0.16666666666666666 + (-0.16666666666666666 * (im_m * (im_m * -0.16666666666666666)))));
double tmp;
if (im_m <= 1.0) {
tmp = (im_m * cos(re)) * ((-1.0 + (t_1 * (t_3 * (t_0 * t_0)))) / (1.0 + (t_1 * (t_1 - -1.0))));
} else if (im_m <= 2.5e+24) {
tmp = (exp(im_m) * -0.5) + (0.5 / exp(im_m));
} else if (im_m <= 4.6e+44) {
tmp = (im_m * (cos(re) * (((-1.0 + (t_3 * ((im_m * im_m) * -0.004629629629629629))) * t_2) + (t_3 * (t_4 * (6.944444444444444e-5 + (t_3 * -3.936759889140842e-8))))))) / (t_2 * t_4);
} else {
tmp = im_m * (cos(re) * (-1.0 + ((im_m * im_m) * (-0.16666666666666666 + (im_m * (im_m * (-0.008333333333333333 + ((im_m * im_m) * -0.0001984126984126984))))))));
}
return im_s * tmp;
}
im\_m = abs(im)
im\_s = copysign(1.0d0, im)
real(8) function code(im_s, re, im_m)
real(8), intent (in) :: im_s
real(8), intent (in) :: re
real(8), intent (in) :: im_m
real(8) :: t_0
real(8) :: t_1
real(8) :: t_2
real(8) :: t_3
real(8) :: t_4
real(8) :: tmp
t_0 = (-0.16666666666666666d0) + (im_m * (im_m * ((-0.008333333333333333d0) + (im_m * (im_m * (-0.0001984126984126984d0))))))
t_1 = im_m * (im_m * t_0)
t_2 = (-0.008333333333333333d0) + ((im_m * im_m) * 0.0001984126984126984d0)
t_3 = (im_m * im_m) * (im_m * im_m)
t_4 = 1.0d0 + ((im_m * im_m) * ((-0.16666666666666666d0) + ((-0.16666666666666666d0) * (im_m * (im_m * (-0.16666666666666666d0))))))
if (im_m <= 1.0d0) then
tmp = (im_m * cos(re)) * (((-1.0d0) + (t_1 * (t_3 * (t_0 * t_0)))) / (1.0d0 + (t_1 * (t_1 - (-1.0d0)))))
else if (im_m <= 2.5d+24) then
tmp = (exp(im_m) * (-0.5d0)) + (0.5d0 / exp(im_m))
else if (im_m <= 4.6d+44) then
tmp = (im_m * (cos(re) * ((((-1.0d0) + (t_3 * ((im_m * im_m) * (-0.004629629629629629d0)))) * t_2) + (t_3 * (t_4 * (6.944444444444444d-5 + (t_3 * (-3.936759889140842d-8)))))))) / (t_2 * t_4)
else
tmp = im_m * (cos(re) * ((-1.0d0) + ((im_m * im_m) * ((-0.16666666666666666d0) + (im_m * (im_m * ((-0.008333333333333333d0) + ((im_m * im_m) * (-0.0001984126984126984d0)))))))))
end if
code = im_s * tmp
end function
im\_m = Math.abs(im);
im\_s = Math.copySign(1.0, im);
public static double code(double im_s, double re, double im_m) {
double t_0 = -0.16666666666666666 + (im_m * (im_m * (-0.008333333333333333 + (im_m * (im_m * -0.0001984126984126984)))));
double t_1 = im_m * (im_m * t_0);
double t_2 = -0.008333333333333333 + ((im_m * im_m) * 0.0001984126984126984);
double t_3 = (im_m * im_m) * (im_m * im_m);
double t_4 = 1.0 + ((im_m * im_m) * (-0.16666666666666666 + (-0.16666666666666666 * (im_m * (im_m * -0.16666666666666666)))));
double tmp;
if (im_m <= 1.0) {
tmp = (im_m * Math.cos(re)) * ((-1.0 + (t_1 * (t_3 * (t_0 * t_0)))) / (1.0 + (t_1 * (t_1 - -1.0))));
} else if (im_m <= 2.5e+24) {
tmp = (Math.exp(im_m) * -0.5) + (0.5 / Math.exp(im_m));
} else if (im_m <= 4.6e+44) {
tmp = (im_m * (Math.cos(re) * (((-1.0 + (t_3 * ((im_m * im_m) * -0.004629629629629629))) * t_2) + (t_3 * (t_4 * (6.944444444444444e-5 + (t_3 * -3.936759889140842e-8))))))) / (t_2 * t_4);
} else {
tmp = im_m * (Math.cos(re) * (-1.0 + ((im_m * im_m) * (-0.16666666666666666 + (im_m * (im_m * (-0.008333333333333333 + ((im_m * im_m) * -0.0001984126984126984))))))));
}
return im_s * tmp;
}
im\_m = math.fabs(im) im\_s = math.copysign(1.0, im) def code(im_s, re, im_m): t_0 = -0.16666666666666666 + (im_m * (im_m * (-0.008333333333333333 + (im_m * (im_m * -0.0001984126984126984))))) t_1 = im_m * (im_m * t_0) t_2 = -0.008333333333333333 + ((im_m * im_m) * 0.0001984126984126984) t_3 = (im_m * im_m) * (im_m * im_m) t_4 = 1.0 + ((im_m * im_m) * (-0.16666666666666666 + (-0.16666666666666666 * (im_m * (im_m * -0.16666666666666666))))) tmp = 0 if im_m <= 1.0: tmp = (im_m * math.cos(re)) * ((-1.0 + (t_1 * (t_3 * (t_0 * t_0)))) / (1.0 + (t_1 * (t_1 - -1.0)))) elif im_m <= 2.5e+24: tmp = (math.exp(im_m) * -0.5) + (0.5 / math.exp(im_m)) elif im_m <= 4.6e+44: tmp = (im_m * (math.cos(re) * (((-1.0 + (t_3 * ((im_m * im_m) * -0.004629629629629629))) * t_2) + (t_3 * (t_4 * (6.944444444444444e-5 + (t_3 * -3.936759889140842e-8))))))) / (t_2 * t_4) else: tmp = im_m * (math.cos(re) * (-1.0 + ((im_m * im_m) * (-0.16666666666666666 + (im_m * (im_m * (-0.008333333333333333 + ((im_m * im_m) * -0.0001984126984126984)))))))) return im_s * tmp
im\_m = abs(im) im\_s = copysign(1.0, im) function code(im_s, re, im_m) t_0 = Float64(-0.16666666666666666 + Float64(im_m * Float64(im_m * Float64(-0.008333333333333333 + Float64(im_m * Float64(im_m * -0.0001984126984126984)))))) t_1 = Float64(im_m * Float64(im_m * t_0)) t_2 = Float64(-0.008333333333333333 + Float64(Float64(im_m * im_m) * 0.0001984126984126984)) t_3 = Float64(Float64(im_m * im_m) * Float64(im_m * im_m)) t_4 = Float64(1.0 + Float64(Float64(im_m * im_m) * Float64(-0.16666666666666666 + Float64(-0.16666666666666666 * Float64(im_m * Float64(im_m * -0.16666666666666666)))))) tmp = 0.0 if (im_m <= 1.0) tmp = Float64(Float64(im_m * cos(re)) * Float64(Float64(-1.0 + Float64(t_1 * Float64(t_3 * Float64(t_0 * t_0)))) / Float64(1.0 + Float64(t_1 * Float64(t_1 - -1.0))))); elseif (im_m <= 2.5e+24) tmp = Float64(Float64(exp(im_m) * -0.5) + Float64(0.5 / exp(im_m))); elseif (im_m <= 4.6e+44) tmp = Float64(Float64(im_m * Float64(cos(re) * Float64(Float64(Float64(-1.0 + Float64(t_3 * Float64(Float64(im_m * im_m) * -0.004629629629629629))) * t_2) + Float64(t_3 * Float64(t_4 * Float64(6.944444444444444e-5 + Float64(t_3 * -3.936759889140842e-8))))))) / Float64(t_2 * t_4)); else tmp = Float64(im_m * Float64(cos(re) * Float64(-1.0 + Float64(Float64(im_m * im_m) * Float64(-0.16666666666666666 + Float64(im_m * Float64(im_m * Float64(-0.008333333333333333 + Float64(Float64(im_m * im_m) * -0.0001984126984126984))))))))); end return Float64(im_s * tmp) end
im\_m = abs(im); im\_s = sign(im) * abs(1.0); function tmp_2 = code(im_s, re, im_m) t_0 = -0.16666666666666666 + (im_m * (im_m * (-0.008333333333333333 + (im_m * (im_m * -0.0001984126984126984))))); t_1 = im_m * (im_m * t_0); t_2 = -0.008333333333333333 + ((im_m * im_m) * 0.0001984126984126984); t_3 = (im_m * im_m) * (im_m * im_m); t_4 = 1.0 + ((im_m * im_m) * (-0.16666666666666666 + (-0.16666666666666666 * (im_m * (im_m * -0.16666666666666666))))); tmp = 0.0; if (im_m <= 1.0) tmp = (im_m * cos(re)) * ((-1.0 + (t_1 * (t_3 * (t_0 * t_0)))) / (1.0 + (t_1 * (t_1 - -1.0)))); elseif (im_m <= 2.5e+24) tmp = (exp(im_m) * -0.5) + (0.5 / exp(im_m)); elseif (im_m <= 4.6e+44) tmp = (im_m * (cos(re) * (((-1.0 + (t_3 * ((im_m * im_m) * -0.004629629629629629))) * t_2) + (t_3 * (t_4 * (6.944444444444444e-5 + (t_3 * -3.936759889140842e-8))))))) / (t_2 * t_4); else tmp = im_m * (cos(re) * (-1.0 + ((im_m * im_m) * (-0.16666666666666666 + (im_m * (im_m * (-0.008333333333333333 + ((im_m * im_m) * -0.0001984126984126984)))))))); end tmp_2 = im_s * tmp; end
im\_m = N[Abs[im], $MachinePrecision]
im\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[im]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[im$95$s_, re_, im$95$m_] := Block[{t$95$0 = N[(-0.16666666666666666 + N[(im$95$m * N[(im$95$m * N[(-0.008333333333333333 + N[(im$95$m * N[(im$95$m * -0.0001984126984126984), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(im$95$m * N[(im$95$m * t$95$0), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(-0.008333333333333333 + N[(N[(im$95$m * im$95$m), $MachinePrecision] * 0.0001984126984126984), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = N[(N[(im$95$m * im$95$m), $MachinePrecision] * N[(im$95$m * im$95$m), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$4 = N[(1.0 + N[(N[(im$95$m * im$95$m), $MachinePrecision] * N[(-0.16666666666666666 + N[(-0.16666666666666666 * N[(im$95$m * N[(im$95$m * -0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, N[(im$95$s * If[LessEqual[im$95$m, 1.0], N[(N[(im$95$m * N[Cos[re], $MachinePrecision]), $MachinePrecision] * N[(N[(-1.0 + N[(t$95$1 * N[(t$95$3 * N[(t$95$0 * t$95$0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(1.0 + N[(t$95$1 * N[(t$95$1 - -1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[im$95$m, 2.5e+24], N[(N[(N[Exp[im$95$m], $MachinePrecision] * -0.5), $MachinePrecision] + N[(0.5 / N[Exp[im$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[im$95$m, 4.6e+44], N[(N[(im$95$m * N[(N[Cos[re], $MachinePrecision] * N[(N[(N[(-1.0 + N[(t$95$3 * N[(N[(im$95$m * im$95$m), $MachinePrecision] * -0.004629629629629629), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * t$95$2), $MachinePrecision] + N[(t$95$3 * N[(t$95$4 * N[(6.944444444444444e-5 + N[(t$95$3 * -3.936759889140842e-8), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(t$95$2 * t$95$4), $MachinePrecision]), $MachinePrecision], N[(im$95$m * N[(N[Cos[re], $MachinePrecision] * N[(-1.0 + N[(N[(im$95$m * im$95$m), $MachinePrecision] * N[(-0.16666666666666666 + N[(im$95$m * N[(im$95$m * N[(-0.008333333333333333 + N[(N[(im$95$m * im$95$m), $MachinePrecision] * -0.0001984126984126984), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]), $MachinePrecision]]]]]]
\begin{array}{l}
im\_m = \left|im\right|
\\
im\_s = \mathsf{copysign}\left(1, im\right)
\\
\begin{array}{l}
t_0 := -0.16666666666666666 + im\_m \cdot \left(im\_m \cdot \left(-0.008333333333333333 + im\_m \cdot \left(im\_m \cdot -0.0001984126984126984\right)\right)\right)\\
t_1 := im\_m \cdot \left(im\_m \cdot t\_0\right)\\
t_2 := -0.008333333333333333 + \left(im\_m \cdot im\_m\right) \cdot 0.0001984126984126984\\
t_3 := \left(im\_m \cdot im\_m\right) \cdot \left(im\_m \cdot im\_m\right)\\
t_4 := 1 + \left(im\_m \cdot im\_m\right) \cdot \left(-0.16666666666666666 + -0.16666666666666666 \cdot \left(im\_m \cdot \left(im\_m \cdot -0.16666666666666666\right)\right)\right)\\
im\_s \cdot \begin{array}{l}
\mathbf{if}\;im\_m \leq 1:\\
\;\;\;\;\left(im\_m \cdot \cos re\right) \cdot \frac{-1 + t\_1 \cdot \left(t\_3 \cdot \left(t\_0 \cdot t\_0\right)\right)}{1 + t\_1 \cdot \left(t\_1 - -1\right)}\\
\mathbf{elif}\;im\_m \leq 2.5 \cdot 10^{+24}:\\
\;\;\;\;e^{im\_m} \cdot -0.5 + \frac{0.5}{e^{im\_m}}\\
\mathbf{elif}\;im\_m \leq 4.6 \cdot 10^{+44}:\\
\;\;\;\;\frac{im\_m \cdot \left(\cos re \cdot \left(\left(-1 + t\_3 \cdot \left(\left(im\_m \cdot im\_m\right) \cdot -0.004629629629629629\right)\right) \cdot t\_2 + t\_3 \cdot \left(t\_4 \cdot \left(6.944444444444444 \cdot 10^{-5} + t\_3 \cdot -3.936759889140842 \cdot 10^{-8}\right)\right)\right)\right)}{t\_2 \cdot t\_4}\\
\mathbf{else}:\\
\;\;\;\;im\_m \cdot \left(\cos re \cdot \left(-1 + \left(im\_m \cdot im\_m\right) \cdot \left(-0.16666666666666666 + im\_m \cdot \left(im\_m \cdot \left(-0.008333333333333333 + \left(im\_m \cdot im\_m\right) \cdot -0.0001984126984126984\right)\right)\right)\right)\right)\\
\end{array}
\end{array}
\end{array}
if im < 1Initial program 46.4%
Taylor expanded in im around 0
*-lowering-*.f64N/A
distribute-rgt-inN/A
associate-+r+N/A
*-commutativeN/A
+-commutativeN/A
+-commutativeN/A
Simplified92.8%
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
Applied egg-rr92.8%
Taylor expanded in im around 0
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/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
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6492.8%
Simplified92.8%
flip3-+N/A
/-lowering-/.f64N/A
Applied egg-rr61.6%
if 1 < im < 2.50000000000000023e24Initial program 100.0%
Taylor expanded in re around 0
sub-negN/A
+-commutativeN/A
distribute-rgt-inN/A
*-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
distribute-lft-neg-outN/A
distribute-rgt-neg-inN/A
metadata-evalN/A
*-lowering-*.f64N/A
exp-lowering-exp.f64N/A
exp-negN/A
associate-*l/N/A
metadata-evalN/A
/-lowering-/.f64N/A
exp-lowering-exp.f64100.0%
Simplified100.0%
if 2.50000000000000023e24 < im < 4.60000000000000009e44Initial program 100.0%
Taylor expanded in im around 0
*-lowering-*.f64N/A
distribute-rgt-inN/A
associate-+r+N/A
*-commutativeN/A
+-commutativeN/A
+-commutativeN/A
Simplified7.0%
Applied egg-rr81.0%
Applied egg-rr100.0%
if 4.60000000000000009e44 < im Initial program 100.0%
Taylor expanded in im around 0
*-lowering-*.f64N/A
distribute-rgt-inN/A
associate-+r+N/A
*-commutativeN/A
+-commutativeN/A
+-commutativeN/A
Simplified100.0%
Taylor expanded in im around 0
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/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
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64100.0%
Simplified100.0%
Final simplification70.9%
im\_m = (fabs.f64 im)
im\_s = (copysign.f64 #s(literal 1 binary64) im)
(FPCore (im_s re im_m)
:precision binary64
(let* ((t_0
(+
-0.16666666666666666
(*
im_m
(*
im_m
(+
-0.008333333333333333
(* im_m (* im_m -0.0001984126984126984)))))))
(t_1 (- (* (* im_m im_m) -0.16666666666666666) -1.0))
(t_2 (* im_m (* im_m (* im_m im_m))))
(t_3 (* im_m (* im_m t_0)))
(t_4
(- -0.008333333333333333 (* (* im_m im_m) -0.0001984126984126984))))
(*
im_s
(if (<= im_m 2.05e+26)
(*
(* im_m (cos re))
(/
(+ -1.0 (* t_3 (* (* (* im_m im_m) (* im_m im_m)) (* t_0 t_0))))
(+ 1.0 (* t_3 (- t_3 -1.0)))))
(if (<= im_m 1.3e+62)
(*
im_m
(*
(cos re)
(/
(+
(* (+ -1.0 (* t_2 0.027777777777777776)) t_4)
(*
t_1
(* t_2 (- 6.944444444444444e-5 (* t_2 3.936759889140842e-8)))))
(* t_4 t_1))))
(*
im_m
(*
(cos re)
(+
-1.0
(*
(* im_m im_m)
(+
-0.16666666666666666
(* (* im_m im_m) -0.008333333333333333)))))))))))im\_m = fabs(im);
im\_s = copysign(1.0, im);
double code(double im_s, double re, double im_m) {
double t_0 = -0.16666666666666666 + (im_m * (im_m * (-0.008333333333333333 + (im_m * (im_m * -0.0001984126984126984)))));
double t_1 = ((im_m * im_m) * -0.16666666666666666) - -1.0;
double t_2 = im_m * (im_m * (im_m * im_m));
double t_3 = im_m * (im_m * t_0);
double t_4 = -0.008333333333333333 - ((im_m * im_m) * -0.0001984126984126984);
double tmp;
if (im_m <= 2.05e+26) {
tmp = (im_m * cos(re)) * ((-1.0 + (t_3 * (((im_m * im_m) * (im_m * im_m)) * (t_0 * t_0)))) / (1.0 + (t_3 * (t_3 - -1.0))));
} else if (im_m <= 1.3e+62) {
tmp = im_m * (cos(re) * ((((-1.0 + (t_2 * 0.027777777777777776)) * t_4) + (t_1 * (t_2 * (6.944444444444444e-5 - (t_2 * 3.936759889140842e-8))))) / (t_4 * t_1)));
} else {
tmp = im_m * (cos(re) * (-1.0 + ((im_m * im_m) * (-0.16666666666666666 + ((im_m * im_m) * -0.008333333333333333)))));
}
return im_s * tmp;
}
im\_m = abs(im)
im\_s = copysign(1.0d0, im)
real(8) function code(im_s, re, im_m)
real(8), intent (in) :: im_s
real(8), intent (in) :: re
real(8), intent (in) :: im_m
real(8) :: t_0
real(8) :: t_1
real(8) :: t_2
real(8) :: t_3
real(8) :: t_4
real(8) :: tmp
t_0 = (-0.16666666666666666d0) + (im_m * (im_m * ((-0.008333333333333333d0) + (im_m * (im_m * (-0.0001984126984126984d0))))))
t_1 = ((im_m * im_m) * (-0.16666666666666666d0)) - (-1.0d0)
t_2 = im_m * (im_m * (im_m * im_m))
t_3 = im_m * (im_m * t_0)
t_4 = (-0.008333333333333333d0) - ((im_m * im_m) * (-0.0001984126984126984d0))
if (im_m <= 2.05d+26) then
tmp = (im_m * cos(re)) * (((-1.0d0) + (t_3 * (((im_m * im_m) * (im_m * im_m)) * (t_0 * t_0)))) / (1.0d0 + (t_3 * (t_3 - (-1.0d0)))))
else if (im_m <= 1.3d+62) then
tmp = im_m * (cos(re) * (((((-1.0d0) + (t_2 * 0.027777777777777776d0)) * t_4) + (t_1 * (t_2 * (6.944444444444444d-5 - (t_2 * 3.936759889140842d-8))))) / (t_4 * t_1)))
else
tmp = im_m * (cos(re) * ((-1.0d0) + ((im_m * im_m) * ((-0.16666666666666666d0) + ((im_m * im_m) * (-0.008333333333333333d0))))))
end if
code = im_s * tmp
end function
im\_m = Math.abs(im);
im\_s = Math.copySign(1.0, im);
public static double code(double im_s, double re, double im_m) {
double t_0 = -0.16666666666666666 + (im_m * (im_m * (-0.008333333333333333 + (im_m * (im_m * -0.0001984126984126984)))));
double t_1 = ((im_m * im_m) * -0.16666666666666666) - -1.0;
double t_2 = im_m * (im_m * (im_m * im_m));
double t_3 = im_m * (im_m * t_0);
double t_4 = -0.008333333333333333 - ((im_m * im_m) * -0.0001984126984126984);
double tmp;
if (im_m <= 2.05e+26) {
tmp = (im_m * Math.cos(re)) * ((-1.0 + (t_3 * (((im_m * im_m) * (im_m * im_m)) * (t_0 * t_0)))) / (1.0 + (t_3 * (t_3 - -1.0))));
} else if (im_m <= 1.3e+62) {
tmp = im_m * (Math.cos(re) * ((((-1.0 + (t_2 * 0.027777777777777776)) * t_4) + (t_1 * (t_2 * (6.944444444444444e-5 - (t_2 * 3.936759889140842e-8))))) / (t_4 * t_1)));
} else {
tmp = im_m * (Math.cos(re) * (-1.0 + ((im_m * im_m) * (-0.16666666666666666 + ((im_m * im_m) * -0.008333333333333333)))));
}
return im_s * tmp;
}
im\_m = math.fabs(im) im\_s = math.copysign(1.0, im) def code(im_s, re, im_m): t_0 = -0.16666666666666666 + (im_m * (im_m * (-0.008333333333333333 + (im_m * (im_m * -0.0001984126984126984))))) t_1 = ((im_m * im_m) * -0.16666666666666666) - -1.0 t_2 = im_m * (im_m * (im_m * im_m)) t_3 = im_m * (im_m * t_0) t_4 = -0.008333333333333333 - ((im_m * im_m) * -0.0001984126984126984) tmp = 0 if im_m <= 2.05e+26: tmp = (im_m * math.cos(re)) * ((-1.0 + (t_3 * (((im_m * im_m) * (im_m * im_m)) * (t_0 * t_0)))) / (1.0 + (t_3 * (t_3 - -1.0)))) elif im_m <= 1.3e+62: tmp = im_m * (math.cos(re) * ((((-1.0 + (t_2 * 0.027777777777777776)) * t_4) + (t_1 * (t_2 * (6.944444444444444e-5 - (t_2 * 3.936759889140842e-8))))) / (t_4 * t_1))) else: tmp = im_m * (math.cos(re) * (-1.0 + ((im_m * im_m) * (-0.16666666666666666 + ((im_m * im_m) * -0.008333333333333333))))) return im_s * tmp
im\_m = abs(im) im\_s = copysign(1.0, im) function code(im_s, re, im_m) t_0 = Float64(-0.16666666666666666 + Float64(im_m * Float64(im_m * Float64(-0.008333333333333333 + Float64(im_m * Float64(im_m * -0.0001984126984126984)))))) t_1 = Float64(Float64(Float64(im_m * im_m) * -0.16666666666666666) - -1.0) t_2 = Float64(im_m * Float64(im_m * Float64(im_m * im_m))) t_3 = Float64(im_m * Float64(im_m * t_0)) t_4 = Float64(-0.008333333333333333 - Float64(Float64(im_m * im_m) * -0.0001984126984126984)) tmp = 0.0 if (im_m <= 2.05e+26) tmp = Float64(Float64(im_m * cos(re)) * Float64(Float64(-1.0 + Float64(t_3 * Float64(Float64(Float64(im_m * im_m) * Float64(im_m * im_m)) * Float64(t_0 * t_0)))) / Float64(1.0 + Float64(t_3 * Float64(t_3 - -1.0))))); elseif (im_m <= 1.3e+62) tmp = Float64(im_m * Float64(cos(re) * Float64(Float64(Float64(Float64(-1.0 + Float64(t_2 * 0.027777777777777776)) * t_4) + Float64(t_1 * Float64(t_2 * Float64(6.944444444444444e-5 - Float64(t_2 * 3.936759889140842e-8))))) / Float64(t_4 * t_1)))); else tmp = Float64(im_m * Float64(cos(re) * Float64(-1.0 + Float64(Float64(im_m * im_m) * Float64(-0.16666666666666666 + Float64(Float64(im_m * im_m) * -0.008333333333333333)))))); end return Float64(im_s * tmp) end
im\_m = abs(im); im\_s = sign(im) * abs(1.0); function tmp_2 = code(im_s, re, im_m) t_0 = -0.16666666666666666 + (im_m * (im_m * (-0.008333333333333333 + (im_m * (im_m * -0.0001984126984126984))))); t_1 = ((im_m * im_m) * -0.16666666666666666) - -1.0; t_2 = im_m * (im_m * (im_m * im_m)); t_3 = im_m * (im_m * t_0); t_4 = -0.008333333333333333 - ((im_m * im_m) * -0.0001984126984126984); tmp = 0.0; if (im_m <= 2.05e+26) tmp = (im_m * cos(re)) * ((-1.0 + (t_3 * (((im_m * im_m) * (im_m * im_m)) * (t_0 * t_0)))) / (1.0 + (t_3 * (t_3 - -1.0)))); elseif (im_m <= 1.3e+62) tmp = im_m * (cos(re) * ((((-1.0 + (t_2 * 0.027777777777777776)) * t_4) + (t_1 * (t_2 * (6.944444444444444e-5 - (t_2 * 3.936759889140842e-8))))) / (t_4 * t_1))); else tmp = im_m * (cos(re) * (-1.0 + ((im_m * im_m) * (-0.16666666666666666 + ((im_m * im_m) * -0.008333333333333333))))); end tmp_2 = im_s * tmp; end
im\_m = N[Abs[im], $MachinePrecision]
im\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[im]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[im$95$s_, re_, im$95$m_] := Block[{t$95$0 = N[(-0.16666666666666666 + N[(im$95$m * N[(im$95$m * N[(-0.008333333333333333 + N[(im$95$m * N[(im$95$m * -0.0001984126984126984), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[(N[(im$95$m * im$95$m), $MachinePrecision] * -0.16666666666666666), $MachinePrecision] - -1.0), $MachinePrecision]}, Block[{t$95$2 = N[(im$95$m * N[(im$95$m * N[(im$95$m * im$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = N[(im$95$m * N[(im$95$m * t$95$0), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$4 = N[(-0.008333333333333333 - N[(N[(im$95$m * im$95$m), $MachinePrecision] * -0.0001984126984126984), $MachinePrecision]), $MachinePrecision]}, N[(im$95$s * If[LessEqual[im$95$m, 2.05e+26], N[(N[(im$95$m * N[Cos[re], $MachinePrecision]), $MachinePrecision] * N[(N[(-1.0 + N[(t$95$3 * N[(N[(N[(im$95$m * im$95$m), $MachinePrecision] * N[(im$95$m * im$95$m), $MachinePrecision]), $MachinePrecision] * N[(t$95$0 * t$95$0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(1.0 + N[(t$95$3 * N[(t$95$3 - -1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[im$95$m, 1.3e+62], N[(im$95$m * N[(N[Cos[re], $MachinePrecision] * N[(N[(N[(N[(-1.0 + N[(t$95$2 * 0.027777777777777776), $MachinePrecision]), $MachinePrecision] * t$95$4), $MachinePrecision] + N[(t$95$1 * N[(t$95$2 * N[(6.944444444444444e-5 - N[(t$95$2 * 3.936759889140842e-8), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(t$95$4 * t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(im$95$m * N[(N[Cos[re], $MachinePrecision] * N[(-1.0 + N[(N[(im$95$m * im$95$m), $MachinePrecision] * N[(-0.16666666666666666 + N[(N[(im$95$m * im$95$m), $MachinePrecision] * -0.008333333333333333), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]]]]]]
\begin{array}{l}
im\_m = \left|im\right|
\\
im\_s = \mathsf{copysign}\left(1, im\right)
\\
\begin{array}{l}
t_0 := -0.16666666666666666 + im\_m \cdot \left(im\_m \cdot \left(-0.008333333333333333 + im\_m \cdot \left(im\_m \cdot -0.0001984126984126984\right)\right)\right)\\
t_1 := \left(im\_m \cdot im\_m\right) \cdot -0.16666666666666666 - -1\\
t_2 := im\_m \cdot \left(im\_m \cdot \left(im\_m \cdot im\_m\right)\right)\\
t_3 := im\_m \cdot \left(im\_m \cdot t\_0\right)\\
t_4 := -0.008333333333333333 - \left(im\_m \cdot im\_m\right) \cdot -0.0001984126984126984\\
im\_s \cdot \begin{array}{l}
\mathbf{if}\;im\_m \leq 2.05 \cdot 10^{+26}:\\
\;\;\;\;\left(im\_m \cdot \cos re\right) \cdot \frac{-1 + t\_3 \cdot \left(\left(\left(im\_m \cdot im\_m\right) \cdot \left(im\_m \cdot im\_m\right)\right) \cdot \left(t\_0 \cdot t\_0\right)\right)}{1 + t\_3 \cdot \left(t\_3 - -1\right)}\\
\mathbf{elif}\;im\_m \leq 1.3 \cdot 10^{+62}:\\
\;\;\;\;im\_m \cdot \left(\cos re \cdot \frac{\left(-1 + t\_2 \cdot 0.027777777777777776\right) \cdot t\_4 + t\_1 \cdot \left(t\_2 \cdot \left(6.944444444444444 \cdot 10^{-5} - t\_2 \cdot 3.936759889140842 \cdot 10^{-8}\right)\right)}{t\_4 \cdot t\_1}\right)\\
\mathbf{else}:\\
\;\;\;\;im\_m \cdot \left(\cos re \cdot \left(-1 + \left(im\_m \cdot im\_m\right) \cdot \left(-0.16666666666666666 + \left(im\_m \cdot im\_m\right) \cdot -0.008333333333333333\right)\right)\right)\\
\end{array}
\end{array}
\end{array}
if im < 2.04999999999999992e26Initial program 48.8%
Taylor expanded in im around 0
*-lowering-*.f64N/A
distribute-rgt-inN/A
associate-+r+N/A
*-commutativeN/A
+-commutativeN/A
+-commutativeN/A
Simplified88.9%
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
Applied egg-rr88.9%
Taylor expanded in im around 0
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/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
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6488.9%
Simplified88.9%
flip3-+N/A
/-lowering-/.f64N/A
Applied egg-rr62.0%
if 2.04999999999999992e26 < im < 1.29999999999999992e62Initial program 100.0%
Taylor expanded in im around 0
*-lowering-*.f64N/A
distribute-rgt-inN/A
associate-+r+N/A
*-commutativeN/A
+-commutativeN/A
+-commutativeN/A
Simplified58.9%
Applied egg-rr100.0%
if 1.29999999999999992e62 < im Initial program 100.0%
Taylor expanded in im around 0
*-lowering-*.f64N/A
*-commutativeN/A
+-commutativeN/A
distribute-lft-inN/A
associate-*r*N/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
distribute-rgt-outN/A
distribute-lft-outN/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
Simplified100.0%
Final simplification69.8%
im\_m = (fabs.f64 im)
im\_s = (copysign.f64 #s(literal 1 binary64) im)
(FPCore (im_s re im_m)
:precision binary64
(let* ((t_0
(+
-1.0
(*
(* im_m im_m)
(+
-0.16666666666666666
(*
im_m
(*
im_m
(+
-0.008333333333333333
(* (* im_m im_m) -0.0001984126984126984)))))))))
(*
im_s
(if (<= im_m 5.2)
(* (* im_m (cos re)) t_0)
(if (<= im_m 4e+44)
(*
(-
(/
1.0
(+
1.0
(* im_m (+ 1.0 (* im_m (+ 0.5 (* im_m 0.16666666666666666)))))))
(exp im_m))
(+ 0.5 (* -0.25 (* re re))))
(* im_m (* (cos re) t_0)))))))im\_m = fabs(im);
im\_s = copysign(1.0, im);
double code(double im_s, double re, double im_m) {
double t_0 = -1.0 + ((im_m * im_m) * (-0.16666666666666666 + (im_m * (im_m * (-0.008333333333333333 + ((im_m * im_m) * -0.0001984126984126984))))));
double tmp;
if (im_m <= 5.2) {
tmp = (im_m * cos(re)) * t_0;
} else if (im_m <= 4e+44) {
tmp = ((1.0 / (1.0 + (im_m * (1.0 + (im_m * (0.5 + (im_m * 0.16666666666666666))))))) - exp(im_m)) * (0.5 + (-0.25 * (re * re)));
} else {
tmp = im_m * (cos(re) * t_0);
}
return im_s * tmp;
}
im\_m = abs(im)
im\_s = copysign(1.0d0, im)
real(8) function code(im_s, re, im_m)
real(8), intent (in) :: im_s
real(8), intent (in) :: re
real(8), intent (in) :: im_m
real(8) :: t_0
real(8) :: tmp
t_0 = (-1.0d0) + ((im_m * im_m) * ((-0.16666666666666666d0) + (im_m * (im_m * ((-0.008333333333333333d0) + ((im_m * im_m) * (-0.0001984126984126984d0)))))))
if (im_m <= 5.2d0) then
tmp = (im_m * cos(re)) * t_0
else if (im_m <= 4d+44) then
tmp = ((1.0d0 / (1.0d0 + (im_m * (1.0d0 + (im_m * (0.5d0 + (im_m * 0.16666666666666666d0))))))) - exp(im_m)) * (0.5d0 + ((-0.25d0) * (re * re)))
else
tmp = im_m * (cos(re) * t_0)
end if
code = im_s * tmp
end function
im\_m = Math.abs(im);
im\_s = Math.copySign(1.0, im);
public static double code(double im_s, double re, double im_m) {
double t_0 = -1.0 + ((im_m * im_m) * (-0.16666666666666666 + (im_m * (im_m * (-0.008333333333333333 + ((im_m * im_m) * -0.0001984126984126984))))));
double tmp;
if (im_m <= 5.2) {
tmp = (im_m * Math.cos(re)) * t_0;
} else if (im_m <= 4e+44) {
tmp = ((1.0 / (1.0 + (im_m * (1.0 + (im_m * (0.5 + (im_m * 0.16666666666666666))))))) - Math.exp(im_m)) * (0.5 + (-0.25 * (re * re)));
} else {
tmp = im_m * (Math.cos(re) * t_0);
}
return im_s * tmp;
}
im\_m = math.fabs(im) im\_s = math.copysign(1.0, im) def code(im_s, re, im_m): t_0 = -1.0 + ((im_m * im_m) * (-0.16666666666666666 + (im_m * (im_m * (-0.008333333333333333 + ((im_m * im_m) * -0.0001984126984126984)))))) tmp = 0 if im_m <= 5.2: tmp = (im_m * math.cos(re)) * t_0 elif im_m <= 4e+44: tmp = ((1.0 / (1.0 + (im_m * (1.0 + (im_m * (0.5 + (im_m * 0.16666666666666666))))))) - math.exp(im_m)) * (0.5 + (-0.25 * (re * re))) else: tmp = im_m * (math.cos(re) * t_0) return im_s * tmp
im\_m = abs(im) im\_s = copysign(1.0, im) function code(im_s, re, im_m) t_0 = Float64(-1.0 + Float64(Float64(im_m * im_m) * Float64(-0.16666666666666666 + Float64(im_m * Float64(im_m * Float64(-0.008333333333333333 + Float64(Float64(im_m * im_m) * -0.0001984126984126984))))))) tmp = 0.0 if (im_m <= 5.2) tmp = Float64(Float64(im_m * cos(re)) * t_0); elseif (im_m <= 4e+44) tmp = Float64(Float64(Float64(1.0 / Float64(1.0 + Float64(im_m * Float64(1.0 + Float64(im_m * Float64(0.5 + Float64(im_m * 0.16666666666666666))))))) - exp(im_m)) * Float64(0.5 + Float64(-0.25 * Float64(re * re)))); else tmp = Float64(im_m * Float64(cos(re) * t_0)); end return Float64(im_s * tmp) end
im\_m = abs(im); im\_s = sign(im) * abs(1.0); function tmp_2 = code(im_s, re, im_m) t_0 = -1.0 + ((im_m * im_m) * (-0.16666666666666666 + (im_m * (im_m * (-0.008333333333333333 + ((im_m * im_m) * -0.0001984126984126984)))))); tmp = 0.0; if (im_m <= 5.2) tmp = (im_m * cos(re)) * t_0; elseif (im_m <= 4e+44) tmp = ((1.0 / (1.0 + (im_m * (1.0 + (im_m * (0.5 + (im_m * 0.16666666666666666))))))) - exp(im_m)) * (0.5 + (-0.25 * (re * re))); else tmp = im_m * (cos(re) * t_0); end tmp_2 = im_s * tmp; end
im\_m = N[Abs[im], $MachinePrecision]
im\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[im]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[im$95$s_, re_, im$95$m_] := Block[{t$95$0 = N[(-1.0 + N[(N[(im$95$m * im$95$m), $MachinePrecision] * N[(-0.16666666666666666 + N[(im$95$m * N[(im$95$m * N[(-0.008333333333333333 + N[(N[(im$95$m * im$95$m), $MachinePrecision] * -0.0001984126984126984), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, N[(im$95$s * If[LessEqual[im$95$m, 5.2], N[(N[(im$95$m * N[Cos[re], $MachinePrecision]), $MachinePrecision] * t$95$0), $MachinePrecision], If[LessEqual[im$95$m, 4e+44], N[(N[(N[(1.0 / N[(1.0 + N[(im$95$m * N[(1.0 + N[(im$95$m * N[(0.5 + N[(im$95$m * 0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[Exp[im$95$m], $MachinePrecision]), $MachinePrecision] * N[(0.5 + N[(-0.25 * N[(re * re), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(im$95$m * N[(N[Cos[re], $MachinePrecision] * t$95$0), $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]]
\begin{array}{l}
im\_m = \left|im\right|
\\
im\_s = \mathsf{copysign}\left(1, im\right)
\\
\begin{array}{l}
t_0 := -1 + \left(im\_m \cdot im\_m\right) \cdot \left(-0.16666666666666666 + im\_m \cdot \left(im\_m \cdot \left(-0.008333333333333333 + \left(im\_m \cdot im\_m\right) \cdot -0.0001984126984126984\right)\right)\right)\\
im\_s \cdot \begin{array}{l}
\mathbf{if}\;im\_m \leq 5.2:\\
\;\;\;\;\left(im\_m \cdot \cos re\right) \cdot t\_0\\
\mathbf{elif}\;im\_m \leq 4 \cdot 10^{+44}:\\
\;\;\;\;\left(\frac{1}{1 + im\_m \cdot \left(1 + im\_m \cdot \left(0.5 + im\_m \cdot 0.16666666666666666\right)\right)} - e^{im\_m}\right) \cdot \left(0.5 + -0.25 \cdot \left(re \cdot re\right)\right)\\
\mathbf{else}:\\
\;\;\;\;im\_m \cdot \left(\cos re \cdot t\_0\right)\\
\end{array}
\end{array}
\end{array}
if im < 5.20000000000000018Initial program 46.4%
Taylor expanded in im around 0
*-lowering-*.f64N/A
distribute-rgt-inN/A
associate-+r+N/A
*-commutativeN/A
+-commutativeN/A
+-commutativeN/A
Simplified92.8%
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
Applied egg-rr92.8%
Taylor expanded in im around 0
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/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
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6492.8%
Simplified92.8%
if 5.20000000000000018 < im < 4.0000000000000004e44Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
associate-*r*N/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
--lowering--.f64N/A
exp-negN/A
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
exp-lowering-exp.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6461.5%
Simplified61.5%
Taylor expanded in im around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6461.5%
Simplified61.5%
if 4.0000000000000004e44 < im Initial program 100.0%
Taylor expanded in im around 0
*-lowering-*.f64N/A
distribute-rgt-inN/A
associate-+r+N/A
*-commutativeN/A
+-commutativeN/A
+-commutativeN/A
Simplified100.0%
Taylor expanded in im around 0
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/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
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64100.0%
Simplified100.0%
Final simplification92.6%
im\_m = (fabs.f64 im)
im\_s = (copysign.f64 #s(literal 1 binary64) im)
(FPCore (im_s re im_m)
:precision binary64
(let* ((t_0
(+
-1.0
(*
(* im_m im_m)
(+
-0.16666666666666666
(*
im_m
(*
im_m
(+
-0.008333333333333333
(* (* im_m im_m) -0.0001984126984126984)))))))))
(*
im_s
(if (<= im_m 5.4)
(* (* im_m (cos re)) t_0)
(if (<= im_m 3.8e+44)
(*
(+ 0.5 (* -0.25 (* re re)))
(- (/ 1.0 (+ 1.0 (* im_m (+ 1.0 (* im_m 0.5))))) (exp im_m)))
(* im_m (* (cos re) t_0)))))))im\_m = fabs(im);
im\_s = copysign(1.0, im);
double code(double im_s, double re, double im_m) {
double t_0 = -1.0 + ((im_m * im_m) * (-0.16666666666666666 + (im_m * (im_m * (-0.008333333333333333 + ((im_m * im_m) * -0.0001984126984126984))))));
double tmp;
if (im_m <= 5.4) {
tmp = (im_m * cos(re)) * t_0;
} else if (im_m <= 3.8e+44) {
tmp = (0.5 + (-0.25 * (re * re))) * ((1.0 / (1.0 + (im_m * (1.0 + (im_m * 0.5))))) - exp(im_m));
} else {
tmp = im_m * (cos(re) * t_0);
}
return im_s * tmp;
}
im\_m = abs(im)
im\_s = copysign(1.0d0, im)
real(8) function code(im_s, re, im_m)
real(8), intent (in) :: im_s
real(8), intent (in) :: re
real(8), intent (in) :: im_m
real(8) :: t_0
real(8) :: tmp
t_0 = (-1.0d0) + ((im_m * im_m) * ((-0.16666666666666666d0) + (im_m * (im_m * ((-0.008333333333333333d0) + ((im_m * im_m) * (-0.0001984126984126984d0)))))))
if (im_m <= 5.4d0) then
tmp = (im_m * cos(re)) * t_0
else if (im_m <= 3.8d+44) then
tmp = (0.5d0 + ((-0.25d0) * (re * re))) * ((1.0d0 / (1.0d0 + (im_m * (1.0d0 + (im_m * 0.5d0))))) - exp(im_m))
else
tmp = im_m * (cos(re) * t_0)
end if
code = im_s * tmp
end function
im\_m = Math.abs(im);
im\_s = Math.copySign(1.0, im);
public static double code(double im_s, double re, double im_m) {
double t_0 = -1.0 + ((im_m * im_m) * (-0.16666666666666666 + (im_m * (im_m * (-0.008333333333333333 + ((im_m * im_m) * -0.0001984126984126984))))));
double tmp;
if (im_m <= 5.4) {
tmp = (im_m * Math.cos(re)) * t_0;
} else if (im_m <= 3.8e+44) {
tmp = (0.5 + (-0.25 * (re * re))) * ((1.0 / (1.0 + (im_m * (1.0 + (im_m * 0.5))))) - Math.exp(im_m));
} else {
tmp = im_m * (Math.cos(re) * t_0);
}
return im_s * tmp;
}
im\_m = math.fabs(im) im\_s = math.copysign(1.0, im) def code(im_s, re, im_m): t_0 = -1.0 + ((im_m * im_m) * (-0.16666666666666666 + (im_m * (im_m * (-0.008333333333333333 + ((im_m * im_m) * -0.0001984126984126984)))))) tmp = 0 if im_m <= 5.4: tmp = (im_m * math.cos(re)) * t_0 elif im_m <= 3.8e+44: tmp = (0.5 + (-0.25 * (re * re))) * ((1.0 / (1.0 + (im_m * (1.0 + (im_m * 0.5))))) - math.exp(im_m)) else: tmp = im_m * (math.cos(re) * t_0) return im_s * tmp
im\_m = abs(im) im\_s = copysign(1.0, im) function code(im_s, re, im_m) t_0 = Float64(-1.0 + Float64(Float64(im_m * im_m) * Float64(-0.16666666666666666 + Float64(im_m * Float64(im_m * Float64(-0.008333333333333333 + Float64(Float64(im_m * im_m) * -0.0001984126984126984))))))) tmp = 0.0 if (im_m <= 5.4) tmp = Float64(Float64(im_m * cos(re)) * t_0); elseif (im_m <= 3.8e+44) tmp = Float64(Float64(0.5 + Float64(-0.25 * Float64(re * re))) * Float64(Float64(1.0 / Float64(1.0 + Float64(im_m * Float64(1.0 + Float64(im_m * 0.5))))) - exp(im_m))); else tmp = Float64(im_m * Float64(cos(re) * t_0)); end return Float64(im_s * tmp) end
im\_m = abs(im); im\_s = sign(im) * abs(1.0); function tmp_2 = code(im_s, re, im_m) t_0 = -1.0 + ((im_m * im_m) * (-0.16666666666666666 + (im_m * (im_m * (-0.008333333333333333 + ((im_m * im_m) * -0.0001984126984126984)))))); tmp = 0.0; if (im_m <= 5.4) tmp = (im_m * cos(re)) * t_0; elseif (im_m <= 3.8e+44) tmp = (0.5 + (-0.25 * (re * re))) * ((1.0 / (1.0 + (im_m * (1.0 + (im_m * 0.5))))) - exp(im_m)); else tmp = im_m * (cos(re) * t_0); end tmp_2 = im_s * tmp; end
im\_m = N[Abs[im], $MachinePrecision]
im\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[im]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[im$95$s_, re_, im$95$m_] := Block[{t$95$0 = N[(-1.0 + N[(N[(im$95$m * im$95$m), $MachinePrecision] * N[(-0.16666666666666666 + N[(im$95$m * N[(im$95$m * N[(-0.008333333333333333 + N[(N[(im$95$m * im$95$m), $MachinePrecision] * -0.0001984126984126984), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, N[(im$95$s * If[LessEqual[im$95$m, 5.4], N[(N[(im$95$m * N[Cos[re], $MachinePrecision]), $MachinePrecision] * t$95$0), $MachinePrecision], If[LessEqual[im$95$m, 3.8e+44], N[(N[(0.5 + N[(-0.25 * N[(re * re), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[(1.0 / N[(1.0 + N[(im$95$m * N[(1.0 + N[(im$95$m * 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[Exp[im$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(im$95$m * N[(N[Cos[re], $MachinePrecision] * t$95$0), $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]]
\begin{array}{l}
im\_m = \left|im\right|
\\
im\_s = \mathsf{copysign}\left(1, im\right)
\\
\begin{array}{l}
t_0 := -1 + \left(im\_m \cdot im\_m\right) \cdot \left(-0.16666666666666666 + im\_m \cdot \left(im\_m \cdot \left(-0.008333333333333333 + \left(im\_m \cdot im\_m\right) \cdot -0.0001984126984126984\right)\right)\right)\\
im\_s \cdot \begin{array}{l}
\mathbf{if}\;im\_m \leq 5.4:\\
\;\;\;\;\left(im\_m \cdot \cos re\right) \cdot t\_0\\
\mathbf{elif}\;im\_m \leq 3.8 \cdot 10^{+44}:\\
\;\;\;\;\left(0.5 + -0.25 \cdot \left(re \cdot re\right)\right) \cdot \left(\frac{1}{1 + im\_m \cdot \left(1 + im\_m \cdot 0.5\right)} - e^{im\_m}\right)\\
\mathbf{else}:\\
\;\;\;\;im\_m \cdot \left(\cos re \cdot t\_0\right)\\
\end{array}
\end{array}
\end{array}
if im < 5.4000000000000004Initial program 46.4%
Taylor expanded in im around 0
*-lowering-*.f64N/A
distribute-rgt-inN/A
associate-+r+N/A
*-commutativeN/A
+-commutativeN/A
+-commutativeN/A
Simplified92.8%
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
Applied egg-rr92.8%
Taylor expanded in im around 0
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/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
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6492.8%
Simplified92.8%
if 5.4000000000000004 < im < 3.8000000000000002e44Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
associate-*r*N/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
--lowering--.f64N/A
exp-negN/A
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
exp-lowering-exp.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6461.5%
Simplified61.5%
Taylor expanded in im around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6461.2%
Simplified61.2%
if 3.8000000000000002e44 < im Initial program 100.0%
Taylor expanded in im around 0
*-lowering-*.f64N/A
distribute-rgt-inN/A
associate-+r+N/A
*-commutativeN/A
+-commutativeN/A
+-commutativeN/A
Simplified100.0%
Taylor expanded in im around 0
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/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
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64100.0%
Simplified100.0%
Final simplification92.6%
im\_m = (fabs.f64 im)
im\_s = (copysign.f64 #s(literal 1 binary64) im)
(FPCore (im_s re im_m)
:precision binary64
(let* ((t_0
(*
im_m
(*
(cos re)
(+
-1.0
(*
(* im_m im_m)
(+
-0.16666666666666666
(*
im_m
(*
im_m
(+
-0.008333333333333333
(* (* im_m im_m) -0.0001984126984126984)))))))))))
(*
im_s
(if (<= im_m 5.2)
t_0
(if (<= im_m 3.6e+44)
(*
(+ 0.5 (* -0.25 (* re re)))
(- (/ 1.0 (+ 1.0 (* im_m (+ 1.0 (* im_m 0.5))))) (exp im_m)))
t_0)))))im\_m = fabs(im);
im\_s = copysign(1.0, im);
double code(double im_s, double re, double im_m) {
double t_0 = im_m * (cos(re) * (-1.0 + ((im_m * im_m) * (-0.16666666666666666 + (im_m * (im_m * (-0.008333333333333333 + ((im_m * im_m) * -0.0001984126984126984))))))));
double tmp;
if (im_m <= 5.2) {
tmp = t_0;
} else if (im_m <= 3.6e+44) {
tmp = (0.5 + (-0.25 * (re * re))) * ((1.0 / (1.0 + (im_m * (1.0 + (im_m * 0.5))))) - exp(im_m));
} else {
tmp = t_0;
}
return im_s * tmp;
}
im\_m = abs(im)
im\_s = copysign(1.0d0, im)
real(8) function code(im_s, re, im_m)
real(8), intent (in) :: im_s
real(8), intent (in) :: re
real(8), intent (in) :: im_m
real(8) :: t_0
real(8) :: tmp
t_0 = im_m * (cos(re) * ((-1.0d0) + ((im_m * im_m) * ((-0.16666666666666666d0) + (im_m * (im_m * ((-0.008333333333333333d0) + ((im_m * im_m) * (-0.0001984126984126984d0)))))))))
if (im_m <= 5.2d0) then
tmp = t_0
else if (im_m <= 3.6d+44) then
tmp = (0.5d0 + ((-0.25d0) * (re * re))) * ((1.0d0 / (1.0d0 + (im_m * (1.0d0 + (im_m * 0.5d0))))) - exp(im_m))
else
tmp = t_0
end if
code = im_s * tmp
end function
im\_m = Math.abs(im);
im\_s = Math.copySign(1.0, im);
public static double code(double im_s, double re, double im_m) {
double t_0 = im_m * (Math.cos(re) * (-1.0 + ((im_m * im_m) * (-0.16666666666666666 + (im_m * (im_m * (-0.008333333333333333 + ((im_m * im_m) * -0.0001984126984126984))))))));
double tmp;
if (im_m <= 5.2) {
tmp = t_0;
} else if (im_m <= 3.6e+44) {
tmp = (0.5 + (-0.25 * (re * re))) * ((1.0 / (1.0 + (im_m * (1.0 + (im_m * 0.5))))) - Math.exp(im_m));
} else {
tmp = t_0;
}
return im_s * tmp;
}
im\_m = math.fabs(im) im\_s = math.copysign(1.0, im) def code(im_s, re, im_m): t_0 = im_m * (math.cos(re) * (-1.0 + ((im_m * im_m) * (-0.16666666666666666 + (im_m * (im_m * (-0.008333333333333333 + ((im_m * im_m) * -0.0001984126984126984)))))))) tmp = 0 if im_m <= 5.2: tmp = t_0 elif im_m <= 3.6e+44: tmp = (0.5 + (-0.25 * (re * re))) * ((1.0 / (1.0 + (im_m * (1.0 + (im_m * 0.5))))) - math.exp(im_m)) else: tmp = t_0 return im_s * tmp
im\_m = abs(im) im\_s = copysign(1.0, im) function code(im_s, re, im_m) t_0 = Float64(im_m * Float64(cos(re) * Float64(-1.0 + Float64(Float64(im_m * im_m) * Float64(-0.16666666666666666 + Float64(im_m * Float64(im_m * Float64(-0.008333333333333333 + Float64(Float64(im_m * im_m) * -0.0001984126984126984))))))))) tmp = 0.0 if (im_m <= 5.2) tmp = t_0; elseif (im_m <= 3.6e+44) tmp = Float64(Float64(0.5 + Float64(-0.25 * Float64(re * re))) * Float64(Float64(1.0 / Float64(1.0 + Float64(im_m * Float64(1.0 + Float64(im_m * 0.5))))) - exp(im_m))); else tmp = t_0; end return Float64(im_s * tmp) end
im\_m = abs(im); im\_s = sign(im) * abs(1.0); function tmp_2 = code(im_s, re, im_m) t_0 = im_m * (cos(re) * (-1.0 + ((im_m * im_m) * (-0.16666666666666666 + (im_m * (im_m * (-0.008333333333333333 + ((im_m * im_m) * -0.0001984126984126984)))))))); tmp = 0.0; if (im_m <= 5.2) tmp = t_0; elseif (im_m <= 3.6e+44) tmp = (0.5 + (-0.25 * (re * re))) * ((1.0 / (1.0 + (im_m * (1.0 + (im_m * 0.5))))) - exp(im_m)); else tmp = t_0; end tmp_2 = im_s * tmp; end
im\_m = N[Abs[im], $MachinePrecision]
im\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[im]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[im$95$s_, re_, im$95$m_] := Block[{t$95$0 = N[(im$95$m * N[(N[Cos[re], $MachinePrecision] * N[(-1.0 + N[(N[(im$95$m * im$95$m), $MachinePrecision] * N[(-0.16666666666666666 + N[(im$95$m * N[(im$95$m * N[(-0.008333333333333333 + N[(N[(im$95$m * im$95$m), $MachinePrecision] * -0.0001984126984126984), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, N[(im$95$s * If[LessEqual[im$95$m, 5.2], t$95$0, If[LessEqual[im$95$m, 3.6e+44], N[(N[(0.5 + N[(-0.25 * N[(re * re), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[(1.0 / N[(1.0 + N[(im$95$m * N[(1.0 + N[(im$95$m * 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[Exp[im$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$0]]), $MachinePrecision]]
\begin{array}{l}
im\_m = \left|im\right|
\\
im\_s = \mathsf{copysign}\left(1, im\right)
\\
\begin{array}{l}
t_0 := im\_m \cdot \left(\cos re \cdot \left(-1 + \left(im\_m \cdot im\_m\right) \cdot \left(-0.16666666666666666 + im\_m \cdot \left(im\_m \cdot \left(-0.008333333333333333 + \left(im\_m \cdot im\_m\right) \cdot -0.0001984126984126984\right)\right)\right)\right)\right)\\
im\_s \cdot \begin{array}{l}
\mathbf{if}\;im\_m \leq 5.2:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;im\_m \leq 3.6 \cdot 10^{+44}:\\
\;\;\;\;\left(0.5 + -0.25 \cdot \left(re \cdot re\right)\right) \cdot \left(\frac{1}{1 + im\_m \cdot \left(1 + im\_m \cdot 0.5\right)} - e^{im\_m}\right)\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
\end{array}
if im < 5.20000000000000018 or 3.6e44 < im Initial program 57.2%
Taylor expanded in im around 0
*-lowering-*.f64N/A
distribute-rgt-inN/A
associate-+r+N/A
*-commutativeN/A
+-commutativeN/A
+-commutativeN/A
Simplified94.3%
Taylor expanded in im around 0
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/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
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6494.3%
Simplified94.3%
if 5.20000000000000018 < im < 3.6e44Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
associate-*r*N/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
--lowering--.f64N/A
exp-negN/A
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
exp-lowering-exp.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6461.5%
Simplified61.5%
Taylor expanded in im around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6461.2%
Simplified61.2%
Final simplification92.6%
im\_m = (fabs.f64 im)
im\_s = (copysign.f64 #s(literal 1 binary64) im)
(FPCore (im_s re im_m)
:precision binary64
(let* ((t_0
(+
-1.0
(*
(* im_m im_m)
(+ -0.16666666666666666 (* (* im_m im_m) -0.008333333333333333))))))
(*
im_s
(if (<= im_m 5.8)
(* (* im_m (cos re)) t_0)
(if (<= im_m 9e+59)
(*
(+ 0.5 (* -0.25 (* re re)))
(- (/ 1.0 (+ 1.0 (* im_m (+ 1.0 (* im_m 0.5))))) (exp im_m)))
(* im_m (* (cos re) t_0)))))))im\_m = fabs(im);
im\_s = copysign(1.0, im);
double code(double im_s, double re, double im_m) {
double t_0 = -1.0 + ((im_m * im_m) * (-0.16666666666666666 + ((im_m * im_m) * -0.008333333333333333)));
double tmp;
if (im_m <= 5.8) {
tmp = (im_m * cos(re)) * t_0;
} else if (im_m <= 9e+59) {
tmp = (0.5 + (-0.25 * (re * re))) * ((1.0 / (1.0 + (im_m * (1.0 + (im_m * 0.5))))) - exp(im_m));
} else {
tmp = im_m * (cos(re) * t_0);
}
return im_s * tmp;
}
im\_m = abs(im)
im\_s = copysign(1.0d0, im)
real(8) function code(im_s, re, im_m)
real(8), intent (in) :: im_s
real(8), intent (in) :: re
real(8), intent (in) :: im_m
real(8) :: t_0
real(8) :: tmp
t_0 = (-1.0d0) + ((im_m * im_m) * ((-0.16666666666666666d0) + ((im_m * im_m) * (-0.008333333333333333d0))))
if (im_m <= 5.8d0) then
tmp = (im_m * cos(re)) * t_0
else if (im_m <= 9d+59) then
tmp = (0.5d0 + ((-0.25d0) * (re * re))) * ((1.0d0 / (1.0d0 + (im_m * (1.0d0 + (im_m * 0.5d0))))) - exp(im_m))
else
tmp = im_m * (cos(re) * t_0)
end if
code = im_s * tmp
end function
im\_m = Math.abs(im);
im\_s = Math.copySign(1.0, im);
public static double code(double im_s, double re, double im_m) {
double t_0 = -1.0 + ((im_m * im_m) * (-0.16666666666666666 + ((im_m * im_m) * -0.008333333333333333)));
double tmp;
if (im_m <= 5.8) {
tmp = (im_m * Math.cos(re)) * t_0;
} else if (im_m <= 9e+59) {
tmp = (0.5 + (-0.25 * (re * re))) * ((1.0 / (1.0 + (im_m * (1.0 + (im_m * 0.5))))) - Math.exp(im_m));
} else {
tmp = im_m * (Math.cos(re) * t_0);
}
return im_s * tmp;
}
im\_m = math.fabs(im) im\_s = math.copysign(1.0, im) def code(im_s, re, im_m): t_0 = -1.0 + ((im_m * im_m) * (-0.16666666666666666 + ((im_m * im_m) * -0.008333333333333333))) tmp = 0 if im_m <= 5.8: tmp = (im_m * math.cos(re)) * t_0 elif im_m <= 9e+59: tmp = (0.5 + (-0.25 * (re * re))) * ((1.0 / (1.0 + (im_m * (1.0 + (im_m * 0.5))))) - math.exp(im_m)) else: tmp = im_m * (math.cos(re) * t_0) return im_s * tmp
im\_m = abs(im) im\_s = copysign(1.0, im) function code(im_s, re, im_m) t_0 = Float64(-1.0 + Float64(Float64(im_m * im_m) * Float64(-0.16666666666666666 + Float64(Float64(im_m * im_m) * -0.008333333333333333)))) tmp = 0.0 if (im_m <= 5.8) tmp = Float64(Float64(im_m * cos(re)) * t_0); elseif (im_m <= 9e+59) tmp = Float64(Float64(0.5 + Float64(-0.25 * Float64(re * re))) * Float64(Float64(1.0 / Float64(1.0 + Float64(im_m * Float64(1.0 + Float64(im_m * 0.5))))) - exp(im_m))); else tmp = Float64(im_m * Float64(cos(re) * t_0)); end return Float64(im_s * tmp) end
im\_m = abs(im); im\_s = sign(im) * abs(1.0); function tmp_2 = code(im_s, re, im_m) t_0 = -1.0 + ((im_m * im_m) * (-0.16666666666666666 + ((im_m * im_m) * -0.008333333333333333))); tmp = 0.0; if (im_m <= 5.8) tmp = (im_m * cos(re)) * t_0; elseif (im_m <= 9e+59) tmp = (0.5 + (-0.25 * (re * re))) * ((1.0 / (1.0 + (im_m * (1.0 + (im_m * 0.5))))) - exp(im_m)); else tmp = im_m * (cos(re) * t_0); end tmp_2 = im_s * tmp; end
im\_m = N[Abs[im], $MachinePrecision]
im\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[im]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[im$95$s_, re_, im$95$m_] := Block[{t$95$0 = N[(-1.0 + N[(N[(im$95$m * im$95$m), $MachinePrecision] * N[(-0.16666666666666666 + N[(N[(im$95$m * im$95$m), $MachinePrecision] * -0.008333333333333333), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, N[(im$95$s * If[LessEqual[im$95$m, 5.8], N[(N[(im$95$m * N[Cos[re], $MachinePrecision]), $MachinePrecision] * t$95$0), $MachinePrecision], If[LessEqual[im$95$m, 9e+59], N[(N[(0.5 + N[(-0.25 * N[(re * re), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[(1.0 / N[(1.0 + N[(im$95$m * N[(1.0 + N[(im$95$m * 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[Exp[im$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(im$95$m * N[(N[Cos[re], $MachinePrecision] * t$95$0), $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]]
\begin{array}{l}
im\_m = \left|im\right|
\\
im\_s = \mathsf{copysign}\left(1, im\right)
\\
\begin{array}{l}
t_0 := -1 + \left(im\_m \cdot im\_m\right) \cdot \left(-0.16666666666666666 + \left(im\_m \cdot im\_m\right) \cdot -0.008333333333333333\right)\\
im\_s \cdot \begin{array}{l}
\mathbf{if}\;im\_m \leq 5.8:\\
\;\;\;\;\left(im\_m \cdot \cos re\right) \cdot t\_0\\
\mathbf{elif}\;im\_m \leq 9 \cdot 10^{+59}:\\
\;\;\;\;\left(0.5 + -0.25 \cdot \left(re \cdot re\right)\right) \cdot \left(\frac{1}{1 + im\_m \cdot \left(1 + im\_m \cdot 0.5\right)} - e^{im\_m}\right)\\
\mathbf{else}:\\
\;\;\;\;im\_m \cdot \left(\cos re \cdot t\_0\right)\\
\end{array}
\end{array}
\end{array}
if im < 5.79999999999999982Initial program 46.4%
Taylor expanded in im around 0
*-lowering-*.f64N/A
*-commutativeN/A
+-commutativeN/A
distribute-lft-inN/A
associate-*r*N/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
distribute-rgt-outN/A
distribute-lft-outN/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
Simplified90.0%
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
cos-lowering-cos.f6490.0%
Applied egg-rr90.0%
if 5.79999999999999982 < im < 8.99999999999999919e59Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
associate-*r*N/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
--lowering--.f64N/A
exp-negN/A
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
exp-lowering-exp.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6466.7%
Simplified66.7%
Taylor expanded in im around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6466.4%
Simplified66.4%
if 8.99999999999999919e59 < im Initial program 100.0%
Taylor expanded in im around 0
*-lowering-*.f64N/A
*-commutativeN/A
+-commutativeN/A
distribute-lft-inN/A
associate-*r*N/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
distribute-rgt-outN/A
distribute-lft-outN/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
Simplified100.0%
Final simplification90.0%
im\_m = (fabs.f64 im)
im\_s = (copysign.f64 #s(literal 1 binary64) im)
(FPCore (im_s re im_m)
:precision binary64
(let* ((t_0
(+
-1.0
(*
(* im_m im_m)
(+ -0.16666666666666666 (* (* im_m im_m) -0.008333333333333333))))))
(*
im_s
(if (<= im_m 5.5)
(* (* im_m (cos re)) t_0)
(if (<= im_m 9e+59)
(* (+ 0.5 (* -0.25 (* re re))) (- (/ 1.0 (+ im_m 1.0)) (exp im_m)))
(* im_m (* (cos re) t_0)))))))im\_m = fabs(im);
im\_s = copysign(1.0, im);
double code(double im_s, double re, double im_m) {
double t_0 = -1.0 + ((im_m * im_m) * (-0.16666666666666666 + ((im_m * im_m) * -0.008333333333333333)));
double tmp;
if (im_m <= 5.5) {
tmp = (im_m * cos(re)) * t_0;
} else if (im_m <= 9e+59) {
tmp = (0.5 + (-0.25 * (re * re))) * ((1.0 / (im_m + 1.0)) - exp(im_m));
} else {
tmp = im_m * (cos(re) * t_0);
}
return im_s * tmp;
}
im\_m = abs(im)
im\_s = copysign(1.0d0, im)
real(8) function code(im_s, re, im_m)
real(8), intent (in) :: im_s
real(8), intent (in) :: re
real(8), intent (in) :: im_m
real(8) :: t_0
real(8) :: tmp
t_0 = (-1.0d0) + ((im_m * im_m) * ((-0.16666666666666666d0) + ((im_m * im_m) * (-0.008333333333333333d0))))
if (im_m <= 5.5d0) then
tmp = (im_m * cos(re)) * t_0
else if (im_m <= 9d+59) then
tmp = (0.5d0 + ((-0.25d0) * (re * re))) * ((1.0d0 / (im_m + 1.0d0)) - exp(im_m))
else
tmp = im_m * (cos(re) * t_0)
end if
code = im_s * tmp
end function
im\_m = Math.abs(im);
im\_s = Math.copySign(1.0, im);
public static double code(double im_s, double re, double im_m) {
double t_0 = -1.0 + ((im_m * im_m) * (-0.16666666666666666 + ((im_m * im_m) * -0.008333333333333333)));
double tmp;
if (im_m <= 5.5) {
tmp = (im_m * Math.cos(re)) * t_0;
} else if (im_m <= 9e+59) {
tmp = (0.5 + (-0.25 * (re * re))) * ((1.0 / (im_m + 1.0)) - Math.exp(im_m));
} else {
tmp = im_m * (Math.cos(re) * t_0);
}
return im_s * tmp;
}
im\_m = math.fabs(im) im\_s = math.copysign(1.0, im) def code(im_s, re, im_m): t_0 = -1.0 + ((im_m * im_m) * (-0.16666666666666666 + ((im_m * im_m) * -0.008333333333333333))) tmp = 0 if im_m <= 5.5: tmp = (im_m * math.cos(re)) * t_0 elif im_m <= 9e+59: tmp = (0.5 + (-0.25 * (re * re))) * ((1.0 / (im_m + 1.0)) - math.exp(im_m)) else: tmp = im_m * (math.cos(re) * t_0) return im_s * tmp
im\_m = abs(im) im\_s = copysign(1.0, im) function code(im_s, re, im_m) t_0 = Float64(-1.0 + Float64(Float64(im_m * im_m) * Float64(-0.16666666666666666 + Float64(Float64(im_m * im_m) * -0.008333333333333333)))) tmp = 0.0 if (im_m <= 5.5) tmp = Float64(Float64(im_m * cos(re)) * t_0); elseif (im_m <= 9e+59) tmp = Float64(Float64(0.5 + Float64(-0.25 * Float64(re * re))) * Float64(Float64(1.0 / Float64(im_m + 1.0)) - exp(im_m))); else tmp = Float64(im_m * Float64(cos(re) * t_0)); end return Float64(im_s * tmp) end
im\_m = abs(im); im\_s = sign(im) * abs(1.0); function tmp_2 = code(im_s, re, im_m) t_0 = -1.0 + ((im_m * im_m) * (-0.16666666666666666 + ((im_m * im_m) * -0.008333333333333333))); tmp = 0.0; if (im_m <= 5.5) tmp = (im_m * cos(re)) * t_0; elseif (im_m <= 9e+59) tmp = (0.5 + (-0.25 * (re * re))) * ((1.0 / (im_m + 1.0)) - exp(im_m)); else tmp = im_m * (cos(re) * t_0); end tmp_2 = im_s * tmp; end
im\_m = N[Abs[im], $MachinePrecision]
im\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[im]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[im$95$s_, re_, im$95$m_] := Block[{t$95$0 = N[(-1.0 + N[(N[(im$95$m * im$95$m), $MachinePrecision] * N[(-0.16666666666666666 + N[(N[(im$95$m * im$95$m), $MachinePrecision] * -0.008333333333333333), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, N[(im$95$s * If[LessEqual[im$95$m, 5.5], N[(N[(im$95$m * N[Cos[re], $MachinePrecision]), $MachinePrecision] * t$95$0), $MachinePrecision], If[LessEqual[im$95$m, 9e+59], N[(N[(0.5 + N[(-0.25 * N[(re * re), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[(1.0 / N[(im$95$m + 1.0), $MachinePrecision]), $MachinePrecision] - N[Exp[im$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(im$95$m * N[(N[Cos[re], $MachinePrecision] * t$95$0), $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]]
\begin{array}{l}
im\_m = \left|im\right|
\\
im\_s = \mathsf{copysign}\left(1, im\right)
\\
\begin{array}{l}
t_0 := -1 + \left(im\_m \cdot im\_m\right) \cdot \left(-0.16666666666666666 + \left(im\_m \cdot im\_m\right) \cdot -0.008333333333333333\right)\\
im\_s \cdot \begin{array}{l}
\mathbf{if}\;im\_m \leq 5.5:\\
\;\;\;\;\left(im\_m \cdot \cos re\right) \cdot t\_0\\
\mathbf{elif}\;im\_m \leq 9 \cdot 10^{+59}:\\
\;\;\;\;\left(0.5 + -0.25 \cdot \left(re \cdot re\right)\right) \cdot \left(\frac{1}{im\_m + 1} - e^{im\_m}\right)\\
\mathbf{else}:\\
\;\;\;\;im\_m \cdot \left(\cos re \cdot t\_0\right)\\
\end{array}
\end{array}
\end{array}
if im < 5.5Initial program 46.4%
Taylor expanded in im around 0
*-lowering-*.f64N/A
*-commutativeN/A
+-commutativeN/A
distribute-lft-inN/A
associate-*r*N/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
distribute-rgt-outN/A
distribute-lft-outN/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
Simplified90.0%
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
cos-lowering-cos.f6490.0%
Applied egg-rr90.0%
if 5.5 < im < 8.99999999999999919e59Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
associate-*r*N/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
--lowering--.f64N/A
exp-negN/A
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
exp-lowering-exp.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6466.7%
Simplified66.7%
Taylor expanded in im around 0
+-commutativeN/A
+-lowering-+.f6466.1%
Simplified66.1%
if 8.99999999999999919e59 < im Initial program 100.0%
Taylor expanded in im around 0
*-lowering-*.f64N/A
*-commutativeN/A
+-commutativeN/A
distribute-lft-inN/A
associate-*r*N/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
distribute-rgt-outN/A
distribute-lft-outN/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
Simplified100.0%
Final simplification90.0%
im\_m = (fabs.f64 im)
im\_s = (copysign.f64 #s(literal 1 binary64) im)
(FPCore (im_s re im_m)
:precision binary64
(let* ((t_0
(*
im_m
(*
(cos re)
(+
-1.0
(*
(* im_m im_m)
(+
-0.16666666666666666
(* (* im_m im_m) -0.008333333333333333))))))))
(*
im_s
(if (<= im_m 5.2)
t_0
(if (<= im_m 9e+59)
(* (+ 0.5 (* -0.25 (* re re))) (- (/ 1.0 (+ im_m 1.0)) (exp im_m)))
t_0)))))im\_m = fabs(im);
im\_s = copysign(1.0, im);
double code(double im_s, double re, double im_m) {
double t_0 = im_m * (cos(re) * (-1.0 + ((im_m * im_m) * (-0.16666666666666666 + ((im_m * im_m) * -0.008333333333333333)))));
double tmp;
if (im_m <= 5.2) {
tmp = t_0;
} else if (im_m <= 9e+59) {
tmp = (0.5 + (-0.25 * (re * re))) * ((1.0 / (im_m + 1.0)) - exp(im_m));
} else {
tmp = t_0;
}
return im_s * tmp;
}
im\_m = abs(im)
im\_s = copysign(1.0d0, im)
real(8) function code(im_s, re, im_m)
real(8), intent (in) :: im_s
real(8), intent (in) :: re
real(8), intent (in) :: im_m
real(8) :: t_0
real(8) :: tmp
t_0 = im_m * (cos(re) * ((-1.0d0) + ((im_m * im_m) * ((-0.16666666666666666d0) + ((im_m * im_m) * (-0.008333333333333333d0))))))
if (im_m <= 5.2d0) then
tmp = t_0
else if (im_m <= 9d+59) then
tmp = (0.5d0 + ((-0.25d0) * (re * re))) * ((1.0d0 / (im_m + 1.0d0)) - exp(im_m))
else
tmp = t_0
end if
code = im_s * tmp
end function
im\_m = Math.abs(im);
im\_s = Math.copySign(1.0, im);
public static double code(double im_s, double re, double im_m) {
double t_0 = im_m * (Math.cos(re) * (-1.0 + ((im_m * im_m) * (-0.16666666666666666 + ((im_m * im_m) * -0.008333333333333333)))));
double tmp;
if (im_m <= 5.2) {
tmp = t_0;
} else if (im_m <= 9e+59) {
tmp = (0.5 + (-0.25 * (re * re))) * ((1.0 / (im_m + 1.0)) - Math.exp(im_m));
} else {
tmp = t_0;
}
return im_s * tmp;
}
im\_m = math.fabs(im) im\_s = math.copysign(1.0, im) def code(im_s, re, im_m): t_0 = im_m * (math.cos(re) * (-1.0 + ((im_m * im_m) * (-0.16666666666666666 + ((im_m * im_m) * -0.008333333333333333))))) tmp = 0 if im_m <= 5.2: tmp = t_0 elif im_m <= 9e+59: tmp = (0.5 + (-0.25 * (re * re))) * ((1.0 / (im_m + 1.0)) - math.exp(im_m)) else: tmp = t_0 return im_s * tmp
im\_m = abs(im) im\_s = copysign(1.0, im) function code(im_s, re, im_m) t_0 = Float64(im_m * Float64(cos(re) * Float64(-1.0 + Float64(Float64(im_m * im_m) * Float64(-0.16666666666666666 + Float64(Float64(im_m * im_m) * -0.008333333333333333)))))) tmp = 0.0 if (im_m <= 5.2) tmp = t_0; elseif (im_m <= 9e+59) tmp = Float64(Float64(0.5 + Float64(-0.25 * Float64(re * re))) * Float64(Float64(1.0 / Float64(im_m + 1.0)) - exp(im_m))); else tmp = t_0; end return Float64(im_s * tmp) end
im\_m = abs(im); im\_s = sign(im) * abs(1.0); function tmp_2 = code(im_s, re, im_m) t_0 = im_m * (cos(re) * (-1.0 + ((im_m * im_m) * (-0.16666666666666666 + ((im_m * im_m) * -0.008333333333333333))))); tmp = 0.0; if (im_m <= 5.2) tmp = t_0; elseif (im_m <= 9e+59) tmp = (0.5 + (-0.25 * (re * re))) * ((1.0 / (im_m + 1.0)) - exp(im_m)); else tmp = t_0; end tmp_2 = im_s * tmp; end
im\_m = N[Abs[im], $MachinePrecision]
im\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[im]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[im$95$s_, re_, im$95$m_] := Block[{t$95$0 = N[(im$95$m * N[(N[Cos[re], $MachinePrecision] * N[(-1.0 + N[(N[(im$95$m * im$95$m), $MachinePrecision] * N[(-0.16666666666666666 + N[(N[(im$95$m * im$95$m), $MachinePrecision] * -0.008333333333333333), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, N[(im$95$s * If[LessEqual[im$95$m, 5.2], t$95$0, If[LessEqual[im$95$m, 9e+59], N[(N[(0.5 + N[(-0.25 * N[(re * re), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[(1.0 / N[(im$95$m + 1.0), $MachinePrecision]), $MachinePrecision] - N[Exp[im$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$0]]), $MachinePrecision]]
\begin{array}{l}
im\_m = \left|im\right|
\\
im\_s = \mathsf{copysign}\left(1, im\right)
\\
\begin{array}{l}
t_0 := im\_m \cdot \left(\cos re \cdot \left(-1 + \left(im\_m \cdot im\_m\right) \cdot \left(-0.16666666666666666 + \left(im\_m \cdot im\_m\right) \cdot -0.008333333333333333\right)\right)\right)\\
im\_s \cdot \begin{array}{l}
\mathbf{if}\;im\_m \leq 5.2:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;im\_m \leq 9 \cdot 10^{+59}:\\
\;\;\;\;\left(0.5 + -0.25 \cdot \left(re \cdot re\right)\right) \cdot \left(\frac{1}{im\_m + 1} - e^{im\_m}\right)\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
\end{array}
if im < 5.20000000000000018 or 8.99999999999999919e59 < im Initial program 56.3%
Taylor expanded in im around 0
*-lowering-*.f64N/A
*-commutativeN/A
+-commutativeN/A
distribute-lft-inN/A
associate-*r*N/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
distribute-rgt-outN/A
distribute-lft-outN/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
Simplified91.8%
if 5.20000000000000018 < im < 8.99999999999999919e59Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
associate-*r*N/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
--lowering--.f64N/A
exp-negN/A
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
exp-lowering-exp.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6466.7%
Simplified66.7%
Taylor expanded in im around 0
+-commutativeN/A
+-lowering-+.f6466.1%
Simplified66.1%
Final simplification90.0%
im\_m = (fabs.f64 im)
im\_s = (copysign.f64 #s(literal 1 binary64) im)
(FPCore (im_s re im_m)
:precision binary64
(let* ((t_0 (+ -1.0 (* (* im_m im_m) -0.16666666666666666))))
(*
im_s
(if (<= im_m 3.45)
(* (cos re) (* im_m t_0))
(if (<= im_m 1.15e+103)
(* (+ 0.5 (* -0.25 (* re re))) (- (/ 1.0 (+ im_m 1.0)) (exp im_m)))
(* im_m (* (cos re) t_0)))))))im\_m = fabs(im);
im\_s = copysign(1.0, im);
double code(double im_s, double re, double im_m) {
double t_0 = -1.0 + ((im_m * im_m) * -0.16666666666666666);
double tmp;
if (im_m <= 3.45) {
tmp = cos(re) * (im_m * t_0);
} else if (im_m <= 1.15e+103) {
tmp = (0.5 + (-0.25 * (re * re))) * ((1.0 / (im_m + 1.0)) - exp(im_m));
} else {
tmp = im_m * (cos(re) * t_0);
}
return im_s * tmp;
}
im\_m = abs(im)
im\_s = copysign(1.0d0, im)
real(8) function code(im_s, re, im_m)
real(8), intent (in) :: im_s
real(8), intent (in) :: re
real(8), intent (in) :: im_m
real(8) :: t_0
real(8) :: tmp
t_0 = (-1.0d0) + ((im_m * im_m) * (-0.16666666666666666d0))
if (im_m <= 3.45d0) then
tmp = cos(re) * (im_m * t_0)
else if (im_m <= 1.15d+103) then
tmp = (0.5d0 + ((-0.25d0) * (re * re))) * ((1.0d0 / (im_m + 1.0d0)) - exp(im_m))
else
tmp = im_m * (cos(re) * t_0)
end if
code = im_s * tmp
end function
im\_m = Math.abs(im);
im\_s = Math.copySign(1.0, im);
public static double code(double im_s, double re, double im_m) {
double t_0 = -1.0 + ((im_m * im_m) * -0.16666666666666666);
double tmp;
if (im_m <= 3.45) {
tmp = Math.cos(re) * (im_m * t_0);
} else if (im_m <= 1.15e+103) {
tmp = (0.5 + (-0.25 * (re * re))) * ((1.0 / (im_m + 1.0)) - Math.exp(im_m));
} else {
tmp = im_m * (Math.cos(re) * t_0);
}
return im_s * tmp;
}
im\_m = math.fabs(im) im\_s = math.copysign(1.0, im) def code(im_s, re, im_m): t_0 = -1.0 + ((im_m * im_m) * -0.16666666666666666) tmp = 0 if im_m <= 3.45: tmp = math.cos(re) * (im_m * t_0) elif im_m <= 1.15e+103: tmp = (0.5 + (-0.25 * (re * re))) * ((1.0 / (im_m + 1.0)) - math.exp(im_m)) else: tmp = im_m * (math.cos(re) * t_0) return im_s * tmp
im\_m = abs(im) im\_s = copysign(1.0, im) function code(im_s, re, im_m) t_0 = Float64(-1.0 + Float64(Float64(im_m * im_m) * -0.16666666666666666)) tmp = 0.0 if (im_m <= 3.45) tmp = Float64(cos(re) * Float64(im_m * t_0)); elseif (im_m <= 1.15e+103) tmp = Float64(Float64(0.5 + Float64(-0.25 * Float64(re * re))) * Float64(Float64(1.0 / Float64(im_m + 1.0)) - exp(im_m))); else tmp = Float64(im_m * Float64(cos(re) * t_0)); end return Float64(im_s * tmp) end
im\_m = abs(im); im\_s = sign(im) * abs(1.0); function tmp_2 = code(im_s, re, im_m) t_0 = -1.0 + ((im_m * im_m) * -0.16666666666666666); tmp = 0.0; if (im_m <= 3.45) tmp = cos(re) * (im_m * t_0); elseif (im_m <= 1.15e+103) tmp = (0.5 + (-0.25 * (re * re))) * ((1.0 / (im_m + 1.0)) - exp(im_m)); else tmp = im_m * (cos(re) * t_0); end tmp_2 = im_s * tmp; end
im\_m = N[Abs[im], $MachinePrecision]
im\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[im]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[im$95$s_, re_, im$95$m_] := Block[{t$95$0 = N[(-1.0 + N[(N[(im$95$m * im$95$m), $MachinePrecision] * -0.16666666666666666), $MachinePrecision]), $MachinePrecision]}, N[(im$95$s * If[LessEqual[im$95$m, 3.45], N[(N[Cos[re], $MachinePrecision] * N[(im$95$m * t$95$0), $MachinePrecision]), $MachinePrecision], If[LessEqual[im$95$m, 1.15e+103], N[(N[(0.5 + N[(-0.25 * N[(re * re), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[(1.0 / N[(im$95$m + 1.0), $MachinePrecision]), $MachinePrecision] - N[Exp[im$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(im$95$m * N[(N[Cos[re], $MachinePrecision] * t$95$0), $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]]
\begin{array}{l}
im\_m = \left|im\right|
\\
im\_s = \mathsf{copysign}\left(1, im\right)
\\
\begin{array}{l}
t_0 := -1 + \left(im\_m \cdot im\_m\right) \cdot -0.16666666666666666\\
im\_s \cdot \begin{array}{l}
\mathbf{if}\;im\_m \leq 3.45:\\
\;\;\;\;\cos re \cdot \left(im\_m \cdot t\_0\right)\\
\mathbf{elif}\;im\_m \leq 1.15 \cdot 10^{+103}:\\
\;\;\;\;\left(0.5 + -0.25 \cdot \left(re \cdot re\right)\right) \cdot \left(\frac{1}{im\_m + 1} - e^{im\_m}\right)\\
\mathbf{else}:\\
\;\;\;\;im\_m \cdot \left(\cos re \cdot t\_0\right)\\
\end{array}
\end{array}
\end{array}
if im < 3.4500000000000002Initial program 46.4%
Taylor expanded in im around 0
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
*-commutativeN/A
associate-*r*N/A
distribute-lft-out--N/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
associate-*r*N/A
distribute-rgt-out--N/A
unsub-negN/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
neg-mul-1N/A
*-commutativeN/A
Simplified86.4%
if 3.4500000000000002 < im < 1.15000000000000004e103Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
associate-*r*N/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
--lowering--.f64N/A
exp-negN/A
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
exp-lowering-exp.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6476.0%
Simplified76.0%
Taylor expanded in im around 0
+-commutativeN/A
+-lowering-+.f6475.6%
Simplified75.6%
if 1.15000000000000004e103 < im Initial program 100.0%
Taylor expanded in im around 0
*-lowering-*.f64N/A
distribute-rgt-inN/A
associate-+r+N/A
*-commutativeN/A
+-commutativeN/A
+-commutativeN/A
Simplified100.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
metadata-evalN/A
metadata-evalN/A
associate-*l*N/A
distribute-rgt-out--N/A
*-commutativeN/A
+-commutativeN/A
*-lowering-*.f64N/A
Simplified100.0%
Final simplification87.3%
im\_m = (fabs.f64 im)
im\_s = (copysign.f64 #s(literal 1 binary64) im)
(FPCore (im_s re im_m)
:precision binary64
(let* ((t_0 (+ -1.0 (* (* im_m im_m) -0.16666666666666666))))
(*
im_s
(if (<= im_m 4.7)
(* (cos re) (* im_m t_0))
(if (<= im_m 1.1e+103)
(* (+ 0.5 (* -0.25 (* re re))) (- 1.0 (exp im_m)))
(* im_m (* (cos re) t_0)))))))im\_m = fabs(im);
im\_s = copysign(1.0, im);
double code(double im_s, double re, double im_m) {
double t_0 = -1.0 + ((im_m * im_m) * -0.16666666666666666);
double tmp;
if (im_m <= 4.7) {
tmp = cos(re) * (im_m * t_0);
} else if (im_m <= 1.1e+103) {
tmp = (0.5 + (-0.25 * (re * re))) * (1.0 - exp(im_m));
} else {
tmp = im_m * (cos(re) * t_0);
}
return im_s * tmp;
}
im\_m = abs(im)
im\_s = copysign(1.0d0, im)
real(8) function code(im_s, re, im_m)
real(8), intent (in) :: im_s
real(8), intent (in) :: re
real(8), intent (in) :: im_m
real(8) :: t_0
real(8) :: tmp
t_0 = (-1.0d0) + ((im_m * im_m) * (-0.16666666666666666d0))
if (im_m <= 4.7d0) then
tmp = cos(re) * (im_m * t_0)
else if (im_m <= 1.1d+103) then
tmp = (0.5d0 + ((-0.25d0) * (re * re))) * (1.0d0 - exp(im_m))
else
tmp = im_m * (cos(re) * t_0)
end if
code = im_s * tmp
end function
im\_m = Math.abs(im);
im\_s = Math.copySign(1.0, im);
public static double code(double im_s, double re, double im_m) {
double t_0 = -1.0 + ((im_m * im_m) * -0.16666666666666666);
double tmp;
if (im_m <= 4.7) {
tmp = Math.cos(re) * (im_m * t_0);
} else if (im_m <= 1.1e+103) {
tmp = (0.5 + (-0.25 * (re * re))) * (1.0 - Math.exp(im_m));
} else {
tmp = im_m * (Math.cos(re) * t_0);
}
return im_s * tmp;
}
im\_m = math.fabs(im) im\_s = math.copysign(1.0, im) def code(im_s, re, im_m): t_0 = -1.0 + ((im_m * im_m) * -0.16666666666666666) tmp = 0 if im_m <= 4.7: tmp = math.cos(re) * (im_m * t_0) elif im_m <= 1.1e+103: tmp = (0.5 + (-0.25 * (re * re))) * (1.0 - math.exp(im_m)) else: tmp = im_m * (math.cos(re) * t_0) return im_s * tmp
im\_m = abs(im) im\_s = copysign(1.0, im) function code(im_s, re, im_m) t_0 = Float64(-1.0 + Float64(Float64(im_m * im_m) * -0.16666666666666666)) tmp = 0.0 if (im_m <= 4.7) tmp = Float64(cos(re) * Float64(im_m * t_0)); elseif (im_m <= 1.1e+103) tmp = Float64(Float64(0.5 + Float64(-0.25 * Float64(re * re))) * Float64(1.0 - exp(im_m))); else tmp = Float64(im_m * Float64(cos(re) * t_0)); end return Float64(im_s * tmp) end
im\_m = abs(im); im\_s = sign(im) * abs(1.0); function tmp_2 = code(im_s, re, im_m) t_0 = -1.0 + ((im_m * im_m) * -0.16666666666666666); tmp = 0.0; if (im_m <= 4.7) tmp = cos(re) * (im_m * t_0); elseif (im_m <= 1.1e+103) tmp = (0.5 + (-0.25 * (re * re))) * (1.0 - exp(im_m)); else tmp = im_m * (cos(re) * t_0); end tmp_2 = im_s * tmp; end
im\_m = N[Abs[im], $MachinePrecision]
im\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[im]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[im$95$s_, re_, im$95$m_] := Block[{t$95$0 = N[(-1.0 + N[(N[(im$95$m * im$95$m), $MachinePrecision] * -0.16666666666666666), $MachinePrecision]), $MachinePrecision]}, N[(im$95$s * If[LessEqual[im$95$m, 4.7], N[(N[Cos[re], $MachinePrecision] * N[(im$95$m * t$95$0), $MachinePrecision]), $MachinePrecision], If[LessEqual[im$95$m, 1.1e+103], N[(N[(0.5 + N[(-0.25 * N[(re * re), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(1.0 - N[Exp[im$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(im$95$m * N[(N[Cos[re], $MachinePrecision] * t$95$0), $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]]
\begin{array}{l}
im\_m = \left|im\right|
\\
im\_s = \mathsf{copysign}\left(1, im\right)
\\
\begin{array}{l}
t_0 := -1 + \left(im\_m \cdot im\_m\right) \cdot -0.16666666666666666\\
im\_s \cdot \begin{array}{l}
\mathbf{if}\;im\_m \leq 4.7:\\
\;\;\;\;\cos re \cdot \left(im\_m \cdot t\_0\right)\\
\mathbf{elif}\;im\_m \leq 1.1 \cdot 10^{+103}:\\
\;\;\;\;\left(0.5 + -0.25 \cdot \left(re \cdot re\right)\right) \cdot \left(1 - e^{im\_m}\right)\\
\mathbf{else}:\\
\;\;\;\;im\_m \cdot \left(\cos re \cdot t\_0\right)\\
\end{array}
\end{array}
\end{array}
if im < 4.70000000000000018Initial program 46.4%
Taylor expanded in im around 0
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
*-commutativeN/A
associate-*r*N/A
distribute-lft-out--N/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
associate-*r*N/A
distribute-rgt-out--N/A
unsub-negN/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
neg-mul-1N/A
*-commutativeN/A
Simplified86.4%
if 4.70000000000000018 < im < 1.09999999999999996e103Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
associate-*r*N/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
--lowering--.f64N/A
exp-negN/A
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
exp-lowering-exp.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6476.0%
Simplified76.0%
Taylor expanded in im around 0
Simplified75.3%
if 1.09999999999999996e103 < im Initial program 100.0%
Taylor expanded in im around 0
*-lowering-*.f64N/A
distribute-rgt-inN/A
associate-+r+N/A
*-commutativeN/A
+-commutativeN/A
+-commutativeN/A
Simplified100.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
metadata-evalN/A
metadata-evalN/A
associate-*l*N/A
distribute-rgt-out--N/A
*-commutativeN/A
+-commutativeN/A
*-lowering-*.f64N/A
Simplified100.0%
Final simplification87.3%
im\_m = (fabs.f64 im)
im\_s = (copysign.f64 #s(literal 1 binary64) im)
(FPCore (im_s re im_m)
:precision binary64
(let* ((t_0 (+ -1.0 (* (* im_m im_m) -0.16666666666666666))))
(*
im_s
(if (<= im_m 580.0)
(* (cos re) (* im_m t_0))
(if (<= im_m 1.25e+103)
(*
(+
-1.0
(*
(* im_m im_m)
(+
-0.16666666666666666
(*
im_m
(*
im_m
(+
-0.008333333333333333
(* (* im_m im_m) -0.0001984126984126984)))))))
(+
im_m
(*
(* re re)
(+ (* im_m -0.5) (* 0.041666666666666664 (* im_m (* re re)))))))
(* im_m (* (cos re) t_0)))))))im\_m = fabs(im);
im\_s = copysign(1.0, im);
double code(double im_s, double re, double im_m) {
double t_0 = -1.0 + ((im_m * im_m) * -0.16666666666666666);
double tmp;
if (im_m <= 580.0) {
tmp = cos(re) * (im_m * t_0);
} else if (im_m <= 1.25e+103) {
tmp = (-1.0 + ((im_m * im_m) * (-0.16666666666666666 + (im_m * (im_m * (-0.008333333333333333 + ((im_m * im_m) * -0.0001984126984126984))))))) * (im_m + ((re * re) * ((im_m * -0.5) + (0.041666666666666664 * (im_m * (re * re))))));
} else {
tmp = im_m * (cos(re) * t_0);
}
return im_s * tmp;
}
im\_m = abs(im)
im\_s = copysign(1.0d0, im)
real(8) function code(im_s, re, im_m)
real(8), intent (in) :: im_s
real(8), intent (in) :: re
real(8), intent (in) :: im_m
real(8) :: t_0
real(8) :: tmp
t_0 = (-1.0d0) + ((im_m * im_m) * (-0.16666666666666666d0))
if (im_m <= 580.0d0) then
tmp = cos(re) * (im_m * t_0)
else if (im_m <= 1.25d+103) then
tmp = ((-1.0d0) + ((im_m * im_m) * ((-0.16666666666666666d0) + (im_m * (im_m * ((-0.008333333333333333d0) + ((im_m * im_m) * (-0.0001984126984126984d0)))))))) * (im_m + ((re * re) * ((im_m * (-0.5d0)) + (0.041666666666666664d0 * (im_m * (re * re))))))
else
tmp = im_m * (cos(re) * t_0)
end if
code = im_s * tmp
end function
im\_m = Math.abs(im);
im\_s = Math.copySign(1.0, im);
public static double code(double im_s, double re, double im_m) {
double t_0 = -1.0 + ((im_m * im_m) * -0.16666666666666666);
double tmp;
if (im_m <= 580.0) {
tmp = Math.cos(re) * (im_m * t_0);
} else if (im_m <= 1.25e+103) {
tmp = (-1.0 + ((im_m * im_m) * (-0.16666666666666666 + (im_m * (im_m * (-0.008333333333333333 + ((im_m * im_m) * -0.0001984126984126984))))))) * (im_m + ((re * re) * ((im_m * -0.5) + (0.041666666666666664 * (im_m * (re * re))))));
} else {
tmp = im_m * (Math.cos(re) * t_0);
}
return im_s * tmp;
}
im\_m = math.fabs(im) im\_s = math.copysign(1.0, im) def code(im_s, re, im_m): t_0 = -1.0 + ((im_m * im_m) * -0.16666666666666666) tmp = 0 if im_m <= 580.0: tmp = math.cos(re) * (im_m * t_0) elif im_m <= 1.25e+103: tmp = (-1.0 + ((im_m * im_m) * (-0.16666666666666666 + (im_m * (im_m * (-0.008333333333333333 + ((im_m * im_m) * -0.0001984126984126984))))))) * (im_m + ((re * re) * ((im_m * -0.5) + (0.041666666666666664 * (im_m * (re * re)))))) else: tmp = im_m * (math.cos(re) * t_0) return im_s * tmp
im\_m = abs(im) im\_s = copysign(1.0, im) function code(im_s, re, im_m) t_0 = Float64(-1.0 + Float64(Float64(im_m * im_m) * -0.16666666666666666)) tmp = 0.0 if (im_m <= 580.0) tmp = Float64(cos(re) * Float64(im_m * t_0)); elseif (im_m <= 1.25e+103) tmp = Float64(Float64(-1.0 + Float64(Float64(im_m * im_m) * Float64(-0.16666666666666666 + Float64(im_m * Float64(im_m * Float64(-0.008333333333333333 + Float64(Float64(im_m * im_m) * -0.0001984126984126984))))))) * Float64(im_m + Float64(Float64(re * re) * Float64(Float64(im_m * -0.5) + Float64(0.041666666666666664 * Float64(im_m * Float64(re * re))))))); else tmp = Float64(im_m * Float64(cos(re) * t_0)); end return Float64(im_s * tmp) end
im\_m = abs(im); im\_s = sign(im) * abs(1.0); function tmp_2 = code(im_s, re, im_m) t_0 = -1.0 + ((im_m * im_m) * -0.16666666666666666); tmp = 0.0; if (im_m <= 580.0) tmp = cos(re) * (im_m * t_0); elseif (im_m <= 1.25e+103) tmp = (-1.0 + ((im_m * im_m) * (-0.16666666666666666 + (im_m * (im_m * (-0.008333333333333333 + ((im_m * im_m) * -0.0001984126984126984))))))) * (im_m + ((re * re) * ((im_m * -0.5) + (0.041666666666666664 * (im_m * (re * re)))))); else tmp = im_m * (cos(re) * t_0); end tmp_2 = im_s * tmp; end
im\_m = N[Abs[im], $MachinePrecision]
im\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[im]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[im$95$s_, re_, im$95$m_] := Block[{t$95$0 = N[(-1.0 + N[(N[(im$95$m * im$95$m), $MachinePrecision] * -0.16666666666666666), $MachinePrecision]), $MachinePrecision]}, N[(im$95$s * If[LessEqual[im$95$m, 580.0], N[(N[Cos[re], $MachinePrecision] * N[(im$95$m * t$95$0), $MachinePrecision]), $MachinePrecision], If[LessEqual[im$95$m, 1.25e+103], N[(N[(-1.0 + N[(N[(im$95$m * im$95$m), $MachinePrecision] * N[(-0.16666666666666666 + N[(im$95$m * N[(im$95$m * N[(-0.008333333333333333 + N[(N[(im$95$m * im$95$m), $MachinePrecision] * -0.0001984126984126984), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(im$95$m + N[(N[(re * re), $MachinePrecision] * N[(N[(im$95$m * -0.5), $MachinePrecision] + N[(0.041666666666666664 * N[(im$95$m * N[(re * re), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(im$95$m * N[(N[Cos[re], $MachinePrecision] * t$95$0), $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]]
\begin{array}{l}
im\_m = \left|im\right|
\\
im\_s = \mathsf{copysign}\left(1, im\right)
\\
\begin{array}{l}
t_0 := -1 + \left(im\_m \cdot im\_m\right) \cdot -0.16666666666666666\\
im\_s \cdot \begin{array}{l}
\mathbf{if}\;im\_m \leq 580:\\
\;\;\;\;\cos re \cdot \left(im\_m \cdot t\_0\right)\\
\mathbf{elif}\;im\_m \leq 1.25 \cdot 10^{+103}:\\
\;\;\;\;\left(-1 + \left(im\_m \cdot im\_m\right) \cdot \left(-0.16666666666666666 + im\_m \cdot \left(im\_m \cdot \left(-0.008333333333333333 + \left(im\_m \cdot im\_m\right) \cdot -0.0001984126984126984\right)\right)\right)\right) \cdot \left(im\_m + \left(re \cdot re\right) \cdot \left(im\_m \cdot -0.5 + 0.041666666666666664 \cdot \left(im\_m \cdot \left(re \cdot re\right)\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;im\_m \cdot \left(\cos re \cdot t\_0\right)\\
\end{array}
\end{array}
\end{array}
if im < 580Initial program 46.7%
Taylor expanded in im around 0
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
*-commutativeN/A
associate-*r*N/A
distribute-lft-out--N/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
associate-*r*N/A
distribute-rgt-out--N/A
unsub-negN/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
neg-mul-1N/A
*-commutativeN/A
Simplified86.0%
if 580 < im < 1.25e103Initial program 100.0%
Taylor expanded in im around 0
*-lowering-*.f64N/A
distribute-rgt-inN/A
associate-+r+N/A
*-commutativeN/A
+-commutativeN/A
+-commutativeN/A
Simplified52.7%
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
Applied egg-rr52.7%
Taylor expanded in im around 0
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/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
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6452.7%
Simplified52.7%
Taylor expanded in re around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6459.3%
Simplified59.3%
if 1.25e103 < im Initial program 100.0%
Taylor expanded in im around 0
*-lowering-*.f64N/A
distribute-rgt-inN/A
associate-+r+N/A
*-commutativeN/A
+-commutativeN/A
+-commutativeN/A
Simplified100.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
metadata-evalN/A
metadata-evalN/A
associate-*l*N/A
distribute-rgt-out--N/A
*-commutativeN/A
+-commutativeN/A
*-lowering-*.f64N/A
Simplified100.0%
Final simplification85.5%
im\_m = (fabs.f64 im)
im\_s = (copysign.f64 #s(literal 1 binary64) im)
(FPCore (im_s re im_m)
:precision binary64
(let* ((t_0
(*
im_m
(* (cos re) (+ -1.0 (* (* im_m im_m) -0.16666666666666666))))))
(*
im_s
(if (<= im_m 460.0)
t_0
(if (<= im_m 1.4e+103)
(*
(+
-1.0
(*
(* im_m im_m)
(+
-0.16666666666666666
(*
im_m
(*
im_m
(+
-0.008333333333333333
(* (* im_m im_m) -0.0001984126984126984)))))))
(+
im_m
(*
(* re re)
(+ (* im_m -0.5) (* 0.041666666666666664 (* im_m (* re re)))))))
t_0)))))im\_m = fabs(im);
im\_s = copysign(1.0, im);
double code(double im_s, double re, double im_m) {
double t_0 = im_m * (cos(re) * (-1.0 + ((im_m * im_m) * -0.16666666666666666)));
double tmp;
if (im_m <= 460.0) {
tmp = t_0;
} else if (im_m <= 1.4e+103) {
tmp = (-1.0 + ((im_m * im_m) * (-0.16666666666666666 + (im_m * (im_m * (-0.008333333333333333 + ((im_m * im_m) * -0.0001984126984126984))))))) * (im_m + ((re * re) * ((im_m * -0.5) + (0.041666666666666664 * (im_m * (re * re))))));
} else {
tmp = t_0;
}
return im_s * tmp;
}
im\_m = abs(im)
im\_s = copysign(1.0d0, im)
real(8) function code(im_s, re, im_m)
real(8), intent (in) :: im_s
real(8), intent (in) :: re
real(8), intent (in) :: im_m
real(8) :: t_0
real(8) :: tmp
t_0 = im_m * (cos(re) * ((-1.0d0) + ((im_m * im_m) * (-0.16666666666666666d0))))
if (im_m <= 460.0d0) then
tmp = t_0
else if (im_m <= 1.4d+103) then
tmp = ((-1.0d0) + ((im_m * im_m) * ((-0.16666666666666666d0) + (im_m * (im_m * ((-0.008333333333333333d0) + ((im_m * im_m) * (-0.0001984126984126984d0)))))))) * (im_m + ((re * re) * ((im_m * (-0.5d0)) + (0.041666666666666664d0 * (im_m * (re * re))))))
else
tmp = t_0
end if
code = im_s * tmp
end function
im\_m = Math.abs(im);
im\_s = Math.copySign(1.0, im);
public static double code(double im_s, double re, double im_m) {
double t_0 = im_m * (Math.cos(re) * (-1.0 + ((im_m * im_m) * -0.16666666666666666)));
double tmp;
if (im_m <= 460.0) {
tmp = t_0;
} else if (im_m <= 1.4e+103) {
tmp = (-1.0 + ((im_m * im_m) * (-0.16666666666666666 + (im_m * (im_m * (-0.008333333333333333 + ((im_m * im_m) * -0.0001984126984126984))))))) * (im_m + ((re * re) * ((im_m * -0.5) + (0.041666666666666664 * (im_m * (re * re))))));
} else {
tmp = t_0;
}
return im_s * tmp;
}
im\_m = math.fabs(im) im\_s = math.copysign(1.0, im) def code(im_s, re, im_m): t_0 = im_m * (math.cos(re) * (-1.0 + ((im_m * im_m) * -0.16666666666666666))) tmp = 0 if im_m <= 460.0: tmp = t_0 elif im_m <= 1.4e+103: tmp = (-1.0 + ((im_m * im_m) * (-0.16666666666666666 + (im_m * (im_m * (-0.008333333333333333 + ((im_m * im_m) * -0.0001984126984126984))))))) * (im_m + ((re * re) * ((im_m * -0.5) + (0.041666666666666664 * (im_m * (re * re)))))) else: tmp = t_0 return im_s * tmp
im\_m = abs(im) im\_s = copysign(1.0, im) function code(im_s, re, im_m) t_0 = Float64(im_m * Float64(cos(re) * Float64(-1.0 + Float64(Float64(im_m * im_m) * -0.16666666666666666)))) tmp = 0.0 if (im_m <= 460.0) tmp = t_0; elseif (im_m <= 1.4e+103) tmp = Float64(Float64(-1.0 + Float64(Float64(im_m * im_m) * Float64(-0.16666666666666666 + Float64(im_m * Float64(im_m * Float64(-0.008333333333333333 + Float64(Float64(im_m * im_m) * -0.0001984126984126984))))))) * Float64(im_m + Float64(Float64(re * re) * Float64(Float64(im_m * -0.5) + Float64(0.041666666666666664 * Float64(im_m * Float64(re * re))))))); else tmp = t_0; end return Float64(im_s * tmp) end
im\_m = abs(im); im\_s = sign(im) * abs(1.0); function tmp_2 = code(im_s, re, im_m) t_0 = im_m * (cos(re) * (-1.0 + ((im_m * im_m) * -0.16666666666666666))); tmp = 0.0; if (im_m <= 460.0) tmp = t_0; elseif (im_m <= 1.4e+103) tmp = (-1.0 + ((im_m * im_m) * (-0.16666666666666666 + (im_m * (im_m * (-0.008333333333333333 + ((im_m * im_m) * -0.0001984126984126984))))))) * (im_m + ((re * re) * ((im_m * -0.5) + (0.041666666666666664 * (im_m * (re * re)))))); else tmp = t_0; end tmp_2 = im_s * tmp; end
im\_m = N[Abs[im], $MachinePrecision]
im\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[im]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[im$95$s_, re_, im$95$m_] := Block[{t$95$0 = N[(im$95$m * N[(N[Cos[re], $MachinePrecision] * N[(-1.0 + N[(N[(im$95$m * im$95$m), $MachinePrecision] * -0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, N[(im$95$s * If[LessEqual[im$95$m, 460.0], t$95$0, If[LessEqual[im$95$m, 1.4e+103], N[(N[(-1.0 + N[(N[(im$95$m * im$95$m), $MachinePrecision] * N[(-0.16666666666666666 + N[(im$95$m * N[(im$95$m * N[(-0.008333333333333333 + N[(N[(im$95$m * im$95$m), $MachinePrecision] * -0.0001984126984126984), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(im$95$m + N[(N[(re * re), $MachinePrecision] * N[(N[(im$95$m * -0.5), $MachinePrecision] + N[(0.041666666666666664 * N[(im$95$m * N[(re * re), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$0]]), $MachinePrecision]]
\begin{array}{l}
im\_m = \left|im\right|
\\
im\_s = \mathsf{copysign}\left(1, im\right)
\\
\begin{array}{l}
t_0 := im\_m \cdot \left(\cos re \cdot \left(-1 + \left(im\_m \cdot im\_m\right) \cdot -0.16666666666666666\right)\right)\\
im\_s \cdot \begin{array}{l}
\mathbf{if}\;im\_m \leq 460:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;im\_m \leq 1.4 \cdot 10^{+103}:\\
\;\;\;\;\left(-1 + \left(im\_m \cdot im\_m\right) \cdot \left(-0.16666666666666666 + im\_m \cdot \left(im\_m \cdot \left(-0.008333333333333333 + \left(im\_m \cdot im\_m\right) \cdot -0.0001984126984126984\right)\right)\right)\right) \cdot \left(im\_m + \left(re \cdot re\right) \cdot \left(im\_m \cdot -0.5 + 0.041666666666666664 \cdot \left(im\_m \cdot \left(re \cdot re\right)\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
\end{array}
if im < 460 or 1.40000000000000004e103 < im Initial program 55.2%
Taylor expanded in im around 0
*-lowering-*.f64N/A
distribute-rgt-inN/A
associate-+r+N/A
*-commutativeN/A
+-commutativeN/A
+-commutativeN/A
Simplified93.6%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
metadata-evalN/A
metadata-evalN/A
associate-*l*N/A
distribute-rgt-out--N/A
*-commutativeN/A
+-commutativeN/A
*-lowering-*.f64N/A
Simplified88.3%
if 460 < im < 1.40000000000000004e103Initial program 100.0%
Taylor expanded in im around 0
*-lowering-*.f64N/A
distribute-rgt-inN/A
associate-+r+N/A
*-commutativeN/A
+-commutativeN/A
+-commutativeN/A
Simplified52.7%
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
Applied egg-rr52.7%
Taylor expanded in im around 0
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/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
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6452.7%
Simplified52.7%
Taylor expanded in re around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6459.3%
Simplified59.3%
Final simplification85.5%
im\_m = (fabs.f64 im)
im\_s = (copysign.f64 #s(literal 1 binary64) im)
(FPCore (im_s re im_m)
:precision binary64
(*
im_s
(if (<= im_m 600.0)
(* (cos re) (- 0.0 im_m))
(if (<= im_m 5e+121)
(*
im_m
(*
(+ 1.0 (* re (* re (+ -0.5 (* (* re re) 0.041666666666666664)))))
(+
(+ -1.0 (* (* im_m im_m) -0.16666666666666666))
(*
(+ -0.008333333333333333 (* (* im_m im_m) -0.0001984126984126984))
(* im_m (* im_m (* im_m im_m)))))))
(*
im_m
(*
(+ 0.5 (* -0.25 (* re re)))
(+ -2.0 (* (* im_m im_m) -0.3333333333333333))))))))im\_m = fabs(im);
im\_s = copysign(1.0, im);
double code(double im_s, double re, double im_m) {
double tmp;
if (im_m <= 600.0) {
tmp = cos(re) * (0.0 - im_m);
} else if (im_m <= 5e+121) {
tmp = im_m * ((1.0 + (re * (re * (-0.5 + ((re * re) * 0.041666666666666664))))) * ((-1.0 + ((im_m * im_m) * -0.16666666666666666)) + ((-0.008333333333333333 + ((im_m * im_m) * -0.0001984126984126984)) * (im_m * (im_m * (im_m * im_m))))));
} else {
tmp = im_m * ((0.5 + (-0.25 * (re * re))) * (-2.0 + ((im_m * im_m) * -0.3333333333333333)));
}
return im_s * tmp;
}
im\_m = abs(im)
im\_s = copysign(1.0d0, im)
real(8) function code(im_s, re, im_m)
real(8), intent (in) :: im_s
real(8), intent (in) :: re
real(8), intent (in) :: im_m
real(8) :: tmp
if (im_m <= 600.0d0) then
tmp = cos(re) * (0.0d0 - im_m)
else if (im_m <= 5d+121) then
tmp = im_m * ((1.0d0 + (re * (re * ((-0.5d0) + ((re * re) * 0.041666666666666664d0))))) * (((-1.0d0) + ((im_m * im_m) * (-0.16666666666666666d0))) + (((-0.008333333333333333d0) + ((im_m * im_m) * (-0.0001984126984126984d0))) * (im_m * (im_m * (im_m * im_m))))))
else
tmp = im_m * ((0.5d0 + ((-0.25d0) * (re * re))) * ((-2.0d0) + ((im_m * im_m) * (-0.3333333333333333d0))))
end if
code = im_s * tmp
end function
im\_m = Math.abs(im);
im\_s = Math.copySign(1.0, im);
public static double code(double im_s, double re, double im_m) {
double tmp;
if (im_m <= 600.0) {
tmp = Math.cos(re) * (0.0 - im_m);
} else if (im_m <= 5e+121) {
tmp = im_m * ((1.0 + (re * (re * (-0.5 + ((re * re) * 0.041666666666666664))))) * ((-1.0 + ((im_m * im_m) * -0.16666666666666666)) + ((-0.008333333333333333 + ((im_m * im_m) * -0.0001984126984126984)) * (im_m * (im_m * (im_m * im_m))))));
} else {
tmp = im_m * ((0.5 + (-0.25 * (re * re))) * (-2.0 + ((im_m * im_m) * -0.3333333333333333)));
}
return im_s * tmp;
}
im\_m = math.fabs(im) im\_s = math.copysign(1.0, im) def code(im_s, re, im_m): tmp = 0 if im_m <= 600.0: tmp = math.cos(re) * (0.0 - im_m) elif im_m <= 5e+121: tmp = im_m * ((1.0 + (re * (re * (-0.5 + ((re * re) * 0.041666666666666664))))) * ((-1.0 + ((im_m * im_m) * -0.16666666666666666)) + ((-0.008333333333333333 + ((im_m * im_m) * -0.0001984126984126984)) * (im_m * (im_m * (im_m * im_m)))))) else: tmp = im_m * ((0.5 + (-0.25 * (re * re))) * (-2.0 + ((im_m * im_m) * -0.3333333333333333))) return im_s * tmp
im\_m = abs(im) im\_s = copysign(1.0, im) function code(im_s, re, im_m) tmp = 0.0 if (im_m <= 600.0) tmp = Float64(cos(re) * Float64(0.0 - im_m)); elseif (im_m <= 5e+121) tmp = Float64(im_m * Float64(Float64(1.0 + Float64(re * Float64(re * Float64(-0.5 + Float64(Float64(re * re) * 0.041666666666666664))))) * Float64(Float64(-1.0 + Float64(Float64(im_m * im_m) * -0.16666666666666666)) + Float64(Float64(-0.008333333333333333 + Float64(Float64(im_m * im_m) * -0.0001984126984126984)) * Float64(im_m * Float64(im_m * Float64(im_m * im_m))))))); else tmp = Float64(im_m * Float64(Float64(0.5 + Float64(-0.25 * Float64(re * re))) * Float64(-2.0 + Float64(Float64(im_m * im_m) * -0.3333333333333333)))); end return Float64(im_s * tmp) end
im\_m = abs(im); im\_s = sign(im) * abs(1.0); function tmp_2 = code(im_s, re, im_m) tmp = 0.0; if (im_m <= 600.0) tmp = cos(re) * (0.0 - im_m); elseif (im_m <= 5e+121) tmp = im_m * ((1.0 + (re * (re * (-0.5 + ((re * re) * 0.041666666666666664))))) * ((-1.0 + ((im_m * im_m) * -0.16666666666666666)) + ((-0.008333333333333333 + ((im_m * im_m) * -0.0001984126984126984)) * (im_m * (im_m * (im_m * im_m)))))); else tmp = im_m * ((0.5 + (-0.25 * (re * re))) * (-2.0 + ((im_m * im_m) * -0.3333333333333333))); end tmp_2 = im_s * tmp; end
im\_m = N[Abs[im], $MachinePrecision]
im\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[im]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[im$95$s_, re_, im$95$m_] := N[(im$95$s * If[LessEqual[im$95$m, 600.0], N[(N[Cos[re], $MachinePrecision] * N[(0.0 - im$95$m), $MachinePrecision]), $MachinePrecision], If[LessEqual[im$95$m, 5e+121], N[(im$95$m * N[(N[(1.0 + N[(re * N[(re * N[(-0.5 + N[(N[(re * re), $MachinePrecision] * 0.041666666666666664), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[(-1.0 + N[(N[(im$95$m * im$95$m), $MachinePrecision] * -0.16666666666666666), $MachinePrecision]), $MachinePrecision] + N[(N[(-0.008333333333333333 + N[(N[(im$95$m * im$95$m), $MachinePrecision] * -0.0001984126984126984), $MachinePrecision]), $MachinePrecision] * N[(im$95$m * N[(im$95$m * N[(im$95$m * im$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(im$95$m * N[(N[(0.5 + N[(-0.25 * N[(re * re), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(-2.0 + N[(N[(im$95$m * im$95$m), $MachinePrecision] * -0.3333333333333333), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]
\begin{array}{l}
im\_m = \left|im\right|
\\
im\_s = \mathsf{copysign}\left(1, im\right)
\\
im\_s \cdot \begin{array}{l}
\mathbf{if}\;im\_m \leq 600:\\
\;\;\;\;\cos re \cdot \left(0 - im\_m\right)\\
\mathbf{elif}\;im\_m \leq 5 \cdot 10^{+121}:\\
\;\;\;\;im\_m \cdot \left(\left(1 + re \cdot \left(re \cdot \left(-0.5 + \left(re \cdot re\right) \cdot 0.041666666666666664\right)\right)\right) \cdot \left(\left(-1 + \left(im\_m \cdot im\_m\right) \cdot -0.16666666666666666\right) + \left(-0.008333333333333333 + \left(im\_m \cdot im\_m\right) \cdot -0.0001984126984126984\right) \cdot \left(im\_m \cdot \left(im\_m \cdot \left(im\_m \cdot im\_m\right)\right)\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;im\_m \cdot \left(\left(0.5 + -0.25 \cdot \left(re \cdot re\right)\right) \cdot \left(-2 + \left(im\_m \cdot im\_m\right) \cdot -0.3333333333333333\right)\right)\\
\end{array}
\end{array}
if im < 600Initial program 46.7%
Taylor expanded in im around 0
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
cos-lowering-cos.f6460.8%
Simplified60.8%
sub0-negN/A
neg-lowering-neg.f64N/A
*-lowering-*.f64N/A
cos-lowering-cos.f6460.8%
Applied egg-rr60.8%
if 600 < im < 5.00000000000000007e121Initial program 100.0%
Taylor expanded in im around 0
*-lowering-*.f64N/A
distribute-rgt-inN/A
associate-+r+N/A
*-commutativeN/A
+-commutativeN/A
+-commutativeN/A
Simplified62.2%
Taylor expanded in re around 0
+-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
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6467.5%
Simplified67.5%
if 5.00000000000000007e121 < im Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
associate-*r*N/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
--lowering--.f64N/A
exp-negN/A
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
exp-lowering-exp.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6480.6%
Simplified80.6%
Taylor expanded in im around 0
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-commutativeN/A
associate-*r*N/A
*-commutativeN/A
distribute-rgt-outN/A
metadata-evalN/A
sub-negN/A
*-lowering-*.f64N/A
Simplified80.6%
Final simplification64.0%
im\_m = (fabs.f64 im)
im\_s = (copysign.f64 #s(literal 1 binary64) im)
(FPCore (im_s re im_m)
:precision binary64
(let* ((t_0
(+
-0.16666666666666666
(*
im_m
(*
im_m
(+
-0.008333333333333333
(* im_m (* im_m -0.0001984126984126984)))))))
(t_1 (* im_m (* im_m t_0))))
(*
im_s
(if (<= im_m 3.4e+24)
(/
(*
im_m
(+ -1.0 (* t_1 (* (* (* im_m im_m) (* im_m im_m)) (* t_0 t_0)))))
(+ 1.0 (* t_1 (- t_1 -1.0))))
(*
(+ 0.5 (* -0.25 (* re re)))
(*
im_m
(+
-2.0
(*
(* im_m im_m)
(+
-0.3333333333333333
(*
im_m
(*
im_m
(+
-0.016666666666666666
(* (* im_m im_m) -0.0003968253968253968)))))))))))))im\_m = fabs(im);
im\_s = copysign(1.0, im);
double code(double im_s, double re, double im_m) {
double t_0 = -0.16666666666666666 + (im_m * (im_m * (-0.008333333333333333 + (im_m * (im_m * -0.0001984126984126984)))));
double t_1 = im_m * (im_m * t_0);
double tmp;
if (im_m <= 3.4e+24) {
tmp = (im_m * (-1.0 + (t_1 * (((im_m * im_m) * (im_m * im_m)) * (t_0 * t_0))))) / (1.0 + (t_1 * (t_1 - -1.0)));
} else {
tmp = (0.5 + (-0.25 * (re * re))) * (im_m * (-2.0 + ((im_m * im_m) * (-0.3333333333333333 + (im_m * (im_m * (-0.016666666666666666 + ((im_m * im_m) * -0.0003968253968253968))))))));
}
return im_s * tmp;
}
im\_m = abs(im)
im\_s = copysign(1.0d0, im)
real(8) function code(im_s, re, im_m)
real(8), intent (in) :: im_s
real(8), intent (in) :: re
real(8), intent (in) :: im_m
real(8) :: t_0
real(8) :: t_1
real(8) :: tmp
t_0 = (-0.16666666666666666d0) + (im_m * (im_m * ((-0.008333333333333333d0) + (im_m * (im_m * (-0.0001984126984126984d0))))))
t_1 = im_m * (im_m * t_0)
if (im_m <= 3.4d+24) then
tmp = (im_m * ((-1.0d0) + (t_1 * (((im_m * im_m) * (im_m * im_m)) * (t_0 * t_0))))) / (1.0d0 + (t_1 * (t_1 - (-1.0d0))))
else
tmp = (0.5d0 + ((-0.25d0) * (re * re))) * (im_m * ((-2.0d0) + ((im_m * im_m) * ((-0.3333333333333333d0) + (im_m * (im_m * ((-0.016666666666666666d0) + ((im_m * im_m) * (-0.0003968253968253968d0)))))))))
end if
code = im_s * tmp
end function
im\_m = Math.abs(im);
im\_s = Math.copySign(1.0, im);
public static double code(double im_s, double re, double im_m) {
double t_0 = -0.16666666666666666 + (im_m * (im_m * (-0.008333333333333333 + (im_m * (im_m * -0.0001984126984126984)))));
double t_1 = im_m * (im_m * t_0);
double tmp;
if (im_m <= 3.4e+24) {
tmp = (im_m * (-1.0 + (t_1 * (((im_m * im_m) * (im_m * im_m)) * (t_0 * t_0))))) / (1.0 + (t_1 * (t_1 - -1.0)));
} else {
tmp = (0.5 + (-0.25 * (re * re))) * (im_m * (-2.0 + ((im_m * im_m) * (-0.3333333333333333 + (im_m * (im_m * (-0.016666666666666666 + ((im_m * im_m) * -0.0003968253968253968))))))));
}
return im_s * tmp;
}
im\_m = math.fabs(im) im\_s = math.copysign(1.0, im) def code(im_s, re, im_m): t_0 = -0.16666666666666666 + (im_m * (im_m * (-0.008333333333333333 + (im_m * (im_m * -0.0001984126984126984))))) t_1 = im_m * (im_m * t_0) tmp = 0 if im_m <= 3.4e+24: tmp = (im_m * (-1.0 + (t_1 * (((im_m * im_m) * (im_m * im_m)) * (t_0 * t_0))))) / (1.0 + (t_1 * (t_1 - -1.0))) else: tmp = (0.5 + (-0.25 * (re * re))) * (im_m * (-2.0 + ((im_m * im_m) * (-0.3333333333333333 + (im_m * (im_m * (-0.016666666666666666 + ((im_m * im_m) * -0.0003968253968253968)))))))) return im_s * tmp
im\_m = abs(im) im\_s = copysign(1.0, im) function code(im_s, re, im_m) t_0 = Float64(-0.16666666666666666 + Float64(im_m * Float64(im_m * Float64(-0.008333333333333333 + Float64(im_m * Float64(im_m * -0.0001984126984126984)))))) t_1 = Float64(im_m * Float64(im_m * t_0)) tmp = 0.0 if (im_m <= 3.4e+24) tmp = Float64(Float64(im_m * Float64(-1.0 + Float64(t_1 * Float64(Float64(Float64(im_m * im_m) * Float64(im_m * im_m)) * Float64(t_0 * t_0))))) / Float64(1.0 + Float64(t_1 * Float64(t_1 - -1.0)))); else tmp = Float64(Float64(0.5 + Float64(-0.25 * Float64(re * re))) * Float64(im_m * Float64(-2.0 + Float64(Float64(im_m * im_m) * Float64(-0.3333333333333333 + Float64(im_m * Float64(im_m * Float64(-0.016666666666666666 + Float64(Float64(im_m * im_m) * -0.0003968253968253968))))))))); end return Float64(im_s * tmp) end
im\_m = abs(im); im\_s = sign(im) * abs(1.0); function tmp_2 = code(im_s, re, im_m) t_0 = -0.16666666666666666 + (im_m * (im_m * (-0.008333333333333333 + (im_m * (im_m * -0.0001984126984126984))))); t_1 = im_m * (im_m * t_0); tmp = 0.0; if (im_m <= 3.4e+24) tmp = (im_m * (-1.0 + (t_1 * (((im_m * im_m) * (im_m * im_m)) * (t_0 * t_0))))) / (1.0 + (t_1 * (t_1 - -1.0))); else tmp = (0.5 + (-0.25 * (re * re))) * (im_m * (-2.0 + ((im_m * im_m) * (-0.3333333333333333 + (im_m * (im_m * (-0.016666666666666666 + ((im_m * im_m) * -0.0003968253968253968)))))))); end tmp_2 = im_s * tmp; end
im\_m = N[Abs[im], $MachinePrecision]
im\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[im]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[im$95$s_, re_, im$95$m_] := Block[{t$95$0 = N[(-0.16666666666666666 + N[(im$95$m * N[(im$95$m * N[(-0.008333333333333333 + N[(im$95$m * N[(im$95$m * -0.0001984126984126984), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(im$95$m * N[(im$95$m * t$95$0), $MachinePrecision]), $MachinePrecision]}, N[(im$95$s * If[LessEqual[im$95$m, 3.4e+24], N[(N[(im$95$m * N[(-1.0 + N[(t$95$1 * N[(N[(N[(im$95$m * im$95$m), $MachinePrecision] * N[(im$95$m * im$95$m), $MachinePrecision]), $MachinePrecision] * N[(t$95$0 * t$95$0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(1.0 + N[(t$95$1 * N[(t$95$1 - -1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(0.5 + N[(-0.25 * N[(re * re), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(im$95$m * N[(-2.0 + N[(N[(im$95$m * im$95$m), $MachinePrecision] * N[(-0.3333333333333333 + N[(im$95$m * N[(im$95$m * N[(-0.016666666666666666 + N[(N[(im$95$m * im$95$m), $MachinePrecision] * -0.0003968253968253968), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]]]
\begin{array}{l}
im\_m = \left|im\right|
\\
im\_s = \mathsf{copysign}\left(1, im\right)
\\
\begin{array}{l}
t_0 := -0.16666666666666666 + im\_m \cdot \left(im\_m \cdot \left(-0.008333333333333333 + im\_m \cdot \left(im\_m \cdot -0.0001984126984126984\right)\right)\right)\\
t_1 := im\_m \cdot \left(im\_m \cdot t\_0\right)\\
im\_s \cdot \begin{array}{l}
\mathbf{if}\;im\_m \leq 3.4 \cdot 10^{+24}:\\
\;\;\;\;\frac{im\_m \cdot \left(-1 + t\_1 \cdot \left(\left(\left(im\_m \cdot im\_m\right) \cdot \left(im\_m \cdot im\_m\right)\right) \cdot \left(t\_0 \cdot t\_0\right)\right)\right)}{1 + t\_1 \cdot \left(t\_1 - -1\right)}\\
\mathbf{else}:\\
\;\;\;\;\left(0.5 + -0.25 \cdot \left(re \cdot re\right)\right) \cdot \left(im\_m \cdot \left(-2 + \left(im\_m \cdot im\_m\right) \cdot \left(-0.3333333333333333 + im\_m \cdot \left(im\_m \cdot \left(-0.016666666666666666 + \left(im\_m \cdot im\_m\right) \cdot -0.0003968253968253968\right)\right)\right)\right)\right)\\
\end{array}
\end{array}
\end{array}
if im < 3.4000000000000001e24Initial program 48.5%
Taylor expanded in im around 0
*-lowering-*.f64N/A
distribute-rgt-inN/A
associate-+r+N/A
*-commutativeN/A
+-commutativeN/A
+-commutativeN/A
Simplified89.3%
Taylor expanded in re around 0
Simplified58.3%
*-lft-identityN/A
flip-+N/A
div-subN/A
--lowering--.f64N/A
Applied egg-rr36.4%
Applied egg-rr36.3%
if 3.4000000000000001e24 < im Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
associate-*r*N/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
--lowering--.f64N/A
exp-negN/A
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
exp-lowering-exp.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6477.8%
Simplified77.8%
Taylor expanded in im around 0
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/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
Simplified76.0%
Final simplification44.6%
im\_m = (fabs.f64 im)
im\_s = (copysign.f64 #s(literal 1 binary64) im)
(FPCore (im_s re im_m)
:precision binary64
(let* ((t_0
(+
-0.16666666666666666
(*
im_m
(*
im_m
(+
-0.008333333333333333
(* im_m (* im_m -0.0001984126984126984))))))))
(*
im_s
(if (<= im_m 3.4e+39)
(/
(* im_m (- 1.0 (* (* (* im_m im_m) (* im_m im_m)) (* t_0 t_0))))
(- -1.0 (* im_m (* im_m t_0))))
(*
(+ 0.5 (* -0.25 (* re re)))
(*
im_m
(+
-2.0
(*
(* im_m im_m)
(+
-0.3333333333333333
(*
im_m
(*
im_m
(+
-0.016666666666666666
(* (* im_m im_m) -0.0003968253968253968)))))))))))))im\_m = fabs(im);
im\_s = copysign(1.0, im);
double code(double im_s, double re, double im_m) {
double t_0 = -0.16666666666666666 + (im_m * (im_m * (-0.008333333333333333 + (im_m * (im_m * -0.0001984126984126984)))));
double tmp;
if (im_m <= 3.4e+39) {
tmp = (im_m * (1.0 - (((im_m * im_m) * (im_m * im_m)) * (t_0 * t_0)))) / (-1.0 - (im_m * (im_m * t_0)));
} else {
tmp = (0.5 + (-0.25 * (re * re))) * (im_m * (-2.0 + ((im_m * im_m) * (-0.3333333333333333 + (im_m * (im_m * (-0.016666666666666666 + ((im_m * im_m) * -0.0003968253968253968))))))));
}
return im_s * tmp;
}
im\_m = abs(im)
im\_s = copysign(1.0d0, im)
real(8) function code(im_s, re, im_m)
real(8), intent (in) :: im_s
real(8), intent (in) :: re
real(8), intent (in) :: im_m
real(8) :: t_0
real(8) :: tmp
t_0 = (-0.16666666666666666d0) + (im_m * (im_m * ((-0.008333333333333333d0) + (im_m * (im_m * (-0.0001984126984126984d0))))))
if (im_m <= 3.4d+39) then
tmp = (im_m * (1.0d0 - (((im_m * im_m) * (im_m * im_m)) * (t_0 * t_0)))) / ((-1.0d0) - (im_m * (im_m * t_0)))
else
tmp = (0.5d0 + ((-0.25d0) * (re * re))) * (im_m * ((-2.0d0) + ((im_m * im_m) * ((-0.3333333333333333d0) + (im_m * (im_m * ((-0.016666666666666666d0) + ((im_m * im_m) * (-0.0003968253968253968d0)))))))))
end if
code = im_s * tmp
end function
im\_m = Math.abs(im);
im\_s = Math.copySign(1.0, im);
public static double code(double im_s, double re, double im_m) {
double t_0 = -0.16666666666666666 + (im_m * (im_m * (-0.008333333333333333 + (im_m * (im_m * -0.0001984126984126984)))));
double tmp;
if (im_m <= 3.4e+39) {
tmp = (im_m * (1.0 - (((im_m * im_m) * (im_m * im_m)) * (t_0 * t_0)))) / (-1.0 - (im_m * (im_m * t_0)));
} else {
tmp = (0.5 + (-0.25 * (re * re))) * (im_m * (-2.0 + ((im_m * im_m) * (-0.3333333333333333 + (im_m * (im_m * (-0.016666666666666666 + ((im_m * im_m) * -0.0003968253968253968))))))));
}
return im_s * tmp;
}
im\_m = math.fabs(im) im\_s = math.copysign(1.0, im) def code(im_s, re, im_m): t_0 = -0.16666666666666666 + (im_m * (im_m * (-0.008333333333333333 + (im_m * (im_m * -0.0001984126984126984))))) tmp = 0 if im_m <= 3.4e+39: tmp = (im_m * (1.0 - (((im_m * im_m) * (im_m * im_m)) * (t_0 * t_0)))) / (-1.0 - (im_m * (im_m * t_0))) else: tmp = (0.5 + (-0.25 * (re * re))) * (im_m * (-2.0 + ((im_m * im_m) * (-0.3333333333333333 + (im_m * (im_m * (-0.016666666666666666 + ((im_m * im_m) * -0.0003968253968253968)))))))) return im_s * tmp
im\_m = abs(im) im\_s = copysign(1.0, im) function code(im_s, re, im_m) t_0 = Float64(-0.16666666666666666 + Float64(im_m * Float64(im_m * Float64(-0.008333333333333333 + Float64(im_m * Float64(im_m * -0.0001984126984126984)))))) tmp = 0.0 if (im_m <= 3.4e+39) tmp = Float64(Float64(im_m * Float64(1.0 - Float64(Float64(Float64(im_m * im_m) * Float64(im_m * im_m)) * Float64(t_0 * t_0)))) / Float64(-1.0 - Float64(im_m * Float64(im_m * t_0)))); else tmp = Float64(Float64(0.5 + Float64(-0.25 * Float64(re * re))) * Float64(im_m * Float64(-2.0 + Float64(Float64(im_m * im_m) * Float64(-0.3333333333333333 + Float64(im_m * Float64(im_m * Float64(-0.016666666666666666 + Float64(Float64(im_m * im_m) * -0.0003968253968253968))))))))); end return Float64(im_s * tmp) end
im\_m = abs(im); im\_s = sign(im) * abs(1.0); function tmp_2 = code(im_s, re, im_m) t_0 = -0.16666666666666666 + (im_m * (im_m * (-0.008333333333333333 + (im_m * (im_m * -0.0001984126984126984))))); tmp = 0.0; if (im_m <= 3.4e+39) tmp = (im_m * (1.0 - (((im_m * im_m) * (im_m * im_m)) * (t_0 * t_0)))) / (-1.0 - (im_m * (im_m * t_0))); else tmp = (0.5 + (-0.25 * (re * re))) * (im_m * (-2.0 + ((im_m * im_m) * (-0.3333333333333333 + (im_m * (im_m * (-0.016666666666666666 + ((im_m * im_m) * -0.0003968253968253968)))))))); end tmp_2 = im_s * tmp; end
im\_m = N[Abs[im], $MachinePrecision]
im\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[im]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[im$95$s_, re_, im$95$m_] := Block[{t$95$0 = N[(-0.16666666666666666 + N[(im$95$m * N[(im$95$m * N[(-0.008333333333333333 + N[(im$95$m * N[(im$95$m * -0.0001984126984126984), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, N[(im$95$s * If[LessEqual[im$95$m, 3.4e+39], N[(N[(im$95$m * N[(1.0 - N[(N[(N[(im$95$m * im$95$m), $MachinePrecision] * N[(im$95$m * im$95$m), $MachinePrecision]), $MachinePrecision] * N[(t$95$0 * t$95$0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(-1.0 - N[(im$95$m * N[(im$95$m * t$95$0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(0.5 + N[(-0.25 * N[(re * re), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(im$95$m * N[(-2.0 + N[(N[(im$95$m * im$95$m), $MachinePrecision] * N[(-0.3333333333333333 + N[(im$95$m * N[(im$95$m * N[(-0.016666666666666666 + N[(N[(im$95$m * im$95$m), $MachinePrecision] * -0.0003968253968253968), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]]
\begin{array}{l}
im\_m = \left|im\right|
\\
im\_s = \mathsf{copysign}\left(1, im\right)
\\
\begin{array}{l}
t_0 := -0.16666666666666666 + im\_m \cdot \left(im\_m \cdot \left(-0.008333333333333333 + im\_m \cdot \left(im\_m \cdot -0.0001984126984126984\right)\right)\right)\\
im\_s \cdot \begin{array}{l}
\mathbf{if}\;im\_m \leq 3.4 \cdot 10^{+39}:\\
\;\;\;\;\frac{im\_m \cdot \left(1 - \left(\left(im\_m \cdot im\_m\right) \cdot \left(im\_m \cdot im\_m\right)\right) \cdot \left(t\_0 \cdot t\_0\right)\right)}{-1 - im\_m \cdot \left(im\_m \cdot t\_0\right)}\\
\mathbf{else}:\\
\;\;\;\;\left(0.5 + -0.25 \cdot \left(re \cdot re\right)\right) \cdot \left(im\_m \cdot \left(-2 + \left(im\_m \cdot im\_m\right) \cdot \left(-0.3333333333333333 + im\_m \cdot \left(im\_m \cdot \left(-0.016666666666666666 + \left(im\_m \cdot im\_m\right) \cdot -0.0003968253968253968\right)\right)\right)\right)\right)\\
\end{array}
\end{array}
\end{array}
if im < 3.3999999999999999e39Initial program 49.5%
Taylor expanded in im around 0
*-lowering-*.f64N/A
distribute-rgt-inN/A
associate-+r+N/A
*-commutativeN/A
+-commutativeN/A
+-commutativeN/A
Simplified87.7%
Taylor expanded in re around 0
Simplified57.2%
*-lft-identityN/A
flip-+N/A
div-subN/A
--lowering--.f64N/A
Applied egg-rr36.7%
Applied egg-rr37.1%
if 3.3999999999999999e39 < im Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
associate-*r*N/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
--lowering--.f64N/A
exp-negN/A
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
exp-lowering-exp.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6478.0%
Simplified78.0%
Taylor expanded in im around 0
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/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
Simplified78.0%
Final simplification45.1%
im\_m = (fabs.f64 im)
im\_s = (copysign.f64 #s(literal 1 binary64) im)
(FPCore (im_s re im_m)
:precision binary64
(let* ((t_0
(+
-1.0
(*
(* im_m im_m)
(+
-0.16666666666666666
(*
im_m
(*
im_m
(+
-0.008333333333333333
(* (* im_m im_m) -0.0001984126984126984)))))))))
(*
im_s
(if (<= im_m 480.0)
(* im_m t_0)
(if (<= im_m 2e+107)
(*
t_0
(+
im_m
(*
(* re re)
(+ (* im_m -0.5) (* 0.041666666666666664 (* im_m (* re re)))))))
(*
im_m
(*
(+ 0.5 (* -0.25 (* re re)))
(+ -2.0 (* (* im_m im_m) -0.3333333333333333)))))))))im\_m = fabs(im);
im\_s = copysign(1.0, im);
double code(double im_s, double re, double im_m) {
double t_0 = -1.0 + ((im_m * im_m) * (-0.16666666666666666 + (im_m * (im_m * (-0.008333333333333333 + ((im_m * im_m) * -0.0001984126984126984))))));
double tmp;
if (im_m <= 480.0) {
tmp = im_m * t_0;
} else if (im_m <= 2e+107) {
tmp = t_0 * (im_m + ((re * re) * ((im_m * -0.5) + (0.041666666666666664 * (im_m * (re * re))))));
} else {
tmp = im_m * ((0.5 + (-0.25 * (re * re))) * (-2.0 + ((im_m * im_m) * -0.3333333333333333)));
}
return im_s * tmp;
}
im\_m = abs(im)
im\_s = copysign(1.0d0, im)
real(8) function code(im_s, re, im_m)
real(8), intent (in) :: im_s
real(8), intent (in) :: re
real(8), intent (in) :: im_m
real(8) :: t_0
real(8) :: tmp
t_0 = (-1.0d0) + ((im_m * im_m) * ((-0.16666666666666666d0) + (im_m * (im_m * ((-0.008333333333333333d0) + ((im_m * im_m) * (-0.0001984126984126984d0)))))))
if (im_m <= 480.0d0) then
tmp = im_m * t_0
else if (im_m <= 2d+107) then
tmp = t_0 * (im_m + ((re * re) * ((im_m * (-0.5d0)) + (0.041666666666666664d0 * (im_m * (re * re))))))
else
tmp = im_m * ((0.5d0 + ((-0.25d0) * (re * re))) * ((-2.0d0) + ((im_m * im_m) * (-0.3333333333333333d0))))
end if
code = im_s * tmp
end function
im\_m = Math.abs(im);
im\_s = Math.copySign(1.0, im);
public static double code(double im_s, double re, double im_m) {
double t_0 = -1.0 + ((im_m * im_m) * (-0.16666666666666666 + (im_m * (im_m * (-0.008333333333333333 + ((im_m * im_m) * -0.0001984126984126984))))));
double tmp;
if (im_m <= 480.0) {
tmp = im_m * t_0;
} else if (im_m <= 2e+107) {
tmp = t_0 * (im_m + ((re * re) * ((im_m * -0.5) + (0.041666666666666664 * (im_m * (re * re))))));
} else {
tmp = im_m * ((0.5 + (-0.25 * (re * re))) * (-2.0 + ((im_m * im_m) * -0.3333333333333333)));
}
return im_s * tmp;
}
im\_m = math.fabs(im) im\_s = math.copysign(1.0, im) def code(im_s, re, im_m): t_0 = -1.0 + ((im_m * im_m) * (-0.16666666666666666 + (im_m * (im_m * (-0.008333333333333333 + ((im_m * im_m) * -0.0001984126984126984)))))) tmp = 0 if im_m <= 480.0: tmp = im_m * t_0 elif im_m <= 2e+107: tmp = t_0 * (im_m + ((re * re) * ((im_m * -0.5) + (0.041666666666666664 * (im_m * (re * re)))))) else: tmp = im_m * ((0.5 + (-0.25 * (re * re))) * (-2.0 + ((im_m * im_m) * -0.3333333333333333))) return im_s * tmp
im\_m = abs(im) im\_s = copysign(1.0, im) function code(im_s, re, im_m) t_0 = Float64(-1.0 + Float64(Float64(im_m * im_m) * Float64(-0.16666666666666666 + Float64(im_m * Float64(im_m * Float64(-0.008333333333333333 + Float64(Float64(im_m * im_m) * -0.0001984126984126984))))))) tmp = 0.0 if (im_m <= 480.0) tmp = Float64(im_m * t_0); elseif (im_m <= 2e+107) tmp = Float64(t_0 * Float64(im_m + Float64(Float64(re * re) * Float64(Float64(im_m * -0.5) + Float64(0.041666666666666664 * Float64(im_m * Float64(re * re))))))); else tmp = Float64(im_m * Float64(Float64(0.5 + Float64(-0.25 * Float64(re * re))) * Float64(-2.0 + Float64(Float64(im_m * im_m) * -0.3333333333333333)))); end return Float64(im_s * tmp) end
im\_m = abs(im); im\_s = sign(im) * abs(1.0); function tmp_2 = code(im_s, re, im_m) t_0 = -1.0 + ((im_m * im_m) * (-0.16666666666666666 + (im_m * (im_m * (-0.008333333333333333 + ((im_m * im_m) * -0.0001984126984126984)))))); tmp = 0.0; if (im_m <= 480.0) tmp = im_m * t_0; elseif (im_m <= 2e+107) tmp = t_0 * (im_m + ((re * re) * ((im_m * -0.5) + (0.041666666666666664 * (im_m * (re * re)))))); else tmp = im_m * ((0.5 + (-0.25 * (re * re))) * (-2.0 + ((im_m * im_m) * -0.3333333333333333))); end tmp_2 = im_s * tmp; end
im\_m = N[Abs[im], $MachinePrecision]
im\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[im]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[im$95$s_, re_, im$95$m_] := Block[{t$95$0 = N[(-1.0 + N[(N[(im$95$m * im$95$m), $MachinePrecision] * N[(-0.16666666666666666 + N[(im$95$m * N[(im$95$m * N[(-0.008333333333333333 + N[(N[(im$95$m * im$95$m), $MachinePrecision] * -0.0001984126984126984), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, N[(im$95$s * If[LessEqual[im$95$m, 480.0], N[(im$95$m * t$95$0), $MachinePrecision], If[LessEqual[im$95$m, 2e+107], N[(t$95$0 * N[(im$95$m + N[(N[(re * re), $MachinePrecision] * N[(N[(im$95$m * -0.5), $MachinePrecision] + N[(0.041666666666666664 * N[(im$95$m * N[(re * re), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(im$95$m * N[(N[(0.5 + N[(-0.25 * N[(re * re), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(-2.0 + N[(N[(im$95$m * im$95$m), $MachinePrecision] * -0.3333333333333333), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]]
\begin{array}{l}
im\_m = \left|im\right|
\\
im\_s = \mathsf{copysign}\left(1, im\right)
\\
\begin{array}{l}
t_0 := -1 + \left(im\_m \cdot im\_m\right) \cdot \left(-0.16666666666666666 + im\_m \cdot \left(im\_m \cdot \left(-0.008333333333333333 + \left(im\_m \cdot im\_m\right) \cdot -0.0001984126984126984\right)\right)\right)\\
im\_s \cdot \begin{array}{l}
\mathbf{if}\;im\_m \leq 480:\\
\;\;\;\;im\_m \cdot t\_0\\
\mathbf{elif}\;im\_m \leq 2 \cdot 10^{+107}:\\
\;\;\;\;t\_0 \cdot \left(im\_m + \left(re \cdot re\right) \cdot \left(im\_m \cdot -0.5 + 0.041666666666666664 \cdot \left(im\_m \cdot \left(re \cdot re\right)\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;im\_m \cdot \left(\left(0.5 + -0.25 \cdot \left(re \cdot re\right)\right) \cdot \left(-2 + \left(im\_m \cdot im\_m\right) \cdot -0.3333333333333333\right)\right)\\
\end{array}
\end{array}
\end{array}
if im < 480Initial program 46.7%
Taylor expanded in im around 0
*-lowering-*.f64N/A
distribute-rgt-inN/A
associate-+r+N/A
*-commutativeN/A
+-commutativeN/A
+-commutativeN/A
Simplified92.4%
Taylor expanded in re around 0
Simplified60.2%
Taylor expanded in im around 0
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/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
Simplified60.2%
if 480 < im < 1.9999999999999999e107Initial program 100.0%
Taylor expanded in im around 0
*-lowering-*.f64N/A
distribute-rgt-inN/A
associate-+r+N/A
*-commutativeN/A
+-commutativeN/A
+-commutativeN/A
Simplified59.5%
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
Applied egg-rr59.5%
Taylor expanded in im around 0
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/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
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6459.5%
Simplified59.5%
Taylor expanded in re around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6465.1%
Simplified65.1%
if 1.9999999999999999e107 < im Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
associate-*r*N/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
--lowering--.f64N/A
exp-negN/A
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
exp-lowering-exp.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6478.8%
Simplified78.8%
Taylor expanded in im around 0
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-commutativeN/A
associate-*r*N/A
*-commutativeN/A
distribute-rgt-outN/A
metadata-evalN/A
sub-negN/A
*-lowering-*.f64N/A
Simplified78.8%
Final simplification63.2%
im\_m = (fabs.f64 im)
im\_s = (copysign.f64 #s(literal 1 binary64) im)
(FPCore (im_s re im_m)
:precision binary64
(*
im_s
(if (<= im_m 430.0)
(*
im_m
(+
-1.0
(*
(* im_m im_m)
(+
-0.16666666666666666
(*
im_m
(*
im_m
(+
-0.008333333333333333
(* (* im_m im_m) -0.0001984126984126984))))))))
(if (<= im_m 2.6e+24)
(*
(+
-1.0
(*
(* im_m im_m)
(+ -0.16666666666666666 (* (* im_m im_m) -0.008333333333333333))))
(+
im_m
(* (* im_m (* re re)) (+ -0.5 (* (* re re) 0.041666666666666664)))))
(*
(+ 0.5 (* -0.25 (* re re)))
(*
im_m
(+
-2.0
(*
(* im_m im_m)
(+
-0.3333333333333333
(*
im_m
(*
im_m
(+
-0.016666666666666666
(* (* im_m im_m) -0.0003968253968253968)))))))))))))im\_m = fabs(im);
im\_s = copysign(1.0, im);
double code(double im_s, double re, double im_m) {
double tmp;
if (im_m <= 430.0) {
tmp = im_m * (-1.0 + ((im_m * im_m) * (-0.16666666666666666 + (im_m * (im_m * (-0.008333333333333333 + ((im_m * im_m) * -0.0001984126984126984)))))));
} else if (im_m <= 2.6e+24) {
tmp = (-1.0 + ((im_m * im_m) * (-0.16666666666666666 + ((im_m * im_m) * -0.008333333333333333)))) * (im_m + ((im_m * (re * re)) * (-0.5 + ((re * re) * 0.041666666666666664))));
} else {
tmp = (0.5 + (-0.25 * (re * re))) * (im_m * (-2.0 + ((im_m * im_m) * (-0.3333333333333333 + (im_m * (im_m * (-0.016666666666666666 + ((im_m * im_m) * -0.0003968253968253968))))))));
}
return im_s * tmp;
}
im\_m = abs(im)
im\_s = copysign(1.0d0, im)
real(8) function code(im_s, re, im_m)
real(8), intent (in) :: im_s
real(8), intent (in) :: re
real(8), intent (in) :: im_m
real(8) :: tmp
if (im_m <= 430.0d0) then
tmp = im_m * ((-1.0d0) + ((im_m * im_m) * ((-0.16666666666666666d0) + (im_m * (im_m * ((-0.008333333333333333d0) + ((im_m * im_m) * (-0.0001984126984126984d0))))))))
else if (im_m <= 2.6d+24) then
tmp = ((-1.0d0) + ((im_m * im_m) * ((-0.16666666666666666d0) + ((im_m * im_m) * (-0.008333333333333333d0))))) * (im_m + ((im_m * (re * re)) * ((-0.5d0) + ((re * re) * 0.041666666666666664d0))))
else
tmp = (0.5d0 + ((-0.25d0) * (re * re))) * (im_m * ((-2.0d0) + ((im_m * im_m) * ((-0.3333333333333333d0) + (im_m * (im_m * ((-0.016666666666666666d0) + ((im_m * im_m) * (-0.0003968253968253968d0)))))))))
end if
code = im_s * tmp
end function
im\_m = Math.abs(im);
im\_s = Math.copySign(1.0, im);
public static double code(double im_s, double re, double im_m) {
double tmp;
if (im_m <= 430.0) {
tmp = im_m * (-1.0 + ((im_m * im_m) * (-0.16666666666666666 + (im_m * (im_m * (-0.008333333333333333 + ((im_m * im_m) * -0.0001984126984126984)))))));
} else if (im_m <= 2.6e+24) {
tmp = (-1.0 + ((im_m * im_m) * (-0.16666666666666666 + ((im_m * im_m) * -0.008333333333333333)))) * (im_m + ((im_m * (re * re)) * (-0.5 + ((re * re) * 0.041666666666666664))));
} else {
tmp = (0.5 + (-0.25 * (re * re))) * (im_m * (-2.0 + ((im_m * im_m) * (-0.3333333333333333 + (im_m * (im_m * (-0.016666666666666666 + ((im_m * im_m) * -0.0003968253968253968))))))));
}
return im_s * tmp;
}
im\_m = math.fabs(im) im\_s = math.copysign(1.0, im) def code(im_s, re, im_m): tmp = 0 if im_m <= 430.0: tmp = im_m * (-1.0 + ((im_m * im_m) * (-0.16666666666666666 + (im_m * (im_m * (-0.008333333333333333 + ((im_m * im_m) * -0.0001984126984126984))))))) elif im_m <= 2.6e+24: tmp = (-1.0 + ((im_m * im_m) * (-0.16666666666666666 + ((im_m * im_m) * -0.008333333333333333)))) * (im_m + ((im_m * (re * re)) * (-0.5 + ((re * re) * 0.041666666666666664)))) else: tmp = (0.5 + (-0.25 * (re * re))) * (im_m * (-2.0 + ((im_m * im_m) * (-0.3333333333333333 + (im_m * (im_m * (-0.016666666666666666 + ((im_m * im_m) * -0.0003968253968253968)))))))) return im_s * tmp
im\_m = abs(im) im\_s = copysign(1.0, im) function code(im_s, re, im_m) tmp = 0.0 if (im_m <= 430.0) tmp = Float64(im_m * Float64(-1.0 + Float64(Float64(im_m * im_m) * Float64(-0.16666666666666666 + Float64(im_m * Float64(im_m * Float64(-0.008333333333333333 + Float64(Float64(im_m * im_m) * -0.0001984126984126984)))))))); elseif (im_m <= 2.6e+24) tmp = Float64(Float64(-1.0 + Float64(Float64(im_m * im_m) * Float64(-0.16666666666666666 + Float64(Float64(im_m * im_m) * -0.008333333333333333)))) * Float64(im_m + Float64(Float64(im_m * Float64(re * re)) * Float64(-0.5 + Float64(Float64(re * re) * 0.041666666666666664))))); else tmp = Float64(Float64(0.5 + Float64(-0.25 * Float64(re * re))) * Float64(im_m * Float64(-2.0 + Float64(Float64(im_m * im_m) * Float64(-0.3333333333333333 + Float64(im_m * Float64(im_m * Float64(-0.016666666666666666 + Float64(Float64(im_m * im_m) * -0.0003968253968253968))))))))); end return Float64(im_s * tmp) end
im\_m = abs(im); im\_s = sign(im) * abs(1.0); function tmp_2 = code(im_s, re, im_m) tmp = 0.0; if (im_m <= 430.0) tmp = im_m * (-1.0 + ((im_m * im_m) * (-0.16666666666666666 + (im_m * (im_m * (-0.008333333333333333 + ((im_m * im_m) * -0.0001984126984126984))))))); elseif (im_m <= 2.6e+24) tmp = (-1.0 + ((im_m * im_m) * (-0.16666666666666666 + ((im_m * im_m) * -0.008333333333333333)))) * (im_m + ((im_m * (re * re)) * (-0.5 + ((re * re) * 0.041666666666666664)))); else tmp = (0.5 + (-0.25 * (re * re))) * (im_m * (-2.0 + ((im_m * im_m) * (-0.3333333333333333 + (im_m * (im_m * (-0.016666666666666666 + ((im_m * im_m) * -0.0003968253968253968)))))))); end tmp_2 = im_s * tmp; end
im\_m = N[Abs[im], $MachinePrecision]
im\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[im]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[im$95$s_, re_, im$95$m_] := N[(im$95$s * If[LessEqual[im$95$m, 430.0], N[(im$95$m * N[(-1.0 + N[(N[(im$95$m * im$95$m), $MachinePrecision] * N[(-0.16666666666666666 + N[(im$95$m * N[(im$95$m * N[(-0.008333333333333333 + N[(N[(im$95$m * im$95$m), $MachinePrecision] * -0.0001984126984126984), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[im$95$m, 2.6e+24], N[(N[(-1.0 + N[(N[(im$95$m * im$95$m), $MachinePrecision] * N[(-0.16666666666666666 + N[(N[(im$95$m * im$95$m), $MachinePrecision] * -0.008333333333333333), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(im$95$m + N[(N[(im$95$m * N[(re * re), $MachinePrecision]), $MachinePrecision] * N[(-0.5 + N[(N[(re * re), $MachinePrecision] * 0.041666666666666664), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(0.5 + N[(-0.25 * N[(re * re), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(im$95$m * N[(-2.0 + N[(N[(im$95$m * im$95$m), $MachinePrecision] * N[(-0.3333333333333333 + N[(im$95$m * N[(im$95$m * N[(-0.016666666666666666 + N[(N[(im$95$m * im$95$m), $MachinePrecision] * -0.0003968253968253968), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]
\begin{array}{l}
im\_m = \left|im\right|
\\
im\_s = \mathsf{copysign}\left(1, im\right)
\\
im\_s \cdot \begin{array}{l}
\mathbf{if}\;im\_m \leq 430:\\
\;\;\;\;im\_m \cdot \left(-1 + \left(im\_m \cdot im\_m\right) \cdot \left(-0.16666666666666666 + im\_m \cdot \left(im\_m \cdot \left(-0.008333333333333333 + \left(im\_m \cdot im\_m\right) \cdot -0.0001984126984126984\right)\right)\right)\right)\\
\mathbf{elif}\;im\_m \leq 2.6 \cdot 10^{+24}:\\
\;\;\;\;\left(-1 + \left(im\_m \cdot im\_m\right) \cdot \left(-0.16666666666666666 + \left(im\_m \cdot im\_m\right) \cdot -0.008333333333333333\right)\right) \cdot \left(im\_m + \left(im\_m \cdot \left(re \cdot re\right)\right) \cdot \left(-0.5 + \left(re \cdot re\right) \cdot 0.041666666666666664\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\left(0.5 + -0.25 \cdot \left(re \cdot re\right)\right) \cdot \left(im\_m \cdot \left(-2 + \left(im\_m \cdot im\_m\right) \cdot \left(-0.3333333333333333 + im\_m \cdot \left(im\_m \cdot \left(-0.016666666666666666 + \left(im\_m \cdot im\_m\right) \cdot -0.0003968253968253968\right)\right)\right)\right)\right)\\
\end{array}
\end{array}
if im < 430Initial program 46.7%
Taylor expanded in im around 0
*-lowering-*.f64N/A
distribute-rgt-inN/A
associate-+r+N/A
*-commutativeN/A
+-commutativeN/A
+-commutativeN/A
Simplified92.4%
Taylor expanded in re around 0
Simplified60.2%
Taylor expanded in im around 0
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/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
Simplified60.2%
if 430 < im < 2.5999999999999998e24Initial program 100.0%
Taylor expanded in im around 0
*-lowering-*.f64N/A
*-commutativeN/A
+-commutativeN/A
distribute-lft-inN/A
associate-*r*N/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
distribute-rgt-outN/A
distribute-lft-outN/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
Simplified4.0%
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
cos-lowering-cos.f644.0%
Applied egg-rr4.0%
Taylor expanded in re around 0
+-lowering-+.f64N/A
distribute-rgt-inN/A
associate-*r*N/A
*-commutativeN/A
*-commutativeN/A
associate-*l*N/A
distribute-lft-outN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6445.1%
Simplified45.1%
if 2.5999999999999998e24 < im Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
associate-*r*N/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
--lowering--.f64N/A
exp-negN/A
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
exp-lowering-exp.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6477.8%
Simplified77.8%
Taylor expanded in im around 0
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/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
Simplified76.0%
Final simplification63.1%
im\_m = (fabs.f64 im)
im\_s = (copysign.f64 #s(literal 1 binary64) im)
(FPCore (im_s re im_m)
:precision binary64
(*
im_s
(if (<= im_m 410.0)
(*
im_m
(+
-1.0
(*
(* im_m im_m)
(+
-0.16666666666666666
(*
im_m
(*
im_m
(+
-0.008333333333333333
(* (* im_m im_m) -0.0001984126984126984))))))))
(if (<= im_m 2.8e+24)
(*
im_m
(*
(+
-1.0
(*
(* im_m im_m)
(+ -0.16666666666666666 (* (* im_m im_m) -0.008333333333333333))))
(+ 1.0 (* re (* re (+ -0.5 (* (* re re) 0.041666666666666664)))))))
(*
(+ 0.5 (* -0.25 (* re re)))
(*
im_m
(+
-2.0
(*
(* im_m im_m)
(+
-0.3333333333333333
(*
im_m
(*
im_m
(+
-0.016666666666666666
(* (* im_m im_m) -0.0003968253968253968)))))))))))))im\_m = fabs(im);
im\_s = copysign(1.0, im);
double code(double im_s, double re, double im_m) {
double tmp;
if (im_m <= 410.0) {
tmp = im_m * (-1.0 + ((im_m * im_m) * (-0.16666666666666666 + (im_m * (im_m * (-0.008333333333333333 + ((im_m * im_m) * -0.0001984126984126984)))))));
} else if (im_m <= 2.8e+24) {
tmp = im_m * ((-1.0 + ((im_m * im_m) * (-0.16666666666666666 + ((im_m * im_m) * -0.008333333333333333)))) * (1.0 + (re * (re * (-0.5 + ((re * re) * 0.041666666666666664))))));
} else {
tmp = (0.5 + (-0.25 * (re * re))) * (im_m * (-2.0 + ((im_m * im_m) * (-0.3333333333333333 + (im_m * (im_m * (-0.016666666666666666 + ((im_m * im_m) * -0.0003968253968253968))))))));
}
return im_s * tmp;
}
im\_m = abs(im)
im\_s = copysign(1.0d0, im)
real(8) function code(im_s, re, im_m)
real(8), intent (in) :: im_s
real(8), intent (in) :: re
real(8), intent (in) :: im_m
real(8) :: tmp
if (im_m <= 410.0d0) then
tmp = im_m * ((-1.0d0) + ((im_m * im_m) * ((-0.16666666666666666d0) + (im_m * (im_m * ((-0.008333333333333333d0) + ((im_m * im_m) * (-0.0001984126984126984d0))))))))
else if (im_m <= 2.8d+24) then
tmp = im_m * (((-1.0d0) + ((im_m * im_m) * ((-0.16666666666666666d0) + ((im_m * im_m) * (-0.008333333333333333d0))))) * (1.0d0 + (re * (re * ((-0.5d0) + ((re * re) * 0.041666666666666664d0))))))
else
tmp = (0.5d0 + ((-0.25d0) * (re * re))) * (im_m * ((-2.0d0) + ((im_m * im_m) * ((-0.3333333333333333d0) + (im_m * (im_m * ((-0.016666666666666666d0) + ((im_m * im_m) * (-0.0003968253968253968d0)))))))))
end if
code = im_s * tmp
end function
im\_m = Math.abs(im);
im\_s = Math.copySign(1.0, im);
public static double code(double im_s, double re, double im_m) {
double tmp;
if (im_m <= 410.0) {
tmp = im_m * (-1.0 + ((im_m * im_m) * (-0.16666666666666666 + (im_m * (im_m * (-0.008333333333333333 + ((im_m * im_m) * -0.0001984126984126984)))))));
} else if (im_m <= 2.8e+24) {
tmp = im_m * ((-1.0 + ((im_m * im_m) * (-0.16666666666666666 + ((im_m * im_m) * -0.008333333333333333)))) * (1.0 + (re * (re * (-0.5 + ((re * re) * 0.041666666666666664))))));
} else {
tmp = (0.5 + (-0.25 * (re * re))) * (im_m * (-2.0 + ((im_m * im_m) * (-0.3333333333333333 + (im_m * (im_m * (-0.016666666666666666 + ((im_m * im_m) * -0.0003968253968253968))))))));
}
return im_s * tmp;
}
im\_m = math.fabs(im) im\_s = math.copysign(1.0, im) def code(im_s, re, im_m): tmp = 0 if im_m <= 410.0: tmp = im_m * (-1.0 + ((im_m * im_m) * (-0.16666666666666666 + (im_m * (im_m * (-0.008333333333333333 + ((im_m * im_m) * -0.0001984126984126984))))))) elif im_m <= 2.8e+24: tmp = im_m * ((-1.0 + ((im_m * im_m) * (-0.16666666666666666 + ((im_m * im_m) * -0.008333333333333333)))) * (1.0 + (re * (re * (-0.5 + ((re * re) * 0.041666666666666664)))))) else: tmp = (0.5 + (-0.25 * (re * re))) * (im_m * (-2.0 + ((im_m * im_m) * (-0.3333333333333333 + (im_m * (im_m * (-0.016666666666666666 + ((im_m * im_m) * -0.0003968253968253968)))))))) return im_s * tmp
im\_m = abs(im) im\_s = copysign(1.0, im) function code(im_s, re, im_m) tmp = 0.0 if (im_m <= 410.0) tmp = Float64(im_m * Float64(-1.0 + Float64(Float64(im_m * im_m) * Float64(-0.16666666666666666 + Float64(im_m * Float64(im_m * Float64(-0.008333333333333333 + Float64(Float64(im_m * im_m) * -0.0001984126984126984)))))))); elseif (im_m <= 2.8e+24) tmp = Float64(im_m * Float64(Float64(-1.0 + Float64(Float64(im_m * im_m) * Float64(-0.16666666666666666 + Float64(Float64(im_m * im_m) * -0.008333333333333333)))) * Float64(1.0 + Float64(re * Float64(re * Float64(-0.5 + Float64(Float64(re * re) * 0.041666666666666664))))))); else tmp = Float64(Float64(0.5 + Float64(-0.25 * Float64(re * re))) * Float64(im_m * Float64(-2.0 + Float64(Float64(im_m * im_m) * Float64(-0.3333333333333333 + Float64(im_m * Float64(im_m * Float64(-0.016666666666666666 + Float64(Float64(im_m * im_m) * -0.0003968253968253968))))))))); end return Float64(im_s * tmp) end
im\_m = abs(im); im\_s = sign(im) * abs(1.0); function tmp_2 = code(im_s, re, im_m) tmp = 0.0; if (im_m <= 410.0) tmp = im_m * (-1.0 + ((im_m * im_m) * (-0.16666666666666666 + (im_m * (im_m * (-0.008333333333333333 + ((im_m * im_m) * -0.0001984126984126984))))))); elseif (im_m <= 2.8e+24) tmp = im_m * ((-1.0 + ((im_m * im_m) * (-0.16666666666666666 + ((im_m * im_m) * -0.008333333333333333)))) * (1.0 + (re * (re * (-0.5 + ((re * re) * 0.041666666666666664)))))); else tmp = (0.5 + (-0.25 * (re * re))) * (im_m * (-2.0 + ((im_m * im_m) * (-0.3333333333333333 + (im_m * (im_m * (-0.016666666666666666 + ((im_m * im_m) * -0.0003968253968253968)))))))); end tmp_2 = im_s * tmp; end
im\_m = N[Abs[im], $MachinePrecision]
im\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[im]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[im$95$s_, re_, im$95$m_] := N[(im$95$s * If[LessEqual[im$95$m, 410.0], N[(im$95$m * N[(-1.0 + N[(N[(im$95$m * im$95$m), $MachinePrecision] * N[(-0.16666666666666666 + N[(im$95$m * N[(im$95$m * N[(-0.008333333333333333 + N[(N[(im$95$m * im$95$m), $MachinePrecision] * -0.0001984126984126984), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[im$95$m, 2.8e+24], N[(im$95$m * N[(N[(-1.0 + N[(N[(im$95$m * im$95$m), $MachinePrecision] * N[(-0.16666666666666666 + N[(N[(im$95$m * im$95$m), $MachinePrecision] * -0.008333333333333333), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(1.0 + N[(re * N[(re * N[(-0.5 + N[(N[(re * re), $MachinePrecision] * 0.041666666666666664), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(0.5 + N[(-0.25 * N[(re * re), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(im$95$m * N[(-2.0 + N[(N[(im$95$m * im$95$m), $MachinePrecision] * N[(-0.3333333333333333 + N[(im$95$m * N[(im$95$m * N[(-0.016666666666666666 + N[(N[(im$95$m * im$95$m), $MachinePrecision] * -0.0003968253968253968), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]
\begin{array}{l}
im\_m = \left|im\right|
\\
im\_s = \mathsf{copysign}\left(1, im\right)
\\
im\_s \cdot \begin{array}{l}
\mathbf{if}\;im\_m \leq 410:\\
\;\;\;\;im\_m \cdot \left(-1 + \left(im\_m \cdot im\_m\right) \cdot \left(-0.16666666666666666 + im\_m \cdot \left(im\_m \cdot \left(-0.008333333333333333 + \left(im\_m \cdot im\_m\right) \cdot -0.0001984126984126984\right)\right)\right)\right)\\
\mathbf{elif}\;im\_m \leq 2.8 \cdot 10^{+24}:\\
\;\;\;\;im\_m \cdot \left(\left(-1 + \left(im\_m \cdot im\_m\right) \cdot \left(-0.16666666666666666 + \left(im\_m \cdot im\_m\right) \cdot -0.008333333333333333\right)\right) \cdot \left(1 + re \cdot \left(re \cdot \left(-0.5 + \left(re \cdot re\right) \cdot 0.041666666666666664\right)\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\left(0.5 + -0.25 \cdot \left(re \cdot re\right)\right) \cdot \left(im\_m \cdot \left(-2 + \left(im\_m \cdot im\_m\right) \cdot \left(-0.3333333333333333 + im\_m \cdot \left(im\_m \cdot \left(-0.016666666666666666 + \left(im\_m \cdot im\_m\right) \cdot -0.0003968253968253968\right)\right)\right)\right)\right)\\
\end{array}
\end{array}
if im < 410Initial program 46.7%
Taylor expanded in im around 0
*-lowering-*.f64N/A
distribute-rgt-inN/A
associate-+r+N/A
*-commutativeN/A
+-commutativeN/A
+-commutativeN/A
Simplified92.4%
Taylor expanded in re around 0
Simplified60.2%
Taylor expanded in im around 0
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/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
Simplified60.2%
if 410 < im < 2.8000000000000002e24Initial program 100.0%
Taylor expanded in im around 0
*-lowering-*.f64N/A
*-commutativeN/A
+-commutativeN/A
distribute-lft-inN/A
associate-*r*N/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
distribute-rgt-outN/A
distribute-lft-outN/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
Simplified4.0%
Taylor expanded in re around 0
+-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
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6445.1%
Simplified45.1%
if 2.8000000000000002e24 < im Initial program 100.0%
Taylor expanded in re around 0
+-commutativeN/A
associate-*r*N/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
--lowering--.f64N/A
exp-negN/A
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
exp-lowering-exp.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6477.8%
Simplified77.8%
Taylor expanded in im around 0
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/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
Simplified76.0%
Final simplification63.1%
im\_m = (fabs.f64 im)
im\_s = (copysign.f64 #s(literal 1 binary64) im)
(FPCore (im_s re im_m)
:precision binary64
(let* ((t_0
(*
im_m
(+
-1.0
(*
(* im_m im_m)
(+
-0.16666666666666666
(*
im_m
(*
im_m
(+
-0.008333333333333333
(* (* im_m im_m) -0.0001984126984126984))))))))))
(*
im_s
(if (<= re 4.4e+101)
t_0
(if (<= re 2.1e+198)
(*
im_m
(*
(+ 0.5 (* -0.25 (* re re)))
(+ -2.0 (* (* im_m im_m) -0.3333333333333333))))
t_0)))))im\_m = fabs(im);
im\_s = copysign(1.0, im);
double code(double im_s, double re, double im_m) {
double t_0 = im_m * (-1.0 + ((im_m * im_m) * (-0.16666666666666666 + (im_m * (im_m * (-0.008333333333333333 + ((im_m * im_m) * -0.0001984126984126984)))))));
double tmp;
if (re <= 4.4e+101) {
tmp = t_0;
} else if (re <= 2.1e+198) {
tmp = im_m * ((0.5 + (-0.25 * (re * re))) * (-2.0 + ((im_m * im_m) * -0.3333333333333333)));
} else {
tmp = t_0;
}
return im_s * tmp;
}
im\_m = abs(im)
im\_s = copysign(1.0d0, im)
real(8) function code(im_s, re, im_m)
real(8), intent (in) :: im_s
real(8), intent (in) :: re
real(8), intent (in) :: im_m
real(8) :: t_0
real(8) :: tmp
t_0 = im_m * ((-1.0d0) + ((im_m * im_m) * ((-0.16666666666666666d0) + (im_m * (im_m * ((-0.008333333333333333d0) + ((im_m * im_m) * (-0.0001984126984126984d0))))))))
if (re <= 4.4d+101) then
tmp = t_0
else if (re <= 2.1d+198) then
tmp = im_m * ((0.5d0 + ((-0.25d0) * (re * re))) * ((-2.0d0) + ((im_m * im_m) * (-0.3333333333333333d0))))
else
tmp = t_0
end if
code = im_s * tmp
end function
im\_m = Math.abs(im);
im\_s = Math.copySign(1.0, im);
public static double code(double im_s, double re, double im_m) {
double t_0 = im_m * (-1.0 + ((im_m * im_m) * (-0.16666666666666666 + (im_m * (im_m * (-0.008333333333333333 + ((im_m * im_m) * -0.0001984126984126984)))))));
double tmp;
if (re <= 4.4e+101) {
tmp = t_0;
} else if (re <= 2.1e+198) {
tmp = im_m * ((0.5 + (-0.25 * (re * re))) * (-2.0 + ((im_m * im_m) * -0.3333333333333333)));
} else {
tmp = t_0;
}
return im_s * tmp;
}
im\_m = math.fabs(im) im\_s = math.copysign(1.0, im) def code(im_s, re, im_m): t_0 = im_m * (-1.0 + ((im_m * im_m) * (-0.16666666666666666 + (im_m * (im_m * (-0.008333333333333333 + ((im_m * im_m) * -0.0001984126984126984))))))) tmp = 0 if re <= 4.4e+101: tmp = t_0 elif re <= 2.1e+198: tmp = im_m * ((0.5 + (-0.25 * (re * re))) * (-2.0 + ((im_m * im_m) * -0.3333333333333333))) else: tmp = t_0 return im_s * tmp
im\_m = abs(im) im\_s = copysign(1.0, im) function code(im_s, re, im_m) t_0 = Float64(im_m * Float64(-1.0 + Float64(Float64(im_m * im_m) * Float64(-0.16666666666666666 + Float64(im_m * Float64(im_m * Float64(-0.008333333333333333 + Float64(Float64(im_m * im_m) * -0.0001984126984126984)))))))) tmp = 0.0 if (re <= 4.4e+101) tmp = t_0; elseif (re <= 2.1e+198) tmp = Float64(im_m * Float64(Float64(0.5 + Float64(-0.25 * Float64(re * re))) * Float64(-2.0 + Float64(Float64(im_m * im_m) * -0.3333333333333333)))); else tmp = t_0; end return Float64(im_s * tmp) end
im\_m = abs(im); im\_s = sign(im) * abs(1.0); function tmp_2 = code(im_s, re, im_m) t_0 = im_m * (-1.0 + ((im_m * im_m) * (-0.16666666666666666 + (im_m * (im_m * (-0.008333333333333333 + ((im_m * im_m) * -0.0001984126984126984))))))); tmp = 0.0; if (re <= 4.4e+101) tmp = t_0; elseif (re <= 2.1e+198) tmp = im_m * ((0.5 + (-0.25 * (re * re))) * (-2.0 + ((im_m * im_m) * -0.3333333333333333))); else tmp = t_0; end tmp_2 = im_s * tmp; end
im\_m = N[Abs[im], $MachinePrecision]
im\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[im]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[im$95$s_, re_, im$95$m_] := Block[{t$95$0 = N[(im$95$m * N[(-1.0 + N[(N[(im$95$m * im$95$m), $MachinePrecision] * N[(-0.16666666666666666 + N[(im$95$m * N[(im$95$m * N[(-0.008333333333333333 + N[(N[(im$95$m * im$95$m), $MachinePrecision] * -0.0001984126984126984), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, N[(im$95$s * If[LessEqual[re, 4.4e+101], t$95$0, If[LessEqual[re, 2.1e+198], N[(im$95$m * N[(N[(0.5 + N[(-0.25 * N[(re * re), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(-2.0 + N[(N[(im$95$m * im$95$m), $MachinePrecision] * -0.3333333333333333), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$0]]), $MachinePrecision]]
\begin{array}{l}
im\_m = \left|im\right|
\\
im\_s = \mathsf{copysign}\left(1, im\right)
\\
\begin{array}{l}
t_0 := im\_m \cdot \left(-1 + \left(im\_m \cdot im\_m\right) \cdot \left(-0.16666666666666666 + im\_m \cdot \left(im\_m \cdot \left(-0.008333333333333333 + \left(im\_m \cdot im\_m\right) \cdot -0.0001984126984126984\right)\right)\right)\right)\\
im\_s \cdot \begin{array}{l}
\mathbf{if}\;re \leq 4.4 \cdot 10^{+101}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;re \leq 2.1 \cdot 10^{+198}:\\
\;\;\;\;im\_m \cdot \left(\left(0.5 + -0.25 \cdot \left(re \cdot re\right)\right) \cdot \left(-2 + \left(im\_m \cdot im\_m\right) \cdot -0.3333333333333333\right)\right)\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
\end{array}
if re < 4.4000000000000001e101 or 2.10000000000000013e198 < re Initial program 57.5%
Taylor expanded in im around 0
*-lowering-*.f64N/A
distribute-rgt-inN/A
associate-+r+N/A
*-commutativeN/A
+-commutativeN/A
+-commutativeN/A
Simplified90.2%
Taylor expanded in re around 0
Simplified63.1%
Taylor expanded in im around 0
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/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
Simplified63.1%
if 4.4000000000000001e101 < re < 2.10000000000000013e198Initial program 83.9%
Taylor expanded in re around 0
+-commutativeN/A
associate-*r*N/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
--lowering--.f64N/A
exp-negN/A
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
exp-lowering-exp.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6439.4%
Simplified39.4%
Taylor expanded in im around 0
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-commutativeN/A
associate-*r*N/A
*-commutativeN/A
distribute-rgt-outN/A
metadata-evalN/A
sub-negN/A
*-lowering-*.f64N/A
Simplified39.3%
Final simplification61.4%
im\_m = (fabs.f64 im)
im\_s = (copysign.f64 #s(literal 1 binary64) im)
(FPCore (im_s re im_m)
:precision binary64
(let* ((t_0
(*
im_m
(+
-1.0
(*
im_m
(*
im_m
(+
-0.16666666666666666
(* (* im_m im_m) -0.008333333333333333))))))))
(*
im_s
(if (<= re 4.4e+101)
t_0
(if (<= re 2.1e+198)
(*
im_m
(*
(+ 0.5 (* -0.25 (* re re)))
(+ -2.0 (* (* im_m im_m) -0.3333333333333333))))
t_0)))))im\_m = fabs(im);
im\_s = copysign(1.0, im);
double code(double im_s, double re, double im_m) {
double t_0 = im_m * (-1.0 + (im_m * (im_m * (-0.16666666666666666 + ((im_m * im_m) * -0.008333333333333333)))));
double tmp;
if (re <= 4.4e+101) {
tmp = t_0;
} else if (re <= 2.1e+198) {
tmp = im_m * ((0.5 + (-0.25 * (re * re))) * (-2.0 + ((im_m * im_m) * -0.3333333333333333)));
} else {
tmp = t_0;
}
return im_s * tmp;
}
im\_m = abs(im)
im\_s = copysign(1.0d0, im)
real(8) function code(im_s, re, im_m)
real(8), intent (in) :: im_s
real(8), intent (in) :: re
real(8), intent (in) :: im_m
real(8) :: t_0
real(8) :: tmp
t_0 = im_m * ((-1.0d0) + (im_m * (im_m * ((-0.16666666666666666d0) + ((im_m * im_m) * (-0.008333333333333333d0))))))
if (re <= 4.4d+101) then
tmp = t_0
else if (re <= 2.1d+198) then
tmp = im_m * ((0.5d0 + ((-0.25d0) * (re * re))) * ((-2.0d0) + ((im_m * im_m) * (-0.3333333333333333d0))))
else
tmp = t_0
end if
code = im_s * tmp
end function
im\_m = Math.abs(im);
im\_s = Math.copySign(1.0, im);
public static double code(double im_s, double re, double im_m) {
double t_0 = im_m * (-1.0 + (im_m * (im_m * (-0.16666666666666666 + ((im_m * im_m) * -0.008333333333333333)))));
double tmp;
if (re <= 4.4e+101) {
tmp = t_0;
} else if (re <= 2.1e+198) {
tmp = im_m * ((0.5 + (-0.25 * (re * re))) * (-2.0 + ((im_m * im_m) * -0.3333333333333333)));
} else {
tmp = t_0;
}
return im_s * tmp;
}
im\_m = math.fabs(im) im\_s = math.copysign(1.0, im) def code(im_s, re, im_m): t_0 = im_m * (-1.0 + (im_m * (im_m * (-0.16666666666666666 + ((im_m * im_m) * -0.008333333333333333))))) tmp = 0 if re <= 4.4e+101: tmp = t_0 elif re <= 2.1e+198: tmp = im_m * ((0.5 + (-0.25 * (re * re))) * (-2.0 + ((im_m * im_m) * -0.3333333333333333))) else: tmp = t_0 return im_s * tmp
im\_m = abs(im) im\_s = copysign(1.0, im) function code(im_s, re, im_m) t_0 = Float64(im_m * Float64(-1.0 + Float64(im_m * Float64(im_m * Float64(-0.16666666666666666 + Float64(Float64(im_m * im_m) * -0.008333333333333333)))))) tmp = 0.0 if (re <= 4.4e+101) tmp = t_0; elseif (re <= 2.1e+198) tmp = Float64(im_m * Float64(Float64(0.5 + Float64(-0.25 * Float64(re * re))) * Float64(-2.0 + Float64(Float64(im_m * im_m) * -0.3333333333333333)))); else tmp = t_0; end return Float64(im_s * tmp) end
im\_m = abs(im); im\_s = sign(im) * abs(1.0); function tmp_2 = code(im_s, re, im_m) t_0 = im_m * (-1.0 + (im_m * (im_m * (-0.16666666666666666 + ((im_m * im_m) * -0.008333333333333333))))); tmp = 0.0; if (re <= 4.4e+101) tmp = t_0; elseif (re <= 2.1e+198) tmp = im_m * ((0.5 + (-0.25 * (re * re))) * (-2.0 + ((im_m * im_m) * -0.3333333333333333))); else tmp = t_0; end tmp_2 = im_s * tmp; end
im\_m = N[Abs[im], $MachinePrecision]
im\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[im]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[im$95$s_, re_, im$95$m_] := Block[{t$95$0 = N[(im$95$m * N[(-1.0 + N[(im$95$m * N[(im$95$m * N[(-0.16666666666666666 + N[(N[(im$95$m * im$95$m), $MachinePrecision] * -0.008333333333333333), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, N[(im$95$s * If[LessEqual[re, 4.4e+101], t$95$0, If[LessEqual[re, 2.1e+198], N[(im$95$m * N[(N[(0.5 + N[(-0.25 * N[(re * re), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(-2.0 + N[(N[(im$95$m * im$95$m), $MachinePrecision] * -0.3333333333333333), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$0]]), $MachinePrecision]]
\begin{array}{l}
im\_m = \left|im\right|
\\
im\_s = \mathsf{copysign}\left(1, im\right)
\\
\begin{array}{l}
t_0 := im\_m \cdot \left(-1 + im\_m \cdot \left(im\_m \cdot \left(-0.16666666666666666 + \left(im\_m \cdot im\_m\right) \cdot -0.008333333333333333\right)\right)\right)\\
im\_s \cdot \begin{array}{l}
\mathbf{if}\;re \leq 4.4 \cdot 10^{+101}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;re \leq 2.1 \cdot 10^{+198}:\\
\;\;\;\;im\_m \cdot \left(\left(0.5 + -0.25 \cdot \left(re \cdot re\right)\right) \cdot \left(-2 + \left(im\_m \cdot im\_m\right) \cdot -0.3333333333333333\right)\right)\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
\end{array}
if re < 4.4000000000000001e101 or 2.10000000000000013e198 < re Initial program 57.5%
Taylor expanded in im around 0
*-lowering-*.f64N/A
*-commutativeN/A
+-commutativeN/A
distribute-lft-inN/A
associate-*r*N/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
distribute-rgt-outN/A
distribute-lft-outN/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
Simplified86.3%
Taylor expanded in re around 0
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/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
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6459.5%
Simplified59.5%
if 4.4000000000000001e101 < re < 2.10000000000000013e198Initial program 83.9%
Taylor expanded in re around 0
+-commutativeN/A
associate-*r*N/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
--lowering--.f64N/A
exp-negN/A
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
exp-lowering-exp.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6439.4%
Simplified39.4%
Taylor expanded in im around 0
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-commutativeN/A
associate-*r*N/A
*-commutativeN/A
distribute-rgt-outN/A
metadata-evalN/A
sub-negN/A
*-lowering-*.f64N/A
Simplified39.3%
Final simplification58.0%
im\_m = (fabs.f64 im)
im\_s = (copysign.f64 #s(literal 1 binary64) im)
(FPCore (im_s re im_m)
:precision binary64
(let* ((t_0
(*
im_m
(+
-1.0
(*
im_m
(*
im_m
(+
-0.16666666666666666
(* (* im_m im_m) -0.008333333333333333))))))))
(*
im_s
(if (<= re 4.4e+101)
t_0
(if (<= re 2.1e+198)
(* (+ 0.5 (* -0.25 (* re re))) (* im_m -2.0))
t_0)))))im\_m = fabs(im);
im\_s = copysign(1.0, im);
double code(double im_s, double re, double im_m) {
double t_0 = im_m * (-1.0 + (im_m * (im_m * (-0.16666666666666666 + ((im_m * im_m) * -0.008333333333333333)))));
double tmp;
if (re <= 4.4e+101) {
tmp = t_0;
} else if (re <= 2.1e+198) {
tmp = (0.5 + (-0.25 * (re * re))) * (im_m * -2.0);
} else {
tmp = t_0;
}
return im_s * tmp;
}
im\_m = abs(im)
im\_s = copysign(1.0d0, im)
real(8) function code(im_s, re, im_m)
real(8), intent (in) :: im_s
real(8), intent (in) :: re
real(8), intent (in) :: im_m
real(8) :: t_0
real(8) :: tmp
t_0 = im_m * ((-1.0d0) + (im_m * (im_m * ((-0.16666666666666666d0) + ((im_m * im_m) * (-0.008333333333333333d0))))))
if (re <= 4.4d+101) then
tmp = t_0
else if (re <= 2.1d+198) then
tmp = (0.5d0 + ((-0.25d0) * (re * re))) * (im_m * (-2.0d0))
else
tmp = t_0
end if
code = im_s * tmp
end function
im\_m = Math.abs(im);
im\_s = Math.copySign(1.0, im);
public static double code(double im_s, double re, double im_m) {
double t_0 = im_m * (-1.0 + (im_m * (im_m * (-0.16666666666666666 + ((im_m * im_m) * -0.008333333333333333)))));
double tmp;
if (re <= 4.4e+101) {
tmp = t_0;
} else if (re <= 2.1e+198) {
tmp = (0.5 + (-0.25 * (re * re))) * (im_m * -2.0);
} else {
tmp = t_0;
}
return im_s * tmp;
}
im\_m = math.fabs(im) im\_s = math.copysign(1.0, im) def code(im_s, re, im_m): t_0 = im_m * (-1.0 + (im_m * (im_m * (-0.16666666666666666 + ((im_m * im_m) * -0.008333333333333333))))) tmp = 0 if re <= 4.4e+101: tmp = t_0 elif re <= 2.1e+198: tmp = (0.5 + (-0.25 * (re * re))) * (im_m * -2.0) else: tmp = t_0 return im_s * tmp
im\_m = abs(im) im\_s = copysign(1.0, im) function code(im_s, re, im_m) t_0 = Float64(im_m * Float64(-1.0 + Float64(im_m * Float64(im_m * Float64(-0.16666666666666666 + Float64(Float64(im_m * im_m) * -0.008333333333333333)))))) tmp = 0.0 if (re <= 4.4e+101) tmp = t_0; elseif (re <= 2.1e+198) tmp = Float64(Float64(0.5 + Float64(-0.25 * Float64(re * re))) * Float64(im_m * -2.0)); else tmp = t_0; end return Float64(im_s * tmp) end
im\_m = abs(im); im\_s = sign(im) * abs(1.0); function tmp_2 = code(im_s, re, im_m) t_0 = im_m * (-1.0 + (im_m * (im_m * (-0.16666666666666666 + ((im_m * im_m) * -0.008333333333333333))))); tmp = 0.0; if (re <= 4.4e+101) tmp = t_0; elseif (re <= 2.1e+198) tmp = (0.5 + (-0.25 * (re * re))) * (im_m * -2.0); else tmp = t_0; end tmp_2 = im_s * tmp; end
im\_m = N[Abs[im], $MachinePrecision]
im\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[im]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[im$95$s_, re_, im$95$m_] := Block[{t$95$0 = N[(im$95$m * N[(-1.0 + N[(im$95$m * N[(im$95$m * N[(-0.16666666666666666 + N[(N[(im$95$m * im$95$m), $MachinePrecision] * -0.008333333333333333), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, N[(im$95$s * If[LessEqual[re, 4.4e+101], t$95$0, If[LessEqual[re, 2.1e+198], N[(N[(0.5 + N[(-0.25 * N[(re * re), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(im$95$m * -2.0), $MachinePrecision]), $MachinePrecision], t$95$0]]), $MachinePrecision]]
\begin{array}{l}
im\_m = \left|im\right|
\\
im\_s = \mathsf{copysign}\left(1, im\right)
\\
\begin{array}{l}
t_0 := im\_m \cdot \left(-1 + im\_m \cdot \left(im\_m \cdot \left(-0.16666666666666666 + \left(im\_m \cdot im\_m\right) \cdot -0.008333333333333333\right)\right)\right)\\
im\_s \cdot \begin{array}{l}
\mathbf{if}\;re \leq 4.4 \cdot 10^{+101}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;re \leq 2.1 \cdot 10^{+198}:\\
\;\;\;\;\left(0.5 + -0.25 \cdot \left(re \cdot re\right)\right) \cdot \left(im\_m \cdot -2\right)\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
\end{array}
if re < 4.4000000000000001e101 or 2.10000000000000013e198 < re Initial program 57.5%
Taylor expanded in im around 0
*-lowering-*.f64N/A
*-commutativeN/A
+-commutativeN/A
distribute-lft-inN/A
associate-*r*N/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
distribute-rgt-outN/A
distribute-lft-outN/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
Simplified86.3%
Taylor expanded in re around 0
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/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
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6459.5%
Simplified59.5%
if 4.4000000000000001e101 < re < 2.10000000000000013e198Initial program 83.9%
Taylor expanded in re around 0
+-commutativeN/A
associate-*r*N/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
--lowering--.f64N/A
exp-negN/A
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
exp-lowering-exp.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6439.4%
Simplified39.4%
Taylor expanded in im around 0
*-commutativeN/A
*-lowering-*.f6439.3%
Simplified39.3%
Final simplification58.0%
im\_m = (fabs.f64 im)
im\_s = (copysign.f64 #s(literal 1 binary64) im)
(FPCore (im_s re im_m)
:precision binary64
(let* ((t_0 (* im_m (+ -1.0 (* (* im_m im_m) -0.16666666666666666)))))
(*
im_s
(if (<= re 4.4e+101)
t_0
(if (<= re 2.1e+198)
(* (+ 0.5 (* -0.25 (* re re))) (* im_m -2.0))
t_0)))))im\_m = fabs(im);
im\_s = copysign(1.0, im);
double code(double im_s, double re, double im_m) {
double t_0 = im_m * (-1.0 + ((im_m * im_m) * -0.16666666666666666));
double tmp;
if (re <= 4.4e+101) {
tmp = t_0;
} else if (re <= 2.1e+198) {
tmp = (0.5 + (-0.25 * (re * re))) * (im_m * -2.0);
} else {
tmp = t_0;
}
return im_s * tmp;
}
im\_m = abs(im)
im\_s = copysign(1.0d0, im)
real(8) function code(im_s, re, im_m)
real(8), intent (in) :: im_s
real(8), intent (in) :: re
real(8), intent (in) :: im_m
real(8) :: t_0
real(8) :: tmp
t_0 = im_m * ((-1.0d0) + ((im_m * im_m) * (-0.16666666666666666d0)))
if (re <= 4.4d+101) then
tmp = t_0
else if (re <= 2.1d+198) then
tmp = (0.5d0 + ((-0.25d0) * (re * re))) * (im_m * (-2.0d0))
else
tmp = t_0
end if
code = im_s * tmp
end function
im\_m = Math.abs(im);
im\_s = Math.copySign(1.0, im);
public static double code(double im_s, double re, double im_m) {
double t_0 = im_m * (-1.0 + ((im_m * im_m) * -0.16666666666666666));
double tmp;
if (re <= 4.4e+101) {
tmp = t_0;
} else if (re <= 2.1e+198) {
tmp = (0.5 + (-0.25 * (re * re))) * (im_m * -2.0);
} else {
tmp = t_0;
}
return im_s * tmp;
}
im\_m = math.fabs(im) im\_s = math.copysign(1.0, im) def code(im_s, re, im_m): t_0 = im_m * (-1.0 + ((im_m * im_m) * -0.16666666666666666)) tmp = 0 if re <= 4.4e+101: tmp = t_0 elif re <= 2.1e+198: tmp = (0.5 + (-0.25 * (re * re))) * (im_m * -2.0) else: tmp = t_0 return im_s * tmp
im\_m = abs(im) im\_s = copysign(1.0, im) function code(im_s, re, im_m) t_0 = Float64(im_m * Float64(-1.0 + Float64(Float64(im_m * im_m) * -0.16666666666666666))) tmp = 0.0 if (re <= 4.4e+101) tmp = t_0; elseif (re <= 2.1e+198) tmp = Float64(Float64(0.5 + Float64(-0.25 * Float64(re * re))) * Float64(im_m * -2.0)); else tmp = t_0; end return Float64(im_s * tmp) end
im\_m = abs(im); im\_s = sign(im) * abs(1.0); function tmp_2 = code(im_s, re, im_m) t_0 = im_m * (-1.0 + ((im_m * im_m) * -0.16666666666666666)); tmp = 0.0; if (re <= 4.4e+101) tmp = t_0; elseif (re <= 2.1e+198) tmp = (0.5 + (-0.25 * (re * re))) * (im_m * -2.0); else tmp = t_0; end tmp_2 = im_s * tmp; end
im\_m = N[Abs[im], $MachinePrecision]
im\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[im]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[im$95$s_, re_, im$95$m_] := Block[{t$95$0 = N[(im$95$m * N[(-1.0 + N[(N[(im$95$m * im$95$m), $MachinePrecision] * -0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, N[(im$95$s * If[LessEqual[re, 4.4e+101], t$95$0, If[LessEqual[re, 2.1e+198], N[(N[(0.5 + N[(-0.25 * N[(re * re), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(im$95$m * -2.0), $MachinePrecision]), $MachinePrecision], t$95$0]]), $MachinePrecision]]
\begin{array}{l}
im\_m = \left|im\right|
\\
im\_s = \mathsf{copysign}\left(1, im\right)
\\
\begin{array}{l}
t_0 := im\_m \cdot \left(-1 + \left(im\_m \cdot im\_m\right) \cdot -0.16666666666666666\right)\\
im\_s \cdot \begin{array}{l}
\mathbf{if}\;re \leq 4.4 \cdot 10^{+101}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;re \leq 2.1 \cdot 10^{+198}:\\
\;\;\;\;\left(0.5 + -0.25 \cdot \left(re \cdot re\right)\right) \cdot \left(im\_m \cdot -2\right)\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
\end{array}
if re < 4.4000000000000001e101 or 2.10000000000000013e198 < re Initial program 57.5%
Taylor expanded in im around 0
*-lowering-*.f64N/A
distribute-rgt-inN/A
associate-+r+N/A
*-commutativeN/A
+-commutativeN/A
+-commutativeN/A
Simplified90.2%
Taylor expanded in re around 0
Simplified63.1%
Taylor expanded in im around 0
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6454.9%
Simplified54.9%
if 4.4000000000000001e101 < re < 2.10000000000000013e198Initial program 83.9%
Taylor expanded in re around 0
+-commutativeN/A
associate-*r*N/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
--lowering--.f64N/A
exp-negN/A
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
exp-lowering-exp.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6439.4%
Simplified39.4%
Taylor expanded in im around 0
*-commutativeN/A
*-lowering-*.f6439.3%
Simplified39.3%
Final simplification53.8%
im\_m = (fabs.f64 im)
im\_s = (copysign.f64 #s(literal 1 binary64) im)
(FPCore (im_s re im_m)
:precision binary64
(let* ((t_0 (* im_m (+ -1.0 (* (* im_m im_m) -0.16666666666666666)))))
(*
im_s
(if (<= re 4.4e+101)
t_0
(if (<= re 2.1e+198) (* im_m (+ -1.0 (* 0.5 (* re re)))) t_0)))))im\_m = fabs(im);
im\_s = copysign(1.0, im);
double code(double im_s, double re, double im_m) {
double t_0 = im_m * (-1.0 + ((im_m * im_m) * -0.16666666666666666));
double tmp;
if (re <= 4.4e+101) {
tmp = t_0;
} else if (re <= 2.1e+198) {
tmp = im_m * (-1.0 + (0.5 * (re * re)));
} else {
tmp = t_0;
}
return im_s * tmp;
}
im\_m = abs(im)
im\_s = copysign(1.0d0, im)
real(8) function code(im_s, re, im_m)
real(8), intent (in) :: im_s
real(8), intent (in) :: re
real(8), intent (in) :: im_m
real(8) :: t_0
real(8) :: tmp
t_0 = im_m * ((-1.0d0) + ((im_m * im_m) * (-0.16666666666666666d0)))
if (re <= 4.4d+101) then
tmp = t_0
else if (re <= 2.1d+198) then
tmp = im_m * ((-1.0d0) + (0.5d0 * (re * re)))
else
tmp = t_0
end if
code = im_s * tmp
end function
im\_m = Math.abs(im);
im\_s = Math.copySign(1.0, im);
public static double code(double im_s, double re, double im_m) {
double t_0 = im_m * (-1.0 + ((im_m * im_m) * -0.16666666666666666));
double tmp;
if (re <= 4.4e+101) {
tmp = t_0;
} else if (re <= 2.1e+198) {
tmp = im_m * (-1.0 + (0.5 * (re * re)));
} else {
tmp = t_0;
}
return im_s * tmp;
}
im\_m = math.fabs(im) im\_s = math.copysign(1.0, im) def code(im_s, re, im_m): t_0 = im_m * (-1.0 + ((im_m * im_m) * -0.16666666666666666)) tmp = 0 if re <= 4.4e+101: tmp = t_0 elif re <= 2.1e+198: tmp = im_m * (-1.0 + (0.5 * (re * re))) else: tmp = t_0 return im_s * tmp
im\_m = abs(im) im\_s = copysign(1.0, im) function code(im_s, re, im_m) t_0 = Float64(im_m * Float64(-1.0 + Float64(Float64(im_m * im_m) * -0.16666666666666666))) tmp = 0.0 if (re <= 4.4e+101) tmp = t_0; elseif (re <= 2.1e+198) tmp = Float64(im_m * Float64(-1.0 + Float64(0.5 * Float64(re * re)))); else tmp = t_0; end return Float64(im_s * tmp) end
im\_m = abs(im); im\_s = sign(im) * abs(1.0); function tmp_2 = code(im_s, re, im_m) t_0 = im_m * (-1.0 + ((im_m * im_m) * -0.16666666666666666)); tmp = 0.0; if (re <= 4.4e+101) tmp = t_0; elseif (re <= 2.1e+198) tmp = im_m * (-1.0 + (0.5 * (re * re))); else tmp = t_0; end tmp_2 = im_s * tmp; end
im\_m = N[Abs[im], $MachinePrecision]
im\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[im]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[im$95$s_, re_, im$95$m_] := Block[{t$95$0 = N[(im$95$m * N[(-1.0 + N[(N[(im$95$m * im$95$m), $MachinePrecision] * -0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, N[(im$95$s * If[LessEqual[re, 4.4e+101], t$95$0, If[LessEqual[re, 2.1e+198], N[(im$95$m * N[(-1.0 + N[(0.5 * N[(re * re), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$0]]), $MachinePrecision]]
\begin{array}{l}
im\_m = \left|im\right|
\\
im\_s = \mathsf{copysign}\left(1, im\right)
\\
\begin{array}{l}
t_0 := im\_m \cdot \left(-1 + \left(im\_m \cdot im\_m\right) \cdot -0.16666666666666666\right)\\
im\_s \cdot \begin{array}{l}
\mathbf{if}\;re \leq 4.4 \cdot 10^{+101}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;re \leq 2.1 \cdot 10^{+198}:\\
\;\;\;\;im\_m \cdot \left(-1 + 0.5 \cdot \left(re \cdot re\right)\right)\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
\end{array}
if re < 4.4000000000000001e101 or 2.10000000000000013e198 < re Initial program 57.5%
Taylor expanded in im around 0
*-lowering-*.f64N/A
distribute-rgt-inN/A
associate-+r+N/A
*-commutativeN/A
+-commutativeN/A
+-commutativeN/A
Simplified90.2%
Taylor expanded in re around 0
Simplified63.1%
Taylor expanded in im around 0
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6454.9%
Simplified54.9%
if 4.4000000000000001e101 < re < 2.10000000000000013e198Initial program 83.9%
Taylor expanded in im around 0
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
cos-lowering-cos.f6422.7%
Simplified22.7%
Taylor expanded in re around 0
sub-negN/A
+-commutativeN/A
mul-1-negN/A
*-commutativeN/A
associate-*r*N/A
metadata-evalN/A
associate-*r*N/A
distribute-rgt-inN/A
metadata-evalN/A
distribute-lft-inN/A
*-lowering-*.f64N/A
distribute-lft-inN/A
metadata-evalN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-commutativeN/A
associate-*l*N/A
metadata-evalN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6439.3%
Simplified39.3%
Final simplification53.8%
im\_m = (fabs.f64 im) im\_s = (copysign.f64 #s(literal 1 binary64) im) (FPCore (im_s re im_m) :precision binary64 (* im_s (if (<= im_m 1.1e+171) (- 0.0 im_m) (- 0.0 (/ (* im_m im_m) im_m)))))
im\_m = fabs(im);
im\_s = copysign(1.0, im);
double code(double im_s, double re, double im_m) {
double tmp;
if (im_m <= 1.1e+171) {
tmp = 0.0 - im_m;
} else {
tmp = 0.0 - ((im_m * im_m) / im_m);
}
return im_s * tmp;
}
im\_m = abs(im)
im\_s = copysign(1.0d0, im)
real(8) function code(im_s, re, im_m)
real(8), intent (in) :: im_s
real(8), intent (in) :: re
real(8), intent (in) :: im_m
real(8) :: tmp
if (im_m <= 1.1d+171) then
tmp = 0.0d0 - im_m
else
tmp = 0.0d0 - ((im_m * im_m) / im_m)
end if
code = im_s * tmp
end function
im\_m = Math.abs(im);
im\_s = Math.copySign(1.0, im);
public static double code(double im_s, double re, double im_m) {
double tmp;
if (im_m <= 1.1e+171) {
tmp = 0.0 - im_m;
} else {
tmp = 0.0 - ((im_m * im_m) / im_m);
}
return im_s * tmp;
}
im\_m = math.fabs(im) im\_s = math.copysign(1.0, im) def code(im_s, re, im_m): tmp = 0 if im_m <= 1.1e+171: tmp = 0.0 - im_m else: tmp = 0.0 - ((im_m * im_m) / im_m) return im_s * tmp
im\_m = abs(im) im\_s = copysign(1.0, im) function code(im_s, re, im_m) tmp = 0.0 if (im_m <= 1.1e+171) tmp = Float64(0.0 - im_m); else tmp = Float64(0.0 - Float64(Float64(im_m * im_m) / im_m)); end return Float64(im_s * tmp) end
im\_m = abs(im); im\_s = sign(im) * abs(1.0); function tmp_2 = code(im_s, re, im_m) tmp = 0.0; if (im_m <= 1.1e+171) tmp = 0.0 - im_m; else tmp = 0.0 - ((im_m * im_m) / im_m); end tmp_2 = im_s * tmp; end
im\_m = N[Abs[im], $MachinePrecision]
im\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[im]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[im$95$s_, re_, im$95$m_] := N[(im$95$s * If[LessEqual[im$95$m, 1.1e+171], N[(0.0 - im$95$m), $MachinePrecision], N[(0.0 - N[(N[(im$95$m * im$95$m), $MachinePrecision] / im$95$m), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
im\_m = \left|im\right|
\\
im\_s = \mathsf{copysign}\left(1, im\right)
\\
im\_s \cdot \begin{array}{l}
\mathbf{if}\;im\_m \leq 1.1 \cdot 10^{+171}:\\
\;\;\;\;0 - im\_m\\
\mathbf{else}:\\
\;\;\;\;0 - \frac{im\_m \cdot im\_m}{im\_m}\\
\end{array}
\end{array}
if im < 1.1e171Initial program 56.1%
Taylor expanded in im around 0
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
cos-lowering-cos.f6450.7%
Simplified50.7%
Taylor expanded in re around 0
mul-1-negN/A
neg-sub0N/A
--lowering--.f6430.0%
Simplified30.0%
sub0-negN/A
neg-lowering-neg.f6430.0%
Applied egg-rr30.0%
if 1.1e171 < im Initial program 100.0%
Taylor expanded in im around 0
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
cos-lowering-cos.f647.5%
Simplified7.5%
Taylor expanded in re around 0
mul-1-negN/A
neg-sub0N/A
--lowering--.f646.3%
Simplified6.3%
flip--N/A
+-lft-identityN/A
/-lowering-/.f64N/A
metadata-evalN/A
--lowering--.f64N/A
*-lowering-*.f6484.2%
Applied egg-rr84.2%
Final simplification34.0%
im\_m = (fabs.f64 im) im\_s = (copysign.f64 #s(literal 1 binary64) im) (FPCore (im_s re im_m) :precision binary64 (* im_s (* im_m (+ -1.0 (* (* im_m im_m) -0.16666666666666666)))))
im\_m = fabs(im);
im\_s = copysign(1.0, im);
double code(double im_s, double re, double im_m) {
return im_s * (im_m * (-1.0 + ((im_m * im_m) * -0.16666666666666666)));
}
im\_m = abs(im)
im\_s = copysign(1.0d0, im)
real(8) function code(im_s, re, im_m)
real(8), intent (in) :: im_s
real(8), intent (in) :: re
real(8), intent (in) :: im_m
code = im_s * (im_m * ((-1.0d0) + ((im_m * im_m) * (-0.16666666666666666d0))))
end function
im\_m = Math.abs(im);
im\_s = Math.copySign(1.0, im);
public static double code(double im_s, double re, double im_m) {
return im_s * (im_m * (-1.0 + ((im_m * im_m) * -0.16666666666666666)));
}
im\_m = math.fabs(im) im\_s = math.copysign(1.0, im) def code(im_s, re, im_m): return im_s * (im_m * (-1.0 + ((im_m * im_m) * -0.16666666666666666)))
im\_m = abs(im) im\_s = copysign(1.0, im) function code(im_s, re, im_m) return Float64(im_s * Float64(im_m * Float64(-1.0 + Float64(Float64(im_m * im_m) * -0.16666666666666666)))) end
im\_m = abs(im); im\_s = sign(im) * abs(1.0); function tmp = code(im_s, re, im_m) tmp = im_s * (im_m * (-1.0 + ((im_m * im_m) * -0.16666666666666666))); end
im\_m = N[Abs[im], $MachinePrecision]
im\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[im]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[im$95$s_, re_, im$95$m_] := N[(im$95$s * N[(im$95$m * N[(-1.0 + N[(N[(im$95$m * im$95$m), $MachinePrecision] * -0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
im\_m = \left|im\right|
\\
im\_s = \mathsf{copysign}\left(1, im\right)
\\
im\_s \cdot \left(im\_m \cdot \left(-1 + \left(im\_m \cdot im\_m\right) \cdot -0.16666666666666666\right)\right)
\end{array}
Initial program 59.4%
Taylor expanded in im around 0
*-lowering-*.f64N/A
distribute-rgt-inN/A
associate-+r+N/A
*-commutativeN/A
+-commutativeN/A
+-commutativeN/A
Simplified89.8%
Taylor expanded in re around 0
Simplified60.5%
Taylor expanded in im around 0
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6452.5%
Simplified52.5%
Final simplification52.5%
im\_m = (fabs.f64 im) im\_s = (copysign.f64 #s(literal 1 binary64) im) (FPCore (im_s re im_m) :precision binary64 (* im_s (- 0.0 im_m)))
im\_m = fabs(im);
im\_s = copysign(1.0, im);
double code(double im_s, double re, double im_m) {
return im_s * (0.0 - im_m);
}
im\_m = abs(im)
im\_s = copysign(1.0d0, im)
real(8) function code(im_s, re, im_m)
real(8), intent (in) :: im_s
real(8), intent (in) :: re
real(8), intent (in) :: im_m
code = im_s * (0.0d0 - im_m)
end function
im\_m = Math.abs(im);
im\_s = Math.copySign(1.0, im);
public static double code(double im_s, double re, double im_m) {
return im_s * (0.0 - im_m);
}
im\_m = math.fabs(im) im\_s = math.copysign(1.0, im) def code(im_s, re, im_m): return im_s * (0.0 - im_m)
im\_m = abs(im) im\_s = copysign(1.0, im) function code(im_s, re, im_m) return Float64(im_s * Float64(0.0 - im_m)) end
im\_m = abs(im); im\_s = sign(im) * abs(1.0); function tmp = code(im_s, re, im_m) tmp = im_s * (0.0 - im_m); end
im\_m = N[Abs[im], $MachinePrecision]
im\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[im]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[im$95$s_, re_, im$95$m_] := N[(im$95$s * N[(0.0 - im$95$m), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
im\_m = \left|im\right|
\\
im\_s = \mathsf{copysign}\left(1, im\right)
\\
im\_s \cdot \left(0 - im\_m\right)
\end{array}
Initial program 59.4%
Taylor expanded in im around 0
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
cos-lowering-cos.f6447.5%
Simplified47.5%
Taylor expanded in re around 0
mul-1-negN/A
neg-sub0N/A
--lowering--.f6428.2%
Simplified28.2%
sub0-negN/A
neg-lowering-neg.f6428.2%
Applied egg-rr28.2%
Final simplification28.2%
(FPCore (re im)
:precision binary64
(if (< (fabs im) 1.0)
(-
(*
(cos re)
(+
(+ im (* (* (* 0.16666666666666666 im) im) im))
(* (* (* (* (* 0.008333333333333333 im) im) im) im) im))))
(* (* 0.5 (cos re)) (- (exp (- 0.0 im)) (exp im)))))
double code(double re, double im) {
double tmp;
if (fabs(im) < 1.0) {
tmp = -(cos(re) * ((im + (((0.16666666666666666 * im) * im) * im)) + (((((0.008333333333333333 * im) * im) * im) * im) * im)));
} else {
tmp = (0.5 * cos(re)) * (exp((0.0 - im)) - exp(im));
}
return tmp;
}
real(8) function code(re, im)
real(8), intent (in) :: re
real(8), intent (in) :: im
real(8) :: tmp
if (abs(im) < 1.0d0) then
tmp = -(cos(re) * ((im + (((0.16666666666666666d0 * im) * im) * im)) + (((((0.008333333333333333d0 * im) * im) * im) * im) * im)))
else
tmp = (0.5d0 * cos(re)) * (exp((0.0d0 - im)) - exp(im))
end if
code = tmp
end function
public static double code(double re, double im) {
double tmp;
if (Math.abs(im) < 1.0) {
tmp = -(Math.cos(re) * ((im + (((0.16666666666666666 * im) * im) * im)) + (((((0.008333333333333333 * im) * im) * im) * im) * im)));
} else {
tmp = (0.5 * Math.cos(re)) * (Math.exp((0.0 - im)) - Math.exp(im));
}
return tmp;
}
def code(re, im): tmp = 0 if math.fabs(im) < 1.0: tmp = -(math.cos(re) * ((im + (((0.16666666666666666 * im) * im) * im)) + (((((0.008333333333333333 * im) * im) * im) * im) * im))) else: tmp = (0.5 * math.cos(re)) * (math.exp((0.0 - im)) - math.exp(im)) return tmp
function code(re, im) tmp = 0.0 if (abs(im) < 1.0) tmp = Float64(-Float64(cos(re) * Float64(Float64(im + Float64(Float64(Float64(0.16666666666666666 * im) * im) * im)) + Float64(Float64(Float64(Float64(Float64(0.008333333333333333 * im) * im) * im) * im) * im)))); else tmp = Float64(Float64(0.5 * cos(re)) * Float64(exp(Float64(0.0 - im)) - exp(im))); end return tmp end
function tmp_2 = code(re, im) tmp = 0.0; if (abs(im) < 1.0) tmp = -(cos(re) * ((im + (((0.16666666666666666 * im) * im) * im)) + (((((0.008333333333333333 * im) * im) * im) * im) * im))); else tmp = (0.5 * cos(re)) * (exp((0.0 - im)) - exp(im)); end tmp_2 = tmp; end
code[re_, im_] := If[Less[N[Abs[im], $MachinePrecision], 1.0], (-N[(N[Cos[re], $MachinePrecision] * N[(N[(im + N[(N[(N[(0.16666666666666666 * im), $MachinePrecision] * im), $MachinePrecision] * im), $MachinePrecision]), $MachinePrecision] + N[(N[(N[(N[(N[(0.008333333333333333 * im), $MachinePrecision] * im), $MachinePrecision] * im), $MachinePrecision] * im), $MachinePrecision] * im), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), N[(N[(0.5 * N[Cos[re], $MachinePrecision]), $MachinePrecision] * N[(N[Exp[N[(0.0 - im), $MachinePrecision]], $MachinePrecision] - N[Exp[im], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\left|im\right| < 1:\\
\;\;\;\;-\cos re \cdot \left(\left(im + \left(\left(0.16666666666666666 \cdot im\right) \cdot im\right) \cdot im\right) + \left(\left(\left(\left(0.008333333333333333 \cdot im\right) \cdot im\right) \cdot im\right) \cdot im\right) \cdot im\right)\\
\mathbf{else}:\\
\;\;\;\;\left(0.5 \cdot \cos re\right) \cdot \left(e^{0 - im} - e^{im}\right)\\
\end{array}
\end{array}
herbie shell --seed 2024155
(FPCore (re im)
:name "math.sin on complex, imaginary part"
:precision binary64
:alt
(! :herbie-platform default (if (< (fabs im) 1) (- (* (cos re) (+ im (* 1/6 im im im) (* 1/120 im im im im im)))) (* (* 1/2 (cos re)) (- (exp (- 0 im)) (exp im)))))
(* (* 0.5 (cos re)) (- (exp (- 0.0 im)) (exp im))))