
(FPCore (t l k) :precision binary64 (/ 2.0 (* (* (* (/ (pow t 3.0) (* l l)) (sin k)) (tan k)) (+ (+ 1.0 (pow (/ k t) 2.0)) 1.0))))
double code(double t, double l, double k) {
return 2.0 / ((((pow(t, 3.0) / (l * l)) * sin(k)) * tan(k)) * ((1.0 + pow((k / t), 2.0)) + 1.0));
}
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(t, l, k)
use fmin_fmax_functions
real(8), intent (in) :: t
real(8), intent (in) :: l
real(8), intent (in) :: k
code = 2.0d0 / (((((t ** 3.0d0) / (l * l)) * sin(k)) * tan(k)) * ((1.0d0 + ((k / t) ** 2.0d0)) + 1.0d0))
end function
public static double code(double t, double l, double k) {
return 2.0 / ((((Math.pow(t, 3.0) / (l * l)) * Math.sin(k)) * Math.tan(k)) * ((1.0 + Math.pow((k / t), 2.0)) + 1.0));
}
def code(t, l, k): return 2.0 / ((((math.pow(t, 3.0) / (l * l)) * math.sin(k)) * math.tan(k)) * ((1.0 + math.pow((k / t), 2.0)) + 1.0))
function code(t, l, k) return Float64(2.0 / Float64(Float64(Float64(Float64((t ^ 3.0) / Float64(l * l)) * sin(k)) * tan(k)) * Float64(Float64(1.0 + (Float64(k / t) ^ 2.0)) + 1.0))) end
function tmp = code(t, l, k) tmp = 2.0 / (((((t ^ 3.0) / (l * l)) * sin(k)) * tan(k)) * ((1.0 + ((k / t) ^ 2.0)) + 1.0)); end
code[t_, l_, k_] := N[(2.0 / N[(N[(N[(N[(N[Power[t, 3.0], $MachinePrecision] / N[(l * l), $MachinePrecision]), $MachinePrecision] * N[Sin[k], $MachinePrecision]), $MachinePrecision] * N[Tan[k], $MachinePrecision]), $MachinePrecision] * N[(N[(1.0 + N[Power[N[(k / t), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{2}{\left(\left(\frac{{t}^{3}}{\ell \cdot \ell} \cdot \sin k\right) \cdot \tan k\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 21 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (t l k) :precision binary64 (/ 2.0 (* (* (* (/ (pow t 3.0) (* l l)) (sin k)) (tan k)) (+ (+ 1.0 (pow (/ k t) 2.0)) 1.0))))
double code(double t, double l, double k) {
return 2.0 / ((((pow(t, 3.0) / (l * l)) * sin(k)) * tan(k)) * ((1.0 + pow((k / t), 2.0)) + 1.0));
}
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(t, l, k)
use fmin_fmax_functions
real(8), intent (in) :: t
real(8), intent (in) :: l
real(8), intent (in) :: k
code = 2.0d0 / (((((t ** 3.0d0) / (l * l)) * sin(k)) * tan(k)) * ((1.0d0 + ((k / t) ** 2.0d0)) + 1.0d0))
end function
public static double code(double t, double l, double k) {
return 2.0 / ((((Math.pow(t, 3.0) / (l * l)) * Math.sin(k)) * Math.tan(k)) * ((1.0 + Math.pow((k / t), 2.0)) + 1.0));
}
def code(t, l, k): return 2.0 / ((((math.pow(t, 3.0) / (l * l)) * math.sin(k)) * math.tan(k)) * ((1.0 + math.pow((k / t), 2.0)) + 1.0))
function code(t, l, k) return Float64(2.0 / Float64(Float64(Float64(Float64((t ^ 3.0) / Float64(l * l)) * sin(k)) * tan(k)) * Float64(Float64(1.0 + (Float64(k / t) ^ 2.0)) + 1.0))) end
function tmp = code(t, l, k) tmp = 2.0 / (((((t ^ 3.0) / (l * l)) * sin(k)) * tan(k)) * ((1.0 + ((k / t) ^ 2.0)) + 1.0)); end
code[t_, l_, k_] := N[(2.0 / N[(N[(N[(N[(N[Power[t, 3.0], $MachinePrecision] / N[(l * l), $MachinePrecision]), $MachinePrecision] * N[Sin[k], $MachinePrecision]), $MachinePrecision] * N[Tan[k], $MachinePrecision]), $MachinePrecision] * N[(N[(1.0 + N[Power[N[(k / t), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{2}{\left(\left(\frac{{t}^{3}}{\ell \cdot \ell} \cdot \sin k\right) \cdot \tan k\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)}
\end{array}
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k)
:precision binary64
(*
t_s
(if (<= t_m 1.02e-82)
(/ 2.0 (/ (/ (* (* (* (pow (sin k) 2.0) t_m) k) k) (* (cos k) l)) l))
(/
2.0
(*
(* (* (/ t_m l) (sin k)) t_m)
(* (/ t_m l) (* (tan k) (+ (pow (/ k t_m) 2.0) 2.0))))))))t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
double tmp;
if (t_m <= 1.02e-82) {
tmp = 2.0 / (((((pow(sin(k), 2.0) * t_m) * k) * k) / (cos(k) * l)) / l);
} else {
tmp = 2.0 / ((((t_m / l) * sin(k)) * t_m) * ((t_m / l) * (tan(k) * (pow((k / t_m), 2.0) + 2.0))));
}
return t_s * tmp;
}
t\_m = private
t\_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(t_s, t_m, l, k)
use fmin_fmax_functions
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k
real(8) :: tmp
if (t_m <= 1.02d-82) then
tmp = 2.0d0 / ((((((sin(k) ** 2.0d0) * t_m) * k) * k) / (cos(k) * l)) / l)
else
tmp = 2.0d0 / ((((t_m / l) * sin(k)) * t_m) * ((t_m / l) * (tan(k) * (((k / t_m) ** 2.0d0) + 2.0d0))))
end if
code = t_s * tmp
end function
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k) {
double tmp;
if (t_m <= 1.02e-82) {
tmp = 2.0 / (((((Math.pow(Math.sin(k), 2.0) * t_m) * k) * k) / (Math.cos(k) * l)) / l);
} else {
tmp = 2.0 / ((((t_m / l) * Math.sin(k)) * t_m) * ((t_m / l) * (Math.tan(k) * (Math.pow((k / t_m), 2.0) + 2.0))));
}
return t_s * tmp;
}
t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k): tmp = 0 if t_m <= 1.02e-82: tmp = 2.0 / (((((math.pow(math.sin(k), 2.0) * t_m) * k) * k) / (math.cos(k) * l)) / l) else: tmp = 2.0 / ((((t_m / l) * math.sin(k)) * t_m) * ((t_m / l) * (math.tan(k) * (math.pow((k / t_m), 2.0) + 2.0)))) return t_s * tmp
t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k) tmp = 0.0 if (t_m <= 1.02e-82) tmp = Float64(2.0 / Float64(Float64(Float64(Float64(Float64((sin(k) ^ 2.0) * t_m) * k) * k) / Float64(cos(k) * l)) / l)); else tmp = Float64(2.0 / Float64(Float64(Float64(Float64(t_m / l) * sin(k)) * t_m) * Float64(Float64(t_m / l) * Float64(tan(k) * Float64((Float64(k / t_m) ^ 2.0) + 2.0))))); end return Float64(t_s * tmp) end
t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k) tmp = 0.0; if (t_m <= 1.02e-82) tmp = 2.0 / ((((((sin(k) ^ 2.0) * t_m) * k) * k) / (cos(k) * l)) / l); else tmp = 2.0 / ((((t_m / l) * sin(k)) * t_m) * ((t_m / l) * (tan(k) * (((k / t_m) ^ 2.0) + 2.0)))); end tmp_2 = t_s * tmp; end
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k_] := N[(t$95$s * If[LessEqual[t$95$m, 1.02e-82], N[(2.0 / N[(N[(N[(N[(N[(N[Power[N[Sin[k], $MachinePrecision], 2.0], $MachinePrecision] * t$95$m), $MachinePrecision] * k), $MachinePrecision] * k), $MachinePrecision] / N[(N[Cos[k], $MachinePrecision] * l), $MachinePrecision]), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(N[(N[(t$95$m / l), $MachinePrecision] * N[Sin[k], $MachinePrecision]), $MachinePrecision] * t$95$m), $MachinePrecision] * N[(N[(t$95$m / l), $MachinePrecision] * N[(N[Tan[k], $MachinePrecision] * N[(N[Power[N[(k / t$95$m), $MachinePrecision], 2.0], $MachinePrecision] + 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_m \leq 1.02 \cdot 10^{-82}:\\
\;\;\;\;\frac{2}{\frac{\frac{\left(\left({\sin k}^{2} \cdot t\_m\right) \cdot k\right) \cdot k}{\cos k \cdot \ell}}{\ell}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\left(\left(\frac{t\_m}{\ell} \cdot \sin k\right) \cdot t\_m\right) \cdot \left(\frac{t\_m}{\ell} \cdot \left(\tan k \cdot \left({\left(\frac{k}{t\_m}\right)}^{2} + 2\right)\right)\right)}\\
\end{array}
\end{array}
if t < 1.02000000000000007e-82Initial program 54.0%
Applied rewrites48.0%
Taylor expanded in t around 0
lower-/.f64N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-pow.f64N/A
lower-sin.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-cos.f6477.8
Applied rewrites77.8%
if 1.02000000000000007e-82 < t Initial program 71.3%
lift-*.f64N/A
*-commutativeN/A
lift-/.f64N/A
lift-pow.f64N/A
cube-multN/A
associate-/l*N/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
lift-*.f64N/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
lower-/.f6487.2
Applied rewrites87.2%
lift-*.f64N/A
lift-*.f64N/A
associate-*l*N/A
lift-*.f64N/A
lift-*.f64N/A
associate-*r*N/A
associate-*l*N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lift-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f64N/A
Applied rewrites94.2%
lift-*.f64N/A
lift-*.f64N/A
*-commutativeN/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f6494.2
Applied rewrites94.2%
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k)
:precision binary64
(let* ((t_2 (* (* k k) t_m)))
(*
t_s
(if (<=
(/
2.0
(*
(* (* (/ (pow t_m 3.0) (* l l)) (sin k)) (tan k))
(+ (+ 1.0 (pow (/ k t_m) 2.0)) 1.0)))
1e+188)
(/ (* l l) (* (* t_m t_m) t_2))
(/ (* (* l l) 2.0) (* t_2 (* k k)))))))t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
double t_2 = (k * k) * t_m;
double tmp;
if ((2.0 / ((((pow(t_m, 3.0) / (l * l)) * sin(k)) * tan(k)) * ((1.0 + pow((k / t_m), 2.0)) + 1.0))) <= 1e+188) {
tmp = (l * l) / ((t_m * t_m) * t_2);
} else {
tmp = ((l * l) * 2.0) / (t_2 * (k * k));
}
return t_s * tmp;
}
t\_m = private
t\_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(t_s, t_m, l, k)
use fmin_fmax_functions
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k
real(8) :: t_2
real(8) :: tmp
t_2 = (k * k) * t_m
if ((2.0d0 / (((((t_m ** 3.0d0) / (l * l)) * sin(k)) * tan(k)) * ((1.0d0 + ((k / t_m) ** 2.0d0)) + 1.0d0))) <= 1d+188) then
tmp = (l * l) / ((t_m * t_m) * t_2)
else
tmp = ((l * l) * 2.0d0) / (t_2 * (k * k))
end if
code = t_s * tmp
end function
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k) {
double t_2 = (k * k) * t_m;
double tmp;
if ((2.0 / ((((Math.pow(t_m, 3.0) / (l * l)) * Math.sin(k)) * Math.tan(k)) * ((1.0 + Math.pow((k / t_m), 2.0)) + 1.0))) <= 1e+188) {
tmp = (l * l) / ((t_m * t_m) * t_2);
} else {
tmp = ((l * l) * 2.0) / (t_2 * (k * k));
}
return t_s * tmp;
}
t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k): t_2 = (k * k) * t_m tmp = 0 if (2.0 / ((((math.pow(t_m, 3.0) / (l * l)) * math.sin(k)) * math.tan(k)) * ((1.0 + math.pow((k / t_m), 2.0)) + 1.0))) <= 1e+188: tmp = (l * l) / ((t_m * t_m) * t_2) else: tmp = ((l * l) * 2.0) / (t_2 * (k * k)) return t_s * tmp
t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k) t_2 = Float64(Float64(k * k) * t_m) tmp = 0.0 if (Float64(2.0 / Float64(Float64(Float64(Float64((t_m ^ 3.0) / Float64(l * l)) * sin(k)) * tan(k)) * Float64(Float64(1.0 + (Float64(k / t_m) ^ 2.0)) + 1.0))) <= 1e+188) tmp = Float64(Float64(l * l) / Float64(Float64(t_m * t_m) * t_2)); else tmp = Float64(Float64(Float64(l * l) * 2.0) / Float64(t_2 * Float64(k * k))); end return Float64(t_s * tmp) end
t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k) t_2 = (k * k) * t_m; tmp = 0.0; if ((2.0 / (((((t_m ^ 3.0) / (l * l)) * sin(k)) * tan(k)) * ((1.0 + ((k / t_m) ^ 2.0)) + 1.0))) <= 1e+188) tmp = (l * l) / ((t_m * t_m) * t_2); else tmp = ((l * l) * 2.0) / (t_2 * (k * k)); end tmp_2 = t_s * tmp; end
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k_] := Block[{t$95$2 = N[(N[(k * k), $MachinePrecision] * t$95$m), $MachinePrecision]}, N[(t$95$s * If[LessEqual[N[(2.0 / N[(N[(N[(N[(N[Power[t$95$m, 3.0], $MachinePrecision] / N[(l * l), $MachinePrecision]), $MachinePrecision] * N[Sin[k], $MachinePrecision]), $MachinePrecision] * N[Tan[k], $MachinePrecision]), $MachinePrecision] * N[(N[(1.0 + N[Power[N[(k / t$95$m), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 1e+188], N[(N[(l * l), $MachinePrecision] / N[(N[(t$95$m * t$95$m), $MachinePrecision] * t$95$2), $MachinePrecision]), $MachinePrecision], N[(N[(N[(l * l), $MachinePrecision] * 2.0), $MachinePrecision] / N[(t$95$2 * N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]]
\begin{array}{l}
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
\begin{array}{l}
t_2 := \left(k \cdot k\right) \cdot t\_m\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;\frac{2}{\left(\left(\frac{{t\_m}^{3}}{\ell \cdot \ell} \cdot \sin k\right) \cdot \tan k\right) \cdot \left(\left(1 + {\left(\frac{k}{t\_m}\right)}^{2}\right) + 1\right)} \leq 10^{+188}:\\
\;\;\;\;\frac{\ell \cdot \ell}{\left(t\_m \cdot t\_m\right) \cdot t\_2}\\
\mathbf{else}:\\
\;\;\;\;\frac{\left(\ell \cdot \ell\right) \cdot 2}{t\_2 \cdot \left(k \cdot k\right)}\\
\end{array}
\end{array}
\end{array}
if (/.f64 #s(literal 2 binary64) (*.f64 (*.f64 (*.f64 (/.f64 (pow.f64 t #s(literal 3 binary64)) (*.f64 l l)) (sin.f64 k)) (tan.f64 k)) (+.f64 (+.f64 #s(literal 1 binary64) (pow.f64 (/.f64 k t) #s(literal 2 binary64))) #s(literal 1 binary64)))) < 1e188Initial program 80.6%
Taylor expanded in k around 0
unpow2N/A
*-commutativeN/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
lower-pow.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f6470.4
Applied rewrites70.4%
Applied rewrites70.4%
Applied rewrites70.4%
Applied rewrites70.3%
if 1e188 < (/.f64 #s(literal 2 binary64) (*.f64 (*.f64 (*.f64 (/.f64 (pow.f64 t #s(literal 3 binary64)) (*.f64 l l)) (sin.f64 k)) (tan.f64 k)) (+.f64 (+.f64 #s(literal 1 binary64) (pow.f64 (/.f64 k t) #s(literal 2 binary64))) #s(literal 1 binary64)))) Initial program 27.6%
Applied rewrites34.1%
Taylor expanded in k around 0
associate-*r/N/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
lower-pow.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f6450.2
Applied rewrites50.2%
Applied rewrites49.2%
Applied rewrites50.2%
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k)
:precision binary64
(let* ((t_2 (* (sin k) t_m)))
(*
t_s
(if (<= t_m 4.9e-60)
(/ 2.0 (/ (/ (* (* (* (pow (sin k) 2.0) t_m) k) k) (* (cos k) l)) l))
(if (<= t_m 3.4e+158)
(*
(/
(* (/ 2.0 (* t_2 t_m)) l)
(* (* t_m (+ (pow (/ k t_m) 2.0) 2.0)) (tan k)))
l)
(/
2.0
(*
(* (* t_2 (* (/ t_m l) (/ t_m l))) (tan k))
(fma (/ k t_m) (/ k t_m) 2.0))))))))t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
double t_2 = sin(k) * t_m;
double tmp;
if (t_m <= 4.9e-60) {
tmp = 2.0 / (((((pow(sin(k), 2.0) * t_m) * k) * k) / (cos(k) * l)) / l);
} else if (t_m <= 3.4e+158) {
tmp = (((2.0 / (t_2 * t_m)) * l) / ((t_m * (pow((k / t_m), 2.0) + 2.0)) * tan(k))) * l;
} else {
tmp = 2.0 / (((t_2 * ((t_m / l) * (t_m / l))) * tan(k)) * fma((k / t_m), (k / t_m), 2.0));
}
return t_s * tmp;
}
t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k) t_2 = Float64(sin(k) * t_m) tmp = 0.0 if (t_m <= 4.9e-60) tmp = Float64(2.0 / Float64(Float64(Float64(Float64(Float64((sin(k) ^ 2.0) * t_m) * k) * k) / Float64(cos(k) * l)) / l)); elseif (t_m <= 3.4e+158) tmp = Float64(Float64(Float64(Float64(2.0 / Float64(t_2 * t_m)) * l) / Float64(Float64(t_m * Float64((Float64(k / t_m) ^ 2.0) + 2.0)) * tan(k))) * l); else tmp = Float64(2.0 / Float64(Float64(Float64(t_2 * Float64(Float64(t_m / l) * Float64(t_m / l))) * tan(k)) * fma(Float64(k / t_m), Float64(k / t_m), 2.0))); end return Float64(t_s * tmp) end
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k_] := Block[{t$95$2 = N[(N[Sin[k], $MachinePrecision] * t$95$m), $MachinePrecision]}, N[(t$95$s * If[LessEqual[t$95$m, 4.9e-60], N[(2.0 / N[(N[(N[(N[(N[(N[Power[N[Sin[k], $MachinePrecision], 2.0], $MachinePrecision] * t$95$m), $MachinePrecision] * k), $MachinePrecision] * k), $MachinePrecision] / N[(N[Cos[k], $MachinePrecision] * l), $MachinePrecision]), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$m, 3.4e+158], N[(N[(N[(N[(2.0 / N[(t$95$2 * t$95$m), $MachinePrecision]), $MachinePrecision] * l), $MachinePrecision] / N[(N[(t$95$m * N[(N[Power[N[(k / t$95$m), $MachinePrecision], 2.0], $MachinePrecision] + 2.0), $MachinePrecision]), $MachinePrecision] * N[Tan[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * l), $MachinePrecision], N[(2.0 / N[(N[(N[(t$95$2 * N[(N[(t$95$m / l), $MachinePrecision] * N[(t$95$m / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[Tan[k], $MachinePrecision]), $MachinePrecision] * N[(N[(k / t$95$m), $MachinePrecision] * N[(k / t$95$m), $MachinePrecision] + 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]]
\begin{array}{l}
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
\begin{array}{l}
t_2 := \sin k \cdot t\_m\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_m \leq 4.9 \cdot 10^{-60}:\\
\;\;\;\;\frac{2}{\frac{\frac{\left(\left({\sin k}^{2} \cdot t\_m\right) \cdot k\right) \cdot k}{\cos k \cdot \ell}}{\ell}}\\
\mathbf{elif}\;t\_m \leq 3.4 \cdot 10^{+158}:\\
\;\;\;\;\frac{\frac{2}{t\_2 \cdot t\_m} \cdot \ell}{\left(t\_m \cdot \left({\left(\frac{k}{t\_m}\right)}^{2} + 2\right)\right) \cdot \tan k} \cdot \ell\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\left(\left(t\_2 \cdot \left(\frac{t\_m}{\ell} \cdot \frac{t\_m}{\ell}\right)\right) \cdot \tan k\right) \cdot \mathsf{fma}\left(\frac{k}{t\_m}, \frac{k}{t\_m}, 2\right)}\\
\end{array}
\end{array}
\end{array}
if t < 4.89999999999999988e-60Initial program 54.2%
Applied rewrites48.4%
Taylor expanded in t around 0
lower-/.f64N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-pow.f64N/A
lower-sin.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-cos.f6477.6
Applied rewrites77.6%
if 4.89999999999999988e-60 < t < 3.3999999999999999e158Initial program 71.2%
lift-*.f64N/A
*-commutativeN/A
lift-/.f64N/A
lift-pow.f64N/A
cube-multN/A
associate-/l*N/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
lift-*.f64N/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
lower-/.f6484.9
Applied rewrites84.9%
lift-*.f64N/A
lift-*.f64N/A
associate-*l*N/A
lift-*.f64N/A
lift-*.f64N/A
associate-*r*N/A
associate-*l*N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lift-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f64N/A
Applied rewrites91.5%
Applied rewrites95.3%
if 3.3999999999999999e158 < t Initial program 74.5%
lift-*.f64N/A
*-commutativeN/A
lift-/.f64N/A
lift-pow.f64N/A
cube-multN/A
associate-/l*N/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
lift-*.f64N/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
lower-/.f6496.8
Applied rewrites96.8%
lift-+.f64N/A
lift-+.f64N/A
+-commutativeN/A
associate-+l+N/A
metadata-evalN/A
lift-pow.f64N/A
unpow2N/A
lower-fma.f6496.8
Applied rewrites96.8%
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k)
:precision binary64
(*
t_s
(if (<= t_m 1e-58)
(/ 2.0 (/ (/ (* (* (* (pow (sin k) 2.0) t_m) k) k) (* (cos k) l)) l))
(/
2.0
(*
(* (* (* (sin k) t_m) (* (/ t_m l) (/ t_m l))) (tan k))
(fma (/ k t_m) (/ k t_m) 2.0))))))t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
double tmp;
if (t_m <= 1e-58) {
tmp = 2.0 / (((((pow(sin(k), 2.0) * t_m) * k) * k) / (cos(k) * l)) / l);
} else {
tmp = 2.0 / ((((sin(k) * t_m) * ((t_m / l) * (t_m / l))) * tan(k)) * fma((k / t_m), (k / t_m), 2.0));
}
return t_s * tmp;
}
t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k) tmp = 0.0 if (t_m <= 1e-58) tmp = Float64(2.0 / Float64(Float64(Float64(Float64(Float64((sin(k) ^ 2.0) * t_m) * k) * k) / Float64(cos(k) * l)) / l)); else tmp = Float64(2.0 / Float64(Float64(Float64(Float64(sin(k) * t_m) * Float64(Float64(t_m / l) * Float64(t_m / l))) * tan(k)) * fma(Float64(k / t_m), Float64(k / t_m), 2.0))); end return Float64(t_s * tmp) end
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k_] := N[(t$95$s * If[LessEqual[t$95$m, 1e-58], N[(2.0 / N[(N[(N[(N[(N[(N[Power[N[Sin[k], $MachinePrecision], 2.0], $MachinePrecision] * t$95$m), $MachinePrecision] * k), $MachinePrecision] * k), $MachinePrecision] / N[(N[Cos[k], $MachinePrecision] * l), $MachinePrecision]), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(N[(N[(N[Sin[k], $MachinePrecision] * t$95$m), $MachinePrecision] * N[(N[(t$95$m / l), $MachinePrecision] * N[(t$95$m / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[Tan[k], $MachinePrecision]), $MachinePrecision] * N[(N[(k / t$95$m), $MachinePrecision] * N[(k / t$95$m), $MachinePrecision] + 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_m \leq 10^{-58}:\\
\;\;\;\;\frac{2}{\frac{\frac{\left(\left({\sin k}^{2} \cdot t\_m\right) \cdot k\right) \cdot k}{\cos k \cdot \ell}}{\ell}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\left(\left(\left(\sin k \cdot t\_m\right) \cdot \left(\frac{t\_m}{\ell} \cdot \frac{t\_m}{\ell}\right)\right) \cdot \tan k\right) \cdot \mathsf{fma}\left(\frac{k}{t\_m}, \frac{k}{t\_m}, 2\right)}\\
\end{array}
\end{array}
if t < 1e-58Initial program 54.2%
Applied rewrites48.4%
Taylor expanded in t around 0
lower-/.f64N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-pow.f64N/A
lower-sin.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-cos.f6477.6
Applied rewrites77.6%
if 1e-58 < t Initial program 72.6%
lift-*.f64N/A
*-commutativeN/A
lift-/.f64N/A
lift-pow.f64N/A
cube-multN/A
associate-/l*N/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
lift-*.f64N/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
lower-/.f6489.9
Applied rewrites89.9%
lift-+.f64N/A
lift-+.f64N/A
+-commutativeN/A
associate-+l+N/A
metadata-evalN/A
lift-pow.f64N/A
unpow2N/A
lower-fma.f6489.9
Applied rewrites89.9%
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k)
:precision binary64
(*
t_s
(if (<= t_m 3.8e-45)
(/ 2.0 (/ (* (* (/ (sin k) l) (tan k)) (* (* k k) t_m)) l))
(/
2.0
(*
(* (* (* (sin k) t_m) (* (/ t_m l) (/ t_m l))) (tan k))
(fma (/ k t_m) (/ k t_m) 2.0))))))t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
double tmp;
if (t_m <= 3.8e-45) {
tmp = 2.0 / ((((sin(k) / l) * tan(k)) * ((k * k) * t_m)) / l);
} else {
tmp = 2.0 / ((((sin(k) * t_m) * ((t_m / l) * (t_m / l))) * tan(k)) * fma((k / t_m), (k / t_m), 2.0));
}
return t_s * tmp;
}
t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k) tmp = 0.0 if (t_m <= 3.8e-45) tmp = Float64(2.0 / Float64(Float64(Float64(Float64(sin(k) / l) * tan(k)) * Float64(Float64(k * k) * t_m)) / l)); else tmp = Float64(2.0 / Float64(Float64(Float64(Float64(sin(k) * t_m) * Float64(Float64(t_m / l) * Float64(t_m / l))) * tan(k)) * fma(Float64(k / t_m), Float64(k / t_m), 2.0))); end return Float64(t_s * tmp) end
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k_] := N[(t$95$s * If[LessEqual[t$95$m, 3.8e-45], N[(2.0 / N[(N[(N[(N[(N[Sin[k], $MachinePrecision] / l), $MachinePrecision] * N[Tan[k], $MachinePrecision]), $MachinePrecision] * N[(N[(k * k), $MachinePrecision] * t$95$m), $MachinePrecision]), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(N[(N[(N[Sin[k], $MachinePrecision] * t$95$m), $MachinePrecision] * N[(N[(t$95$m / l), $MachinePrecision] * N[(t$95$m / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[Tan[k], $MachinePrecision]), $MachinePrecision] * N[(N[(k / t$95$m), $MachinePrecision] * N[(k / t$95$m), $MachinePrecision] + 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_m \leq 3.8 \cdot 10^{-45}:\\
\;\;\;\;\frac{2}{\frac{\left(\frac{\sin k}{\ell} \cdot \tan k\right) \cdot \left(\left(k \cdot k\right) \cdot t\_m\right)}{\ell}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\left(\left(\left(\sin k \cdot t\_m\right) \cdot \left(\frac{t\_m}{\ell} \cdot \frac{t\_m}{\ell}\right)\right) \cdot \tan k\right) \cdot \mathsf{fma}\left(\frac{k}{t\_m}, \frac{k}{t\_m}, 2\right)}\\
\end{array}
\end{array}
if t < 3.79999999999999997e-45Initial program 54.8%
Taylor expanded in t around 0
associate-*r*N/A
associate-/l*N/A
lower-*.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
lower-/.f64N/A
lower-pow.f64N/A
lower-sin.f64N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
lower-cos.f6469.6
Applied rewrites69.6%
Applied rewrites74.6%
if 3.79999999999999997e-45 < t Initial program 72.0%
lift-*.f64N/A
*-commutativeN/A
lift-/.f64N/A
lift-pow.f64N/A
cube-multN/A
associate-/l*N/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
lift-*.f64N/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
lower-/.f6490.3
Applied rewrites90.3%
lift-+.f64N/A
lift-+.f64N/A
+-commutativeN/A
associate-+l+N/A
metadata-evalN/A
lift-pow.f64N/A
unpow2N/A
lower-fma.f6490.3
Applied rewrites90.3%
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k)
:precision binary64
(*
t_s
(if (<= t_m 3.8e-45)
(/ 2.0 (/ (* (* (/ (sin k) l) (tan k)) (* (* k k) t_m)) l))
(/
2.0
(*
(* (* (* (sin k) t_m) (* (/ t_m l) (/ t_m l))) (tan k))
(fma k (/ k (* t_m t_m)) 2.0))))))t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
double tmp;
if (t_m <= 3.8e-45) {
tmp = 2.0 / ((((sin(k) / l) * tan(k)) * ((k * k) * t_m)) / l);
} else {
tmp = 2.0 / ((((sin(k) * t_m) * ((t_m / l) * (t_m / l))) * tan(k)) * fma(k, (k / (t_m * t_m)), 2.0));
}
return t_s * tmp;
}
t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k) tmp = 0.0 if (t_m <= 3.8e-45) tmp = Float64(2.0 / Float64(Float64(Float64(Float64(sin(k) / l) * tan(k)) * Float64(Float64(k * k) * t_m)) / l)); else tmp = Float64(2.0 / Float64(Float64(Float64(Float64(sin(k) * t_m) * Float64(Float64(t_m / l) * Float64(t_m / l))) * tan(k)) * fma(k, Float64(k / Float64(t_m * t_m)), 2.0))); end return Float64(t_s * tmp) end
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k_] := N[(t$95$s * If[LessEqual[t$95$m, 3.8e-45], N[(2.0 / N[(N[(N[(N[(N[Sin[k], $MachinePrecision] / l), $MachinePrecision] * N[Tan[k], $MachinePrecision]), $MachinePrecision] * N[(N[(k * k), $MachinePrecision] * t$95$m), $MachinePrecision]), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(N[(N[(N[Sin[k], $MachinePrecision] * t$95$m), $MachinePrecision] * N[(N[(t$95$m / l), $MachinePrecision] * N[(t$95$m / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[Tan[k], $MachinePrecision]), $MachinePrecision] * N[(k * N[(k / N[(t$95$m * t$95$m), $MachinePrecision]), $MachinePrecision] + 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_m \leq 3.8 \cdot 10^{-45}:\\
\;\;\;\;\frac{2}{\frac{\left(\frac{\sin k}{\ell} \cdot \tan k\right) \cdot \left(\left(k \cdot k\right) \cdot t\_m\right)}{\ell}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\left(\left(\left(\sin k \cdot t\_m\right) \cdot \left(\frac{t\_m}{\ell} \cdot \frac{t\_m}{\ell}\right)\right) \cdot \tan k\right) \cdot \mathsf{fma}\left(k, \frac{k}{t\_m \cdot t\_m}, 2\right)}\\
\end{array}
\end{array}
if t < 3.79999999999999997e-45Initial program 54.8%
Taylor expanded in t around 0
associate-*r*N/A
associate-/l*N/A
lower-*.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
lower-/.f64N/A
lower-pow.f64N/A
lower-sin.f64N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
lower-cos.f6469.6
Applied rewrites69.6%
Applied rewrites74.6%
if 3.79999999999999997e-45 < t Initial program 72.0%
lift-*.f64N/A
*-commutativeN/A
lift-/.f64N/A
lift-pow.f64N/A
cube-multN/A
associate-/l*N/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
lift-*.f64N/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
lower-/.f6490.3
Applied rewrites90.3%
Taylor expanded in t around 0
lower-/.f64N/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6446.0
Applied rewrites46.0%
Applied rewrites87.6%
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k)
:precision binary64
(*
t_s
(if (<= t_m 5.7e-86)
(/
2.0
(*
(/ t_m l)
(*
(/
(fma
(* (fma 0.3333333333333333 (* t_m t_m) 1.0) k)
k
(* (* t_m t_m) 2.0))
l)
(* k k))))
(/
2.0
(*
(* (/ t_m l) (* t_m (sin k)))
(*
(/ t_m l)
(*
(fma (+ 0.6666666666666666 (pow (* t_m t_m) -1.0)) (* k k) 2.0)
k)))))))t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
double tmp;
if (t_m <= 5.7e-86) {
tmp = 2.0 / ((t_m / l) * ((fma((fma(0.3333333333333333, (t_m * t_m), 1.0) * k), k, ((t_m * t_m) * 2.0)) / l) * (k * k)));
} else {
tmp = 2.0 / (((t_m / l) * (t_m * sin(k))) * ((t_m / l) * (fma((0.6666666666666666 + pow((t_m * t_m), -1.0)), (k * k), 2.0) * k)));
}
return t_s * tmp;
}
t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k) tmp = 0.0 if (t_m <= 5.7e-86) tmp = Float64(2.0 / Float64(Float64(t_m / l) * Float64(Float64(fma(Float64(fma(0.3333333333333333, Float64(t_m * t_m), 1.0) * k), k, Float64(Float64(t_m * t_m) * 2.0)) / l) * Float64(k * k)))); else tmp = Float64(2.0 / Float64(Float64(Float64(t_m / l) * Float64(t_m * sin(k))) * Float64(Float64(t_m / l) * Float64(fma(Float64(0.6666666666666666 + (Float64(t_m * t_m) ^ -1.0)), Float64(k * k), 2.0) * k)))); end return Float64(t_s * tmp) end
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k_] := N[(t$95$s * If[LessEqual[t$95$m, 5.7e-86], N[(2.0 / N[(N[(t$95$m / l), $MachinePrecision] * N[(N[(N[(N[(N[(0.3333333333333333 * N[(t$95$m * t$95$m), $MachinePrecision] + 1.0), $MachinePrecision] * k), $MachinePrecision] * k + N[(N[(t$95$m * t$95$m), $MachinePrecision] * 2.0), $MachinePrecision]), $MachinePrecision] / l), $MachinePrecision] * N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(N[(t$95$m / l), $MachinePrecision] * N[(t$95$m * N[Sin[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[(t$95$m / l), $MachinePrecision] * N[(N[(N[(0.6666666666666666 + N[Power[N[(t$95$m * t$95$m), $MachinePrecision], -1.0], $MachinePrecision]), $MachinePrecision] * N[(k * k), $MachinePrecision] + 2.0), $MachinePrecision] * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_m \leq 5.7 \cdot 10^{-86}:\\
\;\;\;\;\frac{2}{\frac{t\_m}{\ell} \cdot \left(\frac{\mathsf{fma}\left(\mathsf{fma}\left(0.3333333333333333, t\_m \cdot t\_m, 1\right) \cdot k, k, \left(t\_m \cdot t\_m\right) \cdot 2\right)}{\ell} \cdot \left(k \cdot k\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\left(\frac{t\_m}{\ell} \cdot \left(t\_m \cdot \sin k\right)\right) \cdot \left(\frac{t\_m}{\ell} \cdot \left(\mathsf{fma}\left(0.6666666666666666 + {\left(t\_m \cdot t\_m\right)}^{-1}, k \cdot k, 2\right) \cdot k\right)\right)}\\
\end{array}
\end{array}
if t < 5.7000000000000004e-86Initial program 54.0%
lift-*.f64N/A
lift-*.f64N/A
associate-*l*N/A
lift-*.f64N/A
associate-*l*N/A
lift-/.f64N/A
lift-pow.f64N/A
cube-multN/A
lift-*.f64N/A
times-fracN/A
associate-*l*N/A
lower-*.f64N/A
lower-/.f64N/A
lower-*.f64N/A
Applied rewrites67.4%
Taylor expanded in k around 0
*-commutativeN/A
lower-*.f64N/A
Applied rewrites69.6%
if 5.7000000000000004e-86 < t Initial program 71.3%
lift-*.f64N/A
*-commutativeN/A
lift-/.f64N/A
lift-pow.f64N/A
cube-multN/A
associate-/l*N/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
lift-*.f64N/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
lower-/.f6487.2
Applied rewrites87.2%
lift-*.f64N/A
lift-*.f64N/A
associate-*l*N/A
lift-*.f64N/A
lift-*.f64N/A
associate-*r*N/A
associate-*l*N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lift-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f64N/A
Applied rewrites94.2%
Taylor expanded in k 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-*.f6486.2
Applied rewrites86.2%
Final simplification74.9%
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k)
:precision binary64
(*
t_s
(if (<= k 6.5e-161)
(/ 2.0 (* (* (* (/ t_m l) (sin k)) t_m) (* (/ (* k t_m) l) 2.0)))
(if (<= k 0.48)
(/
2.0
(*
(/ t_m l)
(*
(/
(fma
(* (fma 0.3333333333333333 (* t_m t_m) 1.0) k)
k
(* (* t_m t_m) 2.0))
l)
(* k k))))
(if (<= k 2.05e+148)
(/ 2.0 (/ (* (* (/ (sin k) l) (tan k)) (* (* k k) t_m)) l))
(/ 2.0 (* k (* (* t_m k) (* (tan k) (/ (sin k) (* l l)))))))))))t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
double tmp;
if (k <= 6.5e-161) {
tmp = 2.0 / ((((t_m / l) * sin(k)) * t_m) * (((k * t_m) / l) * 2.0));
} else if (k <= 0.48) {
tmp = 2.0 / ((t_m / l) * ((fma((fma(0.3333333333333333, (t_m * t_m), 1.0) * k), k, ((t_m * t_m) * 2.0)) / l) * (k * k)));
} else if (k <= 2.05e+148) {
tmp = 2.0 / ((((sin(k) / l) * tan(k)) * ((k * k) * t_m)) / l);
} else {
tmp = 2.0 / (k * ((t_m * k) * (tan(k) * (sin(k) / (l * l)))));
}
return t_s * tmp;
}
t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k) tmp = 0.0 if (k <= 6.5e-161) tmp = Float64(2.0 / Float64(Float64(Float64(Float64(t_m / l) * sin(k)) * t_m) * Float64(Float64(Float64(k * t_m) / l) * 2.0))); elseif (k <= 0.48) tmp = Float64(2.0 / Float64(Float64(t_m / l) * Float64(Float64(fma(Float64(fma(0.3333333333333333, Float64(t_m * t_m), 1.0) * k), k, Float64(Float64(t_m * t_m) * 2.0)) / l) * Float64(k * k)))); elseif (k <= 2.05e+148) tmp = Float64(2.0 / Float64(Float64(Float64(Float64(sin(k) / l) * tan(k)) * Float64(Float64(k * k) * t_m)) / l)); else tmp = Float64(2.0 / Float64(k * Float64(Float64(t_m * k) * Float64(tan(k) * Float64(sin(k) / Float64(l * l)))))); end return Float64(t_s * tmp) end
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k_] := N[(t$95$s * If[LessEqual[k, 6.5e-161], N[(2.0 / N[(N[(N[(N[(t$95$m / l), $MachinePrecision] * N[Sin[k], $MachinePrecision]), $MachinePrecision] * t$95$m), $MachinePrecision] * N[(N[(N[(k * t$95$m), $MachinePrecision] / l), $MachinePrecision] * 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[k, 0.48], N[(2.0 / N[(N[(t$95$m / l), $MachinePrecision] * N[(N[(N[(N[(N[(0.3333333333333333 * N[(t$95$m * t$95$m), $MachinePrecision] + 1.0), $MachinePrecision] * k), $MachinePrecision] * k + N[(N[(t$95$m * t$95$m), $MachinePrecision] * 2.0), $MachinePrecision]), $MachinePrecision] / l), $MachinePrecision] * N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[k, 2.05e+148], N[(2.0 / N[(N[(N[(N[(N[Sin[k], $MachinePrecision] / l), $MachinePrecision] * N[Tan[k], $MachinePrecision]), $MachinePrecision] * N[(N[(k * k), $MachinePrecision] * t$95$m), $MachinePrecision]), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(k * N[(N[(t$95$m * k), $MachinePrecision] * N[(N[Tan[k], $MachinePrecision] * N[(N[Sin[k], $MachinePrecision] / N[(l * l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]), $MachinePrecision]
\begin{array}{l}
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k \leq 6.5 \cdot 10^{-161}:\\
\;\;\;\;\frac{2}{\left(\left(\frac{t\_m}{\ell} \cdot \sin k\right) \cdot t\_m\right) \cdot \left(\frac{k \cdot t\_m}{\ell} \cdot 2\right)}\\
\mathbf{elif}\;k \leq 0.48:\\
\;\;\;\;\frac{2}{\frac{t\_m}{\ell} \cdot \left(\frac{\mathsf{fma}\left(\mathsf{fma}\left(0.3333333333333333, t\_m \cdot t\_m, 1\right) \cdot k, k, \left(t\_m \cdot t\_m\right) \cdot 2\right)}{\ell} \cdot \left(k \cdot k\right)\right)}\\
\mathbf{elif}\;k \leq 2.05 \cdot 10^{+148}:\\
\;\;\;\;\frac{2}{\frac{\left(\frac{\sin k}{\ell} \cdot \tan k\right) \cdot \left(\left(k \cdot k\right) \cdot t\_m\right)}{\ell}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{k \cdot \left(\left(t\_m \cdot k\right) \cdot \left(\tan k \cdot \frac{\sin k}{\ell \cdot \ell}\right)\right)}\\
\end{array}
\end{array}
if k < 6.50000000000000008e-161Initial program 59.6%
lift-*.f64N/A
*-commutativeN/A
lift-/.f64N/A
lift-pow.f64N/A
cube-multN/A
associate-/l*N/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
lift-*.f64N/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
lower-/.f6475.1
Applied rewrites75.1%
lift-*.f64N/A
lift-*.f64N/A
associate-*l*N/A
lift-*.f64N/A
lift-*.f64N/A
associate-*r*N/A
associate-*l*N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lift-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f64N/A
Applied rewrites83.7%
lift-*.f64N/A
lift-*.f64N/A
*-commutativeN/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f6484.5
Applied rewrites84.5%
Taylor expanded in k around 0
*-commutativeN/A
lower-*.f64N/A
lower-/.f64N/A
lower-*.f6472.4
Applied rewrites72.4%
if 6.50000000000000008e-161 < k < 0.47999999999999998Initial program 67.9%
lift-*.f64N/A
lift-*.f64N/A
associate-*l*N/A
lift-*.f64N/A
associate-*l*N/A
lift-/.f64N/A
lift-pow.f64N/A
cube-multN/A
lift-*.f64N/A
times-fracN/A
associate-*l*N/A
lower-*.f64N/A
lower-/.f64N/A
lower-*.f64N/A
Applied rewrites85.2%
Taylor expanded in k around 0
*-commutativeN/A
lower-*.f64N/A
Applied rewrites98.0%
if 0.47999999999999998 < k < 2.0499999999999999e148Initial program 50.3%
Taylor expanded in t around 0
associate-*r*N/A
associate-/l*N/A
lower-*.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
lower-/.f64N/A
lower-pow.f64N/A
lower-sin.f64N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
lower-cos.f6475.7
Applied rewrites75.7%
Applied rewrites84.0%
if 2.0499999999999999e148 < k Initial program 61.8%
Taylor expanded in t around 0
associate-*r*N/A
associate-/l*N/A
lower-*.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
lower-/.f64N/A
lower-pow.f64N/A
lower-sin.f64N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
lower-cos.f6471.6
Applied rewrites71.6%
Applied rewrites84.3%
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k)
:precision binary64
(*
t_s
(if (<= k 6.5e-161)
(/ 2.0 (* (* (* (/ t_m l) (sin k)) t_m) (* (/ (* k t_m) l) 2.0)))
(if (<= k 900000.0)
(/
2.0
(*
(/ t_m l)
(*
(/
(fma
(* (fma 0.3333333333333333 (* t_m t_m) 1.0) k)
k
(* (* t_m t_m) 2.0))
l)
(* k k))))
(/ 2.0 (* k (* (* t_m k) (* (tan k) (/ (sin k) (* l l))))))))))t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
double tmp;
if (k <= 6.5e-161) {
tmp = 2.0 / ((((t_m / l) * sin(k)) * t_m) * (((k * t_m) / l) * 2.0));
} else if (k <= 900000.0) {
tmp = 2.0 / ((t_m / l) * ((fma((fma(0.3333333333333333, (t_m * t_m), 1.0) * k), k, ((t_m * t_m) * 2.0)) / l) * (k * k)));
} else {
tmp = 2.0 / (k * ((t_m * k) * (tan(k) * (sin(k) / (l * l)))));
}
return t_s * tmp;
}
t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k) tmp = 0.0 if (k <= 6.5e-161) tmp = Float64(2.0 / Float64(Float64(Float64(Float64(t_m / l) * sin(k)) * t_m) * Float64(Float64(Float64(k * t_m) / l) * 2.0))); elseif (k <= 900000.0) tmp = Float64(2.0 / Float64(Float64(t_m / l) * Float64(Float64(fma(Float64(fma(0.3333333333333333, Float64(t_m * t_m), 1.0) * k), k, Float64(Float64(t_m * t_m) * 2.0)) / l) * Float64(k * k)))); else tmp = Float64(2.0 / Float64(k * Float64(Float64(t_m * k) * Float64(tan(k) * Float64(sin(k) / Float64(l * l)))))); end return Float64(t_s * tmp) end
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k_] := N[(t$95$s * If[LessEqual[k, 6.5e-161], N[(2.0 / N[(N[(N[(N[(t$95$m / l), $MachinePrecision] * N[Sin[k], $MachinePrecision]), $MachinePrecision] * t$95$m), $MachinePrecision] * N[(N[(N[(k * t$95$m), $MachinePrecision] / l), $MachinePrecision] * 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[k, 900000.0], N[(2.0 / N[(N[(t$95$m / l), $MachinePrecision] * N[(N[(N[(N[(N[(0.3333333333333333 * N[(t$95$m * t$95$m), $MachinePrecision] + 1.0), $MachinePrecision] * k), $MachinePrecision] * k + N[(N[(t$95$m * t$95$m), $MachinePrecision] * 2.0), $MachinePrecision]), $MachinePrecision] / l), $MachinePrecision] * N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(k * N[(N[(t$95$m * k), $MachinePrecision] * N[(N[Tan[k], $MachinePrecision] * N[(N[Sin[k], $MachinePrecision] / N[(l * l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]
\begin{array}{l}
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k \leq 6.5 \cdot 10^{-161}:\\
\;\;\;\;\frac{2}{\left(\left(\frac{t\_m}{\ell} \cdot \sin k\right) \cdot t\_m\right) \cdot \left(\frac{k \cdot t\_m}{\ell} \cdot 2\right)}\\
\mathbf{elif}\;k \leq 900000:\\
\;\;\;\;\frac{2}{\frac{t\_m}{\ell} \cdot \left(\frac{\mathsf{fma}\left(\mathsf{fma}\left(0.3333333333333333, t\_m \cdot t\_m, 1\right) \cdot k, k, \left(t\_m \cdot t\_m\right) \cdot 2\right)}{\ell} \cdot \left(k \cdot k\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{k \cdot \left(\left(t\_m \cdot k\right) \cdot \left(\tan k \cdot \frac{\sin k}{\ell \cdot \ell}\right)\right)}\\
\end{array}
\end{array}
if k < 6.50000000000000008e-161Initial program 59.6%
lift-*.f64N/A
*-commutativeN/A
lift-/.f64N/A
lift-pow.f64N/A
cube-multN/A
associate-/l*N/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
lift-*.f64N/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
lower-/.f6475.1
Applied rewrites75.1%
lift-*.f64N/A
lift-*.f64N/A
associate-*l*N/A
lift-*.f64N/A
lift-*.f64N/A
associate-*r*N/A
associate-*l*N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lift-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f64N/A
Applied rewrites83.7%
lift-*.f64N/A
lift-*.f64N/A
*-commutativeN/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f6484.5
Applied rewrites84.5%
Taylor expanded in k around 0
*-commutativeN/A
lower-*.f64N/A
lower-/.f64N/A
lower-*.f6472.4
Applied rewrites72.4%
if 6.50000000000000008e-161 < k < 9e5Initial program 65.6%
lift-*.f64N/A
lift-*.f64N/A
associate-*l*N/A
lift-*.f64N/A
associate-*l*N/A
lift-/.f64N/A
lift-pow.f64N/A
cube-multN/A
lift-*.f64N/A
times-fracN/A
associate-*l*N/A
lower-*.f64N/A
lower-/.f64N/A
lower-*.f64N/A
Applied rewrites85.7%
Taylor expanded in k around 0
*-commutativeN/A
lower-*.f64N/A
Applied rewrites95.1%
if 9e5 < k Initial program 56.5%
Taylor expanded in t around 0
associate-*r*N/A
associate-/l*N/A
lower-*.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
lower-/.f64N/A
lower-pow.f64N/A
lower-sin.f64N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
lower-cos.f6474.9
Applied rewrites74.9%
Applied rewrites80.9%
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k)
:precision binary64
(*
t_s
(if (<= k 6.5e-161)
(/ 2.0 (* (* (* (/ t_m l) (sin k)) t_m) (* (/ (* k t_m) l) 2.0)))
(if (<= k 900000.0)
(/
2.0
(*
(/ t_m l)
(*
(/
(fma
(* (fma 0.3333333333333333 (* t_m t_m) 1.0) k)
k
(* (* t_m t_m) 2.0))
l)
(* k k))))
(/ 2.0 (* k (* k (* (* (tan k) (/ (sin k) (* l l))) t_m))))))))t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
double tmp;
if (k <= 6.5e-161) {
tmp = 2.0 / ((((t_m / l) * sin(k)) * t_m) * (((k * t_m) / l) * 2.0));
} else if (k <= 900000.0) {
tmp = 2.0 / ((t_m / l) * ((fma((fma(0.3333333333333333, (t_m * t_m), 1.0) * k), k, ((t_m * t_m) * 2.0)) / l) * (k * k)));
} else {
tmp = 2.0 / (k * (k * ((tan(k) * (sin(k) / (l * l))) * t_m)));
}
return t_s * tmp;
}
t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k) tmp = 0.0 if (k <= 6.5e-161) tmp = Float64(2.0 / Float64(Float64(Float64(Float64(t_m / l) * sin(k)) * t_m) * Float64(Float64(Float64(k * t_m) / l) * 2.0))); elseif (k <= 900000.0) tmp = Float64(2.0 / Float64(Float64(t_m / l) * Float64(Float64(fma(Float64(fma(0.3333333333333333, Float64(t_m * t_m), 1.0) * k), k, Float64(Float64(t_m * t_m) * 2.0)) / l) * Float64(k * k)))); else tmp = Float64(2.0 / Float64(k * Float64(k * Float64(Float64(tan(k) * Float64(sin(k) / Float64(l * l))) * t_m)))); end return Float64(t_s * tmp) end
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k_] := N[(t$95$s * If[LessEqual[k, 6.5e-161], N[(2.0 / N[(N[(N[(N[(t$95$m / l), $MachinePrecision] * N[Sin[k], $MachinePrecision]), $MachinePrecision] * t$95$m), $MachinePrecision] * N[(N[(N[(k * t$95$m), $MachinePrecision] / l), $MachinePrecision] * 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[k, 900000.0], N[(2.0 / N[(N[(t$95$m / l), $MachinePrecision] * N[(N[(N[(N[(N[(0.3333333333333333 * N[(t$95$m * t$95$m), $MachinePrecision] + 1.0), $MachinePrecision] * k), $MachinePrecision] * k + N[(N[(t$95$m * t$95$m), $MachinePrecision] * 2.0), $MachinePrecision]), $MachinePrecision] / l), $MachinePrecision] * N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(k * N[(k * N[(N[(N[Tan[k], $MachinePrecision] * N[(N[Sin[k], $MachinePrecision] / N[(l * l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * t$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]
\begin{array}{l}
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k \leq 6.5 \cdot 10^{-161}:\\
\;\;\;\;\frac{2}{\left(\left(\frac{t\_m}{\ell} \cdot \sin k\right) \cdot t\_m\right) \cdot \left(\frac{k \cdot t\_m}{\ell} \cdot 2\right)}\\
\mathbf{elif}\;k \leq 900000:\\
\;\;\;\;\frac{2}{\frac{t\_m}{\ell} \cdot \left(\frac{\mathsf{fma}\left(\mathsf{fma}\left(0.3333333333333333, t\_m \cdot t\_m, 1\right) \cdot k, k, \left(t\_m \cdot t\_m\right) \cdot 2\right)}{\ell} \cdot \left(k \cdot k\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{k \cdot \left(k \cdot \left(\left(\tan k \cdot \frac{\sin k}{\ell \cdot \ell}\right) \cdot t\_m\right)\right)}\\
\end{array}
\end{array}
if k < 6.50000000000000008e-161Initial program 59.6%
lift-*.f64N/A
*-commutativeN/A
lift-/.f64N/A
lift-pow.f64N/A
cube-multN/A
associate-/l*N/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
lift-*.f64N/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
lower-/.f6475.1
Applied rewrites75.1%
lift-*.f64N/A
lift-*.f64N/A
associate-*l*N/A
lift-*.f64N/A
lift-*.f64N/A
associate-*r*N/A
associate-*l*N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lift-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f64N/A
Applied rewrites83.7%
lift-*.f64N/A
lift-*.f64N/A
*-commutativeN/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f6484.5
Applied rewrites84.5%
Taylor expanded in k around 0
*-commutativeN/A
lower-*.f64N/A
lower-/.f64N/A
lower-*.f6472.4
Applied rewrites72.4%
if 6.50000000000000008e-161 < k < 9e5Initial program 65.6%
lift-*.f64N/A
lift-*.f64N/A
associate-*l*N/A
lift-*.f64N/A
associate-*l*N/A
lift-/.f64N/A
lift-pow.f64N/A
cube-multN/A
lift-*.f64N/A
times-fracN/A
associate-*l*N/A
lower-*.f64N/A
lower-/.f64N/A
lower-*.f64N/A
Applied rewrites85.7%
Taylor expanded in k around 0
*-commutativeN/A
lower-*.f64N/A
Applied rewrites95.1%
if 9e5 < k Initial program 56.5%
Taylor expanded in t around 0
associate-*r*N/A
associate-/l*N/A
lower-*.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
lower-/.f64N/A
lower-pow.f64N/A
lower-sin.f64N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
lower-cos.f6474.9
Applied rewrites74.9%
Applied rewrites79.6%
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k)
:precision binary64
(*
t_s
(if (<= t_m 3.6e-13)
(/
2.0
(*
(/ t_m l)
(*
(/
(fma
(* (fma 0.3333333333333333 (* t_m t_m) 1.0) k)
k
(* (* t_m t_m) 2.0))
l)
(* k k))))
(/ 2.0 (* (* (* (/ t_m l) (sin k)) t_m) (* (/ (* k t_m) l) 2.0))))))t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
double tmp;
if (t_m <= 3.6e-13) {
tmp = 2.0 / ((t_m / l) * ((fma((fma(0.3333333333333333, (t_m * t_m), 1.0) * k), k, ((t_m * t_m) * 2.0)) / l) * (k * k)));
} else {
tmp = 2.0 / ((((t_m / l) * sin(k)) * t_m) * (((k * t_m) / l) * 2.0));
}
return t_s * tmp;
}
t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k) tmp = 0.0 if (t_m <= 3.6e-13) tmp = Float64(2.0 / Float64(Float64(t_m / l) * Float64(Float64(fma(Float64(fma(0.3333333333333333, Float64(t_m * t_m), 1.0) * k), k, Float64(Float64(t_m * t_m) * 2.0)) / l) * Float64(k * k)))); else tmp = Float64(2.0 / Float64(Float64(Float64(Float64(t_m / l) * sin(k)) * t_m) * Float64(Float64(Float64(k * t_m) / l) * 2.0))); end return Float64(t_s * tmp) end
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k_] := N[(t$95$s * If[LessEqual[t$95$m, 3.6e-13], N[(2.0 / N[(N[(t$95$m / l), $MachinePrecision] * N[(N[(N[(N[(N[(0.3333333333333333 * N[(t$95$m * t$95$m), $MachinePrecision] + 1.0), $MachinePrecision] * k), $MachinePrecision] * k + N[(N[(t$95$m * t$95$m), $MachinePrecision] * 2.0), $MachinePrecision]), $MachinePrecision] / l), $MachinePrecision] * N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(N[(N[(t$95$m / l), $MachinePrecision] * N[Sin[k], $MachinePrecision]), $MachinePrecision] * t$95$m), $MachinePrecision] * N[(N[(N[(k * t$95$m), $MachinePrecision] / l), $MachinePrecision] * 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_m \leq 3.6 \cdot 10^{-13}:\\
\;\;\;\;\frac{2}{\frac{t\_m}{\ell} \cdot \left(\frac{\mathsf{fma}\left(\mathsf{fma}\left(0.3333333333333333, t\_m \cdot t\_m, 1\right) \cdot k, k, \left(t\_m \cdot t\_m\right) \cdot 2\right)}{\ell} \cdot \left(k \cdot k\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\left(\left(\frac{t\_m}{\ell} \cdot \sin k\right) \cdot t\_m\right) \cdot \left(\frac{k \cdot t\_m}{\ell} \cdot 2\right)}\\
\end{array}
\end{array}
if t < 3.5999999999999998e-13Initial program 56.4%
lift-*.f64N/A
lift-*.f64N/A
associate-*l*N/A
lift-*.f64N/A
associate-*l*N/A
lift-/.f64N/A
lift-pow.f64N/A
cube-multN/A
lift-*.f64N/A
times-fracN/A
associate-*l*N/A
lower-*.f64N/A
lower-/.f64N/A
lower-*.f64N/A
Applied rewrites68.7%
Taylor expanded in k around 0
*-commutativeN/A
lower-*.f64N/A
Applied rewrites70.8%
if 3.5999999999999998e-13 < t Initial program 68.8%
lift-*.f64N/A
*-commutativeN/A
lift-/.f64N/A
lift-pow.f64N/A
cube-multN/A
associate-/l*N/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
lift-*.f64N/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
lower-/.f6489.2
Applied rewrites89.2%
lift-*.f64N/A
lift-*.f64N/A
associate-*l*N/A
lift-*.f64N/A
lift-*.f64N/A
associate-*r*N/A
associate-*l*N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lift-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f64N/A
Applied rewrites95.2%
lift-*.f64N/A
lift-*.f64N/A
*-commutativeN/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f6495.1
Applied rewrites95.1%
Taylor expanded in k around 0
*-commutativeN/A
lower-*.f64N/A
lower-/.f64N/A
lower-*.f6484.9
Applied rewrites84.9%
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k)
:precision binary64
(*
t_s
(if (<= t_m 2200000000000.0)
(/
2.0
(*
(/ t_m l)
(*
(/
(fma
(* (fma 0.3333333333333333 (* t_m t_m) 1.0) k)
k
(* (* t_m t_m) 2.0))
l)
(* k k))))
(/ (/ (pow (/ l t_m) 2.0) (* k t_m)) k))))t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
double tmp;
if (t_m <= 2200000000000.0) {
tmp = 2.0 / ((t_m / l) * ((fma((fma(0.3333333333333333, (t_m * t_m), 1.0) * k), k, ((t_m * t_m) * 2.0)) / l) * (k * k)));
} else {
tmp = (pow((l / t_m), 2.0) / (k * t_m)) / k;
}
return t_s * tmp;
}
t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k) tmp = 0.0 if (t_m <= 2200000000000.0) tmp = Float64(2.0 / Float64(Float64(t_m / l) * Float64(Float64(fma(Float64(fma(0.3333333333333333, Float64(t_m * t_m), 1.0) * k), k, Float64(Float64(t_m * t_m) * 2.0)) / l) * Float64(k * k)))); else tmp = Float64(Float64((Float64(l / t_m) ^ 2.0) / Float64(k * t_m)) / k); end return Float64(t_s * tmp) end
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k_] := N[(t$95$s * If[LessEqual[t$95$m, 2200000000000.0], N[(2.0 / N[(N[(t$95$m / l), $MachinePrecision] * N[(N[(N[(N[(N[(0.3333333333333333 * N[(t$95$m * t$95$m), $MachinePrecision] + 1.0), $MachinePrecision] * k), $MachinePrecision] * k + N[(N[(t$95$m * t$95$m), $MachinePrecision] * 2.0), $MachinePrecision]), $MachinePrecision] / l), $MachinePrecision] * N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[Power[N[(l / t$95$m), $MachinePrecision], 2.0], $MachinePrecision] / N[(k * t$95$m), $MachinePrecision]), $MachinePrecision] / k), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_m \leq 2200000000000:\\
\;\;\;\;\frac{2}{\frac{t\_m}{\ell} \cdot \left(\frac{\mathsf{fma}\left(\mathsf{fma}\left(0.3333333333333333, t\_m \cdot t\_m, 1\right) \cdot k, k, \left(t\_m \cdot t\_m\right) \cdot 2\right)}{\ell} \cdot \left(k \cdot k\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{{\left(\frac{\ell}{t\_m}\right)}^{2}}{k \cdot t\_m}}{k}\\
\end{array}
\end{array}
if t < 2.2e12Initial program 56.5%
lift-*.f64N/A
lift-*.f64N/A
associate-*l*N/A
lift-*.f64N/A
associate-*l*N/A
lift-/.f64N/A
lift-pow.f64N/A
cube-multN/A
lift-*.f64N/A
times-fracN/A
associate-*l*N/A
lower-*.f64N/A
lower-/.f64N/A
lower-*.f64N/A
Applied rewrites68.4%
Taylor expanded in k around 0
*-commutativeN/A
lower-*.f64N/A
Applied rewrites70.4%
if 2.2e12 < t Initial program 70.1%
Taylor expanded in k around 0
unpow2N/A
*-commutativeN/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
lower-pow.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f6454.5
Applied rewrites54.5%
Applied rewrites54.5%
Applied rewrites61.9%
Applied rewrites79.6%
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k)
:precision binary64
(*
t_s
(if (<= k 6.5e-161)
(* (/ (/ l (* t_m t_m)) k) (/ l (* k t_m)))
(/
2.0
(*
(/ t_m l)
(*
(/
(fma
(* (fma 0.3333333333333333 (* t_m t_m) 1.0) k)
k
(* (* t_m t_m) 2.0))
l)
(* k k)))))))t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
double tmp;
if (k <= 6.5e-161) {
tmp = ((l / (t_m * t_m)) / k) * (l / (k * t_m));
} else {
tmp = 2.0 / ((t_m / l) * ((fma((fma(0.3333333333333333, (t_m * t_m), 1.0) * k), k, ((t_m * t_m) * 2.0)) / l) * (k * k)));
}
return t_s * tmp;
}
t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k) tmp = 0.0 if (k <= 6.5e-161) tmp = Float64(Float64(Float64(l / Float64(t_m * t_m)) / k) * Float64(l / Float64(k * t_m))); else tmp = Float64(2.0 / Float64(Float64(t_m / l) * Float64(Float64(fma(Float64(fma(0.3333333333333333, Float64(t_m * t_m), 1.0) * k), k, Float64(Float64(t_m * t_m) * 2.0)) / l) * Float64(k * k)))); end return Float64(t_s * tmp) end
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k_] := N[(t$95$s * If[LessEqual[k, 6.5e-161], N[(N[(N[(l / N[(t$95$m * t$95$m), $MachinePrecision]), $MachinePrecision] / k), $MachinePrecision] * N[(l / N[(k * t$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(t$95$m / l), $MachinePrecision] * N[(N[(N[(N[(N[(0.3333333333333333 * N[(t$95$m * t$95$m), $MachinePrecision] + 1.0), $MachinePrecision] * k), $MachinePrecision] * k + N[(N[(t$95$m * t$95$m), $MachinePrecision] * 2.0), $MachinePrecision]), $MachinePrecision] / l), $MachinePrecision] * N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k \leq 6.5 \cdot 10^{-161}:\\
\;\;\;\;\frac{\frac{\ell}{t\_m \cdot t\_m}}{k} \cdot \frac{\ell}{k \cdot t\_m}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\frac{t\_m}{\ell} \cdot \left(\frac{\mathsf{fma}\left(\mathsf{fma}\left(0.3333333333333333, t\_m \cdot t\_m, 1\right) \cdot k, k, \left(t\_m \cdot t\_m\right) \cdot 2\right)}{\ell} \cdot \left(k \cdot k\right)\right)}\\
\end{array}
\end{array}
if k < 6.50000000000000008e-161Initial program 59.6%
Taylor expanded in k around 0
unpow2N/A
*-commutativeN/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
lower-pow.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f6456.7
Applied rewrites56.7%
Applied rewrites56.7%
Applied rewrites58.3%
Applied rewrites68.3%
if 6.50000000000000008e-161 < k Initial program 59.3%
lift-*.f64N/A
lift-*.f64N/A
associate-*l*N/A
lift-*.f64N/A
associate-*l*N/A
lift-/.f64N/A
lift-pow.f64N/A
cube-multN/A
lift-*.f64N/A
times-fracN/A
associate-*l*N/A
lower-*.f64N/A
lower-/.f64N/A
lower-*.f64N/A
Applied rewrites72.8%
Taylor expanded in k around 0
*-commutativeN/A
lower-*.f64N/A
Applied rewrites74.2%
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k)
:precision binary64
(*
t_s
(if (<= t_m 9.6e-99)
(* (/ 2.0 (* (* k k) (* k k))) (/ (* l l) t_m))
(if (<= t_m 5.5e+76)
(* (/ l (* (* t_m t_m) t_m)) (/ (/ l k) k))
(/ (* (/ l (* t_m t_m)) l) (* (* k t_m) k))))))t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
double tmp;
if (t_m <= 9.6e-99) {
tmp = (2.0 / ((k * k) * (k * k))) * ((l * l) / t_m);
} else if (t_m <= 5.5e+76) {
tmp = (l / ((t_m * t_m) * t_m)) * ((l / k) / k);
} else {
tmp = ((l / (t_m * t_m)) * l) / ((k * t_m) * k);
}
return t_s * tmp;
}
t\_m = private
t\_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(t_s, t_m, l, k)
use fmin_fmax_functions
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k
real(8) :: tmp
if (t_m <= 9.6d-99) then
tmp = (2.0d0 / ((k * k) * (k * k))) * ((l * l) / t_m)
else if (t_m <= 5.5d+76) then
tmp = (l / ((t_m * t_m) * t_m)) * ((l / k) / k)
else
tmp = ((l / (t_m * t_m)) * l) / ((k * t_m) * k)
end if
code = t_s * tmp
end function
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k) {
double tmp;
if (t_m <= 9.6e-99) {
tmp = (2.0 / ((k * k) * (k * k))) * ((l * l) / t_m);
} else if (t_m <= 5.5e+76) {
tmp = (l / ((t_m * t_m) * t_m)) * ((l / k) / k);
} else {
tmp = ((l / (t_m * t_m)) * l) / ((k * t_m) * k);
}
return t_s * tmp;
}
t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k): tmp = 0 if t_m <= 9.6e-99: tmp = (2.0 / ((k * k) * (k * k))) * ((l * l) / t_m) elif t_m <= 5.5e+76: tmp = (l / ((t_m * t_m) * t_m)) * ((l / k) / k) else: tmp = ((l / (t_m * t_m)) * l) / ((k * t_m) * k) return t_s * tmp
t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k) tmp = 0.0 if (t_m <= 9.6e-99) tmp = Float64(Float64(2.0 / Float64(Float64(k * k) * Float64(k * k))) * Float64(Float64(l * l) / t_m)); elseif (t_m <= 5.5e+76) tmp = Float64(Float64(l / Float64(Float64(t_m * t_m) * t_m)) * Float64(Float64(l / k) / k)); else tmp = Float64(Float64(Float64(l / Float64(t_m * t_m)) * l) / Float64(Float64(k * t_m) * k)); end return Float64(t_s * tmp) end
t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k) tmp = 0.0; if (t_m <= 9.6e-99) tmp = (2.0 / ((k * k) * (k * k))) * ((l * l) / t_m); elseif (t_m <= 5.5e+76) tmp = (l / ((t_m * t_m) * t_m)) * ((l / k) / k); else tmp = ((l / (t_m * t_m)) * l) / ((k * t_m) * k); end tmp_2 = t_s * tmp; end
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k_] := N[(t$95$s * If[LessEqual[t$95$m, 9.6e-99], N[(N[(2.0 / N[(N[(k * k), $MachinePrecision] * N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[(l * l), $MachinePrecision] / t$95$m), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$m, 5.5e+76], N[(N[(l / N[(N[(t$95$m * t$95$m), $MachinePrecision] * t$95$m), $MachinePrecision]), $MachinePrecision] * N[(N[(l / k), $MachinePrecision] / k), $MachinePrecision]), $MachinePrecision], N[(N[(N[(l / N[(t$95$m * t$95$m), $MachinePrecision]), $MachinePrecision] * l), $MachinePrecision] / N[(N[(k * t$95$m), $MachinePrecision] * k), $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]
\begin{array}{l}
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_m \leq 9.6 \cdot 10^{-99}:\\
\;\;\;\;\frac{2}{\left(k \cdot k\right) \cdot \left(k \cdot k\right)} \cdot \frac{\ell \cdot \ell}{t\_m}\\
\mathbf{elif}\;t\_m \leq 5.5 \cdot 10^{+76}:\\
\;\;\;\;\frac{\ell}{\left(t\_m \cdot t\_m\right) \cdot t\_m} \cdot \frac{\frac{\ell}{k}}{k}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{\ell}{t\_m \cdot t\_m} \cdot \ell}{\left(k \cdot t\_m\right) \cdot k}\\
\end{array}
\end{array}
if t < 9.6000000000000002e-99Initial program 54.6%
Applied rewrites48.0%
Taylor expanded in k around 0
associate-*r/N/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
lower-pow.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f6457.3
Applied rewrites57.3%
Applied rewrites57.3%
if 9.6000000000000002e-99 < t < 5.5000000000000001e76Initial program 70.9%
Taylor expanded in k around 0
unpow2N/A
*-commutativeN/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
lower-pow.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f6463.7
Applied rewrites63.7%
Applied rewrites63.7%
Applied rewrites70.5%
if 5.5000000000000001e76 < t Initial program 67.9%
Taylor expanded in k around 0
unpow2N/A
*-commutativeN/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
lower-pow.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f6457.2
Applied rewrites57.2%
Applied rewrites57.2%
Applied rewrites66.3%
Applied rewrites79.6%
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k)
:precision binary64
(*
t_s
(if (<= t_m 2.3e-80)
(/ 2.0 (* (* (* k k) t_m) (* (/ k l) (/ k l))))
(* (/ (/ l (* t_m t_m)) k) (/ l (* k t_m))))))t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
double tmp;
if (t_m <= 2.3e-80) {
tmp = 2.0 / (((k * k) * t_m) * ((k / l) * (k / l)));
} else {
tmp = ((l / (t_m * t_m)) / k) * (l / (k * t_m));
}
return t_s * tmp;
}
t\_m = private
t\_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(t_s, t_m, l, k)
use fmin_fmax_functions
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k
real(8) :: tmp
if (t_m <= 2.3d-80) then
tmp = 2.0d0 / (((k * k) * t_m) * ((k / l) * (k / l)))
else
tmp = ((l / (t_m * t_m)) / k) * (l / (k * t_m))
end if
code = t_s * tmp
end function
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k) {
double tmp;
if (t_m <= 2.3e-80) {
tmp = 2.0 / (((k * k) * t_m) * ((k / l) * (k / l)));
} else {
tmp = ((l / (t_m * t_m)) / k) * (l / (k * t_m));
}
return t_s * tmp;
}
t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k): tmp = 0 if t_m <= 2.3e-80: tmp = 2.0 / (((k * k) * t_m) * ((k / l) * (k / l))) else: tmp = ((l / (t_m * t_m)) / k) * (l / (k * t_m)) return t_s * tmp
t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k) tmp = 0.0 if (t_m <= 2.3e-80) tmp = Float64(2.0 / Float64(Float64(Float64(k * k) * t_m) * Float64(Float64(k / l) * Float64(k / l)))); else tmp = Float64(Float64(Float64(l / Float64(t_m * t_m)) / k) * Float64(l / Float64(k * t_m))); end return Float64(t_s * tmp) end
t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k) tmp = 0.0; if (t_m <= 2.3e-80) tmp = 2.0 / (((k * k) * t_m) * ((k / l) * (k / l))); else tmp = ((l / (t_m * t_m)) / k) * (l / (k * t_m)); end tmp_2 = t_s * tmp; end
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k_] := N[(t$95$s * If[LessEqual[t$95$m, 2.3e-80], N[(2.0 / N[(N[(N[(k * k), $MachinePrecision] * t$95$m), $MachinePrecision] * N[(N[(k / l), $MachinePrecision] * N[(k / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(l / N[(t$95$m * t$95$m), $MachinePrecision]), $MachinePrecision] / k), $MachinePrecision] * N[(l / N[(k * t$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_m \leq 2.3 \cdot 10^{-80}:\\
\;\;\;\;\frac{2}{\left(\left(k \cdot k\right) \cdot t\_m\right) \cdot \left(\frac{k}{\ell} \cdot \frac{k}{\ell}\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{\ell}{t\_m \cdot t\_m}}{k} \cdot \frac{\ell}{k \cdot t\_m}\\
\end{array}
\end{array}
if t < 2.2999999999999998e-80Initial program 54.3%
Taylor expanded in t around 0
associate-*r*N/A
associate-/l*N/A
lower-*.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
lower-/.f64N/A
lower-pow.f64N/A
lower-sin.f64N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
lower-cos.f6469.9
Applied rewrites69.9%
Taylor expanded in k around 0
Applied rewrites62.4%
if 2.2999999999999998e-80 < t Initial program 71.0%
Taylor expanded in k around 0
unpow2N/A
*-commutativeN/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
lower-pow.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f6460.0
Applied rewrites60.0%
Applied rewrites60.0%
Applied rewrites64.7%
Applied rewrites79.7%
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k)
:precision binary64
(*
t_s
(if (<= t_m 2.2e-80)
(/ (/ (* (/ (* l l) t_m) 2.0) (* k k)) (* k k))
(* (/ (/ l (* t_m t_m)) k) (/ l (* k t_m))))))t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
double tmp;
if (t_m <= 2.2e-80) {
tmp = ((((l * l) / t_m) * 2.0) / (k * k)) / (k * k);
} else {
tmp = ((l / (t_m * t_m)) / k) * (l / (k * t_m));
}
return t_s * tmp;
}
t\_m = private
t\_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(t_s, t_m, l, k)
use fmin_fmax_functions
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k
real(8) :: tmp
if (t_m <= 2.2d-80) then
tmp = ((((l * l) / t_m) * 2.0d0) / (k * k)) / (k * k)
else
tmp = ((l / (t_m * t_m)) / k) * (l / (k * t_m))
end if
code = t_s * tmp
end function
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k) {
double tmp;
if (t_m <= 2.2e-80) {
tmp = ((((l * l) / t_m) * 2.0) / (k * k)) / (k * k);
} else {
tmp = ((l / (t_m * t_m)) / k) * (l / (k * t_m));
}
return t_s * tmp;
}
t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k): tmp = 0 if t_m <= 2.2e-80: tmp = ((((l * l) / t_m) * 2.0) / (k * k)) / (k * k) else: tmp = ((l / (t_m * t_m)) / k) * (l / (k * t_m)) return t_s * tmp
t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k) tmp = 0.0 if (t_m <= 2.2e-80) tmp = Float64(Float64(Float64(Float64(Float64(l * l) / t_m) * 2.0) / Float64(k * k)) / Float64(k * k)); else tmp = Float64(Float64(Float64(l / Float64(t_m * t_m)) / k) * Float64(l / Float64(k * t_m))); end return Float64(t_s * tmp) end
t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k) tmp = 0.0; if (t_m <= 2.2e-80) tmp = ((((l * l) / t_m) * 2.0) / (k * k)) / (k * k); else tmp = ((l / (t_m * t_m)) / k) * (l / (k * t_m)); end tmp_2 = t_s * tmp; end
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k_] := N[(t$95$s * If[LessEqual[t$95$m, 2.2e-80], N[(N[(N[(N[(N[(l * l), $MachinePrecision] / t$95$m), $MachinePrecision] * 2.0), $MachinePrecision] / N[(k * k), $MachinePrecision]), $MachinePrecision] / N[(k * k), $MachinePrecision]), $MachinePrecision], N[(N[(N[(l / N[(t$95$m * t$95$m), $MachinePrecision]), $MachinePrecision] / k), $MachinePrecision] * N[(l / N[(k * t$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_m \leq 2.2 \cdot 10^{-80}:\\
\;\;\;\;\frac{\frac{\frac{\ell \cdot \ell}{t\_m} \cdot 2}{k \cdot k}}{k \cdot k}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{\ell}{t\_m \cdot t\_m}}{k} \cdot \frac{\ell}{k \cdot t\_m}\\
\end{array}
\end{array}
if t < 2.2000000000000001e-80Initial program 54.3%
Applied rewrites48.3%
Taylor expanded in k around 0
associate-*r/N/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
lower-pow.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f6456.9
Applied rewrites56.9%
Applied rewrites61.0%
if 2.2000000000000001e-80 < t Initial program 71.0%
Taylor expanded in k around 0
unpow2N/A
*-commutativeN/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
lower-pow.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f6460.0
Applied rewrites60.0%
Applied rewrites60.0%
Applied rewrites64.7%
Applied rewrites79.7%
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k)
:precision binary64
(*
t_s
(if (<= t_m 5.2e-88)
(* (/ 2.0 (* (* k k) (* k k))) (/ (* l l) t_m))
(* (/ (/ l (* t_m t_m)) k) (/ l (* k t_m))))))t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
double tmp;
if (t_m <= 5.2e-88) {
tmp = (2.0 / ((k * k) * (k * k))) * ((l * l) / t_m);
} else {
tmp = ((l / (t_m * t_m)) / k) * (l / (k * t_m));
}
return t_s * tmp;
}
t\_m = private
t\_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(t_s, t_m, l, k)
use fmin_fmax_functions
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k
real(8) :: tmp
if (t_m <= 5.2d-88) then
tmp = (2.0d0 / ((k * k) * (k * k))) * ((l * l) / t_m)
else
tmp = ((l / (t_m * t_m)) / k) * (l / (k * t_m))
end if
code = t_s * tmp
end function
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k) {
double tmp;
if (t_m <= 5.2e-88) {
tmp = (2.0 / ((k * k) * (k * k))) * ((l * l) / t_m);
} else {
tmp = ((l / (t_m * t_m)) / k) * (l / (k * t_m));
}
return t_s * tmp;
}
t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k): tmp = 0 if t_m <= 5.2e-88: tmp = (2.0 / ((k * k) * (k * k))) * ((l * l) / t_m) else: tmp = ((l / (t_m * t_m)) / k) * (l / (k * t_m)) return t_s * tmp
t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k) tmp = 0.0 if (t_m <= 5.2e-88) tmp = Float64(Float64(2.0 / Float64(Float64(k * k) * Float64(k * k))) * Float64(Float64(l * l) / t_m)); else tmp = Float64(Float64(Float64(l / Float64(t_m * t_m)) / k) * Float64(l / Float64(k * t_m))); end return Float64(t_s * tmp) end
t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k) tmp = 0.0; if (t_m <= 5.2e-88) tmp = (2.0 / ((k * k) * (k * k))) * ((l * l) / t_m); else tmp = ((l / (t_m * t_m)) / k) * (l / (k * t_m)); end tmp_2 = t_s * tmp; end
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k_] := N[(t$95$s * If[LessEqual[t$95$m, 5.2e-88], N[(N[(2.0 / N[(N[(k * k), $MachinePrecision] * N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[(l * l), $MachinePrecision] / t$95$m), $MachinePrecision]), $MachinePrecision], N[(N[(N[(l / N[(t$95$m * t$95$m), $MachinePrecision]), $MachinePrecision] / k), $MachinePrecision] * N[(l / N[(k * t$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_m \leq 5.2 \cdot 10^{-88}:\\
\;\;\;\;\frac{2}{\left(k \cdot k\right) \cdot \left(k \cdot k\right)} \cdot \frac{\ell \cdot \ell}{t\_m}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{\ell}{t\_m \cdot t\_m}}{k} \cdot \frac{\ell}{k \cdot t\_m}\\
\end{array}
\end{array}
if t < 5.20000000000000027e-88Initial program 54.3%
Applied rewrites48.3%
Taylor expanded in k around 0
associate-*r/N/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
lower-pow.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f6457.0
Applied rewrites57.0%
Applied rewrites57.0%
if 5.20000000000000027e-88 < t Initial program 70.5%
Taylor expanded in k around 0
unpow2N/A
*-commutativeN/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
lower-pow.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f6461.0
Applied rewrites61.0%
Applied rewrites61.0%
Applied rewrites65.6%
Applied rewrites80.2%
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k)
:precision binary64
(*
t_s
(if (<= k 1e-137)
(/ (* (/ l (* t_m t_m)) l) (* (* k t_m) k))
(/ (* (/ l t_m) l) (* t_m (* (* k k) t_m))))))t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
double tmp;
if (k <= 1e-137) {
tmp = ((l / (t_m * t_m)) * l) / ((k * t_m) * k);
} else {
tmp = ((l / t_m) * l) / (t_m * ((k * k) * t_m));
}
return t_s * tmp;
}
t\_m = private
t\_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(t_s, t_m, l, k)
use fmin_fmax_functions
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k
real(8) :: tmp
if (k <= 1d-137) then
tmp = ((l / (t_m * t_m)) * l) / ((k * t_m) * k)
else
tmp = ((l / t_m) * l) / (t_m * ((k * k) * t_m))
end if
code = t_s * tmp
end function
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k) {
double tmp;
if (k <= 1e-137) {
tmp = ((l / (t_m * t_m)) * l) / ((k * t_m) * k);
} else {
tmp = ((l / t_m) * l) / (t_m * ((k * k) * t_m));
}
return t_s * tmp;
}
t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k): tmp = 0 if k <= 1e-137: tmp = ((l / (t_m * t_m)) * l) / ((k * t_m) * k) else: tmp = ((l / t_m) * l) / (t_m * ((k * k) * t_m)) return t_s * tmp
t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k) tmp = 0.0 if (k <= 1e-137) tmp = Float64(Float64(Float64(l / Float64(t_m * t_m)) * l) / Float64(Float64(k * t_m) * k)); else tmp = Float64(Float64(Float64(l / t_m) * l) / Float64(t_m * Float64(Float64(k * k) * t_m))); end return Float64(t_s * tmp) end
t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k) tmp = 0.0; if (k <= 1e-137) tmp = ((l / (t_m * t_m)) * l) / ((k * t_m) * k); else tmp = ((l / t_m) * l) / (t_m * ((k * k) * t_m)); end tmp_2 = t_s * tmp; end
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k_] := N[(t$95$s * If[LessEqual[k, 1e-137], N[(N[(N[(l / N[(t$95$m * t$95$m), $MachinePrecision]), $MachinePrecision] * l), $MachinePrecision] / N[(N[(k * t$95$m), $MachinePrecision] * k), $MachinePrecision]), $MachinePrecision], N[(N[(N[(l / t$95$m), $MachinePrecision] * l), $MachinePrecision] / N[(t$95$m * N[(N[(k * k), $MachinePrecision] * t$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k \leq 10^{-137}:\\
\;\;\;\;\frac{\frac{\ell}{t\_m \cdot t\_m} \cdot \ell}{\left(k \cdot t\_m\right) \cdot k}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{\ell}{t\_m} \cdot \ell}{t\_m \cdot \left(\left(k \cdot k\right) \cdot t\_m\right)}\\
\end{array}
\end{array}
if k < 9.99999999999999978e-138Initial program 60.4%
Taylor expanded in k around 0
unpow2N/A
*-commutativeN/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
lower-pow.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f6457.5
Applied rewrites57.5%
Applied rewrites57.5%
Applied rewrites59.1%
Applied rewrites64.3%
if 9.99999999999999978e-138 < k Initial program 57.9%
Taylor expanded in k around 0
unpow2N/A
*-commutativeN/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
lower-pow.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f6459.1
Applied rewrites59.1%
Applied rewrites59.1%
Applied rewrites63.1%
Applied rewrites66.3%
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k)
:precision binary64
(let* ((t_2 (* (* k k) t_m)))
(*
t_s
(if (<= t_m 4e-46)
(/ (* (* l l) 2.0) (* t_2 (* k k)))
(* l (/ (/ l (* t_m t_m)) t_2))))))t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
double t_2 = (k * k) * t_m;
double tmp;
if (t_m <= 4e-46) {
tmp = ((l * l) * 2.0) / (t_2 * (k * k));
} else {
tmp = l * ((l / (t_m * t_m)) / t_2);
}
return t_s * tmp;
}
t\_m = private
t\_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(t_s, t_m, l, k)
use fmin_fmax_functions
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k
real(8) :: t_2
real(8) :: tmp
t_2 = (k * k) * t_m
if (t_m <= 4d-46) then
tmp = ((l * l) * 2.0d0) / (t_2 * (k * k))
else
tmp = l * ((l / (t_m * t_m)) / t_2)
end if
code = t_s * tmp
end function
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k) {
double t_2 = (k * k) * t_m;
double tmp;
if (t_m <= 4e-46) {
tmp = ((l * l) * 2.0) / (t_2 * (k * k));
} else {
tmp = l * ((l / (t_m * t_m)) / t_2);
}
return t_s * tmp;
}
t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k): t_2 = (k * k) * t_m tmp = 0 if t_m <= 4e-46: tmp = ((l * l) * 2.0) / (t_2 * (k * k)) else: tmp = l * ((l / (t_m * t_m)) / t_2) return t_s * tmp
t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k) t_2 = Float64(Float64(k * k) * t_m) tmp = 0.0 if (t_m <= 4e-46) tmp = Float64(Float64(Float64(l * l) * 2.0) / Float64(t_2 * Float64(k * k))); else tmp = Float64(l * Float64(Float64(l / Float64(t_m * t_m)) / t_2)); end return Float64(t_s * tmp) end
t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k) t_2 = (k * k) * t_m; tmp = 0.0; if (t_m <= 4e-46) tmp = ((l * l) * 2.0) / (t_2 * (k * k)); else tmp = l * ((l / (t_m * t_m)) / t_2); end tmp_2 = t_s * tmp; end
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k_] := Block[{t$95$2 = N[(N[(k * k), $MachinePrecision] * t$95$m), $MachinePrecision]}, N[(t$95$s * If[LessEqual[t$95$m, 4e-46], N[(N[(N[(l * l), $MachinePrecision] * 2.0), $MachinePrecision] / N[(t$95$2 * N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(l * N[(N[(l / N[(t$95$m * t$95$m), $MachinePrecision]), $MachinePrecision] / t$95$2), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]]
\begin{array}{l}
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
\begin{array}{l}
t_2 := \left(k \cdot k\right) \cdot t\_m\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_m \leq 4 \cdot 10^{-46}:\\
\;\;\;\;\frac{\left(\ell \cdot \ell\right) \cdot 2}{t\_2 \cdot \left(k \cdot k\right)}\\
\mathbf{else}:\\
\;\;\;\;\ell \cdot \frac{\frac{\ell}{t\_m \cdot t\_m}}{t\_2}\\
\end{array}
\end{array}
\end{array}
if t < 4.00000000000000009e-46Initial program 54.8%
Applied rewrites49.2%
Taylor expanded in k around 0
associate-*r/N/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
lower-pow.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f6457.3
Applied rewrites57.3%
Applied rewrites57.6%
Applied rewrites59.3%
if 4.00000000000000009e-46 < t Initial program 72.0%
Taylor expanded in k around 0
unpow2N/A
*-commutativeN/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
lower-pow.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f6459.4
Applied rewrites59.4%
Applied rewrites59.4%
Applied rewrites65.3%
Applied rewrites65.4%
t\_m = (fabs.f64 t) t\_s = (copysign.f64 #s(literal 1 binary64) t) (FPCore (t_s t_m l k) :precision binary64 (* t_s (/ (* (/ l t_m) l) (* t_m (* (* k k) t_m)))))
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
return t_s * (((l / t_m) * l) / (t_m * ((k * k) * t_m)));
}
t\_m = private
t\_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(t_s, t_m, l, k)
use fmin_fmax_functions
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k
code = t_s * (((l / t_m) * l) / (t_m * ((k * k) * t_m)))
end function
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k) {
return t_s * (((l / t_m) * l) / (t_m * ((k * k) * t_m)));
}
t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k): return t_s * (((l / t_m) * l) / (t_m * ((k * k) * t_m)))
t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k) return Float64(t_s * Float64(Float64(Float64(l / t_m) * l) / Float64(t_m * Float64(Float64(k * k) * t_m)))) end
t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp = code(t_s, t_m, l, k) tmp = t_s * (((l / t_m) * l) / (t_m * ((k * k) * t_m))); end
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k_] := N[(t$95$s * N[(N[(N[(l / t$95$m), $MachinePrecision] * l), $MachinePrecision] / N[(t$95$m * N[(N[(k * k), $MachinePrecision] * t$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \frac{\frac{\ell}{t\_m} \cdot \ell}{t\_m \cdot \left(\left(k \cdot k\right) \cdot t\_m\right)}
\end{array}
Initial program 59.5%
Taylor expanded in k around 0
unpow2N/A
*-commutativeN/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
lower-pow.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f6458.1
Applied rewrites58.1%
Applied rewrites58.1%
Applied rewrites60.6%
Applied rewrites64.1%
t\_m = (fabs.f64 t) t\_s = (copysign.f64 #s(literal 1 binary64) t) (FPCore (t_s t_m l k) :precision binary64 (* t_s (/ (* l l) (* (* t_m t_m) (* (* k k) t_m)))))
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
return t_s * ((l * l) / ((t_m * t_m) * ((k * k) * t_m)));
}
t\_m = private
t\_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(t_s, t_m, l, k)
use fmin_fmax_functions
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k
code = t_s * ((l * l) / ((t_m * t_m) * ((k * k) * t_m)))
end function
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k) {
return t_s * ((l * l) / ((t_m * t_m) * ((k * k) * t_m)));
}
t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k): return t_s * ((l * l) / ((t_m * t_m) * ((k * k) * t_m)))
t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k) return Float64(t_s * Float64(Float64(l * l) / Float64(Float64(t_m * t_m) * Float64(Float64(k * k) * t_m)))) end
t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp = code(t_s, t_m, l, k) tmp = t_s * ((l * l) / ((t_m * t_m) * ((k * k) * t_m))); end
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k_] := N[(t$95$s * N[(N[(l * l), $MachinePrecision] / N[(N[(t$95$m * t$95$m), $MachinePrecision] * N[(N[(k * k), $MachinePrecision] * t$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \frac{\ell \cdot \ell}{\left(t\_m \cdot t\_m\right) \cdot \left(\left(k \cdot k\right) \cdot t\_m\right)}
\end{array}
Initial program 59.5%
Taylor expanded in k around 0
unpow2N/A
*-commutativeN/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
lower-pow.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f6458.1
Applied rewrites58.1%
Applied rewrites58.1%
Applied rewrites60.6%
Applied rewrites57.3%
herbie shell --seed 2024347
(FPCore (t l k)
:name "Toniolo and Linder, Equation (10+)"
:precision binary64
(/ 2.0 (* (* (* (/ (pow t 3.0) (* l l)) (sin k)) (tan k)) (+ (+ 1.0 (pow (/ k t) 2.0)) 1.0))))