
(FPCore (x eps) :precision binary64 (- (tan (+ x eps)) (tan x)))
double code(double x, double eps) {
return tan((x + eps)) - tan(x);
}
real(8) function code(x, eps)
real(8), intent (in) :: x
real(8), intent (in) :: eps
code = tan((x + eps)) - tan(x)
end function
public static double code(double x, double eps) {
return Math.tan((x + eps)) - Math.tan(x);
}
def code(x, eps): return math.tan((x + eps)) - math.tan(x)
function code(x, eps) return Float64(tan(Float64(x + eps)) - tan(x)) end
function tmp = code(x, eps) tmp = tan((x + eps)) - tan(x); end
code[x_, eps_] := N[(N[Tan[N[(x + eps), $MachinePrecision]], $MachinePrecision] - N[Tan[x], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\tan \left(x + \varepsilon\right) - \tan x
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 16 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x eps) :precision binary64 (- (tan (+ x eps)) (tan x)))
double code(double x, double eps) {
return tan((x + eps)) - tan(x);
}
real(8) function code(x, eps)
real(8), intent (in) :: x
real(8), intent (in) :: eps
code = tan((x + eps)) - tan(x)
end function
public static double code(double x, double eps) {
return Math.tan((x + eps)) - Math.tan(x);
}
def code(x, eps): return math.tan((x + eps)) - math.tan(x)
function code(x, eps) return Float64(tan(Float64(x + eps)) - tan(x)) end
function tmp = code(x, eps) tmp = tan((x + eps)) - tan(x); end
code[x_, eps_] := N[(N[Tan[N[(x + eps), $MachinePrecision]], $MachinePrecision] - N[Tan[x], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\tan \left(x + \varepsilon\right) - \tan x
\end{array}
(FPCore (x eps)
:precision binary64
(let* ((t_0 (pow (sin x) 2.0))
(t_1 (pow (cos x) 2.0))
(t_2
(-
(/ (+ t_0 (/ (pow (sin x) 4.0) t_1)) t_1)
(+
0.16666666666666666
(+
(+ -0.5 (/ (* t_0 -0.5) t_1))
(/ (* t_0 0.16666666666666666) t_1)))))
(t_3 (/ (+ (sin x) (/ (pow (sin x) 3.0) t_1)) (cos x))))
(*
eps
(+
(+
(*
eps
(+
t_3
(*
eps
(+
t_2
(*
eps
(+ (* t_2 (/ (sin x) (cos x))) (* t_3 0.3333333333333333)))))))
(/ t_0 t_1))
1.0))))
double code(double x, double eps) {
double t_0 = pow(sin(x), 2.0);
double t_1 = pow(cos(x), 2.0);
double t_2 = ((t_0 + (pow(sin(x), 4.0) / t_1)) / t_1) - (0.16666666666666666 + ((-0.5 + ((t_0 * -0.5) / t_1)) + ((t_0 * 0.16666666666666666) / t_1)));
double t_3 = (sin(x) + (pow(sin(x), 3.0) / t_1)) / cos(x);
return eps * (((eps * (t_3 + (eps * (t_2 + (eps * ((t_2 * (sin(x) / cos(x))) + (t_3 * 0.3333333333333333))))))) + (t_0 / t_1)) + 1.0);
}
real(8) function code(x, eps)
real(8), intent (in) :: x
real(8), intent (in) :: eps
real(8) :: t_0
real(8) :: t_1
real(8) :: t_2
real(8) :: t_3
t_0 = sin(x) ** 2.0d0
t_1 = cos(x) ** 2.0d0
t_2 = ((t_0 + ((sin(x) ** 4.0d0) / t_1)) / t_1) - (0.16666666666666666d0 + (((-0.5d0) + ((t_0 * (-0.5d0)) / t_1)) + ((t_0 * 0.16666666666666666d0) / t_1)))
t_3 = (sin(x) + ((sin(x) ** 3.0d0) / t_1)) / cos(x)
code = eps * (((eps * (t_3 + (eps * (t_2 + (eps * ((t_2 * (sin(x) / cos(x))) + (t_3 * 0.3333333333333333d0))))))) + (t_0 / t_1)) + 1.0d0)
end function
public static double code(double x, double eps) {
double t_0 = Math.pow(Math.sin(x), 2.0);
double t_1 = Math.pow(Math.cos(x), 2.0);
double t_2 = ((t_0 + (Math.pow(Math.sin(x), 4.0) / t_1)) / t_1) - (0.16666666666666666 + ((-0.5 + ((t_0 * -0.5) / t_1)) + ((t_0 * 0.16666666666666666) / t_1)));
double t_3 = (Math.sin(x) + (Math.pow(Math.sin(x), 3.0) / t_1)) / Math.cos(x);
return eps * (((eps * (t_3 + (eps * (t_2 + (eps * ((t_2 * (Math.sin(x) / Math.cos(x))) + (t_3 * 0.3333333333333333))))))) + (t_0 / t_1)) + 1.0);
}
def code(x, eps): t_0 = math.pow(math.sin(x), 2.0) t_1 = math.pow(math.cos(x), 2.0) t_2 = ((t_0 + (math.pow(math.sin(x), 4.0) / t_1)) / t_1) - (0.16666666666666666 + ((-0.5 + ((t_0 * -0.5) / t_1)) + ((t_0 * 0.16666666666666666) / t_1))) t_3 = (math.sin(x) + (math.pow(math.sin(x), 3.0) / t_1)) / math.cos(x) return eps * (((eps * (t_3 + (eps * (t_2 + (eps * ((t_2 * (math.sin(x) / math.cos(x))) + (t_3 * 0.3333333333333333))))))) + (t_0 / t_1)) + 1.0)
function code(x, eps) t_0 = sin(x) ^ 2.0 t_1 = cos(x) ^ 2.0 t_2 = Float64(Float64(Float64(t_0 + Float64((sin(x) ^ 4.0) / t_1)) / t_1) - Float64(0.16666666666666666 + Float64(Float64(-0.5 + Float64(Float64(t_0 * -0.5) / t_1)) + Float64(Float64(t_0 * 0.16666666666666666) / t_1)))) t_3 = Float64(Float64(sin(x) + Float64((sin(x) ^ 3.0) / t_1)) / cos(x)) return Float64(eps * Float64(Float64(Float64(eps * Float64(t_3 + Float64(eps * Float64(t_2 + Float64(eps * Float64(Float64(t_2 * Float64(sin(x) / cos(x))) + Float64(t_3 * 0.3333333333333333))))))) + Float64(t_0 / t_1)) + 1.0)) end
function tmp = code(x, eps) t_0 = sin(x) ^ 2.0; t_1 = cos(x) ^ 2.0; t_2 = ((t_0 + ((sin(x) ^ 4.0) / t_1)) / t_1) - (0.16666666666666666 + ((-0.5 + ((t_0 * -0.5) / t_1)) + ((t_0 * 0.16666666666666666) / t_1))); t_3 = (sin(x) + ((sin(x) ^ 3.0) / t_1)) / cos(x); tmp = eps * (((eps * (t_3 + (eps * (t_2 + (eps * ((t_2 * (sin(x) / cos(x))) + (t_3 * 0.3333333333333333))))))) + (t_0 / t_1)) + 1.0); end
code[x_, eps_] := Block[{t$95$0 = N[Power[N[Sin[x], $MachinePrecision], 2.0], $MachinePrecision]}, Block[{t$95$1 = N[Power[N[Cos[x], $MachinePrecision], 2.0], $MachinePrecision]}, Block[{t$95$2 = N[(N[(N[(t$95$0 + N[(N[Power[N[Sin[x], $MachinePrecision], 4.0], $MachinePrecision] / t$95$1), $MachinePrecision]), $MachinePrecision] / t$95$1), $MachinePrecision] - N[(0.16666666666666666 + N[(N[(-0.5 + N[(N[(t$95$0 * -0.5), $MachinePrecision] / t$95$1), $MachinePrecision]), $MachinePrecision] + N[(N[(t$95$0 * 0.16666666666666666), $MachinePrecision] / t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = N[(N[(N[Sin[x], $MachinePrecision] + N[(N[Power[N[Sin[x], $MachinePrecision], 3.0], $MachinePrecision] / t$95$1), $MachinePrecision]), $MachinePrecision] / N[Cos[x], $MachinePrecision]), $MachinePrecision]}, N[(eps * N[(N[(N[(eps * N[(t$95$3 + N[(eps * N[(t$95$2 + N[(eps * N[(N[(t$95$2 * N[(N[Sin[x], $MachinePrecision] / N[Cos[x], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(t$95$3 * 0.3333333333333333), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(t$95$0 / t$95$1), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := {\sin x}^{2}\\
t_1 := {\cos x}^{2}\\
t_2 := \frac{t\_0 + \frac{{\sin x}^{4}}{t\_1}}{t\_1} - \left(0.16666666666666666 + \left(\left(-0.5 + \frac{t\_0 \cdot -0.5}{t\_1}\right) + \frac{t\_0 \cdot 0.16666666666666666}{t\_1}\right)\right)\\
t_3 := \frac{\sin x + \frac{{\sin x}^{3}}{t\_1}}{\cos x}\\
\varepsilon \cdot \left(\left(\varepsilon \cdot \left(t\_3 + \varepsilon \cdot \left(t\_2 + \varepsilon \cdot \left(t\_2 \cdot \frac{\sin x}{\cos x} + t\_3 \cdot 0.3333333333333333\right)\right)\right) + \frac{t\_0}{t\_1}\right) + 1\right)
\end{array}
\end{array}
Initial program 61.4%
Taylor expanded in eps around 0
Simplified99.9%
Final simplification99.9%
(FPCore (x eps)
:precision binary64
(let* ((t_0 (pow (sin x) 2.0)) (t_1 (pow (cos x) 2.0)))
(*
eps
(+
(+
(/ t_0 t_1)
(*
eps
(+
(/ (+ (sin x) (/ (pow (sin x) 3.0) t_1)) (cos x))
(*
eps
(+
(-
(/ (+ t_0 (/ (pow (sin x) 4.0) t_1)) t_1)
(+
0.16666666666666666
(+
(+ -0.5 (/ (* t_0 -0.5) t_1))
(/ (* t_0 0.16666666666666666) t_1))))
(*
x
(+
(* eps 0.6666666666666666)
(* 1.8888888888888888 (* x (* eps x))))))))))
1.0))))
double code(double x, double eps) {
double t_0 = pow(sin(x), 2.0);
double t_1 = pow(cos(x), 2.0);
return eps * (((t_0 / t_1) + (eps * (((sin(x) + (pow(sin(x), 3.0) / t_1)) / cos(x)) + (eps * ((((t_0 + (pow(sin(x), 4.0) / t_1)) / t_1) - (0.16666666666666666 + ((-0.5 + ((t_0 * -0.5) / t_1)) + ((t_0 * 0.16666666666666666) / t_1)))) + (x * ((eps * 0.6666666666666666) + (1.8888888888888888 * (x * (eps * x)))))))))) + 1.0);
}
real(8) function code(x, eps)
real(8), intent (in) :: x
real(8), intent (in) :: eps
real(8) :: t_0
real(8) :: t_1
t_0 = sin(x) ** 2.0d0
t_1 = cos(x) ** 2.0d0
code = eps * (((t_0 / t_1) + (eps * (((sin(x) + ((sin(x) ** 3.0d0) / t_1)) / cos(x)) + (eps * ((((t_0 + ((sin(x) ** 4.0d0) / t_1)) / t_1) - (0.16666666666666666d0 + (((-0.5d0) + ((t_0 * (-0.5d0)) / t_1)) + ((t_0 * 0.16666666666666666d0) / t_1)))) + (x * ((eps * 0.6666666666666666d0) + (1.8888888888888888d0 * (x * (eps * x)))))))))) + 1.0d0)
end function
public static double code(double x, double eps) {
double t_0 = Math.pow(Math.sin(x), 2.0);
double t_1 = Math.pow(Math.cos(x), 2.0);
return eps * (((t_0 / t_1) + (eps * (((Math.sin(x) + (Math.pow(Math.sin(x), 3.0) / t_1)) / Math.cos(x)) + (eps * ((((t_0 + (Math.pow(Math.sin(x), 4.0) / t_1)) / t_1) - (0.16666666666666666 + ((-0.5 + ((t_0 * -0.5) / t_1)) + ((t_0 * 0.16666666666666666) / t_1)))) + (x * ((eps * 0.6666666666666666) + (1.8888888888888888 * (x * (eps * x)))))))))) + 1.0);
}
def code(x, eps): t_0 = math.pow(math.sin(x), 2.0) t_1 = math.pow(math.cos(x), 2.0) return eps * (((t_0 / t_1) + (eps * (((math.sin(x) + (math.pow(math.sin(x), 3.0) / t_1)) / math.cos(x)) + (eps * ((((t_0 + (math.pow(math.sin(x), 4.0) / t_1)) / t_1) - (0.16666666666666666 + ((-0.5 + ((t_0 * -0.5) / t_1)) + ((t_0 * 0.16666666666666666) / t_1)))) + (x * ((eps * 0.6666666666666666) + (1.8888888888888888 * (x * (eps * x)))))))))) + 1.0)
function code(x, eps) t_0 = sin(x) ^ 2.0 t_1 = cos(x) ^ 2.0 return Float64(eps * Float64(Float64(Float64(t_0 / t_1) + Float64(eps * Float64(Float64(Float64(sin(x) + Float64((sin(x) ^ 3.0) / t_1)) / cos(x)) + Float64(eps * Float64(Float64(Float64(Float64(t_0 + Float64((sin(x) ^ 4.0) / t_1)) / t_1) - Float64(0.16666666666666666 + Float64(Float64(-0.5 + Float64(Float64(t_0 * -0.5) / t_1)) + Float64(Float64(t_0 * 0.16666666666666666) / t_1)))) + Float64(x * Float64(Float64(eps * 0.6666666666666666) + Float64(1.8888888888888888 * Float64(x * Float64(eps * x)))))))))) + 1.0)) end
function tmp = code(x, eps) t_0 = sin(x) ^ 2.0; t_1 = cos(x) ^ 2.0; tmp = eps * (((t_0 / t_1) + (eps * (((sin(x) + ((sin(x) ^ 3.0) / t_1)) / cos(x)) + (eps * ((((t_0 + ((sin(x) ^ 4.0) / t_1)) / t_1) - (0.16666666666666666 + ((-0.5 + ((t_0 * -0.5) / t_1)) + ((t_0 * 0.16666666666666666) / t_1)))) + (x * ((eps * 0.6666666666666666) + (1.8888888888888888 * (x * (eps * x)))))))))) + 1.0); end
code[x_, eps_] := Block[{t$95$0 = N[Power[N[Sin[x], $MachinePrecision], 2.0], $MachinePrecision]}, Block[{t$95$1 = N[Power[N[Cos[x], $MachinePrecision], 2.0], $MachinePrecision]}, N[(eps * N[(N[(N[(t$95$0 / t$95$1), $MachinePrecision] + N[(eps * N[(N[(N[(N[Sin[x], $MachinePrecision] + N[(N[Power[N[Sin[x], $MachinePrecision], 3.0], $MachinePrecision] / t$95$1), $MachinePrecision]), $MachinePrecision] / N[Cos[x], $MachinePrecision]), $MachinePrecision] + N[(eps * N[(N[(N[(N[(t$95$0 + N[(N[Power[N[Sin[x], $MachinePrecision], 4.0], $MachinePrecision] / t$95$1), $MachinePrecision]), $MachinePrecision] / t$95$1), $MachinePrecision] - N[(0.16666666666666666 + N[(N[(-0.5 + N[(N[(t$95$0 * -0.5), $MachinePrecision] / t$95$1), $MachinePrecision]), $MachinePrecision] + N[(N[(t$95$0 * 0.16666666666666666), $MachinePrecision] / t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(x * N[(N[(eps * 0.6666666666666666), $MachinePrecision] + N[(1.8888888888888888 * N[(x * N[(eps * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := {\sin x}^{2}\\
t_1 := {\cos x}^{2}\\
\varepsilon \cdot \left(\left(\frac{t\_0}{t\_1} + \varepsilon \cdot \left(\frac{\sin x + \frac{{\sin x}^{3}}{t\_1}}{\cos x} + \varepsilon \cdot \left(\left(\frac{t\_0 + \frac{{\sin x}^{4}}{t\_1}}{t\_1} - \left(0.16666666666666666 + \left(\left(-0.5 + \frac{t\_0 \cdot -0.5}{t\_1}\right) + \frac{t\_0 \cdot 0.16666666666666666}{t\_1}\right)\right)\right) + x \cdot \left(\varepsilon \cdot 0.6666666666666666 + 1.8888888888888888 \cdot \left(x \cdot \left(\varepsilon \cdot x\right)\right)\right)\right)\right)\right) + 1\right)
\end{array}
\end{array}
Initial program 61.4%
Taylor expanded in eps around 0
Simplified99.9%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6499.9%
Simplified99.9%
Final simplification99.9%
(FPCore (x eps)
:precision binary64
(let* ((t_0 (pow (tan x) 2.0))
(t_1 (* 0.5 (cos (* x 2.0))))
(t_2 (+ 0.5 t_1)))
(*
eps
(+
(+
t_0
(*
eps
(+
(/ (+ (sin x) (* (sin x) t_0)) (cos x))
(*
eps
(+
(/ (+ (- 0.5 t_1) (/ (pow (sin x) 4.0) t_2)) t_2)
(-
(*
x
(+
(* eps 0.6666666666666666)
(* x (* 1.8888888888888888 (* eps x)))))
(+
(*
(* x x)
(+
(*
(* x x)
(+
(*
x
(* x (+ (* (* x x) -0.0656084656084656) -0.1259259259259259)))
-0.2222222222222222))
-0.3333333333333333))
-0.3333333333333333)))))))
1.0))))
double code(double x, double eps) {
double t_0 = pow(tan(x), 2.0);
double t_1 = 0.5 * cos((x * 2.0));
double t_2 = 0.5 + t_1;
return eps * ((t_0 + (eps * (((sin(x) + (sin(x) * t_0)) / cos(x)) + (eps * ((((0.5 - t_1) + (pow(sin(x), 4.0) / t_2)) / t_2) + ((x * ((eps * 0.6666666666666666) + (x * (1.8888888888888888 * (eps * x))))) - (((x * x) * (((x * x) * ((x * (x * (((x * x) * -0.0656084656084656) + -0.1259259259259259))) + -0.2222222222222222)) + -0.3333333333333333)) + -0.3333333333333333))))))) + 1.0);
}
real(8) function code(x, eps)
real(8), intent (in) :: x
real(8), intent (in) :: eps
real(8) :: t_0
real(8) :: t_1
real(8) :: t_2
t_0 = tan(x) ** 2.0d0
t_1 = 0.5d0 * cos((x * 2.0d0))
t_2 = 0.5d0 + t_1
code = eps * ((t_0 + (eps * (((sin(x) + (sin(x) * t_0)) / cos(x)) + (eps * ((((0.5d0 - t_1) + ((sin(x) ** 4.0d0) / t_2)) / t_2) + ((x * ((eps * 0.6666666666666666d0) + (x * (1.8888888888888888d0 * (eps * x))))) - (((x * x) * (((x * x) * ((x * (x * (((x * x) * (-0.0656084656084656d0)) + (-0.1259259259259259d0)))) + (-0.2222222222222222d0))) + (-0.3333333333333333d0))) + (-0.3333333333333333d0)))))))) + 1.0d0)
end function
public static double code(double x, double eps) {
double t_0 = Math.pow(Math.tan(x), 2.0);
double t_1 = 0.5 * Math.cos((x * 2.0));
double t_2 = 0.5 + t_1;
return eps * ((t_0 + (eps * (((Math.sin(x) + (Math.sin(x) * t_0)) / Math.cos(x)) + (eps * ((((0.5 - t_1) + (Math.pow(Math.sin(x), 4.0) / t_2)) / t_2) + ((x * ((eps * 0.6666666666666666) + (x * (1.8888888888888888 * (eps * x))))) - (((x * x) * (((x * x) * ((x * (x * (((x * x) * -0.0656084656084656) + -0.1259259259259259))) + -0.2222222222222222)) + -0.3333333333333333)) + -0.3333333333333333))))))) + 1.0);
}
def code(x, eps): t_0 = math.pow(math.tan(x), 2.0) t_1 = 0.5 * math.cos((x * 2.0)) t_2 = 0.5 + t_1 return eps * ((t_0 + (eps * (((math.sin(x) + (math.sin(x) * t_0)) / math.cos(x)) + (eps * ((((0.5 - t_1) + (math.pow(math.sin(x), 4.0) / t_2)) / t_2) + ((x * ((eps * 0.6666666666666666) + (x * (1.8888888888888888 * (eps * x))))) - (((x * x) * (((x * x) * ((x * (x * (((x * x) * -0.0656084656084656) + -0.1259259259259259))) + -0.2222222222222222)) + -0.3333333333333333)) + -0.3333333333333333))))))) + 1.0)
function code(x, eps) t_0 = tan(x) ^ 2.0 t_1 = Float64(0.5 * cos(Float64(x * 2.0))) t_2 = Float64(0.5 + t_1) return Float64(eps * Float64(Float64(t_0 + Float64(eps * Float64(Float64(Float64(sin(x) + Float64(sin(x) * t_0)) / cos(x)) + Float64(eps * Float64(Float64(Float64(Float64(0.5 - t_1) + Float64((sin(x) ^ 4.0) / t_2)) / t_2) + Float64(Float64(x * Float64(Float64(eps * 0.6666666666666666) + Float64(x * Float64(1.8888888888888888 * Float64(eps * x))))) - Float64(Float64(Float64(x * x) * Float64(Float64(Float64(x * x) * Float64(Float64(x * Float64(x * Float64(Float64(Float64(x * x) * -0.0656084656084656) + -0.1259259259259259))) + -0.2222222222222222)) + -0.3333333333333333)) + -0.3333333333333333))))))) + 1.0)) end
function tmp = code(x, eps) t_0 = tan(x) ^ 2.0; t_1 = 0.5 * cos((x * 2.0)); t_2 = 0.5 + t_1; tmp = eps * ((t_0 + (eps * (((sin(x) + (sin(x) * t_0)) / cos(x)) + (eps * ((((0.5 - t_1) + ((sin(x) ^ 4.0) / t_2)) / t_2) + ((x * ((eps * 0.6666666666666666) + (x * (1.8888888888888888 * (eps * x))))) - (((x * x) * (((x * x) * ((x * (x * (((x * x) * -0.0656084656084656) + -0.1259259259259259))) + -0.2222222222222222)) + -0.3333333333333333)) + -0.3333333333333333))))))) + 1.0); end
code[x_, eps_] := Block[{t$95$0 = N[Power[N[Tan[x], $MachinePrecision], 2.0], $MachinePrecision]}, Block[{t$95$1 = N[(0.5 * N[Cos[N[(x * 2.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(0.5 + t$95$1), $MachinePrecision]}, N[(eps * N[(N[(t$95$0 + N[(eps * N[(N[(N[(N[Sin[x], $MachinePrecision] + N[(N[Sin[x], $MachinePrecision] * t$95$0), $MachinePrecision]), $MachinePrecision] / N[Cos[x], $MachinePrecision]), $MachinePrecision] + N[(eps * N[(N[(N[(N[(0.5 - t$95$1), $MachinePrecision] + N[(N[Power[N[Sin[x], $MachinePrecision], 4.0], $MachinePrecision] / t$95$2), $MachinePrecision]), $MachinePrecision] / t$95$2), $MachinePrecision] + N[(N[(x * N[(N[(eps * 0.6666666666666666), $MachinePrecision] + N[(x * N[(1.8888888888888888 * N[(eps * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(N[(N[(x * x), $MachinePrecision] * N[(N[(N[(x * x), $MachinePrecision] * N[(N[(x * N[(x * N[(N[(N[(x * x), $MachinePrecision] * -0.0656084656084656), $MachinePrecision] + -0.1259259259259259), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + -0.2222222222222222), $MachinePrecision]), $MachinePrecision] + -0.3333333333333333), $MachinePrecision]), $MachinePrecision] + -0.3333333333333333), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := {\tan x}^{2}\\
t_1 := 0.5 \cdot \cos \left(x \cdot 2\right)\\
t_2 := 0.5 + t\_1\\
\varepsilon \cdot \left(\left(t\_0 + \varepsilon \cdot \left(\frac{\sin x + \sin x \cdot t\_0}{\cos x} + \varepsilon \cdot \left(\frac{\left(0.5 - t\_1\right) + \frac{{\sin x}^{4}}{t\_2}}{t\_2} + \left(x \cdot \left(\varepsilon \cdot 0.6666666666666666 + x \cdot \left(1.8888888888888888 \cdot \left(\varepsilon \cdot x\right)\right)\right) - \left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot \left(x \cdot \left(\left(x \cdot x\right) \cdot -0.0656084656084656 + -0.1259259259259259\right)\right) + -0.2222222222222222\right) + -0.3333333333333333\right) + -0.3333333333333333\right)\right)\right)\right)\right) + 1\right)
\end{array}
\end{array}
Initial program 61.4%
Taylor expanded in eps around 0
Simplified99.9%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6499.9%
Simplified99.9%
Applied egg-rr99.9%
Taylor expanded in x around 0
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
Simplified99.9%
Final simplification99.9%
(FPCore (x eps)
:precision binary64
(let* ((t_0 (pow (tan x) 2.0))
(t_1 (* 0.5 (cos (* x 2.0))))
(t_2 (+ 0.5 t_1)))
(*
eps
(+
(+
t_0
(*
eps
(+
(/ (+ (sin x) (* (sin x) t_0)) (cos x))
(*
eps
(+
(/ (+ (- 0.5 t_1) (/ (pow (sin x) 4.0) t_2)) t_2)
(-
(*
x
(+
(* eps 0.6666666666666666)
(* x (* 1.8888888888888888 (* eps x)))))
(+
(*
(* x x)
(+
(*
x
(* x (+ -0.2222222222222222 (* (* x x) -0.1259259259259259))))
-0.3333333333333333))
-0.3333333333333333)))))))
1.0))))
double code(double x, double eps) {
double t_0 = pow(tan(x), 2.0);
double t_1 = 0.5 * cos((x * 2.0));
double t_2 = 0.5 + t_1;
return eps * ((t_0 + (eps * (((sin(x) + (sin(x) * t_0)) / cos(x)) + (eps * ((((0.5 - t_1) + (pow(sin(x), 4.0) / t_2)) / t_2) + ((x * ((eps * 0.6666666666666666) + (x * (1.8888888888888888 * (eps * x))))) - (((x * x) * ((x * (x * (-0.2222222222222222 + ((x * x) * -0.1259259259259259)))) + -0.3333333333333333)) + -0.3333333333333333))))))) + 1.0);
}
real(8) function code(x, eps)
real(8), intent (in) :: x
real(8), intent (in) :: eps
real(8) :: t_0
real(8) :: t_1
real(8) :: t_2
t_0 = tan(x) ** 2.0d0
t_1 = 0.5d0 * cos((x * 2.0d0))
t_2 = 0.5d0 + t_1
code = eps * ((t_0 + (eps * (((sin(x) + (sin(x) * t_0)) / cos(x)) + (eps * ((((0.5d0 - t_1) + ((sin(x) ** 4.0d0) / t_2)) / t_2) + ((x * ((eps * 0.6666666666666666d0) + (x * (1.8888888888888888d0 * (eps * x))))) - (((x * x) * ((x * (x * ((-0.2222222222222222d0) + ((x * x) * (-0.1259259259259259d0))))) + (-0.3333333333333333d0))) + (-0.3333333333333333d0)))))))) + 1.0d0)
end function
public static double code(double x, double eps) {
double t_0 = Math.pow(Math.tan(x), 2.0);
double t_1 = 0.5 * Math.cos((x * 2.0));
double t_2 = 0.5 + t_1;
return eps * ((t_0 + (eps * (((Math.sin(x) + (Math.sin(x) * t_0)) / Math.cos(x)) + (eps * ((((0.5 - t_1) + (Math.pow(Math.sin(x), 4.0) / t_2)) / t_2) + ((x * ((eps * 0.6666666666666666) + (x * (1.8888888888888888 * (eps * x))))) - (((x * x) * ((x * (x * (-0.2222222222222222 + ((x * x) * -0.1259259259259259)))) + -0.3333333333333333)) + -0.3333333333333333))))))) + 1.0);
}
def code(x, eps): t_0 = math.pow(math.tan(x), 2.0) t_1 = 0.5 * math.cos((x * 2.0)) t_2 = 0.5 + t_1 return eps * ((t_0 + (eps * (((math.sin(x) + (math.sin(x) * t_0)) / math.cos(x)) + (eps * ((((0.5 - t_1) + (math.pow(math.sin(x), 4.0) / t_2)) / t_2) + ((x * ((eps * 0.6666666666666666) + (x * (1.8888888888888888 * (eps * x))))) - (((x * x) * ((x * (x * (-0.2222222222222222 + ((x * x) * -0.1259259259259259)))) + -0.3333333333333333)) + -0.3333333333333333))))))) + 1.0)
function code(x, eps) t_0 = tan(x) ^ 2.0 t_1 = Float64(0.5 * cos(Float64(x * 2.0))) t_2 = Float64(0.5 + t_1) return Float64(eps * Float64(Float64(t_0 + Float64(eps * Float64(Float64(Float64(sin(x) + Float64(sin(x) * t_0)) / cos(x)) + Float64(eps * Float64(Float64(Float64(Float64(0.5 - t_1) + Float64((sin(x) ^ 4.0) / t_2)) / t_2) + Float64(Float64(x * Float64(Float64(eps * 0.6666666666666666) + Float64(x * Float64(1.8888888888888888 * Float64(eps * x))))) - Float64(Float64(Float64(x * x) * Float64(Float64(x * Float64(x * Float64(-0.2222222222222222 + Float64(Float64(x * x) * -0.1259259259259259)))) + -0.3333333333333333)) + -0.3333333333333333))))))) + 1.0)) end
function tmp = code(x, eps) t_0 = tan(x) ^ 2.0; t_1 = 0.5 * cos((x * 2.0)); t_2 = 0.5 + t_1; tmp = eps * ((t_0 + (eps * (((sin(x) + (sin(x) * t_0)) / cos(x)) + (eps * ((((0.5 - t_1) + ((sin(x) ^ 4.0) / t_2)) / t_2) + ((x * ((eps * 0.6666666666666666) + (x * (1.8888888888888888 * (eps * x))))) - (((x * x) * ((x * (x * (-0.2222222222222222 + ((x * x) * -0.1259259259259259)))) + -0.3333333333333333)) + -0.3333333333333333))))))) + 1.0); end
code[x_, eps_] := Block[{t$95$0 = N[Power[N[Tan[x], $MachinePrecision], 2.0], $MachinePrecision]}, Block[{t$95$1 = N[(0.5 * N[Cos[N[(x * 2.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(0.5 + t$95$1), $MachinePrecision]}, N[(eps * N[(N[(t$95$0 + N[(eps * N[(N[(N[(N[Sin[x], $MachinePrecision] + N[(N[Sin[x], $MachinePrecision] * t$95$0), $MachinePrecision]), $MachinePrecision] / N[Cos[x], $MachinePrecision]), $MachinePrecision] + N[(eps * N[(N[(N[(N[(0.5 - t$95$1), $MachinePrecision] + N[(N[Power[N[Sin[x], $MachinePrecision], 4.0], $MachinePrecision] / t$95$2), $MachinePrecision]), $MachinePrecision] / t$95$2), $MachinePrecision] + N[(N[(x * N[(N[(eps * 0.6666666666666666), $MachinePrecision] + N[(x * N[(1.8888888888888888 * N[(eps * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(N[(N[(x * x), $MachinePrecision] * N[(N[(x * N[(x * N[(-0.2222222222222222 + N[(N[(x * x), $MachinePrecision] * -0.1259259259259259), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + -0.3333333333333333), $MachinePrecision]), $MachinePrecision] + -0.3333333333333333), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := {\tan x}^{2}\\
t_1 := 0.5 \cdot \cos \left(x \cdot 2\right)\\
t_2 := 0.5 + t\_1\\
\varepsilon \cdot \left(\left(t\_0 + \varepsilon \cdot \left(\frac{\sin x + \sin x \cdot t\_0}{\cos x} + \varepsilon \cdot \left(\frac{\left(0.5 - t\_1\right) + \frac{{\sin x}^{4}}{t\_2}}{t\_2} + \left(x \cdot \left(\varepsilon \cdot 0.6666666666666666 + x \cdot \left(1.8888888888888888 \cdot \left(\varepsilon \cdot x\right)\right)\right) - \left(\left(x \cdot x\right) \cdot \left(x \cdot \left(x \cdot \left(-0.2222222222222222 + \left(x \cdot x\right) \cdot -0.1259259259259259\right)\right) + -0.3333333333333333\right) + -0.3333333333333333\right)\right)\right)\right)\right) + 1\right)
\end{array}
\end{array}
Initial program 61.4%
Taylor expanded in eps around 0
Simplified99.9%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6499.9%
Simplified99.9%
Applied egg-rr99.9%
Taylor expanded in x around 0
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6499.9%
Simplified99.9%
Final simplification99.9%
(FPCore (x eps)
:precision binary64
(let* ((t_0 (pow (cos x) 2.0)))
(*
eps
(+
(+
(/ (pow (sin x) 2.0) t_0)
(*
eps
(+
(/ (+ (sin x) (/ (pow (sin x) 3.0) t_0)) (cos x))
(* eps 0.3333333333333333))))
1.0))))
double code(double x, double eps) {
double t_0 = pow(cos(x), 2.0);
return eps * (((pow(sin(x), 2.0) / t_0) + (eps * (((sin(x) + (pow(sin(x), 3.0) / t_0)) / cos(x)) + (eps * 0.3333333333333333)))) + 1.0);
}
real(8) function code(x, eps)
real(8), intent (in) :: x
real(8), intent (in) :: eps
real(8) :: t_0
t_0 = cos(x) ** 2.0d0
code = eps * ((((sin(x) ** 2.0d0) / t_0) + (eps * (((sin(x) + ((sin(x) ** 3.0d0) / t_0)) / cos(x)) + (eps * 0.3333333333333333d0)))) + 1.0d0)
end function
public static double code(double x, double eps) {
double t_0 = Math.pow(Math.cos(x), 2.0);
return eps * (((Math.pow(Math.sin(x), 2.0) / t_0) + (eps * (((Math.sin(x) + (Math.pow(Math.sin(x), 3.0) / t_0)) / Math.cos(x)) + (eps * 0.3333333333333333)))) + 1.0);
}
def code(x, eps): t_0 = math.pow(math.cos(x), 2.0) return eps * (((math.pow(math.sin(x), 2.0) / t_0) + (eps * (((math.sin(x) + (math.pow(math.sin(x), 3.0) / t_0)) / math.cos(x)) + (eps * 0.3333333333333333)))) + 1.0)
function code(x, eps) t_0 = cos(x) ^ 2.0 return Float64(eps * Float64(Float64(Float64((sin(x) ^ 2.0) / t_0) + Float64(eps * Float64(Float64(Float64(sin(x) + Float64((sin(x) ^ 3.0) / t_0)) / cos(x)) + Float64(eps * 0.3333333333333333)))) + 1.0)) end
function tmp = code(x, eps) t_0 = cos(x) ^ 2.0; tmp = eps * ((((sin(x) ^ 2.0) / t_0) + (eps * (((sin(x) + ((sin(x) ^ 3.0) / t_0)) / cos(x)) + (eps * 0.3333333333333333)))) + 1.0); end
code[x_, eps_] := Block[{t$95$0 = N[Power[N[Cos[x], $MachinePrecision], 2.0], $MachinePrecision]}, N[(eps * N[(N[(N[(N[Power[N[Sin[x], $MachinePrecision], 2.0], $MachinePrecision] / t$95$0), $MachinePrecision] + N[(eps * N[(N[(N[(N[Sin[x], $MachinePrecision] + N[(N[Power[N[Sin[x], $MachinePrecision], 3.0], $MachinePrecision] / t$95$0), $MachinePrecision]), $MachinePrecision] / N[Cos[x], $MachinePrecision]), $MachinePrecision] + N[(eps * 0.3333333333333333), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := {\cos x}^{2}\\
\varepsilon \cdot \left(\left(\frac{{\sin x}^{2}}{t\_0} + \varepsilon \cdot \left(\frac{\sin x + \frac{{\sin x}^{3}}{t\_0}}{\cos x} + \varepsilon \cdot 0.3333333333333333\right)\right) + 1\right)
\end{array}
\end{array}
Initial program 61.4%
Taylor expanded in eps around 0
Simplified99.9%
Taylor expanded in x around 0
Simplified99.7%
Final simplification99.7%
(FPCore (x eps) :precision binary64 (let* ((t_0 (pow (tan x) 2.0))) (+ eps (* eps (+ t_0 (/ eps (/ (cos x) (+ (sin x) (* (sin x) t_0)))))))))
double code(double x, double eps) {
double t_0 = pow(tan(x), 2.0);
return eps + (eps * (t_0 + (eps / (cos(x) / (sin(x) + (sin(x) * t_0))))));
}
real(8) function code(x, eps)
real(8), intent (in) :: x
real(8), intent (in) :: eps
real(8) :: t_0
t_0 = tan(x) ** 2.0d0
code = eps + (eps * (t_0 + (eps / (cos(x) / (sin(x) + (sin(x) * t_0))))))
end function
public static double code(double x, double eps) {
double t_0 = Math.pow(Math.tan(x), 2.0);
return eps + (eps * (t_0 + (eps / (Math.cos(x) / (Math.sin(x) + (Math.sin(x) * t_0))))));
}
def code(x, eps): t_0 = math.pow(math.tan(x), 2.0) return eps + (eps * (t_0 + (eps / (math.cos(x) / (math.sin(x) + (math.sin(x) * t_0))))))
function code(x, eps) t_0 = tan(x) ^ 2.0 return Float64(eps + Float64(eps * Float64(t_0 + Float64(eps / Float64(cos(x) / Float64(sin(x) + Float64(sin(x) * t_0))))))) end
function tmp = code(x, eps) t_0 = tan(x) ^ 2.0; tmp = eps + (eps * (t_0 + (eps / (cos(x) / (sin(x) + (sin(x) * t_0)))))); end
code[x_, eps_] := Block[{t$95$0 = N[Power[N[Tan[x], $MachinePrecision], 2.0], $MachinePrecision]}, N[(eps + N[(eps * N[(t$95$0 + N[(eps / N[(N[Cos[x], $MachinePrecision] / N[(N[Sin[x], $MachinePrecision] + N[(N[Sin[x], $MachinePrecision] * t$95$0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := {\tan x}^{2}\\
\varepsilon + \varepsilon \cdot \left(t\_0 + \frac{\varepsilon}{\frac{\cos x}{\sin x + \sin x \cdot t\_0}}\right)
\end{array}
\end{array}
Initial program 61.4%
Taylor expanded in eps around 0
*-lowering-*.f64N/A
+-commutativeN/A
associate--l+N/A
+-commutativeN/A
+-lowering-+.f64N/A
Simplified99.7%
associate-+l+N/A
distribute-rgt-inN/A
*-lft-identityN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
Applied egg-rr99.7%
Final simplification99.7%
(FPCore (x eps) :precision binary64 (* eps (+ (pow (tan x) 2.0) (+ (* eps x) 1.0))))
double code(double x, double eps) {
return eps * (pow(tan(x), 2.0) + ((eps * x) + 1.0));
}
real(8) function code(x, eps)
real(8), intent (in) :: x
real(8), intent (in) :: eps
code = eps * ((tan(x) ** 2.0d0) + ((eps * x) + 1.0d0))
end function
public static double code(double x, double eps) {
return eps * (Math.pow(Math.tan(x), 2.0) + ((eps * x) + 1.0));
}
def code(x, eps): return eps * (math.pow(math.tan(x), 2.0) + ((eps * x) + 1.0))
function code(x, eps) return Float64(eps * Float64((tan(x) ^ 2.0) + Float64(Float64(eps * x) + 1.0))) end
function tmp = code(x, eps) tmp = eps * ((tan(x) ^ 2.0) + ((eps * x) + 1.0)); end
code[x_, eps_] := N[(eps * N[(N[Power[N[Tan[x], $MachinePrecision], 2.0], $MachinePrecision] + N[(N[(eps * x), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\varepsilon \cdot \left({\tan x}^{2} + \left(\varepsilon \cdot x + 1\right)\right)
\end{array}
Initial program 61.4%
Taylor expanded in eps around 0
*-lowering-*.f64N/A
+-commutativeN/A
associate--l+N/A
+-commutativeN/A
+-lowering-+.f64N/A
Simplified99.7%
Taylor expanded in x around 0
Simplified99.4%
*-commutativeN/A
*-lowering-*.f64N/A
+-commutativeN/A
associate-+l+N/A
unpow2N/A
unpow2N/A
frac-timesN/A
tan-quotN/A
tan-quotN/A
unpow2N/A
+-lowering-+.f64N/A
pow-lowering-pow.f64N/A
tan-lowering-tan.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6499.4%
Applied egg-rr99.4%
Final simplification99.4%
(FPCore (x eps)
:precision binary64
(*
eps
(+
(+
(*
(* x x)
(+
(*
(* x x)
(+
0.6666666666666666
(* (* x x) (+ 0.37777777777777777 (* (* x x) 0.19682539682539682)))))
1.0))
1.0)
(* eps x))))
double code(double x, double eps) {
return eps * ((((x * x) * (((x * x) * (0.6666666666666666 + ((x * x) * (0.37777777777777777 + ((x * x) * 0.19682539682539682))))) + 1.0)) + 1.0) + (eps * x));
}
real(8) function code(x, eps)
real(8), intent (in) :: x
real(8), intent (in) :: eps
code = eps * ((((x * x) * (((x * x) * (0.6666666666666666d0 + ((x * x) * (0.37777777777777777d0 + ((x * x) * 0.19682539682539682d0))))) + 1.0d0)) + 1.0d0) + (eps * x))
end function
public static double code(double x, double eps) {
return eps * ((((x * x) * (((x * x) * (0.6666666666666666 + ((x * x) * (0.37777777777777777 + ((x * x) * 0.19682539682539682))))) + 1.0)) + 1.0) + (eps * x));
}
def code(x, eps): return eps * ((((x * x) * (((x * x) * (0.6666666666666666 + ((x * x) * (0.37777777777777777 + ((x * x) * 0.19682539682539682))))) + 1.0)) + 1.0) + (eps * x))
function code(x, eps) return Float64(eps * Float64(Float64(Float64(Float64(x * x) * Float64(Float64(Float64(x * x) * Float64(0.6666666666666666 + Float64(Float64(x * x) * Float64(0.37777777777777777 + Float64(Float64(x * x) * 0.19682539682539682))))) + 1.0)) + 1.0) + Float64(eps * x))) end
function tmp = code(x, eps) tmp = eps * ((((x * x) * (((x * x) * (0.6666666666666666 + ((x * x) * (0.37777777777777777 + ((x * x) * 0.19682539682539682))))) + 1.0)) + 1.0) + (eps * x)); end
code[x_, eps_] := N[(eps * N[(N[(N[(N[(x * x), $MachinePrecision] * N[(N[(N[(x * x), $MachinePrecision] * N[(0.6666666666666666 + N[(N[(x * x), $MachinePrecision] * N[(0.37777777777777777 + N[(N[(x * x), $MachinePrecision] * 0.19682539682539682), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision] + N[(eps * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\varepsilon \cdot \left(\left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(0.6666666666666666 + \left(x \cdot x\right) \cdot \left(0.37777777777777777 + \left(x \cdot x\right) \cdot 0.19682539682539682\right)\right) + 1\right) + 1\right) + \varepsilon \cdot x\right)
\end{array}
Initial program 61.4%
Taylor expanded in eps around 0
*-lowering-*.f64N/A
+-commutativeN/A
associate--l+N/A
+-commutativeN/A
+-lowering-+.f64N/A
Simplified99.7%
Taylor expanded in x around 0
Simplified99.4%
Taylor expanded in x around 0
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6498.9%
Simplified98.9%
Final simplification98.9%
(FPCore (x eps)
:precision binary64
(*
eps
(+
(* eps x)
(+
(*
(* x x)
(+
(* (* x x) (+ 0.6666666666666666 (* (* x x) 0.37777777777777777)))
1.0))
1.0))))
double code(double x, double eps) {
return eps * ((eps * x) + (((x * x) * (((x * x) * (0.6666666666666666 + ((x * x) * 0.37777777777777777))) + 1.0)) + 1.0));
}
real(8) function code(x, eps)
real(8), intent (in) :: x
real(8), intent (in) :: eps
code = eps * ((eps * x) + (((x * x) * (((x * x) * (0.6666666666666666d0 + ((x * x) * 0.37777777777777777d0))) + 1.0d0)) + 1.0d0))
end function
public static double code(double x, double eps) {
return eps * ((eps * x) + (((x * x) * (((x * x) * (0.6666666666666666 + ((x * x) * 0.37777777777777777))) + 1.0)) + 1.0));
}
def code(x, eps): return eps * ((eps * x) + (((x * x) * (((x * x) * (0.6666666666666666 + ((x * x) * 0.37777777777777777))) + 1.0)) + 1.0))
function code(x, eps) return Float64(eps * Float64(Float64(eps * x) + Float64(Float64(Float64(x * x) * Float64(Float64(Float64(x * x) * Float64(0.6666666666666666 + Float64(Float64(x * x) * 0.37777777777777777))) + 1.0)) + 1.0))) end
function tmp = code(x, eps) tmp = eps * ((eps * x) + (((x * x) * (((x * x) * (0.6666666666666666 + ((x * x) * 0.37777777777777777))) + 1.0)) + 1.0)); end
code[x_, eps_] := N[(eps * N[(N[(eps * x), $MachinePrecision] + N[(N[(N[(x * x), $MachinePrecision] * N[(N[(N[(x * x), $MachinePrecision] * N[(0.6666666666666666 + N[(N[(x * x), $MachinePrecision] * 0.37777777777777777), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\varepsilon \cdot \left(\varepsilon \cdot x + \left(\left(x \cdot x\right) \cdot \left(\left(x \cdot x\right) \cdot \left(0.6666666666666666 + \left(x \cdot x\right) \cdot 0.37777777777777777\right) + 1\right) + 1\right)\right)
\end{array}
Initial program 61.4%
Taylor expanded in eps around 0
*-lowering-*.f64N/A
+-commutativeN/A
associate--l+N/A
+-commutativeN/A
+-lowering-+.f64N/A
Simplified99.7%
Taylor expanded in x around 0
Simplified99.4%
Taylor expanded in x around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6498.9%
Simplified98.9%
Final simplification98.9%
(FPCore (x eps) :precision binary64 (* eps (+ (* x (+ eps (* x (+ (* x (* x 0.6666666666666666)) 1.0)))) 1.0)))
double code(double x, double eps) {
return eps * ((x * (eps + (x * ((x * (x * 0.6666666666666666)) + 1.0)))) + 1.0);
}
real(8) function code(x, eps)
real(8), intent (in) :: x
real(8), intent (in) :: eps
code = eps * ((x * (eps + (x * ((x * (x * 0.6666666666666666d0)) + 1.0d0)))) + 1.0d0)
end function
public static double code(double x, double eps) {
return eps * ((x * (eps + (x * ((x * (x * 0.6666666666666666)) + 1.0)))) + 1.0);
}
def code(x, eps): return eps * ((x * (eps + (x * ((x * (x * 0.6666666666666666)) + 1.0)))) + 1.0)
function code(x, eps) return Float64(eps * Float64(Float64(x * Float64(eps + Float64(x * Float64(Float64(x * Float64(x * 0.6666666666666666)) + 1.0)))) + 1.0)) end
function tmp = code(x, eps) tmp = eps * ((x * (eps + (x * ((x * (x * 0.6666666666666666)) + 1.0)))) + 1.0); end
code[x_, eps_] := N[(eps * N[(N[(x * N[(eps + N[(x * N[(N[(x * N[(x * 0.6666666666666666), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\varepsilon \cdot \left(x \cdot \left(\varepsilon + x \cdot \left(x \cdot \left(x \cdot 0.6666666666666666\right) + 1\right)\right) + 1\right)
\end{array}
Initial program 61.4%
Taylor expanded in eps around 0
*-lowering-*.f64N/A
+-commutativeN/A
associate--l+N/A
+-commutativeN/A
+-lowering-+.f64N/A
Simplified99.7%
Taylor expanded in x around 0
Simplified99.4%
Taylor expanded in x around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f6498.9%
Simplified98.9%
Final simplification98.9%
(FPCore (x eps) :precision binary64 (+ eps (* x (* eps (+ eps x)))))
double code(double x, double eps) {
return eps + (x * (eps * (eps + x)));
}
real(8) function code(x, eps)
real(8), intent (in) :: x
real(8), intent (in) :: eps
code = eps + (x * (eps * (eps + x)))
end function
public static double code(double x, double eps) {
return eps + (x * (eps * (eps + x)));
}
def code(x, eps): return eps + (x * (eps * (eps + x)))
function code(x, eps) return Float64(eps + Float64(x * Float64(eps * Float64(eps + x)))) end
function tmp = code(x, eps) tmp = eps + (x * (eps * (eps + x))); end
code[x_, eps_] := N[(eps + N[(x * N[(eps * N[(eps + x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\varepsilon + x \cdot \left(\varepsilon \cdot \left(\varepsilon + x\right)\right)
\end{array}
Initial program 61.4%
Taylor expanded in eps around 0
*-lowering-*.f64N/A
+-commutativeN/A
associate--l+N/A
+-commutativeN/A
+-lowering-+.f64N/A
Simplified99.7%
Taylor expanded in x around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
unpow2N/A
distribute-lft-outN/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f6498.8%
Simplified98.8%
Final simplification98.8%
(FPCore (x eps) :precision binary64 (+ eps (* x (* eps x))))
double code(double x, double eps) {
return eps + (x * (eps * x));
}
real(8) function code(x, eps)
real(8), intent (in) :: x
real(8), intent (in) :: eps
code = eps + (x * (eps * x))
end function
public static double code(double x, double eps) {
return eps + (x * (eps * x));
}
def code(x, eps): return eps + (x * (eps * x))
function code(x, eps) return Float64(eps + Float64(x * Float64(eps * x))) end
function tmp = code(x, eps) tmp = eps + (x * (eps * x)); end
code[x_, eps_] := N[(eps + N[(x * N[(eps * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\varepsilon + x \cdot \left(\varepsilon \cdot x\right)
\end{array}
Initial program 61.4%
Taylor expanded in eps around 0
*-lowering-*.f64N/A
+-commutativeN/A
associate--l+N/A
+-commutativeN/A
+-lowering-+.f64N/A
Simplified99.7%
Taylor expanded in x around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
distribute-rgt1-inN/A
+-commutativeN/A
*-lowering-*.f64N/A
Simplified98.8%
Taylor expanded in eps around 0
*-commutativeN/A
*-lowering-*.f6498.8%
Simplified98.8%
Final simplification98.8%
(FPCore (x eps) :precision binary64 (+ eps (* x (* eps eps))))
double code(double x, double eps) {
return eps + (x * (eps * eps));
}
real(8) function code(x, eps)
real(8), intent (in) :: x
real(8), intent (in) :: eps
code = eps + (x * (eps * eps))
end function
public static double code(double x, double eps) {
return eps + (x * (eps * eps));
}
def code(x, eps): return eps + (x * (eps * eps))
function code(x, eps) return Float64(eps + Float64(x * Float64(eps * eps))) end
function tmp = code(x, eps) tmp = eps + (x * (eps * eps)); end
code[x_, eps_] := N[(eps + N[(x * N[(eps * eps), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\varepsilon + x \cdot \left(\varepsilon \cdot \varepsilon\right)
\end{array}
Initial program 61.4%
Taylor expanded in eps around 0
*-lowering-*.f64N/A
+-commutativeN/A
associate--l+N/A
+-commutativeN/A
+-lowering-+.f64N/A
Simplified99.7%
Taylor expanded in x around 0
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6498.1%
Simplified98.1%
(FPCore (x eps) :precision binary64 (* eps (+ (* eps x) 1.0)))
double code(double x, double eps) {
return eps * ((eps * x) + 1.0);
}
real(8) function code(x, eps)
real(8), intent (in) :: x
real(8), intent (in) :: eps
code = eps * ((eps * x) + 1.0d0)
end function
public static double code(double x, double eps) {
return eps * ((eps * x) + 1.0);
}
def code(x, eps): return eps * ((eps * x) + 1.0)
function code(x, eps) return Float64(eps * Float64(Float64(eps * x) + 1.0)) end
function tmp = code(x, eps) tmp = eps * ((eps * x) + 1.0); end
code[x_, eps_] := N[(eps * N[(N[(eps * x), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\varepsilon \cdot \left(\varepsilon \cdot x + 1\right)
\end{array}
Initial program 61.4%
Taylor expanded in eps around 0
*-lowering-*.f64N/A
+-commutativeN/A
associate--l+N/A
+-commutativeN/A
+-lowering-+.f64N/A
Simplified99.7%
Taylor expanded in x around 0
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6498.1%
Simplified98.1%
Final simplification98.1%
(FPCore (x eps) :precision binary64 eps)
double code(double x, double eps) {
return eps;
}
real(8) function code(x, eps)
real(8), intent (in) :: x
real(8), intent (in) :: eps
code = eps
end function
public static double code(double x, double eps) {
return eps;
}
def code(x, eps): return eps
function code(x, eps) return eps end
function tmp = code(x, eps) tmp = eps; end
code[x_, eps_] := eps
\begin{array}{l}
\\
\varepsilon
\end{array}
Initial program 61.4%
Taylor expanded in eps around 0
*-lowering-*.f64N/A
+-commutativeN/A
associate--l+N/A
+-commutativeN/A
+-lowering-+.f64N/A
Simplified99.7%
Taylor expanded in x around 0
Simplified98.1%
(FPCore (x eps) :precision binary64 0.0)
double code(double x, double eps) {
return 0.0;
}
real(8) function code(x, eps)
real(8), intent (in) :: x
real(8), intent (in) :: eps
code = 0.0d0
end function
public static double code(double x, double eps) {
return 0.0;
}
def code(x, eps): return 0.0
function code(x, eps) return 0.0 end
function tmp = code(x, eps) tmp = 0.0; end
code[x_, eps_] := 0.0
\begin{array}{l}
\\
0
\end{array}
Initial program 61.4%
sub-negN/A
+-commutativeN/A
tan-quotN/A
div-invN/A
distribute-rgt-neg-inN/A
accelerator-lowering-fma.f64N/A
sin-lowering-sin.f64N/A
neg-lowering-neg.f64N/A
/-lowering-/.f64N/A
cos-lowering-cos.f64N/A
tan-lowering-tan.f64N/A
+-lowering-+.f6461.4%
Applied egg-rr61.4%
Taylor expanded in eps around 0
distribute-lft1-inN/A
metadata-evalN/A
mul0-lft5.2%
Simplified5.2%
(FPCore (x eps) :precision binary64 (/ (sin eps) (* (cos x) (cos (+ x eps)))))
double code(double x, double eps) {
return sin(eps) / (cos(x) * cos((x + eps)));
}
real(8) function code(x, eps)
real(8), intent (in) :: x
real(8), intent (in) :: eps
code = sin(eps) / (cos(x) * cos((x + eps)))
end function
public static double code(double x, double eps) {
return Math.sin(eps) / (Math.cos(x) * Math.cos((x + eps)));
}
def code(x, eps): return math.sin(eps) / (math.cos(x) * math.cos((x + eps)))
function code(x, eps) return Float64(sin(eps) / Float64(cos(x) * cos(Float64(x + eps)))) end
function tmp = code(x, eps) tmp = sin(eps) / (cos(x) * cos((x + eps))); end
code[x_, eps_] := N[(N[Sin[eps], $MachinePrecision] / N[(N[Cos[x], $MachinePrecision] * N[Cos[N[(x + eps), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\sin \varepsilon}{\cos x \cdot \cos \left(x + \varepsilon\right)}
\end{array}
(FPCore (x eps) :precision binary64 (- (/ (+ (tan x) (tan eps)) (- 1.0 (* (tan x) (tan eps)))) (tan x)))
double code(double x, double eps) {
return ((tan(x) + tan(eps)) / (1.0 - (tan(x) * tan(eps)))) - tan(x);
}
real(8) function code(x, eps)
real(8), intent (in) :: x
real(8), intent (in) :: eps
code = ((tan(x) + tan(eps)) / (1.0d0 - (tan(x) * tan(eps)))) - tan(x)
end function
public static double code(double x, double eps) {
return ((Math.tan(x) + Math.tan(eps)) / (1.0 - (Math.tan(x) * Math.tan(eps)))) - Math.tan(x);
}
def code(x, eps): return ((math.tan(x) + math.tan(eps)) / (1.0 - (math.tan(x) * math.tan(eps)))) - math.tan(x)
function code(x, eps) return Float64(Float64(Float64(tan(x) + tan(eps)) / Float64(1.0 - Float64(tan(x) * tan(eps)))) - tan(x)) end
function tmp = code(x, eps) tmp = ((tan(x) + tan(eps)) / (1.0 - (tan(x) * tan(eps)))) - tan(x); end
code[x_, eps_] := N[(N[(N[(N[Tan[x], $MachinePrecision] + N[Tan[eps], $MachinePrecision]), $MachinePrecision] / N[(1.0 - N[(N[Tan[x], $MachinePrecision] * N[Tan[eps], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[Tan[x], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\tan x + \tan \varepsilon}{1 - \tan x \cdot \tan \varepsilon} - \tan x
\end{array}
(FPCore (x eps) :precision binary64 (+ eps (* (* eps (tan x)) (tan x))))
double code(double x, double eps) {
return eps + ((eps * tan(x)) * tan(x));
}
real(8) function code(x, eps)
real(8), intent (in) :: x
real(8), intent (in) :: eps
code = eps + ((eps * tan(x)) * tan(x))
end function
public static double code(double x, double eps) {
return eps + ((eps * Math.tan(x)) * Math.tan(x));
}
def code(x, eps): return eps + ((eps * math.tan(x)) * math.tan(x))
function code(x, eps) return Float64(eps + Float64(Float64(eps * tan(x)) * tan(x))) end
function tmp = code(x, eps) tmp = eps + ((eps * tan(x)) * tan(x)); end
code[x_, eps_] := N[(eps + N[(N[(eps * N[Tan[x], $MachinePrecision]), $MachinePrecision] * N[Tan[x], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\varepsilon + \left(\varepsilon \cdot \tan x\right) \cdot \tan x
\end{array}
herbie shell --seed 2024191
(FPCore (x eps)
:name "2tan (problem 3.3.2)"
:precision binary64
:pre (and (and (and (<= -10000.0 x) (<= x 10000.0)) (< (* 1e-16 (fabs x)) eps)) (< eps (fabs x)))
:alt
(! :herbie-platform default (/ (sin eps) (* (cos x) (cos (+ x eps)))))
:alt
(! :herbie-platform default (- (/ (+ (tan x) (tan eps)) (- 1 (* (tan x) (tan eps)))) (tan x)))
:alt
(! :herbie-platform default (+ eps (* eps (tan x) (tan x))))
(- (tan (+ x eps)) (tan x)))