
(FPCore (eps) :precision binary64 (/ (- eps (sin eps)) (- eps (tan eps))))
double code(double eps) {
return (eps - sin(eps)) / (eps - tan(eps));
}
real(8) function code(eps)
real(8), intent (in) :: eps
code = (eps - sin(eps)) / (eps - tan(eps))
end function
public static double code(double eps) {
return (eps - Math.sin(eps)) / (eps - Math.tan(eps));
}
def code(eps): return (eps - math.sin(eps)) / (eps - math.tan(eps))
function code(eps) return Float64(Float64(eps - sin(eps)) / Float64(eps - tan(eps))) end
function tmp = code(eps) tmp = (eps - sin(eps)) / (eps - tan(eps)); end
code[eps_] := N[(N[(eps - N[Sin[eps], $MachinePrecision]), $MachinePrecision] / N[(eps - N[Tan[eps], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\varepsilon - \sin \varepsilon}{\varepsilon - \tan \varepsilon}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 6 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (eps) :precision binary64 (/ (- eps (sin eps)) (- eps (tan eps))))
double code(double eps) {
return (eps - sin(eps)) / (eps - tan(eps));
}
real(8) function code(eps)
real(8), intent (in) :: eps
code = (eps - sin(eps)) / (eps - tan(eps))
end function
public static double code(double eps) {
return (eps - Math.sin(eps)) / (eps - Math.tan(eps));
}
def code(eps): return (eps - math.sin(eps)) / (eps - math.tan(eps))
function code(eps) return Float64(Float64(eps - sin(eps)) / Float64(eps - tan(eps))) end
function tmp = code(eps) tmp = (eps - sin(eps)) / (eps - tan(eps)); end
code[eps_] := N[(N[(eps - N[Sin[eps], $MachinePrecision]), $MachinePrecision] / N[(eps - N[Tan[eps], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\varepsilon - \sin \varepsilon}{\varepsilon - \tan \varepsilon}
\end{array}
(FPCore (eps)
:precision binary64
(let* ((t_0
(+
0.225
(*
(* eps eps)
(+ -0.009642857142857142 (* (* eps eps) 0.00024107142857142857)))))
(t_1 (* eps t_0)))
(/
(+
(*
(* eps (* eps eps))
(*
t_1
(*
(* eps eps)
(+
0.050625
(*
eps
(*
eps
(+
-0.004339285714285714
(* (* eps eps) 0.00014722576530612244))))))))
-0.125)
(+ (* (* eps eps) (* t_1 t_1)) (- 0.25 (* t_0 (* (* eps eps) -0.5)))))))
double code(double eps) {
double t_0 = 0.225 + ((eps * eps) * (-0.009642857142857142 + ((eps * eps) * 0.00024107142857142857)));
double t_1 = eps * t_0;
return (((eps * (eps * eps)) * (t_1 * ((eps * eps) * (0.050625 + (eps * (eps * (-0.004339285714285714 + ((eps * eps) * 0.00014722576530612244)))))))) + -0.125) / (((eps * eps) * (t_1 * t_1)) + (0.25 - (t_0 * ((eps * eps) * -0.5))));
}
real(8) function code(eps)
real(8), intent (in) :: eps
real(8) :: t_0
real(8) :: t_1
t_0 = 0.225d0 + ((eps * eps) * ((-0.009642857142857142d0) + ((eps * eps) * 0.00024107142857142857d0)))
t_1 = eps * t_0
code = (((eps * (eps * eps)) * (t_1 * ((eps * eps) * (0.050625d0 + (eps * (eps * ((-0.004339285714285714d0) + ((eps * eps) * 0.00014722576530612244d0)))))))) + (-0.125d0)) / (((eps * eps) * (t_1 * t_1)) + (0.25d0 - (t_0 * ((eps * eps) * (-0.5d0)))))
end function
public static double code(double eps) {
double t_0 = 0.225 + ((eps * eps) * (-0.009642857142857142 + ((eps * eps) * 0.00024107142857142857)));
double t_1 = eps * t_0;
return (((eps * (eps * eps)) * (t_1 * ((eps * eps) * (0.050625 + (eps * (eps * (-0.004339285714285714 + ((eps * eps) * 0.00014722576530612244)))))))) + -0.125) / (((eps * eps) * (t_1 * t_1)) + (0.25 - (t_0 * ((eps * eps) * -0.5))));
}
def code(eps): t_0 = 0.225 + ((eps * eps) * (-0.009642857142857142 + ((eps * eps) * 0.00024107142857142857))) t_1 = eps * t_0 return (((eps * (eps * eps)) * (t_1 * ((eps * eps) * (0.050625 + (eps * (eps * (-0.004339285714285714 + ((eps * eps) * 0.00014722576530612244)))))))) + -0.125) / (((eps * eps) * (t_1 * t_1)) + (0.25 - (t_0 * ((eps * eps) * -0.5))))
function code(eps) t_0 = Float64(0.225 + Float64(Float64(eps * eps) * Float64(-0.009642857142857142 + Float64(Float64(eps * eps) * 0.00024107142857142857)))) t_1 = Float64(eps * t_0) return Float64(Float64(Float64(Float64(eps * Float64(eps * eps)) * Float64(t_1 * Float64(Float64(eps * eps) * Float64(0.050625 + Float64(eps * Float64(eps * Float64(-0.004339285714285714 + Float64(Float64(eps * eps) * 0.00014722576530612244)))))))) + -0.125) / Float64(Float64(Float64(eps * eps) * Float64(t_1 * t_1)) + Float64(0.25 - Float64(t_0 * Float64(Float64(eps * eps) * -0.5))))) end
function tmp = code(eps) t_0 = 0.225 + ((eps * eps) * (-0.009642857142857142 + ((eps * eps) * 0.00024107142857142857))); t_1 = eps * t_0; tmp = (((eps * (eps * eps)) * (t_1 * ((eps * eps) * (0.050625 + (eps * (eps * (-0.004339285714285714 + ((eps * eps) * 0.00014722576530612244)))))))) + -0.125) / (((eps * eps) * (t_1 * t_1)) + (0.25 - (t_0 * ((eps * eps) * -0.5)))); end
code[eps_] := Block[{t$95$0 = N[(0.225 + N[(N[(eps * eps), $MachinePrecision] * N[(-0.009642857142857142 + N[(N[(eps * eps), $MachinePrecision] * 0.00024107142857142857), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(eps * t$95$0), $MachinePrecision]}, N[(N[(N[(N[(eps * N[(eps * eps), $MachinePrecision]), $MachinePrecision] * N[(t$95$1 * N[(N[(eps * eps), $MachinePrecision] * N[(0.050625 + N[(eps * N[(eps * N[(-0.004339285714285714 + N[(N[(eps * eps), $MachinePrecision] * 0.00014722576530612244), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + -0.125), $MachinePrecision] / N[(N[(N[(eps * eps), $MachinePrecision] * N[(t$95$1 * t$95$1), $MachinePrecision]), $MachinePrecision] + N[(0.25 - N[(t$95$0 * N[(N[(eps * eps), $MachinePrecision] * -0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := 0.225 + \left(\varepsilon \cdot \varepsilon\right) \cdot \left(-0.009642857142857142 + \left(\varepsilon \cdot \varepsilon\right) \cdot 0.00024107142857142857\right)\\
t_1 := \varepsilon \cdot t\_0\\
\frac{\left(\varepsilon \cdot \left(\varepsilon \cdot \varepsilon\right)\right) \cdot \left(t\_1 \cdot \left(\left(\varepsilon \cdot \varepsilon\right) \cdot \left(0.050625 + \varepsilon \cdot \left(\varepsilon \cdot \left(-0.004339285714285714 + \left(\varepsilon \cdot \varepsilon\right) \cdot 0.00014722576530612244\right)\right)\right)\right)\right) + -0.125}{\left(\varepsilon \cdot \varepsilon\right) \cdot \left(t\_1 \cdot t\_1\right) + \left(0.25 - t\_0 \cdot \left(\left(\varepsilon \cdot \varepsilon\right) \cdot -0.5\right)\right)}
\end{array}
\end{array}
Initial program 1.4%
Taylor expanded in eps around 0
sub-negN/A
+-lowering-+.f64N/A
Simplified99.9%
flip3-+N/A
/-lowering-/.f64N/A
Applied egg-rr99.9%
Taylor expanded in eps around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6499.9%
Simplified99.9%
Taylor expanded in eps around 0
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6499.9%
Simplified99.9%
(FPCore (eps)
:precision binary64
(let* ((t_0
(+
0.225
(*
(* eps eps)
(+ -0.009642857142857142 (* (* eps eps) 0.00024107142857142857)))))
(t_1 (* eps t_0)))
(/
(+
-0.125
(*
(* eps (* eps eps))
(* t_1 (* t_1 (* eps (+ 0.225 (* (* eps eps) -0.009642857142857142)))))))
(+
(- 0.25 (* t_0 (* (* eps eps) -0.5)))
(*
(* (* eps eps) (* eps eps))
(+
0.050625
(*
(* eps eps)
(+ -0.004339285714285714 (* (* eps eps) 0.00020146683673469387)))))))))
double code(double eps) {
double t_0 = 0.225 + ((eps * eps) * (-0.009642857142857142 + ((eps * eps) * 0.00024107142857142857)));
double t_1 = eps * t_0;
return (-0.125 + ((eps * (eps * eps)) * (t_1 * (t_1 * (eps * (0.225 + ((eps * eps) * -0.009642857142857142))))))) / ((0.25 - (t_0 * ((eps * eps) * -0.5))) + (((eps * eps) * (eps * eps)) * (0.050625 + ((eps * eps) * (-0.004339285714285714 + ((eps * eps) * 0.00020146683673469387))))));
}
real(8) function code(eps)
real(8), intent (in) :: eps
real(8) :: t_0
real(8) :: t_1
t_0 = 0.225d0 + ((eps * eps) * ((-0.009642857142857142d0) + ((eps * eps) * 0.00024107142857142857d0)))
t_1 = eps * t_0
code = ((-0.125d0) + ((eps * (eps * eps)) * (t_1 * (t_1 * (eps * (0.225d0 + ((eps * eps) * (-0.009642857142857142d0)))))))) / ((0.25d0 - (t_0 * ((eps * eps) * (-0.5d0)))) + (((eps * eps) * (eps * eps)) * (0.050625d0 + ((eps * eps) * ((-0.004339285714285714d0) + ((eps * eps) * 0.00020146683673469387d0))))))
end function
public static double code(double eps) {
double t_0 = 0.225 + ((eps * eps) * (-0.009642857142857142 + ((eps * eps) * 0.00024107142857142857)));
double t_1 = eps * t_0;
return (-0.125 + ((eps * (eps * eps)) * (t_1 * (t_1 * (eps * (0.225 + ((eps * eps) * -0.009642857142857142))))))) / ((0.25 - (t_0 * ((eps * eps) * -0.5))) + (((eps * eps) * (eps * eps)) * (0.050625 + ((eps * eps) * (-0.004339285714285714 + ((eps * eps) * 0.00020146683673469387))))));
}
def code(eps): t_0 = 0.225 + ((eps * eps) * (-0.009642857142857142 + ((eps * eps) * 0.00024107142857142857))) t_1 = eps * t_0 return (-0.125 + ((eps * (eps * eps)) * (t_1 * (t_1 * (eps * (0.225 + ((eps * eps) * -0.009642857142857142))))))) / ((0.25 - (t_0 * ((eps * eps) * -0.5))) + (((eps * eps) * (eps * eps)) * (0.050625 + ((eps * eps) * (-0.004339285714285714 + ((eps * eps) * 0.00020146683673469387))))))
function code(eps) t_0 = Float64(0.225 + Float64(Float64(eps * eps) * Float64(-0.009642857142857142 + Float64(Float64(eps * eps) * 0.00024107142857142857)))) t_1 = Float64(eps * t_0) return Float64(Float64(-0.125 + Float64(Float64(eps * Float64(eps * eps)) * Float64(t_1 * Float64(t_1 * Float64(eps * Float64(0.225 + Float64(Float64(eps * eps) * -0.009642857142857142))))))) / Float64(Float64(0.25 - Float64(t_0 * Float64(Float64(eps * eps) * -0.5))) + Float64(Float64(Float64(eps * eps) * Float64(eps * eps)) * Float64(0.050625 + Float64(Float64(eps * eps) * Float64(-0.004339285714285714 + Float64(Float64(eps * eps) * 0.00020146683673469387))))))) end
function tmp = code(eps) t_0 = 0.225 + ((eps * eps) * (-0.009642857142857142 + ((eps * eps) * 0.00024107142857142857))); t_1 = eps * t_0; tmp = (-0.125 + ((eps * (eps * eps)) * (t_1 * (t_1 * (eps * (0.225 + ((eps * eps) * -0.009642857142857142))))))) / ((0.25 - (t_0 * ((eps * eps) * -0.5))) + (((eps * eps) * (eps * eps)) * (0.050625 + ((eps * eps) * (-0.004339285714285714 + ((eps * eps) * 0.00020146683673469387)))))); end
code[eps_] := Block[{t$95$0 = N[(0.225 + N[(N[(eps * eps), $MachinePrecision] * N[(-0.009642857142857142 + N[(N[(eps * eps), $MachinePrecision] * 0.00024107142857142857), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(eps * t$95$0), $MachinePrecision]}, N[(N[(-0.125 + N[(N[(eps * N[(eps * eps), $MachinePrecision]), $MachinePrecision] * N[(t$95$1 * N[(t$95$1 * N[(eps * N[(0.225 + N[(N[(eps * eps), $MachinePrecision] * -0.009642857142857142), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(N[(0.25 - N[(t$95$0 * N[(N[(eps * eps), $MachinePrecision] * -0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(N[(eps * eps), $MachinePrecision] * N[(eps * eps), $MachinePrecision]), $MachinePrecision] * N[(0.050625 + N[(N[(eps * eps), $MachinePrecision] * N[(-0.004339285714285714 + N[(N[(eps * eps), $MachinePrecision] * 0.00020146683673469387), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := 0.225 + \left(\varepsilon \cdot \varepsilon\right) \cdot \left(-0.009642857142857142 + \left(\varepsilon \cdot \varepsilon\right) \cdot 0.00024107142857142857\right)\\
t_1 := \varepsilon \cdot t\_0\\
\frac{-0.125 + \left(\varepsilon \cdot \left(\varepsilon \cdot \varepsilon\right)\right) \cdot \left(t\_1 \cdot \left(t\_1 \cdot \left(\varepsilon \cdot \left(0.225 + \left(\varepsilon \cdot \varepsilon\right) \cdot -0.009642857142857142\right)\right)\right)\right)}{\left(0.25 - t\_0 \cdot \left(\left(\varepsilon \cdot \varepsilon\right) \cdot -0.5\right)\right) + \left(\left(\varepsilon \cdot \varepsilon\right) \cdot \left(\varepsilon \cdot \varepsilon\right)\right) \cdot \left(0.050625 + \left(\varepsilon \cdot \varepsilon\right) \cdot \left(-0.004339285714285714 + \left(\varepsilon \cdot \varepsilon\right) \cdot 0.00020146683673469387\right)\right)}
\end{array}
\end{array}
Initial program 1.4%
Taylor expanded in eps around 0
sub-negN/A
+-lowering-+.f64N/A
Simplified99.9%
flip3-+N/A
/-lowering-/.f64N/A
Applied egg-rr99.9%
Taylor expanded in eps around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6499.9%
Simplified99.9%
Taylor expanded in eps around 0
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6499.9%
Simplified99.9%
Final simplification99.9%
(FPCore (eps)
:precision binary64
(+
-0.5
(*
(* eps eps)
(+
0.225
(*
eps
(*
eps
(+ -0.009642857142857142 (* (* eps eps) 0.00024107142857142857))))))))
double code(double eps) {
return -0.5 + ((eps * eps) * (0.225 + (eps * (eps * (-0.009642857142857142 + ((eps * eps) * 0.00024107142857142857))))));
}
real(8) function code(eps)
real(8), intent (in) :: eps
code = (-0.5d0) + ((eps * eps) * (0.225d0 + (eps * (eps * ((-0.009642857142857142d0) + ((eps * eps) * 0.00024107142857142857d0))))))
end function
public static double code(double eps) {
return -0.5 + ((eps * eps) * (0.225 + (eps * (eps * (-0.009642857142857142 + ((eps * eps) * 0.00024107142857142857))))));
}
def code(eps): return -0.5 + ((eps * eps) * (0.225 + (eps * (eps * (-0.009642857142857142 + ((eps * eps) * 0.00024107142857142857))))))
function code(eps) return Float64(-0.5 + Float64(Float64(eps * eps) * Float64(0.225 + Float64(eps * Float64(eps * Float64(-0.009642857142857142 + Float64(Float64(eps * eps) * 0.00024107142857142857))))))) end
function tmp = code(eps) tmp = -0.5 + ((eps * eps) * (0.225 + (eps * (eps * (-0.009642857142857142 + ((eps * eps) * 0.00024107142857142857)))))); end
code[eps_] := N[(-0.5 + N[(N[(eps * eps), $MachinePrecision] * N[(0.225 + N[(eps * N[(eps * N[(-0.009642857142857142 + N[(N[(eps * eps), $MachinePrecision] * 0.00024107142857142857), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
-0.5 + \left(\varepsilon \cdot \varepsilon\right) \cdot \left(0.225 + \varepsilon \cdot \left(\varepsilon \cdot \left(-0.009642857142857142 + \left(\varepsilon \cdot \varepsilon\right) \cdot 0.00024107142857142857\right)\right)\right)
\end{array}
Initial program 1.4%
Taylor expanded in eps around 0
sub-negN/A
+-lowering-+.f64N/A
Simplified99.9%
Final simplification99.9%
(FPCore (eps) :precision binary64 (+ -0.5 (* eps (* eps (+ 0.225 (* (* eps eps) -0.009642857142857142))))))
double code(double eps) {
return -0.5 + (eps * (eps * (0.225 + ((eps * eps) * -0.009642857142857142))));
}
real(8) function code(eps)
real(8), intent (in) :: eps
code = (-0.5d0) + (eps * (eps * (0.225d0 + ((eps * eps) * (-0.009642857142857142d0)))))
end function
public static double code(double eps) {
return -0.5 + (eps * (eps * (0.225 + ((eps * eps) * -0.009642857142857142))));
}
def code(eps): return -0.5 + (eps * (eps * (0.225 + ((eps * eps) * -0.009642857142857142))))
function code(eps) return Float64(-0.5 + Float64(eps * Float64(eps * Float64(0.225 + Float64(Float64(eps * eps) * -0.009642857142857142))))) end
function tmp = code(eps) tmp = -0.5 + (eps * (eps * (0.225 + ((eps * eps) * -0.009642857142857142)))); end
code[eps_] := N[(-0.5 + N[(eps * N[(eps * N[(0.225 + N[(N[(eps * eps), $MachinePrecision] * -0.009642857142857142), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
-0.5 + \varepsilon \cdot \left(\varepsilon \cdot \left(0.225 + \left(\varepsilon \cdot \varepsilon\right) \cdot -0.009642857142857142\right)\right)
\end{array}
Initial program 1.4%
Taylor expanded in eps around 0
sub-negN/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
metadata-eval99.8%
Simplified99.8%
Final simplification99.8%
(FPCore (eps) :precision binary64 (+ -0.5 (* (* eps eps) 0.225)))
double code(double eps) {
return -0.5 + ((eps * eps) * 0.225);
}
real(8) function code(eps)
real(8), intent (in) :: eps
code = (-0.5d0) + ((eps * eps) * 0.225d0)
end function
public static double code(double eps) {
return -0.5 + ((eps * eps) * 0.225);
}
def code(eps): return -0.5 + ((eps * eps) * 0.225)
function code(eps) return Float64(-0.5 + Float64(Float64(eps * eps) * 0.225)) end
function tmp = code(eps) tmp = -0.5 + ((eps * eps) * 0.225); end
code[eps_] := N[(-0.5 + N[(N[(eps * eps), $MachinePrecision] * 0.225), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
-0.5 + \left(\varepsilon \cdot \varepsilon\right) \cdot 0.225
\end{array}
Initial program 1.4%
Taylor expanded in eps around 0
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
metadata-eval99.7%
Simplified99.7%
Final simplification99.7%
(FPCore (eps) :precision binary64 -0.5)
double code(double eps) {
return -0.5;
}
real(8) function code(eps)
real(8), intent (in) :: eps
code = -0.5d0
end function
public static double code(double eps) {
return -0.5;
}
def code(eps): return -0.5
function code(eps) return -0.5 end
function tmp = code(eps) tmp = -0.5; end
code[eps_] := -0.5
\begin{array}{l}
\\
-0.5
\end{array}
Initial program 1.4%
Taylor expanded in eps around 0
Simplified99.5%
(FPCore (eps)
:precision binary64
(let* ((t_0 (* (* (* eps eps) eps) eps)))
(+
(+ (+ -0.5 (/ (* 9.0 (* eps eps)) 40.0)) (/ (* -27.0 t_0) 2800.0))
(/ (* 27.0 (* (* t_0 eps) eps)) 112000.0))))
double code(double eps) {
double t_0 = ((eps * eps) * eps) * eps;
return ((-0.5 + ((9.0 * (eps * eps)) / 40.0)) + ((-27.0 * t_0) / 2800.0)) + ((27.0 * ((t_0 * eps) * eps)) / 112000.0);
}
real(8) function code(eps)
real(8), intent (in) :: eps
real(8) :: t_0
t_0 = ((eps * eps) * eps) * eps
code = (((-0.5d0) + ((9.0d0 * (eps * eps)) / 40.0d0)) + (((-27.0d0) * t_0) / 2800.0d0)) + ((27.0d0 * ((t_0 * eps) * eps)) / 112000.0d0)
end function
public static double code(double eps) {
double t_0 = ((eps * eps) * eps) * eps;
return ((-0.5 + ((9.0 * (eps * eps)) / 40.0)) + ((-27.0 * t_0) / 2800.0)) + ((27.0 * ((t_0 * eps) * eps)) / 112000.0);
}
def code(eps): t_0 = ((eps * eps) * eps) * eps return ((-0.5 + ((9.0 * (eps * eps)) / 40.0)) + ((-27.0 * t_0) / 2800.0)) + ((27.0 * ((t_0 * eps) * eps)) / 112000.0)
function code(eps) t_0 = Float64(Float64(Float64(eps * eps) * eps) * eps) return Float64(Float64(Float64(-0.5 + Float64(Float64(9.0 * Float64(eps * eps)) / 40.0)) + Float64(Float64(-27.0 * t_0) / 2800.0)) + Float64(Float64(27.0 * Float64(Float64(t_0 * eps) * eps)) / 112000.0)) end
function tmp = code(eps) t_0 = ((eps * eps) * eps) * eps; tmp = ((-0.5 + ((9.0 * (eps * eps)) / 40.0)) + ((-27.0 * t_0) / 2800.0)) + ((27.0 * ((t_0 * eps) * eps)) / 112000.0); end
code[eps_] := Block[{t$95$0 = N[(N[(N[(eps * eps), $MachinePrecision] * eps), $MachinePrecision] * eps), $MachinePrecision]}, N[(N[(N[(-0.5 + N[(N[(9.0 * N[(eps * eps), $MachinePrecision]), $MachinePrecision] / 40.0), $MachinePrecision]), $MachinePrecision] + N[(N[(-27.0 * t$95$0), $MachinePrecision] / 2800.0), $MachinePrecision]), $MachinePrecision] + N[(N[(27.0 * N[(N[(t$95$0 * eps), $MachinePrecision] * eps), $MachinePrecision]), $MachinePrecision] / 112000.0), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(\left(\varepsilon \cdot \varepsilon\right) \cdot \varepsilon\right) \cdot \varepsilon\\
\left(\left(-0.5 + \frac{9 \cdot \left(\varepsilon \cdot \varepsilon\right)}{40}\right) + \frac{-27 \cdot t\_0}{2800}\right) + \frac{27 \cdot \left(\left(t\_0 \cdot \varepsilon\right) \cdot \varepsilon\right)}{112000}
\end{array}
\end{array}
herbie shell --seed 2024158
(FPCore (eps)
:name "sintan (problem 3.4.5)"
:precision binary64
:pre (and (<= -0.4 eps) (<= eps 0.4))
:alt
(! :herbie-platform default (+ -1/2 (/ (* 9 (* eps eps)) 40) (/ (* -27 (* eps eps eps eps)) 2800) (/ (* 27 (* eps eps eps eps eps eps)) 112000)))
(/ (- eps (sin eps)) (- eps (tan eps))))