
(FPCore (K m n M l) :precision binary64 (* (cos (- (/ (* K (+ m n)) 2.0) M)) (exp (- (- (pow (- (/ (+ m n) 2.0) M) 2.0)) (- l (fabs (- m n)))))))
double code(double K, double m, double n, double M, double l) {
return cos((((K * (m + n)) / 2.0) - M)) * exp((-pow((((m + n) / 2.0) - M), 2.0) - (l - fabs((m - n)))));
}
real(8) function code(k, m, n, m_1, l)
real(8), intent (in) :: k
real(8), intent (in) :: m
real(8), intent (in) :: n
real(8), intent (in) :: m_1
real(8), intent (in) :: l
code = cos((((k * (m + n)) / 2.0d0) - m_1)) * exp((-((((m + n) / 2.0d0) - m_1) ** 2.0d0) - (l - abs((m - n)))))
end function
public static double code(double K, double m, double n, double M, double l) {
return Math.cos((((K * (m + n)) / 2.0) - M)) * Math.exp((-Math.pow((((m + n) / 2.0) - M), 2.0) - (l - Math.abs((m - n)))));
}
def code(K, m, n, M, l): return math.cos((((K * (m + n)) / 2.0) - M)) * math.exp((-math.pow((((m + n) / 2.0) - M), 2.0) - (l - math.fabs((m - n)))))
function code(K, m, n, M, l) return Float64(cos(Float64(Float64(Float64(K * Float64(m + n)) / 2.0) - M)) * exp(Float64(Float64(-(Float64(Float64(Float64(m + n) / 2.0) - M) ^ 2.0)) - Float64(l - abs(Float64(m - n)))))) end
function tmp = code(K, m, n, M, l) tmp = cos((((K * (m + n)) / 2.0) - M)) * exp((-((((m + n) / 2.0) - M) ^ 2.0) - (l - abs((m - n))))); end
code[K_, m_, n_, M_, l_] := N[(N[Cos[N[(N[(N[(K * N[(m + n), $MachinePrecision]), $MachinePrecision] / 2.0), $MachinePrecision] - M), $MachinePrecision]], $MachinePrecision] * N[Exp[N[((-N[Power[N[(N[(N[(m + n), $MachinePrecision] / 2.0), $MachinePrecision] - M), $MachinePrecision], 2.0], $MachinePrecision]) - N[(l - N[Abs[N[(m - n), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\cos \left(\frac{K \cdot \left(m + n\right)}{2} - M\right) \cdot e^{\left(-{\left(\frac{m + n}{2} - M\right)}^{2}\right) - \left(\ell - \left|m - n\right|\right)}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 11 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (K m n M l) :precision binary64 (* (cos (- (/ (* K (+ m n)) 2.0) M)) (exp (- (- (pow (- (/ (+ m n) 2.0) M) 2.0)) (- l (fabs (- m n)))))))
double code(double K, double m, double n, double M, double l) {
return cos((((K * (m + n)) / 2.0) - M)) * exp((-pow((((m + n) / 2.0) - M), 2.0) - (l - fabs((m - n)))));
}
real(8) function code(k, m, n, m_1, l)
real(8), intent (in) :: k
real(8), intent (in) :: m
real(8), intent (in) :: n
real(8), intent (in) :: m_1
real(8), intent (in) :: l
code = cos((((k * (m + n)) / 2.0d0) - m_1)) * exp((-((((m + n) / 2.0d0) - m_1) ** 2.0d0) - (l - abs((m - n)))))
end function
public static double code(double K, double m, double n, double M, double l) {
return Math.cos((((K * (m + n)) / 2.0) - M)) * Math.exp((-Math.pow((((m + n) / 2.0) - M), 2.0) - (l - Math.abs((m - n)))));
}
def code(K, m, n, M, l): return math.cos((((K * (m + n)) / 2.0) - M)) * math.exp((-math.pow((((m + n) / 2.0) - M), 2.0) - (l - math.fabs((m - n)))))
function code(K, m, n, M, l) return Float64(cos(Float64(Float64(Float64(K * Float64(m + n)) / 2.0) - M)) * exp(Float64(Float64(-(Float64(Float64(Float64(m + n) / 2.0) - M) ^ 2.0)) - Float64(l - abs(Float64(m - n)))))) end
function tmp = code(K, m, n, M, l) tmp = cos((((K * (m + n)) / 2.0) - M)) * exp((-((((m + n) / 2.0) - M) ^ 2.0) - (l - abs((m - n))))); end
code[K_, m_, n_, M_, l_] := N[(N[Cos[N[(N[(N[(K * N[(m + n), $MachinePrecision]), $MachinePrecision] / 2.0), $MachinePrecision] - M), $MachinePrecision]], $MachinePrecision] * N[Exp[N[((-N[Power[N[(N[(N[(m + n), $MachinePrecision] / 2.0), $MachinePrecision] - M), $MachinePrecision], 2.0], $MachinePrecision]) - N[(l - N[Abs[N[(m - n), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\cos \left(\frac{K \cdot \left(m + n\right)}{2} - M\right) \cdot e^{\left(-{\left(\frac{m + n}{2} - M\right)}^{2}\right) - \left(\ell - \left|m - n\right|\right)}
\end{array}
(FPCore (K m n M l) :precision binary64 (* (exp (+ (- (* (+ M (* -0.5 (+ n m))) (- (* 0.5 (+ n m)) M)) l) (fabs (- n m)))) (cos M)))
double code(double K, double m, double n, double M, double l) {
return exp(((((M + (-0.5 * (n + m))) * ((0.5 * (n + m)) - M)) - l) + fabs((n - m)))) * cos(M);
}
real(8) function code(k, m, n, m_1, l)
real(8), intent (in) :: k
real(8), intent (in) :: m
real(8), intent (in) :: n
real(8), intent (in) :: m_1
real(8), intent (in) :: l
code = exp(((((m_1 + ((-0.5d0) * (n + m))) * ((0.5d0 * (n + m)) - m_1)) - l) + abs((n - m)))) * cos(m_1)
end function
public static double code(double K, double m, double n, double M, double l) {
return Math.exp(((((M + (-0.5 * (n + m))) * ((0.5 * (n + m)) - M)) - l) + Math.abs((n - m)))) * Math.cos(M);
}
def code(K, m, n, M, l): return math.exp(((((M + (-0.5 * (n + m))) * ((0.5 * (n + m)) - M)) - l) + math.fabs((n - m)))) * math.cos(M)
function code(K, m, n, M, l) return Float64(exp(Float64(Float64(Float64(Float64(M + Float64(-0.5 * Float64(n + m))) * Float64(Float64(0.5 * Float64(n + m)) - M)) - l) + abs(Float64(n - m)))) * cos(M)) end
function tmp = code(K, m, n, M, l) tmp = exp(((((M + (-0.5 * (n + m))) * ((0.5 * (n + m)) - M)) - l) + abs((n - m)))) * cos(M); end
code[K_, m_, n_, M_, l_] := N[(N[Exp[N[(N[(N[(N[(M + N[(-0.5 * N[(n + m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[(0.5 * N[(n + m), $MachinePrecision]), $MachinePrecision] - M), $MachinePrecision]), $MachinePrecision] - l), $MachinePrecision] + N[Abs[N[(n - m), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * N[Cos[M], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
e^{\left(\left(M + -0.5 \cdot \left(n + m\right)\right) \cdot \left(0.5 \cdot \left(n + m\right) - M\right) - \ell\right) + \left|n - m\right|} \cdot \cos M
\end{array}
Initial program 78.2%
*-lowering-*.f64N/A
sub-negN/A
+-commutativeN/A
neg-sub0N/A
associate-+l-N/A
sub0-negN/A
cos-negN/A
cos-lowering-cos.f64N/A
sub-negN/A
+-lowering-+.f64N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
metadata-evalN/A
exp-lowering-exp.f64N/A
Simplified78.2%
Taylor expanded in K around 0
*-commutativeN/A
*-lowering-*.f64N/A
Simplified98.9%
Final simplification98.9%
(FPCore (K m n M l)
:precision binary64
(let* ((t_0 (* (cos M) (exp (- 0.0 (* M M))))))
(if (<= M -125000.0)
t_0
(if (<= M 13500.0)
(exp (+ (fabs (- n m)) (- (* -0.25 (* (+ n m) (+ n m))) l)))
t_0))))
double code(double K, double m, double n, double M, double l) {
double t_0 = cos(M) * exp((0.0 - (M * M)));
double tmp;
if (M <= -125000.0) {
tmp = t_0;
} else if (M <= 13500.0) {
tmp = exp((fabs((n - m)) + ((-0.25 * ((n + m) * (n + m))) - l)));
} else {
tmp = t_0;
}
return tmp;
}
real(8) function code(k, m, n, m_1, l)
real(8), intent (in) :: k
real(8), intent (in) :: m
real(8), intent (in) :: n
real(8), intent (in) :: m_1
real(8), intent (in) :: l
real(8) :: t_0
real(8) :: tmp
t_0 = cos(m_1) * exp((0.0d0 - (m_1 * m_1)))
if (m_1 <= (-125000.0d0)) then
tmp = t_0
else if (m_1 <= 13500.0d0) then
tmp = exp((abs((n - m)) + (((-0.25d0) * ((n + m) * (n + m))) - l)))
else
tmp = t_0
end if
code = tmp
end function
public static double code(double K, double m, double n, double M, double l) {
double t_0 = Math.cos(M) * Math.exp((0.0 - (M * M)));
double tmp;
if (M <= -125000.0) {
tmp = t_0;
} else if (M <= 13500.0) {
tmp = Math.exp((Math.abs((n - m)) + ((-0.25 * ((n + m) * (n + m))) - l)));
} else {
tmp = t_0;
}
return tmp;
}
def code(K, m, n, M, l): t_0 = math.cos(M) * math.exp((0.0 - (M * M))) tmp = 0 if M <= -125000.0: tmp = t_0 elif M <= 13500.0: tmp = math.exp((math.fabs((n - m)) + ((-0.25 * ((n + m) * (n + m))) - l))) else: tmp = t_0 return tmp
function code(K, m, n, M, l) t_0 = Float64(cos(M) * exp(Float64(0.0 - Float64(M * M)))) tmp = 0.0 if (M <= -125000.0) tmp = t_0; elseif (M <= 13500.0) tmp = exp(Float64(abs(Float64(n - m)) + Float64(Float64(-0.25 * Float64(Float64(n + m) * Float64(n + m))) - l))); else tmp = t_0; end return tmp end
function tmp_2 = code(K, m, n, M, l) t_0 = cos(M) * exp((0.0 - (M * M))); tmp = 0.0; if (M <= -125000.0) tmp = t_0; elseif (M <= 13500.0) tmp = exp((abs((n - m)) + ((-0.25 * ((n + m) * (n + m))) - l))); else tmp = t_0; end tmp_2 = tmp; end
code[K_, m_, n_, M_, l_] := Block[{t$95$0 = N[(N[Cos[M], $MachinePrecision] * N[Exp[N[(0.0 - N[(M * M), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[M, -125000.0], t$95$0, If[LessEqual[M, 13500.0], N[Exp[N[(N[Abs[N[(n - m), $MachinePrecision]], $MachinePrecision] + N[(N[(-0.25 * N[(N[(n + m), $MachinePrecision] * N[(n + m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - l), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], t$95$0]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \cos M \cdot e^{0 - M \cdot M}\\
\mathbf{if}\;M \leq -125000:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;M \leq 13500:\\
\;\;\;\;e^{\left|n - m\right| + \left(-0.25 \cdot \left(\left(n + m\right) \cdot \left(n + m\right)\right) - \ell\right)}\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if M < -125000 or 13500 < M Initial program 77.0%
*-lowering-*.f64N/A
sub-negN/A
+-commutativeN/A
neg-sub0N/A
associate-+l-N/A
sub0-negN/A
cos-negN/A
cos-lowering-cos.f64N/A
sub-negN/A
+-lowering-+.f64N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
metadata-evalN/A
exp-lowering-exp.f64N/A
Simplified77.0%
Taylor expanded in K around 0
*-commutativeN/A
*-lowering-*.f64N/A
Simplified99.2%
Taylor expanded in M around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
unpow2N/A
*-lowering-*.f6498.4%
Simplified98.4%
if -125000 < M < 13500Initial program 79.4%
*-lowering-*.f64N/A
sub-negN/A
+-commutativeN/A
neg-sub0N/A
associate-+l-N/A
sub0-negN/A
cos-negN/A
cos-lowering-cos.f64N/A
sub-negN/A
+-lowering-+.f64N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
metadata-evalN/A
exp-lowering-exp.f64N/A
Simplified79.4%
Taylor expanded in K around 0
*-commutativeN/A
*-lowering-*.f64N/A
Simplified98.6%
Taylor expanded in M around 0
exp-lowering-exp.f64N/A
associate--l+N/A
+-lowering-+.f64N/A
fabs-subN/A
sub-negN/A
mul-1-negN/A
fabs-lowering-fabs.f64N/A
mul-1-negN/A
sub-negN/A
--lowering--.f64N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
+-commutativeN/A
+-lowering-+.f6498.1%
Simplified98.1%
Final simplification98.3%
(FPCore (K m n M l)
:precision binary64
(let* ((t_0 (* -0.25 (* n n))))
(if (<= n -3.6e-158)
(exp (+ (fabs (- n m)) (* -0.25 (* m m))))
(if (<= n 1.45e-178)
(* (cos M) t_0)
(if (<= n 54.0) (* (cos M) (exp (- 0.0 (* M M)))) (exp t_0))))))
double code(double K, double m, double n, double M, double l) {
double t_0 = -0.25 * (n * n);
double tmp;
if (n <= -3.6e-158) {
tmp = exp((fabs((n - m)) + (-0.25 * (m * m))));
} else if (n <= 1.45e-178) {
tmp = cos(M) * t_0;
} else if (n <= 54.0) {
tmp = cos(M) * exp((0.0 - (M * M)));
} else {
tmp = exp(t_0);
}
return tmp;
}
real(8) function code(k, m, n, m_1, l)
real(8), intent (in) :: k
real(8), intent (in) :: m
real(8), intent (in) :: n
real(8), intent (in) :: m_1
real(8), intent (in) :: l
real(8) :: t_0
real(8) :: tmp
t_0 = (-0.25d0) * (n * n)
if (n <= (-3.6d-158)) then
tmp = exp((abs((n - m)) + ((-0.25d0) * (m * m))))
else if (n <= 1.45d-178) then
tmp = cos(m_1) * t_0
else if (n <= 54.0d0) then
tmp = cos(m_1) * exp((0.0d0 - (m_1 * m_1)))
else
tmp = exp(t_0)
end if
code = tmp
end function
public static double code(double K, double m, double n, double M, double l) {
double t_0 = -0.25 * (n * n);
double tmp;
if (n <= -3.6e-158) {
tmp = Math.exp((Math.abs((n - m)) + (-0.25 * (m * m))));
} else if (n <= 1.45e-178) {
tmp = Math.cos(M) * t_0;
} else if (n <= 54.0) {
tmp = Math.cos(M) * Math.exp((0.0 - (M * M)));
} else {
tmp = Math.exp(t_0);
}
return tmp;
}
def code(K, m, n, M, l): t_0 = -0.25 * (n * n) tmp = 0 if n <= -3.6e-158: tmp = math.exp((math.fabs((n - m)) + (-0.25 * (m * m)))) elif n <= 1.45e-178: tmp = math.cos(M) * t_0 elif n <= 54.0: tmp = math.cos(M) * math.exp((0.0 - (M * M))) else: tmp = math.exp(t_0) return tmp
function code(K, m, n, M, l) t_0 = Float64(-0.25 * Float64(n * n)) tmp = 0.0 if (n <= -3.6e-158) tmp = exp(Float64(abs(Float64(n - m)) + Float64(-0.25 * Float64(m * m)))); elseif (n <= 1.45e-178) tmp = Float64(cos(M) * t_0); elseif (n <= 54.0) tmp = Float64(cos(M) * exp(Float64(0.0 - Float64(M * M)))); else tmp = exp(t_0); end return tmp end
function tmp_2 = code(K, m, n, M, l) t_0 = -0.25 * (n * n); tmp = 0.0; if (n <= -3.6e-158) tmp = exp((abs((n - m)) + (-0.25 * (m * m)))); elseif (n <= 1.45e-178) tmp = cos(M) * t_0; elseif (n <= 54.0) tmp = cos(M) * exp((0.0 - (M * M))); else tmp = exp(t_0); end tmp_2 = tmp; end
code[K_, m_, n_, M_, l_] := Block[{t$95$0 = N[(-0.25 * N[(n * n), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[n, -3.6e-158], N[Exp[N[(N[Abs[N[(n - m), $MachinePrecision]], $MachinePrecision] + N[(-0.25 * N[(m * m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], If[LessEqual[n, 1.45e-178], N[(N[Cos[M], $MachinePrecision] * t$95$0), $MachinePrecision], If[LessEqual[n, 54.0], N[(N[Cos[M], $MachinePrecision] * N[Exp[N[(0.0 - N[(M * M), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], N[Exp[t$95$0], $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := -0.25 \cdot \left(n \cdot n\right)\\
\mathbf{if}\;n \leq -3.6 \cdot 10^{-158}:\\
\;\;\;\;e^{\left|n - m\right| + -0.25 \cdot \left(m \cdot m\right)}\\
\mathbf{elif}\;n \leq 1.45 \cdot 10^{-178}:\\
\;\;\;\;\cos M \cdot t\_0\\
\mathbf{elif}\;n \leq 54:\\
\;\;\;\;\cos M \cdot e^{0 - M \cdot M}\\
\mathbf{else}:\\
\;\;\;\;e^{t\_0}\\
\end{array}
\end{array}
if n < -3.59999999999999991e-158Initial program 74.4%
*-lowering-*.f64N/A
sub-negN/A
+-commutativeN/A
neg-sub0N/A
associate-+l-N/A
sub0-negN/A
cos-negN/A
cos-lowering-cos.f64N/A
sub-negN/A
+-lowering-+.f64N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
metadata-evalN/A
exp-lowering-exp.f64N/A
Simplified74.4%
Taylor expanded in K around 0
*-commutativeN/A
*-lowering-*.f64N/A
Simplified97.6%
Taylor expanded in M around 0
exp-lowering-exp.f64N/A
associate--l+N/A
+-lowering-+.f64N/A
fabs-subN/A
sub-negN/A
mul-1-negN/A
fabs-lowering-fabs.f64N/A
mul-1-negN/A
sub-negN/A
--lowering--.f64N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
+-commutativeN/A
+-lowering-+.f6490.4%
Simplified90.4%
Taylor expanded in m around inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6442.5%
Simplified42.5%
if -3.59999999999999991e-158 < n < 1.4499999999999999e-178Initial program 82.3%
*-lowering-*.f64N/A
sub-negN/A
+-commutativeN/A
neg-sub0N/A
associate-+l-N/A
sub0-negN/A
cos-negN/A
cos-lowering-cos.f64N/A
sub-negN/A
+-lowering-+.f64N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
metadata-evalN/A
exp-lowering-exp.f64N/A
Simplified82.3%
Taylor expanded in K around 0
*-commutativeN/A
*-lowering-*.f64N/A
Simplified100.0%
Taylor expanded in n around inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f648.3%
Simplified8.3%
Taylor expanded in n around 0
associate-*r*N/A
distribute-rgt1-inN/A
+-commutativeN/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
cos-lowering-cos.f648.3%
Simplified8.3%
Taylor expanded in n around inf
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
cos-lowering-cos.f6480.3%
Simplified80.3%
if 1.4499999999999999e-178 < n < 54Initial program 89.5%
*-lowering-*.f64N/A
sub-negN/A
+-commutativeN/A
neg-sub0N/A
associate-+l-N/A
sub0-negN/A
cos-negN/A
cos-lowering-cos.f64N/A
sub-negN/A
+-lowering-+.f64N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
metadata-evalN/A
exp-lowering-exp.f64N/A
Simplified89.5%
Taylor expanded in K around 0
*-commutativeN/A
*-lowering-*.f64N/A
Simplified97.7%
Taylor expanded in M around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
unpow2N/A
*-lowering-*.f6457.4%
Simplified57.4%
if 54 < n Initial program 73.7%
*-lowering-*.f64N/A
sub-negN/A
+-commutativeN/A
neg-sub0N/A
associate-+l-N/A
sub0-negN/A
cos-negN/A
cos-lowering-cos.f64N/A
sub-negN/A
+-lowering-+.f64N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
metadata-evalN/A
exp-lowering-exp.f64N/A
Simplified73.7%
Taylor expanded in K around 0
*-commutativeN/A
*-lowering-*.f64N/A
Simplified100.0%
Taylor expanded in M around 0
exp-lowering-exp.f64N/A
associate--l+N/A
+-lowering-+.f64N/A
fabs-subN/A
sub-negN/A
mul-1-negN/A
fabs-lowering-fabs.f64N/A
mul-1-negN/A
sub-negN/A
--lowering--.f64N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
+-commutativeN/A
+-lowering-+.f6498.7%
Simplified98.7%
Taylor expanded in n around inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6498.7%
Simplified98.7%
Final simplification70.4%
(FPCore (K m n M l)
:precision binary64
(let* ((t_0 (* -0.25 (* n n))))
(if (<= n -1.02e-157)
(* (cos M) (exp (* -0.25 (* m m))))
(if (<= n 1.4e-178)
(* (cos M) t_0)
(if (<= n 54.0) (* (cos M) (exp (- 0.0 (* M M)))) (exp t_0))))))
double code(double K, double m, double n, double M, double l) {
double t_0 = -0.25 * (n * n);
double tmp;
if (n <= -1.02e-157) {
tmp = cos(M) * exp((-0.25 * (m * m)));
} else if (n <= 1.4e-178) {
tmp = cos(M) * t_0;
} else if (n <= 54.0) {
tmp = cos(M) * exp((0.0 - (M * M)));
} else {
tmp = exp(t_0);
}
return tmp;
}
real(8) function code(k, m, n, m_1, l)
real(8), intent (in) :: k
real(8), intent (in) :: m
real(8), intent (in) :: n
real(8), intent (in) :: m_1
real(8), intent (in) :: l
real(8) :: t_0
real(8) :: tmp
t_0 = (-0.25d0) * (n * n)
if (n <= (-1.02d-157)) then
tmp = cos(m_1) * exp(((-0.25d0) * (m * m)))
else if (n <= 1.4d-178) then
tmp = cos(m_1) * t_0
else if (n <= 54.0d0) then
tmp = cos(m_1) * exp((0.0d0 - (m_1 * m_1)))
else
tmp = exp(t_0)
end if
code = tmp
end function
public static double code(double K, double m, double n, double M, double l) {
double t_0 = -0.25 * (n * n);
double tmp;
if (n <= -1.02e-157) {
tmp = Math.cos(M) * Math.exp((-0.25 * (m * m)));
} else if (n <= 1.4e-178) {
tmp = Math.cos(M) * t_0;
} else if (n <= 54.0) {
tmp = Math.cos(M) * Math.exp((0.0 - (M * M)));
} else {
tmp = Math.exp(t_0);
}
return tmp;
}
def code(K, m, n, M, l): t_0 = -0.25 * (n * n) tmp = 0 if n <= -1.02e-157: tmp = math.cos(M) * math.exp((-0.25 * (m * m))) elif n <= 1.4e-178: tmp = math.cos(M) * t_0 elif n <= 54.0: tmp = math.cos(M) * math.exp((0.0 - (M * M))) else: tmp = math.exp(t_0) return tmp
function code(K, m, n, M, l) t_0 = Float64(-0.25 * Float64(n * n)) tmp = 0.0 if (n <= -1.02e-157) tmp = Float64(cos(M) * exp(Float64(-0.25 * Float64(m * m)))); elseif (n <= 1.4e-178) tmp = Float64(cos(M) * t_0); elseif (n <= 54.0) tmp = Float64(cos(M) * exp(Float64(0.0 - Float64(M * M)))); else tmp = exp(t_0); end return tmp end
function tmp_2 = code(K, m, n, M, l) t_0 = -0.25 * (n * n); tmp = 0.0; if (n <= -1.02e-157) tmp = cos(M) * exp((-0.25 * (m * m))); elseif (n <= 1.4e-178) tmp = cos(M) * t_0; elseif (n <= 54.0) tmp = cos(M) * exp((0.0 - (M * M))); else tmp = exp(t_0); end tmp_2 = tmp; end
code[K_, m_, n_, M_, l_] := Block[{t$95$0 = N[(-0.25 * N[(n * n), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[n, -1.02e-157], N[(N[Cos[M], $MachinePrecision] * N[Exp[N[(-0.25 * N[(m * m), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], If[LessEqual[n, 1.4e-178], N[(N[Cos[M], $MachinePrecision] * t$95$0), $MachinePrecision], If[LessEqual[n, 54.0], N[(N[Cos[M], $MachinePrecision] * N[Exp[N[(0.0 - N[(M * M), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], N[Exp[t$95$0], $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := -0.25 \cdot \left(n \cdot n\right)\\
\mathbf{if}\;n \leq -1.02 \cdot 10^{-157}:\\
\;\;\;\;\cos M \cdot e^{-0.25 \cdot \left(m \cdot m\right)}\\
\mathbf{elif}\;n \leq 1.4 \cdot 10^{-178}:\\
\;\;\;\;\cos M \cdot t\_0\\
\mathbf{elif}\;n \leq 54:\\
\;\;\;\;\cos M \cdot e^{0 - M \cdot M}\\
\mathbf{else}:\\
\;\;\;\;e^{t\_0}\\
\end{array}
\end{array}
if n < -1.0200000000000001e-157Initial program 74.4%
*-lowering-*.f64N/A
sub-negN/A
+-commutativeN/A
neg-sub0N/A
associate-+l-N/A
sub0-negN/A
cos-negN/A
cos-lowering-cos.f64N/A
sub-negN/A
+-lowering-+.f64N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
metadata-evalN/A
exp-lowering-exp.f64N/A
Simplified74.4%
Taylor expanded in K around 0
*-commutativeN/A
*-lowering-*.f64N/A
Simplified97.6%
Taylor expanded in m around inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6443.7%
Simplified43.7%
if -1.0200000000000001e-157 < n < 1.4000000000000001e-178Initial program 82.3%
*-lowering-*.f64N/A
sub-negN/A
+-commutativeN/A
neg-sub0N/A
associate-+l-N/A
sub0-negN/A
cos-negN/A
cos-lowering-cos.f64N/A
sub-negN/A
+-lowering-+.f64N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
metadata-evalN/A
exp-lowering-exp.f64N/A
Simplified82.3%
Taylor expanded in K around 0
*-commutativeN/A
*-lowering-*.f64N/A
Simplified100.0%
Taylor expanded in n around inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f648.3%
Simplified8.3%
Taylor expanded in n around 0
associate-*r*N/A
distribute-rgt1-inN/A
+-commutativeN/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
cos-lowering-cos.f648.3%
Simplified8.3%
Taylor expanded in n around inf
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
cos-lowering-cos.f6480.3%
Simplified80.3%
if 1.4000000000000001e-178 < n < 54Initial program 89.5%
*-lowering-*.f64N/A
sub-negN/A
+-commutativeN/A
neg-sub0N/A
associate-+l-N/A
sub0-negN/A
cos-negN/A
cos-lowering-cos.f64N/A
sub-negN/A
+-lowering-+.f64N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
metadata-evalN/A
exp-lowering-exp.f64N/A
Simplified89.5%
Taylor expanded in K around 0
*-commutativeN/A
*-lowering-*.f64N/A
Simplified97.7%
Taylor expanded in M around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
unpow2N/A
*-lowering-*.f6457.4%
Simplified57.4%
if 54 < n Initial program 73.7%
*-lowering-*.f64N/A
sub-negN/A
+-commutativeN/A
neg-sub0N/A
associate-+l-N/A
sub0-negN/A
cos-negN/A
cos-lowering-cos.f64N/A
sub-negN/A
+-lowering-+.f64N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
metadata-evalN/A
exp-lowering-exp.f64N/A
Simplified73.7%
Taylor expanded in K around 0
*-commutativeN/A
*-lowering-*.f64N/A
Simplified100.0%
Taylor expanded in M around 0
exp-lowering-exp.f64N/A
associate--l+N/A
+-lowering-+.f64N/A
fabs-subN/A
sub-negN/A
mul-1-negN/A
fabs-lowering-fabs.f64N/A
mul-1-negN/A
sub-negN/A
--lowering--.f64N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
+-commutativeN/A
+-lowering-+.f6498.7%
Simplified98.7%
Taylor expanded in n around inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6498.7%
Simplified98.7%
Final simplification70.8%
(FPCore (K m n M l)
:precision binary64
(let* ((t_0 (* -0.25 (* n n))) (t_1 (* (cos M) (exp (* -0.25 (* m m))))))
(if (<= n -3.3e-158)
t_1
(if (<= n 2.25e-147) (* (cos M) t_0) (if (<= n 54.0) t_1 (exp t_0))))))
double code(double K, double m, double n, double M, double l) {
double t_0 = -0.25 * (n * n);
double t_1 = cos(M) * exp((-0.25 * (m * m)));
double tmp;
if (n <= -3.3e-158) {
tmp = t_1;
} else if (n <= 2.25e-147) {
tmp = cos(M) * t_0;
} else if (n <= 54.0) {
tmp = t_1;
} else {
tmp = exp(t_0);
}
return tmp;
}
real(8) function code(k, m, n, m_1, l)
real(8), intent (in) :: k
real(8), intent (in) :: m
real(8), intent (in) :: n
real(8), intent (in) :: m_1
real(8), intent (in) :: l
real(8) :: t_0
real(8) :: t_1
real(8) :: tmp
t_0 = (-0.25d0) * (n * n)
t_1 = cos(m_1) * exp(((-0.25d0) * (m * m)))
if (n <= (-3.3d-158)) then
tmp = t_1
else if (n <= 2.25d-147) then
tmp = cos(m_1) * t_0
else if (n <= 54.0d0) then
tmp = t_1
else
tmp = exp(t_0)
end if
code = tmp
end function
public static double code(double K, double m, double n, double M, double l) {
double t_0 = -0.25 * (n * n);
double t_1 = Math.cos(M) * Math.exp((-0.25 * (m * m)));
double tmp;
if (n <= -3.3e-158) {
tmp = t_1;
} else if (n <= 2.25e-147) {
tmp = Math.cos(M) * t_0;
} else if (n <= 54.0) {
tmp = t_1;
} else {
tmp = Math.exp(t_0);
}
return tmp;
}
def code(K, m, n, M, l): t_0 = -0.25 * (n * n) t_1 = math.cos(M) * math.exp((-0.25 * (m * m))) tmp = 0 if n <= -3.3e-158: tmp = t_1 elif n <= 2.25e-147: tmp = math.cos(M) * t_0 elif n <= 54.0: tmp = t_1 else: tmp = math.exp(t_0) return tmp
function code(K, m, n, M, l) t_0 = Float64(-0.25 * Float64(n * n)) t_1 = Float64(cos(M) * exp(Float64(-0.25 * Float64(m * m)))) tmp = 0.0 if (n <= -3.3e-158) tmp = t_1; elseif (n <= 2.25e-147) tmp = Float64(cos(M) * t_0); elseif (n <= 54.0) tmp = t_1; else tmp = exp(t_0); end return tmp end
function tmp_2 = code(K, m, n, M, l) t_0 = -0.25 * (n * n); t_1 = cos(M) * exp((-0.25 * (m * m))); tmp = 0.0; if (n <= -3.3e-158) tmp = t_1; elseif (n <= 2.25e-147) tmp = cos(M) * t_0; elseif (n <= 54.0) tmp = t_1; else tmp = exp(t_0); end tmp_2 = tmp; end
code[K_, m_, n_, M_, l_] := Block[{t$95$0 = N[(-0.25 * N[(n * n), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[Cos[M], $MachinePrecision] * N[Exp[N[(-0.25 * N[(m * m), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[n, -3.3e-158], t$95$1, If[LessEqual[n, 2.25e-147], N[(N[Cos[M], $MachinePrecision] * t$95$0), $MachinePrecision], If[LessEqual[n, 54.0], t$95$1, N[Exp[t$95$0], $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := -0.25 \cdot \left(n \cdot n\right)\\
t_1 := \cos M \cdot e^{-0.25 \cdot \left(m \cdot m\right)}\\
\mathbf{if}\;n \leq -3.3 \cdot 10^{-158}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;n \leq 2.25 \cdot 10^{-147}:\\
\;\;\;\;\cos M \cdot t\_0\\
\mathbf{elif}\;n \leq 54:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;e^{t\_0}\\
\end{array}
\end{array}
if n < -3.3000000000000002e-158 or 2.24999999999999986e-147 < n < 54Initial program 78.8%
*-lowering-*.f64N/A
sub-negN/A
+-commutativeN/A
neg-sub0N/A
associate-+l-N/A
sub0-negN/A
cos-negN/A
cos-lowering-cos.f64N/A
sub-negN/A
+-lowering-+.f64N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
metadata-evalN/A
exp-lowering-exp.f64N/A
Simplified78.8%
Taylor expanded in K around 0
*-commutativeN/A
*-lowering-*.f64N/A
Simplified98.2%
Taylor expanded in m around inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6446.3%
Simplified46.3%
if -3.3000000000000002e-158 < n < 2.24999999999999986e-147Initial program 82.4%
*-lowering-*.f64N/A
sub-negN/A
+-commutativeN/A
neg-sub0N/A
associate-+l-N/A
sub0-negN/A
cos-negN/A
cos-lowering-cos.f64N/A
sub-negN/A
+-lowering-+.f64N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
metadata-evalN/A
exp-lowering-exp.f64N/A
Simplified82.4%
Taylor expanded in K around 0
*-commutativeN/A
*-lowering-*.f64N/A
Simplified98.8%
Taylor expanded in n around inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f648.1%
Simplified8.1%
Taylor expanded in n around 0
associate-*r*N/A
distribute-rgt1-inN/A
+-commutativeN/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
cos-lowering-cos.f648.1%
Simplified8.1%
Taylor expanded in n around inf
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
cos-lowering-cos.f6478.6%
Simplified78.6%
if 54 < n Initial program 73.7%
*-lowering-*.f64N/A
sub-negN/A
+-commutativeN/A
neg-sub0N/A
associate-+l-N/A
sub0-negN/A
cos-negN/A
cos-lowering-cos.f64N/A
sub-negN/A
+-lowering-+.f64N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
metadata-evalN/A
exp-lowering-exp.f64N/A
Simplified73.7%
Taylor expanded in K around 0
*-commutativeN/A
*-lowering-*.f64N/A
Simplified100.0%
Taylor expanded in M around 0
exp-lowering-exp.f64N/A
associate--l+N/A
+-lowering-+.f64N/A
fabs-subN/A
sub-negN/A
mul-1-negN/A
fabs-lowering-fabs.f64N/A
mul-1-negN/A
sub-negN/A
--lowering--.f64N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
+-commutativeN/A
+-lowering-+.f6498.7%
Simplified98.7%
Taylor expanded in n around inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6498.7%
Simplified98.7%
Final simplification70.3%
(FPCore (K m n M l)
:precision binary64
(let* ((t_0 (* -0.25 (* n n))) (t_1 (exp (* -0.25 (* m m)))))
(if (<= n -5.6e-160)
t_1
(if (<= n 2.25e-147) (* (cos M) t_0) (if (<= n 54.0) t_1 (exp t_0))))))
double code(double K, double m, double n, double M, double l) {
double t_0 = -0.25 * (n * n);
double t_1 = exp((-0.25 * (m * m)));
double tmp;
if (n <= -5.6e-160) {
tmp = t_1;
} else if (n <= 2.25e-147) {
tmp = cos(M) * t_0;
} else if (n <= 54.0) {
tmp = t_1;
} else {
tmp = exp(t_0);
}
return tmp;
}
real(8) function code(k, m, n, m_1, l)
real(8), intent (in) :: k
real(8), intent (in) :: m
real(8), intent (in) :: n
real(8), intent (in) :: m_1
real(8), intent (in) :: l
real(8) :: t_0
real(8) :: t_1
real(8) :: tmp
t_0 = (-0.25d0) * (n * n)
t_1 = exp(((-0.25d0) * (m * m)))
if (n <= (-5.6d-160)) then
tmp = t_1
else if (n <= 2.25d-147) then
tmp = cos(m_1) * t_0
else if (n <= 54.0d0) then
tmp = t_1
else
tmp = exp(t_0)
end if
code = tmp
end function
public static double code(double K, double m, double n, double M, double l) {
double t_0 = -0.25 * (n * n);
double t_1 = Math.exp((-0.25 * (m * m)));
double tmp;
if (n <= -5.6e-160) {
tmp = t_1;
} else if (n <= 2.25e-147) {
tmp = Math.cos(M) * t_0;
} else if (n <= 54.0) {
tmp = t_1;
} else {
tmp = Math.exp(t_0);
}
return tmp;
}
def code(K, m, n, M, l): t_0 = -0.25 * (n * n) t_1 = math.exp((-0.25 * (m * m))) tmp = 0 if n <= -5.6e-160: tmp = t_1 elif n <= 2.25e-147: tmp = math.cos(M) * t_0 elif n <= 54.0: tmp = t_1 else: tmp = math.exp(t_0) return tmp
function code(K, m, n, M, l) t_0 = Float64(-0.25 * Float64(n * n)) t_1 = exp(Float64(-0.25 * Float64(m * m))) tmp = 0.0 if (n <= -5.6e-160) tmp = t_1; elseif (n <= 2.25e-147) tmp = Float64(cos(M) * t_0); elseif (n <= 54.0) tmp = t_1; else tmp = exp(t_0); end return tmp end
function tmp_2 = code(K, m, n, M, l) t_0 = -0.25 * (n * n); t_1 = exp((-0.25 * (m * m))); tmp = 0.0; if (n <= -5.6e-160) tmp = t_1; elseif (n <= 2.25e-147) tmp = cos(M) * t_0; elseif (n <= 54.0) tmp = t_1; else tmp = exp(t_0); end tmp_2 = tmp; end
code[K_, m_, n_, M_, l_] := Block[{t$95$0 = N[(-0.25 * N[(n * n), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[Exp[N[(-0.25 * N[(m * m), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]}, If[LessEqual[n, -5.6e-160], t$95$1, If[LessEqual[n, 2.25e-147], N[(N[Cos[M], $MachinePrecision] * t$95$0), $MachinePrecision], If[LessEqual[n, 54.0], t$95$1, N[Exp[t$95$0], $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := -0.25 \cdot \left(n \cdot n\right)\\
t_1 := e^{-0.25 \cdot \left(m \cdot m\right)}\\
\mathbf{if}\;n \leq -5.6 \cdot 10^{-160}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;n \leq 2.25 \cdot 10^{-147}:\\
\;\;\;\;\cos M \cdot t\_0\\
\mathbf{elif}\;n \leq 54:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;e^{t\_0}\\
\end{array}
\end{array}
if n < -5.60000000000000032e-160 or 2.24999999999999986e-147 < n < 54Initial program 79.0%
*-lowering-*.f64N/A
sub-negN/A
+-commutativeN/A
neg-sub0N/A
associate-+l-N/A
sub0-negN/A
cos-negN/A
cos-lowering-cos.f64N/A
sub-negN/A
+-lowering-+.f64N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
metadata-evalN/A
exp-lowering-exp.f64N/A
Simplified79.0%
Taylor expanded in K around 0
*-commutativeN/A
*-lowering-*.f64N/A
Simplified98.3%
Taylor expanded in M around 0
exp-lowering-exp.f64N/A
associate--l+N/A
+-lowering-+.f64N/A
fabs-subN/A
sub-negN/A
mul-1-negN/A
fabs-lowering-fabs.f64N/A
mul-1-negN/A
sub-negN/A
--lowering--.f64N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
+-commutativeN/A
+-lowering-+.f6484.5%
Simplified84.5%
Taylor expanded in m around inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6445.9%
Simplified45.9%
if -5.60000000000000032e-160 < n < 2.24999999999999986e-147Initial program 82.1%
*-lowering-*.f64N/A
sub-negN/A
+-commutativeN/A
neg-sub0N/A
associate-+l-N/A
sub0-negN/A
cos-negN/A
cos-lowering-cos.f64N/A
sub-negN/A
+-lowering-+.f64N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
metadata-evalN/A
exp-lowering-exp.f64N/A
Simplified82.1%
Taylor expanded in K around 0
*-commutativeN/A
*-lowering-*.f64N/A
Simplified98.8%
Taylor expanded in n around inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f648.2%
Simplified8.2%
Taylor expanded in n around 0
associate-*r*N/A
distribute-rgt1-inN/A
+-commutativeN/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
cos-lowering-cos.f648.2%
Simplified8.2%
Taylor expanded in n around inf
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
cos-lowering-cos.f6478.7%
Simplified78.7%
if 54 < n Initial program 73.7%
*-lowering-*.f64N/A
sub-negN/A
+-commutativeN/A
neg-sub0N/A
associate-+l-N/A
sub0-negN/A
cos-negN/A
cos-lowering-cos.f64N/A
sub-negN/A
+-lowering-+.f64N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
metadata-evalN/A
exp-lowering-exp.f64N/A
Simplified73.7%
Taylor expanded in K around 0
*-commutativeN/A
*-lowering-*.f64N/A
Simplified100.0%
Taylor expanded in M around 0
exp-lowering-exp.f64N/A
associate--l+N/A
+-lowering-+.f64N/A
fabs-subN/A
sub-negN/A
mul-1-negN/A
fabs-lowering-fabs.f64N/A
mul-1-negN/A
sub-negN/A
--lowering--.f64N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
+-commutativeN/A
+-lowering-+.f6498.7%
Simplified98.7%
Taylor expanded in n around inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6498.7%
Simplified98.7%
Final simplification70.0%
(FPCore (K m n M l) :precision binary64 (let* ((t_0 (exp (* -0.25 (* m m))))) (if (<= m -3.55) t_0 (if (<= m 3.3e-78) (exp (- 0.0 l)) t_0))))
double code(double K, double m, double n, double M, double l) {
double t_0 = exp((-0.25 * (m * m)));
double tmp;
if (m <= -3.55) {
tmp = t_0;
} else if (m <= 3.3e-78) {
tmp = exp((0.0 - l));
} else {
tmp = t_0;
}
return tmp;
}
real(8) function code(k, m, n, m_1, l)
real(8), intent (in) :: k
real(8), intent (in) :: m
real(8), intent (in) :: n
real(8), intent (in) :: m_1
real(8), intent (in) :: l
real(8) :: t_0
real(8) :: tmp
t_0 = exp(((-0.25d0) * (m * m)))
if (m <= (-3.55d0)) then
tmp = t_0
else if (m <= 3.3d-78) then
tmp = exp((0.0d0 - l))
else
tmp = t_0
end if
code = tmp
end function
public static double code(double K, double m, double n, double M, double l) {
double t_0 = Math.exp((-0.25 * (m * m)));
double tmp;
if (m <= -3.55) {
tmp = t_0;
} else if (m <= 3.3e-78) {
tmp = Math.exp((0.0 - l));
} else {
tmp = t_0;
}
return tmp;
}
def code(K, m, n, M, l): t_0 = math.exp((-0.25 * (m * m))) tmp = 0 if m <= -3.55: tmp = t_0 elif m <= 3.3e-78: tmp = math.exp((0.0 - l)) else: tmp = t_0 return tmp
function code(K, m, n, M, l) t_0 = exp(Float64(-0.25 * Float64(m * m))) tmp = 0.0 if (m <= -3.55) tmp = t_0; elseif (m <= 3.3e-78) tmp = exp(Float64(0.0 - l)); else tmp = t_0; end return tmp end
function tmp_2 = code(K, m, n, M, l) t_0 = exp((-0.25 * (m * m))); tmp = 0.0; if (m <= -3.55) tmp = t_0; elseif (m <= 3.3e-78) tmp = exp((0.0 - l)); else tmp = t_0; end tmp_2 = tmp; end
code[K_, m_, n_, M_, l_] := Block[{t$95$0 = N[Exp[N[(-0.25 * N[(m * m), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]}, If[LessEqual[m, -3.55], t$95$0, If[LessEqual[m, 3.3e-78], N[Exp[N[(0.0 - l), $MachinePrecision]], $MachinePrecision], t$95$0]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := e^{-0.25 \cdot \left(m \cdot m\right)}\\
\mathbf{if}\;m \leq -3.55:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;m \leq 3.3 \cdot 10^{-78}:\\
\;\;\;\;e^{0 - \ell}\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if m < -3.5499999999999998 or 3.29999999999999982e-78 < m Initial program 73.9%
*-lowering-*.f64N/A
sub-negN/A
+-commutativeN/A
neg-sub0N/A
associate-+l-N/A
sub0-negN/A
cos-negN/A
cos-lowering-cos.f64N/A
sub-negN/A
+-lowering-+.f64N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
metadata-evalN/A
exp-lowering-exp.f64N/A
Simplified73.9%
Taylor expanded in K around 0
*-commutativeN/A
*-lowering-*.f64N/A
Simplified99.3%
Taylor expanded in M around 0
exp-lowering-exp.f64N/A
associate--l+N/A
+-lowering-+.f64N/A
fabs-subN/A
sub-negN/A
mul-1-negN/A
fabs-lowering-fabs.f64N/A
mul-1-negN/A
sub-negN/A
--lowering--.f64N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
+-commutativeN/A
+-lowering-+.f6495.6%
Simplified95.6%
Taylor expanded in m around inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6490.6%
Simplified90.6%
if -3.5499999999999998 < m < 3.29999999999999982e-78Initial program 83.0%
*-lowering-*.f64N/A
sub-negN/A
+-commutativeN/A
neg-sub0N/A
associate-+l-N/A
sub0-negN/A
cos-negN/A
cos-lowering-cos.f64N/A
sub-negN/A
+-lowering-+.f64N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
metadata-evalN/A
exp-lowering-exp.f64N/A
Simplified83.0%
Taylor expanded in K around 0
*-commutativeN/A
*-lowering-*.f64N/A
Simplified98.5%
Taylor expanded in M around 0
exp-lowering-exp.f64N/A
associate--l+N/A
+-lowering-+.f64N/A
fabs-subN/A
sub-negN/A
mul-1-negN/A
fabs-lowering-fabs.f64N/A
mul-1-negN/A
sub-negN/A
--lowering--.f64N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
+-commutativeN/A
+-lowering-+.f6480.3%
Simplified80.3%
Taylor expanded in l around inf
neg-mul-1N/A
neg-sub0N/A
--lowering--.f6438.3%
Simplified38.3%
(FPCore (K m n M l) :precision binary64 (if (<= m -0.0106) (exp (* -0.25 (* m m))) (exp (* -0.25 (* n n)))))
double code(double K, double m, double n, double M, double l) {
double tmp;
if (m <= -0.0106) {
tmp = exp((-0.25 * (m * m)));
} else {
tmp = exp((-0.25 * (n * n)));
}
return tmp;
}
real(8) function code(k, m, n, m_1, l)
real(8), intent (in) :: k
real(8), intent (in) :: m
real(8), intent (in) :: n
real(8), intent (in) :: m_1
real(8), intent (in) :: l
real(8) :: tmp
if (m <= (-0.0106d0)) then
tmp = exp(((-0.25d0) * (m * m)))
else
tmp = exp(((-0.25d0) * (n * n)))
end if
code = tmp
end function
public static double code(double K, double m, double n, double M, double l) {
double tmp;
if (m <= -0.0106) {
tmp = Math.exp((-0.25 * (m * m)));
} else {
tmp = Math.exp((-0.25 * (n * n)));
}
return tmp;
}
def code(K, m, n, M, l): tmp = 0 if m <= -0.0106: tmp = math.exp((-0.25 * (m * m))) else: tmp = math.exp((-0.25 * (n * n))) return tmp
function code(K, m, n, M, l) tmp = 0.0 if (m <= -0.0106) tmp = exp(Float64(-0.25 * Float64(m * m))); else tmp = exp(Float64(-0.25 * Float64(n * n))); end return tmp end
function tmp_2 = code(K, m, n, M, l) tmp = 0.0; if (m <= -0.0106) tmp = exp((-0.25 * (m * m))); else tmp = exp((-0.25 * (n * n))); end tmp_2 = tmp; end
code[K_, m_, n_, M_, l_] := If[LessEqual[m, -0.0106], N[Exp[N[(-0.25 * N[(m * m), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], N[Exp[N[(-0.25 * N[(n * n), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;m \leq -0.0106:\\
\;\;\;\;e^{-0.25 \cdot \left(m \cdot m\right)}\\
\mathbf{else}:\\
\;\;\;\;e^{-0.25 \cdot \left(n \cdot n\right)}\\
\end{array}
\end{array}
if m < -0.0106Initial program 77.3%
*-lowering-*.f64N/A
sub-negN/A
+-commutativeN/A
neg-sub0N/A
associate-+l-N/A
sub0-negN/A
cos-negN/A
cos-lowering-cos.f64N/A
sub-negN/A
+-lowering-+.f64N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
metadata-evalN/A
exp-lowering-exp.f64N/A
Simplified77.3%
Taylor expanded in K around 0
*-commutativeN/A
*-lowering-*.f64N/A
Simplified98.5%
Taylor expanded in M around 0
exp-lowering-exp.f64N/A
associate--l+N/A
+-lowering-+.f64N/A
fabs-subN/A
sub-negN/A
mul-1-negN/A
fabs-lowering-fabs.f64N/A
mul-1-negN/A
sub-negN/A
--lowering--.f64N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
+-commutativeN/A
+-lowering-+.f6495.5%
Simplified95.5%
Taylor expanded in m around inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6494.1%
Simplified94.1%
if -0.0106 < m Initial program 78.5%
*-lowering-*.f64N/A
sub-negN/A
+-commutativeN/A
neg-sub0N/A
associate-+l-N/A
sub0-negN/A
cos-negN/A
cos-lowering-cos.f64N/A
sub-negN/A
+-lowering-+.f64N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
metadata-evalN/A
exp-lowering-exp.f64N/A
Simplified78.5%
Taylor expanded in K around 0
*-commutativeN/A
*-lowering-*.f64N/A
Simplified99.1%
Taylor expanded in M around 0
exp-lowering-exp.f64N/A
associate--l+N/A
+-lowering-+.f64N/A
fabs-subN/A
sub-negN/A
mul-1-negN/A
fabs-lowering-fabs.f64N/A
mul-1-negN/A
sub-negN/A
--lowering--.f64N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
+-commutativeN/A
+-lowering-+.f6485.8%
Simplified85.8%
Taylor expanded in n around inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6458.1%
Simplified58.1%
(FPCore (K m n M l) :precision binary64 (exp (- 0.0 l)))
double code(double K, double m, double n, double M, double l) {
return exp((0.0 - l));
}
real(8) function code(k, m, n, m_1, l)
real(8), intent (in) :: k
real(8), intent (in) :: m
real(8), intent (in) :: n
real(8), intent (in) :: m_1
real(8), intent (in) :: l
code = exp((0.0d0 - l))
end function
public static double code(double K, double m, double n, double M, double l) {
return Math.exp((0.0 - l));
}
def code(K, m, n, M, l): return math.exp((0.0 - l))
function code(K, m, n, M, l) return exp(Float64(0.0 - l)) end
function tmp = code(K, m, n, M, l) tmp = exp((0.0 - l)); end
code[K_, m_, n_, M_, l_] := N[Exp[N[(0.0 - l), $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
\\
e^{0 - \ell}
\end{array}
Initial program 78.2%
*-lowering-*.f64N/A
sub-negN/A
+-commutativeN/A
neg-sub0N/A
associate-+l-N/A
sub0-negN/A
cos-negN/A
cos-lowering-cos.f64N/A
sub-negN/A
+-lowering-+.f64N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
metadata-evalN/A
exp-lowering-exp.f64N/A
Simplified78.2%
Taylor expanded in K around 0
*-commutativeN/A
*-lowering-*.f64N/A
Simplified98.9%
Taylor expanded in M around 0
exp-lowering-exp.f64N/A
associate--l+N/A
+-lowering-+.f64N/A
fabs-subN/A
sub-negN/A
mul-1-negN/A
fabs-lowering-fabs.f64N/A
mul-1-negN/A
sub-negN/A
--lowering--.f64N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
+-commutativeN/A
+-lowering-+.f6488.3%
Simplified88.3%
Taylor expanded in l around inf
neg-mul-1N/A
neg-sub0N/A
--lowering--.f6435.6%
Simplified35.6%
(FPCore (K m n M l) :precision binary64 (cos M))
double code(double K, double m, double n, double M, double l) {
return cos(M);
}
real(8) function code(k, m, n, m_1, l)
real(8), intent (in) :: k
real(8), intent (in) :: m
real(8), intent (in) :: n
real(8), intent (in) :: m_1
real(8), intent (in) :: l
code = cos(m_1)
end function
public static double code(double K, double m, double n, double M, double l) {
return Math.cos(M);
}
def code(K, m, n, M, l): return math.cos(M)
function code(K, m, n, M, l) return cos(M) end
function tmp = code(K, m, n, M, l) tmp = cos(M); end
code[K_, m_, n_, M_, l_] := N[Cos[M], $MachinePrecision]
\begin{array}{l}
\\
\cos M
\end{array}
Initial program 78.2%
*-lowering-*.f64N/A
sub-negN/A
+-commutativeN/A
neg-sub0N/A
associate-+l-N/A
sub0-negN/A
cos-negN/A
cos-lowering-cos.f64N/A
sub-negN/A
+-lowering-+.f64N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
metadata-evalN/A
exp-lowering-exp.f64N/A
Simplified78.2%
Taylor expanded in K around 0
*-commutativeN/A
*-lowering-*.f64N/A
Simplified98.9%
Taylor expanded in l around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f6435.6%
Simplified35.6%
Taylor expanded in l around 0
cos-lowering-cos.f646.4%
Simplified6.4%
(FPCore (K m n M l) :precision binary64 1.0)
double code(double K, double m, double n, double M, double l) {
return 1.0;
}
real(8) function code(k, m, n, m_1, l)
real(8), intent (in) :: k
real(8), intent (in) :: m
real(8), intent (in) :: n
real(8), intent (in) :: m_1
real(8), intent (in) :: l
code = 1.0d0
end function
public static double code(double K, double m, double n, double M, double l) {
return 1.0;
}
def code(K, m, n, M, l): return 1.0
function code(K, m, n, M, l) return 1.0 end
function tmp = code(K, m, n, M, l) tmp = 1.0; end
code[K_, m_, n_, M_, l_] := 1.0
\begin{array}{l}
\\
1
\end{array}
Initial program 78.2%
*-lowering-*.f64N/A
sub-negN/A
+-commutativeN/A
neg-sub0N/A
associate-+l-N/A
sub0-negN/A
cos-negN/A
cos-lowering-cos.f64N/A
sub-negN/A
+-lowering-+.f64N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
metadata-evalN/A
exp-lowering-exp.f64N/A
Simplified78.2%
Taylor expanded in K around 0
*-commutativeN/A
*-lowering-*.f64N/A
Simplified98.9%
Taylor expanded in M around 0
exp-lowering-exp.f64N/A
associate--l+N/A
+-lowering-+.f64N/A
fabs-subN/A
sub-negN/A
mul-1-negN/A
fabs-lowering-fabs.f64N/A
mul-1-negN/A
sub-negN/A
--lowering--.f64N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
+-commutativeN/A
+-lowering-+.f6488.3%
Simplified88.3%
Taylor expanded in n around inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6457.1%
Simplified57.1%
Taylor expanded in n around 0
Simplified6.4%
herbie shell --seed 2024161
(FPCore (K m n M l)
:name "Maksimov and Kolovsky, Equation (32)"
:precision binary64
(* (cos (- (/ (* K (+ m n)) 2.0) M)) (exp (- (- (pow (- (/ (+ m n) 2.0) M) 2.0)) (- l (fabs (- m n)))))))