
(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 17 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)) (+ (exp (- im)) (exp im)))))
(*
re_s
(if (<= t_1 (- INFINITY))
(* (* t_0 (* (* re_m re_m) -0.08333333333333333)) re_m)
(if (<= t_1 1.0)
(fma
(* (sin re_m) (fma 0.041666666666666664 (* im im) 0.5))
(* im im)
(sin re_m))
(* (* 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)) * (exp(-im) + exp(im));
double tmp;
if (t_1 <= -((double) INFINITY)) {
tmp = (t_0 * ((re_m * re_m) * -0.08333333333333333)) * re_m;
} else if (t_1 <= 1.0) {
tmp = fma((sin(re_m) * fma(0.041666666666666664, (im * im), 0.5)), (im * im), sin(re_m));
} 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(Float64(0.5 * sin(re_m)) * Float64(exp(Float64(-im)) + exp(im))) tmp = 0.0 if (t_1 <= Float64(-Inf)) tmp = Float64(Float64(t_0 * Float64(Float64(re_m * re_m) * -0.08333333333333333)) * re_m); elseif (t_1 <= 1.0) tmp = fma(Float64(sin(re_m) * fma(0.041666666666666664, Float64(im * im), 0.5)), Float64(im * im), sin(re_m)); 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[(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$1, (-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$1, 1.0], N[(N[(N[Sin[re$95$m], $MachinePrecision] * N[(0.041666666666666664 * N[(im * im), $MachinePrecision] + 0.5), $MachinePrecision]), $MachinePrecision] * N[(im * im), $MachinePrecision] + N[Sin[re$95$m], $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 := \left(0.5 \cdot \sin re\_m\right) \cdot \left(e^{-im} + e^{im}\right)\\
re\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_1 \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\_1 \leq 1:\\
\;\;\;\;\mathsf{fma}\left(\sin re\_m \cdot \mathsf{fma}\left(0.041666666666666664, im \cdot im, 0.5\right), im \cdot im, \sin re\_m\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 rewrites77.6%
Taylor expanded in re around inf
*-commutativeN/A
lower-*.f64N/A
pow2N/A
lift-*.f6414.9
Applied rewrites14.9%
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
*-commutativeN/A
lower-fma.f64N/A
associate-*r*N/A
distribute-rgt-outN/A
lower-*.f64N/A
lift-sin.f64N/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
lift-sin.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.f6473.1
Applied rewrites73.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 (* 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 (fma (* im im) 0.08333333333333333 1.0) (* 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(fma((im * im), 0.08333333333333333, 1.0), (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(fma(Float64(im * im), 0.08333333333333333, 1.0), Float64(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[(N[(N[(im * im), $MachinePrecision] * 0.08333333333333333 + 1.0), $MachinePrecision] * N[(im * im), $MachinePrecision] + 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(\mathsf{fma}\left(im \cdot im, 0.08333333333333333, 1\right), im \cdot 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 rewrites77.6%
Taylor expanded in re around inf
*-commutativeN/A
lower-*.f64N/A
pow2N/A
lift-*.f6414.9
Applied rewrites14.9%
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
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.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.f6473.1
Applied rewrites73.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 (* 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 rewrites77.6%
Taylor expanded in re around inf
*-commutativeN/A
lower-*.f64N/A
pow2N/A
lift-*.f6414.9
Applied rewrites14.9%
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.f6499.9
Applied rewrites99.9%
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.f6473.1
Applied rewrites73.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 -9.92063492063492e-5 (* re_m re_m) 0.004166666666666667)
(* re_m re_m))
0.08333333333333333)
(* re_m re_m)
0.5)
re_m)
(fma (fma (* im im) 0.08333333333333333 1.0) (* im im) 2.0))
(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(-9.92063492063492e-5, (re_m * re_m), 0.004166666666666667) * (re_m * re_m)) - 0.08333333333333333), (re_m * re_m), 0.5) * re_m) * fma(fma((im * im), 0.08333333333333333, 1.0), (im * im), 2.0);
} 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(Float64(Float64(fma(-9.92063492063492e-5, Float64(re_m * re_m), 0.004166666666666667) * Float64(re_m * re_m)) - 0.08333333333333333), Float64(re_m * re_m), 0.5) * re_m) * fma(fma(Float64(im * im), 0.08333333333333333, 1.0), Float64(im * im), 2.0)); 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[(N[(-9.92063492063492e-5 * N[(re$95$m * re$95$m), $MachinePrecision] + 0.004166666666666667), $MachinePrecision] * N[(re$95$m * re$95$m), $MachinePrecision]), $MachinePrecision] - 0.08333333333333333), $MachinePrecision] * N[(re$95$m * re$95$m), $MachinePrecision] + 0.5), $MachinePrecision] * re$95$m), $MachinePrecision] * N[(N[(N[(im * im), $MachinePrecision] * 0.08333333333333333 + 1.0), $MachinePrecision] * N[(im * im), $MachinePrecision] + 2.0), $MachinePrecision]), $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(-9.92063492063492 \cdot 10^{-5}, re\_m \cdot re\_m, 0.004166666666666667\right) \cdot \left(re\_m \cdot re\_m\right) - 0.08333333333333333, re\_m \cdot re\_m, 0.5\right) \cdot re\_m\right) \cdot \mathsf{fma}\left(\mathsf{fma}\left(im \cdot im, 0.08333333333333333, 1\right), im \cdot im, 2\right)\\
\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 im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6463.9
Applied rewrites63.9%
Taylor expanded in re around 0
*-commutativeN/A
lower-*.f64N/A
Applied rewrites51.6%
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.f6499.9
Applied rewrites99.9%
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.f6473.1
Applied rewrites73.1%
Final simplification80.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 -9.92063492063492e-5 (* re_m re_m) 0.004166666666666667)
(* re_m re_m))
0.08333333333333333)
(* re_m re_m)
0.5)
re_m)
(fma (fma (* im im) 0.08333333333333333 1.0) (* im im) 2.0))
(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(-9.92063492063492e-5, (re_m * re_m), 0.004166666666666667) * (re_m * re_m)) - 0.08333333333333333), (re_m * re_m), 0.5) * re_m) * fma(fma((im * im), 0.08333333333333333, 1.0), (im * im), 2.0);
} 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(Float64(Float64(fma(-9.92063492063492e-5, Float64(re_m * re_m), 0.004166666666666667) * Float64(re_m * re_m)) - 0.08333333333333333), Float64(re_m * re_m), 0.5) * re_m) * fma(fma(Float64(im * im), 0.08333333333333333, 1.0), Float64(im * im), 2.0)); 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[(N[(-9.92063492063492e-5 * N[(re$95$m * re$95$m), $MachinePrecision] + 0.004166666666666667), $MachinePrecision] * N[(re$95$m * re$95$m), $MachinePrecision]), $MachinePrecision] - 0.08333333333333333), $MachinePrecision] * N[(re$95$m * re$95$m), $MachinePrecision] + 0.5), $MachinePrecision] * re$95$m), $MachinePrecision] * N[(N[(N[(im * im), $MachinePrecision] * 0.08333333333333333 + 1.0), $MachinePrecision] * N[(im * im), $MachinePrecision] + 2.0), $MachinePrecision]), $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(-9.92063492063492 \cdot 10^{-5}, re\_m \cdot re\_m, 0.004166666666666667\right) \cdot \left(re\_m \cdot re\_m\right) - 0.08333333333333333, re\_m \cdot re\_m, 0.5\right) \cdot re\_m\right) \cdot \mathsf{fma}\left(\mathsf{fma}\left(im \cdot im, 0.08333333333333333, 1\right), im \cdot im, 2\right)\\
\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 im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6463.9
Applied rewrites63.9%
Taylor expanded in re around 0
*-commutativeN/A
lower-*.f64N/A
Applied rewrites51.6%
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.2
Applied rewrites99.2%
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.f6473.1
Applied rewrites73.1%
Final simplification79.9%
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 -9.92063492063492e-5 (* re_m re_m) 0.004166666666666667)
(* re_m re_m))
0.08333333333333333)
(* re_m re_m)
0.5)
re_m)
(fma (fma (* im im) 0.08333333333333333 1.0) (* im im) 2.0))
(if (<= t_0 1.0)
(sin re_m)
(*
(* re_m 0.5)
(fma
(*
(fma
(fma (* im im) 0.002777777777777778 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(-9.92063492063492e-5, (re_m * re_m), 0.004166666666666667) * (re_m * re_m)) - 0.08333333333333333), (re_m * re_m), 0.5) * re_m) * fma(fma((im * im), 0.08333333333333333, 1.0), (im * im), 2.0);
} else if (t_0 <= 1.0) {
tmp = sin(re_m);
} else {
tmp = (re_m * 0.5) * fma((fma(fma((im * im), 0.002777777777777778, 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(Float64(Float64(fma(-9.92063492063492e-5, Float64(re_m * re_m), 0.004166666666666667) * Float64(re_m * re_m)) - 0.08333333333333333), Float64(re_m * re_m), 0.5) * re_m) * fma(fma(Float64(im * im), 0.08333333333333333, 1.0), Float64(im * im), 2.0)); elseif (t_0 <= 1.0) tmp = sin(re_m); else tmp = Float64(Float64(re_m * 0.5) * fma(Float64(fma(fma(Float64(im * im), 0.002777777777777778, 0.08333333333333333), 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_] := 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[(N[(-9.92063492063492e-5 * N[(re$95$m * re$95$m), $MachinePrecision] + 0.004166666666666667), $MachinePrecision] * N[(re$95$m * re$95$m), $MachinePrecision]), $MachinePrecision] - 0.08333333333333333), $MachinePrecision] * N[(re$95$m * re$95$m), $MachinePrecision] + 0.5), $MachinePrecision] * re$95$m), $MachinePrecision] * N[(N[(N[(im * im), $MachinePrecision] * 0.08333333333333333 + 1.0), $MachinePrecision] * N[(im * im), $MachinePrecision] + 2.0), $MachinePrecision]), $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[(im * im), $MachinePrecision] * 0.002777777777777778 + 0.08333333333333333), $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)
\\
\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(-9.92063492063492 \cdot 10^{-5}, re\_m \cdot re\_m, 0.004166666666666667\right) \cdot \left(re\_m \cdot re\_m\right) - 0.08333333333333333, re\_m \cdot re\_m, 0.5\right) \cdot re\_m\right) \cdot \mathsf{fma}\left(\mathsf{fma}\left(im \cdot im, 0.08333333333333333, 1\right), im \cdot im, 2\right)\\
\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(im \cdot im, 0.002777777777777778, 0.08333333333333333\right), im \cdot 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 im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6463.9
Applied rewrites63.9%
Taylor expanded in re around 0
*-commutativeN/A
lower-*.f64N/A
Applied rewrites51.6%
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.2
Applied rewrites99.2%
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.f6473.1
Applied rewrites73.1%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
pow2N/A
lift-*.f64N/A
pow2N/A
lift-*.f64N/A
pow2N/A
lift-*.f6464.5
Applied rewrites64.5%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
pow2N/A
lift-*.f64N/A
pow2N/A
lift-*.f64N/A
pow2N/A
lift-*.f6464.5
Applied rewrites64.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
lower-*.f64N/A
lower-fma.f64N/A
*-commutativeN/A
lower-fma.f64N/A
lift-*.f64N/A
lift-*.f6464.5
Applied rewrites64.5%
Final simplification77.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))) 4e-5)
(*
(* (fma (* re_m re_m) -0.08333333333333333 0.5) re_m)
(fma (fma (* im im) 0.08333333333333333 1.0) (* im im) 2.0))
(*
(* re_m 0.5)
(fma
(*
(fma
(fma (* im im) 0.002777777777777778 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))) <= 4e-5) {
tmp = (fma((re_m * re_m), -0.08333333333333333, 0.5) * re_m) * fma(fma((im * im), 0.08333333333333333, 1.0), (im * im), 2.0);
} else {
tmp = (re_m * 0.5) * fma((fma(fma((im * im), 0.002777777777777778, 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))) <= 4e-5) tmp = Float64(Float64(fma(Float64(re_m * re_m), -0.08333333333333333, 0.5) * re_m) * fma(fma(Float64(im * im), 0.08333333333333333, 1.0), Float64(im * im), 2.0)); else tmp = Float64(Float64(re_m * 0.5) * fma(Float64(fma(fma(Float64(im * im), 0.002777777777777778, 0.08333333333333333), 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], 4e-5], N[(N[(N[(N[(re$95$m * re$95$m), $MachinePrecision] * -0.08333333333333333 + 0.5), $MachinePrecision] * re$95$m), $MachinePrecision] * N[(N[(N[(im * im), $MachinePrecision] * 0.08333333333333333 + 1.0), $MachinePrecision] * N[(im * im), $MachinePrecision] + 2.0), $MachinePrecision]), $MachinePrecision], N[(N[(re$95$m * 0.5), $MachinePrecision] * N[(N[(N[(N[(N[(im * im), $MachinePrecision] * 0.002777777777777778 + 0.08333333333333333), $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 4 \cdot 10^{-5}:\\
\;\;\;\;\left(\mathsf{fma}\left(re\_m \cdot re\_m, -0.08333333333333333, 0.5\right) \cdot re\_m\right) \cdot \mathsf{fma}\left(\mathsf{fma}\left(im \cdot im, 0.08333333333333333, 1\right), im \cdot im, 2\right)\\
\mathbf{else}:\\
\;\;\;\;\left(re\_m \cdot 0.5\right) \cdot \mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(im \cdot im, 0.002777777777777778, 0.08333333333333333\right), 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))) < 4.00000000000000033e-5Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6484.9
Applied rewrites84.9%
Taylor expanded in re around 0
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
pow2N/A
+-commutativeN/A
lift-fma.f64N/A
lift-*.f6459.4
Applied rewrites59.4%
if 4.00000000000000033e-5 < (*.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.f6452.0
Applied rewrites52.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
pow2N/A
lift-*.f64N/A
pow2N/A
lift-*.f64N/A
pow2N/A
lift-*.f6446.0
Applied rewrites46.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
pow2N/A
lift-*.f64N/A
pow2N/A
lift-*.f64N/A
pow2N/A
lift-*.f6446.0
Applied rewrites46.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
lower-*.f64N/A
lower-fma.f64N/A
*-commutativeN/A
lower-fma.f64N/A
lift-*.f64N/A
lift-*.f6446.0
Applied rewrites46.0%
Final simplification54.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))) 4e-5)
(*
(* (fma (* re_m re_m) -0.08333333333333333 0.5) re_m)
(fma (* (* im im) 0.08333333333333333) (* im im) 2.0))
(*
(* re_m 0.5)
(fma
(*
(fma
(fma (* im im) 0.002777777777777778 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))) <= 4e-5) {
tmp = (fma((re_m * re_m), -0.08333333333333333, 0.5) * re_m) * fma(((im * im) * 0.08333333333333333), (im * im), 2.0);
} else {
tmp = (re_m * 0.5) * fma((fma(fma((im * im), 0.002777777777777778, 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))) <= 4e-5) tmp = Float64(Float64(fma(Float64(re_m * re_m), -0.08333333333333333, 0.5) * re_m) * fma(Float64(Float64(im * im) * 0.08333333333333333), Float64(im * im), 2.0)); else tmp = Float64(Float64(re_m * 0.5) * fma(Float64(fma(fma(Float64(im * im), 0.002777777777777778, 0.08333333333333333), 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], 4e-5], N[(N[(N[(N[(re$95$m * re$95$m), $MachinePrecision] * -0.08333333333333333 + 0.5), $MachinePrecision] * re$95$m), $MachinePrecision] * N[(N[(N[(im * im), $MachinePrecision] * 0.08333333333333333), $MachinePrecision] * N[(im * im), $MachinePrecision] + 2.0), $MachinePrecision]), $MachinePrecision], N[(N[(re$95$m * 0.5), $MachinePrecision] * N[(N[(N[(N[(N[(im * im), $MachinePrecision] * 0.002777777777777778 + 0.08333333333333333), $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 4 \cdot 10^{-5}:\\
\;\;\;\;\left(\mathsf{fma}\left(re\_m \cdot re\_m, -0.08333333333333333, 0.5\right) \cdot re\_m\right) \cdot \mathsf{fma}\left(\left(im \cdot im\right) \cdot 0.08333333333333333, im \cdot im, 2\right)\\
\mathbf{else}:\\
\;\;\;\;\left(re\_m \cdot 0.5\right) \cdot \mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(im \cdot im, 0.002777777777777778, 0.08333333333333333\right), 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))) < 4.00000000000000033e-5Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6484.9
Applied rewrites84.9%
Taylor expanded in re around 0
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
pow2N/A
+-commutativeN/A
lift-fma.f64N/A
lift-*.f6459.4
Applied rewrites59.4%
Taylor expanded in im around inf
pow2N/A
*-commutativeN/A
lower-*.f64N/A
lift-*.f6459.2
Applied rewrites59.2%
if 4.00000000000000033e-5 < (*.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.f6452.0
Applied rewrites52.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
pow2N/A
lift-*.f64N/A
pow2N/A
lift-*.f64N/A
pow2N/A
lift-*.f6446.0
Applied rewrites46.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
pow2N/A
lift-*.f64N/A
pow2N/A
lift-*.f64N/A
pow2N/A
lift-*.f6446.0
Applied rewrites46.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
lower-*.f64N/A
lower-fma.f64N/A
*-commutativeN/A
lower-fma.f64N/A
lift-*.f64N/A
lift-*.f6446.0
Applied rewrites46.0%
Final simplification54.2%
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))) 4e-5)
(*
(* (fma (* re_m re_m) -0.08333333333333333 0.5) re_m)
(fma (* (* im im) 0.08333333333333333) (* im im) 2.0))
(*
(* re_m 0.5)
(fma
(fma (* (* im im) 0.002777777777777778) (* 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))) <= 4e-5) {
tmp = (fma((re_m * re_m), -0.08333333333333333, 0.5) * re_m) * fma(((im * im) * 0.08333333333333333), (im * im), 2.0);
} else {
tmp = (re_m * 0.5) * fma(fma(((im * im) * 0.002777777777777778), (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))) <= 4e-5) tmp = Float64(Float64(fma(Float64(re_m * re_m), -0.08333333333333333, 0.5) * re_m) * fma(Float64(Float64(im * im) * 0.08333333333333333), Float64(im * im), 2.0)); else tmp = Float64(Float64(re_m * 0.5) * fma(fma(Float64(Float64(im * im) * 0.002777777777777778), Float64(im * im), 1.0), 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], 4e-5], N[(N[(N[(N[(re$95$m * re$95$m), $MachinePrecision] * -0.08333333333333333 + 0.5), $MachinePrecision] * re$95$m), $MachinePrecision] * N[(N[(N[(im * im), $MachinePrecision] * 0.08333333333333333), $MachinePrecision] * N[(im * im), $MachinePrecision] + 2.0), $MachinePrecision]), $MachinePrecision], N[(N[(re$95$m * 0.5), $MachinePrecision] * N[(N[(N[(N[(im * im), $MachinePrecision] * 0.002777777777777778), $MachinePrecision] * N[(im * im), $MachinePrecision] + 1.0), $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 4 \cdot 10^{-5}:\\
\;\;\;\;\left(\mathsf{fma}\left(re\_m \cdot re\_m, -0.08333333333333333, 0.5\right) \cdot re\_m\right) \cdot \mathsf{fma}\left(\left(im \cdot im\right) \cdot 0.08333333333333333, im \cdot im, 2\right)\\
\mathbf{else}:\\
\;\;\;\;\left(re\_m \cdot 0.5\right) \cdot \mathsf{fma}\left(\mathsf{fma}\left(\left(im \cdot im\right) \cdot 0.002777777777777778, im \cdot im, 1\right), 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))) < 4.00000000000000033e-5Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6484.9
Applied rewrites84.9%
Taylor expanded in re around 0
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
pow2N/A
+-commutativeN/A
lift-fma.f64N/A
lift-*.f6459.4
Applied rewrites59.4%
Taylor expanded in im around inf
pow2N/A
*-commutativeN/A
lower-*.f64N/A
lift-*.f6459.2
Applied rewrites59.2%
if 4.00000000000000033e-5 < (*.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.f6452.0
Applied rewrites52.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
pow2N/A
lift-*.f64N/A
pow2N/A
lift-*.f64N/A
pow2N/A
lift-*.f6446.0
Applied rewrites46.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
pow2N/A
lift-*.f64N/A
pow2N/A
lift-*.f64N/A
pow2N/A
lift-*.f6446.0
Applied rewrites46.0%
Taylor expanded in im around inf
pow2N/A
*-commutativeN/A
lower-*.f64N/A
lift-*.f6446.0
Applied rewrites46.0%
Final simplification54.2%
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))) 4e-5)
(* (* (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.002777777777777778) (* 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))) <= 4e-5) {
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.002777777777777778), (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))) <= 4e-5) 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(fma(Float64(Float64(im * im) * 0.002777777777777778), Float64(im * im), 1.0), 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], 4e-5], 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.002777777777777778), $MachinePrecision] * N[(im * im), $MachinePrecision] + 1.0), $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 4 \cdot 10^{-5}:\\
\;\;\;\;\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(im \cdot im\right) \cdot 0.002777777777777778, im \cdot im, 1\right), 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))) < 4.00000000000000033e-5Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
unpow2N/A
lower-fma.f6476.4
Applied rewrites76.4%
Taylor expanded in re around 0
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
pow2N/A
+-commutativeN/A
lift-fma.f64N/A
lift-*.f6454.6
Applied rewrites54.6%
if 4.00000000000000033e-5 < (*.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.f6452.0
Applied rewrites52.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
pow2N/A
lift-*.f64N/A
pow2N/A
lift-*.f64N/A
pow2N/A
lift-*.f6446.0
Applied rewrites46.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
pow2N/A
lift-*.f64N/A
pow2N/A
lift-*.f64N/A
pow2N/A
lift-*.f6446.0
Applied rewrites46.0%
Taylor expanded in im around inf
pow2N/A
*-commutativeN/A
lower-*.f64N/A
lift-*.f6446.0
Applied rewrites46.0%
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))) 4e-5)
(* (* (fma (* re_m re_m) -0.08333333333333333 0.5) re_m) (fma im im 2.0))
(* (fma (fma (* im im) 0.041666666666666664 0.5) (* im im) 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)) * (exp(-im) + exp(im))) <= 4e-5) {
tmp = (fma((re_m * re_m), -0.08333333333333333, 0.5) * re_m) * fma(im, im, 2.0);
} else {
tmp = fma(fma((im * im), 0.041666666666666664, 0.5), (im * im), 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(Float64(0.5 * sin(re_m)) * Float64(exp(Float64(-im)) + exp(im))) <= 4e-5) tmp = Float64(Float64(fma(Float64(re_m * re_m), -0.08333333333333333, 0.5) * re_m) * fma(im, im, 2.0)); else tmp = Float64(fma(fma(Float64(im * im), 0.041666666666666664, 0.5), Float64(im * im), 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[(N[(0.5 * N[Sin[re$95$m], $MachinePrecision]), $MachinePrecision] * N[(N[Exp[(-im)], $MachinePrecision] + N[Exp[im], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 4e-5], 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[(im * im), $MachinePrecision] * 0.041666666666666664 + 0.5), $MachinePrecision] * N[(im * im), $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}\;\left(0.5 \cdot \sin re\_m\right) \cdot \left(e^{-im} + e^{im}\right) \leq 4 \cdot 10^{-5}:\\
\;\;\;\;\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(\mathsf{fma}\left(im \cdot im, 0.041666666666666664, 0.5\right), im \cdot im, 1\right) \cdot 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))) < 4.00000000000000033e-5Initial program 100.0%
Taylor expanded in im around 0
+-commutativeN/A
unpow2N/A
lower-fma.f6476.4
Applied rewrites76.4%
Taylor expanded in re around 0
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
pow2N/A
+-commutativeN/A
lift-fma.f64N/A
lift-*.f6454.6
Applied rewrites54.6%
if 4.00000000000000033e-5 < (*.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.f6452.0
Applied rewrites52.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-fma.f64N/A
lower-*.f64N/A
pow2N/A
lift-*.f64N/A
*-commutativeN/A
lift-*.f64N/A
pow2N/A
lift-*.f6440.1
Applied rewrites40.1%
Taylor expanded in re around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
pow2N/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
pow2N/A
*-commutativeN/A
lower-fma.f64N/A
lift-*.f64N/A
lift-*.f6444.0
Applied rewrites44.0%
Final simplification50.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
(if (<= (* (* 0.5 (sin re_m)) (+ (exp (- im)) (exp im))) 4e-5)
(* (fma -0.16666666666666666 (* re_m re_m) 1.0) re_m)
(* (fma (fma (* im im) 0.041666666666666664 0.5) (* im im) 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)) * (exp(-im) + exp(im))) <= 4e-5) {
tmp = fma(-0.16666666666666666, (re_m * re_m), 1.0) * re_m;
} else {
tmp = fma(fma((im * im), 0.041666666666666664, 0.5), (im * im), 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(Float64(0.5 * sin(re_m)) * Float64(exp(Float64(-im)) + exp(im))) <= 4e-5) tmp = Float64(fma(-0.16666666666666666, Float64(re_m * re_m), 1.0) * re_m); else tmp = Float64(fma(fma(Float64(im * im), 0.041666666666666664, 0.5), Float64(im * im), 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[(N[(0.5 * N[Sin[re$95$m], $MachinePrecision]), $MachinePrecision] * N[(N[Exp[(-im)], $MachinePrecision] + N[Exp[im], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 4e-5], N[(N[(-0.16666666666666666 * N[(re$95$m * re$95$m), $MachinePrecision] + 1.0), $MachinePrecision] * re$95$m), $MachinePrecision], N[(N[(N[(N[(im * im), $MachinePrecision] * 0.041666666666666664 + 0.5), $MachinePrecision] * N[(im * im), $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}\;\left(0.5 \cdot \sin re\_m\right) \cdot \left(e^{-im} + e^{im}\right) \leq 4 \cdot 10^{-5}:\\
\;\;\;\;\mathsf{fma}\left(-0.16666666666666666, re\_m \cdot re\_m, 1\right) \cdot re\_m\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(im \cdot im, 0.041666666666666664, 0.5\right), im \cdot im, 1\right) \cdot 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))) < 4.00000000000000033e-5Initial program 100.0%
Taylor expanded in im around 0
lift-sin.f6458.6
Applied rewrites58.6%
Taylor expanded in re around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
lower-fma.f64N/A
pow2N/A
lift-*.f6442.1
Applied rewrites42.1%
if 4.00000000000000033e-5 < (*.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.f6452.0
Applied rewrites52.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-fma.f64N/A
lower-*.f64N/A
pow2N/A
lift-*.f64N/A
*-commutativeN/A
lift-*.f64N/A
pow2N/A
lift-*.f6440.1
Applied rewrites40.1%
Taylor expanded in re around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
pow2N/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
pow2N/A
*-commutativeN/A
lower-fma.f64N/A
lift-*.f64N/A
lift-*.f6444.0
Applied rewrites44.0%
Final simplification42.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
(if (<= (* (* 0.5 (sin re_m)) (+ (exp (- im)) (exp im))) 4e-5)
(* (fma -0.16666666666666666 (* re_m re_m) 1.0) re_m)
(fma (* (* (* im im) re_m) 0.041666666666666664) (* im im) 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))) <= 4e-5) {
tmp = fma(-0.16666666666666666, (re_m * re_m), 1.0) * re_m;
} else {
tmp = fma((((im * im) * re_m) * 0.041666666666666664), (im * im), 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))) <= 4e-5) tmp = Float64(fma(-0.16666666666666666, Float64(re_m * re_m), 1.0) * re_m); else tmp = fma(Float64(Float64(Float64(im * im) * re_m) * 0.041666666666666664), Float64(im * im), 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[(N[(0.5 * N[Sin[re$95$m], $MachinePrecision]), $MachinePrecision] * N[(N[Exp[(-im)], $MachinePrecision] + N[Exp[im], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 4e-5], N[(N[(-0.16666666666666666 * N[(re$95$m * re$95$m), $MachinePrecision] + 1.0), $MachinePrecision] * re$95$m), $MachinePrecision], N[(N[(N[(N[(im * im), $MachinePrecision] * re$95$m), $MachinePrecision] * 0.041666666666666664), $MachinePrecision] * N[(im * im), $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}\;\left(0.5 \cdot \sin re\_m\right) \cdot \left(e^{-im} + e^{im}\right) \leq 4 \cdot 10^{-5}:\\
\;\;\;\;\mathsf{fma}\left(-0.16666666666666666, re\_m \cdot re\_m, 1\right) \cdot re\_m\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\left(\left(im \cdot im\right) \cdot re\_m\right) \cdot 0.041666666666666664, im \cdot im, re\_m\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))) < 4.00000000000000033e-5Initial program 100.0%
Taylor expanded in im around 0
lift-sin.f6458.6
Applied rewrites58.6%
Taylor expanded in re around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
lower-fma.f64N/A
pow2N/A
lift-*.f6442.1
Applied rewrites42.1%
if 4.00000000000000033e-5 < (*.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.f6452.0
Applied rewrites52.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-fma.f64N/A
lower-*.f64N/A
pow2N/A
lift-*.f64N/A
*-commutativeN/A
lift-*.f64N/A
pow2N/A
lift-*.f6440.1
Applied rewrites40.1%
Taylor expanded in im around inf
*-commutativeN/A
lower-*.f64N/A
pow2N/A
lift-*.f64N/A
lift-*.f6440.1
Applied rewrites40.1%
Final simplification41.3%
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))) 4e-5)
(* (fma -0.16666666666666666 (* re_m re_m) 1.0) re_m)
(fma (* (* im im) re_m) 0.5 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))) <= 4e-5) {
tmp = fma(-0.16666666666666666, (re_m * re_m), 1.0) * re_m;
} else {
tmp = fma(((im * im) * re_m), 0.5, 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))) <= 4e-5) tmp = Float64(fma(-0.16666666666666666, Float64(re_m * re_m), 1.0) * re_m); else tmp = fma(Float64(Float64(im * im) * re_m), 0.5, 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[(N[(0.5 * N[Sin[re$95$m], $MachinePrecision]), $MachinePrecision] * N[(N[Exp[(-im)], $MachinePrecision] + N[Exp[im], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 4e-5], N[(N[(-0.16666666666666666 * N[(re$95$m * re$95$m), $MachinePrecision] + 1.0), $MachinePrecision] * re$95$m), $MachinePrecision], N[(N[(N[(im * im), $MachinePrecision] * re$95$m), $MachinePrecision] * 0.5 + 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}\;\left(0.5 \cdot \sin re\_m\right) \cdot \left(e^{-im} + e^{im}\right) \leq 4 \cdot 10^{-5}:\\
\;\;\;\;\mathsf{fma}\left(-0.16666666666666666, re\_m \cdot re\_m, 1\right) \cdot re\_m\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\left(im \cdot im\right) \cdot re\_m, 0.5, re\_m\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))) < 4.00000000000000033e-5Initial program 100.0%
Taylor expanded in im around 0
lift-sin.f6458.6
Applied rewrites58.6%
Taylor expanded in re around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
lower-fma.f64N/A
pow2N/A
lift-*.f6442.1
Applied rewrites42.1%
if 4.00000000000000033e-5 < (*.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.f6452.0
Applied rewrites52.0%
Taylor expanded in im around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-*.f64N/A
pow2N/A
lift-*.f6435.1
Applied rewrites35.1%
Final simplification39.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 (* (* 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 (* 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.f6448.7
Applied rewrites48.7%
Taylor expanded in re around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
lower-fma.f64N/A
pow2N/A
lift-*.f6431.3
Applied rewrites31.3%
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.f6465.5
Applied rewrites65.5%
Taylor expanded in im around 0
Applied rewrites25.3%
herbie shell --seed 2025083
(FPCore (re im)
:name "math.sin on complex, real part"
:precision binary64
(* (* 0.5 (sin re)) (+ (exp (- 0.0 im)) (exp im))))