
(FPCore (re im) :precision binary64 (* (* 0.5 (sin re)) (- (exp (- im)) (exp im))))
double code(double re, double im) {
return (0.5 * sin(re)) * (exp(-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(-im) - exp(im))
end function
public static double code(double re, double im) {
return (0.5 * Math.sin(re)) * (Math.exp(-im) - Math.exp(im));
}
def code(re, im): return (0.5 * math.sin(re)) * (math.exp(-im) - math.exp(im))
function code(re, im) return Float64(Float64(0.5 * sin(re)) * Float64(exp(Float64(-im)) - exp(im))) end
function tmp = code(re, im) tmp = (0.5 * sin(re)) * (exp(-im) - exp(im)); end
code[re_, im_] := N[(N[(0.5 * N[Sin[re], $MachinePrecision]), $MachinePrecision] * N[(N[Exp[(-im)], $MachinePrecision] - N[Exp[im], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(0.5 \cdot \sin re\right) \cdot \left(e^{-im} - e^{im}\right)
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 20 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (re im) :precision binary64 (* (* 0.5 (sin re)) (- (exp (- im)) (exp im))))
double code(double re, double im) {
return (0.5 * sin(re)) * (exp(-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(-im) - exp(im))
end function
public static double code(double re, double im) {
return (0.5 * Math.sin(re)) * (Math.exp(-im) - Math.exp(im));
}
def code(re, im): return (0.5 * math.sin(re)) * (math.exp(-im) - math.exp(im))
function code(re, im) return Float64(Float64(0.5 * sin(re)) * Float64(exp(Float64(-im)) - exp(im))) end
function tmp = code(re, im) tmp = (0.5 * sin(re)) * (exp(-im) - exp(im)); end
code[re_, im_] := N[(N[(0.5 * N[Sin[re], $MachinePrecision]), $MachinePrecision] * N[(N[Exp[(-im)], $MachinePrecision] - N[Exp[im], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(0.5 \cdot \sin re\right) \cdot \left(e^{-im} - e^{im}\right)
\end{array}
(FPCore (re im)
:precision binary64
(let* ((t_0 (* 0.5 (sin re))))
(if (<= im -3.8)
(* t_0 (- (exp (- im)) 1.0))
(if (<= im 3.8)
(*
t_0
(*
(-
(*
(-
(*
(*
(- (* -0.0003968253968253968 (* im im)) 0.016666666666666666)
im)
im)
0.3333333333333333)
(* im im))
2.0)
im))
(* t_0 (- 1.0 (exp im)))))))
double code(double re, double im) {
double t_0 = 0.5 * sin(re);
double tmp;
if (im <= -3.8) {
tmp = t_0 * (exp(-im) - 1.0);
} else if (im <= 3.8) {
tmp = t_0 * ((((((((-0.0003968253968253968 * (im * im)) - 0.016666666666666666) * im) * im) - 0.3333333333333333) * (im * im)) - 2.0) * im);
} else {
tmp = t_0 * (1.0 - exp(im));
}
return tmp;
}
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
real(8) :: t_0
real(8) :: tmp
t_0 = 0.5d0 * sin(re)
if (im <= (-3.8d0)) then
tmp = t_0 * (exp(-im) - 1.0d0)
else if (im <= 3.8d0) then
tmp = t_0 * (((((((((-0.0003968253968253968d0) * (im * im)) - 0.016666666666666666d0) * im) * im) - 0.3333333333333333d0) * (im * im)) - 2.0d0) * im)
else
tmp = t_0 * (1.0d0 - exp(im))
end if
code = tmp
end function
public static double code(double re, double im) {
double t_0 = 0.5 * Math.sin(re);
double tmp;
if (im <= -3.8) {
tmp = t_0 * (Math.exp(-im) - 1.0);
} else if (im <= 3.8) {
tmp = t_0 * ((((((((-0.0003968253968253968 * (im * im)) - 0.016666666666666666) * im) * im) - 0.3333333333333333) * (im * im)) - 2.0) * im);
} else {
tmp = t_0 * (1.0 - Math.exp(im));
}
return tmp;
}
def code(re, im): t_0 = 0.5 * math.sin(re) tmp = 0 if im <= -3.8: tmp = t_0 * (math.exp(-im) - 1.0) elif im <= 3.8: tmp = t_0 * ((((((((-0.0003968253968253968 * (im * im)) - 0.016666666666666666) * im) * im) - 0.3333333333333333) * (im * im)) - 2.0) * im) else: tmp = t_0 * (1.0 - math.exp(im)) return tmp
function code(re, im) t_0 = Float64(0.5 * sin(re)) tmp = 0.0 if (im <= -3.8) tmp = Float64(t_0 * Float64(exp(Float64(-im)) - 1.0)); elseif (im <= 3.8) tmp = Float64(t_0 * Float64(Float64(Float64(Float64(Float64(Float64(Float64(Float64(-0.0003968253968253968 * Float64(im * im)) - 0.016666666666666666) * im) * im) - 0.3333333333333333) * Float64(im * im)) - 2.0) * im)); else tmp = Float64(t_0 * Float64(1.0 - exp(im))); end return tmp end
function tmp_2 = code(re, im) t_0 = 0.5 * sin(re); tmp = 0.0; if (im <= -3.8) tmp = t_0 * (exp(-im) - 1.0); elseif (im <= 3.8) tmp = t_0 * ((((((((-0.0003968253968253968 * (im * im)) - 0.016666666666666666) * im) * im) - 0.3333333333333333) * (im * im)) - 2.0) * im); else tmp = t_0 * (1.0 - exp(im)); end tmp_2 = tmp; end
code[re_, im_] := Block[{t$95$0 = N[(0.5 * N[Sin[re], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[im, -3.8], N[(t$95$0 * N[(N[Exp[(-im)], $MachinePrecision] - 1.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[im, 3.8], N[(t$95$0 * N[(N[(N[(N[(N[(N[(N[(N[(-0.0003968253968253968 * N[(im * im), $MachinePrecision]), $MachinePrecision] - 0.016666666666666666), $MachinePrecision] * im), $MachinePrecision] * im), $MachinePrecision] - 0.3333333333333333), $MachinePrecision] * N[(im * im), $MachinePrecision]), $MachinePrecision] - 2.0), $MachinePrecision] * im), $MachinePrecision]), $MachinePrecision], N[(t$95$0 * N[(1.0 - N[Exp[im], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := 0.5 \cdot \sin re\\
\mathbf{if}\;im \leq -3.8:\\
\;\;\;\;t\_0 \cdot \left(e^{-im} - 1\right)\\
\mathbf{elif}\;im \leq 3.8:\\
\;\;\;\;t\_0 \cdot \left(\left(\left(\left(\left(-0.0003968253968253968 \cdot \left(im \cdot im\right) - 0.016666666666666666\right) \cdot im\right) \cdot im - 0.3333333333333333\right) \cdot \left(im \cdot im\right) - 2\right) \cdot im\right)\\
\mathbf{else}:\\
\;\;\;\;t\_0 \cdot \left(1 - e^{im}\right)\\
\end{array}
\end{array}
if im < -3.7999999999999998Initial program 100.0%
Taylor expanded in im around 0
Applied rewrites100.0%
if -3.7999999999999998 < im < 3.7999999999999998Initial program 35.2%
Taylor expanded in im around 0
*-commutativeN/A
lower-*.f64N/A
Applied rewrites99.9%
if 3.7999999999999998 < im Initial program 100.0%
Taylor expanded in im around 0
Applied rewrites100.0%
(FPCore (re im)
:precision binary64
(let* ((t_0 (* 0.5 (sin re))))
(if (<= im 3.8)
(*
t_0
(*
(-
(*
(-
(*
(* (- (* -0.0003968253968253968 (* im im)) 0.016666666666666666) im)
im)
0.3333333333333333)
(* im im))
2.0)
im))
(* t_0 (- 1.0 (exp im))))))
double code(double re, double im) {
double t_0 = 0.5 * sin(re);
double tmp;
if (im <= 3.8) {
tmp = t_0 * ((((((((-0.0003968253968253968 * (im * im)) - 0.016666666666666666) * im) * im) - 0.3333333333333333) * (im * im)) - 2.0) * im);
} else {
tmp = t_0 * (1.0 - exp(im));
}
return tmp;
}
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
real(8) :: t_0
real(8) :: tmp
t_0 = 0.5d0 * sin(re)
if (im <= 3.8d0) then
tmp = t_0 * (((((((((-0.0003968253968253968d0) * (im * im)) - 0.016666666666666666d0) * im) * im) - 0.3333333333333333d0) * (im * im)) - 2.0d0) * im)
else
tmp = t_0 * (1.0d0 - exp(im))
end if
code = tmp
end function
public static double code(double re, double im) {
double t_0 = 0.5 * Math.sin(re);
double tmp;
if (im <= 3.8) {
tmp = t_0 * ((((((((-0.0003968253968253968 * (im * im)) - 0.016666666666666666) * im) * im) - 0.3333333333333333) * (im * im)) - 2.0) * im);
} else {
tmp = t_0 * (1.0 - Math.exp(im));
}
return tmp;
}
def code(re, im): t_0 = 0.5 * math.sin(re) tmp = 0 if im <= 3.8: tmp = t_0 * ((((((((-0.0003968253968253968 * (im * im)) - 0.016666666666666666) * im) * im) - 0.3333333333333333) * (im * im)) - 2.0) * im) else: tmp = t_0 * (1.0 - math.exp(im)) return tmp
function code(re, im) t_0 = Float64(0.5 * sin(re)) tmp = 0.0 if (im <= 3.8) tmp = Float64(t_0 * Float64(Float64(Float64(Float64(Float64(Float64(Float64(Float64(-0.0003968253968253968 * Float64(im * im)) - 0.016666666666666666) * im) * im) - 0.3333333333333333) * Float64(im * im)) - 2.0) * im)); else tmp = Float64(t_0 * Float64(1.0 - exp(im))); end return tmp end
function tmp_2 = code(re, im) t_0 = 0.5 * sin(re); tmp = 0.0; if (im <= 3.8) tmp = t_0 * ((((((((-0.0003968253968253968 * (im * im)) - 0.016666666666666666) * im) * im) - 0.3333333333333333) * (im * im)) - 2.0) * im); else tmp = t_0 * (1.0 - exp(im)); end tmp_2 = tmp; end
code[re_, im_] := Block[{t$95$0 = N[(0.5 * N[Sin[re], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[im, 3.8], N[(t$95$0 * N[(N[(N[(N[(N[(N[(N[(N[(-0.0003968253968253968 * N[(im * im), $MachinePrecision]), $MachinePrecision] - 0.016666666666666666), $MachinePrecision] * im), $MachinePrecision] * im), $MachinePrecision] - 0.3333333333333333), $MachinePrecision] * N[(im * im), $MachinePrecision]), $MachinePrecision] - 2.0), $MachinePrecision] * im), $MachinePrecision]), $MachinePrecision], N[(t$95$0 * N[(1.0 - N[Exp[im], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := 0.5 \cdot \sin re\\
\mathbf{if}\;im \leq 3.8:\\
\;\;\;\;t\_0 \cdot \left(\left(\left(\left(\left(-0.0003968253968253968 \cdot \left(im \cdot im\right) - 0.016666666666666666\right) \cdot im\right) \cdot im - 0.3333333333333333\right) \cdot \left(im \cdot im\right) - 2\right) \cdot im\right)\\
\mathbf{else}:\\
\;\;\;\;t\_0 \cdot \left(1 - e^{im}\right)\\
\end{array}
\end{array}
if im < 3.7999999999999998Initial program 52.9%
Taylor expanded in im around 0
*-commutativeN/A
lower-*.f64N/A
Applied rewrites94.8%
if 3.7999999999999998 < im Initial program 100.0%
Taylor expanded in im around 0
Applied rewrites100.0%
(FPCore (re im)
:precision binary64
(if (<= (* 0.5 (sin re)) -0.05)
(*
(*
(* (fma (* re re) -0.16666666666666666 1.0) re)
(fma (* -0.16666666666666666 im) im -1.0))
im)
(*
(*
(fma
(- (* 0.004166666666666667 (* re re)) 0.08333333333333333)
(* re re)
0.5)
re)
(*
(-
(*
(-
(*
(* (- (* -0.0003968253968253968 (* im im)) 0.016666666666666666) im)
im)
0.3333333333333333)
(* im im))
2.0)
im))))
double code(double re, double im) {
double tmp;
if ((0.5 * sin(re)) <= -0.05) {
tmp = ((fma((re * re), -0.16666666666666666, 1.0) * re) * fma((-0.16666666666666666 * im), im, -1.0)) * im;
} else {
tmp = (fma(((0.004166666666666667 * (re * re)) - 0.08333333333333333), (re * re), 0.5) * re) * ((((((((-0.0003968253968253968 * (im * im)) - 0.016666666666666666) * im) * im) - 0.3333333333333333) * (im * im)) - 2.0) * im);
}
return tmp;
}
function code(re, im) tmp = 0.0 if (Float64(0.5 * sin(re)) <= -0.05) tmp = Float64(Float64(Float64(fma(Float64(re * re), -0.16666666666666666, 1.0) * re) * fma(Float64(-0.16666666666666666 * im), im, -1.0)) * im); else tmp = Float64(Float64(fma(Float64(Float64(0.004166666666666667 * Float64(re * re)) - 0.08333333333333333), Float64(re * re), 0.5) * re) * Float64(Float64(Float64(Float64(Float64(Float64(Float64(Float64(-0.0003968253968253968 * Float64(im * im)) - 0.016666666666666666) * im) * im) - 0.3333333333333333) * Float64(im * im)) - 2.0) * im)); end return tmp end
code[re_, im_] := If[LessEqual[N[(0.5 * N[Sin[re], $MachinePrecision]), $MachinePrecision], -0.05], N[(N[(N[(N[(N[(re * re), $MachinePrecision] * -0.16666666666666666 + 1.0), $MachinePrecision] * re), $MachinePrecision] * N[(N[(-0.16666666666666666 * im), $MachinePrecision] * im + -1.0), $MachinePrecision]), $MachinePrecision] * im), $MachinePrecision], N[(N[(N[(N[(N[(0.004166666666666667 * N[(re * re), $MachinePrecision]), $MachinePrecision] - 0.08333333333333333), $MachinePrecision] * N[(re * re), $MachinePrecision] + 0.5), $MachinePrecision] * re), $MachinePrecision] * N[(N[(N[(N[(N[(N[(N[(N[(-0.0003968253968253968 * N[(im * im), $MachinePrecision]), $MachinePrecision] - 0.016666666666666666), $MachinePrecision] * im), $MachinePrecision] * im), $MachinePrecision] - 0.3333333333333333), $MachinePrecision] * N[(im * im), $MachinePrecision]), $MachinePrecision] - 2.0), $MachinePrecision] * im), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;0.5 \cdot \sin re \leq -0.05:\\
\;\;\;\;\left(\left(\mathsf{fma}\left(re \cdot re, -0.16666666666666666, 1\right) \cdot re\right) \cdot \mathsf{fma}\left(-0.16666666666666666 \cdot im, im, -1\right)\right) \cdot im\\
\mathbf{else}:\\
\;\;\;\;\left(\mathsf{fma}\left(0.004166666666666667 \cdot \left(re \cdot re\right) - 0.08333333333333333, re \cdot re, 0.5\right) \cdot re\right) \cdot \left(\left(\left(\left(\left(-0.0003968253968253968 \cdot \left(im \cdot im\right) - 0.016666666666666666\right) \cdot im\right) \cdot im - 0.3333333333333333\right) \cdot \left(im \cdot im\right) - 2\right) \cdot im\right)\\
\end{array}
\end{array}
if (*.f64 #s(literal 1/2 binary64) (sin.f64 re)) < -0.050000000000000003Initial program 51.3%
Taylor expanded in im around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
associate-*r*N/A
distribute-rgt-outN/A
lower-*.f64N/A
lift-sin.f64N/A
unpow2N/A
associate-*r*N/A
lower-fma.f64N/A
lower-*.f6486.3
Applied rewrites86.3%
Taylor expanded in re around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6428.3
Applied rewrites28.3%
if -0.050000000000000003 < (*.f64 #s(literal 1/2 binary64) (sin.f64 re)) Initial program 70.2%
Taylor expanded in im around 0
*-commutativeN/A
lower-*.f64N/A
Applied rewrites93.0%
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
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6475.0
Applied rewrites75.0%
(FPCore (re im)
:precision binary64
(let* ((t_0
(*
(* 0.5 (sin re))
(*
(-
(*
(* (- (* -0.016666666666666666 (* im im)) 0.3333333333333333) im)
im)
2.0)
im))))
(if (<= im -1.2e+64)
t_0
(if (<= im -8.2)
(*
(*
0.5
(*
(fma
(- (* 0.008333333333333333 (* re re)) 0.16666666666666666)
(* re re)
1.0)
re))
(- (exp (- im)) 1.0))
(if (or (<= im 7.6) (not (<= im 1.01e+62)))
t_0
(*
(* 0.5 (* (fma (* re re) -0.16666666666666666 1.0) re))
(- 1.0 (exp im))))))))
double code(double re, double im) {
double t_0 = (0.5 * sin(re)) * ((((((-0.016666666666666666 * (im * im)) - 0.3333333333333333) * im) * im) - 2.0) * im);
double tmp;
if (im <= -1.2e+64) {
tmp = t_0;
} else if (im <= -8.2) {
tmp = (0.5 * (fma(((0.008333333333333333 * (re * re)) - 0.16666666666666666), (re * re), 1.0) * re)) * (exp(-im) - 1.0);
} else if ((im <= 7.6) || !(im <= 1.01e+62)) {
tmp = t_0;
} else {
tmp = (0.5 * (fma((re * re), -0.16666666666666666, 1.0) * re)) * (1.0 - exp(im));
}
return tmp;
}
function code(re, im) t_0 = Float64(Float64(0.5 * sin(re)) * Float64(Float64(Float64(Float64(Float64(Float64(-0.016666666666666666 * Float64(im * im)) - 0.3333333333333333) * im) * im) - 2.0) * im)) tmp = 0.0 if (im <= -1.2e+64) tmp = t_0; elseif (im <= -8.2) tmp = Float64(Float64(0.5 * Float64(fma(Float64(Float64(0.008333333333333333 * Float64(re * re)) - 0.16666666666666666), Float64(re * re), 1.0) * re)) * Float64(exp(Float64(-im)) - 1.0)); elseif ((im <= 7.6) || !(im <= 1.01e+62)) tmp = t_0; else tmp = Float64(Float64(0.5 * Float64(fma(Float64(re * re), -0.16666666666666666, 1.0) * re)) * Float64(1.0 - exp(im))); end return tmp end
code[re_, im_] := Block[{t$95$0 = N[(N[(0.5 * N[Sin[re], $MachinePrecision]), $MachinePrecision] * N[(N[(N[(N[(N[(N[(-0.016666666666666666 * N[(im * im), $MachinePrecision]), $MachinePrecision] - 0.3333333333333333), $MachinePrecision] * im), $MachinePrecision] * im), $MachinePrecision] - 2.0), $MachinePrecision] * im), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[im, -1.2e+64], t$95$0, If[LessEqual[im, -8.2], N[(N[(0.5 * N[(N[(N[(N[(0.008333333333333333 * N[(re * re), $MachinePrecision]), $MachinePrecision] - 0.16666666666666666), $MachinePrecision] * N[(re * re), $MachinePrecision] + 1.0), $MachinePrecision] * re), $MachinePrecision]), $MachinePrecision] * N[(N[Exp[(-im)], $MachinePrecision] - 1.0), $MachinePrecision]), $MachinePrecision], If[Or[LessEqual[im, 7.6], N[Not[LessEqual[im, 1.01e+62]], $MachinePrecision]], t$95$0, N[(N[(0.5 * N[(N[(N[(re * re), $MachinePrecision] * -0.16666666666666666 + 1.0), $MachinePrecision] * re), $MachinePrecision]), $MachinePrecision] * N[(1.0 - N[Exp[im], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(0.5 \cdot \sin re\right) \cdot \left(\left(\left(\left(-0.016666666666666666 \cdot \left(im \cdot im\right) - 0.3333333333333333\right) \cdot im\right) \cdot im - 2\right) \cdot im\right)\\
\mathbf{if}\;im \leq -1.2 \cdot 10^{+64}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;im \leq -8.2:\\
\;\;\;\;\left(0.5 \cdot \left(\mathsf{fma}\left(0.008333333333333333 \cdot \left(re \cdot re\right) - 0.16666666666666666, re \cdot re, 1\right) \cdot re\right)\right) \cdot \left(e^{-im} - 1\right)\\
\mathbf{elif}\;im \leq 7.6 \lor \neg \left(im \leq 1.01 \cdot 10^{+62}\right):\\
\;\;\;\;t\_0\\
\mathbf{else}:\\
\;\;\;\;\left(0.5 \cdot \left(\mathsf{fma}\left(re \cdot re, -0.16666666666666666, 1\right) \cdot re\right)\right) \cdot \left(1 - e^{im}\right)\\
\end{array}
\end{array}
if im < -1.2e64 or -8.1999999999999993 < im < 7.5999999999999996 or 1.01e62 < im Initial program 61.7%
Taylor expanded in im around 0
*-commutativeN/A
lower-*.f64N/A
lower--.f64N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
lower--.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6499.8
Applied rewrites99.8%
if -1.2e64 < im < -8.1999999999999993Initial program 100.0%
Taylor expanded in im around 0
Applied rewrites100.0%
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-*.f6483.3
Applied rewrites83.3%
if 7.5999999999999996 < im < 1.01e62Initial program 100.0%
Taylor expanded in im around 0
Applied rewrites100.0%
Taylor expanded in re around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6485.7
Applied rewrites85.7%
Final simplification98.3%
(FPCore (re im)
:precision binary64
(if (<= (* 0.5 (sin re)) -0.05)
(*
(*
(* (fma (* re re) -0.16666666666666666 1.0) re)
(fma (* -0.16666666666666666 im) im -1.0))
im)
(*
(* 0.5 re)
(*
(-
(*
(*
(-
(*
(* (- (* (* -0.0003968253968253968 im) im) 0.016666666666666666) im)
im)
0.3333333333333333)
im)
im)
2.0)
im))))
double code(double re, double im) {
double tmp;
if ((0.5 * sin(re)) <= -0.05) {
tmp = ((fma((re * re), -0.16666666666666666, 1.0) * re) * fma((-0.16666666666666666 * im), im, -1.0)) * im;
} else {
tmp = (0.5 * re) * ((((((((((-0.0003968253968253968 * im) * im) - 0.016666666666666666) * im) * im) - 0.3333333333333333) * im) * im) - 2.0) * im);
}
return tmp;
}
function code(re, im) tmp = 0.0 if (Float64(0.5 * sin(re)) <= -0.05) tmp = Float64(Float64(Float64(fma(Float64(re * re), -0.16666666666666666, 1.0) * re) * fma(Float64(-0.16666666666666666 * im), im, -1.0)) * im); else tmp = Float64(Float64(0.5 * re) * Float64(Float64(Float64(Float64(Float64(Float64(Float64(Float64(Float64(Float64(-0.0003968253968253968 * im) * im) - 0.016666666666666666) * im) * im) - 0.3333333333333333) * im) * im) - 2.0) * im)); end return tmp end
code[re_, im_] := If[LessEqual[N[(0.5 * N[Sin[re], $MachinePrecision]), $MachinePrecision], -0.05], N[(N[(N[(N[(N[(re * re), $MachinePrecision] * -0.16666666666666666 + 1.0), $MachinePrecision] * re), $MachinePrecision] * N[(N[(-0.16666666666666666 * im), $MachinePrecision] * im + -1.0), $MachinePrecision]), $MachinePrecision] * im), $MachinePrecision], N[(N[(0.5 * re), $MachinePrecision] * N[(N[(N[(N[(N[(N[(N[(N[(N[(N[(-0.0003968253968253968 * im), $MachinePrecision] * im), $MachinePrecision] - 0.016666666666666666), $MachinePrecision] * im), $MachinePrecision] * im), $MachinePrecision] - 0.3333333333333333), $MachinePrecision] * im), $MachinePrecision] * im), $MachinePrecision] - 2.0), $MachinePrecision] * im), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;0.5 \cdot \sin re \leq -0.05:\\
\;\;\;\;\left(\left(\mathsf{fma}\left(re \cdot re, -0.16666666666666666, 1\right) \cdot re\right) \cdot \mathsf{fma}\left(-0.16666666666666666 \cdot im, im, -1\right)\right) \cdot im\\
\mathbf{else}:\\
\;\;\;\;\left(0.5 \cdot re\right) \cdot \left(\left(\left(\left(\left(\left(\left(-0.0003968253968253968 \cdot im\right) \cdot im - 0.016666666666666666\right) \cdot im\right) \cdot im - 0.3333333333333333\right) \cdot im\right) \cdot im - 2\right) \cdot im\right)\\
\end{array}
\end{array}
if (*.f64 #s(literal 1/2 binary64) (sin.f64 re)) < -0.050000000000000003Initial program 51.3%
Taylor expanded in im around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
associate-*r*N/A
distribute-rgt-outN/A
lower-*.f64N/A
lift-sin.f64N/A
unpow2N/A
associate-*r*N/A
lower-fma.f64N/A
lower-*.f6486.3
Applied rewrites86.3%
Taylor expanded in re around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6428.3
Applied rewrites28.3%
if -0.050000000000000003 < (*.f64 #s(literal 1/2 binary64) (sin.f64 re)) Initial program 70.2%
Taylor expanded in im around 0
*-commutativeN/A
lower-*.f64N/A
Applied rewrites93.0%
lift-*.f64N/A
lift-sin.f64N/A
*-commutativeN/A
lower-*.f64N/A
lift-sin.f6493.0
Applied rewrites93.0%
Taylor expanded in re around 0
*-commutativeN/A
lower-*.f6474.3
Applied rewrites74.3%
lift-*.f64N/A
lift-*.f64N/A
lift--.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift--.f64N/A
lift-*.f64N/A
lift-*.f64N/A
associate-*r*N/A
lower-*.f64N/A
Applied rewrites74.3%
(FPCore (re im)
:precision binary64
(if (<= (* 0.5 (sin re)) -0.05)
(*
(*
(* (fma (* re re) -0.16666666666666666 1.0) re)
(fma (* -0.16666666666666666 im) im -1.0))
im)
(*
(* 0.5 re)
(*
(-
(*
(-
(* (* (* -0.0003968253968253968 im) im) (* im im))
0.3333333333333333)
(* im im))
2.0)
im))))
double code(double re, double im) {
double tmp;
if ((0.5 * sin(re)) <= -0.05) {
tmp = ((fma((re * re), -0.16666666666666666, 1.0) * re) * fma((-0.16666666666666666 * im), im, -1.0)) * im;
} else {
tmp = (0.5 * re) * (((((((-0.0003968253968253968 * im) * im) * (im * im)) - 0.3333333333333333) * (im * im)) - 2.0) * im);
}
return tmp;
}
function code(re, im) tmp = 0.0 if (Float64(0.5 * sin(re)) <= -0.05) tmp = Float64(Float64(Float64(fma(Float64(re * re), -0.16666666666666666, 1.0) * re) * fma(Float64(-0.16666666666666666 * im), im, -1.0)) * im); else tmp = Float64(Float64(0.5 * re) * Float64(Float64(Float64(Float64(Float64(Float64(Float64(-0.0003968253968253968 * im) * im) * Float64(im * im)) - 0.3333333333333333) * Float64(im * im)) - 2.0) * im)); end return tmp end
code[re_, im_] := If[LessEqual[N[(0.5 * N[Sin[re], $MachinePrecision]), $MachinePrecision], -0.05], N[(N[(N[(N[(N[(re * re), $MachinePrecision] * -0.16666666666666666 + 1.0), $MachinePrecision] * re), $MachinePrecision] * N[(N[(-0.16666666666666666 * im), $MachinePrecision] * im + -1.0), $MachinePrecision]), $MachinePrecision] * im), $MachinePrecision], N[(N[(0.5 * re), $MachinePrecision] * N[(N[(N[(N[(N[(N[(N[(-0.0003968253968253968 * im), $MachinePrecision] * im), $MachinePrecision] * N[(im * im), $MachinePrecision]), $MachinePrecision] - 0.3333333333333333), $MachinePrecision] * N[(im * im), $MachinePrecision]), $MachinePrecision] - 2.0), $MachinePrecision] * im), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;0.5 \cdot \sin re \leq -0.05:\\
\;\;\;\;\left(\left(\mathsf{fma}\left(re \cdot re, -0.16666666666666666, 1\right) \cdot re\right) \cdot \mathsf{fma}\left(-0.16666666666666666 \cdot im, im, -1\right)\right) \cdot im\\
\mathbf{else}:\\
\;\;\;\;\left(0.5 \cdot re\right) \cdot \left(\left(\left(\left(\left(-0.0003968253968253968 \cdot im\right) \cdot im\right) \cdot \left(im \cdot im\right) - 0.3333333333333333\right) \cdot \left(im \cdot im\right) - 2\right) \cdot im\right)\\
\end{array}
\end{array}
if (*.f64 #s(literal 1/2 binary64) (sin.f64 re)) < -0.050000000000000003Initial program 51.3%
Taylor expanded in im around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
associate-*r*N/A
distribute-rgt-outN/A
lower-*.f64N/A
lift-sin.f64N/A
unpow2N/A
associate-*r*N/A
lower-fma.f64N/A
lower-*.f6486.3
Applied rewrites86.3%
Taylor expanded in re around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6428.3
Applied rewrites28.3%
if -0.050000000000000003 < (*.f64 #s(literal 1/2 binary64) (sin.f64 re)) Initial program 70.2%
Taylor expanded in im around 0
*-commutativeN/A
lower-*.f64N/A
Applied rewrites93.0%
lift-*.f64N/A
lift-sin.f64N/A
*-commutativeN/A
lower-*.f64N/A
lift-sin.f6493.0
Applied rewrites93.0%
Taylor expanded in re around 0
*-commutativeN/A
lower-*.f6474.3
Applied rewrites74.3%
Taylor expanded in im around inf
pow2N/A
associate-*l*N/A
lift-*.f64N/A
lift-*.f6474.3
Applied rewrites74.3%
(FPCore (re im)
:precision binary64
(*
(* 0.5 (sin re))
(*
(-
(*
(-
(*
(* (- (* -0.0003968253968253968 (* im im)) 0.016666666666666666) im)
im)
0.3333333333333333)
(* im im))
2.0)
im)))
double code(double re, double im) {
return (0.5 * sin(re)) * ((((((((-0.0003968253968253968 * (im * im)) - 0.016666666666666666) * im) * im) - 0.3333333333333333) * (im * im)) - 2.0) * 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)) * (((((((((-0.0003968253968253968d0) * (im * im)) - 0.016666666666666666d0) * im) * im) - 0.3333333333333333d0) * (im * im)) - 2.0d0) * im)
end function
public static double code(double re, double im) {
return (0.5 * Math.sin(re)) * ((((((((-0.0003968253968253968 * (im * im)) - 0.016666666666666666) * im) * im) - 0.3333333333333333) * (im * im)) - 2.0) * im);
}
def code(re, im): return (0.5 * math.sin(re)) * ((((((((-0.0003968253968253968 * (im * im)) - 0.016666666666666666) * im) * im) - 0.3333333333333333) * (im * im)) - 2.0) * im)
function code(re, im) return Float64(Float64(0.5 * sin(re)) * Float64(Float64(Float64(Float64(Float64(Float64(Float64(Float64(-0.0003968253968253968 * Float64(im * im)) - 0.016666666666666666) * im) * im) - 0.3333333333333333) * Float64(im * im)) - 2.0) * im)) end
function tmp = code(re, im) tmp = (0.5 * sin(re)) * ((((((((-0.0003968253968253968 * (im * im)) - 0.016666666666666666) * im) * im) - 0.3333333333333333) * (im * im)) - 2.0) * im); end
code[re_, im_] := N[(N[(0.5 * N[Sin[re], $MachinePrecision]), $MachinePrecision] * N[(N[(N[(N[(N[(N[(N[(N[(-0.0003968253968253968 * N[(im * im), $MachinePrecision]), $MachinePrecision] - 0.016666666666666666), $MachinePrecision] * im), $MachinePrecision] * im), $MachinePrecision] - 0.3333333333333333), $MachinePrecision] * N[(im * im), $MachinePrecision]), $MachinePrecision] - 2.0), $MachinePrecision] * im), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(0.5 \cdot \sin re\right) \cdot \left(\left(\left(\left(\left(-0.0003968253968253968 \cdot \left(im \cdot im\right) - 0.016666666666666666\right) \cdot im\right) \cdot im - 0.3333333333333333\right) \cdot \left(im \cdot im\right) - 2\right) \cdot im\right)
\end{array}
Initial program 65.6%
Taylor expanded in im around 0
*-commutativeN/A
lower-*.f64N/A
Applied rewrites93.3%
(FPCore (re im)
:precision binary64
(let* ((t_0
(*
(* 0.5 (sin re))
(* (- (* -0.3333333333333333 (* im im)) 2.0) im))))
(if (<= im -8e+102)
t_0
(if (<= im -8.2)
(*
(* 0.5 (* (fma -0.16666666666666666 (* re re) 1.0) re))
(- (exp (- im)) 1.0))
(if (<= im 7.4)
(* (* (sin re) im) (- (* (* im im) -0.16666666666666666) 1.0))
(if (<= im 8e+102)
(*
(* 0.5 (* (fma (* re re) -0.16666666666666666 1.0) re))
(- 1.0 (exp im)))
t_0))))))
double code(double re, double im) {
double t_0 = (0.5 * sin(re)) * (((-0.3333333333333333 * (im * im)) - 2.0) * im);
double tmp;
if (im <= -8e+102) {
tmp = t_0;
} else if (im <= -8.2) {
tmp = (0.5 * (fma(-0.16666666666666666, (re * re), 1.0) * re)) * (exp(-im) - 1.0);
} else if (im <= 7.4) {
tmp = (sin(re) * im) * (((im * im) * -0.16666666666666666) - 1.0);
} else if (im <= 8e+102) {
tmp = (0.5 * (fma((re * re), -0.16666666666666666, 1.0) * re)) * (1.0 - exp(im));
} else {
tmp = t_0;
}
return tmp;
}
function code(re, im) t_0 = Float64(Float64(0.5 * sin(re)) * Float64(Float64(Float64(-0.3333333333333333 * Float64(im * im)) - 2.0) * im)) tmp = 0.0 if (im <= -8e+102) tmp = t_0; elseif (im <= -8.2) tmp = Float64(Float64(0.5 * Float64(fma(-0.16666666666666666, Float64(re * re), 1.0) * re)) * Float64(exp(Float64(-im)) - 1.0)); elseif (im <= 7.4) tmp = Float64(Float64(sin(re) * im) * Float64(Float64(Float64(im * im) * -0.16666666666666666) - 1.0)); elseif (im <= 8e+102) tmp = Float64(Float64(0.5 * Float64(fma(Float64(re * re), -0.16666666666666666, 1.0) * re)) * Float64(1.0 - exp(im))); else tmp = t_0; end return tmp end
code[re_, im_] := Block[{t$95$0 = N[(N[(0.5 * N[Sin[re], $MachinePrecision]), $MachinePrecision] * N[(N[(N[(-0.3333333333333333 * N[(im * im), $MachinePrecision]), $MachinePrecision] - 2.0), $MachinePrecision] * im), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[im, -8e+102], t$95$0, If[LessEqual[im, -8.2], N[(N[(0.5 * N[(N[(-0.16666666666666666 * N[(re * re), $MachinePrecision] + 1.0), $MachinePrecision] * re), $MachinePrecision]), $MachinePrecision] * N[(N[Exp[(-im)], $MachinePrecision] - 1.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[im, 7.4], N[(N[(N[Sin[re], $MachinePrecision] * im), $MachinePrecision] * N[(N[(N[(im * im), $MachinePrecision] * -0.16666666666666666), $MachinePrecision] - 1.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[im, 8e+102], N[(N[(0.5 * N[(N[(N[(re * re), $MachinePrecision] * -0.16666666666666666 + 1.0), $MachinePrecision] * re), $MachinePrecision]), $MachinePrecision] * N[(1.0 - N[Exp[im], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$0]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(0.5 \cdot \sin re\right) \cdot \left(\left(-0.3333333333333333 \cdot \left(im \cdot im\right) - 2\right) \cdot im\right)\\
\mathbf{if}\;im \leq -8 \cdot 10^{+102}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;im \leq -8.2:\\
\;\;\;\;\left(0.5 \cdot \left(\mathsf{fma}\left(-0.16666666666666666, re \cdot re, 1\right) \cdot re\right)\right) \cdot \left(e^{-im} - 1\right)\\
\mathbf{elif}\;im \leq 7.4:\\
\;\;\;\;\left(\sin re \cdot im\right) \cdot \left(\left(im \cdot im\right) \cdot -0.16666666666666666 - 1\right)\\
\mathbf{elif}\;im \leq 8 \cdot 10^{+102}:\\
\;\;\;\;\left(0.5 \cdot \left(\mathsf{fma}\left(re \cdot re, -0.16666666666666666, 1\right) \cdot re\right)\right) \cdot \left(1 - e^{im}\right)\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if im < -7.99999999999999982e102 or 7.99999999999999982e102 < im Initial program 100.0%
Taylor expanded in im around 0
*-commutativeN/A
lower-*.f64N/A
lower--.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64100.0
Applied rewrites100.0%
if -7.99999999999999982e102 < im < -8.1999999999999993Initial program 100.0%
Taylor expanded in im around 0
Applied rewrites100.0%
Taylor expanded in re around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
lower-fma.f64N/A
pow2N/A
lift-*.f6458.8
Applied rewrites58.8%
if -8.1999999999999993 < im < 7.4000000000000004Initial program 35.2%
Taylor expanded in im around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
associate-*r*N/A
distribute-rgt-outN/A
lower-*.f64N/A
lift-sin.f64N/A
unpow2N/A
associate-*r*N/A
lower-fma.f64N/A
lower-*.f6499.6
Applied rewrites99.6%
Taylor expanded in re around inf
associate-*r*N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lift-sin.f64N/A
lower--.f64N/A
*-commutativeN/A
lower-*.f64N/A
pow2N/A
lift-*.f6499.6
Applied rewrites99.6%
if 7.4000000000000004 < im < 7.99999999999999982e102Initial program 100.0%
Taylor expanded in im around 0
Applied rewrites100.0%
Taylor expanded in re around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6487.5
Applied rewrites87.5%
(FPCore (re im)
:precision binary64
(if (<= im -2.1e+137)
(* (* (sin re) (fma (* -0.16666666666666666 im) im -1.0)) im)
(if (<= im -0.048)
(*
(*
(fma
(- (* 0.004166666666666667 (* re re)) 0.08333333333333333)
(* re re)
0.5)
re)
(*
(-
(*
(-
(*
(* (- (* -0.0003968253968253968 (* im im)) 0.016666666666666666) im)
im)
0.3333333333333333)
(* im im))
2.0)
im))
(if (<= im 7.4)
(* (* (sin re) im) (- (* (* im im) -0.16666666666666666) 1.0))
(if (<= im 8e+102)
(*
(* 0.5 (* (fma (* re re) -0.16666666666666666 1.0) re))
(- 1.0 (exp im)))
(*
(* 0.5 (sin re))
(* (- (* -0.3333333333333333 (* im im)) 2.0) im)))))))
double code(double re, double im) {
double tmp;
if (im <= -2.1e+137) {
tmp = (sin(re) * fma((-0.16666666666666666 * im), im, -1.0)) * im;
} else if (im <= -0.048) {
tmp = (fma(((0.004166666666666667 * (re * re)) - 0.08333333333333333), (re * re), 0.5) * re) * ((((((((-0.0003968253968253968 * (im * im)) - 0.016666666666666666) * im) * im) - 0.3333333333333333) * (im * im)) - 2.0) * im);
} else if (im <= 7.4) {
tmp = (sin(re) * im) * (((im * im) * -0.16666666666666666) - 1.0);
} else if (im <= 8e+102) {
tmp = (0.5 * (fma((re * re), -0.16666666666666666, 1.0) * re)) * (1.0 - exp(im));
} else {
tmp = (0.5 * sin(re)) * (((-0.3333333333333333 * (im * im)) - 2.0) * im);
}
return tmp;
}
function code(re, im) tmp = 0.0 if (im <= -2.1e+137) tmp = Float64(Float64(sin(re) * fma(Float64(-0.16666666666666666 * im), im, -1.0)) * im); elseif (im <= -0.048) tmp = Float64(Float64(fma(Float64(Float64(0.004166666666666667 * Float64(re * re)) - 0.08333333333333333), Float64(re * re), 0.5) * re) * Float64(Float64(Float64(Float64(Float64(Float64(Float64(Float64(-0.0003968253968253968 * Float64(im * im)) - 0.016666666666666666) * im) * im) - 0.3333333333333333) * Float64(im * im)) - 2.0) * im)); elseif (im <= 7.4) tmp = Float64(Float64(sin(re) * im) * Float64(Float64(Float64(im * im) * -0.16666666666666666) - 1.0)); elseif (im <= 8e+102) tmp = Float64(Float64(0.5 * Float64(fma(Float64(re * re), -0.16666666666666666, 1.0) * re)) * Float64(1.0 - exp(im))); else tmp = Float64(Float64(0.5 * sin(re)) * Float64(Float64(Float64(-0.3333333333333333 * Float64(im * im)) - 2.0) * im)); end return tmp end
code[re_, im_] := If[LessEqual[im, -2.1e+137], N[(N[(N[Sin[re], $MachinePrecision] * N[(N[(-0.16666666666666666 * im), $MachinePrecision] * im + -1.0), $MachinePrecision]), $MachinePrecision] * im), $MachinePrecision], If[LessEqual[im, -0.048], N[(N[(N[(N[(N[(0.004166666666666667 * N[(re * re), $MachinePrecision]), $MachinePrecision] - 0.08333333333333333), $MachinePrecision] * N[(re * re), $MachinePrecision] + 0.5), $MachinePrecision] * re), $MachinePrecision] * N[(N[(N[(N[(N[(N[(N[(N[(-0.0003968253968253968 * N[(im * im), $MachinePrecision]), $MachinePrecision] - 0.016666666666666666), $MachinePrecision] * im), $MachinePrecision] * im), $MachinePrecision] - 0.3333333333333333), $MachinePrecision] * N[(im * im), $MachinePrecision]), $MachinePrecision] - 2.0), $MachinePrecision] * im), $MachinePrecision]), $MachinePrecision], If[LessEqual[im, 7.4], N[(N[(N[Sin[re], $MachinePrecision] * im), $MachinePrecision] * N[(N[(N[(im * im), $MachinePrecision] * -0.16666666666666666), $MachinePrecision] - 1.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[im, 8e+102], N[(N[(0.5 * N[(N[(N[(re * re), $MachinePrecision] * -0.16666666666666666 + 1.0), $MachinePrecision] * re), $MachinePrecision]), $MachinePrecision] * N[(1.0 - N[Exp[im], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(0.5 * N[Sin[re], $MachinePrecision]), $MachinePrecision] * N[(N[(N[(-0.3333333333333333 * N[(im * im), $MachinePrecision]), $MachinePrecision] - 2.0), $MachinePrecision] * im), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;im \leq -2.1 \cdot 10^{+137}:\\
\;\;\;\;\left(\sin re \cdot \mathsf{fma}\left(-0.16666666666666666 \cdot im, im, -1\right)\right) \cdot im\\
\mathbf{elif}\;im \leq -0.048:\\
\;\;\;\;\left(\mathsf{fma}\left(0.004166666666666667 \cdot \left(re \cdot re\right) - 0.08333333333333333, re \cdot re, 0.5\right) \cdot re\right) \cdot \left(\left(\left(\left(\left(-0.0003968253968253968 \cdot \left(im \cdot im\right) - 0.016666666666666666\right) \cdot im\right) \cdot im - 0.3333333333333333\right) \cdot \left(im \cdot im\right) - 2\right) \cdot im\right)\\
\mathbf{elif}\;im \leq 7.4:\\
\;\;\;\;\left(\sin re \cdot im\right) \cdot \left(\left(im \cdot im\right) \cdot -0.16666666666666666 - 1\right)\\
\mathbf{elif}\;im \leq 8 \cdot 10^{+102}:\\
\;\;\;\;\left(0.5 \cdot \left(\mathsf{fma}\left(re \cdot re, -0.16666666666666666, 1\right) \cdot re\right)\right) \cdot \left(1 - e^{im}\right)\\
\mathbf{else}:\\
\;\;\;\;\left(0.5 \cdot \sin re\right) \cdot \left(\left(-0.3333333333333333 \cdot \left(im \cdot im\right) - 2\right) \cdot im\right)\\
\end{array}
\end{array}
if im < -2.0999999999999999e137Initial program 100.0%
Taylor expanded in im around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
associate-*r*N/A
distribute-rgt-outN/A
lower-*.f64N/A
lift-sin.f64N/A
unpow2N/A
associate-*r*N/A
lower-fma.f64N/A
lower-*.f64100.0
Applied rewrites100.0%
if -2.0999999999999999e137 < im < -0.048000000000000001Initial program 100.0%
Taylor expanded in im around 0
*-commutativeN/A
lower-*.f64N/A
Applied rewrites54.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
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6462.5
Applied rewrites62.5%
if -0.048000000000000001 < im < 7.4000000000000004Initial program 35.2%
Taylor expanded in im around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
associate-*r*N/A
distribute-rgt-outN/A
lower-*.f64N/A
lift-sin.f64N/A
unpow2N/A
associate-*r*N/A
lower-fma.f64N/A
lower-*.f6499.6
Applied rewrites99.6%
Taylor expanded in re around inf
associate-*r*N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lift-sin.f64N/A
lower--.f64N/A
*-commutativeN/A
lower-*.f64N/A
pow2N/A
lift-*.f6499.6
Applied rewrites99.6%
if 7.4000000000000004 < im < 7.99999999999999982e102Initial program 100.0%
Taylor expanded in im around 0
Applied rewrites100.0%
Taylor expanded in re around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6487.5
Applied rewrites87.5%
if 7.99999999999999982e102 < im Initial program 100.0%
Taylor expanded in im around 0
*-commutativeN/A
lower-*.f64N/A
lower--.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64100.0
Applied rewrites100.0%
(FPCore (re im)
:precision binary64
(if (or (<= im 7.6) (not (<= im 1.01e+62)))
(*
(* 0.5 (sin re))
(*
(-
(* (* (- (* -0.016666666666666666 (* im im)) 0.3333333333333333) im) im)
2.0)
im))
(*
(* 0.5 (* (fma (* re re) -0.16666666666666666 1.0) re))
(- 1.0 (exp im)))))
double code(double re, double im) {
double tmp;
if ((im <= 7.6) || !(im <= 1.01e+62)) {
tmp = (0.5 * sin(re)) * ((((((-0.016666666666666666 * (im * im)) - 0.3333333333333333) * im) * im) - 2.0) * im);
} else {
tmp = (0.5 * (fma((re * re), -0.16666666666666666, 1.0) * re)) * (1.0 - exp(im));
}
return tmp;
}
function code(re, im) tmp = 0.0 if ((im <= 7.6) || !(im <= 1.01e+62)) tmp = Float64(Float64(0.5 * sin(re)) * Float64(Float64(Float64(Float64(Float64(Float64(-0.016666666666666666 * Float64(im * im)) - 0.3333333333333333) * im) * im) - 2.0) * im)); else tmp = Float64(Float64(0.5 * Float64(fma(Float64(re * re), -0.16666666666666666, 1.0) * re)) * Float64(1.0 - exp(im))); end return tmp end
code[re_, im_] := If[Or[LessEqual[im, 7.6], N[Not[LessEqual[im, 1.01e+62]], $MachinePrecision]], N[(N[(0.5 * N[Sin[re], $MachinePrecision]), $MachinePrecision] * N[(N[(N[(N[(N[(N[(-0.016666666666666666 * N[(im * im), $MachinePrecision]), $MachinePrecision] - 0.3333333333333333), $MachinePrecision] * im), $MachinePrecision] * im), $MachinePrecision] - 2.0), $MachinePrecision] * im), $MachinePrecision]), $MachinePrecision], N[(N[(0.5 * N[(N[(N[(re * re), $MachinePrecision] * -0.16666666666666666 + 1.0), $MachinePrecision] * re), $MachinePrecision]), $MachinePrecision] * N[(1.0 - N[Exp[im], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;im \leq 7.6 \lor \neg \left(im \leq 1.01 \cdot 10^{+62}\right):\\
\;\;\;\;\left(0.5 \cdot \sin re\right) \cdot \left(\left(\left(\left(-0.016666666666666666 \cdot \left(im \cdot im\right) - 0.3333333333333333\right) \cdot im\right) \cdot im - 2\right) \cdot im\right)\\
\mathbf{else}:\\
\;\;\;\;\left(0.5 \cdot \left(\mathsf{fma}\left(re \cdot re, -0.16666666666666666, 1\right) \cdot re\right)\right) \cdot \left(1 - e^{im}\right)\\
\end{array}
\end{array}
if im < 7.5999999999999996 or 1.01e62 < im Initial program 63.6%
Taylor expanded in im around 0
*-commutativeN/A
lower-*.f64N/A
lower--.f64N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
lower--.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6495.5
Applied rewrites95.5%
if 7.5999999999999996 < im < 1.01e62Initial program 100.0%
Taylor expanded in im around 0
Applied rewrites100.0%
Taylor expanded in re around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6485.7
Applied rewrites85.7%
Final simplification94.9%
(FPCore (re im)
:precision binary64
(if (<= (* 0.5 (sin re)) -0.05)
(*
(*
(* (fma (* re re) -0.16666666666666666 1.0) re)
(fma (* -0.16666666666666666 im) im -1.0))
im)
(*
(* 0.5 re)
(*
(-
(* (* (- (* -0.016666666666666666 (* im im)) 0.3333333333333333) im) im)
2.0)
im))))
double code(double re, double im) {
double tmp;
if ((0.5 * sin(re)) <= -0.05) {
tmp = ((fma((re * re), -0.16666666666666666, 1.0) * re) * fma((-0.16666666666666666 * im), im, -1.0)) * im;
} else {
tmp = (0.5 * re) * ((((((-0.016666666666666666 * (im * im)) - 0.3333333333333333) * im) * im) - 2.0) * im);
}
return tmp;
}
function code(re, im) tmp = 0.0 if (Float64(0.5 * sin(re)) <= -0.05) tmp = Float64(Float64(Float64(fma(Float64(re * re), -0.16666666666666666, 1.0) * re) * fma(Float64(-0.16666666666666666 * im), im, -1.0)) * im); else tmp = Float64(Float64(0.5 * re) * Float64(Float64(Float64(Float64(Float64(Float64(-0.016666666666666666 * Float64(im * im)) - 0.3333333333333333) * im) * im) - 2.0) * im)); end return tmp end
code[re_, im_] := If[LessEqual[N[(0.5 * N[Sin[re], $MachinePrecision]), $MachinePrecision], -0.05], N[(N[(N[(N[(N[(re * re), $MachinePrecision] * -0.16666666666666666 + 1.0), $MachinePrecision] * re), $MachinePrecision] * N[(N[(-0.16666666666666666 * im), $MachinePrecision] * im + -1.0), $MachinePrecision]), $MachinePrecision] * im), $MachinePrecision], N[(N[(0.5 * re), $MachinePrecision] * N[(N[(N[(N[(N[(N[(-0.016666666666666666 * N[(im * im), $MachinePrecision]), $MachinePrecision] - 0.3333333333333333), $MachinePrecision] * im), $MachinePrecision] * im), $MachinePrecision] - 2.0), $MachinePrecision] * im), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;0.5 \cdot \sin re \leq -0.05:\\
\;\;\;\;\left(\left(\mathsf{fma}\left(re \cdot re, -0.16666666666666666, 1\right) \cdot re\right) \cdot \mathsf{fma}\left(-0.16666666666666666 \cdot im, im, -1\right)\right) \cdot im\\
\mathbf{else}:\\
\;\;\;\;\left(0.5 \cdot re\right) \cdot \left(\left(\left(\left(-0.016666666666666666 \cdot \left(im \cdot im\right) - 0.3333333333333333\right) \cdot im\right) \cdot im - 2\right) \cdot im\right)\\
\end{array}
\end{array}
if (*.f64 #s(literal 1/2 binary64) (sin.f64 re)) < -0.050000000000000003Initial program 51.3%
Taylor expanded in im around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
associate-*r*N/A
distribute-rgt-outN/A
lower-*.f64N/A
lift-sin.f64N/A
unpow2N/A
associate-*r*N/A
lower-fma.f64N/A
lower-*.f6486.3
Applied rewrites86.3%
Taylor expanded in re around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6428.3
Applied rewrites28.3%
if -0.050000000000000003 < (*.f64 #s(literal 1/2 binary64) (sin.f64 re)) Initial program 70.2%
Taylor expanded in im around 0
*-commutativeN/A
lower-*.f64N/A
lower--.f64N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
lower--.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6489.4
Applied rewrites89.4%
Taylor expanded in re around 0
Applied rewrites70.8%
(FPCore (re im)
:precision binary64
(let* ((t_0 (fma (* -0.16666666666666666 im) im -1.0)))
(if (<= (* 0.5 (sin re)) -0.05)
(* (* (* (fma (* re re) -0.16666666666666666 1.0) re) t_0) im)
(* re (* t_0 im)))))
double code(double re, double im) {
double t_0 = fma((-0.16666666666666666 * im), im, -1.0);
double tmp;
if ((0.5 * sin(re)) <= -0.05) {
tmp = ((fma((re * re), -0.16666666666666666, 1.0) * re) * t_0) * im;
} else {
tmp = re * (t_0 * im);
}
return tmp;
}
function code(re, im) t_0 = fma(Float64(-0.16666666666666666 * im), im, -1.0) tmp = 0.0 if (Float64(0.5 * sin(re)) <= -0.05) tmp = Float64(Float64(Float64(fma(Float64(re * re), -0.16666666666666666, 1.0) * re) * t_0) * im); else tmp = Float64(re * Float64(t_0 * im)); end return tmp end
code[re_, im_] := Block[{t$95$0 = N[(N[(-0.16666666666666666 * im), $MachinePrecision] * im + -1.0), $MachinePrecision]}, If[LessEqual[N[(0.5 * N[Sin[re], $MachinePrecision]), $MachinePrecision], -0.05], N[(N[(N[(N[(N[(re * re), $MachinePrecision] * -0.16666666666666666 + 1.0), $MachinePrecision] * re), $MachinePrecision] * t$95$0), $MachinePrecision] * im), $MachinePrecision], N[(re * N[(t$95$0 * im), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \mathsf{fma}\left(-0.16666666666666666 \cdot im, im, -1\right)\\
\mathbf{if}\;0.5 \cdot \sin re \leq -0.05:\\
\;\;\;\;\left(\left(\mathsf{fma}\left(re \cdot re, -0.16666666666666666, 1\right) \cdot re\right) \cdot t\_0\right) \cdot im\\
\mathbf{else}:\\
\;\;\;\;re \cdot \left(t\_0 \cdot im\right)\\
\end{array}
\end{array}
if (*.f64 #s(literal 1/2 binary64) (sin.f64 re)) < -0.050000000000000003Initial program 51.3%
Taylor expanded in im around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
associate-*r*N/A
distribute-rgt-outN/A
lower-*.f64N/A
lift-sin.f64N/A
unpow2N/A
associate-*r*N/A
lower-fma.f64N/A
lower-*.f6486.3
Applied rewrites86.3%
Taylor expanded in re around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6428.3
Applied rewrites28.3%
if -0.050000000000000003 < (*.f64 #s(literal 1/2 binary64) (sin.f64 re)) Initial program 70.2%
Taylor expanded in im around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
associate-*r*N/A
distribute-rgt-outN/A
lower-*.f64N/A
lift-sin.f64N/A
unpow2N/A
associate-*r*N/A
lower-fma.f64N/A
lower-*.f6480.4
Applied rewrites80.4%
Taylor expanded in re around 0
Applied rewrites62.9%
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-fma.f64N/A
associate-*l*N/A
lower-*.f64N/A
lower-*.f64N/A
lift-fma.f64N/A
lift-*.f6466.3
Applied rewrites66.3%
(FPCore (re im)
:precision binary64
(if (<= im -2.8e+140)
(*
(*
(* (fma (* re re) -0.16666666666666666 1.0) re)
(fma (* -0.16666666666666666 im) im -1.0))
im)
(if (<= im -9.8e-8)
(*
(*
(fma
(- (* 0.004166666666666667 (* re re)) 0.08333333333333333)
(* re re)
0.5)
re)
(*
(-
(*
(-
(*
(* (- (* -0.0003968253968253968 (* im im)) 0.016666666666666666) im)
im)
0.3333333333333333)
(* im im))
2.0)
im))
(if (<= im 1.06e+32)
(* (- (sin re)) im)
(* (* 0.5 re) (* (pow im 7.0) -0.0003968253968253968))))))
double code(double re, double im) {
double tmp;
if (im <= -2.8e+140) {
tmp = ((fma((re * re), -0.16666666666666666, 1.0) * re) * fma((-0.16666666666666666 * im), im, -1.0)) * im;
} else if (im <= -9.8e-8) {
tmp = (fma(((0.004166666666666667 * (re * re)) - 0.08333333333333333), (re * re), 0.5) * re) * ((((((((-0.0003968253968253968 * (im * im)) - 0.016666666666666666) * im) * im) - 0.3333333333333333) * (im * im)) - 2.0) * im);
} else if (im <= 1.06e+32) {
tmp = -sin(re) * im;
} else {
tmp = (0.5 * re) * (pow(im, 7.0) * -0.0003968253968253968);
}
return tmp;
}
function code(re, im) tmp = 0.0 if (im <= -2.8e+140) tmp = Float64(Float64(Float64(fma(Float64(re * re), -0.16666666666666666, 1.0) * re) * fma(Float64(-0.16666666666666666 * im), im, -1.0)) * im); elseif (im <= -9.8e-8) tmp = Float64(Float64(fma(Float64(Float64(0.004166666666666667 * Float64(re * re)) - 0.08333333333333333), Float64(re * re), 0.5) * re) * Float64(Float64(Float64(Float64(Float64(Float64(Float64(Float64(-0.0003968253968253968 * Float64(im * im)) - 0.016666666666666666) * im) * im) - 0.3333333333333333) * Float64(im * im)) - 2.0) * im)); elseif (im <= 1.06e+32) tmp = Float64(Float64(-sin(re)) * im); else tmp = Float64(Float64(0.5 * re) * Float64((im ^ 7.0) * -0.0003968253968253968)); end return tmp end
code[re_, im_] := If[LessEqual[im, -2.8e+140], N[(N[(N[(N[(N[(re * re), $MachinePrecision] * -0.16666666666666666 + 1.0), $MachinePrecision] * re), $MachinePrecision] * N[(N[(-0.16666666666666666 * im), $MachinePrecision] * im + -1.0), $MachinePrecision]), $MachinePrecision] * im), $MachinePrecision], If[LessEqual[im, -9.8e-8], N[(N[(N[(N[(N[(0.004166666666666667 * N[(re * re), $MachinePrecision]), $MachinePrecision] - 0.08333333333333333), $MachinePrecision] * N[(re * re), $MachinePrecision] + 0.5), $MachinePrecision] * re), $MachinePrecision] * N[(N[(N[(N[(N[(N[(N[(N[(-0.0003968253968253968 * N[(im * im), $MachinePrecision]), $MachinePrecision] - 0.016666666666666666), $MachinePrecision] * im), $MachinePrecision] * im), $MachinePrecision] - 0.3333333333333333), $MachinePrecision] * N[(im * im), $MachinePrecision]), $MachinePrecision] - 2.0), $MachinePrecision] * im), $MachinePrecision]), $MachinePrecision], If[LessEqual[im, 1.06e+32], N[((-N[Sin[re], $MachinePrecision]) * im), $MachinePrecision], N[(N[(0.5 * re), $MachinePrecision] * N[(N[Power[im, 7.0], $MachinePrecision] * -0.0003968253968253968), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;im \leq -2.8 \cdot 10^{+140}:\\
\;\;\;\;\left(\left(\mathsf{fma}\left(re \cdot re, -0.16666666666666666, 1\right) \cdot re\right) \cdot \mathsf{fma}\left(-0.16666666666666666 \cdot im, im, -1\right)\right) \cdot im\\
\mathbf{elif}\;im \leq -9.8 \cdot 10^{-8}:\\
\;\;\;\;\left(\mathsf{fma}\left(0.004166666666666667 \cdot \left(re \cdot re\right) - 0.08333333333333333, re \cdot re, 0.5\right) \cdot re\right) \cdot \left(\left(\left(\left(\left(-0.0003968253968253968 \cdot \left(im \cdot im\right) - 0.016666666666666666\right) \cdot im\right) \cdot im - 0.3333333333333333\right) \cdot \left(im \cdot im\right) - 2\right) \cdot im\right)\\
\mathbf{elif}\;im \leq 1.06 \cdot 10^{+32}:\\
\;\;\;\;\left(-\sin re\right) \cdot im\\
\mathbf{else}:\\
\;\;\;\;\left(0.5 \cdot re\right) \cdot \left({im}^{7} \cdot -0.0003968253968253968\right)\\
\end{array}
\end{array}
if im < -2.79999999999999983e140Initial program 100.0%
Taylor expanded in im around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
associate-*r*N/A
distribute-rgt-outN/A
lower-*.f64N/A
lift-sin.f64N/A
unpow2N/A
associate-*r*N/A
lower-fma.f64N/A
lower-*.f64100.0
Applied rewrites100.0%
Taylor expanded in re around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6485.7
Applied rewrites85.7%
if -2.79999999999999983e140 < im < -9.8000000000000004e-8Initial program 98.7%
Taylor expanded in im around 0
*-commutativeN/A
lower-*.f64N/A
Applied rewrites60.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
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6463.0
Applied rewrites63.0%
if -9.8000000000000004e-8 < im < 1.0600000000000001e32Initial program 36.8%
Taylor expanded in im around 0
*-commutativeN/A
associate-*r*N/A
lower-*.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lift-sin.f6496.7
Applied rewrites96.7%
if 1.0600000000000001e32 < im Initial program 100.0%
Taylor expanded in im around 0
*-commutativeN/A
lower-*.f64N/A
Applied rewrites94.3%
lift-*.f64N/A
lift-sin.f64N/A
*-commutativeN/A
lower-*.f64N/A
lift-sin.f6494.3
Applied rewrites94.3%
Taylor expanded in re around 0
*-commutativeN/A
lower-*.f6475.5
Applied rewrites75.5%
Taylor expanded in im around inf
*-commutativeN/A
lower-*.f64N/A
lower-pow.f6477.0
Applied rewrites77.0%
(FPCore (re im) :precision binary64 (if (<= (* 0.5 (sin re)) -0.05) (* (fma (* (* re re) im) 0.16666666666666666 (- im)) re) (* re (* (fma (* -0.16666666666666666 im) im -1.0) im))))
double code(double re, double im) {
double tmp;
if ((0.5 * sin(re)) <= -0.05) {
tmp = fma(((re * re) * im), 0.16666666666666666, -im) * re;
} else {
tmp = re * (fma((-0.16666666666666666 * im), im, -1.0) * im);
}
return tmp;
}
function code(re, im) tmp = 0.0 if (Float64(0.5 * sin(re)) <= -0.05) tmp = Float64(fma(Float64(Float64(re * re) * im), 0.16666666666666666, Float64(-im)) * re); else tmp = Float64(re * Float64(fma(Float64(-0.16666666666666666 * im), im, -1.0) * im)); end return tmp end
code[re_, im_] := If[LessEqual[N[(0.5 * N[Sin[re], $MachinePrecision]), $MachinePrecision], -0.05], N[(N[(N[(N[(re * re), $MachinePrecision] * im), $MachinePrecision] * 0.16666666666666666 + (-im)), $MachinePrecision] * re), $MachinePrecision], N[(re * N[(N[(N[(-0.16666666666666666 * im), $MachinePrecision] * im + -1.0), $MachinePrecision] * im), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;0.5 \cdot \sin re \leq -0.05:\\
\;\;\;\;\mathsf{fma}\left(\left(re \cdot re\right) \cdot im, 0.16666666666666666, -im\right) \cdot re\\
\mathbf{else}:\\
\;\;\;\;re \cdot \left(\mathsf{fma}\left(-0.16666666666666666 \cdot im, im, -1\right) \cdot im\right)\\
\end{array}
\end{array}
if (*.f64 #s(literal 1/2 binary64) (sin.f64 re)) < -0.050000000000000003Initial program 51.3%
Taylor expanded in im around 0
*-commutativeN/A
associate-*r*N/A
lower-*.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lift-sin.f6454.8
Applied rewrites54.8%
Taylor expanded in re around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
mul-1-negN/A
lower-neg.f6422.5
Applied rewrites22.5%
if -0.050000000000000003 < (*.f64 #s(literal 1/2 binary64) (sin.f64 re)) Initial program 70.2%
Taylor expanded in im around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
associate-*r*N/A
distribute-rgt-outN/A
lower-*.f64N/A
lift-sin.f64N/A
unpow2N/A
associate-*r*N/A
lower-fma.f64N/A
lower-*.f6480.4
Applied rewrites80.4%
Taylor expanded in re around 0
Applied rewrites62.9%
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-fma.f64N/A
associate-*l*N/A
lower-*.f64N/A
lower-*.f64N/A
lift-fma.f64N/A
lift-*.f6466.3
Applied rewrites66.3%
(FPCore (re im)
:precision binary64
(if (<= re 7.5e-21)
(*
(* 0.5 re)
(*
(-
(*
(*
(-
(*
(* (- (* (* -0.0003968253968253968 im) im) 0.016666666666666666) im)
im)
0.3333333333333333)
im)
im)
2.0)
im))
(* (* (sin re) im) (- (* (* im im) -0.16666666666666666) 1.0))))
double code(double re, double im) {
double tmp;
if (re <= 7.5e-21) {
tmp = (0.5 * re) * ((((((((((-0.0003968253968253968 * im) * im) - 0.016666666666666666) * im) * im) - 0.3333333333333333) * im) * im) - 2.0) * im);
} else {
tmp = (sin(re) * im) * (((im * im) * -0.16666666666666666) - 1.0);
}
return tmp;
}
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
real(8) :: tmp
if (re <= 7.5d-21) then
tmp = (0.5d0 * re) * (((((((((((-0.0003968253968253968d0) * im) * im) - 0.016666666666666666d0) * im) * im) - 0.3333333333333333d0) * im) * im) - 2.0d0) * im)
else
tmp = (sin(re) * im) * (((im * im) * (-0.16666666666666666d0)) - 1.0d0)
end if
code = tmp
end function
public static double code(double re, double im) {
double tmp;
if (re <= 7.5e-21) {
tmp = (0.5 * re) * ((((((((((-0.0003968253968253968 * im) * im) - 0.016666666666666666) * im) * im) - 0.3333333333333333) * im) * im) - 2.0) * im);
} else {
tmp = (Math.sin(re) * im) * (((im * im) * -0.16666666666666666) - 1.0);
}
return tmp;
}
def code(re, im): tmp = 0 if re <= 7.5e-21: tmp = (0.5 * re) * ((((((((((-0.0003968253968253968 * im) * im) - 0.016666666666666666) * im) * im) - 0.3333333333333333) * im) * im) - 2.0) * im) else: tmp = (math.sin(re) * im) * (((im * im) * -0.16666666666666666) - 1.0) return tmp
function code(re, im) tmp = 0.0 if (re <= 7.5e-21) tmp = Float64(Float64(0.5 * re) * Float64(Float64(Float64(Float64(Float64(Float64(Float64(Float64(Float64(Float64(-0.0003968253968253968 * im) * im) - 0.016666666666666666) * im) * im) - 0.3333333333333333) * im) * im) - 2.0) * im)); else tmp = Float64(Float64(sin(re) * im) * Float64(Float64(Float64(im * im) * -0.16666666666666666) - 1.0)); end return tmp end
function tmp_2 = code(re, im) tmp = 0.0; if (re <= 7.5e-21) tmp = (0.5 * re) * ((((((((((-0.0003968253968253968 * im) * im) - 0.016666666666666666) * im) * im) - 0.3333333333333333) * im) * im) - 2.0) * im); else tmp = (sin(re) * im) * (((im * im) * -0.16666666666666666) - 1.0); end tmp_2 = tmp; end
code[re_, im_] := If[LessEqual[re, 7.5e-21], N[(N[(0.5 * re), $MachinePrecision] * N[(N[(N[(N[(N[(N[(N[(N[(N[(N[(-0.0003968253968253968 * im), $MachinePrecision] * im), $MachinePrecision] - 0.016666666666666666), $MachinePrecision] * im), $MachinePrecision] * im), $MachinePrecision] - 0.3333333333333333), $MachinePrecision] * im), $MachinePrecision] * im), $MachinePrecision] - 2.0), $MachinePrecision] * im), $MachinePrecision]), $MachinePrecision], N[(N[(N[Sin[re], $MachinePrecision] * im), $MachinePrecision] * N[(N[(N[(im * im), $MachinePrecision] * -0.16666666666666666), $MachinePrecision] - 1.0), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;re \leq 7.5 \cdot 10^{-21}:\\
\;\;\;\;\left(0.5 \cdot re\right) \cdot \left(\left(\left(\left(\left(\left(\left(-0.0003968253968253968 \cdot im\right) \cdot im - 0.016666666666666666\right) \cdot im\right) \cdot im - 0.3333333333333333\right) \cdot im\right) \cdot im - 2\right) \cdot im\right)\\
\mathbf{else}:\\
\;\;\;\;\left(\sin re \cdot im\right) \cdot \left(\left(im \cdot im\right) \cdot -0.16666666666666666 - 1\right)\\
\end{array}
\end{array}
if re < 7.50000000000000072e-21Initial program 69.4%
Taylor expanded in im around 0
*-commutativeN/A
lower-*.f64N/A
Applied rewrites94.5%
lift-*.f64N/A
lift-sin.f64N/A
*-commutativeN/A
lower-*.f64N/A
lift-sin.f6494.5
Applied rewrites94.5%
Taylor expanded in re around 0
*-commutativeN/A
lower-*.f6474.0
Applied rewrites74.0%
lift-*.f64N/A
lift-*.f64N/A
lift--.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift--.f64N/A
lift-*.f64N/A
lift-*.f64N/A
associate-*r*N/A
lower-*.f64N/A
Applied rewrites74.0%
if 7.50000000000000072e-21 < re Initial program 54.3%
Taylor expanded in im around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
associate-*r*N/A
distribute-rgt-outN/A
lower-*.f64N/A
lift-sin.f64N/A
unpow2N/A
associate-*r*N/A
lower-fma.f64N/A
lower-*.f6486.6
Applied rewrites86.6%
Taylor expanded in re around inf
associate-*r*N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lift-sin.f64N/A
lower--.f64N/A
*-commutativeN/A
lower-*.f64N/A
pow2N/A
lift-*.f6486.6
Applied rewrites86.6%
(FPCore (re im)
:precision binary64
(if (<= re 7.5e-21)
(*
(* 0.5 re)
(*
(-
(*
(*
(-
(*
(* (- (* (* -0.0003968253968253968 im) im) 0.016666666666666666) im)
im)
0.3333333333333333)
im)
im)
2.0)
im))
(* (* (sin re) (fma (* -0.16666666666666666 im) im -1.0)) im)))
double code(double re, double im) {
double tmp;
if (re <= 7.5e-21) {
tmp = (0.5 * re) * ((((((((((-0.0003968253968253968 * im) * im) - 0.016666666666666666) * im) * im) - 0.3333333333333333) * im) * im) - 2.0) * im);
} else {
tmp = (sin(re) * fma((-0.16666666666666666 * im), im, -1.0)) * im;
}
return tmp;
}
function code(re, im) tmp = 0.0 if (re <= 7.5e-21) tmp = Float64(Float64(0.5 * re) * Float64(Float64(Float64(Float64(Float64(Float64(Float64(Float64(Float64(Float64(-0.0003968253968253968 * im) * im) - 0.016666666666666666) * im) * im) - 0.3333333333333333) * im) * im) - 2.0) * im)); else tmp = Float64(Float64(sin(re) * fma(Float64(-0.16666666666666666 * im), im, -1.0)) * im); end return tmp end
code[re_, im_] := If[LessEqual[re, 7.5e-21], N[(N[(0.5 * re), $MachinePrecision] * N[(N[(N[(N[(N[(N[(N[(N[(N[(N[(-0.0003968253968253968 * im), $MachinePrecision] * im), $MachinePrecision] - 0.016666666666666666), $MachinePrecision] * im), $MachinePrecision] * im), $MachinePrecision] - 0.3333333333333333), $MachinePrecision] * im), $MachinePrecision] * im), $MachinePrecision] - 2.0), $MachinePrecision] * im), $MachinePrecision]), $MachinePrecision], N[(N[(N[Sin[re], $MachinePrecision] * N[(N[(-0.16666666666666666 * im), $MachinePrecision] * im + -1.0), $MachinePrecision]), $MachinePrecision] * im), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;re \leq 7.5 \cdot 10^{-21}:\\
\;\;\;\;\left(0.5 \cdot re\right) \cdot \left(\left(\left(\left(\left(\left(\left(-0.0003968253968253968 \cdot im\right) \cdot im - 0.016666666666666666\right) \cdot im\right) \cdot im - 0.3333333333333333\right) \cdot im\right) \cdot im - 2\right) \cdot im\right)\\
\mathbf{else}:\\
\;\;\;\;\left(\sin re \cdot \mathsf{fma}\left(-0.16666666666666666 \cdot im, im, -1\right)\right) \cdot im\\
\end{array}
\end{array}
if re < 7.50000000000000072e-21Initial program 69.4%
Taylor expanded in im around 0
*-commutativeN/A
lower-*.f64N/A
Applied rewrites94.5%
lift-*.f64N/A
lift-sin.f64N/A
*-commutativeN/A
lower-*.f64N/A
lift-sin.f6494.5
Applied rewrites94.5%
Taylor expanded in re around 0
*-commutativeN/A
lower-*.f6474.0
Applied rewrites74.0%
lift-*.f64N/A
lift-*.f64N/A
lift--.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift--.f64N/A
lift-*.f64N/A
lift-*.f64N/A
associate-*r*N/A
lower-*.f64N/A
Applied rewrites74.0%
if 7.50000000000000072e-21 < re Initial program 54.3%
Taylor expanded in im around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
associate-*r*N/A
distribute-rgt-outN/A
lower-*.f64N/A
lift-sin.f64N/A
unpow2N/A
associate-*r*N/A
lower-fma.f64N/A
lower-*.f6486.6
Applied rewrites86.6%
(FPCore (re im)
:precision binary64
(if (<= im -2.8e+140)
(*
(*
(* (fma (* re re) -0.16666666666666666 1.0) re)
(fma (* -0.16666666666666666 im) im -1.0))
im)
(if (or (<= im -9.8e-8) (not (<= im 1.06e+32)))
(*
(*
(fma
(- (* 0.004166666666666667 (* re re)) 0.08333333333333333)
(* re re)
0.5)
re)
(*
(-
(*
(-
(*
(* (- (* -0.0003968253968253968 (* im im)) 0.016666666666666666) im)
im)
0.3333333333333333)
(* im im))
2.0)
im))
(* (- (sin re)) im))))
double code(double re, double im) {
double tmp;
if (im <= -2.8e+140) {
tmp = ((fma((re * re), -0.16666666666666666, 1.0) * re) * fma((-0.16666666666666666 * im), im, -1.0)) * im;
} else if ((im <= -9.8e-8) || !(im <= 1.06e+32)) {
tmp = (fma(((0.004166666666666667 * (re * re)) - 0.08333333333333333), (re * re), 0.5) * re) * ((((((((-0.0003968253968253968 * (im * im)) - 0.016666666666666666) * im) * im) - 0.3333333333333333) * (im * im)) - 2.0) * im);
} else {
tmp = -sin(re) * im;
}
return tmp;
}
function code(re, im) tmp = 0.0 if (im <= -2.8e+140) tmp = Float64(Float64(Float64(fma(Float64(re * re), -0.16666666666666666, 1.0) * re) * fma(Float64(-0.16666666666666666 * im), im, -1.0)) * im); elseif ((im <= -9.8e-8) || !(im <= 1.06e+32)) tmp = Float64(Float64(fma(Float64(Float64(0.004166666666666667 * Float64(re * re)) - 0.08333333333333333), Float64(re * re), 0.5) * re) * Float64(Float64(Float64(Float64(Float64(Float64(Float64(Float64(-0.0003968253968253968 * Float64(im * im)) - 0.016666666666666666) * im) * im) - 0.3333333333333333) * Float64(im * im)) - 2.0) * im)); else tmp = Float64(Float64(-sin(re)) * im); end return tmp end
code[re_, im_] := If[LessEqual[im, -2.8e+140], N[(N[(N[(N[(N[(re * re), $MachinePrecision] * -0.16666666666666666 + 1.0), $MachinePrecision] * re), $MachinePrecision] * N[(N[(-0.16666666666666666 * im), $MachinePrecision] * im + -1.0), $MachinePrecision]), $MachinePrecision] * im), $MachinePrecision], If[Or[LessEqual[im, -9.8e-8], N[Not[LessEqual[im, 1.06e+32]], $MachinePrecision]], N[(N[(N[(N[(N[(0.004166666666666667 * N[(re * re), $MachinePrecision]), $MachinePrecision] - 0.08333333333333333), $MachinePrecision] * N[(re * re), $MachinePrecision] + 0.5), $MachinePrecision] * re), $MachinePrecision] * N[(N[(N[(N[(N[(N[(N[(N[(-0.0003968253968253968 * N[(im * im), $MachinePrecision]), $MachinePrecision] - 0.016666666666666666), $MachinePrecision] * im), $MachinePrecision] * im), $MachinePrecision] - 0.3333333333333333), $MachinePrecision] * N[(im * im), $MachinePrecision]), $MachinePrecision] - 2.0), $MachinePrecision] * im), $MachinePrecision]), $MachinePrecision], N[((-N[Sin[re], $MachinePrecision]) * im), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;im \leq -2.8 \cdot 10^{+140}:\\
\;\;\;\;\left(\left(\mathsf{fma}\left(re \cdot re, -0.16666666666666666, 1\right) \cdot re\right) \cdot \mathsf{fma}\left(-0.16666666666666666 \cdot im, im, -1\right)\right) \cdot im\\
\mathbf{elif}\;im \leq -9.8 \cdot 10^{-8} \lor \neg \left(im \leq 1.06 \cdot 10^{+32}\right):\\
\;\;\;\;\left(\mathsf{fma}\left(0.004166666666666667 \cdot \left(re \cdot re\right) - 0.08333333333333333, re \cdot re, 0.5\right) \cdot re\right) \cdot \left(\left(\left(\left(\left(-0.0003968253968253968 \cdot \left(im \cdot im\right) - 0.016666666666666666\right) \cdot im\right) \cdot im - 0.3333333333333333\right) \cdot \left(im \cdot im\right) - 2\right) \cdot im\right)\\
\mathbf{else}:\\
\;\;\;\;\left(-\sin re\right) \cdot im\\
\end{array}
\end{array}
if im < -2.79999999999999983e140Initial program 100.0%
Taylor expanded in im around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
associate-*r*N/A
distribute-rgt-outN/A
lower-*.f64N/A
lift-sin.f64N/A
unpow2N/A
associate-*r*N/A
lower-fma.f64N/A
lower-*.f64100.0
Applied rewrites100.0%
Taylor expanded in re around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6485.7
Applied rewrites85.7%
if -2.79999999999999983e140 < im < -9.8000000000000004e-8 or 1.0600000000000001e32 < im Initial program 99.7%
Taylor expanded in im around 0
*-commutativeN/A
lower-*.f64N/A
Applied rewrites85.1%
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
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6472.2
Applied rewrites72.2%
if -9.8000000000000004e-8 < im < 1.0600000000000001e32Initial program 36.8%
Taylor expanded in im around 0
*-commutativeN/A
associate-*r*N/A
lower-*.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lift-sin.f6496.7
Applied rewrites96.7%
Final simplification87.0%
(FPCore (re im) :precision binary64 (if (or (<= im -2.45) (not (<= im 2.05e-9))) (* (* re (* (* im im) -0.16666666666666666)) im) (* (- re) im)))
double code(double re, double im) {
double tmp;
if ((im <= -2.45) || !(im <= 2.05e-9)) {
tmp = (re * ((im * im) * -0.16666666666666666)) * im;
} else {
tmp = -re * im;
}
return tmp;
}
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
real(8) :: tmp
if ((im <= (-2.45d0)) .or. (.not. (im <= 2.05d-9))) then
tmp = (re * ((im * im) * (-0.16666666666666666d0))) * im
else
tmp = -re * im
end if
code = tmp
end function
public static double code(double re, double im) {
double tmp;
if ((im <= -2.45) || !(im <= 2.05e-9)) {
tmp = (re * ((im * im) * -0.16666666666666666)) * im;
} else {
tmp = -re * im;
}
return tmp;
}
def code(re, im): tmp = 0 if (im <= -2.45) or not (im <= 2.05e-9): tmp = (re * ((im * im) * -0.16666666666666666)) * im else: tmp = -re * im return tmp
function code(re, im) tmp = 0.0 if ((im <= -2.45) || !(im <= 2.05e-9)) tmp = Float64(Float64(re * Float64(Float64(im * im) * -0.16666666666666666)) * im); else tmp = Float64(Float64(-re) * im); end return tmp end
function tmp_2 = code(re, im) tmp = 0.0; if ((im <= -2.45) || ~((im <= 2.05e-9))) tmp = (re * ((im * im) * -0.16666666666666666)) * im; else tmp = -re * im; end tmp_2 = tmp; end
code[re_, im_] := If[Or[LessEqual[im, -2.45], N[Not[LessEqual[im, 2.05e-9]], $MachinePrecision]], N[(N[(re * N[(N[(im * im), $MachinePrecision] * -0.16666666666666666), $MachinePrecision]), $MachinePrecision] * im), $MachinePrecision], N[((-re) * im), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;im \leq -2.45 \lor \neg \left(im \leq 2.05 \cdot 10^{-9}\right):\\
\;\;\;\;\left(re \cdot \left(\left(im \cdot im\right) \cdot -0.16666666666666666\right)\right) \cdot im\\
\mathbf{else}:\\
\;\;\;\;\left(-re\right) \cdot im\\
\end{array}
\end{array}
if im < -2.4500000000000002 or 2.0500000000000002e-9 < im Initial program 99.6%
Taylor expanded in im around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
associate-*r*N/A
distribute-rgt-outN/A
lower-*.f64N/A
lift-sin.f64N/A
unpow2N/A
associate-*r*N/A
lower-fma.f64N/A
lower-*.f6462.0
Applied rewrites62.0%
Taylor expanded in re around 0
Applied rewrites46.4%
Taylor expanded in im around inf
*-commutativeN/A
pow2N/A
lift-*.f64N/A
lift-*.f6446.4
Applied rewrites46.4%
if -2.4500000000000002 < im < 2.0500000000000002e-9Initial program 34.6%
Taylor expanded in im around 0
*-commutativeN/A
associate-*r*N/A
lower-*.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lift-sin.f6499.8
Applied rewrites99.8%
Taylor expanded in re around 0
Applied rewrites58.2%
Final simplification52.6%
(FPCore (re im) :precision binary64 (* re (* (fma (* -0.16666666666666666 im) im -1.0) im)))
double code(double re, double im) {
return re * (fma((-0.16666666666666666 * im), im, -1.0) * im);
}
function code(re, im) return Float64(re * Float64(fma(Float64(-0.16666666666666666 * im), im, -1.0) * im)) end
code[re_, im_] := N[(re * N[(N[(N[(-0.16666666666666666 * im), $MachinePrecision] * im + -1.0), $MachinePrecision] * im), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
re \cdot \left(\mathsf{fma}\left(-0.16666666666666666 \cdot im, im, -1\right) \cdot im\right)
\end{array}
Initial program 65.6%
Taylor expanded in im around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
associate-*r*N/A
distribute-rgt-outN/A
lower-*.f64N/A
lift-sin.f64N/A
unpow2N/A
associate-*r*N/A
lower-fma.f64N/A
lower-*.f6481.8
Applied rewrites81.8%
Taylor expanded in re around 0
Applied rewrites52.6%
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-fma.f64N/A
associate-*l*N/A
lower-*.f64N/A
lower-*.f64N/A
lift-fma.f64N/A
lift-*.f6455.2
Applied rewrites55.2%
(FPCore (re im) :precision binary64 (* (- re) im))
double code(double re, double im) {
return -re * 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 = -re * im
end function
public static double code(double re, double im) {
return -re * im;
}
def code(re, im): return -re * im
function code(re, im) return Float64(Float64(-re) * im) end
function tmp = code(re, im) tmp = -re * im; end
code[re_, im_] := N[((-re) * im), $MachinePrecision]
\begin{array}{l}
\\
\left(-re\right) \cdot im
\end{array}
Initial program 65.6%
Taylor expanded in im around 0
*-commutativeN/A
associate-*r*N/A
lower-*.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lift-sin.f6454.8
Applied rewrites54.8%
Taylor expanded in re around 0
Applied rewrites38.2%
(FPCore (re im)
:precision binary64
(if (< (fabs im) 1.0)
(-
(*
(sin re)
(+
(+ im (* (* (* 0.16666666666666666 im) im) im))
(* (* (* (* (* 0.008333333333333333 im) im) im) im) im))))
(* (* 0.5 (sin re)) (- (exp (- im)) (exp im)))))
double code(double re, double im) {
double tmp;
if (fabs(im) < 1.0) {
tmp = -(sin(re) * ((im + (((0.16666666666666666 * im) * im) * im)) + (((((0.008333333333333333 * im) * im) * im) * im) * im)));
} else {
tmp = (0.5 * sin(re)) * (exp(-im) - exp(im));
}
return tmp;
}
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
real(8) :: tmp
if (abs(im) < 1.0d0) then
tmp = -(sin(re) * ((im + (((0.16666666666666666d0 * im) * im) * im)) + (((((0.008333333333333333d0 * im) * im) * im) * im) * im)))
else
tmp = (0.5d0 * sin(re)) * (exp(-im) - exp(im))
end if
code = tmp
end function
public static double code(double re, double im) {
double tmp;
if (Math.abs(im) < 1.0) {
tmp = -(Math.sin(re) * ((im + (((0.16666666666666666 * im) * im) * im)) + (((((0.008333333333333333 * im) * im) * im) * im) * im)));
} else {
tmp = (0.5 * Math.sin(re)) * (Math.exp(-im) - Math.exp(im));
}
return tmp;
}
def code(re, im): tmp = 0 if math.fabs(im) < 1.0: tmp = -(math.sin(re) * ((im + (((0.16666666666666666 * im) * im) * im)) + (((((0.008333333333333333 * im) * im) * im) * im) * im))) else: tmp = (0.5 * math.sin(re)) * (math.exp(-im) - math.exp(im)) return tmp
function code(re, im) tmp = 0.0 if (abs(im) < 1.0) tmp = Float64(-Float64(sin(re) * Float64(Float64(im + Float64(Float64(Float64(0.16666666666666666 * im) * im) * im)) + Float64(Float64(Float64(Float64(Float64(0.008333333333333333 * im) * im) * im) * im) * im)))); else tmp = Float64(Float64(0.5 * sin(re)) * Float64(exp(Float64(-im)) - exp(im))); end return tmp end
function tmp_2 = code(re, im) tmp = 0.0; if (abs(im) < 1.0) tmp = -(sin(re) * ((im + (((0.16666666666666666 * im) * im) * im)) + (((((0.008333333333333333 * im) * im) * im) * im) * im))); else tmp = (0.5 * sin(re)) * (exp(-im) - exp(im)); end tmp_2 = tmp; end
code[re_, im_] := If[Less[N[Abs[im], $MachinePrecision], 1.0], (-N[(N[Sin[re], $MachinePrecision] * N[(N[(im + N[(N[(N[(0.16666666666666666 * im), $MachinePrecision] * im), $MachinePrecision] * im), $MachinePrecision]), $MachinePrecision] + N[(N[(N[(N[(N[(0.008333333333333333 * im), $MachinePrecision] * im), $MachinePrecision] * im), $MachinePrecision] * im), $MachinePrecision] * im), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), N[(N[(0.5 * N[Sin[re], $MachinePrecision]), $MachinePrecision] * N[(N[Exp[(-im)], $MachinePrecision] - N[Exp[im], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\left|im\right| < 1:\\
\;\;\;\;-\sin re \cdot \left(\left(im + \left(\left(0.16666666666666666 \cdot im\right) \cdot im\right) \cdot im\right) + \left(\left(\left(\left(0.008333333333333333 \cdot im\right) \cdot im\right) \cdot im\right) \cdot im\right) \cdot im\right)\\
\mathbf{else}:\\
\;\;\;\;\left(0.5 \cdot \sin re\right) \cdot \left(e^{-im} - e^{im}\right)\\
\end{array}
\end{array}
herbie shell --seed 2025037
(FPCore (re im)
:name "math.cos on complex, imaginary part"
:precision binary64
:alt
(! :herbie-platform default (if (< (fabs im) 1) (- (* (sin re) (+ im (* 1/6 im im im) (* 1/120 im im im im im)))) (* (* 1/2 (sin re)) (- (exp (- im)) (exp im)))))
(* (* 0.5 (sin re)) (- (exp (- im)) (exp im))))