
(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 11 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 (* (cos eps) (cos x))) (t_1 (* (sin eps) (sin x))))
(/
(/ (sin eps) (cos x))
(/
(- (pow t_0 3.0) (pow t_1 3.0))
(+ (+ (* t_1 t_0) (pow t_1 2.0)) (pow t_0 2.0))))))
double code(double x, double eps) {
double t_0 = cos(eps) * cos(x);
double t_1 = sin(eps) * sin(x);
return (sin(eps) / cos(x)) / ((pow(t_0, 3.0) - pow(t_1, 3.0)) / (((t_1 * t_0) + pow(t_1, 2.0)) + pow(t_0, 2.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 = cos(eps) * cos(x)
t_1 = sin(eps) * sin(x)
code = (sin(eps) / cos(x)) / (((t_0 ** 3.0d0) - (t_1 ** 3.0d0)) / (((t_1 * t_0) + (t_1 ** 2.0d0)) + (t_0 ** 2.0d0)))
end function
public static double code(double x, double eps) {
double t_0 = Math.cos(eps) * Math.cos(x);
double t_1 = Math.sin(eps) * Math.sin(x);
return (Math.sin(eps) / Math.cos(x)) / ((Math.pow(t_0, 3.0) - Math.pow(t_1, 3.0)) / (((t_1 * t_0) + Math.pow(t_1, 2.0)) + Math.pow(t_0, 2.0)));
}
def code(x, eps): t_0 = math.cos(eps) * math.cos(x) t_1 = math.sin(eps) * math.sin(x) return (math.sin(eps) / math.cos(x)) / ((math.pow(t_0, 3.0) - math.pow(t_1, 3.0)) / (((t_1 * t_0) + math.pow(t_1, 2.0)) + math.pow(t_0, 2.0)))
function code(x, eps) t_0 = Float64(cos(eps) * cos(x)) t_1 = Float64(sin(eps) * sin(x)) return Float64(Float64(sin(eps) / cos(x)) / Float64(Float64((t_0 ^ 3.0) - (t_1 ^ 3.0)) / Float64(Float64(Float64(t_1 * t_0) + (t_1 ^ 2.0)) + (t_0 ^ 2.0)))) end
function tmp = code(x, eps) t_0 = cos(eps) * cos(x); t_1 = sin(eps) * sin(x); tmp = (sin(eps) / cos(x)) / (((t_0 ^ 3.0) - (t_1 ^ 3.0)) / (((t_1 * t_0) + (t_1 ^ 2.0)) + (t_0 ^ 2.0))); end
code[x_, eps_] := Block[{t$95$0 = N[(N[Cos[eps], $MachinePrecision] * N[Cos[x], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[Sin[eps], $MachinePrecision] * N[Sin[x], $MachinePrecision]), $MachinePrecision]}, N[(N[(N[Sin[eps], $MachinePrecision] / N[Cos[x], $MachinePrecision]), $MachinePrecision] / N[(N[(N[Power[t$95$0, 3.0], $MachinePrecision] - N[Power[t$95$1, 3.0], $MachinePrecision]), $MachinePrecision] / N[(N[(N[(t$95$1 * t$95$0), $MachinePrecision] + N[Power[t$95$1, 2.0], $MachinePrecision]), $MachinePrecision] + N[Power[t$95$0, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \cos \varepsilon \cdot \cos x\\
t_1 := \sin \varepsilon \cdot \sin x\\
\frac{\frac{\sin \varepsilon}{\cos x}}{\frac{{t\_0}^{3} - {t\_1}^{3}}{\left(t\_1 \cdot t\_0 + {t\_1}^{2}\right) + {t\_0}^{2}}}
\end{array}
\end{array}
Initial program 60.4%
lift--.f64N/A
lift-tan.f64N/A
tan-quotN/A
lift-tan.f64N/A
tan-quotN/A
frac-subN/A
clear-numN/A
lower-/.f64N/A
lower-/.f64N/A
lower-*.f64N/A
lower-cos.f64N/A
lift-+.f64N/A
+-commutativeN/A
lower-+.f64N/A
lower-cos.f64N/A
sin-diffN/A
lower-sin.f64N/A
lower--.f6460.4
lift-+.f64N/A
+-commutativeN/A
lower-+.f6460.4
Applied rewrites60.4%
lift-/.f64N/A
lift-/.f64N/A
clear-numN/A
lift-*.f64N/A
*-commutativeN/A
associate-/r*N/A
lower-/.f64N/A
lower-/.f6460.4
lift--.f64N/A
lift-+.f64N/A
associate--l+N/A
+-inversesN/A
+-commutativeN/A
lower-+.f6499.7
Applied rewrites99.7%
lift-cos.f64N/A
lift-+.f64N/A
cos-sumN/A
flip3--N/A
lower-/.f64N/A
Applied rewrites100.0%
Final simplification100.0%
(FPCore (x eps) :precision binary64 (/ (sin eps) (* (- (* (cos eps) (cos x)) (* (sin eps) (sin x))) (cos x))))
double code(double x, double eps) {
return sin(eps) / (((cos(eps) * cos(x)) - (sin(eps) * sin(x))) * cos(x));
}
real(8) function code(x, eps)
real(8), intent (in) :: x
real(8), intent (in) :: eps
code = sin(eps) / (((cos(eps) * cos(x)) - (sin(eps) * sin(x))) * cos(x))
end function
public static double code(double x, double eps) {
return Math.sin(eps) / (((Math.cos(eps) * Math.cos(x)) - (Math.sin(eps) * Math.sin(x))) * Math.cos(x));
}
def code(x, eps): return math.sin(eps) / (((math.cos(eps) * math.cos(x)) - (math.sin(eps) * math.sin(x))) * math.cos(x))
function code(x, eps) return Float64(sin(eps) / Float64(Float64(Float64(cos(eps) * cos(x)) - Float64(sin(eps) * sin(x))) * cos(x))) end
function tmp = code(x, eps) tmp = sin(eps) / (((cos(eps) * cos(x)) - (sin(eps) * sin(x))) * cos(x)); end
code[x_, eps_] := N[(N[Sin[eps], $MachinePrecision] / N[(N[(N[(N[Cos[eps], $MachinePrecision] * N[Cos[x], $MachinePrecision]), $MachinePrecision] - N[(N[Sin[eps], $MachinePrecision] * N[Sin[x], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[Cos[x], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\sin \varepsilon}{\left(\cos \varepsilon \cdot \cos x - \sin \varepsilon \cdot \sin x\right) \cdot \cos x}
\end{array}
Initial program 60.4%
lift--.f64N/A
lift-tan.f64N/A
tan-quotN/A
lift-tan.f64N/A
tan-quotN/A
frac-subN/A
clear-numN/A
lower-/.f64N/A
lower-/.f64N/A
lower-*.f64N/A
lower-cos.f64N/A
lift-+.f64N/A
+-commutativeN/A
lower-+.f64N/A
lower-cos.f64N/A
sin-diffN/A
lower-sin.f64N/A
lower--.f6460.4
lift-+.f64N/A
+-commutativeN/A
lower-+.f6460.4
Applied rewrites60.4%
lift-/.f64N/A
lift-/.f64N/A
clear-numN/A
lower-/.f6460.4
lift--.f64N/A
lift-+.f64N/A
associate--l+N/A
+-inversesN/A
+-commutativeN/A
lower-+.f6499.7
Applied rewrites99.7%
lift-cos.f64N/A
lift-+.f64N/A
cos-sumN/A
lower--.f64N/A
lift-cos.f64N/A
lower-*.f64N/A
lower-cos.f64N/A
+-lft-identityN/A
lift-+.f64N/A
lift-sin.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-sin.f64100.0
lift-+.f64N/A
+-lft-identity100.0
Applied rewrites100.0%
Final simplification100.0%
(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(Float64(sin(eps) / 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[(N[Sin[eps], $MachinePrecision] / N[Cos[x], $MachinePrecision]), $MachinePrecision] / N[Cos[N[(x + eps), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\frac{\sin \varepsilon}{\cos x}}{\cos \left(x + \varepsilon\right)}
\end{array}
Initial program 60.4%
lift--.f64N/A
lift-tan.f64N/A
tan-quotN/A
lift-tan.f64N/A
tan-quotN/A
frac-subN/A
clear-numN/A
lower-/.f64N/A
lower-/.f64N/A
lower-*.f64N/A
lower-cos.f64N/A
lift-+.f64N/A
+-commutativeN/A
lower-+.f64N/A
lower-cos.f64N/A
sin-diffN/A
lower-sin.f64N/A
lower--.f6460.4
lift-+.f64N/A
+-commutativeN/A
lower-+.f6460.4
Applied rewrites60.4%
lift-/.f64N/A
lift-/.f64N/A
clear-numN/A
lift-*.f64N/A
*-commutativeN/A
associate-/r*N/A
lower-/.f64N/A
lower-/.f6460.4
lift--.f64N/A
lift-+.f64N/A
associate--l+N/A
+-inversesN/A
+-commutativeN/A
lower-+.f6499.7
Applied rewrites99.7%
lift-+.f64N/A
+-lft-identity99.7
lift-+.f64N/A
+-commutativeN/A
lower-+.f6499.7
Applied rewrites99.7%
(FPCore (x eps) :precision binary64 (/ (/ (sin eps) (cos (+ x eps))) (cos x)))
double code(double x, double eps) {
return (sin(eps) / cos((x + eps))) / cos(x);
}
real(8) function code(x, eps)
real(8), intent (in) :: x
real(8), intent (in) :: eps
code = (sin(eps) / cos((x + eps))) / cos(x)
end function
public static double code(double x, double eps) {
return (Math.sin(eps) / Math.cos((x + eps))) / Math.cos(x);
}
def code(x, eps): return (math.sin(eps) / math.cos((x + eps))) / math.cos(x)
function code(x, eps) return Float64(Float64(sin(eps) / cos(Float64(x + eps))) / cos(x)) end
function tmp = code(x, eps) tmp = (sin(eps) / cos((x + eps))) / cos(x); end
code[x_, eps_] := N[(N[(N[Sin[eps], $MachinePrecision] / N[Cos[N[(x + eps), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / N[Cos[x], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\frac{\sin \varepsilon}{\cos \left(x + \varepsilon\right)}}{\cos x}
\end{array}
Initial program 60.4%
lift--.f64N/A
lift-tan.f64N/A
tan-quotN/A
lift-tan.f64N/A
tan-quotN/A
frac-subN/A
clear-numN/A
lower-/.f64N/A
lower-/.f64N/A
lower-*.f64N/A
lower-cos.f64N/A
lift-+.f64N/A
+-commutativeN/A
lower-+.f64N/A
lower-cos.f64N/A
sin-diffN/A
lower-sin.f64N/A
lower--.f6460.4
lift-+.f64N/A
+-commutativeN/A
lower-+.f6460.4
Applied rewrites60.4%
Taylor expanded in x around inf
*-commutativeN/A
associate-/r*N/A
lower-/.f64N/A
lower-/.f64N/A
lower-sin.f64N/A
lower-cos.f64N/A
lower-+.f64N/A
lower-cos.f6499.7
Applied rewrites99.7%
Final simplification99.7%
(FPCore (x eps) :precision binary64 (/ (sin eps) (* (cos (+ x eps)) (cos x))))
double code(double x, double eps) {
return sin(eps) / (cos((x + eps)) * cos(x));
}
real(8) function code(x, eps)
real(8), intent (in) :: x
real(8), intent (in) :: eps
code = sin(eps) / (cos((x + eps)) * cos(x))
end function
public static double code(double x, double eps) {
return Math.sin(eps) / (Math.cos((x + eps)) * Math.cos(x));
}
def code(x, eps): return math.sin(eps) / (math.cos((x + eps)) * math.cos(x))
function code(x, eps) return Float64(sin(eps) / Float64(cos(Float64(x + eps)) * cos(x))) end
function tmp = code(x, eps) tmp = sin(eps) / (cos((x + eps)) * cos(x)); end
code[x_, eps_] := N[(N[Sin[eps], $MachinePrecision] / N[(N[Cos[N[(x + eps), $MachinePrecision]], $MachinePrecision] * N[Cos[x], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\sin \varepsilon}{\cos \left(x + \varepsilon\right) \cdot \cos x}
\end{array}
Initial program 60.4%
lift--.f64N/A
lift-tan.f64N/A
tan-quotN/A
lift-tan.f64N/A
tan-quotN/A
frac-subN/A
clear-numN/A
lower-/.f64N/A
lower-/.f64N/A
lower-*.f64N/A
lower-cos.f64N/A
lift-+.f64N/A
+-commutativeN/A
lower-+.f64N/A
lower-cos.f64N/A
sin-diffN/A
lower-sin.f64N/A
lower--.f6460.4
lift-+.f64N/A
+-commutativeN/A
lower-+.f6460.4
Applied rewrites60.4%
lift-/.f64N/A
lift-/.f64N/A
clear-numN/A
lift-*.f64N/A
*-commutativeN/A
associate-/r*N/A
lower-/.f64N/A
lower-/.f6460.4
lift--.f64N/A
lift-+.f64N/A
associate--l+N/A
+-inversesN/A
+-commutativeN/A
lower-+.f6499.7
Applied rewrites99.7%
lift-/.f64N/A
lift-/.f64N/A
associate-/l/N/A
lift-*.f64N/A
lift-/.f6499.7
lift-+.f64N/A
+-lft-identity99.7
lift-+.f64N/A
+-commutativeN/A
lower-+.f6499.7
Applied rewrites99.7%
(FPCore (x eps) :precision binary64 (/ (/ eps (cos x)) (cos (+ x eps))))
double code(double x, double eps) {
return (eps / cos(x)) / cos((x + eps));
}
real(8) function code(x, eps)
real(8), intent (in) :: x
real(8), intent (in) :: eps
code = (eps / cos(x)) / cos((x + eps))
end function
public static double code(double x, double eps) {
return (eps / Math.cos(x)) / Math.cos((x + eps));
}
def code(x, eps): return (eps / math.cos(x)) / math.cos((x + eps))
function code(x, eps) return Float64(Float64(eps / cos(x)) / cos(Float64(x + eps))) end
function tmp = code(x, eps) tmp = (eps / cos(x)) / cos((x + eps)); end
code[x_, eps_] := N[(N[(eps / N[Cos[x], $MachinePrecision]), $MachinePrecision] / N[Cos[N[(x + eps), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\frac{\varepsilon}{\cos x}}{\cos \left(x + \varepsilon\right)}
\end{array}
Initial program 60.4%
lift--.f64N/A
lift-tan.f64N/A
tan-quotN/A
lift-tan.f64N/A
tan-quotN/A
frac-subN/A
clear-numN/A
lower-/.f64N/A
lower-/.f64N/A
lower-*.f64N/A
lower-cos.f64N/A
lift-+.f64N/A
+-commutativeN/A
lower-+.f64N/A
lower-cos.f64N/A
sin-diffN/A
lower-sin.f64N/A
lower--.f6460.4
lift-+.f64N/A
+-commutativeN/A
lower-+.f6460.4
Applied rewrites60.4%
lift-/.f64N/A
lift-/.f64N/A
clear-numN/A
lift-*.f64N/A
*-commutativeN/A
associate-/r*N/A
lower-/.f64N/A
lower-/.f6460.4
lift--.f64N/A
lift-+.f64N/A
associate--l+N/A
+-inversesN/A
+-commutativeN/A
lower-+.f6499.7
Applied rewrites99.7%
Taylor expanded in eps around 0
lower-/.f64N/A
lower-cos.f6499.2
Applied rewrites99.2%
Final simplification99.2%
(FPCore (x eps) :precision binary64 (/ eps (pow (cos x) 2.0)))
double code(double x, double eps) {
return eps / pow(cos(x), 2.0);
}
real(8) function code(x, eps)
real(8), intent (in) :: x
real(8), intent (in) :: eps
code = eps / (cos(x) ** 2.0d0)
end function
public static double code(double x, double eps) {
return eps / Math.pow(Math.cos(x), 2.0);
}
def code(x, eps): return eps / math.pow(math.cos(x), 2.0)
function code(x, eps) return Float64(eps / (cos(x) ^ 2.0)) end
function tmp = code(x, eps) tmp = eps / (cos(x) ^ 2.0); end
code[x_, eps_] := N[(eps / N[Power[N[Cos[x], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\varepsilon}{{\cos x}^{2}}
\end{array}
Initial program 60.4%
lift--.f64N/A
lift-tan.f64N/A
tan-quotN/A
lift-tan.f64N/A
tan-quotN/A
frac-subN/A
clear-numN/A
lower-/.f64N/A
lower-/.f64N/A
lower-*.f64N/A
lower-cos.f64N/A
lift-+.f64N/A
+-commutativeN/A
lower-+.f64N/A
lower-cos.f64N/A
sin-diffN/A
lower-sin.f64N/A
lower--.f6460.4
lift-+.f64N/A
+-commutativeN/A
lower-+.f6460.4
Applied rewrites60.4%
Taylor expanded in eps around 0
lower-/.f64N/A
lower-pow.f64N/A
lower-cos.f6498.6
Applied rewrites98.6%
(FPCore (x eps) :precision binary64 (* (/ eps (+ (cos (* 2.0 x)) 1.0)) 2.0))
double code(double x, double eps) {
return (eps / (cos((2.0 * x)) + 1.0)) * 2.0;
}
real(8) function code(x, eps)
real(8), intent (in) :: x
real(8), intent (in) :: eps
code = (eps / (cos((2.0d0 * x)) + 1.0d0)) * 2.0d0
end function
public static double code(double x, double eps) {
return (eps / (Math.cos((2.0 * x)) + 1.0)) * 2.0;
}
def code(x, eps): return (eps / (math.cos((2.0 * x)) + 1.0)) * 2.0
function code(x, eps) return Float64(Float64(eps / Float64(cos(Float64(2.0 * x)) + 1.0)) * 2.0) end
function tmp = code(x, eps) tmp = (eps / (cos((2.0 * x)) + 1.0)) * 2.0; end
code[x_, eps_] := N[(N[(eps / N[(N[Cos[N[(2.0 * x), $MachinePrecision]], $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision] * 2.0), $MachinePrecision]
\begin{array}{l}
\\
\frac{\varepsilon}{\cos \left(2 \cdot x\right) + 1} \cdot 2
\end{array}
Initial program 60.4%
lift--.f64N/A
lift-tan.f64N/A
tan-quotN/A
lift-tan.f64N/A
tan-quotN/A
frac-subN/A
clear-numN/A
lower-/.f64N/A
lower-/.f64N/A
lower-*.f64N/A
lower-cos.f64N/A
lift-+.f64N/A
+-commutativeN/A
lower-+.f64N/A
lower-cos.f64N/A
sin-diffN/A
lower-sin.f64N/A
lower--.f6460.4
lift-+.f64N/A
+-commutativeN/A
lower-+.f6460.4
Applied rewrites60.4%
lift-/.f64N/A
lift-/.f64N/A
clear-numN/A
lift-*.f64N/A
*-commutativeN/A
associate-/r*N/A
lower-/.f64N/A
lower-/.f6460.4
lift--.f64N/A
lift-+.f64N/A
associate--l+N/A
+-inversesN/A
+-commutativeN/A
lower-+.f6499.7
Applied rewrites99.7%
lift-/.f64N/A
lift-/.f64N/A
associate-/l/N/A
lift-cos.f64N/A
lift-cos.f64N/A
cos-multN/A
associate-/r/N/A
lower-*.f64N/A
Applied rewrites99.7%
Taylor expanded in eps around 0
lower-/.f64N/A
+-commutativeN/A
lower-+.f64N/A
lower-cos.f64N/A
lower-*.f6498.6
Applied rewrites98.6%
(FPCore (x eps) :precision binary64 (fma (fma (* 0.6666666666666666 eps) (* x x) eps) (* x x) eps))
double code(double x, double eps) {
return fma(fma((0.6666666666666666 * eps), (x * x), eps), (x * x), eps);
}
function code(x, eps) return fma(fma(Float64(0.6666666666666666 * eps), Float64(x * x), eps), Float64(x * x), eps) end
code[x_, eps_] := N[(N[(N[(0.6666666666666666 * eps), $MachinePrecision] * N[(x * x), $MachinePrecision] + eps), $MachinePrecision] * N[(x * x), $MachinePrecision] + eps), $MachinePrecision]
\begin{array}{l}
\\
\mathsf{fma}\left(\mathsf{fma}\left(0.6666666666666666 \cdot \varepsilon, x \cdot x, \varepsilon\right), x \cdot x, \varepsilon\right)
\end{array}
Initial program 60.4%
lift--.f64N/A
lift-tan.f64N/A
tan-quotN/A
lift-tan.f64N/A
tan-quotN/A
frac-subN/A
clear-numN/A
lower-/.f64N/A
lower-/.f64N/A
lower-*.f64N/A
lower-cos.f64N/A
lift-+.f64N/A
+-commutativeN/A
lower-+.f64N/A
lower-cos.f64N/A
sin-diffN/A
lower-sin.f64N/A
lower--.f6460.4
lift-+.f64N/A
+-commutativeN/A
lower-+.f6460.4
Applied rewrites60.4%
lift-/.f64N/A
lift-/.f64N/A
clear-numN/A
lift-*.f64N/A
*-commutativeN/A
associate-/r*N/A
lower-/.f64N/A
lower-/.f6460.4
lift--.f64N/A
lift-+.f64N/A
associate--l+N/A
+-inversesN/A
+-commutativeN/A
lower-+.f6499.7
Applied rewrites99.7%
Taylor expanded in eps around 0
lower-/.f64N/A
lower-pow.f64N/A
lower-cos.f6498.6
Applied rewrites98.6%
Taylor expanded in x around 0
Applied rewrites98.1%
(FPCore (x eps) :precision binary64 (fma (* x x) eps eps))
double code(double x, double eps) {
return fma((x * x), eps, eps);
}
function code(x, eps) return fma(Float64(x * x), eps, eps) end
code[x_, eps_] := N[(N[(x * x), $MachinePrecision] * eps + eps), $MachinePrecision]
\begin{array}{l}
\\
\mathsf{fma}\left(x \cdot x, \varepsilon, \varepsilon\right)
\end{array}
Initial program 60.4%
lift--.f64N/A
lift-tan.f64N/A
tan-quotN/A
lift-tan.f64N/A
tan-quotN/A
frac-subN/A
clear-numN/A
lower-/.f64N/A
lower-/.f64N/A
lower-*.f64N/A
lower-cos.f64N/A
lift-+.f64N/A
+-commutativeN/A
lower-+.f64N/A
lower-cos.f64N/A
sin-diffN/A
lower-sin.f64N/A
lower--.f6460.4
lift-+.f64N/A
+-commutativeN/A
lower-+.f6460.4
Applied rewrites60.4%
lift-/.f64N/A
lift-/.f64N/A
clear-numN/A
lift-*.f64N/A
*-commutativeN/A
associate-/r*N/A
lower-/.f64N/A
lower-/.f6460.4
lift--.f64N/A
lift-+.f64N/A
associate--l+N/A
+-inversesN/A
+-commutativeN/A
lower-+.f6499.7
Applied rewrites99.7%
Taylor expanded in eps around 0
lower-/.f64N/A
lower-pow.f64N/A
lower-cos.f6498.6
Applied rewrites98.6%
Taylor expanded in x around 0
Applied rewrites98.1%
(FPCore (x eps) :precision binary64 (fma (* x eps) eps eps))
double code(double x, double eps) {
return fma((x * eps), eps, eps);
}
function code(x, eps) return fma(Float64(x * eps), eps, eps) end
code[x_, eps_] := N[(N[(x * eps), $MachinePrecision] * eps + eps), $MachinePrecision]
\begin{array}{l}
\\
\mathsf{fma}\left(x \cdot \varepsilon, \varepsilon, \varepsilon\right)
\end{array}
Initial program 60.4%
Taylor expanded in eps around 0
Applied rewrites99.2%
Taylor expanded in x around 0
Applied rewrites97.8%
Taylor expanded in x around 0
Applied rewrites97.8%
Taylor expanded in x around inf
Applied rewrites97.8%
Final simplification97.8%
(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 2024331
(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 (+ eps (* eps (tan x) (tan x))))
(- (tan (+ x eps)) (tan x)))