
(FPCore (re im) :precision binary64 (* (* 0.5 (sin re)) (+ (exp (- 0.0 im)) (exp im))))
double code(double re, double im) {
return (0.5 * sin(re)) * (exp((0.0 - im)) + exp(im));
}
module fmin_fmax_functions
implicit none
private
public fmax
public fmin
interface fmax
module procedure fmax88
module procedure fmax44
module procedure fmax84
module procedure fmax48
end interface
interface fmin
module procedure fmin88
module procedure fmin44
module procedure fmin84
module procedure fmin48
end interface
contains
real(8) function fmax88(x, y) result (res)
real(8), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(x, max(x, y), y /= y), x /= x)
end function
real(4) function fmax44(x, y) result (res)
real(4), intent (in) :: x
real(4), intent (in) :: y
res = merge(y, merge(x, max(x, y), y /= y), x /= x)
end function
real(8) function fmax84(x, y) result(res)
real(8), intent (in) :: x
real(4), intent (in) :: y
res = merge(dble(y), merge(x, max(x, dble(y)), y /= y), x /= x)
end function
real(8) function fmax48(x, y) result(res)
real(4), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(dble(x), max(dble(x), y), y /= y), x /= x)
end function
real(8) function fmin88(x, y) result (res)
real(8), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(x, min(x, y), y /= y), x /= x)
end function
real(4) function fmin44(x, y) result (res)
real(4), intent (in) :: x
real(4), intent (in) :: y
res = merge(y, merge(x, min(x, y), y /= y), x /= x)
end function
real(8) function fmin84(x, y) result(res)
real(8), intent (in) :: x
real(4), intent (in) :: y
res = merge(dble(y), merge(x, min(x, dble(y)), y /= y), x /= x)
end function
real(8) function fmin48(x, y) result(res)
real(4), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(dble(x), min(dble(x), y), y /= y), x /= x)
end function
end module
real(8) function code(re, im)
use fmin_fmax_functions
real(8), intent (in) :: re
real(8), intent (in) :: im
code = (0.5d0 * sin(re)) * (exp((0.0d0 - im)) + exp(im))
end function
public static double code(double re, double im) {
return (0.5 * Math.sin(re)) * (Math.exp((0.0 - im)) + Math.exp(im));
}
def code(re, im): return (0.5 * math.sin(re)) * (math.exp((0.0 - im)) + math.exp(im))
function code(re, im) return Float64(Float64(0.5 * sin(re)) * Float64(exp(Float64(0.0 - im)) + exp(im))) end
function tmp = code(re, im) tmp = (0.5 * sin(re)) * (exp((0.0 - im)) + exp(im)); end
code[re_, im_] := N[(N[(0.5 * N[Sin[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 \sin re\right) \cdot \left(e^{0 - im} + e^{im}\right)
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 19 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (re im) :precision binary64 (* (* 0.5 (sin re)) (+ (exp (- 0.0 im)) (exp im))))
double code(double re, double im) {
return (0.5 * sin(re)) * (exp((0.0 - im)) + exp(im));
}
module fmin_fmax_functions
implicit none
private
public fmax
public fmin
interface fmax
module procedure fmax88
module procedure fmax44
module procedure fmax84
module procedure fmax48
end interface
interface fmin
module procedure fmin88
module procedure fmin44
module procedure fmin84
module procedure fmin48
end interface
contains
real(8) function fmax88(x, y) result (res)
real(8), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(x, max(x, y), y /= y), x /= x)
end function
real(4) function fmax44(x, y) result (res)
real(4), intent (in) :: x
real(4), intent (in) :: y
res = merge(y, merge(x, max(x, y), y /= y), x /= x)
end function
real(8) function fmax84(x, y) result(res)
real(8), intent (in) :: x
real(4), intent (in) :: y
res = merge(dble(y), merge(x, max(x, dble(y)), y /= y), x /= x)
end function
real(8) function fmax48(x, y) result(res)
real(4), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(dble(x), max(dble(x), y), y /= y), x /= x)
end function
real(8) function fmin88(x, y) result (res)
real(8), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(x, min(x, y), y /= y), x /= x)
end function
real(4) function fmin44(x, y) result (res)
real(4), intent (in) :: x
real(4), intent (in) :: y
res = merge(y, merge(x, min(x, y), y /= y), x /= x)
end function
real(8) function fmin84(x, y) result(res)
real(8), intent (in) :: x
real(4), intent (in) :: y
res = merge(dble(y), merge(x, min(x, dble(y)), y /= y), x /= x)
end function
real(8) function fmin48(x, y) result(res)
real(4), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(dble(x), min(dble(x), y), y /= y), x /= x)
end function
end module
real(8) function code(re, im)
use fmin_fmax_functions
real(8), intent (in) :: re
real(8), intent (in) :: im
code = (0.5d0 * sin(re)) * (exp((0.0d0 - im)) + exp(im))
end function
public static double code(double re, double im) {
return (0.5 * Math.sin(re)) * (Math.exp((0.0 - im)) + Math.exp(im));
}
def code(re, im): return (0.5 * math.sin(re)) * (math.exp((0.0 - im)) + math.exp(im))
function code(re, im) return Float64(Float64(0.5 * sin(re)) * Float64(exp(Float64(0.0 - im)) + exp(im))) end
function tmp = code(re, im) tmp = (0.5 * sin(re)) * (exp((0.0 - im)) + exp(im)); end
code[re_, im_] := N[(N[(0.5 * N[Sin[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 \sin re\right) \cdot \left(e^{0 - im} + e^{im}\right)
\end{array}
re\_m = (fabs.f64 re)
re\_s = (copysign.f64 #s(literal 1 binary64) re)
(FPCore (re_s re_m im)
:precision binary64
(let* ((t_0 (* 2.0 (cosh im)))
(t_1 (* 0.5 (sin re_m)))
(t_2 (* t_1 (+ (exp (- im)) (exp im)))))
(*
re_s
(if (<= t_2 (- INFINITY))
(* (* t_0 (* (* re_m re_m) -0.08333333333333333)) re_m)
(if (<= t_2 1.0) (* t_1 (fma im im 2.0)) (* (* re_m 0.5) t_0))))))re\_m = fabs(re);
re\_s = copysign(1.0, re);
double code(double re_s, double re_m, double im) {
double t_0 = 2.0 * cosh(im);
double t_1 = 0.5 * sin(re_m);
double t_2 = t_1 * (exp(-im) + exp(im));
double tmp;
if (t_2 <= -((double) INFINITY)) {
tmp = (t_0 * ((re_m * re_m) * -0.08333333333333333)) * re_m;
} else if (t_2 <= 1.0) {
tmp = t_1 * fma(im, im, 2.0);
} else {
tmp = (re_m * 0.5) * t_0;
}
return re_s * tmp;
}
re\_m = abs(re) re\_s = copysign(1.0, re) function code(re_s, re_m, im) t_0 = Float64(2.0 * cosh(im)) t_1 = Float64(0.5 * sin(re_m)) t_2 = Float64(t_1 * Float64(exp(Float64(-im)) + exp(im))) tmp = 0.0 if (t_2 <= Float64(-Inf)) tmp = Float64(Float64(t_0 * Float64(Float64(re_m * re_m) * -0.08333333333333333)) * re_m); elseif (t_2 <= 1.0) tmp = Float64(t_1 * fma(im, im, 2.0)); else tmp = Float64(Float64(re_m * 0.5) * t_0); end return Float64(re_s * tmp) end
re\_m = N[Abs[re], $MachinePrecision]
re\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[re]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[re$95$s_, re$95$m_, im_] := Block[{t$95$0 = N[(2.0 * N[Cosh[im], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(0.5 * N[Sin[re$95$m], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(t$95$1 * N[(N[Exp[(-im)], $MachinePrecision] + N[Exp[im], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, N[(re$95$s * If[LessEqual[t$95$2, (-Infinity)], N[(N[(t$95$0 * N[(N[(re$95$m * re$95$m), $MachinePrecision] * -0.08333333333333333), $MachinePrecision]), $MachinePrecision] * re$95$m), $MachinePrecision], If[LessEqual[t$95$2, 1.0], N[(t$95$1 * N[(im * im + 2.0), $MachinePrecision]), $MachinePrecision], N[(N[(re$95$m * 0.5), $MachinePrecision] * t$95$0), $MachinePrecision]]]), $MachinePrecision]]]]
\begin{array}{l}
re\_m = \left|re\right|
\\
re\_s = \mathsf{copysign}\left(1, re\right)
\\
\begin{array}{l}
t_0 := 2 \cdot \cosh im\\
t_1 := 0.5 \cdot \sin re\_m\\
t_2 := t\_1 \cdot \left(e^{-im} + e^{im}\right)\\
re\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_2 \leq -\infty:\\
\;\;\;\;\left(t\_0 \cdot \left(\left(re\_m \cdot re\_m\right) \cdot -0.08333333333333333\right)\right) \cdot re\_m\\
\mathbf{elif}\;t\_2 \leq 1:\\
\;\;\;\;t\_1 \cdot \mathsf{fma}\left(im, im, 2\right)\\
\mathbf{else}:\\
\;\;\;\;\left(re\_m \cdot 0.5\right) \cdot t\_0\\
\end{array}
\end{array}
\end{array}
if (*.f64 (*.f64 #s(literal 1/2 binary64) (sin.f64 re)) (+.f64 (exp.f64 (-.f64 #s(literal 0 binary64) im)) (exp.f64 im))) < -inf.0Initial program 100.0%
Taylor expanded in re around 0
*-commutativeN/A
lower-*.f64N/A
Applied rewrites74.0%
Taylor expanded in re around inf
*-commutativeN/A
lower-*.f64N/A
pow2N/A
lift-*.f6423.4
Applied rewrites23.4%
if -inf.0 < (*.f64 (*.f64 #s(literal 1/2 binary64) (sin.f64 re)) (+.f64 (exp.f64 (-.f64 #s(literal 0 binary64) im)) (exp.f64 im))) < 1Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
unpow2N/A
lower-fma.f64100.0
Applied rewrites100.0%
if 1 < (*.f64 (*.f64 #s(literal 1/2 binary64) (sin.f64 re)) (+.f64 (exp.f64 (-.f64 #s(literal 0 binary64) im)) (exp.f64 im))) Initial program 100.0%
Taylor expanded in re around 0
associate-*r*N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
cosh-undefN/A
lower-*.f64N/A
lower-cosh.f6477.1
Applied rewrites77.1%
Final simplification70.7%
re\_m = (fabs.f64 re)
re\_s = (copysign.f64 #s(literal 1 binary64) re)
(FPCore (re_s re_m im)
:precision binary64
(let* ((t_0 (* 0.5 (sin re_m))) (t_1 (* t_0 (+ (exp (- im)) (exp im)))))
(*
re_s
(if (<= t_1 (- INFINITY))
(*
(*
(fma
(fma
(fma 0.002777777777777778 (* im im) 0.08333333333333333)
(* im im)
1.0)
(* im im)
2.0)
(fma (* re_m re_m) -0.08333333333333333 0.5))
re_m)
(if (<= t_1 1.0)
(* t_0 (fma im im 2.0))
(* (* re_m 0.5) (* 2.0 (cosh im))))))))re\_m = fabs(re);
re\_s = copysign(1.0, re);
double code(double re_s, double re_m, double im) {
double t_0 = 0.5 * sin(re_m);
double t_1 = t_0 * (exp(-im) + exp(im));
double tmp;
if (t_1 <= -((double) INFINITY)) {
tmp = (fma(fma(fma(0.002777777777777778, (im * im), 0.08333333333333333), (im * im), 1.0), (im * im), 2.0) * fma((re_m * re_m), -0.08333333333333333, 0.5)) * re_m;
} else if (t_1 <= 1.0) {
tmp = t_0 * fma(im, im, 2.0);
} else {
tmp = (re_m * 0.5) * (2.0 * cosh(im));
}
return re_s * tmp;
}
re\_m = abs(re) re\_s = copysign(1.0, re) function code(re_s, re_m, im) t_0 = Float64(0.5 * sin(re_m)) t_1 = Float64(t_0 * Float64(exp(Float64(-im)) + exp(im))) tmp = 0.0 if (t_1 <= Float64(-Inf)) tmp = Float64(Float64(fma(fma(fma(0.002777777777777778, Float64(im * im), 0.08333333333333333), Float64(im * im), 1.0), Float64(im * im), 2.0) * fma(Float64(re_m * re_m), -0.08333333333333333, 0.5)) * re_m); elseif (t_1 <= 1.0) tmp = Float64(t_0 * fma(im, im, 2.0)); else tmp = Float64(Float64(re_m * 0.5) * Float64(2.0 * cosh(im))); end return Float64(re_s * tmp) end
re\_m = N[Abs[re], $MachinePrecision]
re\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[re]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[re$95$s_, re$95$m_, im_] := Block[{t$95$0 = N[(0.5 * N[Sin[re$95$m], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(t$95$0 * N[(N[Exp[(-im)], $MachinePrecision] + N[Exp[im], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, N[(re$95$s * If[LessEqual[t$95$1, (-Infinity)], N[(N[(N[(N[(N[(0.002777777777777778 * N[(im * im), $MachinePrecision] + 0.08333333333333333), $MachinePrecision] * N[(im * im), $MachinePrecision] + 1.0), $MachinePrecision] * N[(im * im), $MachinePrecision] + 2.0), $MachinePrecision] * N[(N[(re$95$m * re$95$m), $MachinePrecision] * -0.08333333333333333 + 0.5), $MachinePrecision]), $MachinePrecision] * re$95$m), $MachinePrecision], If[LessEqual[t$95$1, 1.0], N[(t$95$0 * N[(im * im + 2.0), $MachinePrecision]), $MachinePrecision], N[(N[(re$95$m * 0.5), $MachinePrecision] * N[(2.0 * N[Cosh[im], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]]]
\begin{array}{l}
re\_m = \left|re\right|
\\
re\_s = \mathsf{copysign}\left(1, re\right)
\\
\begin{array}{l}
t_0 := 0.5 \cdot \sin re\_m\\
t_1 := t\_0 \cdot \left(e^{-im} + e^{im}\right)\\
re\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_1 \leq -\infty:\\
\;\;\;\;\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(0.002777777777777778, im \cdot im, 0.08333333333333333\right), im \cdot im, 1\right), im \cdot im, 2\right) \cdot \mathsf{fma}\left(re\_m \cdot re\_m, -0.08333333333333333, 0.5\right)\right) \cdot re\_m\\
\mathbf{elif}\;t\_1 \leq 1:\\
\;\;\;\;t\_0 \cdot \mathsf{fma}\left(im, im, 2\right)\\
\mathbf{else}:\\
\;\;\;\;\left(re\_m \cdot 0.5\right) \cdot \left(2 \cdot \cosh im\right)\\
\end{array}
\end{array}
\end{array}
if (*.f64 (*.f64 #s(literal 1/2 binary64) (sin.f64 re)) (+.f64 (exp.f64 (-.f64 #s(literal 0 binary64) im)) (exp.f64 im))) < -inf.0Initial program 100.0%
Taylor expanded in re around 0
*-commutativeN/A
lower-*.f64N/A
Applied rewrites74.0%
Taylor expanded in im around 0
cosh-undef-revN/A
sub0-negN/A
+-commutativeN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
pow2N/A
lower-*.f64N/A
pow2N/A
lower-*.f64N/A
pow2N/A
lower-*.f6465.4
Applied rewrites65.4%
if -inf.0 < (*.f64 (*.f64 #s(literal 1/2 binary64) (sin.f64 re)) (+.f64 (exp.f64 (-.f64 #s(literal 0 binary64) im)) (exp.f64 im))) < 1Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
unpow2N/A
lower-fma.f64100.0
Applied rewrites100.0%
if 1 < (*.f64 (*.f64 #s(literal 1/2 binary64) (sin.f64 re)) (+.f64 (exp.f64 (-.f64 #s(literal 0 binary64) im)) (exp.f64 im))) Initial program 100.0%
Taylor expanded in re around 0
associate-*r*N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
cosh-undefN/A
lower-*.f64N/A
lower-cosh.f6477.1
Applied rewrites77.1%
Final simplification83.3%
re\_m = (fabs.f64 re)
re\_s = (copysign.f64 #s(literal 1 binary64) re)
(FPCore (re_s re_m im)
:precision binary64
(let* ((t_0 (* (* 0.5 (sin re_m)) (+ (exp (- im)) (exp im)))))
(*
re_s
(if (<= t_0 (- INFINITY))
(*
(*
(fma
(fma
(fma 0.002777777777777778 (* im im) 0.08333333333333333)
(* im im)
1.0)
(* im im)
2.0)
(fma (* re_m re_m) -0.08333333333333333 0.5))
re_m)
(if (<= t_0 1.0) (sin re_m) (* (* re_m 0.5) (* 2.0 (cosh im))))))))re\_m = fabs(re);
re\_s = copysign(1.0, re);
double code(double re_s, double re_m, double im) {
double t_0 = (0.5 * sin(re_m)) * (exp(-im) + exp(im));
double tmp;
if (t_0 <= -((double) INFINITY)) {
tmp = (fma(fma(fma(0.002777777777777778, (im * im), 0.08333333333333333), (im * im), 1.0), (im * im), 2.0) * fma((re_m * re_m), -0.08333333333333333, 0.5)) * re_m;
} else if (t_0 <= 1.0) {
tmp = sin(re_m);
} else {
tmp = (re_m * 0.5) * (2.0 * cosh(im));
}
return re_s * tmp;
}
re\_m = abs(re) re\_s = copysign(1.0, re) function code(re_s, re_m, im) t_0 = Float64(Float64(0.5 * sin(re_m)) * Float64(exp(Float64(-im)) + exp(im))) tmp = 0.0 if (t_0 <= Float64(-Inf)) tmp = Float64(Float64(fma(fma(fma(0.002777777777777778, Float64(im * im), 0.08333333333333333), Float64(im * im), 1.0), Float64(im * im), 2.0) * fma(Float64(re_m * re_m), -0.08333333333333333, 0.5)) * re_m); elseif (t_0 <= 1.0) tmp = sin(re_m); else tmp = Float64(Float64(re_m * 0.5) * Float64(2.0 * cosh(im))); end return Float64(re_s * tmp) end
re\_m = N[Abs[re], $MachinePrecision]
re\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[re]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[re$95$s_, re$95$m_, im_] := Block[{t$95$0 = N[(N[(0.5 * N[Sin[re$95$m], $MachinePrecision]), $MachinePrecision] * N[(N[Exp[(-im)], $MachinePrecision] + N[Exp[im], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, N[(re$95$s * If[LessEqual[t$95$0, (-Infinity)], N[(N[(N[(N[(N[(0.002777777777777778 * N[(im * im), $MachinePrecision] + 0.08333333333333333), $MachinePrecision] * N[(im * im), $MachinePrecision] + 1.0), $MachinePrecision] * N[(im * im), $MachinePrecision] + 2.0), $MachinePrecision] * N[(N[(re$95$m * re$95$m), $MachinePrecision] * -0.08333333333333333 + 0.5), $MachinePrecision]), $MachinePrecision] * re$95$m), $MachinePrecision], If[LessEqual[t$95$0, 1.0], N[Sin[re$95$m], $MachinePrecision], N[(N[(re$95$m * 0.5), $MachinePrecision] * N[(2.0 * N[Cosh[im], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]]
\begin{array}{l}
re\_m = \left|re\right|
\\
re\_s = \mathsf{copysign}\left(1, re\right)
\\
\begin{array}{l}
t_0 := \left(0.5 \cdot \sin re\_m\right) \cdot \left(e^{-im} + e^{im}\right)\\
re\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_0 \leq -\infty:\\
\;\;\;\;\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(0.002777777777777778, im \cdot im, 0.08333333333333333\right), im \cdot im, 1\right), im \cdot im, 2\right) \cdot \mathsf{fma}\left(re\_m \cdot re\_m, -0.08333333333333333, 0.5\right)\right) \cdot re\_m\\
\mathbf{elif}\;t\_0 \leq 1:\\
\;\;\;\;\sin re\_m\\
\mathbf{else}:\\
\;\;\;\;\left(re\_m \cdot 0.5\right) \cdot \left(2 \cdot \cosh im\right)\\
\end{array}
\end{array}
\end{array}
if (*.f64 (*.f64 #s(literal 1/2 binary64) (sin.f64 re)) (+.f64 (exp.f64 (-.f64 #s(literal 0 binary64) im)) (exp.f64 im))) < -inf.0Initial program 100.0%
Taylor expanded in re around 0
*-commutativeN/A
lower-*.f64N/A
Applied rewrites74.0%
Taylor expanded in im around 0
cosh-undef-revN/A
sub0-negN/A
+-commutativeN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
pow2N/A
lower-*.f64N/A
pow2N/A
lower-*.f64N/A
pow2N/A
lower-*.f6465.4
Applied rewrites65.4%
if -inf.0 < (*.f64 (*.f64 #s(literal 1/2 binary64) (sin.f64 re)) (+.f64 (exp.f64 (-.f64 #s(literal 0 binary64) im)) (exp.f64 im))) < 1Initial program 100.0%
Taylor expanded in im around 0
lift-sin.f6499.3
Applied rewrites99.3%
if 1 < (*.f64 (*.f64 #s(literal 1/2 binary64) (sin.f64 re)) (+.f64 (exp.f64 (-.f64 #s(literal 0 binary64) im)) (exp.f64 im))) Initial program 100.0%
Taylor expanded in re around 0
associate-*r*N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
cosh-undefN/A
lower-*.f64N/A
lower-cosh.f6477.1
Applied rewrites77.1%
Final simplification83.0%
re\_m = (fabs.f64 re)
re\_s = (copysign.f64 #s(literal 1 binary64) re)
(FPCore (re_s re_m im)
:precision binary64
(let* ((t_0 (* (* 0.5 (sin re_m)) (+ (exp (- im)) (exp im)))))
(*
re_s
(if (<= t_0 (- INFINITY))
(*
(*
(fma
(fma
(fma 0.002777777777777778 (* im im) 0.08333333333333333)
(* im im)
1.0)
(* im im)
2.0)
(fma (* re_m re_m) -0.08333333333333333 0.5))
re_m)
(if (<= t_0 1.0)
(sin re_m)
(*
(* re_m 0.5)
(fma
(*
(fma
(* (fma (* 0.002777777777777778 im) im 0.08333333333333333) im)
im
1.0)
im)
im
2.0)))))))re\_m = fabs(re);
re\_s = copysign(1.0, re);
double code(double re_s, double re_m, double im) {
double t_0 = (0.5 * sin(re_m)) * (exp(-im) + exp(im));
double tmp;
if (t_0 <= -((double) INFINITY)) {
tmp = (fma(fma(fma(0.002777777777777778, (im * im), 0.08333333333333333), (im * im), 1.0), (im * im), 2.0) * fma((re_m * re_m), -0.08333333333333333, 0.5)) * re_m;
} else if (t_0 <= 1.0) {
tmp = sin(re_m);
} else {
tmp = (re_m * 0.5) * fma((fma((fma((0.002777777777777778 * im), im, 0.08333333333333333) * im), im, 1.0) * im), im, 2.0);
}
return re_s * tmp;
}
re\_m = abs(re) re\_s = copysign(1.0, re) function code(re_s, re_m, im) t_0 = Float64(Float64(0.5 * sin(re_m)) * Float64(exp(Float64(-im)) + exp(im))) tmp = 0.0 if (t_0 <= Float64(-Inf)) tmp = Float64(Float64(fma(fma(fma(0.002777777777777778, Float64(im * im), 0.08333333333333333), Float64(im * im), 1.0), Float64(im * im), 2.0) * fma(Float64(re_m * re_m), -0.08333333333333333, 0.5)) * re_m); elseif (t_0 <= 1.0) tmp = sin(re_m); else tmp = Float64(Float64(re_m * 0.5) * fma(Float64(fma(Float64(fma(Float64(0.002777777777777778 * im), im, 0.08333333333333333) * im), im, 1.0) * im), im, 2.0)); end return Float64(re_s * tmp) end
re\_m = N[Abs[re], $MachinePrecision]
re\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[re]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[re$95$s_, re$95$m_, im_] := Block[{t$95$0 = N[(N[(0.5 * N[Sin[re$95$m], $MachinePrecision]), $MachinePrecision] * N[(N[Exp[(-im)], $MachinePrecision] + N[Exp[im], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, N[(re$95$s * If[LessEqual[t$95$0, (-Infinity)], N[(N[(N[(N[(N[(0.002777777777777778 * N[(im * im), $MachinePrecision] + 0.08333333333333333), $MachinePrecision] * N[(im * im), $MachinePrecision] + 1.0), $MachinePrecision] * N[(im * im), $MachinePrecision] + 2.0), $MachinePrecision] * N[(N[(re$95$m * re$95$m), $MachinePrecision] * -0.08333333333333333 + 0.5), $MachinePrecision]), $MachinePrecision] * re$95$m), $MachinePrecision], If[LessEqual[t$95$0, 1.0], N[Sin[re$95$m], $MachinePrecision], N[(N[(re$95$m * 0.5), $MachinePrecision] * N[(N[(N[(N[(N[(N[(0.002777777777777778 * im), $MachinePrecision] * im + 0.08333333333333333), $MachinePrecision] * im), $MachinePrecision] * im + 1.0), $MachinePrecision] * im), $MachinePrecision] * im + 2.0), $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]]
\begin{array}{l}
re\_m = \left|re\right|
\\
re\_s = \mathsf{copysign}\left(1, re\right)
\\
\begin{array}{l}
t_0 := \left(0.5 \cdot \sin re\_m\right) \cdot \left(e^{-im} + e^{im}\right)\\
re\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_0 \leq -\infty:\\
\;\;\;\;\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(0.002777777777777778, im \cdot im, 0.08333333333333333\right), im \cdot im, 1\right), im \cdot im, 2\right) \cdot \mathsf{fma}\left(re\_m \cdot re\_m, -0.08333333333333333, 0.5\right)\right) \cdot re\_m\\
\mathbf{elif}\;t\_0 \leq 1:\\
\;\;\;\;\sin re\_m\\
\mathbf{else}:\\
\;\;\;\;\left(re\_m \cdot 0.5\right) \cdot \mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(0.002777777777777778 \cdot im, im, 0.08333333333333333\right) \cdot im, im, 1\right) \cdot im, im, 2\right)\\
\end{array}
\end{array}
\end{array}
if (*.f64 (*.f64 #s(literal 1/2 binary64) (sin.f64 re)) (+.f64 (exp.f64 (-.f64 #s(literal 0 binary64) im)) (exp.f64 im))) < -inf.0Initial program 100.0%
Taylor expanded in re around 0
*-commutativeN/A
lower-*.f64N/A
Applied rewrites74.0%
Taylor expanded in im around 0
cosh-undef-revN/A
sub0-negN/A
+-commutativeN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
pow2N/A
lower-*.f64N/A
pow2N/A
lower-*.f64N/A
pow2N/A
lower-*.f6465.4
Applied rewrites65.4%
if -inf.0 < (*.f64 (*.f64 #s(literal 1/2 binary64) (sin.f64 re)) (+.f64 (exp.f64 (-.f64 #s(literal 0 binary64) im)) (exp.f64 im))) < 1Initial program 100.0%
Taylor expanded in im around 0
lift-sin.f6499.3
Applied rewrites99.3%
if 1 < (*.f64 (*.f64 #s(literal 1/2 binary64) (sin.f64 re)) (+.f64 (exp.f64 (-.f64 #s(literal 0 binary64) im)) (exp.f64 im))) Initial program 100.0%
Taylor expanded in re around 0
associate-*r*N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
cosh-undefN/A
lower-*.f64N/A
lower-cosh.f6477.1
Applied rewrites77.1%
Taylor expanded in im around 0
cosh-undef-revN/A
sub0-negN/A
+-commutativeN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
pow2N/A
lower-*.f64N/A
pow2N/A
lower-*.f64N/A
pow2N/A
lower-*.f6467.5
Applied rewrites67.5%
lift-*.f64N/A
lift-fma.f64N/A
lift-*.f64N/A
lift-fma.f64N/A
lift-*.f64N/A
lift-fma.f64N/A
associate-*r*N/A
lower-fma.f64N/A
Applied rewrites67.5%
lift-*.f64N/A
lift-*.f64N/A
lift-fma.f64N/A
*-commutativeN/A
associate-*l*N/A
lower-fma.f64N/A
associate-*r*N/A
lower-fma.f64N/A
lower-*.f64N/A
lift-fma.f64N/A
lift-*.f6467.5
Applied rewrites67.5%
Final simplification80.4%
re\_m = (fabs.f64 re)
re\_s = (copysign.f64 #s(literal 1 binary64) re)
(FPCore (re_s re_m im)
:precision binary64
(*
re_s
(if (<= (* (* 0.5 (sin re_m)) (+ (exp (- im)) (exp im))) 0.0002)
(*
(*
(fma (fma im (* im 0.08333333333333333) 1.0) (* im im) 2.0)
(fma (* re_m re_m) -0.08333333333333333 0.5))
re_m)
(*
(* re_m 0.5)
(fma
(*
(fma
(* (fma (* 0.002777777777777778 im) im 0.08333333333333333) im)
im
1.0)
im)
im
2.0)))))re\_m = fabs(re);
re\_s = copysign(1.0, re);
double code(double re_s, double re_m, double im) {
double tmp;
if (((0.5 * sin(re_m)) * (exp(-im) + exp(im))) <= 0.0002) {
tmp = (fma(fma(im, (im * 0.08333333333333333), 1.0), (im * im), 2.0) * fma((re_m * re_m), -0.08333333333333333, 0.5)) * re_m;
} else {
tmp = (re_m * 0.5) * fma((fma((fma((0.002777777777777778 * im), im, 0.08333333333333333) * im), im, 1.0) * im), im, 2.0);
}
return re_s * tmp;
}
re\_m = abs(re) re\_s = copysign(1.0, re) function code(re_s, re_m, im) tmp = 0.0 if (Float64(Float64(0.5 * sin(re_m)) * Float64(exp(Float64(-im)) + exp(im))) <= 0.0002) tmp = Float64(Float64(fma(fma(im, Float64(im * 0.08333333333333333), 1.0), Float64(im * im), 2.0) * fma(Float64(re_m * re_m), -0.08333333333333333, 0.5)) * re_m); else tmp = Float64(Float64(re_m * 0.5) * fma(Float64(fma(Float64(fma(Float64(0.002777777777777778 * im), im, 0.08333333333333333) * im), im, 1.0) * im), im, 2.0)); end return Float64(re_s * tmp) end
re\_m = N[Abs[re], $MachinePrecision]
re\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[re]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[re$95$s_, re$95$m_, im_] := N[(re$95$s * If[LessEqual[N[(N[(0.5 * N[Sin[re$95$m], $MachinePrecision]), $MachinePrecision] * N[(N[Exp[(-im)], $MachinePrecision] + N[Exp[im], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 0.0002], N[(N[(N[(N[(im * N[(im * 0.08333333333333333), $MachinePrecision] + 1.0), $MachinePrecision] * N[(im * im), $MachinePrecision] + 2.0), $MachinePrecision] * N[(N[(re$95$m * re$95$m), $MachinePrecision] * -0.08333333333333333 + 0.5), $MachinePrecision]), $MachinePrecision] * re$95$m), $MachinePrecision], N[(N[(re$95$m * 0.5), $MachinePrecision] * N[(N[(N[(N[(N[(N[(0.002777777777777778 * im), $MachinePrecision] * im + 0.08333333333333333), $MachinePrecision] * im), $MachinePrecision] * im + 1.0), $MachinePrecision] * im), $MachinePrecision] * im + 2.0), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
re\_m = \left|re\right|
\\
re\_s = \mathsf{copysign}\left(1, re\right)
\\
re\_s \cdot \begin{array}{l}
\mathbf{if}\;\left(0.5 \cdot \sin re\_m\right) \cdot \left(e^{-im} + e^{im}\right) \leq 0.0002:\\
\;\;\;\;\left(\mathsf{fma}\left(\mathsf{fma}\left(im, im \cdot 0.08333333333333333, 1\right), im \cdot im, 2\right) \cdot \mathsf{fma}\left(re\_m \cdot re\_m, -0.08333333333333333, 0.5\right)\right) \cdot re\_m\\
\mathbf{else}:\\
\;\;\;\;\left(re\_m \cdot 0.5\right) \cdot \mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(0.002777777777777778 \cdot im, im, 0.08333333333333333\right) \cdot im, im, 1\right) \cdot im, im, 2\right)\\
\end{array}
\end{array}
if (*.f64 (*.f64 #s(literal 1/2 binary64) (sin.f64 re)) (+.f64 (exp.f64 (-.f64 #s(literal 0 binary64) im)) (exp.f64 im))) < 2.0000000000000001e-4Initial program 100.0%
Taylor expanded in re around 0
*-commutativeN/A
lower-*.f64N/A
Applied rewrites68.4%
Taylor expanded in im around 0
cosh-undef-revN/A
sub0-negN/A
+-commutativeN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
pow2N/A
lower-*.f64N/A
pow2N/A
lower-*.f6462.2
Applied rewrites62.2%
lift-*.f64N/A
lift-fma.f64N/A
associate-*l*N/A
lower-fma.f64N/A
lower-*.f6462.2
Applied rewrites62.2%
if 2.0000000000000001e-4 < (*.f64 (*.f64 #s(literal 1/2 binary64) (sin.f64 re)) (+.f64 (exp.f64 (-.f64 #s(literal 0 binary64) im)) (exp.f64 im))) Initial program 100.0%
Taylor expanded in re around 0
associate-*r*N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
cosh-undefN/A
lower-*.f64N/A
lower-cosh.f6454.8
Applied rewrites54.8%
Taylor expanded in im around 0
cosh-undef-revN/A
sub0-negN/A
+-commutativeN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
pow2N/A
lower-*.f64N/A
pow2N/A
lower-*.f64N/A
pow2N/A
lower-*.f6448.0
Applied rewrites48.0%
lift-*.f64N/A
lift-fma.f64N/A
lift-*.f64N/A
lift-fma.f64N/A
lift-*.f64N/A
lift-fma.f64N/A
associate-*r*N/A
lower-fma.f64N/A
Applied rewrites48.0%
lift-*.f64N/A
lift-*.f64N/A
lift-fma.f64N/A
*-commutativeN/A
associate-*l*N/A
lower-fma.f64N/A
associate-*r*N/A
lower-fma.f64N/A
lower-*.f64N/A
lift-fma.f64N/A
lift-*.f6448.0
Applied rewrites48.0%
Final simplification56.7%
re\_m = (fabs.f64 re)
re\_s = (copysign.f64 #s(literal 1 binary64) re)
(FPCore (re_s re_m im)
:precision binary64
(*
re_s
(if (<= (* (* 0.5 (sin re_m)) (+ (exp (- im)) (exp im))) 0.0002)
(* (* (fma (* re_m re_m) -0.08333333333333333 0.5) re_m) (fma im im 2.0))
(*
(* re_m 0.5)
(fma
(*
(fma
(* (fma (* 0.002777777777777778 im) im 0.08333333333333333) im)
im
1.0)
im)
im
2.0)))))re\_m = fabs(re);
re\_s = copysign(1.0, re);
double code(double re_s, double re_m, double im) {
double tmp;
if (((0.5 * sin(re_m)) * (exp(-im) + exp(im))) <= 0.0002) {
tmp = (fma((re_m * re_m), -0.08333333333333333, 0.5) * re_m) * fma(im, im, 2.0);
} else {
tmp = (re_m * 0.5) * fma((fma((fma((0.002777777777777778 * im), im, 0.08333333333333333) * im), im, 1.0) * im), im, 2.0);
}
return re_s * tmp;
}
re\_m = abs(re) re\_s = copysign(1.0, re) function code(re_s, re_m, im) tmp = 0.0 if (Float64(Float64(0.5 * sin(re_m)) * Float64(exp(Float64(-im)) + exp(im))) <= 0.0002) tmp = Float64(Float64(fma(Float64(re_m * re_m), -0.08333333333333333, 0.5) * re_m) * fma(im, im, 2.0)); else tmp = Float64(Float64(re_m * 0.5) * fma(Float64(fma(Float64(fma(Float64(0.002777777777777778 * im), im, 0.08333333333333333) * im), im, 1.0) * im), im, 2.0)); end return Float64(re_s * tmp) end
re\_m = N[Abs[re], $MachinePrecision]
re\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[re]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[re$95$s_, re$95$m_, im_] := N[(re$95$s * If[LessEqual[N[(N[(0.5 * N[Sin[re$95$m], $MachinePrecision]), $MachinePrecision] * N[(N[Exp[(-im)], $MachinePrecision] + N[Exp[im], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 0.0002], N[(N[(N[(N[(re$95$m * re$95$m), $MachinePrecision] * -0.08333333333333333 + 0.5), $MachinePrecision] * re$95$m), $MachinePrecision] * N[(im * im + 2.0), $MachinePrecision]), $MachinePrecision], N[(N[(re$95$m * 0.5), $MachinePrecision] * N[(N[(N[(N[(N[(N[(0.002777777777777778 * im), $MachinePrecision] * im + 0.08333333333333333), $MachinePrecision] * im), $MachinePrecision] * im + 1.0), $MachinePrecision] * im), $MachinePrecision] * im + 2.0), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
re\_m = \left|re\right|
\\
re\_s = \mathsf{copysign}\left(1, re\right)
\\
re\_s \cdot \begin{array}{l}
\mathbf{if}\;\left(0.5 \cdot \sin re\_m\right) \cdot \left(e^{-im} + e^{im}\right) \leq 0.0002:\\
\;\;\;\;\left(\mathsf{fma}\left(re\_m \cdot re\_m, -0.08333333333333333, 0.5\right) \cdot re\_m\right) \cdot \mathsf{fma}\left(im, im, 2\right)\\
\mathbf{else}:\\
\;\;\;\;\left(re\_m \cdot 0.5\right) \cdot \mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(0.002777777777777778 \cdot im, im, 0.08333333333333333\right) \cdot im, im, 1\right) \cdot im, im, 2\right)\\
\end{array}
\end{array}
if (*.f64 (*.f64 #s(literal 1/2 binary64) (sin.f64 re)) (+.f64 (exp.f64 (-.f64 #s(literal 0 binary64) im)) (exp.f64 im))) < 2.0000000000000001e-4Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
unpow2N/A
lower-fma.f6474.8
Applied rewrites74.8%
Taylor expanded in re around 0
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
pow2N/A
+-commutativeN/A
lift-fma.f64N/A
lift-*.f6455.4
Applied rewrites55.4%
if 2.0000000000000001e-4 < (*.f64 (*.f64 #s(literal 1/2 binary64) (sin.f64 re)) (+.f64 (exp.f64 (-.f64 #s(literal 0 binary64) im)) (exp.f64 im))) Initial program 100.0%
Taylor expanded in re around 0
associate-*r*N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
cosh-undefN/A
lower-*.f64N/A
lower-cosh.f6454.8
Applied rewrites54.8%
Taylor expanded in im around 0
cosh-undef-revN/A
sub0-negN/A
+-commutativeN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
pow2N/A
lower-*.f64N/A
pow2N/A
lower-*.f64N/A
pow2N/A
lower-*.f6448.0
Applied rewrites48.0%
lift-*.f64N/A
lift-fma.f64N/A
lift-*.f64N/A
lift-fma.f64N/A
lift-*.f64N/A
lift-fma.f64N/A
associate-*r*N/A
lower-fma.f64N/A
Applied rewrites48.0%
lift-*.f64N/A
lift-*.f64N/A
lift-fma.f64N/A
*-commutativeN/A
associate-*l*N/A
lower-fma.f64N/A
associate-*r*N/A
lower-fma.f64N/A
lower-*.f64N/A
lift-fma.f64N/A
lift-*.f6448.0
Applied rewrites48.0%
Final simplification52.5%
re\_m = (fabs.f64 re)
re\_s = (copysign.f64 #s(literal 1 binary64) re)
(FPCore (re_s re_m im)
:precision binary64
(*
re_s
(if (<= (* (* 0.5 (sin re_m)) (+ (exp (- im)) (exp im))) 0.0002)
(* (* (fma (* re_m re_m) -0.08333333333333333 0.5) re_m) (fma im im 2.0))
(*
(* re_m 0.5)
(fma
(* (fma (* (* 0.002777777777777778 im) im) (* im im) 1.0) im)
im
2.0)))))re\_m = fabs(re);
re\_s = copysign(1.0, re);
double code(double re_s, double re_m, double im) {
double tmp;
if (((0.5 * sin(re_m)) * (exp(-im) + exp(im))) <= 0.0002) {
tmp = (fma((re_m * re_m), -0.08333333333333333, 0.5) * re_m) * fma(im, im, 2.0);
} else {
tmp = (re_m * 0.5) * fma((fma(((0.002777777777777778 * im) * im), (im * im), 1.0) * im), im, 2.0);
}
return re_s * tmp;
}
re\_m = abs(re) re\_s = copysign(1.0, re) function code(re_s, re_m, im) tmp = 0.0 if (Float64(Float64(0.5 * sin(re_m)) * Float64(exp(Float64(-im)) + exp(im))) <= 0.0002) tmp = Float64(Float64(fma(Float64(re_m * re_m), -0.08333333333333333, 0.5) * re_m) * fma(im, im, 2.0)); else tmp = Float64(Float64(re_m * 0.5) * fma(Float64(fma(Float64(Float64(0.002777777777777778 * im) * im), Float64(im * im), 1.0) * im), im, 2.0)); end return Float64(re_s * tmp) end
re\_m = N[Abs[re], $MachinePrecision]
re\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[re]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[re$95$s_, re$95$m_, im_] := N[(re$95$s * If[LessEqual[N[(N[(0.5 * N[Sin[re$95$m], $MachinePrecision]), $MachinePrecision] * N[(N[Exp[(-im)], $MachinePrecision] + N[Exp[im], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 0.0002], N[(N[(N[(N[(re$95$m * re$95$m), $MachinePrecision] * -0.08333333333333333 + 0.5), $MachinePrecision] * re$95$m), $MachinePrecision] * N[(im * im + 2.0), $MachinePrecision]), $MachinePrecision], N[(N[(re$95$m * 0.5), $MachinePrecision] * N[(N[(N[(N[(N[(0.002777777777777778 * im), $MachinePrecision] * im), $MachinePrecision] * N[(im * im), $MachinePrecision] + 1.0), $MachinePrecision] * im), $MachinePrecision] * im + 2.0), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
re\_m = \left|re\right|
\\
re\_s = \mathsf{copysign}\left(1, re\right)
\\
re\_s \cdot \begin{array}{l}
\mathbf{if}\;\left(0.5 \cdot \sin re\_m\right) \cdot \left(e^{-im} + e^{im}\right) \leq 0.0002:\\
\;\;\;\;\left(\mathsf{fma}\left(re\_m \cdot re\_m, -0.08333333333333333, 0.5\right) \cdot re\_m\right) \cdot \mathsf{fma}\left(im, im, 2\right)\\
\mathbf{else}:\\
\;\;\;\;\left(re\_m \cdot 0.5\right) \cdot \mathsf{fma}\left(\mathsf{fma}\left(\left(0.002777777777777778 \cdot im\right) \cdot im, im \cdot im, 1\right) \cdot im, im, 2\right)\\
\end{array}
\end{array}
if (*.f64 (*.f64 #s(literal 1/2 binary64) (sin.f64 re)) (+.f64 (exp.f64 (-.f64 #s(literal 0 binary64) im)) (exp.f64 im))) < 2.0000000000000001e-4Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
unpow2N/A
lower-fma.f6474.8
Applied rewrites74.8%
Taylor expanded in re around 0
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
pow2N/A
+-commutativeN/A
lift-fma.f64N/A
lift-*.f6455.4
Applied rewrites55.4%
if 2.0000000000000001e-4 < (*.f64 (*.f64 #s(literal 1/2 binary64) (sin.f64 re)) (+.f64 (exp.f64 (-.f64 #s(literal 0 binary64) im)) (exp.f64 im))) Initial program 100.0%
Taylor expanded in re around 0
associate-*r*N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
cosh-undefN/A
lower-*.f64N/A
lower-cosh.f6454.8
Applied rewrites54.8%
Taylor expanded in im around 0
cosh-undef-revN/A
sub0-negN/A
+-commutativeN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
pow2N/A
lower-*.f64N/A
pow2N/A
lower-*.f64N/A
pow2N/A
lower-*.f6448.0
Applied rewrites48.0%
lift-*.f64N/A
lift-fma.f64N/A
lift-*.f64N/A
lift-fma.f64N/A
lift-*.f64N/A
lift-fma.f64N/A
associate-*r*N/A
lower-fma.f64N/A
Applied rewrites48.0%
Taylor expanded in im around inf
*-commutativeN/A
associate-*l*N/A
pow2N/A
associate-*l*N/A
lower-*.f64N/A
lift-*.f6448.0
Applied rewrites48.0%
Final simplification52.5%
re\_m = (fabs.f64 re)
re\_s = (copysign.f64 #s(literal 1 binary64) re)
(FPCore (re_s re_m im)
:precision binary64
(*
re_s
(if (<= (* (* 0.5 (sin re_m)) (+ (exp (- im)) (exp im))) 0.0002)
(* (* (fma (* re_m re_m) -0.08333333333333333 0.5) re_m) (fma im im 2.0))
(*
(* re_m 0.5)
(fma (* (fma (* im im) 0.08333333333333333 1.0) im) im 2.0)))))re\_m = fabs(re);
re\_s = copysign(1.0, re);
double code(double re_s, double re_m, double im) {
double tmp;
if (((0.5 * sin(re_m)) * (exp(-im) + exp(im))) <= 0.0002) {
tmp = (fma((re_m * re_m), -0.08333333333333333, 0.5) * re_m) * fma(im, im, 2.0);
} else {
tmp = (re_m * 0.5) * fma((fma((im * im), 0.08333333333333333, 1.0) * im), im, 2.0);
}
return re_s * tmp;
}
re\_m = abs(re) re\_s = copysign(1.0, re) function code(re_s, re_m, im) tmp = 0.0 if (Float64(Float64(0.5 * sin(re_m)) * Float64(exp(Float64(-im)) + exp(im))) <= 0.0002) tmp = Float64(Float64(fma(Float64(re_m * re_m), -0.08333333333333333, 0.5) * re_m) * fma(im, im, 2.0)); else tmp = Float64(Float64(re_m * 0.5) * fma(Float64(fma(Float64(im * im), 0.08333333333333333, 1.0) * im), im, 2.0)); end return Float64(re_s * tmp) end
re\_m = N[Abs[re], $MachinePrecision]
re\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[re]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[re$95$s_, re$95$m_, im_] := N[(re$95$s * If[LessEqual[N[(N[(0.5 * N[Sin[re$95$m], $MachinePrecision]), $MachinePrecision] * N[(N[Exp[(-im)], $MachinePrecision] + N[Exp[im], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 0.0002], N[(N[(N[(N[(re$95$m * re$95$m), $MachinePrecision] * -0.08333333333333333 + 0.5), $MachinePrecision] * re$95$m), $MachinePrecision] * N[(im * im + 2.0), $MachinePrecision]), $MachinePrecision], N[(N[(re$95$m * 0.5), $MachinePrecision] * N[(N[(N[(N[(im * im), $MachinePrecision] * 0.08333333333333333 + 1.0), $MachinePrecision] * im), $MachinePrecision] * im + 2.0), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
re\_m = \left|re\right|
\\
re\_s = \mathsf{copysign}\left(1, re\right)
\\
re\_s \cdot \begin{array}{l}
\mathbf{if}\;\left(0.5 \cdot \sin re\_m\right) \cdot \left(e^{-im} + e^{im}\right) \leq 0.0002:\\
\;\;\;\;\left(\mathsf{fma}\left(re\_m \cdot re\_m, -0.08333333333333333, 0.5\right) \cdot re\_m\right) \cdot \mathsf{fma}\left(im, im, 2\right)\\
\mathbf{else}:\\
\;\;\;\;\left(re\_m \cdot 0.5\right) \cdot \mathsf{fma}\left(\mathsf{fma}\left(im \cdot im, 0.08333333333333333, 1\right) \cdot im, im, 2\right)\\
\end{array}
\end{array}
if (*.f64 (*.f64 #s(literal 1/2 binary64) (sin.f64 re)) (+.f64 (exp.f64 (-.f64 #s(literal 0 binary64) im)) (exp.f64 im))) < 2.0000000000000001e-4Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
unpow2N/A
lower-fma.f6474.8
Applied rewrites74.8%
Taylor expanded in re around 0
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
pow2N/A
+-commutativeN/A
lift-fma.f64N/A
lift-*.f6455.4
Applied rewrites55.4%
if 2.0000000000000001e-4 < (*.f64 (*.f64 #s(literal 1/2 binary64) (sin.f64 re)) (+.f64 (exp.f64 (-.f64 #s(literal 0 binary64) im)) (exp.f64 im))) Initial program 100.0%
Taylor expanded in re around 0
associate-*r*N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
cosh-undefN/A
lower-*.f64N/A
lower-cosh.f6454.8
Applied rewrites54.8%
Taylor expanded in im around 0
cosh-undef-revN/A
sub0-negN/A
+-commutativeN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
pow2N/A
lower-*.f64N/A
pow2N/A
lower-*.f6445.2
Applied rewrites45.2%
lift-*.f64N/A
lift-fma.f64N/A
lift-*.f64N/A
lift-fma.f64N/A
associate-*r*N/A
lower-fma.f64N/A
lower-*.f64N/A
lift-fma.f64N/A
lift-*.f6445.2
Applied rewrites45.2%
Final simplification51.4%
re\_m = (fabs.f64 re)
re\_s = (copysign.f64 #s(literal 1 binary64) re)
(FPCore (re_s re_m im)
:precision binary64
(*
re_s
(if (<= (* (* 0.5 (sin re_m)) (+ (exp (- im)) (exp im))) 0.0002)
(* (* (fma (* re_m re_m) -0.08333333333333333 0.5) re_m) (fma im im 2.0))
(* (* re_m 0.5) (fma (* (* im im) 0.08333333333333333) (* im im) 2.0)))))re\_m = fabs(re);
re\_s = copysign(1.0, re);
double code(double re_s, double re_m, double im) {
double tmp;
if (((0.5 * sin(re_m)) * (exp(-im) + exp(im))) <= 0.0002) {
tmp = (fma((re_m * re_m), -0.08333333333333333, 0.5) * re_m) * fma(im, im, 2.0);
} else {
tmp = (re_m * 0.5) * fma(((im * im) * 0.08333333333333333), (im * im), 2.0);
}
return re_s * tmp;
}
re\_m = abs(re) re\_s = copysign(1.0, re) function code(re_s, re_m, im) tmp = 0.0 if (Float64(Float64(0.5 * sin(re_m)) * Float64(exp(Float64(-im)) + exp(im))) <= 0.0002) tmp = Float64(Float64(fma(Float64(re_m * re_m), -0.08333333333333333, 0.5) * re_m) * fma(im, im, 2.0)); else tmp = Float64(Float64(re_m * 0.5) * fma(Float64(Float64(im * im) * 0.08333333333333333), Float64(im * im), 2.0)); end return Float64(re_s * tmp) end
re\_m = N[Abs[re], $MachinePrecision]
re\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[re]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[re$95$s_, re$95$m_, im_] := N[(re$95$s * If[LessEqual[N[(N[(0.5 * N[Sin[re$95$m], $MachinePrecision]), $MachinePrecision] * N[(N[Exp[(-im)], $MachinePrecision] + N[Exp[im], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 0.0002], N[(N[(N[(N[(re$95$m * re$95$m), $MachinePrecision] * -0.08333333333333333 + 0.5), $MachinePrecision] * re$95$m), $MachinePrecision] * N[(im * im + 2.0), $MachinePrecision]), $MachinePrecision], N[(N[(re$95$m * 0.5), $MachinePrecision] * N[(N[(N[(im * im), $MachinePrecision] * 0.08333333333333333), $MachinePrecision] * N[(im * im), $MachinePrecision] + 2.0), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
re\_m = \left|re\right|
\\
re\_s = \mathsf{copysign}\left(1, re\right)
\\
re\_s \cdot \begin{array}{l}
\mathbf{if}\;\left(0.5 \cdot \sin re\_m\right) \cdot \left(e^{-im} + e^{im}\right) \leq 0.0002:\\
\;\;\;\;\left(\mathsf{fma}\left(re\_m \cdot re\_m, -0.08333333333333333, 0.5\right) \cdot re\_m\right) \cdot \mathsf{fma}\left(im, im, 2\right)\\
\mathbf{else}:\\
\;\;\;\;\left(re\_m \cdot 0.5\right) \cdot \mathsf{fma}\left(\left(im \cdot im\right) \cdot 0.08333333333333333, im \cdot im, 2\right)\\
\end{array}
\end{array}
if (*.f64 (*.f64 #s(literal 1/2 binary64) (sin.f64 re)) (+.f64 (exp.f64 (-.f64 #s(literal 0 binary64) im)) (exp.f64 im))) < 2.0000000000000001e-4Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
unpow2N/A
lower-fma.f6474.8
Applied rewrites74.8%
Taylor expanded in re around 0
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
pow2N/A
+-commutativeN/A
lift-fma.f64N/A
lift-*.f6455.4
Applied rewrites55.4%
if 2.0000000000000001e-4 < (*.f64 (*.f64 #s(literal 1/2 binary64) (sin.f64 re)) (+.f64 (exp.f64 (-.f64 #s(literal 0 binary64) im)) (exp.f64 im))) Initial program 100.0%
Taylor expanded in re around 0
associate-*r*N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
cosh-undefN/A
lower-*.f64N/A
lower-cosh.f6454.8
Applied rewrites54.8%
Taylor expanded in im around 0
cosh-undef-revN/A
sub0-negN/A
+-commutativeN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
pow2N/A
lower-*.f64N/A
pow2N/A
lower-*.f6445.2
Applied rewrites45.2%
Taylor expanded in im around inf
*-commutativeN/A
lower-*.f64N/A
pow2N/A
lift-*.f6445.2
Applied rewrites45.2%
Final simplification51.4%
re\_m = (fabs.f64 re)
re\_s = (copysign.f64 #s(literal 1 binary64) re)
(FPCore (re_s re_m im)
:precision binary64
(*
re_s
(if (<= (* (* 0.5 (sin re_m)) (+ (exp (- im)) (exp im))) 0.0002)
(* (* (fma (* re_m re_m) -0.08333333333333333 0.5) re_m) (fma im im 2.0))
(* (* 0.5 re_m) (fma im im 2.0)))))re\_m = fabs(re);
re\_s = copysign(1.0, re);
double code(double re_s, double re_m, double im) {
double tmp;
if (((0.5 * sin(re_m)) * (exp(-im) + exp(im))) <= 0.0002) {
tmp = (fma((re_m * re_m), -0.08333333333333333, 0.5) * re_m) * fma(im, im, 2.0);
} else {
tmp = (0.5 * re_m) * fma(im, im, 2.0);
}
return re_s * tmp;
}
re\_m = abs(re) re\_s = copysign(1.0, re) function code(re_s, re_m, im) tmp = 0.0 if (Float64(Float64(0.5 * sin(re_m)) * Float64(exp(Float64(-im)) + exp(im))) <= 0.0002) tmp = Float64(Float64(fma(Float64(re_m * re_m), -0.08333333333333333, 0.5) * re_m) * fma(im, im, 2.0)); else tmp = Float64(Float64(0.5 * re_m) * fma(im, im, 2.0)); end return Float64(re_s * tmp) end
re\_m = N[Abs[re], $MachinePrecision]
re\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[re]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[re$95$s_, re$95$m_, im_] := N[(re$95$s * If[LessEqual[N[(N[(0.5 * N[Sin[re$95$m], $MachinePrecision]), $MachinePrecision] * N[(N[Exp[(-im)], $MachinePrecision] + N[Exp[im], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 0.0002], N[(N[(N[(N[(re$95$m * re$95$m), $MachinePrecision] * -0.08333333333333333 + 0.5), $MachinePrecision] * re$95$m), $MachinePrecision] * N[(im * im + 2.0), $MachinePrecision]), $MachinePrecision], N[(N[(0.5 * re$95$m), $MachinePrecision] * N[(im * im + 2.0), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
re\_m = \left|re\right|
\\
re\_s = \mathsf{copysign}\left(1, re\right)
\\
re\_s \cdot \begin{array}{l}
\mathbf{if}\;\left(0.5 \cdot \sin re\_m\right) \cdot \left(e^{-im} + e^{im}\right) \leq 0.0002:\\
\;\;\;\;\left(\mathsf{fma}\left(re\_m \cdot re\_m, -0.08333333333333333, 0.5\right) \cdot re\_m\right) \cdot \mathsf{fma}\left(im, im, 2\right)\\
\mathbf{else}:\\
\;\;\;\;\left(0.5 \cdot re\_m\right) \cdot \mathsf{fma}\left(im, im, 2\right)\\
\end{array}
\end{array}
if (*.f64 (*.f64 #s(literal 1/2 binary64) (sin.f64 re)) (+.f64 (exp.f64 (-.f64 #s(literal 0 binary64) im)) (exp.f64 im))) < 2.0000000000000001e-4Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
unpow2N/A
lower-fma.f6474.8
Applied rewrites74.8%
Taylor expanded in re around 0
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
pow2N/A
+-commutativeN/A
lift-fma.f64N/A
lift-*.f6455.4
Applied rewrites55.4%
if 2.0000000000000001e-4 < (*.f64 (*.f64 #s(literal 1/2 binary64) (sin.f64 re)) (+.f64 (exp.f64 (-.f64 #s(literal 0 binary64) im)) (exp.f64 im))) Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
unpow2N/A
lower-fma.f6471.4
Applied rewrites71.4%
Taylor expanded in re around 0
Applied rewrites33.0%
Final simplification46.7%
re\_m = (fabs.f64 re)
re\_s = (copysign.f64 #s(literal 1 binary64) re)
(FPCore (re_s re_m im)
:precision binary64
(*
re_s
(if (<= (* (* 0.5 (sin re_m)) (+ (exp (- im)) (exp im))) -0.05)
(* (* (fma (* re_m re_m) -0.08333333333333333 0.5) re_m) (* im im))
(* (* 0.5 re_m) (fma im im 2.0)))))re\_m = fabs(re);
re\_s = copysign(1.0, re);
double code(double re_s, double re_m, double im) {
double tmp;
if (((0.5 * sin(re_m)) * (exp(-im) + exp(im))) <= -0.05) {
tmp = (fma((re_m * re_m), -0.08333333333333333, 0.5) * re_m) * (im * im);
} else {
tmp = (0.5 * re_m) * fma(im, im, 2.0);
}
return re_s * tmp;
}
re\_m = abs(re) re\_s = copysign(1.0, re) function code(re_s, re_m, im) tmp = 0.0 if (Float64(Float64(0.5 * sin(re_m)) * Float64(exp(Float64(-im)) + exp(im))) <= -0.05) tmp = Float64(Float64(fma(Float64(re_m * re_m), -0.08333333333333333, 0.5) * re_m) * Float64(im * im)); else tmp = Float64(Float64(0.5 * re_m) * fma(im, im, 2.0)); end return Float64(re_s * tmp) end
re\_m = N[Abs[re], $MachinePrecision]
re\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[re]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[re$95$s_, re$95$m_, im_] := N[(re$95$s * If[LessEqual[N[(N[(0.5 * N[Sin[re$95$m], $MachinePrecision]), $MachinePrecision] * N[(N[Exp[(-im)], $MachinePrecision] + N[Exp[im], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], -0.05], N[(N[(N[(N[(re$95$m * re$95$m), $MachinePrecision] * -0.08333333333333333 + 0.5), $MachinePrecision] * re$95$m), $MachinePrecision] * N[(im * im), $MachinePrecision]), $MachinePrecision], N[(N[(0.5 * re$95$m), $MachinePrecision] * N[(im * im + 2.0), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
re\_m = \left|re\right|
\\
re\_s = \mathsf{copysign}\left(1, re\right)
\\
re\_s \cdot \begin{array}{l}
\mathbf{if}\;\left(0.5 \cdot \sin re\_m\right) \cdot \left(e^{-im} + e^{im}\right) \leq -0.05:\\
\;\;\;\;\left(\mathsf{fma}\left(re\_m \cdot re\_m, -0.08333333333333333, 0.5\right) \cdot re\_m\right) \cdot \left(im \cdot im\right)\\
\mathbf{else}:\\
\;\;\;\;\left(0.5 \cdot re\_m\right) \cdot \mathsf{fma}\left(im, im, 2\right)\\
\end{array}
\end{array}
if (*.f64 (*.f64 #s(literal 1/2 binary64) (sin.f64 re)) (+.f64 (exp.f64 (-.f64 #s(literal 0 binary64) im)) (exp.f64 im))) < -0.050000000000000003Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
unpow2N/A
lower-fma.f6463.3
Applied rewrites63.3%
Taylor expanded in re around 0
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
pow2N/A
+-commutativeN/A
lift-fma.f64N/A
lift-*.f6435.0
Applied rewrites35.0%
Taylor expanded in im around inf
pow2N/A
lift-*.f6435.0
Applied rewrites35.0%
if -0.050000000000000003 < (*.f64 (*.f64 #s(literal 1/2 binary64) (sin.f64 re)) (+.f64 (exp.f64 (-.f64 #s(literal 0 binary64) im)) (exp.f64 im))) Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
unpow2N/A
lower-fma.f6480.8
Applied rewrites80.8%
Taylor expanded in re around 0
Applied rewrites54.8%
Final simplification46.5%
re\_m = (fabs.f64 re)
re\_s = (copysign.f64 #s(literal 1 binary64) re)
(FPCore (re_s re_m im)
:precision binary64
(*
re_s
(if (<= (* (* 0.5 (sin re_m)) (+ (exp (- im)) (exp im))) 0.0002)
(* (fma -0.16666666666666666 (* re_m re_m) 1.0) re_m)
(* (* 0.5 re_m) (fma im im 2.0)))))re\_m = fabs(re);
re\_s = copysign(1.0, re);
double code(double re_s, double re_m, double im) {
double tmp;
if (((0.5 * sin(re_m)) * (exp(-im) + exp(im))) <= 0.0002) {
tmp = fma(-0.16666666666666666, (re_m * re_m), 1.0) * re_m;
} else {
tmp = (0.5 * re_m) * fma(im, im, 2.0);
}
return re_s * tmp;
}
re\_m = abs(re) re\_s = copysign(1.0, re) function code(re_s, re_m, im) tmp = 0.0 if (Float64(Float64(0.5 * sin(re_m)) * Float64(exp(Float64(-im)) + exp(im))) <= 0.0002) tmp = Float64(fma(-0.16666666666666666, Float64(re_m * re_m), 1.0) * re_m); else tmp = Float64(Float64(0.5 * re_m) * fma(im, im, 2.0)); end return Float64(re_s * tmp) end
re\_m = N[Abs[re], $MachinePrecision]
re\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[re]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[re$95$s_, re$95$m_, im_] := N[(re$95$s * If[LessEqual[N[(N[(0.5 * N[Sin[re$95$m], $MachinePrecision]), $MachinePrecision] * N[(N[Exp[(-im)], $MachinePrecision] + N[Exp[im], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 0.0002], N[(N[(-0.16666666666666666 * N[(re$95$m * re$95$m), $MachinePrecision] + 1.0), $MachinePrecision] * re$95$m), $MachinePrecision], N[(N[(0.5 * re$95$m), $MachinePrecision] * N[(im * im + 2.0), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
re\_m = \left|re\right|
\\
re\_s = \mathsf{copysign}\left(1, re\right)
\\
re\_s \cdot \begin{array}{l}
\mathbf{if}\;\left(0.5 \cdot \sin re\_m\right) \cdot \left(e^{-im} + e^{im}\right) \leq 0.0002:\\
\;\;\;\;\mathsf{fma}\left(-0.16666666666666666, re\_m \cdot re\_m, 1\right) \cdot re\_m\\
\mathbf{else}:\\
\;\;\;\;\left(0.5 \cdot re\_m\right) \cdot \mathsf{fma}\left(im, im, 2\right)\\
\end{array}
\end{array}
if (*.f64 (*.f64 #s(literal 1/2 binary64) (sin.f64 re)) (+.f64 (exp.f64 (-.f64 #s(literal 0 binary64) im)) (exp.f64 im))) < 2.0000000000000001e-4Initial program 100.0%
Taylor expanded in im around 0
lift-sin.f6451.8
Applied rewrites51.8%
Taylor expanded in re around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
lower-fma.f64N/A
pow2N/A
lift-*.f6439.1
Applied rewrites39.1%
if 2.0000000000000001e-4 < (*.f64 (*.f64 #s(literal 1/2 binary64) (sin.f64 re)) (+.f64 (exp.f64 (-.f64 #s(literal 0 binary64) im)) (exp.f64 im))) Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
unpow2N/A
lower-fma.f6471.4
Applied rewrites71.4%
Taylor expanded in re around 0
Applied rewrites33.0%
Final simplification36.7%
re\_m = (fabs.f64 re)
re\_s = (copysign.f64 #s(literal 1 binary64) re)
(FPCore (re_s re_m im)
:precision binary64
(*
re_s
(if (<= (* (* 0.5 (sin re_m)) (+ (exp (- im)) (exp im))) -0.002)
(* (* (* re_m re_m) -0.16666666666666666) re_m)
re_m)))re\_m = fabs(re);
re\_s = copysign(1.0, re);
double code(double re_s, double re_m, double im) {
double tmp;
if (((0.5 * sin(re_m)) * (exp(-im) + exp(im))) <= -0.002) {
tmp = ((re_m * re_m) * -0.16666666666666666) * re_m;
} else {
tmp = re_m;
}
return re_s * tmp;
}
re\_m = private
re\_s = private
module fmin_fmax_functions
implicit none
private
public fmax
public fmin
interface fmax
module procedure fmax88
module procedure fmax44
module procedure fmax84
module procedure fmax48
end interface
interface fmin
module procedure fmin88
module procedure fmin44
module procedure fmin84
module procedure fmin48
end interface
contains
real(8) function fmax88(x, y) result (res)
real(8), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(x, max(x, y), y /= y), x /= x)
end function
real(4) function fmax44(x, y) result (res)
real(4), intent (in) :: x
real(4), intent (in) :: y
res = merge(y, merge(x, max(x, y), y /= y), x /= x)
end function
real(8) function fmax84(x, y) result(res)
real(8), intent (in) :: x
real(4), intent (in) :: y
res = merge(dble(y), merge(x, max(x, dble(y)), y /= y), x /= x)
end function
real(8) function fmax48(x, y) result(res)
real(4), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(dble(x), max(dble(x), y), y /= y), x /= x)
end function
real(8) function fmin88(x, y) result (res)
real(8), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(x, min(x, y), y /= y), x /= x)
end function
real(4) function fmin44(x, y) result (res)
real(4), intent (in) :: x
real(4), intent (in) :: y
res = merge(y, merge(x, min(x, y), y /= y), x /= x)
end function
real(8) function fmin84(x, y) result(res)
real(8), intent (in) :: x
real(4), intent (in) :: y
res = merge(dble(y), merge(x, min(x, dble(y)), y /= y), x /= x)
end function
real(8) function fmin48(x, y) result(res)
real(4), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(dble(x), min(dble(x), y), y /= y), x /= x)
end function
end module
real(8) function code(re_s, re_m, im)
use fmin_fmax_functions
real(8), intent (in) :: re_s
real(8), intent (in) :: re_m
real(8), intent (in) :: im
real(8) :: tmp
if (((0.5d0 * sin(re_m)) * (exp(-im) + exp(im))) <= (-0.002d0)) then
tmp = ((re_m * re_m) * (-0.16666666666666666d0)) * re_m
else
tmp = re_m
end if
code = re_s * tmp
end function
re\_m = Math.abs(re);
re\_s = Math.copySign(1.0, re);
public static double code(double re_s, double re_m, double im) {
double tmp;
if (((0.5 * Math.sin(re_m)) * (Math.exp(-im) + Math.exp(im))) <= -0.002) {
tmp = ((re_m * re_m) * -0.16666666666666666) * re_m;
} else {
tmp = re_m;
}
return re_s * tmp;
}
re\_m = math.fabs(re) re\_s = math.copysign(1.0, re) def code(re_s, re_m, im): tmp = 0 if ((0.5 * math.sin(re_m)) * (math.exp(-im) + math.exp(im))) <= -0.002: tmp = ((re_m * re_m) * -0.16666666666666666) * re_m else: tmp = re_m return re_s * tmp
re\_m = abs(re) re\_s = copysign(1.0, re) function code(re_s, re_m, im) tmp = 0.0 if (Float64(Float64(0.5 * sin(re_m)) * Float64(exp(Float64(-im)) + exp(im))) <= -0.002) tmp = Float64(Float64(Float64(re_m * re_m) * -0.16666666666666666) * re_m); else tmp = re_m; end return Float64(re_s * tmp) end
re\_m = abs(re); re\_s = sign(re) * abs(1.0); function tmp_2 = code(re_s, re_m, im) tmp = 0.0; if (((0.5 * sin(re_m)) * (exp(-im) + exp(im))) <= -0.002) tmp = ((re_m * re_m) * -0.16666666666666666) * re_m; else tmp = re_m; end tmp_2 = re_s * tmp; end
re\_m = N[Abs[re], $MachinePrecision]
re\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[re]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[re$95$s_, re$95$m_, im_] := N[(re$95$s * If[LessEqual[N[(N[(0.5 * N[Sin[re$95$m], $MachinePrecision]), $MachinePrecision] * N[(N[Exp[(-im)], $MachinePrecision] + N[Exp[im], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], -0.002], N[(N[(N[(re$95$m * re$95$m), $MachinePrecision] * -0.16666666666666666), $MachinePrecision] * re$95$m), $MachinePrecision], re$95$m]), $MachinePrecision]
\begin{array}{l}
re\_m = \left|re\right|
\\
re\_s = \mathsf{copysign}\left(1, re\right)
\\
re\_s \cdot \begin{array}{l}
\mathbf{if}\;\left(0.5 \cdot \sin re\_m\right) \cdot \left(e^{-im} + e^{im}\right) \leq -0.002:\\
\;\;\;\;\left(\left(re\_m \cdot re\_m\right) \cdot -0.16666666666666666\right) \cdot re\_m\\
\mathbf{else}:\\
\;\;\;\;re\_m\\
\end{array}
\end{array}
if (*.f64 (*.f64 #s(literal 1/2 binary64) (sin.f64 re)) (+.f64 (exp.f64 (-.f64 #s(literal 0 binary64) im)) (exp.f64 im))) < -2e-3Initial program 100.0%
Taylor expanded in im around 0
lift-sin.f6429.8
Applied rewrites29.8%
Taylor expanded in re around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
lower-fma.f64N/A
pow2N/A
lift-*.f6411.2
Applied rewrites11.2%
Taylor expanded in re around inf
*-commutativeN/A
lower-*.f64N/A
pow2N/A
lift-*.f6410.9
Applied rewrites10.9%
if -2e-3 < (*.f64 (*.f64 #s(literal 1/2 binary64) (sin.f64 re)) (+.f64 (exp.f64 (-.f64 #s(literal 0 binary64) im)) (exp.f64 im))) Initial program 100.0%
Taylor expanded in re around 0
associate-*r*N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
cosh-undefN/A
lower-*.f64N/A
lower-cosh.f6469.4
Applied rewrites69.4%
Taylor expanded in im around 0
Applied rewrites34.4%
Final simplification24.6%
re\_m = (fabs.f64 re) re\_s = (copysign.f64 #s(literal 1 binary64) re) (FPCore (re_s re_m im) :precision binary64 (* re_s (* (* 0.5 (sin re_m)) (+ (exp (- im)) (exp im)))))
re\_m = fabs(re);
re\_s = copysign(1.0, re);
double code(double re_s, double re_m, double im) {
return re_s * ((0.5 * sin(re_m)) * (exp(-im) + exp(im)));
}
re\_m = private
re\_s = private
module fmin_fmax_functions
implicit none
private
public fmax
public fmin
interface fmax
module procedure fmax88
module procedure fmax44
module procedure fmax84
module procedure fmax48
end interface
interface fmin
module procedure fmin88
module procedure fmin44
module procedure fmin84
module procedure fmin48
end interface
contains
real(8) function fmax88(x, y) result (res)
real(8), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(x, max(x, y), y /= y), x /= x)
end function
real(4) function fmax44(x, y) result (res)
real(4), intent (in) :: x
real(4), intent (in) :: y
res = merge(y, merge(x, max(x, y), y /= y), x /= x)
end function
real(8) function fmax84(x, y) result(res)
real(8), intent (in) :: x
real(4), intent (in) :: y
res = merge(dble(y), merge(x, max(x, dble(y)), y /= y), x /= x)
end function
real(8) function fmax48(x, y) result(res)
real(4), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(dble(x), max(dble(x), y), y /= y), x /= x)
end function
real(8) function fmin88(x, y) result (res)
real(8), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(x, min(x, y), y /= y), x /= x)
end function
real(4) function fmin44(x, y) result (res)
real(4), intent (in) :: x
real(4), intent (in) :: y
res = merge(y, merge(x, min(x, y), y /= y), x /= x)
end function
real(8) function fmin84(x, y) result(res)
real(8), intent (in) :: x
real(4), intent (in) :: y
res = merge(dble(y), merge(x, min(x, dble(y)), y /= y), x /= x)
end function
real(8) function fmin48(x, y) result(res)
real(4), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(dble(x), min(dble(x), y), y /= y), x /= x)
end function
end module
real(8) function code(re_s, re_m, im)
use fmin_fmax_functions
real(8), intent (in) :: re_s
real(8), intent (in) :: re_m
real(8), intent (in) :: im
code = re_s * ((0.5d0 * sin(re_m)) * (exp(-im) + exp(im)))
end function
re\_m = Math.abs(re);
re\_s = Math.copySign(1.0, re);
public static double code(double re_s, double re_m, double im) {
return re_s * ((0.5 * Math.sin(re_m)) * (Math.exp(-im) + Math.exp(im)));
}
re\_m = math.fabs(re) re\_s = math.copysign(1.0, re) def code(re_s, re_m, im): return re_s * ((0.5 * math.sin(re_m)) * (math.exp(-im) + math.exp(im)))
re\_m = abs(re) re\_s = copysign(1.0, re) function code(re_s, re_m, im) return Float64(re_s * Float64(Float64(0.5 * sin(re_m)) * Float64(exp(Float64(-im)) + exp(im)))) end
re\_m = abs(re); re\_s = sign(re) * abs(1.0); function tmp = code(re_s, re_m, im) tmp = re_s * ((0.5 * sin(re_m)) * (exp(-im) + exp(im))); end
re\_m = N[Abs[re], $MachinePrecision]
re\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[re]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[re$95$s_, re$95$m_, im_] := N[(re$95$s * N[(N[(0.5 * N[Sin[re$95$m], $MachinePrecision]), $MachinePrecision] * N[(N[Exp[(-im)], $MachinePrecision] + N[Exp[im], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
re\_m = \left|re\right|
\\
re\_s = \mathsf{copysign}\left(1, re\right)
\\
re\_s \cdot \left(\left(0.5 \cdot \sin re\_m\right) \cdot \left(e^{-im} + e^{im}\right)\right)
\end{array}
Initial program 100.0%
Final simplification100.0%
re\_m = (fabs.f64 re)
re\_s = (copysign.f64 #s(literal 1 binary64) re)
(FPCore (re_s re_m im)
:precision binary64
(let* ((t_0 (* 0.5 (sin re_m))))
(*
re_s
(if (<= t_0 -0.001)
(* (* (fma (* re_m re_m) -0.08333333333333333 0.5) re_m) (fma im im 2.0))
(if (<= t_0 1e-8)
(* (* re_m 0.5) (fma (* (* im im) 0.08333333333333333) (* im im) 2.0))
(*
(*
(fma
(*
(- (* (* re_m re_m) 0.004166666666666667) 0.08333333333333333)
re_m)
re_m
0.5)
re_m)
(fma im im 2.0)))))))re\_m = fabs(re);
re\_s = copysign(1.0, re);
double code(double re_s, double re_m, double im) {
double t_0 = 0.5 * sin(re_m);
double tmp;
if (t_0 <= -0.001) {
tmp = (fma((re_m * re_m), -0.08333333333333333, 0.5) * re_m) * fma(im, im, 2.0);
} else if (t_0 <= 1e-8) {
tmp = (re_m * 0.5) * fma(((im * im) * 0.08333333333333333), (im * im), 2.0);
} else {
tmp = (fma(((((re_m * re_m) * 0.004166666666666667) - 0.08333333333333333) * re_m), re_m, 0.5) * re_m) * fma(im, im, 2.0);
}
return re_s * tmp;
}
re\_m = abs(re) re\_s = copysign(1.0, re) function code(re_s, re_m, im) t_0 = Float64(0.5 * sin(re_m)) tmp = 0.0 if (t_0 <= -0.001) tmp = Float64(Float64(fma(Float64(re_m * re_m), -0.08333333333333333, 0.5) * re_m) * fma(im, im, 2.0)); elseif (t_0 <= 1e-8) tmp = Float64(Float64(re_m * 0.5) * fma(Float64(Float64(im * im) * 0.08333333333333333), Float64(im * im), 2.0)); else tmp = Float64(Float64(fma(Float64(Float64(Float64(Float64(re_m * re_m) * 0.004166666666666667) - 0.08333333333333333) * re_m), re_m, 0.5) * re_m) * fma(im, im, 2.0)); end return Float64(re_s * tmp) end
re\_m = N[Abs[re], $MachinePrecision]
re\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[re]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[re$95$s_, re$95$m_, im_] := Block[{t$95$0 = N[(0.5 * N[Sin[re$95$m], $MachinePrecision]), $MachinePrecision]}, N[(re$95$s * If[LessEqual[t$95$0, -0.001], N[(N[(N[(N[(re$95$m * re$95$m), $MachinePrecision] * -0.08333333333333333 + 0.5), $MachinePrecision] * re$95$m), $MachinePrecision] * N[(im * im + 2.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, 1e-8], N[(N[(re$95$m * 0.5), $MachinePrecision] * N[(N[(N[(im * im), $MachinePrecision] * 0.08333333333333333), $MachinePrecision] * N[(im * im), $MachinePrecision] + 2.0), $MachinePrecision]), $MachinePrecision], N[(N[(N[(N[(N[(N[(N[(re$95$m * re$95$m), $MachinePrecision] * 0.004166666666666667), $MachinePrecision] - 0.08333333333333333), $MachinePrecision] * re$95$m), $MachinePrecision] * re$95$m + 0.5), $MachinePrecision] * re$95$m), $MachinePrecision] * N[(im * im + 2.0), $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]]
\begin{array}{l}
re\_m = \left|re\right|
\\
re\_s = \mathsf{copysign}\left(1, re\right)
\\
\begin{array}{l}
t_0 := 0.5 \cdot \sin re\_m\\
re\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_0 \leq -0.001:\\
\;\;\;\;\left(\mathsf{fma}\left(re\_m \cdot re\_m, -0.08333333333333333, 0.5\right) \cdot re\_m\right) \cdot \mathsf{fma}\left(im, im, 2\right)\\
\mathbf{elif}\;t\_0 \leq 10^{-8}:\\
\;\;\;\;\left(re\_m \cdot 0.5\right) \cdot \mathsf{fma}\left(\left(im \cdot im\right) \cdot 0.08333333333333333, im \cdot im, 2\right)\\
\mathbf{else}:\\
\;\;\;\;\left(\mathsf{fma}\left(\left(\left(re\_m \cdot re\_m\right) \cdot 0.004166666666666667 - 0.08333333333333333\right) \cdot re\_m, re\_m, 0.5\right) \cdot re\_m\right) \cdot \mathsf{fma}\left(im, im, 2\right)\\
\end{array}
\end{array}
\end{array}
if (*.f64 #s(literal 1/2 binary64) (sin.f64 re)) < -1e-3Initial program 99.9%
Taylor expanded in im around 0
+-commutativeN/A
unpow2N/A
lower-fma.f6469.2
Applied rewrites69.2%
Taylor expanded in re around 0
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
pow2N/A
+-commutativeN/A
lift-fma.f64N/A
lift-*.f6424.7
Applied rewrites24.7%
if -1e-3 < (*.f64 #s(literal 1/2 binary64) (sin.f64 re)) < 1e-8Initial program 100.0%
Taylor expanded in re around 0
associate-*r*N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
cosh-undefN/A
lower-*.f64N/A
lower-cosh.f64100.0
Applied rewrites100.0%
Taylor expanded in im around 0
cosh-undef-revN/A
sub0-negN/A
+-commutativeN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
pow2N/A
lower-*.f64N/A
pow2N/A
lower-*.f6487.0
Applied rewrites87.0%
Taylor expanded in im around inf
*-commutativeN/A
lower-*.f64N/A
pow2N/A
lift-*.f6487.0
Applied rewrites87.0%
if 1e-8 < (*.f64 #s(literal 1/2 binary64) (sin.f64 re)) Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
unpow2N/A
lower-fma.f6480.9
Applied rewrites80.9%
Taylor expanded in re around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f64N/A
lower-*.f64N/A
pow2N/A
lift-*.f64N/A
pow2N/A
lift-*.f6432.7
Applied rewrites32.7%
lift-fma.f64N/A
lift--.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
associate-*r*N/A
lower-fma.f64N/A
lower-*.f64N/A
pow2N/A
lower--.f64N/A
*-commutativeN/A
lower-*.f64N/A
pow2N/A
lift-*.f6432.7
Applied rewrites32.7%
re\_m = (fabs.f64 re)
re\_s = (copysign.f64 #s(literal 1 binary64) re)
(FPCore (re_s re_m im)
:precision binary64
(*
re_s
(if (<= (* 0.5 (sin re_m)) 0.0005)
(*
(*
(fma
(fma
(fma 0.002777777777777778 (* im im) 0.08333333333333333)
(* im im)
1.0)
(* im im)
2.0)
(fma (* re_m re_m) -0.08333333333333333 0.5))
re_m)
(*
(* (fma (* (* re_m re_m) 0.004166666666666667) (* re_m re_m) 0.5) re_m)
(fma im im 2.0)))))re\_m = fabs(re);
re\_s = copysign(1.0, re);
double code(double re_s, double re_m, double im) {
double tmp;
if ((0.5 * sin(re_m)) <= 0.0005) {
tmp = (fma(fma(fma(0.002777777777777778, (im * im), 0.08333333333333333), (im * im), 1.0), (im * im), 2.0) * fma((re_m * re_m), -0.08333333333333333, 0.5)) * re_m;
} else {
tmp = (fma(((re_m * re_m) * 0.004166666666666667), (re_m * re_m), 0.5) * re_m) * fma(im, im, 2.0);
}
return re_s * tmp;
}
re\_m = abs(re) re\_s = copysign(1.0, re) function code(re_s, re_m, im) tmp = 0.0 if (Float64(0.5 * sin(re_m)) <= 0.0005) tmp = Float64(Float64(fma(fma(fma(0.002777777777777778, Float64(im * im), 0.08333333333333333), Float64(im * im), 1.0), Float64(im * im), 2.0) * fma(Float64(re_m * re_m), -0.08333333333333333, 0.5)) * re_m); else tmp = Float64(Float64(fma(Float64(Float64(re_m * re_m) * 0.004166666666666667), Float64(re_m * re_m), 0.5) * re_m) * fma(im, im, 2.0)); end return Float64(re_s * tmp) end
re\_m = N[Abs[re], $MachinePrecision]
re\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[re]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[re$95$s_, re$95$m_, im_] := N[(re$95$s * If[LessEqual[N[(0.5 * N[Sin[re$95$m], $MachinePrecision]), $MachinePrecision], 0.0005], N[(N[(N[(N[(N[(0.002777777777777778 * N[(im * im), $MachinePrecision] + 0.08333333333333333), $MachinePrecision] * N[(im * im), $MachinePrecision] + 1.0), $MachinePrecision] * N[(im * im), $MachinePrecision] + 2.0), $MachinePrecision] * N[(N[(re$95$m * re$95$m), $MachinePrecision] * -0.08333333333333333 + 0.5), $MachinePrecision]), $MachinePrecision] * re$95$m), $MachinePrecision], N[(N[(N[(N[(N[(re$95$m * re$95$m), $MachinePrecision] * 0.004166666666666667), $MachinePrecision] * N[(re$95$m * re$95$m), $MachinePrecision] + 0.5), $MachinePrecision] * re$95$m), $MachinePrecision] * N[(im * im + 2.0), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
re\_m = \left|re\right|
\\
re\_s = \mathsf{copysign}\left(1, re\right)
\\
re\_s \cdot \begin{array}{l}
\mathbf{if}\;0.5 \cdot \sin re\_m \leq 0.0005:\\
\;\;\;\;\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(0.002777777777777778, im \cdot im, 0.08333333333333333\right), im \cdot im, 1\right), im \cdot im, 2\right) \cdot \mathsf{fma}\left(re\_m \cdot re\_m, -0.08333333333333333, 0.5\right)\right) \cdot re\_m\\
\mathbf{else}:\\
\;\;\;\;\left(\mathsf{fma}\left(\left(re\_m \cdot re\_m\right) \cdot 0.004166666666666667, re\_m \cdot re\_m, 0.5\right) \cdot re\_m\right) \cdot \mathsf{fma}\left(im, im, 2\right)\\
\end{array}
\end{array}
if (*.f64 #s(literal 1/2 binary64) (sin.f64 re)) < 5.0000000000000001e-4Initial program 100.0%
Taylor expanded in re around 0
*-commutativeN/A
lower-*.f64N/A
Applied rewrites73.9%
Taylor expanded in im around 0
cosh-undef-revN/A
sub0-negN/A
+-commutativeN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
pow2N/A
lower-*.f64N/A
pow2N/A
lower-*.f64N/A
pow2N/A
lower-*.f6467.3
Applied rewrites67.3%
if 5.0000000000000001e-4 < (*.f64 #s(literal 1/2 binary64) (sin.f64 re)) Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
unpow2N/A
lower-fma.f6480.3
Applied rewrites80.3%
Taylor expanded in re around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f64N/A
lower-*.f64N/A
pow2N/A
lift-*.f64N/A
pow2N/A
lift-*.f6430.7
Applied rewrites30.7%
Taylor expanded in re around inf
*-commutativeN/A
lower-*.f64N/A
pow2N/A
lift-*.f6430.7
Applied rewrites30.7%
re\_m = (fabs.f64 re)
re\_s = (copysign.f64 #s(literal 1 binary64) re)
(FPCore (re_s re_m im)
:precision binary64
(*
re_s
(if (<= (* 0.5 (sin re_m)) 0.0005)
(* (* (fma (* re_m re_m) -0.08333333333333333 0.5) re_m) (fma im im 2.0))
(*
(fma
(- (* 0.008333333333333333 (* re_m re_m)) 0.16666666666666666)
(* re_m re_m)
1.0)
re_m))))re\_m = fabs(re);
re\_s = copysign(1.0, re);
double code(double re_s, double re_m, double im) {
double tmp;
if ((0.5 * sin(re_m)) <= 0.0005) {
tmp = (fma((re_m * re_m), -0.08333333333333333, 0.5) * re_m) * fma(im, im, 2.0);
} else {
tmp = fma(((0.008333333333333333 * (re_m * re_m)) - 0.16666666666666666), (re_m * re_m), 1.0) * re_m;
}
return re_s * tmp;
}
re\_m = abs(re) re\_s = copysign(1.0, re) function code(re_s, re_m, im) tmp = 0.0 if (Float64(0.5 * sin(re_m)) <= 0.0005) tmp = Float64(Float64(fma(Float64(re_m * re_m), -0.08333333333333333, 0.5) * re_m) * fma(im, im, 2.0)); else tmp = Float64(fma(Float64(Float64(0.008333333333333333 * Float64(re_m * re_m)) - 0.16666666666666666), Float64(re_m * re_m), 1.0) * re_m); end return Float64(re_s * tmp) end
re\_m = N[Abs[re], $MachinePrecision]
re\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[re]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[re$95$s_, re$95$m_, im_] := N[(re$95$s * If[LessEqual[N[(0.5 * N[Sin[re$95$m], $MachinePrecision]), $MachinePrecision], 0.0005], N[(N[(N[(N[(re$95$m * re$95$m), $MachinePrecision] * -0.08333333333333333 + 0.5), $MachinePrecision] * re$95$m), $MachinePrecision] * N[(im * im + 2.0), $MachinePrecision]), $MachinePrecision], N[(N[(N[(N[(0.008333333333333333 * N[(re$95$m * re$95$m), $MachinePrecision]), $MachinePrecision] - 0.16666666666666666), $MachinePrecision] * N[(re$95$m * re$95$m), $MachinePrecision] + 1.0), $MachinePrecision] * re$95$m), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
re\_m = \left|re\right|
\\
re\_s = \mathsf{copysign}\left(1, re\right)
\\
re\_s \cdot \begin{array}{l}
\mathbf{if}\;0.5 \cdot \sin re\_m \leq 0.0005:\\
\;\;\;\;\left(\mathsf{fma}\left(re\_m \cdot re\_m, -0.08333333333333333, 0.5\right) \cdot re\_m\right) \cdot \mathsf{fma}\left(im, im, 2\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(0.008333333333333333 \cdot \left(re\_m \cdot re\_m\right) - 0.16666666666666666, re\_m \cdot re\_m, 1\right) \cdot re\_m\\
\end{array}
\end{array}
if (*.f64 #s(literal 1/2 binary64) (sin.f64 re)) < 5.0000000000000001e-4Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
unpow2N/A
lower-fma.f6471.0
Applied rewrites71.0%
Taylor expanded in re around 0
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
pow2N/A
+-commutativeN/A
lift-fma.f64N/A
lift-*.f6455.0
Applied rewrites55.0%
if 5.0000000000000001e-4 < (*.f64 #s(literal 1/2 binary64) (sin.f64 re)) Initial program 100.0%
Taylor expanded in im around 0
lift-sin.f6445.6
Applied rewrites45.6%
Taylor expanded in re around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f64N/A
lower-*.f64N/A
pow2N/A
lift-*.f64N/A
pow2N/A
lift-*.f6425.0
Applied rewrites25.0%
re\_m = (fabs.f64 re) re\_s = (copysign.f64 #s(literal 1 binary64) re) (FPCore (re_s re_m im) :precision binary64 (* re_s (* (fma -0.16666666666666666 (* re_m re_m) 1.0) re_m)))
re\_m = fabs(re);
re\_s = copysign(1.0, re);
double code(double re_s, double re_m, double im) {
return re_s * (fma(-0.16666666666666666, (re_m * re_m), 1.0) * re_m);
}
re\_m = abs(re) re\_s = copysign(1.0, re) function code(re_s, re_m, im) return Float64(re_s * Float64(fma(-0.16666666666666666, Float64(re_m * re_m), 1.0) * re_m)) end
re\_m = N[Abs[re], $MachinePrecision]
re\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[re]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[re$95$s_, re$95$m_, im_] := N[(re$95$s * N[(N[(-0.16666666666666666 * N[(re$95$m * re$95$m), $MachinePrecision] + 1.0), $MachinePrecision] * re$95$m), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
re\_m = \left|re\right|
\\
re\_s = \mathsf{copysign}\left(1, re\right)
\\
re\_s \cdot \left(\mathsf{fma}\left(-0.16666666666666666, re\_m \cdot re\_m, 1\right) \cdot re\_m\right)
\end{array}
Initial program 100.0%
Taylor expanded in im around 0
lift-sin.f6443.8
Applied rewrites43.8%
Taylor expanded in re around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
lower-fma.f64N/A
pow2N/A
lift-*.f6428.8
Applied rewrites28.8%
re\_m = (fabs.f64 re) re\_s = (copysign.f64 #s(literal 1 binary64) re) (FPCore (re_s re_m im) :precision binary64 (* re_s re_m))
re\_m = fabs(re);
re\_s = copysign(1.0, re);
double code(double re_s, double re_m, double im) {
return re_s * re_m;
}
re\_m = private
re\_s = private
module fmin_fmax_functions
implicit none
private
public fmax
public fmin
interface fmax
module procedure fmax88
module procedure fmax44
module procedure fmax84
module procedure fmax48
end interface
interface fmin
module procedure fmin88
module procedure fmin44
module procedure fmin84
module procedure fmin48
end interface
contains
real(8) function fmax88(x, y) result (res)
real(8), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(x, max(x, y), y /= y), x /= x)
end function
real(4) function fmax44(x, y) result (res)
real(4), intent (in) :: x
real(4), intent (in) :: y
res = merge(y, merge(x, max(x, y), y /= y), x /= x)
end function
real(8) function fmax84(x, y) result(res)
real(8), intent (in) :: x
real(4), intent (in) :: y
res = merge(dble(y), merge(x, max(x, dble(y)), y /= y), x /= x)
end function
real(8) function fmax48(x, y) result(res)
real(4), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(dble(x), max(dble(x), y), y /= y), x /= x)
end function
real(8) function fmin88(x, y) result (res)
real(8), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(x, min(x, y), y /= y), x /= x)
end function
real(4) function fmin44(x, y) result (res)
real(4), intent (in) :: x
real(4), intent (in) :: y
res = merge(y, merge(x, min(x, y), y /= y), x /= x)
end function
real(8) function fmin84(x, y) result(res)
real(8), intent (in) :: x
real(4), intent (in) :: y
res = merge(dble(y), merge(x, min(x, dble(y)), y /= y), x /= x)
end function
real(8) function fmin48(x, y) result(res)
real(4), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(dble(x), min(dble(x), y), y /= y), x /= x)
end function
end module
real(8) function code(re_s, re_m, im)
use fmin_fmax_functions
real(8), intent (in) :: re_s
real(8), intent (in) :: re_m
real(8), intent (in) :: im
code = re_s * re_m
end function
re\_m = Math.abs(re);
re\_s = Math.copySign(1.0, re);
public static double code(double re_s, double re_m, double im) {
return re_s * re_m;
}
re\_m = math.fabs(re) re\_s = math.copysign(1.0, re) def code(re_s, re_m, im): return re_s * re_m
re\_m = abs(re) re\_s = copysign(1.0, re) function code(re_s, re_m, im) return Float64(re_s * re_m) end
re\_m = abs(re); re\_s = sign(re) * abs(1.0); function tmp = code(re_s, re_m, im) tmp = re_s * re_m; end
re\_m = N[Abs[re], $MachinePrecision]
re\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[re]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[re$95$s_, re$95$m_, im_] := N[(re$95$s * re$95$m), $MachinePrecision]
\begin{array}{l}
re\_m = \left|re\right|
\\
re\_s = \mathsf{copysign}\left(1, re\right)
\\
re\_s \cdot re\_m
\end{array}
Initial program 100.0%
Taylor expanded in re around 0
associate-*r*N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
cosh-undefN/A
lower-*.f64N/A
lower-cosh.f6463.8
Applied rewrites63.8%
Taylor expanded in im around 0
Applied rewrites21.1%
herbie shell --seed 2025064
(FPCore (re im)
:name "math.sin on complex, real part"
:precision binary64
(* (* 0.5 (sin re)) (+ (exp (- 0.0 im)) (exp im))))