
(FPCore (x) :precision binary64 (/ (* 6.0 (- x 1.0)) (+ (+ x 1.0) (* 4.0 (sqrt x)))))
double code(double x) {
return (6.0 * (x - 1.0)) / ((x + 1.0) + (4.0 * sqrt(x)));
}
real(8) function code(x)
real(8), intent (in) :: x
code = (6.0d0 * (x - 1.0d0)) / ((x + 1.0d0) + (4.0d0 * sqrt(x)))
end function
public static double code(double x) {
return (6.0 * (x - 1.0)) / ((x + 1.0) + (4.0 * Math.sqrt(x)));
}
def code(x): return (6.0 * (x - 1.0)) / ((x + 1.0) + (4.0 * math.sqrt(x)))
function code(x) return Float64(Float64(6.0 * Float64(x - 1.0)) / Float64(Float64(x + 1.0) + Float64(4.0 * sqrt(x)))) end
function tmp = code(x) tmp = (6.0 * (x - 1.0)) / ((x + 1.0) + (4.0 * sqrt(x))); end
code[x_] := N[(N[(6.0 * N[(x - 1.0), $MachinePrecision]), $MachinePrecision] / N[(N[(x + 1.0), $MachinePrecision] + N[(4.0 * N[Sqrt[x], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{6 \cdot \left(x - 1\right)}{\left(x + 1\right) + 4 \cdot \sqrt{x}}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 14 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x) :precision binary64 (/ (* 6.0 (- x 1.0)) (+ (+ x 1.0) (* 4.0 (sqrt x)))))
double code(double x) {
return (6.0 * (x - 1.0)) / ((x + 1.0) + (4.0 * sqrt(x)));
}
real(8) function code(x)
real(8), intent (in) :: x
code = (6.0d0 * (x - 1.0d0)) / ((x + 1.0d0) + (4.0d0 * sqrt(x)))
end function
public static double code(double x) {
return (6.0 * (x - 1.0)) / ((x + 1.0) + (4.0 * Math.sqrt(x)));
}
def code(x): return (6.0 * (x - 1.0)) / ((x + 1.0) + (4.0 * math.sqrt(x)))
function code(x) return Float64(Float64(6.0 * Float64(x - 1.0)) / Float64(Float64(x + 1.0) + Float64(4.0 * sqrt(x)))) end
function tmp = code(x) tmp = (6.0 * (x - 1.0)) / ((x + 1.0) + (4.0 * sqrt(x))); end
code[x_] := N[(N[(6.0 * N[(x - 1.0), $MachinePrecision]), $MachinePrecision] / N[(N[(x + 1.0), $MachinePrecision] + N[(4.0 * N[Sqrt[x], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{6 \cdot \left(x - 1\right)}{\left(x + 1\right) + 4 \cdot \sqrt{x}}
\end{array}
(FPCore (x) :precision binary64 (/ 6.0 (/ (+ x (/ (- 1.0 (* x 16.0)) (- 1.0 (* 4.0 (sqrt x))))) (+ x -1.0))))
double code(double x) {
return 6.0 / ((x + ((1.0 - (x * 16.0)) / (1.0 - (4.0 * sqrt(x))))) / (x + -1.0));
}
real(8) function code(x)
real(8), intent (in) :: x
code = 6.0d0 / ((x + ((1.0d0 - (x * 16.0d0)) / (1.0d0 - (4.0d0 * sqrt(x))))) / (x + (-1.0d0)))
end function
public static double code(double x) {
return 6.0 / ((x + ((1.0 - (x * 16.0)) / (1.0 - (4.0 * Math.sqrt(x))))) / (x + -1.0));
}
def code(x): return 6.0 / ((x + ((1.0 - (x * 16.0)) / (1.0 - (4.0 * math.sqrt(x))))) / (x + -1.0))
function code(x) return Float64(6.0 / Float64(Float64(x + Float64(Float64(1.0 - Float64(x * 16.0)) / Float64(1.0 - Float64(4.0 * sqrt(x))))) / Float64(x + -1.0))) end
function tmp = code(x) tmp = 6.0 / ((x + ((1.0 - (x * 16.0)) / (1.0 - (4.0 * sqrt(x))))) / (x + -1.0)); end
code[x_] := N[(6.0 / N[(N[(x + N[(N[(1.0 - N[(x * 16.0), $MachinePrecision]), $MachinePrecision] / N[(1.0 - N[(4.0 * N[Sqrt[x], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(x + -1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{6}{\frac{x + \frac{1 - x \cdot 16}{1 - 4 \cdot \sqrt{x}}}{x + -1}}
\end{array}
Initial program 99.8%
associate-/l*N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
associate-+l+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-eval99.9%
Applied egg-rr99.9%
flip-+N/A
/-lowering-/.f64N/A
metadata-evalN/A
--lowering--.f64N/A
*-commutativeN/A
*-commutativeN/A
swap-sqrN/A
metadata-evalN/A
metadata-evalN/A
rem-square-sqrtN/A
*-lowering-*.f64N/A
metadata-evalN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f6499.9%
Applied egg-rr99.9%
(FPCore (x)
:precision binary64
(let* ((t_0 (* 4.0 (sqrt x))))
(if (<= x 1.0)
(/ 6.0 (/ (+ 1.0 t_0) (+ x -1.0)))
(* 6.0 (/ (+ x -1.0) (+ x t_0))))))
double code(double x) {
double t_0 = 4.0 * sqrt(x);
double tmp;
if (x <= 1.0) {
tmp = 6.0 / ((1.0 + t_0) / (x + -1.0));
} else {
tmp = 6.0 * ((x + -1.0) / (x + t_0));
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: t_0
real(8) :: tmp
t_0 = 4.0d0 * sqrt(x)
if (x <= 1.0d0) then
tmp = 6.0d0 / ((1.0d0 + t_0) / (x + (-1.0d0)))
else
tmp = 6.0d0 * ((x + (-1.0d0)) / (x + t_0))
end if
code = tmp
end function
public static double code(double x) {
double t_0 = 4.0 * Math.sqrt(x);
double tmp;
if (x <= 1.0) {
tmp = 6.0 / ((1.0 + t_0) / (x + -1.0));
} else {
tmp = 6.0 * ((x + -1.0) / (x + t_0));
}
return tmp;
}
def code(x): t_0 = 4.0 * math.sqrt(x) tmp = 0 if x <= 1.0: tmp = 6.0 / ((1.0 + t_0) / (x + -1.0)) else: tmp = 6.0 * ((x + -1.0) / (x + t_0)) return tmp
function code(x) t_0 = Float64(4.0 * sqrt(x)) tmp = 0.0 if (x <= 1.0) tmp = Float64(6.0 / Float64(Float64(1.0 + t_0) / Float64(x + -1.0))); else tmp = Float64(6.0 * Float64(Float64(x + -1.0) / Float64(x + t_0))); end return tmp end
function tmp_2 = code(x) t_0 = 4.0 * sqrt(x); tmp = 0.0; if (x <= 1.0) tmp = 6.0 / ((1.0 + t_0) / (x + -1.0)); else tmp = 6.0 * ((x + -1.0) / (x + t_0)); end tmp_2 = tmp; end
code[x_] := Block[{t$95$0 = N[(4.0 * N[Sqrt[x], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[x, 1.0], N[(6.0 / N[(N[(1.0 + t$95$0), $MachinePrecision] / N[(x + -1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(6.0 * N[(N[(x + -1.0), $MachinePrecision] / N[(x + t$95$0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := 4 \cdot \sqrt{x}\\
\mathbf{if}\;x \leq 1:\\
\;\;\;\;\frac{6}{\frac{1 + t\_0}{x + -1}}\\
\mathbf{else}:\\
\;\;\;\;6 \cdot \frac{x + -1}{x + t\_0}\\
\end{array}
\end{array}
if x < 1Initial program 99.9%
associate-/l*N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
associate-+l+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-eval99.9%
Applied egg-rr99.9%
Taylor expanded in x around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f6499.3%
Simplified99.3%
if 1 < x Initial program 99.7%
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
associate-+l+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f6499.9%
Applied egg-rr99.9%
Taylor expanded in x around inf
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f6498.0%
Simplified98.0%
Final simplification98.5%
(FPCore (x) :precision binary64 (let* ((t_0 (+ x (* 4.0 (sqrt x))))) (if (<= x 0.29) (/ -6.0 (+ 1.0 t_0)) (* 6.0 (/ (+ x -1.0) t_0)))))
double code(double x) {
double t_0 = x + (4.0 * sqrt(x));
double tmp;
if (x <= 0.29) {
tmp = -6.0 / (1.0 + t_0);
} else {
tmp = 6.0 * ((x + -1.0) / t_0);
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: t_0
real(8) :: tmp
t_0 = x + (4.0d0 * sqrt(x))
if (x <= 0.29d0) then
tmp = (-6.0d0) / (1.0d0 + t_0)
else
tmp = 6.0d0 * ((x + (-1.0d0)) / t_0)
end if
code = tmp
end function
public static double code(double x) {
double t_0 = x + (4.0 * Math.sqrt(x));
double tmp;
if (x <= 0.29) {
tmp = -6.0 / (1.0 + t_0);
} else {
tmp = 6.0 * ((x + -1.0) / t_0);
}
return tmp;
}
def code(x): t_0 = x + (4.0 * math.sqrt(x)) tmp = 0 if x <= 0.29: tmp = -6.0 / (1.0 + t_0) else: tmp = 6.0 * ((x + -1.0) / t_0) return tmp
function code(x) t_0 = Float64(x + Float64(4.0 * sqrt(x))) tmp = 0.0 if (x <= 0.29) tmp = Float64(-6.0 / Float64(1.0 + t_0)); else tmp = Float64(6.0 * Float64(Float64(x + -1.0) / t_0)); end return tmp end
function tmp_2 = code(x) t_0 = x + (4.0 * sqrt(x)); tmp = 0.0; if (x <= 0.29) tmp = -6.0 / (1.0 + t_0); else tmp = 6.0 * ((x + -1.0) / t_0); end tmp_2 = tmp; end
code[x_] := Block[{t$95$0 = N[(x + N[(4.0 * N[Sqrt[x], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[x, 0.29], N[(-6.0 / N[(1.0 + t$95$0), $MachinePrecision]), $MachinePrecision], N[(6.0 * N[(N[(x + -1.0), $MachinePrecision] / t$95$0), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := x + 4 \cdot \sqrt{x}\\
\mathbf{if}\;x \leq 0.29:\\
\;\;\;\;\frac{-6}{1 + t\_0}\\
\mathbf{else}:\\
\;\;\;\;6 \cdot \frac{x + -1}{t\_0}\\
\end{array}
\end{array}
if x < 0.28999999999999998Initial program 99.9%
+-commutativeN/A
associate-+r+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f6499.9%
Applied egg-rr99.9%
Taylor expanded in x around 0
Simplified99.3%
if 0.28999999999999998 < x Initial program 99.7%
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
associate-+l+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f6499.9%
Applied egg-rr99.9%
Taylor expanded in x around inf
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f6498.0%
Simplified98.0%
Final simplification98.5%
(FPCore (x) :precision binary64 (let* ((t_0 (+ x (* 4.0 (sqrt x))))) (if (<= x 1.0) (/ -6.0 (+ 1.0 t_0)) (* 6.0 (/ x t_0)))))
double code(double x) {
double t_0 = x + (4.0 * sqrt(x));
double tmp;
if (x <= 1.0) {
tmp = -6.0 / (1.0 + t_0);
} else {
tmp = 6.0 * (x / t_0);
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: t_0
real(8) :: tmp
t_0 = x + (4.0d0 * sqrt(x))
if (x <= 1.0d0) then
tmp = (-6.0d0) / (1.0d0 + t_0)
else
tmp = 6.0d0 * (x / t_0)
end if
code = tmp
end function
public static double code(double x) {
double t_0 = x + (4.0 * Math.sqrt(x));
double tmp;
if (x <= 1.0) {
tmp = -6.0 / (1.0 + t_0);
} else {
tmp = 6.0 * (x / t_0);
}
return tmp;
}
def code(x): t_0 = x + (4.0 * math.sqrt(x)) tmp = 0 if x <= 1.0: tmp = -6.0 / (1.0 + t_0) else: tmp = 6.0 * (x / t_0) return tmp
function code(x) t_0 = Float64(x + Float64(4.0 * sqrt(x))) tmp = 0.0 if (x <= 1.0) tmp = Float64(-6.0 / Float64(1.0 + t_0)); else tmp = Float64(6.0 * Float64(x / t_0)); end return tmp end
function tmp_2 = code(x) t_0 = x + (4.0 * sqrt(x)); tmp = 0.0; if (x <= 1.0) tmp = -6.0 / (1.0 + t_0); else tmp = 6.0 * (x / t_0); end tmp_2 = tmp; end
code[x_] := Block[{t$95$0 = N[(x + N[(4.0 * N[Sqrt[x], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[x, 1.0], N[(-6.0 / N[(1.0 + t$95$0), $MachinePrecision]), $MachinePrecision], N[(6.0 * N[(x / t$95$0), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := x + 4 \cdot \sqrt{x}\\
\mathbf{if}\;x \leq 1:\\
\;\;\;\;\frac{-6}{1 + t\_0}\\
\mathbf{else}:\\
\;\;\;\;6 \cdot \frac{x}{t\_0}\\
\end{array}
\end{array}
if x < 1Initial program 99.9%
+-commutativeN/A
associate-+r+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f6499.9%
Applied egg-rr99.9%
Taylor expanded in x around 0
Simplified99.3%
if 1 < x Initial program 99.7%
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
associate-+l+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f6499.9%
Applied egg-rr99.9%
Taylor expanded in x around inf
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f6498.0%
Simplified98.0%
Taylor expanded in x around inf
Simplified97.8%
Final simplification98.4%
(FPCore (x) :precision binary64 (let* ((t_0 (* 4.0 (sqrt x)))) (if (<= x 1.0) (/ -6.0 (+ t_0 (+ x 1.0))) (* 6.0 (/ x (+ x t_0))))))
double code(double x) {
double t_0 = 4.0 * sqrt(x);
double tmp;
if (x <= 1.0) {
tmp = -6.0 / (t_0 + (x + 1.0));
} else {
tmp = 6.0 * (x / (x + t_0));
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: t_0
real(8) :: tmp
t_0 = 4.0d0 * sqrt(x)
if (x <= 1.0d0) then
tmp = (-6.0d0) / (t_0 + (x + 1.0d0))
else
tmp = 6.0d0 * (x / (x + t_0))
end if
code = tmp
end function
public static double code(double x) {
double t_0 = 4.0 * Math.sqrt(x);
double tmp;
if (x <= 1.0) {
tmp = -6.0 / (t_0 + (x + 1.0));
} else {
tmp = 6.0 * (x / (x + t_0));
}
return tmp;
}
def code(x): t_0 = 4.0 * math.sqrt(x) tmp = 0 if x <= 1.0: tmp = -6.0 / (t_0 + (x + 1.0)) else: tmp = 6.0 * (x / (x + t_0)) return tmp
function code(x) t_0 = Float64(4.0 * sqrt(x)) tmp = 0.0 if (x <= 1.0) tmp = Float64(-6.0 / Float64(t_0 + Float64(x + 1.0))); else tmp = Float64(6.0 * Float64(x / Float64(x + t_0))); end return tmp end
function tmp_2 = code(x) t_0 = 4.0 * sqrt(x); tmp = 0.0; if (x <= 1.0) tmp = -6.0 / (t_0 + (x + 1.0)); else tmp = 6.0 * (x / (x + t_0)); end tmp_2 = tmp; end
code[x_] := Block[{t$95$0 = N[(4.0 * N[Sqrt[x], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[x, 1.0], N[(-6.0 / N[(t$95$0 + N[(x + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(6.0 * N[(x / N[(x + t$95$0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := 4 \cdot \sqrt{x}\\
\mathbf{if}\;x \leq 1:\\
\;\;\;\;\frac{-6}{t\_0 + \left(x + 1\right)}\\
\mathbf{else}:\\
\;\;\;\;6 \cdot \frac{x}{x + t\_0}\\
\end{array}
\end{array}
if x < 1Initial program 99.9%
Taylor expanded in x around 0
Simplified99.3%
if 1 < x Initial program 99.7%
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
associate-+l+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f6499.9%
Applied egg-rr99.9%
Taylor expanded in x around inf
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f6498.0%
Simplified98.0%
Taylor expanded in x around inf
Simplified97.8%
Final simplification98.4%
(FPCore (x) :precision binary64 (if (<= x 1.0) (* 6.0 (/ 1.0 (+ -1.0 (* (sqrt x) -4.0)))) (* 6.0 (/ x (+ x (* 4.0 (sqrt x)))))))
double code(double x) {
double tmp;
if (x <= 1.0) {
tmp = 6.0 * (1.0 / (-1.0 + (sqrt(x) * -4.0)));
} else {
tmp = 6.0 * (x / (x + (4.0 * sqrt(x))));
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: tmp
if (x <= 1.0d0) then
tmp = 6.0d0 * (1.0d0 / ((-1.0d0) + (sqrt(x) * (-4.0d0))))
else
tmp = 6.0d0 * (x / (x + (4.0d0 * sqrt(x))))
end if
code = tmp
end function
public static double code(double x) {
double tmp;
if (x <= 1.0) {
tmp = 6.0 * (1.0 / (-1.0 + (Math.sqrt(x) * -4.0)));
} else {
tmp = 6.0 * (x / (x + (4.0 * Math.sqrt(x))));
}
return tmp;
}
def code(x): tmp = 0 if x <= 1.0: tmp = 6.0 * (1.0 / (-1.0 + (math.sqrt(x) * -4.0))) else: tmp = 6.0 * (x / (x + (4.0 * math.sqrt(x)))) return tmp
function code(x) tmp = 0.0 if (x <= 1.0) tmp = Float64(6.0 * Float64(1.0 / Float64(-1.0 + Float64(sqrt(x) * -4.0)))); else tmp = Float64(6.0 * Float64(x / Float64(x + Float64(4.0 * sqrt(x))))); end return tmp end
function tmp_2 = code(x) tmp = 0.0; if (x <= 1.0) tmp = 6.0 * (1.0 / (-1.0 + (sqrt(x) * -4.0))); else tmp = 6.0 * (x / (x + (4.0 * sqrt(x)))); end tmp_2 = tmp; end
code[x_] := If[LessEqual[x, 1.0], N[(6.0 * N[(1.0 / N[(-1.0 + N[(N[Sqrt[x], $MachinePrecision] * -4.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(6.0 * N[(x / N[(x + N[(4.0 * N[Sqrt[x], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 1:\\
\;\;\;\;6 \cdot \frac{1}{-1 + \sqrt{x} \cdot -4}\\
\mathbf{else}:\\
\;\;\;\;6 \cdot \frac{x}{x + 4 \cdot \sqrt{x}}\\
\end{array}
\end{array}
if x < 1Initial program 99.9%
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
associate-+l+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f6499.9%
Applied egg-rr99.9%
Taylor expanded in x around 0
metadata-evalN/A
distribute-neg-fracN/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
distribute-neg-inN/A
metadata-evalN/A
distribute-rgt-neg-outN/A
mul-1-negN/A
associate-*r*N/A
metadata-evalN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f6499.3%
Simplified99.3%
if 1 < x Initial program 99.7%
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
associate-+l+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f6499.9%
Applied egg-rr99.9%
Taylor expanded in x around inf
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f6498.0%
Simplified98.0%
Taylor expanded in x around inf
Simplified97.8%
Final simplification98.4%
(FPCore (x) :precision binary64 (if (<= x 1.0) (* 6.0 (/ 1.0 (+ -1.0 (* (sqrt x) -4.0)))) (* -6.0 (/ 1.0 (+ -1.0 (/ -4.0 (sqrt x)))))))
double code(double x) {
double tmp;
if (x <= 1.0) {
tmp = 6.0 * (1.0 / (-1.0 + (sqrt(x) * -4.0)));
} else {
tmp = -6.0 * (1.0 / (-1.0 + (-4.0 / sqrt(x))));
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: tmp
if (x <= 1.0d0) then
tmp = 6.0d0 * (1.0d0 / ((-1.0d0) + (sqrt(x) * (-4.0d0))))
else
tmp = (-6.0d0) * (1.0d0 / ((-1.0d0) + ((-4.0d0) / sqrt(x))))
end if
code = tmp
end function
public static double code(double x) {
double tmp;
if (x <= 1.0) {
tmp = 6.0 * (1.0 / (-1.0 + (Math.sqrt(x) * -4.0)));
} else {
tmp = -6.0 * (1.0 / (-1.0 + (-4.0 / Math.sqrt(x))));
}
return tmp;
}
def code(x): tmp = 0 if x <= 1.0: tmp = 6.0 * (1.0 / (-1.0 + (math.sqrt(x) * -4.0))) else: tmp = -6.0 * (1.0 / (-1.0 + (-4.0 / math.sqrt(x)))) return tmp
function code(x) tmp = 0.0 if (x <= 1.0) tmp = Float64(6.0 * Float64(1.0 / Float64(-1.0 + Float64(sqrt(x) * -4.0)))); else tmp = Float64(-6.0 * Float64(1.0 / Float64(-1.0 + Float64(-4.0 / sqrt(x))))); end return tmp end
function tmp_2 = code(x) tmp = 0.0; if (x <= 1.0) tmp = 6.0 * (1.0 / (-1.0 + (sqrt(x) * -4.0))); else tmp = -6.0 * (1.0 / (-1.0 + (-4.0 / sqrt(x)))); end tmp_2 = tmp; end
code[x_] := If[LessEqual[x, 1.0], N[(6.0 * N[(1.0 / N[(-1.0 + N[(N[Sqrt[x], $MachinePrecision] * -4.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(-6.0 * N[(1.0 / N[(-1.0 + N[(-4.0 / N[Sqrt[x], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 1:\\
\;\;\;\;6 \cdot \frac{1}{-1 + \sqrt{x} \cdot -4}\\
\mathbf{else}:\\
\;\;\;\;-6 \cdot \frac{1}{-1 + \frac{-4}{\sqrt{x}}}\\
\end{array}
\end{array}
if x < 1Initial program 99.9%
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
associate-+l+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f6499.9%
Applied egg-rr99.9%
Taylor expanded in x around 0
metadata-evalN/A
distribute-neg-fracN/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
distribute-neg-inN/A
metadata-evalN/A
distribute-rgt-neg-outN/A
mul-1-negN/A
associate-*r*N/A
metadata-evalN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f6499.3%
Simplified99.3%
if 1 < x Initial program 99.7%
Taylor expanded in x around inf
+-commutativeN/A
remove-double-negN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
rem-square-sqrtN/A
unpow2N/A
*-commutativeN/A
metadata-evalN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
distribute-neg-frac2N/A
distribute-neg-fracN/A
metadata-evalN/A
/-lowering-/.f64N/A
Simplified97.8%
clear-numN/A
associate-/r/N/A
flip3-+N/A
clear-numN/A
*-lowering-*.f64N/A
Applied egg-rr97.8%
Final simplification98.4%
(FPCore (x) :precision binary64 (if (<= x 1.0) (* 6.0 (/ 1.0 (+ -1.0 (* (sqrt x) -4.0)))) (/ 6.0 (+ 1.0 (/ 4.0 (sqrt x))))))
double code(double x) {
double tmp;
if (x <= 1.0) {
tmp = 6.0 * (1.0 / (-1.0 + (sqrt(x) * -4.0)));
} else {
tmp = 6.0 / (1.0 + (4.0 / sqrt(x)));
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: tmp
if (x <= 1.0d0) then
tmp = 6.0d0 * (1.0d0 / ((-1.0d0) + (sqrt(x) * (-4.0d0))))
else
tmp = 6.0d0 / (1.0d0 + (4.0d0 / sqrt(x)))
end if
code = tmp
end function
public static double code(double x) {
double tmp;
if (x <= 1.0) {
tmp = 6.0 * (1.0 / (-1.0 + (Math.sqrt(x) * -4.0)));
} else {
tmp = 6.0 / (1.0 + (4.0 / Math.sqrt(x)));
}
return tmp;
}
def code(x): tmp = 0 if x <= 1.0: tmp = 6.0 * (1.0 / (-1.0 + (math.sqrt(x) * -4.0))) else: tmp = 6.0 / (1.0 + (4.0 / math.sqrt(x))) return tmp
function code(x) tmp = 0.0 if (x <= 1.0) tmp = Float64(6.0 * Float64(1.0 / Float64(-1.0 + Float64(sqrt(x) * -4.0)))); else tmp = Float64(6.0 / Float64(1.0 + Float64(4.0 / sqrt(x)))); end return tmp end
function tmp_2 = code(x) tmp = 0.0; if (x <= 1.0) tmp = 6.0 * (1.0 / (-1.0 + (sqrt(x) * -4.0))); else tmp = 6.0 / (1.0 + (4.0 / sqrt(x))); end tmp_2 = tmp; end
code[x_] := If[LessEqual[x, 1.0], N[(6.0 * N[(1.0 / N[(-1.0 + N[(N[Sqrt[x], $MachinePrecision] * -4.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(6.0 / N[(1.0 + N[(4.0 / N[Sqrt[x], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 1:\\
\;\;\;\;6 \cdot \frac{1}{-1 + \sqrt{x} \cdot -4}\\
\mathbf{else}:\\
\;\;\;\;\frac{6}{1 + \frac{4}{\sqrt{x}}}\\
\end{array}
\end{array}
if x < 1Initial program 99.9%
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
associate-+l+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f6499.9%
Applied egg-rr99.9%
Taylor expanded in x around 0
metadata-evalN/A
distribute-neg-fracN/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
distribute-neg-inN/A
metadata-evalN/A
distribute-rgt-neg-outN/A
mul-1-negN/A
associate-*r*N/A
metadata-evalN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f6499.3%
Simplified99.3%
if 1 < x Initial program 99.7%
Taylor expanded in x around inf
+-commutativeN/A
remove-double-negN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
rem-square-sqrtN/A
unpow2N/A
*-commutativeN/A
metadata-evalN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
distribute-neg-frac2N/A
distribute-neg-fracN/A
metadata-evalN/A
/-lowering-/.f64N/A
Simplified97.8%
frac-2negN/A
metadata-evalN/A
/-lowering-/.f64N/A
distribute-neg-inN/A
metadata-evalN/A
mul-1-negN/A
+-lowering-+.f64N/A
*-commutativeN/A
associate-*r*N/A
metadata-evalN/A
sqrt-divN/A
metadata-evalN/A
un-div-invN/A
/-lowering-/.f64N/A
sqrt-lowering-sqrt.f6497.8%
Applied egg-rr97.8%
Final simplification98.4%
(FPCore (x) :precision binary64 (* 6.0 (/ (+ x -1.0) (+ x (+ 1.0 (* 4.0 (sqrt x)))))))
double code(double x) {
return 6.0 * ((x + -1.0) / (x + (1.0 + (4.0 * sqrt(x)))));
}
real(8) function code(x)
real(8), intent (in) :: x
code = 6.0d0 * ((x + (-1.0d0)) / (x + (1.0d0 + (4.0d0 * sqrt(x)))))
end function
public static double code(double x) {
return 6.0 * ((x + -1.0) / (x + (1.0 + (4.0 * Math.sqrt(x)))));
}
def code(x): return 6.0 * ((x + -1.0) / (x + (1.0 + (4.0 * math.sqrt(x)))))
function code(x) return Float64(6.0 * Float64(Float64(x + -1.0) / Float64(x + Float64(1.0 + Float64(4.0 * sqrt(x)))))) end
function tmp = code(x) tmp = 6.0 * ((x + -1.0) / (x + (1.0 + (4.0 * sqrt(x))))); end
code[x_] := N[(6.0 * N[(N[(x + -1.0), $MachinePrecision] / N[(x + N[(1.0 + N[(4.0 * N[Sqrt[x], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
6 \cdot \frac{x + -1}{x + \left(1 + 4 \cdot \sqrt{x}\right)}
\end{array}
Initial program 99.8%
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-evalN/A
associate-+l+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f6499.9%
Applied egg-rr99.9%
Final simplification99.9%
(FPCore (x) :precision binary64 (if (<= x 1.0) (/ 1.0 (+ -0.16666666666666666 (* (sqrt x) -0.6666666666666666))) (/ 6.0 (+ 1.0 (/ 4.0 (sqrt x))))))
double code(double x) {
double tmp;
if (x <= 1.0) {
tmp = 1.0 / (-0.16666666666666666 + (sqrt(x) * -0.6666666666666666));
} else {
tmp = 6.0 / (1.0 + (4.0 / sqrt(x)));
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: tmp
if (x <= 1.0d0) then
tmp = 1.0d0 / ((-0.16666666666666666d0) + (sqrt(x) * (-0.6666666666666666d0)))
else
tmp = 6.0d0 / (1.0d0 + (4.0d0 / sqrt(x)))
end if
code = tmp
end function
public static double code(double x) {
double tmp;
if (x <= 1.0) {
tmp = 1.0 / (-0.16666666666666666 + (Math.sqrt(x) * -0.6666666666666666));
} else {
tmp = 6.0 / (1.0 + (4.0 / Math.sqrt(x)));
}
return tmp;
}
def code(x): tmp = 0 if x <= 1.0: tmp = 1.0 / (-0.16666666666666666 + (math.sqrt(x) * -0.6666666666666666)) else: tmp = 6.0 / (1.0 + (4.0 / math.sqrt(x))) return tmp
function code(x) tmp = 0.0 if (x <= 1.0) tmp = Float64(1.0 / Float64(-0.16666666666666666 + Float64(sqrt(x) * -0.6666666666666666))); else tmp = Float64(6.0 / Float64(1.0 + Float64(4.0 / sqrt(x)))); end return tmp end
function tmp_2 = code(x) tmp = 0.0; if (x <= 1.0) tmp = 1.0 / (-0.16666666666666666 + (sqrt(x) * -0.6666666666666666)); else tmp = 6.0 / (1.0 + (4.0 / sqrt(x))); end tmp_2 = tmp; end
code[x_] := If[LessEqual[x, 1.0], N[(1.0 / N[(-0.16666666666666666 + N[(N[Sqrt[x], $MachinePrecision] * -0.6666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(6.0 / N[(1.0 + N[(4.0 / N[Sqrt[x], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 1:\\
\;\;\;\;\frac{1}{-0.16666666666666666 + \sqrt{x} \cdot -0.6666666666666666}\\
\mathbf{else}:\\
\;\;\;\;\frac{6}{1 + \frac{4}{\sqrt{x}}}\\
\end{array}
\end{array}
if x < 1Initial program 99.9%
Taylor expanded in x around 0
metadata-evalN/A
distribute-neg-fracN/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
distribute-neg-inN/A
metadata-evalN/A
+-lowering-+.f64N/A
*-commutativeN/A
distribute-rgt-neg-inN/A
metadata-evalN/A
metadata-evalN/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f64N/A
metadata-eval99.2%
Simplified99.2%
clear-numN/A
/-lowering-/.f64N/A
frac-2negN/A
distribute-neg-inN/A
metadata-evalN/A
mul-1-negN/A
*-commutativeN/A
associate-*r*N/A
metadata-evalN/A
metadata-evalN/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f6499.3%
Applied egg-rr99.3%
Taylor expanded in x around 0
distribute-rgt-inN/A
metadata-evalN/A
+-lowering-+.f64N/A
*-commutativeN/A
associate-*l*N/A
metadata-evalN/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f6499.2%
Simplified99.2%
if 1 < x Initial program 99.7%
Taylor expanded in x around inf
+-commutativeN/A
remove-double-negN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
rem-square-sqrtN/A
unpow2N/A
*-commutativeN/A
metadata-evalN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
distribute-neg-frac2N/A
distribute-neg-fracN/A
metadata-evalN/A
/-lowering-/.f64N/A
Simplified97.8%
frac-2negN/A
metadata-evalN/A
/-lowering-/.f64N/A
distribute-neg-inN/A
metadata-evalN/A
mul-1-negN/A
+-lowering-+.f64N/A
*-commutativeN/A
associate-*r*N/A
metadata-evalN/A
sqrt-divN/A
metadata-evalN/A
un-div-invN/A
/-lowering-/.f64N/A
sqrt-lowering-sqrt.f6497.8%
Applied egg-rr97.8%
(FPCore (x) :precision binary64 (if (<= x 0.83) (/ 1.0 (+ -0.16666666666666666 (* (sqrt x) -0.6666666666666666))) (/ 6.0 (/ x (+ x -1.0)))))
double code(double x) {
double tmp;
if (x <= 0.83) {
tmp = 1.0 / (-0.16666666666666666 + (sqrt(x) * -0.6666666666666666));
} else {
tmp = 6.0 / (x / (x + -1.0));
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: tmp
if (x <= 0.83d0) then
tmp = 1.0d0 / ((-0.16666666666666666d0) + (sqrt(x) * (-0.6666666666666666d0)))
else
tmp = 6.0d0 / (x / (x + (-1.0d0)))
end if
code = tmp
end function
public static double code(double x) {
double tmp;
if (x <= 0.83) {
tmp = 1.0 / (-0.16666666666666666 + (Math.sqrt(x) * -0.6666666666666666));
} else {
tmp = 6.0 / (x / (x + -1.0));
}
return tmp;
}
def code(x): tmp = 0 if x <= 0.83: tmp = 1.0 / (-0.16666666666666666 + (math.sqrt(x) * -0.6666666666666666)) else: tmp = 6.0 / (x / (x + -1.0)) return tmp
function code(x) tmp = 0.0 if (x <= 0.83) tmp = Float64(1.0 / Float64(-0.16666666666666666 + Float64(sqrt(x) * -0.6666666666666666))); else tmp = Float64(6.0 / Float64(x / Float64(x + -1.0))); end return tmp end
function tmp_2 = code(x) tmp = 0.0; if (x <= 0.83) tmp = 1.0 / (-0.16666666666666666 + (sqrt(x) * -0.6666666666666666)); else tmp = 6.0 / (x / (x + -1.0)); end tmp_2 = tmp; end
code[x_] := If[LessEqual[x, 0.83], N[(1.0 / N[(-0.16666666666666666 + N[(N[Sqrt[x], $MachinePrecision] * -0.6666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(6.0 / N[(x / N[(x + -1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 0.83:\\
\;\;\;\;\frac{1}{-0.16666666666666666 + \sqrt{x} \cdot -0.6666666666666666}\\
\mathbf{else}:\\
\;\;\;\;\frac{6}{\frac{x}{x + -1}}\\
\end{array}
\end{array}
if x < 0.82999999999999996Initial program 99.9%
Taylor expanded in x around 0
metadata-evalN/A
distribute-neg-fracN/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
distribute-neg-inN/A
metadata-evalN/A
+-lowering-+.f64N/A
*-commutativeN/A
distribute-rgt-neg-inN/A
metadata-evalN/A
metadata-evalN/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f64N/A
metadata-eval99.2%
Simplified99.2%
clear-numN/A
/-lowering-/.f64N/A
frac-2negN/A
distribute-neg-inN/A
metadata-evalN/A
mul-1-negN/A
*-commutativeN/A
associate-*r*N/A
metadata-evalN/A
metadata-evalN/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f6499.3%
Applied egg-rr99.3%
Taylor expanded in x around 0
distribute-rgt-inN/A
metadata-evalN/A
+-lowering-+.f64N/A
*-commutativeN/A
associate-*l*N/A
metadata-evalN/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f6499.2%
Simplified99.2%
if 0.82999999999999996 < x Initial program 99.7%
associate-/l*N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
associate-+l+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-eval99.9%
Applied egg-rr99.9%
flip-+N/A
/-lowering-/.f64N/A
metadata-evalN/A
--lowering--.f64N/A
*-commutativeN/A
*-commutativeN/A
swap-sqrN/A
metadata-evalN/A
metadata-evalN/A
rem-square-sqrtN/A
*-lowering-*.f64N/A
metadata-evalN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f6499.9%
Applied egg-rr99.9%
Taylor expanded in x around inf
Simplified94.6%
(FPCore (x) :precision binary64 (if (<= x 0.061) (+ -6.0 (* -6.0 (* (sqrt x) -4.0))) (/ 6.0 (/ x (+ x -1.0)))))
double code(double x) {
double tmp;
if (x <= 0.061) {
tmp = -6.0 + (-6.0 * (sqrt(x) * -4.0));
} else {
tmp = 6.0 / (x / (x + -1.0));
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: tmp
if (x <= 0.061d0) then
tmp = (-6.0d0) + ((-6.0d0) * (sqrt(x) * (-4.0d0)))
else
tmp = 6.0d0 / (x / (x + (-1.0d0)))
end if
code = tmp
end function
public static double code(double x) {
double tmp;
if (x <= 0.061) {
tmp = -6.0 + (-6.0 * (Math.sqrt(x) * -4.0));
} else {
tmp = 6.0 / (x / (x + -1.0));
}
return tmp;
}
def code(x): tmp = 0 if x <= 0.061: tmp = -6.0 + (-6.0 * (math.sqrt(x) * -4.0)) else: tmp = 6.0 / (x / (x + -1.0)) return tmp
function code(x) tmp = 0.0 if (x <= 0.061) tmp = Float64(-6.0 + Float64(-6.0 * Float64(sqrt(x) * -4.0))); else tmp = Float64(6.0 / Float64(x / Float64(x + -1.0))); end return tmp end
function tmp_2 = code(x) tmp = 0.0; if (x <= 0.061) tmp = -6.0 + (-6.0 * (sqrt(x) * -4.0)); else tmp = 6.0 / (x / (x + -1.0)); end tmp_2 = tmp; end
code[x_] := If[LessEqual[x, 0.061], N[(-6.0 + N[(-6.0 * N[(N[Sqrt[x], $MachinePrecision] * -4.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(6.0 / N[(x / N[(x + -1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 0.061:\\
\;\;\;\;-6 + -6 \cdot \left(\sqrt{x} \cdot -4\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{6}{\frac{x}{x + -1}}\\
\end{array}
\end{array}
if x < 0.060999999999999999Initial program 99.9%
associate-/l*N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
associate-+l+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-eval99.9%
Applied egg-rr99.9%
flip-+N/A
/-lowering-/.f64N/A
metadata-evalN/A
--lowering--.f64N/A
*-commutativeN/A
*-commutativeN/A
swap-sqrN/A
metadata-evalN/A
metadata-evalN/A
rem-square-sqrtN/A
*-lowering-*.f64N/A
metadata-evalN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f6499.9%
Applied egg-rr99.9%
Taylor expanded in x around 0
cancel-sign-sub-invN/A
metadata-evalN/A
distribute-rgt-inN/A
metadata-evalN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f6499.1%
Simplified99.1%
if 0.060999999999999999 < x Initial program 99.7%
associate-/l*N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
associate-+l+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-eval99.9%
Applied egg-rr99.9%
flip-+N/A
/-lowering-/.f64N/A
metadata-evalN/A
--lowering--.f64N/A
*-commutativeN/A
*-commutativeN/A
swap-sqrN/A
metadata-evalN/A
metadata-evalN/A
rem-square-sqrtN/A
*-lowering-*.f64N/A
metadata-evalN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f6499.9%
Applied egg-rr99.9%
Taylor expanded in x around inf
Simplified94.6%
Final simplification96.5%
(FPCore (x) :precision binary64 (/ 6.0 (/ x (+ x -1.0))))
double code(double x) {
return 6.0 / (x / (x + -1.0));
}
real(8) function code(x)
real(8), intent (in) :: x
code = 6.0d0 / (x / (x + (-1.0d0)))
end function
public static double code(double x) {
return 6.0 / (x / (x + -1.0));
}
def code(x): return 6.0 / (x / (x + -1.0))
function code(x) return Float64(6.0 / Float64(x / Float64(x + -1.0))) end
function tmp = code(x) tmp = 6.0 / (x / (x + -1.0)); end
code[x_] := N[(6.0 / N[(x / N[(x + -1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{6}{\frac{x}{x + -1}}
\end{array}
Initial program 99.8%
associate-/l*N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
associate-+l+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-eval99.9%
Applied egg-rr99.9%
flip-+N/A
/-lowering-/.f64N/A
metadata-evalN/A
--lowering--.f64N/A
*-commutativeN/A
*-commutativeN/A
swap-sqrN/A
metadata-evalN/A
metadata-evalN/A
rem-square-sqrtN/A
*-lowering-*.f64N/A
metadata-evalN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f6499.9%
Applied egg-rr99.9%
Taylor expanded in x around inf
Simplified56.9%
(FPCore (x) :precision binary64 6.0)
double code(double x) {
return 6.0;
}
real(8) function code(x)
real(8), intent (in) :: x
code = 6.0d0
end function
public static double code(double x) {
return 6.0;
}
def code(x): return 6.0
function code(x) return 6.0 end
function tmp = code(x) tmp = 6.0; end
code[x_] := 6.0
\begin{array}{l}
\\
6
\end{array}
Initial program 99.8%
associate-/l*N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
associate-+l+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-eval99.9%
Applied egg-rr99.9%
flip-+N/A
/-lowering-/.f64N/A
metadata-evalN/A
--lowering--.f64N/A
*-commutativeN/A
*-commutativeN/A
swap-sqrN/A
metadata-evalN/A
metadata-evalN/A
rem-square-sqrtN/A
*-lowering-*.f64N/A
metadata-evalN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f6499.9%
Applied egg-rr99.9%
Taylor expanded in x around inf
Simplified55.4%
(FPCore (x) :precision binary64 (/ 6.0 (/ (+ (+ x 1.0) (* 4.0 (sqrt x))) (- x 1.0))))
double code(double x) {
return 6.0 / (((x + 1.0) + (4.0 * sqrt(x))) / (x - 1.0));
}
real(8) function code(x)
real(8), intent (in) :: x
code = 6.0d0 / (((x + 1.0d0) + (4.0d0 * sqrt(x))) / (x - 1.0d0))
end function
public static double code(double x) {
return 6.0 / (((x + 1.0) + (4.0 * Math.sqrt(x))) / (x - 1.0));
}
def code(x): return 6.0 / (((x + 1.0) + (4.0 * math.sqrt(x))) / (x - 1.0))
function code(x) return Float64(6.0 / Float64(Float64(Float64(x + 1.0) + Float64(4.0 * sqrt(x))) / Float64(x - 1.0))) end
function tmp = code(x) tmp = 6.0 / (((x + 1.0) + (4.0 * sqrt(x))) / (x - 1.0)); end
code[x_] := N[(6.0 / N[(N[(N[(x + 1.0), $MachinePrecision] + N[(4.0 * N[Sqrt[x], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(x - 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{6}{\frac{\left(x + 1\right) + 4 \cdot \sqrt{x}}{x - 1}}
\end{array}
herbie shell --seed 2024160
(FPCore (x)
:name "Data.Approximate.Numerics:blog from approximate-0.2.2.1"
:precision binary64
:alt
(! :herbie-platform default (/ 6 (/ (+ (+ x 1) (* 4 (sqrt x))) (- x 1))))
(/ (* 6.0 (- x 1.0)) (+ (+ x 1.0) (* 4.0 (sqrt x)))))