
(FPCore (x eps) :precision binary64 (- (cos (+ x eps)) (cos x)))
double code(double x, double eps) {
return cos((x + eps)) - cos(x);
}
real(8) function code(x, eps)
real(8), intent (in) :: x
real(8), intent (in) :: eps
code = cos((x + eps)) - cos(x)
end function
public static double code(double x, double eps) {
return Math.cos((x + eps)) - Math.cos(x);
}
def code(x, eps): return math.cos((x + eps)) - math.cos(x)
function code(x, eps) return Float64(cos(Float64(x + eps)) - cos(x)) end
function tmp = code(x, eps) tmp = cos((x + eps)) - cos(x); end
code[x_, eps_] := N[(N[Cos[N[(x + eps), $MachinePrecision]], $MachinePrecision] - N[Cos[x], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\cos \left(x + \varepsilon\right) - \cos x
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 18 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x eps) :precision binary64 (- (cos (+ x eps)) (cos x)))
double code(double x, double eps) {
return cos((x + eps)) - cos(x);
}
real(8) function code(x, eps)
real(8), intent (in) :: x
real(8), intent (in) :: eps
code = cos((x + eps)) - cos(x)
end function
public static double code(double x, double eps) {
return Math.cos((x + eps)) - Math.cos(x);
}
def code(x, eps): return math.cos((x + eps)) - math.cos(x)
function code(x, eps) return Float64(cos(Float64(x + eps)) - cos(x)) end
function tmp = code(x, eps) tmp = cos((x + eps)) - cos(x); end
code[x_, eps_] := N[(N[Cos[N[(x + eps), $MachinePrecision]], $MachinePrecision] - N[Cos[x], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\cos \left(x + \varepsilon\right) - \cos x
\end{array}
(FPCore (x eps) :precision binary64 (let* ((t_0 (sin (* eps 0.5)))) (* (fma (* (sin x) (cos (* eps 0.5))) t_0 (* t_0 (* t_0 (cos x)))) -2.0)))
double code(double x, double eps) {
double t_0 = sin((eps * 0.5));
return fma((sin(x) * cos((eps * 0.5))), t_0, (t_0 * (t_0 * cos(x)))) * -2.0;
}
function code(x, eps) t_0 = sin(Float64(eps * 0.5)) return Float64(fma(Float64(sin(x) * cos(Float64(eps * 0.5))), t_0, Float64(t_0 * Float64(t_0 * cos(x)))) * -2.0) end
code[x_, eps_] := Block[{t$95$0 = N[Sin[N[(eps * 0.5), $MachinePrecision]], $MachinePrecision]}, N[(N[(N[(N[Sin[x], $MachinePrecision] * N[Cos[N[(eps * 0.5), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] * t$95$0 + N[(t$95$0 * N[(t$95$0 * N[Cos[x], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * -2.0), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sin \left(\varepsilon \cdot 0.5\right)\\
\mathsf{fma}\left(\sin x \cdot \cos \left(\varepsilon \cdot 0.5\right), t\_0, t\_0 \cdot \left(t\_0 \cdot \cos x\right)\right) \cdot -2
\end{array}
\end{array}
Initial program 49.1%
diff-cosN/A
*-commutativeN/A
*-lowering-*.f64N/A
Applied egg-rr99.5%
Taylor expanded in eps around inf
metadata-evalN/A
cancel-sign-sub-invN/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
cancel-sign-sub-invN/A
metadata-evalN/A
+-commutativeN/A
distribute-rgt-inN/A
*-commutativeN/A
associate-*l*N/A
metadata-evalN/A
*-rgt-identityN/A
*-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6499.5%
Simplified99.5%
sin-sumN/A
distribute-rgt-inN/A
fma-defineN/A
fma-lowering-fma.f64N/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
cos-lowering-cos.f64N/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f6499.9%
Applied egg-rr99.9%
Final simplification99.9%
(FPCore (x eps) :precision binary64 (let* ((t_0 (sin (* eps 0.5)))) (* -2.0 (* t_0 (+ (* (sin x) (cos (* eps 0.5))) (* t_0 (cos x)))))))
double code(double x, double eps) {
double t_0 = sin((eps * 0.5));
return -2.0 * (t_0 * ((sin(x) * cos((eps * 0.5))) + (t_0 * cos(x))));
}
real(8) function code(x, eps)
real(8), intent (in) :: x
real(8), intent (in) :: eps
real(8) :: t_0
t_0 = sin((eps * 0.5d0))
code = (-2.0d0) * (t_0 * ((sin(x) * cos((eps * 0.5d0))) + (t_0 * cos(x))))
end function
public static double code(double x, double eps) {
double t_0 = Math.sin((eps * 0.5));
return -2.0 * (t_0 * ((Math.sin(x) * Math.cos((eps * 0.5))) + (t_0 * Math.cos(x))));
}
def code(x, eps): t_0 = math.sin((eps * 0.5)) return -2.0 * (t_0 * ((math.sin(x) * math.cos((eps * 0.5))) + (t_0 * math.cos(x))))
function code(x, eps) t_0 = sin(Float64(eps * 0.5)) return Float64(-2.0 * Float64(t_0 * Float64(Float64(sin(x) * cos(Float64(eps * 0.5))) + Float64(t_0 * cos(x))))) end
function tmp = code(x, eps) t_0 = sin((eps * 0.5)); tmp = -2.0 * (t_0 * ((sin(x) * cos((eps * 0.5))) + (t_0 * cos(x)))); end
code[x_, eps_] := Block[{t$95$0 = N[Sin[N[(eps * 0.5), $MachinePrecision]], $MachinePrecision]}, N[(-2.0 * N[(t$95$0 * N[(N[(N[Sin[x], $MachinePrecision] * N[Cos[N[(eps * 0.5), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] + N[(t$95$0 * N[Cos[x], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sin \left(\varepsilon \cdot 0.5\right)\\
-2 \cdot \left(t\_0 \cdot \left(\sin x \cdot \cos \left(\varepsilon \cdot 0.5\right) + t\_0 \cdot \cos x\right)\right)
\end{array}
\end{array}
Initial program 49.1%
diff-cosN/A
*-commutativeN/A
*-lowering-*.f64N/A
Applied egg-rr99.5%
Taylor expanded in eps around inf
metadata-evalN/A
cancel-sign-sub-invN/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
cancel-sign-sub-invN/A
metadata-evalN/A
+-commutativeN/A
distribute-rgt-inN/A
*-commutativeN/A
associate-*l*N/A
metadata-evalN/A
*-rgt-identityN/A
*-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6499.5%
Simplified99.5%
sin-sumN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
cos-lowering-cos.f64N/A
*-lowering-*.f6499.7%
Applied egg-rr99.7%
Final simplification99.7%
(FPCore (x eps)
:precision binary64
(let* ((t_0
(*
eps
(+
0.5
(*
(* eps eps)
(+
-0.020833333333333332
(*
(* eps eps)
(+
0.00026041666666666666
(* (* eps eps) -1.5500992063492063e-6)))))))))
(*
-2.0
(+
(* (* (sin x) (cos (* eps 0.5))) t_0)
(* (* (sin (* eps 0.5)) (cos x)) t_0)))))
double code(double x, double eps) {
double t_0 = eps * (0.5 + ((eps * eps) * (-0.020833333333333332 + ((eps * eps) * (0.00026041666666666666 + ((eps * eps) * -1.5500992063492063e-6))))));
return -2.0 * (((sin(x) * cos((eps * 0.5))) * t_0) + ((sin((eps * 0.5)) * cos(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 = eps * (0.5d0 + ((eps * eps) * ((-0.020833333333333332d0) + ((eps * eps) * (0.00026041666666666666d0 + ((eps * eps) * (-1.5500992063492063d-6)))))))
code = (-2.0d0) * (((sin(x) * cos((eps * 0.5d0))) * t_0) + ((sin((eps * 0.5d0)) * cos(x)) * t_0))
end function
public static double code(double x, double eps) {
double t_0 = eps * (0.5 + ((eps * eps) * (-0.020833333333333332 + ((eps * eps) * (0.00026041666666666666 + ((eps * eps) * -1.5500992063492063e-6))))));
return -2.0 * (((Math.sin(x) * Math.cos((eps * 0.5))) * t_0) + ((Math.sin((eps * 0.5)) * Math.cos(x)) * t_0));
}
def code(x, eps): t_0 = eps * (0.5 + ((eps * eps) * (-0.020833333333333332 + ((eps * eps) * (0.00026041666666666666 + ((eps * eps) * -1.5500992063492063e-6)))))) return -2.0 * (((math.sin(x) * math.cos((eps * 0.5))) * t_0) + ((math.sin((eps * 0.5)) * math.cos(x)) * t_0))
function code(x, eps) t_0 = Float64(eps * Float64(0.5 + Float64(Float64(eps * eps) * Float64(-0.020833333333333332 + Float64(Float64(eps * eps) * Float64(0.00026041666666666666 + Float64(Float64(eps * eps) * -1.5500992063492063e-6))))))) return Float64(-2.0 * Float64(Float64(Float64(sin(x) * cos(Float64(eps * 0.5))) * t_0) + Float64(Float64(sin(Float64(eps * 0.5)) * cos(x)) * t_0))) end
function tmp = code(x, eps) t_0 = eps * (0.5 + ((eps * eps) * (-0.020833333333333332 + ((eps * eps) * (0.00026041666666666666 + ((eps * eps) * -1.5500992063492063e-6)))))); tmp = -2.0 * (((sin(x) * cos((eps * 0.5))) * t_0) + ((sin((eps * 0.5)) * cos(x)) * t_0)); end
code[x_, eps_] := Block[{t$95$0 = N[(eps * N[(0.5 + N[(N[(eps * eps), $MachinePrecision] * N[(-0.020833333333333332 + N[(N[(eps * eps), $MachinePrecision] * N[(0.00026041666666666666 + N[(N[(eps * eps), $MachinePrecision] * -1.5500992063492063e-6), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, N[(-2.0 * N[(N[(N[(N[Sin[x], $MachinePrecision] * N[Cos[N[(eps * 0.5), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] * t$95$0), $MachinePrecision] + N[(N[(N[Sin[N[(eps * 0.5), $MachinePrecision]], $MachinePrecision] * N[Cos[x], $MachinePrecision]), $MachinePrecision] * t$95$0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \varepsilon \cdot \left(0.5 + \left(\varepsilon \cdot \varepsilon\right) \cdot \left(-0.020833333333333332 + \left(\varepsilon \cdot \varepsilon\right) \cdot \left(0.00026041666666666666 + \left(\varepsilon \cdot \varepsilon\right) \cdot -1.5500992063492063 \cdot 10^{-6}\right)\right)\right)\\
-2 \cdot \left(\left(\sin x \cdot \cos \left(\varepsilon \cdot 0.5\right)\right) \cdot t\_0 + \left(\sin \left(\varepsilon \cdot 0.5\right) \cdot \cos x\right) \cdot t\_0\right)
\end{array}
\end{array}
Initial program 49.1%
diff-cosN/A
*-commutativeN/A
*-lowering-*.f64N/A
Applied egg-rr99.5%
Taylor expanded in eps around inf
metadata-evalN/A
cancel-sign-sub-invN/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
cancel-sign-sub-invN/A
metadata-evalN/A
+-commutativeN/A
distribute-rgt-inN/A
*-commutativeN/A
associate-*l*N/A
metadata-evalN/A
*-rgt-identityN/A
*-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6499.5%
Simplified99.5%
Taylor expanded in eps around 0
*-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
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6499.5%
Simplified99.5%
sin-sumN/A
*-commutativeN/A
distribute-rgt-inN/A
+-lowering-+.f64N/A
Applied egg-rr99.6%
Final simplification99.6%
(FPCore (x eps)
:precision binary64
(*
-2.0
(*
(+ (* (sin x) (cos (* eps 0.5))) (* (sin (* eps 0.5)) (cos x)))
(*
eps
(+
0.5
(*
eps
(*
eps
(+
-0.020833333333333332
(*
(* eps eps)
(+
0.00026041666666666666
(* (* eps eps) -1.5500992063492063e-6)))))))))))
double code(double x, double eps) {
return -2.0 * (((sin(x) * cos((eps * 0.5))) + (sin((eps * 0.5)) * cos(x))) * (eps * (0.5 + (eps * (eps * (-0.020833333333333332 + ((eps * eps) * (0.00026041666666666666 + ((eps * eps) * -1.5500992063492063e-6)))))))));
}
real(8) function code(x, eps)
real(8), intent (in) :: x
real(8), intent (in) :: eps
code = (-2.0d0) * (((sin(x) * cos((eps * 0.5d0))) + (sin((eps * 0.5d0)) * cos(x))) * (eps * (0.5d0 + (eps * (eps * ((-0.020833333333333332d0) + ((eps * eps) * (0.00026041666666666666d0 + ((eps * eps) * (-1.5500992063492063d-6))))))))))
end function
public static double code(double x, double eps) {
return -2.0 * (((Math.sin(x) * Math.cos((eps * 0.5))) + (Math.sin((eps * 0.5)) * Math.cos(x))) * (eps * (0.5 + (eps * (eps * (-0.020833333333333332 + ((eps * eps) * (0.00026041666666666666 + ((eps * eps) * -1.5500992063492063e-6)))))))));
}
def code(x, eps): return -2.0 * (((math.sin(x) * math.cos((eps * 0.5))) + (math.sin((eps * 0.5)) * math.cos(x))) * (eps * (0.5 + (eps * (eps * (-0.020833333333333332 + ((eps * eps) * (0.00026041666666666666 + ((eps * eps) * -1.5500992063492063e-6)))))))))
function code(x, eps) return Float64(-2.0 * Float64(Float64(Float64(sin(x) * cos(Float64(eps * 0.5))) + Float64(sin(Float64(eps * 0.5)) * cos(x))) * Float64(eps * Float64(0.5 + Float64(eps * Float64(eps * Float64(-0.020833333333333332 + Float64(Float64(eps * eps) * Float64(0.00026041666666666666 + Float64(Float64(eps * eps) * -1.5500992063492063e-6)))))))))) end
function tmp = code(x, eps) tmp = -2.0 * (((sin(x) * cos((eps * 0.5))) + (sin((eps * 0.5)) * cos(x))) * (eps * (0.5 + (eps * (eps * (-0.020833333333333332 + ((eps * eps) * (0.00026041666666666666 + ((eps * eps) * -1.5500992063492063e-6))))))))); end
code[x_, eps_] := N[(-2.0 * N[(N[(N[(N[Sin[x], $MachinePrecision] * N[Cos[N[(eps * 0.5), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] + N[(N[Sin[N[(eps * 0.5), $MachinePrecision]], $MachinePrecision] * N[Cos[x], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(eps * N[(0.5 + N[(eps * N[(eps * N[(-0.020833333333333332 + N[(N[(eps * eps), $MachinePrecision] * N[(0.00026041666666666666 + N[(N[(eps * eps), $MachinePrecision] * -1.5500992063492063e-6), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
-2 \cdot \left(\left(\sin x \cdot \cos \left(\varepsilon \cdot 0.5\right) + \sin \left(\varepsilon \cdot 0.5\right) \cdot \cos x\right) \cdot \left(\varepsilon \cdot \left(0.5 + \varepsilon \cdot \left(\varepsilon \cdot \left(-0.020833333333333332 + \left(\varepsilon \cdot \varepsilon\right) \cdot \left(0.00026041666666666666 + \left(\varepsilon \cdot \varepsilon\right) \cdot -1.5500992063492063 \cdot 10^{-6}\right)\right)\right)\right)\right)\right)
\end{array}
Initial program 49.1%
diff-cosN/A
*-commutativeN/A
*-lowering-*.f64N/A
Applied egg-rr99.5%
Taylor expanded in eps around inf
metadata-evalN/A
cancel-sign-sub-invN/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
cancel-sign-sub-invN/A
metadata-evalN/A
+-commutativeN/A
distribute-rgt-inN/A
*-commutativeN/A
associate-*l*N/A
metadata-evalN/A
*-rgt-identityN/A
*-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6499.5%
Simplified99.5%
Taylor expanded in eps around 0
*-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
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6499.5%
Simplified99.5%
+-commutativeN/A
sin-sumN/A
*-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
cos-lowering-cos.f64N/A
*-lowering-*.f6499.6%
Applied egg-rr99.6%
Final simplification99.6%
(FPCore (x eps)
:precision binary64
(*
(+
0.5
(*
(* eps eps)
(+
-0.020833333333333332
(*
(* eps eps)
(+ 0.00026041666666666666 (* (* eps eps) -1.5500992063492063e-6))))))
(* eps (* -2.0 (sin (+ x (* eps 0.5)))))))
double code(double x, double eps) {
return (0.5 + ((eps * eps) * (-0.020833333333333332 + ((eps * eps) * (0.00026041666666666666 + ((eps * eps) * -1.5500992063492063e-6)))))) * (eps * (-2.0 * sin((x + (eps * 0.5)))));
}
real(8) function code(x, eps)
real(8), intent (in) :: x
real(8), intent (in) :: eps
code = (0.5d0 + ((eps * eps) * ((-0.020833333333333332d0) + ((eps * eps) * (0.00026041666666666666d0 + ((eps * eps) * (-1.5500992063492063d-6))))))) * (eps * ((-2.0d0) * sin((x + (eps * 0.5d0)))))
end function
public static double code(double x, double eps) {
return (0.5 + ((eps * eps) * (-0.020833333333333332 + ((eps * eps) * (0.00026041666666666666 + ((eps * eps) * -1.5500992063492063e-6)))))) * (eps * (-2.0 * Math.sin((x + (eps * 0.5)))));
}
def code(x, eps): return (0.5 + ((eps * eps) * (-0.020833333333333332 + ((eps * eps) * (0.00026041666666666666 + ((eps * eps) * -1.5500992063492063e-6)))))) * (eps * (-2.0 * math.sin((x + (eps * 0.5)))))
function code(x, eps) return Float64(Float64(0.5 + Float64(Float64(eps * eps) * Float64(-0.020833333333333332 + Float64(Float64(eps * eps) * Float64(0.00026041666666666666 + Float64(Float64(eps * eps) * -1.5500992063492063e-6)))))) * Float64(eps * Float64(-2.0 * sin(Float64(x + Float64(eps * 0.5)))))) end
function tmp = code(x, eps) tmp = (0.5 + ((eps * eps) * (-0.020833333333333332 + ((eps * eps) * (0.00026041666666666666 + ((eps * eps) * -1.5500992063492063e-6)))))) * (eps * (-2.0 * sin((x + (eps * 0.5))))); end
code[x_, eps_] := N[(N[(0.5 + N[(N[(eps * eps), $MachinePrecision] * N[(-0.020833333333333332 + N[(N[(eps * eps), $MachinePrecision] * N[(0.00026041666666666666 + N[(N[(eps * eps), $MachinePrecision] * -1.5500992063492063e-6), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(eps * N[(-2.0 * N[Sin[N[(x + N[(eps * 0.5), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(0.5 + \left(\varepsilon \cdot \varepsilon\right) \cdot \left(-0.020833333333333332 + \left(\varepsilon \cdot \varepsilon\right) \cdot \left(0.00026041666666666666 + \left(\varepsilon \cdot \varepsilon\right) \cdot -1.5500992063492063 \cdot 10^{-6}\right)\right)\right) \cdot \left(\varepsilon \cdot \left(-2 \cdot \sin \left(x + \varepsilon \cdot 0.5\right)\right)\right)
\end{array}
Initial program 49.1%
diff-cosN/A
*-commutativeN/A
*-lowering-*.f64N/A
Applied egg-rr99.5%
Taylor expanded in eps around inf
metadata-evalN/A
cancel-sign-sub-invN/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
cancel-sign-sub-invN/A
metadata-evalN/A
+-commutativeN/A
distribute-rgt-inN/A
*-commutativeN/A
associate-*l*N/A
metadata-evalN/A
*-rgt-identityN/A
*-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6499.5%
Simplified99.5%
Taylor expanded in eps around 0
*-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
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6499.5%
Simplified99.5%
associate-*l*N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
Applied egg-rr99.5%
(FPCore (x eps)
:precision binary64
(*
-2.0
(*
(*
eps
(+
0.5
(*
eps
(*
eps
(+
-0.020833333333333332
(*
(* eps eps)
(+
0.00026041666666666666
(* (* eps eps) -1.5500992063492063e-6))))))))
(sin (+ x (* eps 0.5))))))
double code(double x, double eps) {
return -2.0 * ((eps * (0.5 + (eps * (eps * (-0.020833333333333332 + ((eps * eps) * (0.00026041666666666666 + ((eps * eps) * -1.5500992063492063e-6)))))))) * sin((x + (eps * 0.5))));
}
real(8) function code(x, eps)
real(8), intent (in) :: x
real(8), intent (in) :: eps
code = (-2.0d0) * ((eps * (0.5d0 + (eps * (eps * ((-0.020833333333333332d0) + ((eps * eps) * (0.00026041666666666666d0 + ((eps * eps) * (-1.5500992063492063d-6))))))))) * sin((x + (eps * 0.5d0))))
end function
public static double code(double x, double eps) {
return -2.0 * ((eps * (0.5 + (eps * (eps * (-0.020833333333333332 + ((eps * eps) * (0.00026041666666666666 + ((eps * eps) * -1.5500992063492063e-6)))))))) * Math.sin((x + (eps * 0.5))));
}
def code(x, eps): return -2.0 * ((eps * (0.5 + (eps * (eps * (-0.020833333333333332 + ((eps * eps) * (0.00026041666666666666 + ((eps * eps) * -1.5500992063492063e-6)))))))) * math.sin((x + (eps * 0.5))))
function code(x, eps) return Float64(-2.0 * Float64(Float64(eps * Float64(0.5 + Float64(eps * Float64(eps * Float64(-0.020833333333333332 + Float64(Float64(eps * eps) * Float64(0.00026041666666666666 + Float64(Float64(eps * eps) * -1.5500992063492063e-6)))))))) * sin(Float64(x + Float64(eps * 0.5))))) end
function tmp = code(x, eps) tmp = -2.0 * ((eps * (0.5 + (eps * (eps * (-0.020833333333333332 + ((eps * eps) * (0.00026041666666666666 + ((eps * eps) * -1.5500992063492063e-6)))))))) * sin((x + (eps * 0.5)))); end
code[x_, eps_] := N[(-2.0 * N[(N[(eps * N[(0.5 + N[(eps * N[(eps * N[(-0.020833333333333332 + N[(N[(eps * eps), $MachinePrecision] * N[(0.00026041666666666666 + N[(N[(eps * eps), $MachinePrecision] * -1.5500992063492063e-6), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[Sin[N[(x + N[(eps * 0.5), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
-2 \cdot \left(\left(\varepsilon \cdot \left(0.5 + \varepsilon \cdot \left(\varepsilon \cdot \left(-0.020833333333333332 + \left(\varepsilon \cdot \varepsilon\right) \cdot \left(0.00026041666666666666 + \left(\varepsilon \cdot \varepsilon\right) \cdot -1.5500992063492063 \cdot 10^{-6}\right)\right)\right)\right)\right) \cdot \sin \left(x + \varepsilon \cdot 0.5\right)\right)
\end{array}
Initial program 49.1%
diff-cosN/A
*-commutativeN/A
*-lowering-*.f64N/A
Applied egg-rr99.5%
Taylor expanded in eps around inf
metadata-evalN/A
cancel-sign-sub-invN/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
cancel-sign-sub-invN/A
metadata-evalN/A
+-commutativeN/A
distribute-rgt-inN/A
*-commutativeN/A
associate-*l*N/A
metadata-evalN/A
*-rgt-identityN/A
*-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6499.5%
Simplified99.5%
Taylor expanded in eps around 0
*-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
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6499.5%
Simplified99.5%
Final simplification99.5%
(FPCore (x eps)
:precision binary64
(*
-2.0
(*
(sin (+ x (* eps 0.5)))
(*
eps
(+
0.5
(*
eps
(*
eps
(+ -0.020833333333333332 (* (* eps eps) 0.00026041666666666666)))))))))
double code(double x, double eps) {
return -2.0 * (sin((x + (eps * 0.5))) * (eps * (0.5 + (eps * (eps * (-0.020833333333333332 + ((eps * eps) * 0.00026041666666666666)))))));
}
real(8) function code(x, eps)
real(8), intent (in) :: x
real(8), intent (in) :: eps
code = (-2.0d0) * (sin((x + (eps * 0.5d0))) * (eps * (0.5d0 + (eps * (eps * ((-0.020833333333333332d0) + ((eps * eps) * 0.00026041666666666666d0)))))))
end function
public static double code(double x, double eps) {
return -2.0 * (Math.sin((x + (eps * 0.5))) * (eps * (0.5 + (eps * (eps * (-0.020833333333333332 + ((eps * eps) * 0.00026041666666666666)))))));
}
def code(x, eps): return -2.0 * (math.sin((x + (eps * 0.5))) * (eps * (0.5 + (eps * (eps * (-0.020833333333333332 + ((eps * eps) * 0.00026041666666666666)))))))
function code(x, eps) return Float64(-2.0 * Float64(sin(Float64(x + Float64(eps * 0.5))) * Float64(eps * Float64(0.5 + Float64(eps * Float64(eps * Float64(-0.020833333333333332 + Float64(Float64(eps * eps) * 0.00026041666666666666)))))))) end
function tmp = code(x, eps) tmp = -2.0 * (sin((x + (eps * 0.5))) * (eps * (0.5 + (eps * (eps * (-0.020833333333333332 + ((eps * eps) * 0.00026041666666666666))))))); end
code[x_, eps_] := N[(-2.0 * N[(N[Sin[N[(x + N[(eps * 0.5), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * N[(eps * N[(0.5 + N[(eps * N[(eps * N[(-0.020833333333333332 + N[(N[(eps * eps), $MachinePrecision] * 0.00026041666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
-2 \cdot \left(\sin \left(x + \varepsilon \cdot 0.5\right) \cdot \left(\varepsilon \cdot \left(0.5 + \varepsilon \cdot \left(\varepsilon \cdot \left(-0.020833333333333332 + \left(\varepsilon \cdot \varepsilon\right) \cdot 0.00026041666666666666\right)\right)\right)\right)\right)
\end{array}
Initial program 49.1%
diff-cosN/A
*-commutativeN/A
*-lowering-*.f64N/A
Applied egg-rr99.5%
Taylor expanded in eps around inf
metadata-evalN/A
cancel-sign-sub-invN/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
cancel-sign-sub-invN/A
metadata-evalN/A
+-commutativeN/A
distribute-rgt-inN/A
*-commutativeN/A
associate-*l*N/A
metadata-evalN/A
*-rgt-identityN/A
*-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6499.5%
Simplified99.5%
Taylor expanded in eps around 0
*-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.4%
Simplified99.4%
Final simplification99.4%
(FPCore (x eps) :precision binary64 (* -2.0 (* (sin (+ x (* eps 0.5))) (* eps (+ 0.5 (* eps (* eps -0.020833333333333332)))))))
double code(double x, double eps) {
return -2.0 * (sin((x + (eps * 0.5))) * (eps * (0.5 + (eps * (eps * -0.020833333333333332)))));
}
real(8) function code(x, eps)
real(8), intent (in) :: x
real(8), intent (in) :: eps
code = (-2.0d0) * (sin((x + (eps * 0.5d0))) * (eps * (0.5d0 + (eps * (eps * (-0.020833333333333332d0))))))
end function
public static double code(double x, double eps) {
return -2.0 * (Math.sin((x + (eps * 0.5))) * (eps * (0.5 + (eps * (eps * -0.020833333333333332)))));
}
def code(x, eps): return -2.0 * (math.sin((x + (eps * 0.5))) * (eps * (0.5 + (eps * (eps * -0.020833333333333332)))))
function code(x, eps) return Float64(-2.0 * Float64(sin(Float64(x + Float64(eps * 0.5))) * Float64(eps * Float64(0.5 + Float64(eps * Float64(eps * -0.020833333333333332)))))) end
function tmp = code(x, eps) tmp = -2.0 * (sin((x + (eps * 0.5))) * (eps * (0.5 + (eps * (eps * -0.020833333333333332))))); end
code[x_, eps_] := N[(-2.0 * N[(N[Sin[N[(x + N[(eps * 0.5), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * N[(eps * N[(0.5 + N[(eps * N[(eps * -0.020833333333333332), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
-2 \cdot \left(\sin \left(x + \varepsilon \cdot 0.5\right) \cdot \left(\varepsilon \cdot \left(0.5 + \varepsilon \cdot \left(\varepsilon \cdot -0.020833333333333332\right)\right)\right)\right)
\end{array}
Initial program 49.1%
diff-cosN/A
*-commutativeN/A
*-lowering-*.f64N/A
Applied egg-rr99.5%
Taylor expanded in eps around inf
metadata-evalN/A
cancel-sign-sub-invN/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
cancel-sign-sub-invN/A
metadata-evalN/A
+-commutativeN/A
distribute-rgt-inN/A
*-commutativeN/A
associate-*l*N/A
metadata-evalN/A
*-rgt-identityN/A
*-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6499.5%
Simplified99.5%
Taylor expanded in eps around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6499.2%
Simplified99.2%
Final simplification99.2%
(FPCore (x eps) :precision binary64 (* -2.0 (* (* eps 0.5) (sin (+ x (* eps 0.5))))))
double code(double x, double eps) {
return -2.0 * ((eps * 0.5) * sin((x + (eps * 0.5))));
}
real(8) function code(x, eps)
real(8), intent (in) :: x
real(8), intent (in) :: eps
code = (-2.0d0) * ((eps * 0.5d0) * sin((x + (eps * 0.5d0))))
end function
public static double code(double x, double eps) {
return -2.0 * ((eps * 0.5) * Math.sin((x + (eps * 0.5))));
}
def code(x, eps): return -2.0 * ((eps * 0.5) * math.sin((x + (eps * 0.5))))
function code(x, eps) return Float64(-2.0 * Float64(Float64(eps * 0.5) * sin(Float64(x + Float64(eps * 0.5))))) end
function tmp = code(x, eps) tmp = -2.0 * ((eps * 0.5) * sin((x + (eps * 0.5)))); end
code[x_, eps_] := N[(-2.0 * N[(N[(eps * 0.5), $MachinePrecision] * N[Sin[N[(x + N[(eps * 0.5), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
-2 \cdot \left(\left(\varepsilon \cdot 0.5\right) \cdot \sin \left(x + \varepsilon \cdot 0.5\right)\right)
\end{array}
Initial program 49.1%
diff-cosN/A
*-commutativeN/A
*-lowering-*.f64N/A
Applied egg-rr99.5%
Taylor expanded in eps around inf
metadata-evalN/A
cancel-sign-sub-invN/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
cancel-sign-sub-invN/A
metadata-evalN/A
+-commutativeN/A
distribute-rgt-inN/A
*-commutativeN/A
associate-*l*N/A
metadata-evalN/A
*-rgt-identityN/A
*-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6499.5%
Simplified99.5%
Taylor expanded in eps around 0
*-commutativeN/A
*-lowering-*.f6498.8%
Simplified98.8%
Final simplification98.8%
(FPCore (x eps) :precision binary64 (* eps (- (* eps -0.5) (sin x))))
double code(double x, double eps) {
return eps * ((eps * -0.5) - sin(x));
}
real(8) function code(x, eps)
real(8), intent (in) :: x
real(8), intent (in) :: eps
code = eps * ((eps * (-0.5d0)) - sin(x))
end function
public static double code(double x, double eps) {
return eps * ((eps * -0.5) - Math.sin(x));
}
def code(x, eps): return eps * ((eps * -0.5) - math.sin(x))
function code(x, eps) return Float64(eps * Float64(Float64(eps * -0.5) - sin(x))) end
function tmp = code(x, eps) tmp = eps * ((eps * -0.5) - sin(x)); end
code[x_, eps_] := N[(eps * N[(N[(eps * -0.5), $MachinePrecision] - N[Sin[x], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\varepsilon \cdot \left(\varepsilon \cdot -0.5 - \sin x\right)
\end{array}
Initial program 49.1%
Taylor expanded in eps around 0
*-lowering-*.f64N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
sin-lowering-sin.f6498.8%
Simplified98.8%
Taylor expanded in x around 0
*-commutativeN/A
*-lowering-*.f6498.2%
Simplified98.2%
(FPCore (x eps) :precision binary64 (+ (* x (* x (* eps (+ (* x 0.16666666666666666) (* eps 0.25))))) (* eps (- (* eps -0.5) x))))
double code(double x, double eps) {
return (x * (x * (eps * ((x * 0.16666666666666666) + (eps * 0.25))))) + (eps * ((eps * -0.5) - x));
}
real(8) function code(x, eps)
real(8), intent (in) :: x
real(8), intent (in) :: eps
code = (x * (x * (eps * ((x * 0.16666666666666666d0) + (eps * 0.25d0))))) + (eps * ((eps * (-0.5d0)) - x))
end function
public static double code(double x, double eps) {
return (x * (x * (eps * ((x * 0.16666666666666666) + (eps * 0.25))))) + (eps * ((eps * -0.5) - x));
}
def code(x, eps): return (x * (x * (eps * ((x * 0.16666666666666666) + (eps * 0.25))))) + (eps * ((eps * -0.5) - x))
function code(x, eps) return Float64(Float64(x * Float64(x * Float64(eps * Float64(Float64(x * 0.16666666666666666) + Float64(eps * 0.25))))) + Float64(eps * Float64(Float64(eps * -0.5) - x))) end
function tmp = code(x, eps) tmp = (x * (x * (eps * ((x * 0.16666666666666666) + (eps * 0.25))))) + (eps * ((eps * -0.5) - x)); end
code[x_, eps_] := N[(N[(x * N[(x * N[(eps * N[(N[(x * 0.16666666666666666), $MachinePrecision] + N[(eps * 0.25), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(eps * N[(N[(eps * -0.5), $MachinePrecision] - x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x \cdot \left(x \cdot \left(\varepsilon \cdot \left(x \cdot 0.16666666666666666 + \varepsilon \cdot 0.25\right)\right)\right) + \varepsilon \cdot \left(\varepsilon \cdot -0.5 - x\right)
\end{array}
Initial program 49.1%
Taylor expanded in eps around 0
*-lowering-*.f64N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
sin-lowering-sin.f6498.8%
Simplified98.8%
Taylor expanded in x around 0
+-commutativeN/A
+-commutativeN/A
distribute-rgt-inN/A
associate-+l+N/A
+-lowering-+.f64N/A
Simplified96.9%
(FPCore (x eps) :precision binary64 (* eps (+ (* eps -0.5) (* x (+ -1.0 (* x (+ (* x 0.16666666666666666) (* eps 0.25))))))))
double code(double x, double eps) {
return eps * ((eps * -0.5) + (x * (-1.0 + (x * ((x * 0.16666666666666666) + (eps * 0.25))))));
}
real(8) function code(x, eps)
real(8), intent (in) :: x
real(8), intent (in) :: eps
code = eps * ((eps * (-0.5d0)) + (x * ((-1.0d0) + (x * ((x * 0.16666666666666666d0) + (eps * 0.25d0))))))
end function
public static double code(double x, double eps) {
return eps * ((eps * -0.5) + (x * (-1.0 + (x * ((x * 0.16666666666666666) + (eps * 0.25))))));
}
def code(x, eps): return eps * ((eps * -0.5) + (x * (-1.0 + (x * ((x * 0.16666666666666666) + (eps * 0.25))))))
function code(x, eps) return Float64(eps * Float64(Float64(eps * -0.5) + Float64(x * Float64(-1.0 + Float64(x * Float64(Float64(x * 0.16666666666666666) + Float64(eps * 0.25))))))) end
function tmp = code(x, eps) tmp = eps * ((eps * -0.5) + (x * (-1.0 + (x * ((x * 0.16666666666666666) + (eps * 0.25)))))); end
code[x_, eps_] := N[(eps * N[(N[(eps * -0.5), $MachinePrecision] + N[(x * N[(-1.0 + N[(x * N[(N[(x * 0.16666666666666666), $MachinePrecision] + N[(eps * 0.25), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\varepsilon \cdot \left(\varepsilon \cdot -0.5 + x \cdot \left(-1 + x \cdot \left(x \cdot 0.16666666666666666 + \varepsilon \cdot 0.25\right)\right)\right)
\end{array}
Initial program 49.1%
Taylor expanded in eps around 0
*-lowering-*.f64N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
sin-lowering-sin.f6498.8%
Simplified98.8%
Taylor expanded in x around 0
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6496.9%
Simplified96.9%
(FPCore (x eps) :precision binary64 (* eps (- (* eps (+ -0.5 (* eps (* x 0.16666666666666666)))) x)))
double code(double x, double eps) {
return eps * ((eps * (-0.5 + (eps * (x * 0.16666666666666666)))) - x);
}
real(8) function code(x, eps)
real(8), intent (in) :: x
real(8), intent (in) :: eps
code = eps * ((eps * ((-0.5d0) + (eps * (x * 0.16666666666666666d0)))) - x)
end function
public static double code(double x, double eps) {
return eps * ((eps * (-0.5 + (eps * (x * 0.16666666666666666)))) - x);
}
def code(x, eps): return eps * ((eps * (-0.5 + (eps * (x * 0.16666666666666666)))) - x)
function code(x, eps) return Float64(eps * Float64(Float64(eps * Float64(-0.5 + Float64(eps * Float64(x * 0.16666666666666666)))) - x)) end
function tmp = code(x, eps) tmp = eps * ((eps * (-0.5 + (eps * (x * 0.16666666666666666)))) - x); end
code[x_, eps_] := N[(eps * N[(N[(eps * N[(-0.5 + N[(eps * N[(x * 0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - x), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\varepsilon \cdot \left(\varepsilon \cdot \left(-0.5 + \varepsilon \cdot \left(x \cdot 0.16666666666666666\right)\right) - x\right)
\end{array}
Initial program 49.1%
Taylor expanded in x around 0
+-commutativeN/A
associate--l+N/A
+-commutativeN/A
mul-1-negN/A
sub-negN/A
--lowering--.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
cos-lowering-cos.f64N/A
*-lowering-*.f64N/A
sin-lowering-sin.f6475.6%
Simplified75.6%
Taylor expanded in eps around 0
*-lowering-*.f64N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
metadata-evalN/A
*-lowering-*.f64N/A
metadata-evalN/A
*-commutativeN/A
*-lowering-*.f6496.4%
Simplified96.4%
(FPCore (x eps) :precision binary64 (if (<= x -3.3e-147) (* 2.0 (* x x)) (* eps (* eps -0.5))))
double code(double x, double eps) {
double tmp;
if (x <= -3.3e-147) {
tmp = 2.0 * (x * x);
} else {
tmp = eps * (eps * -0.5);
}
return tmp;
}
real(8) function code(x, eps)
real(8), intent (in) :: x
real(8), intent (in) :: eps
real(8) :: tmp
if (x <= (-3.3d-147)) then
tmp = 2.0d0 * (x * x)
else
tmp = eps * (eps * (-0.5d0))
end if
code = tmp
end function
public static double code(double x, double eps) {
double tmp;
if (x <= -3.3e-147) {
tmp = 2.0 * (x * x);
} else {
tmp = eps * (eps * -0.5);
}
return tmp;
}
def code(x, eps): tmp = 0 if x <= -3.3e-147: tmp = 2.0 * (x * x) else: tmp = eps * (eps * -0.5) return tmp
function code(x, eps) tmp = 0.0 if (x <= -3.3e-147) tmp = Float64(2.0 * Float64(x * x)); else tmp = Float64(eps * Float64(eps * -0.5)); end return tmp end
function tmp_2 = code(x, eps) tmp = 0.0; if (x <= -3.3e-147) tmp = 2.0 * (x * x); else tmp = eps * (eps * -0.5); end tmp_2 = tmp; end
code[x_, eps_] := If[LessEqual[x, -3.3e-147], N[(2.0 * N[(x * x), $MachinePrecision]), $MachinePrecision], N[(eps * N[(eps * -0.5), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -3.3 \cdot 10^{-147}:\\
\;\;\;\;2 \cdot \left(x \cdot x\right)\\
\mathbf{else}:\\
\;\;\;\;\varepsilon \cdot \left(\varepsilon \cdot -0.5\right)\\
\end{array}
\end{array}
if x < -3.29999999999999987e-147Initial program 9.9%
flip--N/A
clear-numN/A
/-lowering-/.f64N/A
clear-numN/A
flip--N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
cos-lowering-cos.f64N/A
+-lowering-+.f64N/A
cos-lowering-cos.f649.9%
Applied egg-rr9.9%
Taylor expanded in eps around 0
div-subN/A
sub-negN/A
+-lowering-+.f64N/A
Simplified77.4%
Taylor expanded in x around 0
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6411.7%
Simplified11.7%
if -3.29999999999999987e-147 < x Initial program 63.6%
Taylor expanded in x around 0
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
cos-lowering-cos.f6462.7%
Simplified62.7%
Taylor expanded in eps around 0
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6465.0%
Simplified65.0%
(FPCore (x eps) :precision binary64 (* eps (- (* eps -0.5) x)))
double code(double x, double eps) {
return eps * ((eps * -0.5) - x);
}
real(8) function code(x, eps)
real(8), intent (in) :: x
real(8), intent (in) :: eps
code = eps * ((eps * (-0.5d0)) - x)
end function
public static double code(double x, double eps) {
return eps * ((eps * -0.5) - x);
}
def code(x, eps): return eps * ((eps * -0.5) - x)
function code(x, eps) return Float64(eps * Float64(Float64(eps * -0.5) - x)) end
function tmp = code(x, eps) tmp = eps * ((eps * -0.5) - x); end
code[x_, eps_] := N[(eps * N[(N[(eps * -0.5), $MachinePrecision] - x), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\varepsilon \cdot \left(\varepsilon \cdot -0.5 - x\right)
\end{array}
Initial program 49.1%
Taylor expanded in x around 0
+-commutativeN/A
associate--l+N/A
+-commutativeN/A
mul-1-negN/A
sub-negN/A
--lowering--.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
cos-lowering-cos.f64N/A
*-lowering-*.f64N/A
sin-lowering-sin.f6475.6%
Simplified75.6%
Taylor expanded in eps around 0
*-lowering-*.f64N/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6496.4%
Simplified96.4%
(FPCore (x eps) :precision binary64 (- 0.0 (* x eps)))
double code(double x, double eps) {
return 0.0 - (x * eps);
}
real(8) function code(x, eps)
real(8), intent (in) :: x
real(8), intent (in) :: eps
code = 0.0d0 - (x * eps)
end function
public static double code(double x, double eps) {
return 0.0 - (x * eps);
}
def code(x, eps): return 0.0 - (x * eps)
function code(x, eps) return Float64(0.0 - Float64(x * eps)) end
function tmp = code(x, eps) tmp = 0.0 - (x * eps); end
code[x_, eps_] := N[(0.0 - N[(x * eps), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
0 - x \cdot \varepsilon
\end{array}
Initial program 49.1%
Taylor expanded in x around 0
+-commutativeN/A
associate--l+N/A
+-commutativeN/A
mul-1-negN/A
sub-negN/A
--lowering--.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
cos-lowering-cos.f64N/A
*-lowering-*.f64N/A
sin-lowering-sin.f6475.6%
Simplified75.6%
Taylor expanded in eps around 0
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f6475.6%
Simplified75.6%
Final simplification75.6%
(FPCore (x eps) :precision binary64 (* 2.0 (* x x)))
double code(double x, double eps) {
return 2.0 * (x * x);
}
real(8) function code(x, eps)
real(8), intent (in) :: x
real(8), intent (in) :: eps
code = 2.0d0 * (x * x)
end function
public static double code(double x, double eps) {
return 2.0 * (x * x);
}
def code(x, eps): return 2.0 * (x * x)
function code(x, eps) return Float64(2.0 * Float64(x * x)) end
function tmp = code(x, eps) tmp = 2.0 * (x * x); end
code[x_, eps_] := N[(2.0 * N[(x * x), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
2 \cdot \left(x \cdot x\right)
\end{array}
Initial program 49.1%
flip--N/A
clear-numN/A
/-lowering-/.f64N/A
clear-numN/A
flip--N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
cos-lowering-cos.f64N/A
+-lowering-+.f64N/A
cos-lowering-cos.f6449.1%
Applied egg-rr49.1%
Taylor expanded in eps around 0
div-subN/A
sub-negN/A
+-lowering-+.f64N/A
Simplified65.6%
Taylor expanded in x around 0
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6447.9%
Simplified47.9%
(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 49.1%
Taylor expanded in x around 0
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
cos-lowering-cos.f6447.2%
Simplified47.2%
Taylor expanded in eps around 0
Simplified47.2%
metadata-eval47.2%
Applied egg-rr47.2%
(FPCore (x eps) :precision binary64 (* (* -2.0 (sin (+ x (/ eps 2.0)))) (sin (/ eps 2.0))))
double code(double x, double eps) {
return (-2.0 * sin((x + (eps / 2.0)))) * sin((eps / 2.0));
}
real(8) function code(x, eps)
real(8), intent (in) :: x
real(8), intent (in) :: eps
code = ((-2.0d0) * sin((x + (eps / 2.0d0)))) * sin((eps / 2.0d0))
end function
public static double code(double x, double eps) {
return (-2.0 * Math.sin((x + (eps / 2.0)))) * Math.sin((eps / 2.0));
}
def code(x, eps): return (-2.0 * math.sin((x + (eps / 2.0)))) * math.sin((eps / 2.0))
function code(x, eps) return Float64(Float64(-2.0 * sin(Float64(x + Float64(eps / 2.0)))) * sin(Float64(eps / 2.0))) end
function tmp = code(x, eps) tmp = (-2.0 * sin((x + (eps / 2.0)))) * sin((eps / 2.0)); end
code[x_, eps_] := N[(N[(-2.0 * N[Sin[N[(x + N[(eps / 2.0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] * N[Sin[N[(eps / 2.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(-2 \cdot \sin \left(x + \frac{\varepsilon}{2}\right)\right) \cdot \sin \left(\frac{\varepsilon}{2}\right)
\end{array}
herbie shell --seed 2024140
(FPCore (x eps)
:name "2cos (problem 3.3.5)"
:precision binary64
:pre (and (and (and (<= -10000.0 x) (<= x 10000.0)) (< (* 1e-16 (fabs x)) eps)) (< eps (fabs x)))
:alt
(! :herbie-platform default (* -2 (sin (+ x (/ eps 2))) (sin (/ eps 2))))
(- (cos (+ x eps)) (cos x)))