
(FPCore (i n) :precision binary64 (* 100.0 (/ (- (pow (+ 1.0 (/ i n)) n) 1.0) (/ i n))))
double code(double i, double n) {
return 100.0 * ((pow((1.0 + (i / n)), n) - 1.0) / (i / n));
}
real(8) function code(i, n)
real(8), intent (in) :: i
real(8), intent (in) :: n
code = 100.0d0 * ((((1.0d0 + (i / n)) ** n) - 1.0d0) / (i / n))
end function
public static double code(double i, double n) {
return 100.0 * ((Math.pow((1.0 + (i / n)), n) - 1.0) / (i / n));
}
def code(i, n): return 100.0 * ((math.pow((1.0 + (i / n)), n) - 1.0) / (i / n))
function code(i, n) return Float64(100.0 * Float64(Float64((Float64(1.0 + Float64(i / n)) ^ n) - 1.0) / Float64(i / n))) end
function tmp = code(i, n) tmp = 100.0 * ((((1.0 + (i / n)) ^ n) - 1.0) / (i / n)); end
code[i_, n_] := N[(100.0 * N[(N[(N[Power[N[(1.0 + N[(i / n), $MachinePrecision]), $MachinePrecision], n], $MachinePrecision] - 1.0), $MachinePrecision] / N[(i / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
100 \cdot \frac{{\left(1 + \frac{i}{n}\right)}^{n} - 1}{\frac{i}{n}}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 19 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (i n) :precision binary64 (* 100.0 (/ (- (pow (+ 1.0 (/ i n)) n) 1.0) (/ i n))))
double code(double i, double n) {
return 100.0 * ((pow((1.0 + (i / n)), n) - 1.0) / (i / n));
}
real(8) function code(i, n)
real(8), intent (in) :: i
real(8), intent (in) :: n
code = 100.0d0 * ((((1.0d0 + (i / n)) ** n) - 1.0d0) / (i / n))
end function
public static double code(double i, double n) {
return 100.0 * ((Math.pow((1.0 + (i / n)), n) - 1.0) / (i / n));
}
def code(i, n): return 100.0 * ((math.pow((1.0 + (i / n)), n) - 1.0) / (i / n))
function code(i, n) return Float64(100.0 * Float64(Float64((Float64(1.0 + Float64(i / n)) ^ n) - 1.0) / Float64(i / n))) end
function tmp = code(i, n) tmp = 100.0 * ((((1.0 + (i / n)) ^ n) - 1.0) / (i / n)); end
code[i_, n_] := N[(100.0 * N[(N[(N[Power[N[(1.0 + N[(i / n), $MachinePrecision]), $MachinePrecision], n], $MachinePrecision] - 1.0), $MachinePrecision] / N[(i / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
100 \cdot \frac{{\left(1 + \frac{i}{n}\right)}^{n} - 1}{\frac{i}{n}}
\end{array}
(FPCore (i n)
:precision binary64
(let* ((t_0 (/ (- (pow (+ 1.0 (/ i n)) n) 1.0) (/ i n))))
(if (<= t_0 0.0)
(* 100.0 (* (/ (expm1 (* (log1p (/ i n)) n)) i) n))
(if (<= t_0 INFINITY)
(+ (/ (* (pow (+ (/ i n) 1.0) n) 100.0) (/ i n)) (* (/ -100.0 i) n))
(* 100.0 n)))))
double code(double i, double n) {
double t_0 = (pow((1.0 + (i / n)), n) - 1.0) / (i / n);
double tmp;
if (t_0 <= 0.0) {
tmp = 100.0 * ((expm1((log1p((i / n)) * n)) / i) * n);
} else if (t_0 <= ((double) INFINITY)) {
tmp = ((pow(((i / n) + 1.0), n) * 100.0) / (i / n)) + ((-100.0 / i) * n);
} else {
tmp = 100.0 * n;
}
return tmp;
}
public static double code(double i, double n) {
double t_0 = (Math.pow((1.0 + (i / n)), n) - 1.0) / (i / n);
double tmp;
if (t_0 <= 0.0) {
tmp = 100.0 * ((Math.expm1((Math.log1p((i / n)) * n)) / i) * n);
} else if (t_0 <= Double.POSITIVE_INFINITY) {
tmp = ((Math.pow(((i / n) + 1.0), n) * 100.0) / (i / n)) + ((-100.0 / i) * n);
} else {
tmp = 100.0 * n;
}
return tmp;
}
def code(i, n): t_0 = (math.pow((1.0 + (i / n)), n) - 1.0) / (i / n) tmp = 0 if t_0 <= 0.0: tmp = 100.0 * ((math.expm1((math.log1p((i / n)) * n)) / i) * n) elif t_0 <= math.inf: tmp = ((math.pow(((i / n) + 1.0), n) * 100.0) / (i / n)) + ((-100.0 / i) * n) else: tmp = 100.0 * n return tmp
function code(i, n) t_0 = Float64(Float64((Float64(1.0 + Float64(i / n)) ^ n) - 1.0) / Float64(i / n)) tmp = 0.0 if (t_0 <= 0.0) tmp = Float64(100.0 * Float64(Float64(expm1(Float64(log1p(Float64(i / n)) * n)) / i) * n)); elseif (t_0 <= Inf) tmp = Float64(Float64(Float64((Float64(Float64(i / n) + 1.0) ^ n) * 100.0) / Float64(i / n)) + Float64(Float64(-100.0 / i) * n)); else tmp = Float64(100.0 * n); end return tmp end
code[i_, n_] := Block[{t$95$0 = N[(N[(N[Power[N[(1.0 + N[(i / n), $MachinePrecision]), $MachinePrecision], n], $MachinePrecision] - 1.0), $MachinePrecision] / N[(i / n), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, 0.0], N[(100.0 * N[(N[(N[(Exp[N[(N[Log[1 + N[(i / n), $MachinePrecision]], $MachinePrecision] * n), $MachinePrecision]] - 1), $MachinePrecision] / i), $MachinePrecision] * n), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, Infinity], N[(N[(N[(N[Power[N[(N[(i / n), $MachinePrecision] + 1.0), $MachinePrecision], n], $MachinePrecision] * 100.0), $MachinePrecision] / N[(i / n), $MachinePrecision]), $MachinePrecision] + N[(N[(-100.0 / i), $MachinePrecision] * n), $MachinePrecision]), $MachinePrecision], N[(100.0 * n), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{{\left(1 + \frac{i}{n}\right)}^{n} - 1}{\frac{i}{n}}\\
\mathbf{if}\;t\_0 \leq 0:\\
\;\;\;\;100 \cdot \left(\frac{\mathsf{expm1}\left(\mathsf{log1p}\left(\frac{i}{n}\right) \cdot n\right)}{i} \cdot n\right)\\
\mathbf{elif}\;t\_0 \leq \infty:\\
\;\;\;\;\frac{{\left(\frac{i}{n} + 1\right)}^{n} \cdot 100}{\frac{i}{n}} + \frac{-100}{i} \cdot n\\
\mathbf{else}:\\
\;\;\;\;100 \cdot n\\
\end{array}
\end{array}
if (/.f64 (-.f64 (pow.f64 (+.f64 #s(literal 1 binary64) (/.f64 i n)) n) #s(literal 1 binary64)) (/.f64 i n)) < -0.0Initial program 24.3%
lift-/.f64N/A
lift-/.f64N/A
associate-/r/N/A
lower-*.f64N/A
lower-/.f6424.1
lift--.f64N/A
lift-pow.f64N/A
pow-to-expN/A
lower-expm1.f64N/A
lower-*.f64N/A
lift-+.f64N/A
lower-log1p.f6497.4
Applied rewrites97.4%
if -0.0 < (/.f64 (-.f64 (pow.f64 (+.f64 #s(literal 1 binary64) (/.f64 i n)) n) #s(literal 1 binary64)) (/.f64 i n)) < +inf.0Initial program 96.3%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lift--.f64N/A
sub-negN/A
distribute-lft-inN/A
div-addN/A
*-commutativeN/A
lower-+.f64N/A
lower-/.f64N/A
lower-*.f64N/A
lift-+.f64N/A
+-commutativeN/A
lower-+.f64N/A
Applied rewrites96.4%
lift-/.f64N/A
lift-/.f64N/A
associate-/r/N/A
lower-*.f64N/A
lower-/.f6496.5
Applied rewrites96.5%
if +inf.0 < (/.f64 (-.f64 (pow.f64 (+.f64 #s(literal 1 binary64) (/.f64 i n)) n) #s(literal 1 binary64)) (/.f64 i n)) Initial program 0.0%
Taylor expanded in i around 0
lower-*.f6485.5
Applied rewrites85.5%
(FPCore (i n) :precision binary64 (if (or (<= n -5e-311) (not (<= n 2.6e-67))) (* (* (/ (expm1 i) i) 100.0) n) (* (* (* (- (log i) (log n)) n) (/ 100.0 i)) n)))
double code(double i, double n) {
double tmp;
if ((n <= -5e-311) || !(n <= 2.6e-67)) {
tmp = ((expm1(i) / i) * 100.0) * n;
} else {
tmp = (((log(i) - log(n)) * n) * (100.0 / i)) * n;
}
return tmp;
}
public static double code(double i, double n) {
double tmp;
if ((n <= -5e-311) || !(n <= 2.6e-67)) {
tmp = ((Math.expm1(i) / i) * 100.0) * n;
} else {
tmp = (((Math.log(i) - Math.log(n)) * n) * (100.0 / i)) * n;
}
return tmp;
}
def code(i, n): tmp = 0 if (n <= -5e-311) or not (n <= 2.6e-67): tmp = ((math.expm1(i) / i) * 100.0) * n else: tmp = (((math.log(i) - math.log(n)) * n) * (100.0 / i)) * n return tmp
function code(i, n) tmp = 0.0 if ((n <= -5e-311) || !(n <= 2.6e-67)) tmp = Float64(Float64(Float64(expm1(i) / i) * 100.0) * n); else tmp = Float64(Float64(Float64(Float64(log(i) - log(n)) * n) * Float64(100.0 / i)) * n); end return tmp end
code[i_, n_] := If[Or[LessEqual[n, -5e-311], N[Not[LessEqual[n, 2.6e-67]], $MachinePrecision]], N[(N[(N[(N[(Exp[i] - 1), $MachinePrecision] / i), $MachinePrecision] * 100.0), $MachinePrecision] * n), $MachinePrecision], N[(N[(N[(N[(N[Log[i], $MachinePrecision] - N[Log[n], $MachinePrecision]), $MachinePrecision] * n), $MachinePrecision] * N[(100.0 / i), $MachinePrecision]), $MachinePrecision] * n), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;n \leq -5 \cdot 10^{-311} \lor \neg \left(n \leq 2.6 \cdot 10^{-67}\right):\\
\;\;\;\;\left(\frac{\mathsf{expm1}\left(i\right)}{i} \cdot 100\right) \cdot n\\
\mathbf{else}:\\
\;\;\;\;\left(\left(\left(\log i - \log n\right) \cdot n\right) \cdot \frac{100}{i}\right) \cdot n\\
\end{array}
\end{array}
if n < -5.00000000000023e-311 or 2.5999999999999999e-67 < n Initial program 26.8%
Taylor expanded in n around inf
associate-*r/N/A
*-commutativeN/A
associate-*r*N/A
associate-*l/N/A
associate-*r/N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f64N/A
lower-expm1.f6481.8
Applied rewrites81.8%
if -5.00000000000023e-311 < n < 2.5999999999999999e-67Initial program 31.4%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lift-/.f64N/A
associate-/r/N/A
lower-*.f64N/A
*-commutativeN/A
associate-/l*N/A
lower-*.f64N/A
lift--.f64N/A
lift-pow.f64N/A
pow-to-expN/A
lower-expm1.f64N/A
lower-*.f64N/A
lift-+.f64N/A
lower-log1p.f64N/A
lower-/.f6478.1
Applied rewrites78.1%
Taylor expanded in n around 0
*-commutativeN/A
lower-*.f64N/A
mul-1-negN/A
unsub-negN/A
lower--.f64N/A
lower-log.f64N/A
lower-log.f6475.2
Applied rewrites75.2%
Final simplification81.0%
(FPCore (i n) :precision binary64 (if (or (<= n -5e-311) (not (<= n 2.6e-67))) (* (* (/ (expm1 i) i) 100.0) n) (* (* (* n (/ (- (log i) (log n)) i)) 100.0) n)))
double code(double i, double n) {
double tmp;
if ((n <= -5e-311) || !(n <= 2.6e-67)) {
tmp = ((expm1(i) / i) * 100.0) * n;
} else {
tmp = ((n * ((log(i) - log(n)) / i)) * 100.0) * n;
}
return tmp;
}
public static double code(double i, double n) {
double tmp;
if ((n <= -5e-311) || !(n <= 2.6e-67)) {
tmp = ((Math.expm1(i) / i) * 100.0) * n;
} else {
tmp = ((n * ((Math.log(i) - Math.log(n)) / i)) * 100.0) * n;
}
return tmp;
}
def code(i, n): tmp = 0 if (n <= -5e-311) or not (n <= 2.6e-67): tmp = ((math.expm1(i) / i) * 100.0) * n else: tmp = ((n * ((math.log(i) - math.log(n)) / i)) * 100.0) * n return tmp
function code(i, n) tmp = 0.0 if ((n <= -5e-311) || !(n <= 2.6e-67)) tmp = Float64(Float64(Float64(expm1(i) / i) * 100.0) * n); else tmp = Float64(Float64(Float64(n * Float64(Float64(log(i) - log(n)) / i)) * 100.0) * n); end return tmp end
code[i_, n_] := If[Or[LessEqual[n, -5e-311], N[Not[LessEqual[n, 2.6e-67]], $MachinePrecision]], N[(N[(N[(N[(Exp[i] - 1), $MachinePrecision] / i), $MachinePrecision] * 100.0), $MachinePrecision] * n), $MachinePrecision], N[(N[(N[(n * N[(N[(N[Log[i], $MachinePrecision] - N[Log[n], $MachinePrecision]), $MachinePrecision] / i), $MachinePrecision]), $MachinePrecision] * 100.0), $MachinePrecision] * n), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;n \leq -5 \cdot 10^{-311} \lor \neg \left(n \leq 2.6 \cdot 10^{-67}\right):\\
\;\;\;\;\left(\frac{\mathsf{expm1}\left(i\right)}{i} \cdot 100\right) \cdot n\\
\mathbf{else}:\\
\;\;\;\;\left(\left(n \cdot \frac{\log i - \log n}{i}\right) \cdot 100\right) \cdot n\\
\end{array}
\end{array}
if n < -5.00000000000023e-311 or 2.5999999999999999e-67 < n Initial program 26.8%
Taylor expanded in n around inf
associate-*r/N/A
*-commutativeN/A
associate-*r*N/A
associate-*l/N/A
associate-*r/N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f64N/A
lower-expm1.f6481.8
Applied rewrites81.8%
if -5.00000000000023e-311 < n < 2.5999999999999999e-67Initial program 31.4%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lift-/.f64N/A
associate-/r/N/A
lower-*.f64N/A
*-commutativeN/A
associate-/l*N/A
lower-*.f64N/A
lift--.f64N/A
lift-pow.f64N/A
pow-to-expN/A
lower-expm1.f64N/A
lower-*.f64N/A
lift-+.f64N/A
lower-log1p.f64N/A
lower-/.f6478.1
Applied rewrites78.1%
Taylor expanded in n around 0
*-commutativeN/A
lower-*.f64N/A
associate-/l*N/A
div-addN/A
associate-*r/N/A
+-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
associate-*r/N/A
div-addN/A
lower-/.f64N/A
mul-1-negN/A
unsub-negN/A
lower--.f64N/A
lower-log.f64N/A
lower-log.f6475.2
Applied rewrites75.2%
Final simplification81.0%
(FPCore (i n) :precision binary64 (if (or (<= n -5e-311) (not (<= n 2.6e-67))) (* (* (/ (expm1 i) i) 100.0) n) (* 100.0 (* (/ (* (- (log i) (log n)) n) i) n))))
double code(double i, double n) {
double tmp;
if ((n <= -5e-311) || !(n <= 2.6e-67)) {
tmp = ((expm1(i) / i) * 100.0) * n;
} else {
tmp = 100.0 * ((((log(i) - log(n)) * n) / i) * n);
}
return tmp;
}
public static double code(double i, double n) {
double tmp;
if ((n <= -5e-311) || !(n <= 2.6e-67)) {
tmp = ((Math.expm1(i) / i) * 100.0) * n;
} else {
tmp = 100.0 * ((((Math.log(i) - Math.log(n)) * n) / i) * n);
}
return tmp;
}
def code(i, n): tmp = 0 if (n <= -5e-311) or not (n <= 2.6e-67): tmp = ((math.expm1(i) / i) * 100.0) * n else: tmp = 100.0 * ((((math.log(i) - math.log(n)) * n) / i) * n) return tmp
function code(i, n) tmp = 0.0 if ((n <= -5e-311) || !(n <= 2.6e-67)) tmp = Float64(Float64(Float64(expm1(i) / i) * 100.0) * n); else tmp = Float64(100.0 * Float64(Float64(Float64(Float64(log(i) - log(n)) * n) / i) * n)); end return tmp end
code[i_, n_] := If[Or[LessEqual[n, -5e-311], N[Not[LessEqual[n, 2.6e-67]], $MachinePrecision]], N[(N[(N[(N[(Exp[i] - 1), $MachinePrecision] / i), $MachinePrecision] * 100.0), $MachinePrecision] * n), $MachinePrecision], N[(100.0 * N[(N[(N[(N[(N[Log[i], $MachinePrecision] - N[Log[n], $MachinePrecision]), $MachinePrecision] * n), $MachinePrecision] / i), $MachinePrecision] * n), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;n \leq -5 \cdot 10^{-311} \lor \neg \left(n \leq 2.6 \cdot 10^{-67}\right):\\
\;\;\;\;\left(\frac{\mathsf{expm1}\left(i\right)}{i} \cdot 100\right) \cdot n\\
\mathbf{else}:\\
\;\;\;\;100 \cdot \left(\frac{\left(\log i - \log n\right) \cdot n}{i} \cdot n\right)\\
\end{array}
\end{array}
if n < -5.00000000000023e-311 or 2.5999999999999999e-67 < n Initial program 26.8%
Taylor expanded in n around inf
associate-*r/N/A
*-commutativeN/A
associate-*r*N/A
associate-*l/N/A
associate-*r/N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f64N/A
lower-expm1.f6481.8
Applied rewrites81.8%
if -5.00000000000023e-311 < n < 2.5999999999999999e-67Initial program 31.4%
lift-/.f64N/A
lift-/.f64N/A
associate-/r/N/A
lower-*.f64N/A
lower-/.f6431.8
lift--.f64N/A
lift-pow.f64N/A
pow-to-expN/A
lower-expm1.f64N/A
lower-*.f64N/A
lift-+.f64N/A
lower-log1p.f6478.0
Applied rewrites78.0%
Taylor expanded in n around 0
*-commutativeN/A
lower-*.f64N/A
mul-1-negN/A
unsub-negN/A
lower--.f64N/A
lower-log.f64N/A
lower-log.f6475.1
Applied rewrites75.1%
Final simplification81.0%
(FPCore (i n) :precision binary64 (if (or (<= n -5e-311) (not (<= n 2.6e-67))) (* (* (/ (expm1 i) i) 100.0) n) (* 100.0 (* (* (/ (- (log i) (log n)) i) n) n))))
double code(double i, double n) {
double tmp;
if ((n <= -5e-311) || !(n <= 2.6e-67)) {
tmp = ((expm1(i) / i) * 100.0) * n;
} else {
tmp = 100.0 * ((((log(i) - log(n)) / i) * n) * n);
}
return tmp;
}
public static double code(double i, double n) {
double tmp;
if ((n <= -5e-311) || !(n <= 2.6e-67)) {
tmp = ((Math.expm1(i) / i) * 100.0) * n;
} else {
tmp = 100.0 * ((((Math.log(i) - Math.log(n)) / i) * n) * n);
}
return tmp;
}
def code(i, n): tmp = 0 if (n <= -5e-311) or not (n <= 2.6e-67): tmp = ((math.expm1(i) / i) * 100.0) * n else: tmp = 100.0 * ((((math.log(i) - math.log(n)) / i) * n) * n) return tmp
function code(i, n) tmp = 0.0 if ((n <= -5e-311) || !(n <= 2.6e-67)) tmp = Float64(Float64(Float64(expm1(i) / i) * 100.0) * n); else tmp = Float64(100.0 * Float64(Float64(Float64(Float64(log(i) - log(n)) / i) * n) * n)); end return tmp end
code[i_, n_] := If[Or[LessEqual[n, -5e-311], N[Not[LessEqual[n, 2.6e-67]], $MachinePrecision]], N[(N[(N[(N[(Exp[i] - 1), $MachinePrecision] / i), $MachinePrecision] * 100.0), $MachinePrecision] * n), $MachinePrecision], N[(100.0 * N[(N[(N[(N[(N[Log[i], $MachinePrecision] - N[Log[n], $MachinePrecision]), $MachinePrecision] / i), $MachinePrecision] * n), $MachinePrecision] * n), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;n \leq -5 \cdot 10^{-311} \lor \neg \left(n \leq 2.6 \cdot 10^{-67}\right):\\
\;\;\;\;\left(\frac{\mathsf{expm1}\left(i\right)}{i} \cdot 100\right) \cdot n\\
\mathbf{else}:\\
\;\;\;\;100 \cdot \left(\left(\frac{\log i - \log n}{i} \cdot n\right) \cdot n\right)\\
\end{array}
\end{array}
if n < -5.00000000000023e-311 or 2.5999999999999999e-67 < n Initial program 26.8%
Taylor expanded in n around inf
associate-*r/N/A
*-commutativeN/A
associate-*r*N/A
associate-*l/N/A
associate-*r/N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f64N/A
lower-expm1.f6481.8
Applied rewrites81.8%
if -5.00000000000023e-311 < n < 2.5999999999999999e-67Initial program 31.4%
lift-/.f64N/A
lift-/.f64N/A
associate-/r/N/A
lower-*.f64N/A
lower-/.f6431.8
lift--.f64N/A
lift-pow.f64N/A
pow-to-expN/A
lower-expm1.f64N/A
lower-*.f64N/A
lift-+.f64N/A
lower-log1p.f6478.0
Applied rewrites78.0%
Taylor expanded in n around 0
associate-/l*N/A
div-addN/A
associate-*r/N/A
+-commutativeN/A
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
associate-*r/N/A
div-addN/A
lower-/.f64N/A
mul-1-negN/A
unsub-negN/A
lower--.f64N/A
lower-log.f64N/A
lower-log.f6475.0
Applied rewrites75.0%
Final simplification81.0%
(FPCore (i n) :precision binary64 (if (or (<= n -5e-311) (not (<= n 2.6e-67))) (* (* (/ (expm1 i) i) 100.0) n) (* (* (* n n) 100.0) (/ (- (log i) (log n)) i))))
double code(double i, double n) {
double tmp;
if ((n <= -5e-311) || !(n <= 2.6e-67)) {
tmp = ((expm1(i) / i) * 100.0) * n;
} else {
tmp = ((n * n) * 100.0) * ((log(i) - log(n)) / i);
}
return tmp;
}
public static double code(double i, double n) {
double tmp;
if ((n <= -5e-311) || !(n <= 2.6e-67)) {
tmp = ((Math.expm1(i) / i) * 100.0) * n;
} else {
tmp = ((n * n) * 100.0) * ((Math.log(i) - Math.log(n)) / i);
}
return tmp;
}
def code(i, n): tmp = 0 if (n <= -5e-311) or not (n <= 2.6e-67): tmp = ((math.expm1(i) / i) * 100.0) * n else: tmp = ((n * n) * 100.0) * ((math.log(i) - math.log(n)) / i) return tmp
function code(i, n) tmp = 0.0 if ((n <= -5e-311) || !(n <= 2.6e-67)) tmp = Float64(Float64(Float64(expm1(i) / i) * 100.0) * n); else tmp = Float64(Float64(Float64(n * n) * 100.0) * Float64(Float64(log(i) - log(n)) / i)); end return tmp end
code[i_, n_] := If[Or[LessEqual[n, -5e-311], N[Not[LessEqual[n, 2.6e-67]], $MachinePrecision]], N[(N[(N[(N[(Exp[i] - 1), $MachinePrecision] / i), $MachinePrecision] * 100.0), $MachinePrecision] * n), $MachinePrecision], N[(N[(N[(n * n), $MachinePrecision] * 100.0), $MachinePrecision] * N[(N[(N[Log[i], $MachinePrecision] - N[Log[n], $MachinePrecision]), $MachinePrecision] / i), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;n \leq -5 \cdot 10^{-311} \lor \neg \left(n \leq 2.6 \cdot 10^{-67}\right):\\
\;\;\;\;\left(\frac{\mathsf{expm1}\left(i\right)}{i} \cdot 100\right) \cdot n\\
\mathbf{else}:\\
\;\;\;\;\left(\left(n \cdot n\right) \cdot 100\right) \cdot \frac{\log i - \log n}{i}\\
\end{array}
\end{array}
if n < -5.00000000000023e-311 or 2.5999999999999999e-67 < n Initial program 26.8%
Taylor expanded in n around inf
associate-*r/N/A
*-commutativeN/A
associate-*r*N/A
associate-*l/N/A
associate-*r/N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f64N/A
lower-expm1.f6481.8
Applied rewrites81.8%
if -5.00000000000023e-311 < n < 2.5999999999999999e-67Initial program 31.4%
Taylor expanded in n around 0
*-commutativeN/A
associate-/l*N/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
lower-/.f64N/A
mul-1-negN/A
unsub-negN/A
lower--.f64N/A
lower-log.f64N/A
lower-log.f6469.5
Applied rewrites69.5%
Final simplification80.3%
(FPCore (i n) :precision binary64 (if (or (<= n -2.1e-233) (not (<= n 1.95e-184))) (* (* (/ (expm1 i) i) 100.0) n) (* 100.0 (/ (- 1.0 1.0) (/ i n)))))
double code(double i, double n) {
double tmp;
if ((n <= -2.1e-233) || !(n <= 1.95e-184)) {
tmp = ((expm1(i) / i) * 100.0) * n;
} else {
tmp = 100.0 * ((1.0 - 1.0) / (i / n));
}
return tmp;
}
public static double code(double i, double n) {
double tmp;
if ((n <= -2.1e-233) || !(n <= 1.95e-184)) {
tmp = ((Math.expm1(i) / i) * 100.0) * n;
} else {
tmp = 100.0 * ((1.0 - 1.0) / (i / n));
}
return tmp;
}
def code(i, n): tmp = 0 if (n <= -2.1e-233) or not (n <= 1.95e-184): tmp = ((math.expm1(i) / i) * 100.0) * n else: tmp = 100.0 * ((1.0 - 1.0) / (i / n)) return tmp
function code(i, n) tmp = 0.0 if ((n <= -2.1e-233) || !(n <= 1.95e-184)) tmp = Float64(Float64(Float64(expm1(i) / i) * 100.0) * n); else tmp = Float64(100.0 * Float64(Float64(1.0 - 1.0) / Float64(i / n))); end return tmp end
code[i_, n_] := If[Or[LessEqual[n, -2.1e-233], N[Not[LessEqual[n, 1.95e-184]], $MachinePrecision]], N[(N[(N[(N[(Exp[i] - 1), $MachinePrecision] / i), $MachinePrecision] * 100.0), $MachinePrecision] * n), $MachinePrecision], N[(100.0 * N[(N[(1.0 - 1.0), $MachinePrecision] / N[(i / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;n \leq -2.1 \cdot 10^{-233} \lor \neg \left(n \leq 1.95 \cdot 10^{-184}\right):\\
\;\;\;\;\left(\frac{\mathsf{expm1}\left(i\right)}{i} \cdot 100\right) \cdot n\\
\mathbf{else}:\\
\;\;\;\;100 \cdot \frac{1 - 1}{\frac{i}{n}}\\
\end{array}
\end{array}
if n < -2.0999999999999999e-233 or 1.94999999999999997e-184 < n Initial program 23.8%
Taylor expanded in n around inf
associate-*r/N/A
*-commutativeN/A
associate-*r*N/A
associate-*l/N/A
associate-*r/N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f64N/A
lower-expm1.f6479.9
Applied rewrites79.9%
if -2.0999999999999999e-233 < n < 1.94999999999999997e-184Initial program 57.7%
Taylor expanded in i around 0
Applied rewrites76.3%
Final simplification79.6%
(FPCore (i n)
:precision binary64
(if (<= n -1.1e-193)
(*
(fma
(*
100.0
(+
(- 0.5 (/ (+ (fma -0.3333333333333333 (/ i n) (* 0.5 i)) 0.5) n))
(* 0.16666666666666666 i)))
i
100.0)
n)
(if (<= n 1.95e-184)
0.0
(*
(*
(/
(*
(fma
(fma (fma 0.041666666666666664 i 0.16666666666666666) i 0.5)
i
1.0)
i)
i)
100.0)
n))))
double code(double i, double n) {
double tmp;
if (n <= -1.1e-193) {
tmp = fma((100.0 * ((0.5 - ((fma(-0.3333333333333333, (i / n), (0.5 * i)) + 0.5) / n)) + (0.16666666666666666 * i))), i, 100.0) * n;
} else if (n <= 1.95e-184) {
tmp = 0.0;
} else {
tmp = (((fma(fma(fma(0.041666666666666664, i, 0.16666666666666666), i, 0.5), i, 1.0) * i) / i) * 100.0) * n;
}
return tmp;
}
function code(i, n) tmp = 0.0 if (n <= -1.1e-193) tmp = Float64(fma(Float64(100.0 * Float64(Float64(0.5 - Float64(Float64(fma(-0.3333333333333333, Float64(i / n), Float64(0.5 * i)) + 0.5) / n)) + Float64(0.16666666666666666 * i))), i, 100.0) * n); elseif (n <= 1.95e-184) tmp = 0.0; else tmp = Float64(Float64(Float64(Float64(fma(fma(fma(0.041666666666666664, i, 0.16666666666666666), i, 0.5), i, 1.0) * i) / i) * 100.0) * n); end return tmp end
code[i_, n_] := If[LessEqual[n, -1.1e-193], N[(N[(N[(100.0 * N[(N[(0.5 - N[(N[(N[(-0.3333333333333333 * N[(i / n), $MachinePrecision] + N[(0.5 * i), $MachinePrecision]), $MachinePrecision] + 0.5), $MachinePrecision] / n), $MachinePrecision]), $MachinePrecision] + N[(0.16666666666666666 * i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * i + 100.0), $MachinePrecision] * n), $MachinePrecision], If[LessEqual[n, 1.95e-184], 0.0, N[(N[(N[(N[(N[(N[(N[(0.041666666666666664 * i + 0.16666666666666666), $MachinePrecision] * i + 0.5), $MachinePrecision] * i + 1.0), $MachinePrecision] * i), $MachinePrecision] / i), $MachinePrecision] * 100.0), $MachinePrecision] * n), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;n \leq -1.1 \cdot 10^{-193}:\\
\;\;\;\;\mathsf{fma}\left(100 \cdot \left(\left(0.5 - \frac{\mathsf{fma}\left(-0.3333333333333333, \frac{i}{n}, 0.5 \cdot i\right) + 0.5}{n}\right) + 0.16666666666666666 \cdot i\right), i, 100\right) \cdot n\\
\mathbf{elif}\;n \leq 1.95 \cdot 10^{-184}:\\
\;\;\;\;0\\
\mathbf{else}:\\
\;\;\;\;\left(\frac{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(0.041666666666666664, i, 0.16666666666666666\right), i, 0.5\right), i, 1\right) \cdot i}{i} \cdot 100\right) \cdot n\\
\end{array}
\end{array}
if n < -1.09999999999999988e-193Initial program 24.9%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lift-/.f64N/A
associate-/r/N/A
lower-*.f64N/A
*-commutativeN/A
associate-/l*N/A
lower-*.f64N/A
lift--.f64N/A
lift-pow.f64N/A
pow-to-expN/A
lower-expm1.f64N/A
lower-*.f64N/A
lift-+.f64N/A
lower-log1p.f64N/A
lower-/.f6474.7
Applied rewrites74.7%
Taylor expanded in i around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
Applied rewrites54.8%
Taylor expanded in n around -inf
Applied rewrites58.2%
if -1.09999999999999988e-193 < n < 1.94999999999999997e-184Initial program 56.4%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lift--.f64N/A
sub-negN/A
distribute-lft-inN/A
div-addN/A
*-commutativeN/A
lower-+.f64N/A
lower-/.f64N/A
lower-*.f64N/A
lift-+.f64N/A
+-commutativeN/A
lower-+.f64N/A
Applied rewrites56.4%
Taylor expanded in i around 0
cancel-sign-sub-invN/A
metadata-evalN/A
+-inversesN/A
lower-/.f6470.0
Applied rewrites70.0%
Taylor expanded in i around 0
Applied rewrites70.0%
if 1.94999999999999997e-184 < n Initial program 19.6%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lift--.f64N/A
sub-negN/A
distribute-lft-inN/A
div-addN/A
*-commutativeN/A
lower-+.f64N/A
lower-/.f64N/A
lower-*.f64N/A
lift-+.f64N/A
+-commutativeN/A
lower-+.f64N/A
Applied rewrites19.6%
Taylor expanded in n around inf
distribute-lft-out--N/A
div-subN/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f64N/A
lower-expm1.f6484.6
Applied rewrites84.6%
Taylor expanded in i around 0
Applied rewrites70.3%
(FPCore (i n)
:precision binary64
(if (<= n -1.1e-193)
(*
100.0
(*
(fma
(+
(- 0.5 (/ (fma (+ 1.0 i) 0.5 (* -0.3333333333333333 (/ i n))) n))
(* 0.16666666666666666 i))
i
1.0)
n))
(if (<= n 1.95e-184)
0.0
(*
(*
(/
(*
(fma
(fma (fma 0.041666666666666664 i 0.16666666666666666) i 0.5)
i
1.0)
i)
i)
100.0)
n))))
double code(double i, double n) {
double tmp;
if (n <= -1.1e-193) {
tmp = 100.0 * (fma(((0.5 - (fma((1.0 + i), 0.5, (-0.3333333333333333 * (i / n))) / n)) + (0.16666666666666666 * i)), i, 1.0) * n);
} else if (n <= 1.95e-184) {
tmp = 0.0;
} else {
tmp = (((fma(fma(fma(0.041666666666666664, i, 0.16666666666666666), i, 0.5), i, 1.0) * i) / i) * 100.0) * n;
}
return tmp;
}
function code(i, n) tmp = 0.0 if (n <= -1.1e-193) tmp = Float64(100.0 * Float64(fma(Float64(Float64(0.5 - Float64(fma(Float64(1.0 + i), 0.5, Float64(-0.3333333333333333 * Float64(i / n))) / n)) + Float64(0.16666666666666666 * i)), i, 1.0) * n)); elseif (n <= 1.95e-184) tmp = 0.0; else tmp = Float64(Float64(Float64(Float64(fma(fma(fma(0.041666666666666664, i, 0.16666666666666666), i, 0.5), i, 1.0) * i) / i) * 100.0) * n); end return tmp end
code[i_, n_] := If[LessEqual[n, -1.1e-193], N[(100.0 * N[(N[(N[(N[(0.5 - N[(N[(N[(1.0 + i), $MachinePrecision] * 0.5 + N[(-0.3333333333333333 * N[(i / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / n), $MachinePrecision]), $MachinePrecision] + N[(0.16666666666666666 * i), $MachinePrecision]), $MachinePrecision] * i + 1.0), $MachinePrecision] * n), $MachinePrecision]), $MachinePrecision], If[LessEqual[n, 1.95e-184], 0.0, N[(N[(N[(N[(N[(N[(N[(0.041666666666666664 * i + 0.16666666666666666), $MachinePrecision] * i + 0.5), $MachinePrecision] * i + 1.0), $MachinePrecision] * i), $MachinePrecision] / i), $MachinePrecision] * 100.0), $MachinePrecision] * n), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;n \leq -1.1 \cdot 10^{-193}:\\
\;\;\;\;100 \cdot \left(\mathsf{fma}\left(\left(0.5 - \frac{\mathsf{fma}\left(1 + i, 0.5, -0.3333333333333333 \cdot \frac{i}{n}\right)}{n}\right) + 0.16666666666666666 \cdot i, i, 1\right) \cdot n\right)\\
\mathbf{elif}\;n \leq 1.95 \cdot 10^{-184}:\\
\;\;\;\;0\\
\mathbf{else}:\\
\;\;\;\;\left(\frac{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(0.041666666666666664, i, 0.16666666666666666\right), i, 0.5\right), i, 1\right) \cdot i}{i} \cdot 100\right) \cdot n\\
\end{array}
\end{array}
if n < -1.09999999999999988e-193Initial program 24.9%
lift-/.f64N/A
lift-/.f64N/A
associate-/r/N/A
lower-*.f64N/A
lower-/.f6424.8
lift--.f64N/A
lift-pow.f64N/A
pow-to-expN/A
lower-expm1.f64N/A
lower-*.f64N/A
lift-+.f64N/A
lower-log1p.f6475.4
Applied rewrites75.4%
Taylor expanded in i around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
Applied rewrites54.8%
Taylor expanded in n around -inf
Applied rewrites58.2%
if -1.09999999999999988e-193 < n < 1.94999999999999997e-184Initial program 56.4%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lift--.f64N/A
sub-negN/A
distribute-lft-inN/A
div-addN/A
*-commutativeN/A
lower-+.f64N/A
lower-/.f64N/A
lower-*.f64N/A
lift-+.f64N/A
+-commutativeN/A
lower-+.f64N/A
Applied rewrites56.4%
Taylor expanded in i around 0
cancel-sign-sub-invN/A
metadata-evalN/A
+-inversesN/A
lower-/.f6470.0
Applied rewrites70.0%
Taylor expanded in i around 0
Applied rewrites70.0%
if 1.94999999999999997e-184 < n Initial program 19.6%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lift--.f64N/A
sub-negN/A
distribute-lft-inN/A
div-addN/A
*-commutativeN/A
lower-+.f64N/A
lower-/.f64N/A
lower-*.f64N/A
lift-+.f64N/A
+-commutativeN/A
lower-+.f64N/A
Applied rewrites19.6%
Taylor expanded in n around inf
distribute-lft-out--N/A
div-subN/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f64N/A
lower-expm1.f6484.6
Applied rewrites84.6%
Taylor expanded in i around 0
Applied rewrites70.3%
(FPCore (i n)
:precision binary64
(if (<= n -1.1e-193)
(fma (* n (fma 16.666666666666668 i 50.0)) i (* 100.0 n))
(if (<= n 1.95e-184)
0.0
(*
(*
(/
(*
(fma
(fma (fma 0.041666666666666664 i 0.16666666666666666) i 0.5)
i
1.0)
i)
i)
100.0)
n))))
double code(double i, double n) {
double tmp;
if (n <= -1.1e-193) {
tmp = fma((n * fma(16.666666666666668, i, 50.0)), i, (100.0 * n));
} else if (n <= 1.95e-184) {
tmp = 0.0;
} else {
tmp = (((fma(fma(fma(0.041666666666666664, i, 0.16666666666666666), i, 0.5), i, 1.0) * i) / i) * 100.0) * n;
}
return tmp;
}
function code(i, n) tmp = 0.0 if (n <= -1.1e-193) tmp = fma(Float64(n * fma(16.666666666666668, i, 50.0)), i, Float64(100.0 * n)); elseif (n <= 1.95e-184) tmp = 0.0; else tmp = Float64(Float64(Float64(Float64(fma(fma(fma(0.041666666666666664, i, 0.16666666666666666), i, 0.5), i, 1.0) * i) / i) * 100.0) * n); end return tmp end
code[i_, n_] := If[LessEqual[n, -1.1e-193], N[(N[(n * N[(16.666666666666668 * i + 50.0), $MachinePrecision]), $MachinePrecision] * i + N[(100.0 * n), $MachinePrecision]), $MachinePrecision], If[LessEqual[n, 1.95e-184], 0.0, N[(N[(N[(N[(N[(N[(N[(0.041666666666666664 * i + 0.16666666666666666), $MachinePrecision] * i + 0.5), $MachinePrecision] * i + 1.0), $MachinePrecision] * i), $MachinePrecision] / i), $MachinePrecision] * 100.0), $MachinePrecision] * n), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;n \leq -1.1 \cdot 10^{-193}:\\
\;\;\;\;\mathsf{fma}\left(n \cdot \mathsf{fma}\left(16.666666666666668, i, 50\right), i, 100 \cdot n\right)\\
\mathbf{elif}\;n \leq 1.95 \cdot 10^{-184}:\\
\;\;\;\;0\\
\mathbf{else}:\\
\;\;\;\;\left(\frac{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(0.041666666666666664, i, 0.16666666666666666\right), i, 0.5\right), i, 1\right) \cdot i}{i} \cdot 100\right) \cdot n\\
\end{array}
\end{array}
if n < -1.09999999999999988e-193Initial program 24.9%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lift--.f64N/A
sub-negN/A
distribute-lft-inN/A
div-addN/A
*-commutativeN/A
lower-+.f64N/A
lower-/.f64N/A
lower-*.f64N/A
lift-+.f64N/A
+-commutativeN/A
lower-+.f64N/A
Applied rewrites24.9%
Taylor expanded in n around inf
distribute-lft-out--N/A
div-subN/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f64N/A
lower-expm1.f6476.8
Applied rewrites76.8%
Taylor expanded in i around 0
Applied rewrites58.0%
if -1.09999999999999988e-193 < n < 1.94999999999999997e-184Initial program 56.4%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lift--.f64N/A
sub-negN/A
distribute-lft-inN/A
div-addN/A
*-commutativeN/A
lower-+.f64N/A
lower-/.f64N/A
lower-*.f64N/A
lift-+.f64N/A
+-commutativeN/A
lower-+.f64N/A
Applied rewrites56.4%
Taylor expanded in i around 0
cancel-sign-sub-invN/A
metadata-evalN/A
+-inversesN/A
lower-/.f6470.0
Applied rewrites70.0%
Taylor expanded in i around 0
Applied rewrites70.0%
if 1.94999999999999997e-184 < n Initial program 19.6%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lift--.f64N/A
sub-negN/A
distribute-lft-inN/A
div-addN/A
*-commutativeN/A
lower-+.f64N/A
lower-/.f64N/A
lower-*.f64N/A
lift-+.f64N/A
+-commutativeN/A
lower-+.f64N/A
Applied rewrites19.6%
Taylor expanded in n around inf
distribute-lft-out--N/A
div-subN/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f64N/A
lower-expm1.f6484.6
Applied rewrites84.6%
Taylor expanded in i around 0
Applied rewrites70.3%
(FPCore (i n)
:precision binary64
(if (<= n -1.1e-193)
(fma (* n (fma 16.666666666666668 i 50.0)) i (* 100.0 n))
(if (<= n 1.95e-184)
0.0
(*
100.0
(*
(fma (fma (fma 0.041666666666666664 i 0.16666666666666666) i 0.5) i 1.0)
n)))))
double code(double i, double n) {
double tmp;
if (n <= -1.1e-193) {
tmp = fma((n * fma(16.666666666666668, i, 50.0)), i, (100.0 * n));
} else if (n <= 1.95e-184) {
tmp = 0.0;
} else {
tmp = 100.0 * (fma(fma(fma(0.041666666666666664, i, 0.16666666666666666), i, 0.5), i, 1.0) * n);
}
return tmp;
}
function code(i, n) tmp = 0.0 if (n <= -1.1e-193) tmp = fma(Float64(n * fma(16.666666666666668, i, 50.0)), i, Float64(100.0 * n)); elseif (n <= 1.95e-184) tmp = 0.0; else tmp = Float64(100.0 * Float64(fma(fma(fma(0.041666666666666664, i, 0.16666666666666666), i, 0.5), i, 1.0) * n)); end return tmp end
code[i_, n_] := If[LessEqual[n, -1.1e-193], N[(N[(n * N[(16.666666666666668 * i + 50.0), $MachinePrecision]), $MachinePrecision] * i + N[(100.0 * n), $MachinePrecision]), $MachinePrecision], If[LessEqual[n, 1.95e-184], 0.0, N[(100.0 * N[(N[(N[(N[(0.041666666666666664 * i + 0.16666666666666666), $MachinePrecision] * i + 0.5), $MachinePrecision] * i + 1.0), $MachinePrecision] * n), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;n \leq -1.1 \cdot 10^{-193}:\\
\;\;\;\;\mathsf{fma}\left(n \cdot \mathsf{fma}\left(16.666666666666668, i, 50\right), i, 100 \cdot n\right)\\
\mathbf{elif}\;n \leq 1.95 \cdot 10^{-184}:\\
\;\;\;\;0\\
\mathbf{else}:\\
\;\;\;\;100 \cdot \left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(0.041666666666666664, i, 0.16666666666666666\right), i, 0.5\right), i, 1\right) \cdot n\right)\\
\end{array}
\end{array}
if n < -1.09999999999999988e-193Initial program 24.9%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lift--.f64N/A
sub-negN/A
distribute-lft-inN/A
div-addN/A
*-commutativeN/A
lower-+.f64N/A
lower-/.f64N/A
lower-*.f64N/A
lift-+.f64N/A
+-commutativeN/A
lower-+.f64N/A
Applied rewrites24.9%
Taylor expanded in n around inf
distribute-lft-out--N/A
div-subN/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f64N/A
lower-expm1.f6476.8
Applied rewrites76.8%
Taylor expanded in i around 0
Applied rewrites58.0%
if -1.09999999999999988e-193 < n < 1.94999999999999997e-184Initial program 56.4%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lift--.f64N/A
sub-negN/A
distribute-lft-inN/A
div-addN/A
*-commutativeN/A
lower-+.f64N/A
lower-/.f64N/A
lower-*.f64N/A
lift-+.f64N/A
+-commutativeN/A
lower-+.f64N/A
Applied rewrites56.4%
Taylor expanded in i around 0
cancel-sign-sub-invN/A
metadata-evalN/A
+-inversesN/A
lower-/.f6470.0
Applied rewrites70.0%
Taylor expanded in i around 0
Applied rewrites70.0%
if 1.94999999999999997e-184 < n Initial program 19.6%
lift-/.f64N/A
lift-/.f64N/A
associate-/r/N/A
lower-*.f64N/A
lower-/.f6420.1
lift--.f64N/A
lift-pow.f64N/A
pow-to-expN/A
lower-expm1.f64N/A
lower-*.f64N/A
lift-+.f64N/A
lower-log1p.f6474.4
Applied rewrites74.4%
Taylor expanded in i around 0
Applied rewrites64.3%
Taylor expanded in n around inf
Applied rewrites69.4%
(FPCore (i n)
:precision binary64
(if (<= n -1.1e-193)
(fma (* n (fma 16.666666666666668 i 50.0)) i (* 100.0 n))
(if (<= n 1.95e-184)
0.0
(*
(fma (fma (fma 4.166666666666667 i 16.666666666666668) i 50.0) i 100.0)
n))))
double code(double i, double n) {
double tmp;
if (n <= -1.1e-193) {
tmp = fma((n * fma(16.666666666666668, i, 50.0)), i, (100.0 * n));
} else if (n <= 1.95e-184) {
tmp = 0.0;
} else {
tmp = fma(fma(fma(4.166666666666667, i, 16.666666666666668), i, 50.0), i, 100.0) * n;
}
return tmp;
}
function code(i, n) tmp = 0.0 if (n <= -1.1e-193) tmp = fma(Float64(n * fma(16.666666666666668, i, 50.0)), i, Float64(100.0 * n)); elseif (n <= 1.95e-184) tmp = 0.0; else tmp = Float64(fma(fma(fma(4.166666666666667, i, 16.666666666666668), i, 50.0), i, 100.0) * n); end return tmp end
code[i_, n_] := If[LessEqual[n, -1.1e-193], N[(N[(n * N[(16.666666666666668 * i + 50.0), $MachinePrecision]), $MachinePrecision] * i + N[(100.0 * n), $MachinePrecision]), $MachinePrecision], If[LessEqual[n, 1.95e-184], 0.0, N[(N[(N[(N[(4.166666666666667 * i + 16.666666666666668), $MachinePrecision] * i + 50.0), $MachinePrecision] * i + 100.0), $MachinePrecision] * n), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;n \leq -1.1 \cdot 10^{-193}:\\
\;\;\;\;\mathsf{fma}\left(n \cdot \mathsf{fma}\left(16.666666666666668, i, 50\right), i, 100 \cdot n\right)\\
\mathbf{elif}\;n \leq 1.95 \cdot 10^{-184}:\\
\;\;\;\;0\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(4.166666666666667, i, 16.666666666666668\right), i, 50\right), i, 100\right) \cdot n\\
\end{array}
\end{array}
if n < -1.09999999999999988e-193Initial program 24.9%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lift--.f64N/A
sub-negN/A
distribute-lft-inN/A
div-addN/A
*-commutativeN/A
lower-+.f64N/A
lower-/.f64N/A
lower-*.f64N/A
lift-+.f64N/A
+-commutativeN/A
lower-+.f64N/A
Applied rewrites24.9%
Taylor expanded in n around inf
distribute-lft-out--N/A
div-subN/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f64N/A
lower-expm1.f6476.8
Applied rewrites76.8%
Taylor expanded in i around 0
Applied rewrites58.0%
if -1.09999999999999988e-193 < n < 1.94999999999999997e-184Initial program 56.4%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lift--.f64N/A
sub-negN/A
distribute-lft-inN/A
div-addN/A
*-commutativeN/A
lower-+.f64N/A
lower-/.f64N/A
lower-*.f64N/A
lift-+.f64N/A
+-commutativeN/A
lower-+.f64N/A
Applied rewrites56.4%
Taylor expanded in i around 0
cancel-sign-sub-invN/A
metadata-evalN/A
+-inversesN/A
lower-/.f6470.0
Applied rewrites70.0%
Taylor expanded in i around 0
Applied rewrites70.0%
if 1.94999999999999997e-184 < n Initial program 19.6%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lift--.f64N/A
sub-negN/A
distribute-lft-inN/A
div-addN/A
*-commutativeN/A
lower-+.f64N/A
lower-/.f64N/A
lower-*.f64N/A
lift-+.f64N/A
+-commutativeN/A
lower-+.f64N/A
Applied rewrites19.6%
Taylor expanded in n around inf
distribute-lft-out--N/A
div-subN/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f64N/A
lower-expm1.f6484.6
Applied rewrites84.6%
Taylor expanded in i around 0
Applied rewrites69.4%
(FPCore (i n)
:precision binary64
(if (<= n -1.1e-193)
(fma (* n (fma 16.666666666666668 i 50.0)) i (* 100.0 n))
(if (<= n 1.95e-184)
0.0
(* 100.0 (* (fma (fma 0.16666666666666666 i 0.5) i 1.0) n)))))
double code(double i, double n) {
double tmp;
if (n <= -1.1e-193) {
tmp = fma((n * fma(16.666666666666668, i, 50.0)), i, (100.0 * n));
} else if (n <= 1.95e-184) {
tmp = 0.0;
} else {
tmp = 100.0 * (fma(fma(0.16666666666666666, i, 0.5), i, 1.0) * n);
}
return tmp;
}
function code(i, n) tmp = 0.0 if (n <= -1.1e-193) tmp = fma(Float64(n * fma(16.666666666666668, i, 50.0)), i, Float64(100.0 * n)); elseif (n <= 1.95e-184) tmp = 0.0; else tmp = Float64(100.0 * Float64(fma(fma(0.16666666666666666, i, 0.5), i, 1.0) * n)); end return tmp end
code[i_, n_] := If[LessEqual[n, -1.1e-193], N[(N[(n * N[(16.666666666666668 * i + 50.0), $MachinePrecision]), $MachinePrecision] * i + N[(100.0 * n), $MachinePrecision]), $MachinePrecision], If[LessEqual[n, 1.95e-184], 0.0, N[(100.0 * N[(N[(N[(0.16666666666666666 * i + 0.5), $MachinePrecision] * i + 1.0), $MachinePrecision] * n), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;n \leq -1.1 \cdot 10^{-193}:\\
\;\;\;\;\mathsf{fma}\left(n \cdot \mathsf{fma}\left(16.666666666666668, i, 50\right), i, 100 \cdot n\right)\\
\mathbf{elif}\;n \leq 1.95 \cdot 10^{-184}:\\
\;\;\;\;0\\
\mathbf{else}:\\
\;\;\;\;100 \cdot \left(\mathsf{fma}\left(\mathsf{fma}\left(0.16666666666666666, i, 0.5\right), i, 1\right) \cdot n\right)\\
\end{array}
\end{array}
if n < -1.09999999999999988e-193Initial program 24.9%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lift--.f64N/A
sub-negN/A
distribute-lft-inN/A
div-addN/A
*-commutativeN/A
lower-+.f64N/A
lower-/.f64N/A
lower-*.f64N/A
lift-+.f64N/A
+-commutativeN/A
lower-+.f64N/A
Applied rewrites24.9%
Taylor expanded in n around inf
distribute-lft-out--N/A
div-subN/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f64N/A
lower-expm1.f6476.8
Applied rewrites76.8%
Taylor expanded in i around 0
Applied rewrites58.0%
if -1.09999999999999988e-193 < n < 1.94999999999999997e-184Initial program 56.4%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lift--.f64N/A
sub-negN/A
distribute-lft-inN/A
div-addN/A
*-commutativeN/A
lower-+.f64N/A
lower-/.f64N/A
lower-*.f64N/A
lift-+.f64N/A
+-commutativeN/A
lower-+.f64N/A
Applied rewrites56.4%
Taylor expanded in i around 0
cancel-sign-sub-invN/A
metadata-evalN/A
+-inversesN/A
lower-/.f6470.0
Applied rewrites70.0%
Taylor expanded in i around 0
Applied rewrites70.0%
if 1.94999999999999997e-184 < n Initial program 19.6%
lift-/.f64N/A
lift-/.f64N/A
associate-/r/N/A
lower-*.f64N/A
lower-/.f6420.1
lift--.f64N/A
lift-pow.f64N/A
pow-to-expN/A
lower-expm1.f64N/A
lower-*.f64N/A
lift-+.f64N/A
lower-log1p.f6474.4
Applied rewrites74.4%
Taylor expanded in i around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
Applied rewrites64.8%
Taylor expanded in n around inf
Applied rewrites66.9%
(FPCore (i n) :precision binary64 (if (or (<= n -1.1e-193) (not (<= n 1.95e-184))) (* (fma (fma 16.666666666666668 i 50.0) i 100.0) n) 0.0))
double code(double i, double n) {
double tmp;
if ((n <= -1.1e-193) || !(n <= 1.95e-184)) {
tmp = fma(fma(16.666666666666668, i, 50.0), i, 100.0) * n;
} else {
tmp = 0.0;
}
return tmp;
}
function code(i, n) tmp = 0.0 if ((n <= -1.1e-193) || !(n <= 1.95e-184)) tmp = Float64(fma(fma(16.666666666666668, i, 50.0), i, 100.0) * n); else tmp = 0.0; end return tmp end
code[i_, n_] := If[Or[LessEqual[n, -1.1e-193], N[Not[LessEqual[n, 1.95e-184]], $MachinePrecision]], N[(N[(N[(16.666666666666668 * i + 50.0), $MachinePrecision] * i + 100.0), $MachinePrecision] * n), $MachinePrecision], 0.0]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;n \leq -1.1 \cdot 10^{-193} \lor \neg \left(n \leq 1.95 \cdot 10^{-184}\right):\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(16.666666666666668, i, 50\right), i, 100\right) \cdot n\\
\mathbf{else}:\\
\;\;\;\;0\\
\end{array}
\end{array}
if n < -1.09999999999999988e-193 or 1.94999999999999997e-184 < n Initial program 22.5%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lift--.f64N/A
sub-negN/A
distribute-lft-inN/A
div-addN/A
*-commutativeN/A
lower-+.f64N/A
lower-/.f64N/A
lower-*.f64N/A
lift-+.f64N/A
+-commutativeN/A
lower-+.f64N/A
Applied rewrites22.4%
Taylor expanded in n around inf
distribute-lft-out--N/A
div-subN/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f64N/A
lower-expm1.f6480.4
Applied rewrites80.4%
Taylor expanded in i around 0
Applied rewrites62.1%
if -1.09999999999999988e-193 < n < 1.94999999999999997e-184Initial program 56.4%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lift--.f64N/A
sub-negN/A
distribute-lft-inN/A
div-addN/A
*-commutativeN/A
lower-+.f64N/A
lower-/.f64N/A
lower-*.f64N/A
lift-+.f64N/A
+-commutativeN/A
lower-+.f64N/A
Applied rewrites56.4%
Taylor expanded in i around 0
cancel-sign-sub-invN/A
metadata-evalN/A
+-inversesN/A
lower-/.f6470.0
Applied rewrites70.0%
Taylor expanded in i around 0
Applied rewrites70.0%
Final simplification63.3%
(FPCore (i n)
:precision binary64
(if (<= n -1.1e-193)
(fma (* n (fma 16.666666666666668 i 50.0)) i (* 100.0 n))
(if (<= n 1.95e-184)
0.0
(* (fma (fma 16.666666666666668 i 50.0) i 100.0) n))))
double code(double i, double n) {
double tmp;
if (n <= -1.1e-193) {
tmp = fma((n * fma(16.666666666666668, i, 50.0)), i, (100.0 * n));
} else if (n <= 1.95e-184) {
tmp = 0.0;
} else {
tmp = fma(fma(16.666666666666668, i, 50.0), i, 100.0) * n;
}
return tmp;
}
function code(i, n) tmp = 0.0 if (n <= -1.1e-193) tmp = fma(Float64(n * fma(16.666666666666668, i, 50.0)), i, Float64(100.0 * n)); elseif (n <= 1.95e-184) tmp = 0.0; else tmp = Float64(fma(fma(16.666666666666668, i, 50.0), i, 100.0) * n); end return tmp end
code[i_, n_] := If[LessEqual[n, -1.1e-193], N[(N[(n * N[(16.666666666666668 * i + 50.0), $MachinePrecision]), $MachinePrecision] * i + N[(100.0 * n), $MachinePrecision]), $MachinePrecision], If[LessEqual[n, 1.95e-184], 0.0, N[(N[(N[(16.666666666666668 * i + 50.0), $MachinePrecision] * i + 100.0), $MachinePrecision] * n), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;n \leq -1.1 \cdot 10^{-193}:\\
\;\;\;\;\mathsf{fma}\left(n \cdot \mathsf{fma}\left(16.666666666666668, i, 50\right), i, 100 \cdot n\right)\\
\mathbf{elif}\;n \leq 1.95 \cdot 10^{-184}:\\
\;\;\;\;0\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(16.666666666666668, i, 50\right), i, 100\right) \cdot n\\
\end{array}
\end{array}
if n < -1.09999999999999988e-193Initial program 24.9%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lift--.f64N/A
sub-negN/A
distribute-lft-inN/A
div-addN/A
*-commutativeN/A
lower-+.f64N/A
lower-/.f64N/A
lower-*.f64N/A
lift-+.f64N/A
+-commutativeN/A
lower-+.f64N/A
Applied rewrites24.9%
Taylor expanded in n around inf
distribute-lft-out--N/A
div-subN/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f64N/A
lower-expm1.f6476.8
Applied rewrites76.8%
Taylor expanded in i around 0
Applied rewrites58.0%
if -1.09999999999999988e-193 < n < 1.94999999999999997e-184Initial program 56.4%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lift--.f64N/A
sub-negN/A
distribute-lft-inN/A
div-addN/A
*-commutativeN/A
lower-+.f64N/A
lower-/.f64N/A
lower-*.f64N/A
lift-+.f64N/A
+-commutativeN/A
lower-+.f64N/A
Applied rewrites56.4%
Taylor expanded in i around 0
cancel-sign-sub-invN/A
metadata-evalN/A
+-inversesN/A
lower-/.f6470.0
Applied rewrites70.0%
Taylor expanded in i around 0
Applied rewrites70.0%
if 1.94999999999999997e-184 < n Initial program 19.6%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lift--.f64N/A
sub-negN/A
distribute-lft-inN/A
div-addN/A
*-commutativeN/A
lower-+.f64N/A
lower-/.f64N/A
lower-*.f64N/A
lift-+.f64N/A
+-commutativeN/A
lower-+.f64N/A
Applied rewrites19.6%
Taylor expanded in n around inf
distribute-lft-out--N/A
div-subN/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f64N/A
lower-expm1.f6484.6
Applied rewrites84.6%
Taylor expanded in i around 0
Applied rewrites66.9%
(FPCore (i n) :precision binary64 (if (<= n -1.1e-193) (* (fma 50.0 i 100.0) n) (if (<= n 1.95e-184) 0.0 (* (* (fma 0.5 i 1.0) 100.0) n))))
double code(double i, double n) {
double tmp;
if (n <= -1.1e-193) {
tmp = fma(50.0, i, 100.0) * n;
} else if (n <= 1.95e-184) {
tmp = 0.0;
} else {
tmp = (fma(0.5, i, 1.0) * 100.0) * n;
}
return tmp;
}
function code(i, n) tmp = 0.0 if (n <= -1.1e-193) tmp = Float64(fma(50.0, i, 100.0) * n); elseif (n <= 1.95e-184) tmp = 0.0; else tmp = Float64(Float64(fma(0.5, i, 1.0) * 100.0) * n); end return tmp end
code[i_, n_] := If[LessEqual[n, -1.1e-193], N[(N[(50.0 * i + 100.0), $MachinePrecision] * n), $MachinePrecision], If[LessEqual[n, 1.95e-184], 0.0, N[(N[(N[(0.5 * i + 1.0), $MachinePrecision] * 100.0), $MachinePrecision] * n), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;n \leq -1.1 \cdot 10^{-193}:\\
\;\;\;\;\mathsf{fma}\left(50, i, 100\right) \cdot n\\
\mathbf{elif}\;n \leq 1.95 \cdot 10^{-184}:\\
\;\;\;\;0\\
\mathbf{else}:\\
\;\;\;\;\left(\mathsf{fma}\left(0.5, i, 1\right) \cdot 100\right) \cdot n\\
\end{array}
\end{array}
if n < -1.09999999999999988e-193Initial program 24.9%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lift--.f64N/A
sub-negN/A
distribute-lft-inN/A
div-addN/A
*-commutativeN/A
lower-+.f64N/A
lower-/.f64N/A
lower-*.f64N/A
lift-+.f64N/A
+-commutativeN/A
lower-+.f64N/A
Applied rewrites24.9%
Taylor expanded in n around inf
distribute-lft-out--N/A
div-subN/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f64N/A
lower-expm1.f6476.8
Applied rewrites76.8%
Taylor expanded in i around 0
Applied rewrites56.6%
if -1.09999999999999988e-193 < n < 1.94999999999999997e-184Initial program 56.4%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lift--.f64N/A
sub-negN/A
distribute-lft-inN/A
div-addN/A
*-commutativeN/A
lower-+.f64N/A
lower-/.f64N/A
lower-*.f64N/A
lift-+.f64N/A
+-commutativeN/A
lower-+.f64N/A
Applied rewrites56.4%
Taylor expanded in i around 0
cancel-sign-sub-invN/A
metadata-evalN/A
+-inversesN/A
lower-/.f6470.0
Applied rewrites70.0%
Taylor expanded in i around 0
Applied rewrites70.0%
if 1.94999999999999997e-184 < n Initial program 19.6%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lift--.f64N/A
sub-negN/A
distribute-lft-inN/A
div-addN/A
*-commutativeN/A
lower-+.f64N/A
lower-/.f64N/A
lower-*.f64N/A
lift-+.f64N/A
+-commutativeN/A
lower-+.f64N/A
Applied rewrites19.6%
Taylor expanded in n around inf
distribute-lft-out--N/A
div-subN/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f64N/A
lower-expm1.f6484.6
Applied rewrites84.6%
Taylor expanded in i around 0
Applied rewrites65.5%
(FPCore (i n) :precision binary64 (if (or (<= n -1.1e-193) (not (<= n 1.95e-184))) (* (fma 50.0 i 100.0) n) 0.0))
double code(double i, double n) {
double tmp;
if ((n <= -1.1e-193) || !(n <= 1.95e-184)) {
tmp = fma(50.0, i, 100.0) * n;
} else {
tmp = 0.0;
}
return tmp;
}
function code(i, n) tmp = 0.0 if ((n <= -1.1e-193) || !(n <= 1.95e-184)) tmp = Float64(fma(50.0, i, 100.0) * n); else tmp = 0.0; end return tmp end
code[i_, n_] := If[Or[LessEqual[n, -1.1e-193], N[Not[LessEqual[n, 1.95e-184]], $MachinePrecision]], N[(N[(50.0 * i + 100.0), $MachinePrecision] * n), $MachinePrecision], 0.0]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;n \leq -1.1 \cdot 10^{-193} \lor \neg \left(n \leq 1.95 \cdot 10^{-184}\right):\\
\;\;\;\;\mathsf{fma}\left(50, i, 100\right) \cdot n\\
\mathbf{else}:\\
\;\;\;\;0\\
\end{array}
\end{array}
if n < -1.09999999999999988e-193 or 1.94999999999999997e-184 < n Initial program 22.5%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lift--.f64N/A
sub-negN/A
distribute-lft-inN/A
div-addN/A
*-commutativeN/A
lower-+.f64N/A
lower-/.f64N/A
lower-*.f64N/A
lift-+.f64N/A
+-commutativeN/A
lower-+.f64N/A
Applied rewrites22.4%
Taylor expanded in n around inf
distribute-lft-out--N/A
div-subN/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f64N/A
lower-expm1.f6480.4
Applied rewrites80.4%
Taylor expanded in i around 0
Applied rewrites60.7%
if -1.09999999999999988e-193 < n < 1.94999999999999997e-184Initial program 56.4%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lift--.f64N/A
sub-negN/A
distribute-lft-inN/A
div-addN/A
*-commutativeN/A
lower-+.f64N/A
lower-/.f64N/A
lower-*.f64N/A
lift-+.f64N/A
+-commutativeN/A
lower-+.f64N/A
Applied rewrites56.4%
Taylor expanded in i around 0
cancel-sign-sub-invN/A
metadata-evalN/A
+-inversesN/A
lower-/.f6470.0
Applied rewrites70.0%
Taylor expanded in i around 0
Applied rewrites70.0%
Final simplification62.1%
(FPCore (i n) :precision binary64 (if (<= i -8.8e+21) 0.0 (if (<= i 1.25e-15) (* 100.0 n) 0.0)))
double code(double i, double n) {
double tmp;
if (i <= -8.8e+21) {
tmp = 0.0;
} else if (i <= 1.25e-15) {
tmp = 100.0 * n;
} else {
tmp = 0.0;
}
return tmp;
}
real(8) function code(i, n)
real(8), intent (in) :: i
real(8), intent (in) :: n
real(8) :: tmp
if (i <= (-8.8d+21)) then
tmp = 0.0d0
else if (i <= 1.25d-15) then
tmp = 100.0d0 * n
else
tmp = 0.0d0
end if
code = tmp
end function
public static double code(double i, double n) {
double tmp;
if (i <= -8.8e+21) {
tmp = 0.0;
} else if (i <= 1.25e-15) {
tmp = 100.0 * n;
} else {
tmp = 0.0;
}
return tmp;
}
def code(i, n): tmp = 0 if i <= -8.8e+21: tmp = 0.0 elif i <= 1.25e-15: tmp = 100.0 * n else: tmp = 0.0 return tmp
function code(i, n) tmp = 0.0 if (i <= -8.8e+21) tmp = 0.0; elseif (i <= 1.25e-15) tmp = Float64(100.0 * n); else tmp = 0.0; end return tmp end
function tmp_2 = code(i, n) tmp = 0.0; if (i <= -8.8e+21) tmp = 0.0; elseif (i <= 1.25e-15) tmp = 100.0 * n; else tmp = 0.0; end tmp_2 = tmp; end
code[i_, n_] := If[LessEqual[i, -8.8e+21], 0.0, If[LessEqual[i, 1.25e-15], N[(100.0 * n), $MachinePrecision], 0.0]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;i \leq -8.8 \cdot 10^{+21}:\\
\;\;\;\;0\\
\mathbf{elif}\;i \leq 1.25 \cdot 10^{-15}:\\
\;\;\;\;100 \cdot n\\
\mathbf{else}:\\
\;\;\;\;0\\
\end{array}
\end{array}
if i < -8.8e21 or 1.25e-15 < i Initial program 53.1%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lift--.f64N/A
sub-negN/A
distribute-lft-inN/A
div-addN/A
*-commutativeN/A
lower-+.f64N/A
lower-/.f64N/A
lower-*.f64N/A
lift-+.f64N/A
+-commutativeN/A
lower-+.f64N/A
Applied rewrites53.1%
Taylor expanded in i around 0
cancel-sign-sub-invN/A
metadata-evalN/A
+-inversesN/A
lower-/.f6426.4
Applied rewrites26.4%
Taylor expanded in i around 0
Applied rewrites26.4%
if -8.8e21 < i < 1.25e-15Initial program 8.0%
Taylor expanded in i around 0
lower-*.f6480.3
Applied rewrites80.3%
(FPCore (i n) :precision binary64 0.0)
double code(double i, double n) {
return 0.0;
}
real(8) function code(i, n)
real(8), intent (in) :: i
real(8), intent (in) :: n
code = 0.0d0
end function
public static double code(double i, double n) {
return 0.0;
}
def code(i, n): return 0.0
function code(i, n) return 0.0 end
function tmp = code(i, n) tmp = 0.0; end
code[i_, n_] := 0.0
\begin{array}{l}
\\
0
\end{array}
Initial program 27.4%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lift--.f64N/A
sub-negN/A
distribute-lft-inN/A
div-addN/A
*-commutativeN/A
lower-+.f64N/A
lower-/.f64N/A
lower-*.f64N/A
lift-+.f64N/A
+-commutativeN/A
lower-+.f64N/A
Applied rewrites27.3%
Taylor expanded in i around 0
cancel-sign-sub-invN/A
metadata-evalN/A
+-inversesN/A
lower-/.f6415.9
Applied rewrites15.9%
Taylor expanded in i around 0
Applied rewrites15.9%
(FPCore (i n)
:precision binary64
(let* ((t_0 (+ 1.0 (/ i n))))
(*
100.0
(/
(-
(exp
(*
n
(if (== t_0 1.0)
(/ i n)
(/ (* (/ i n) (log t_0)) (- (+ (/ i n) 1.0) 1.0)))))
1.0)
(/ i n)))))
double code(double i, double n) {
double t_0 = 1.0 + (i / n);
double tmp;
if (t_0 == 1.0) {
tmp = i / n;
} else {
tmp = ((i / n) * log(t_0)) / (((i / n) + 1.0) - 1.0);
}
return 100.0 * ((exp((n * tmp)) - 1.0) / (i / n));
}
real(8) function code(i, n)
real(8), intent (in) :: i
real(8), intent (in) :: n
real(8) :: t_0
real(8) :: tmp
t_0 = 1.0d0 + (i / n)
if (t_0 == 1.0d0) then
tmp = i / n
else
tmp = ((i / n) * log(t_0)) / (((i / n) + 1.0d0) - 1.0d0)
end if
code = 100.0d0 * ((exp((n * tmp)) - 1.0d0) / (i / n))
end function
public static double code(double i, double n) {
double t_0 = 1.0 + (i / n);
double tmp;
if (t_0 == 1.0) {
tmp = i / n;
} else {
tmp = ((i / n) * Math.log(t_0)) / (((i / n) + 1.0) - 1.0);
}
return 100.0 * ((Math.exp((n * tmp)) - 1.0) / (i / n));
}
def code(i, n): t_0 = 1.0 + (i / n) tmp = 0 if t_0 == 1.0: tmp = i / n else: tmp = ((i / n) * math.log(t_0)) / (((i / n) + 1.0) - 1.0) return 100.0 * ((math.exp((n * tmp)) - 1.0) / (i / n))
function code(i, n) t_0 = Float64(1.0 + Float64(i / n)) tmp = 0.0 if (t_0 == 1.0) tmp = Float64(i / n); else tmp = Float64(Float64(Float64(i / n) * log(t_0)) / Float64(Float64(Float64(i / n) + 1.0) - 1.0)); end return Float64(100.0 * Float64(Float64(exp(Float64(n * tmp)) - 1.0) / Float64(i / n))) end
function tmp_2 = code(i, n) t_0 = 1.0 + (i / n); tmp = 0.0; if (t_0 == 1.0) tmp = i / n; else tmp = ((i / n) * log(t_0)) / (((i / n) + 1.0) - 1.0); end tmp_2 = 100.0 * ((exp((n * tmp)) - 1.0) / (i / n)); end
code[i_, n_] := Block[{t$95$0 = N[(1.0 + N[(i / n), $MachinePrecision]), $MachinePrecision]}, N[(100.0 * N[(N[(N[Exp[N[(n * If[Equal[t$95$0, 1.0], N[(i / n), $MachinePrecision], N[(N[(N[(i / n), $MachinePrecision] * N[Log[t$95$0], $MachinePrecision]), $MachinePrecision] / N[(N[(N[(i / n), $MachinePrecision] + 1.0), $MachinePrecision] - 1.0), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]], $MachinePrecision] - 1.0), $MachinePrecision] / N[(i / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := 1 + \frac{i}{n}\\
100 \cdot \frac{e^{n \cdot \begin{array}{l}
\mathbf{if}\;t\_0 = 1:\\
\;\;\;\;\frac{i}{n}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{i}{n} \cdot \log t\_0}{\left(\frac{i}{n} + 1\right) - 1}\\
\end{array}} - 1}{\frac{i}{n}}
\end{array}
\end{array}
herbie shell --seed 2024305
(FPCore (i n)
:name "Compound Interest"
:precision binary64
:alt
(! :herbie-platform default (let ((lnbase (if (== (+ 1 (/ i n)) 1) (/ i n) (/ (* (/ i n) (log (+ 1 (/ i n)))) (- (+ (/ i n) 1) 1))))) (* 100 (/ (- (exp (* n lnbase)) 1) (/ i n)))))
(* 100.0 (/ (- (pow (+ 1.0 (/ i n)) n) 1.0) (/ i n))))