
(FPCore (N) :precision binary64 (- (atan (+ N 1.0)) (atan N)))
double code(double N) {
return atan((N + 1.0)) - atan(N);
}
real(8) function code(n)
real(8), intent (in) :: n
code = atan((n + 1.0d0)) - atan(n)
end function
public static double code(double N) {
return Math.atan((N + 1.0)) - Math.atan(N);
}
def code(N): return math.atan((N + 1.0)) - math.atan(N)
function code(N) return Float64(atan(Float64(N + 1.0)) - atan(N)) end
function tmp = code(N) tmp = atan((N + 1.0)) - atan(N); end
code[N_] := N[(N[ArcTan[N[(N + 1.0), $MachinePrecision]], $MachinePrecision] - N[ArcTan[N], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\tan^{-1} \left(N + 1\right) - \tan^{-1} N
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 5 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (N) :precision binary64 (- (atan (+ N 1.0)) (atan N)))
double code(double N) {
return atan((N + 1.0)) - atan(N);
}
real(8) function code(n)
real(8), intent (in) :: n
code = atan((n + 1.0d0)) - atan(n)
end function
public static double code(double N) {
return Math.atan((N + 1.0)) - Math.atan(N);
}
def code(N): return math.atan((N + 1.0)) - math.atan(N)
function code(N) return Float64(atan(Float64(N + 1.0)) - atan(N)) end
function tmp = code(N) tmp = atan((N + 1.0)) - atan(N); end
code[N_] := N[(N[ArcTan[N[(N + 1.0), $MachinePrecision]], $MachinePrecision] - N[ArcTan[N], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\tan^{-1} \left(N + 1\right) - \tan^{-1} N
\end{array}
(FPCore (N) :precision binary64 (atan2 1.0 (+ 1.0 (fma N N N))))
double code(double N) {
return atan2(1.0, (1.0 + fma(N, N, N)));
}
function code(N) return atan(1.0, Float64(1.0 + fma(N, N, N))) end
code[N_] := N[ArcTan[1.0 / N[(1.0 + N[(N * N + N), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
\\
\tan^{-1}_* \frac{1}{1 + \mathsf{fma}\left(N, N, N\right)}
\end{array}
Initial program 78.1%
diff-atan78.1%
associate--l+78.1%
+-commutative78.1%
*-commutative78.1%
fma-def78.1%
Applied egg-rr78.1%
+-commutative78.1%
associate-+l-99.0%
+-inverses99.0%
metadata-eval99.0%
+-commutative99.0%
Simplified99.0%
Taylor expanded in N around -inf 99.0%
fma-udef99.0%
+-commutative99.0%
*-commutative99.0%
cancel-sign-sub-inv99.0%
metadata-eval99.0%
*-lft-identity99.0%
+-commutative99.0%
distribute-lft1-in99.0%
fma-udef99.0%
Simplified99.0%
Final simplification99.0%
(FPCore (N) :precision binary64 (if (or (<= N -1.0) (not (<= N 1.0))) (atan2 1.0 (* N N)) (atan2 1.0 1.0)))
double code(double N) {
double tmp;
if ((N <= -1.0) || !(N <= 1.0)) {
tmp = atan2(1.0, (N * N));
} else {
tmp = atan2(1.0, 1.0);
}
return tmp;
}
real(8) function code(n)
real(8), intent (in) :: n
real(8) :: tmp
if ((n <= (-1.0d0)) .or. (.not. (n <= 1.0d0))) then
tmp = atan2(1.0d0, (n * n))
else
tmp = atan2(1.0d0, 1.0d0)
end if
code = tmp
end function
public static double code(double N) {
double tmp;
if ((N <= -1.0) || !(N <= 1.0)) {
tmp = Math.atan2(1.0, (N * N));
} else {
tmp = Math.atan2(1.0, 1.0);
}
return tmp;
}
def code(N): tmp = 0 if (N <= -1.0) or not (N <= 1.0): tmp = math.atan2(1.0, (N * N)) else: tmp = math.atan2(1.0, 1.0) return tmp
function code(N) tmp = 0.0 if ((N <= -1.0) || !(N <= 1.0)) tmp = atan(1.0, Float64(N * N)); else tmp = atan(1.0, 1.0); end return tmp end
function tmp_2 = code(N) tmp = 0.0; if ((N <= -1.0) || ~((N <= 1.0))) tmp = atan2(1.0, (N * N)); else tmp = atan2(1.0, 1.0); end tmp_2 = tmp; end
code[N_] := If[Or[LessEqual[N, -1.0], N[Not[LessEqual[N, 1.0]], $MachinePrecision]], N[ArcTan[1.0 / N[(N * N), $MachinePrecision]], $MachinePrecision], N[ArcTan[1.0 / 1.0], $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;N \leq -1 \lor \neg \left(N \leq 1\right):\\
\;\;\;\;\tan^{-1}_* \frac{1}{N \cdot N}\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1}_* \frac{1}{1}\\
\end{array}
\end{array}
if N < -1 or 1 < N Initial program 49.5%
diff-atan49.5%
associate--l+49.5%
+-commutative49.5%
*-commutative49.5%
fma-def49.5%
Applied egg-rr49.5%
+-commutative49.5%
associate-+l-97.7%
+-inverses97.7%
metadata-eval97.7%
+-commutative97.7%
Simplified97.7%
Taylor expanded in N around -inf 97.7%
fma-udef97.7%
+-commutative97.7%
*-commutative97.7%
cancel-sign-sub-inv97.7%
metadata-eval97.7%
*-lft-identity97.7%
+-commutative97.7%
distribute-lft1-in97.7%
fma-udef97.7%
Simplified97.7%
Taylor expanded in N around inf 97.7%
unpow297.7%
Simplified97.7%
if -1 < N < 1Initial program 100.0%
diff-atan100.0%
associate--l+100.0%
+-commutative100.0%
*-commutative100.0%
fma-def100.0%
Applied egg-rr100.0%
+-commutative100.0%
associate-+l-100.0%
+-inverses100.0%
metadata-eval100.0%
+-commutative100.0%
Simplified100.0%
Taylor expanded in N around -inf 100.0%
fma-udef100.0%
+-commutative100.0%
*-commutative100.0%
cancel-sign-sub-inv100.0%
metadata-eval100.0%
*-lft-identity100.0%
+-commutative100.0%
distribute-lft1-in100.0%
fma-udef100.0%
Simplified100.0%
Taylor expanded in N around 0 97.2%
Final simplification97.4%
(FPCore (N) :precision binary64 (if (or (<= N -0.62) (not (<= N 1.65))) (atan2 1.0 (* N N)) (atan2 1.0 (+ 1.0 N))))
double code(double N) {
double tmp;
if ((N <= -0.62) || !(N <= 1.65)) {
tmp = atan2(1.0, (N * N));
} else {
tmp = atan2(1.0, (1.0 + N));
}
return tmp;
}
real(8) function code(n)
real(8), intent (in) :: n
real(8) :: tmp
if ((n <= (-0.62d0)) .or. (.not. (n <= 1.65d0))) then
tmp = atan2(1.0d0, (n * n))
else
tmp = atan2(1.0d0, (1.0d0 + n))
end if
code = tmp
end function
public static double code(double N) {
double tmp;
if ((N <= -0.62) || !(N <= 1.65)) {
tmp = Math.atan2(1.0, (N * N));
} else {
tmp = Math.atan2(1.0, (1.0 + N));
}
return tmp;
}
def code(N): tmp = 0 if (N <= -0.62) or not (N <= 1.65): tmp = math.atan2(1.0, (N * N)) else: tmp = math.atan2(1.0, (1.0 + N)) return tmp
function code(N) tmp = 0.0 if ((N <= -0.62) || !(N <= 1.65)) tmp = atan(1.0, Float64(N * N)); else tmp = atan(1.0, Float64(1.0 + N)); end return tmp end
function tmp_2 = code(N) tmp = 0.0; if ((N <= -0.62) || ~((N <= 1.65))) tmp = atan2(1.0, (N * N)); else tmp = atan2(1.0, (1.0 + N)); end tmp_2 = tmp; end
code[N_] := If[Or[LessEqual[N, -0.62], N[Not[LessEqual[N, 1.65]], $MachinePrecision]], N[ArcTan[1.0 / N[(N * N), $MachinePrecision]], $MachinePrecision], N[ArcTan[1.0 / N[(1.0 + N), $MachinePrecision]], $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;N \leq -0.62 \lor \neg \left(N \leq 1.65\right):\\
\;\;\;\;\tan^{-1}_* \frac{1}{N \cdot N}\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1}_* \frac{1}{1 + N}\\
\end{array}
\end{array}
if N < -0.619999999999999996 or 1.6499999999999999 < N Initial program 49.5%
diff-atan49.5%
associate--l+49.5%
+-commutative49.5%
*-commutative49.5%
fma-def49.5%
Applied egg-rr49.5%
+-commutative49.5%
associate-+l-97.7%
+-inverses97.7%
metadata-eval97.7%
+-commutative97.7%
Simplified97.7%
Taylor expanded in N around -inf 97.7%
fma-udef97.7%
+-commutative97.7%
*-commutative97.7%
cancel-sign-sub-inv97.7%
metadata-eval97.7%
*-lft-identity97.7%
+-commutative97.7%
distribute-lft1-in97.7%
fma-udef97.7%
Simplified97.7%
Taylor expanded in N around inf 97.7%
unpow297.7%
Simplified97.7%
if -0.619999999999999996 < N < 1.6499999999999999Initial program 100.0%
diff-atan100.0%
associate--l+100.0%
+-commutative100.0%
*-commutative100.0%
fma-def100.0%
Applied egg-rr100.0%
+-commutative100.0%
associate-+l-100.0%
+-inverses100.0%
metadata-eval100.0%
+-commutative100.0%
Simplified100.0%
Taylor expanded in N around -inf 100.0%
fma-udef100.0%
+-commutative100.0%
*-commutative100.0%
cancel-sign-sub-inv100.0%
metadata-eval100.0%
*-lft-identity100.0%
+-commutative100.0%
distribute-lft1-in100.0%
fma-udef100.0%
Simplified100.0%
Taylor expanded in N around 0 98.6%
Final simplification98.2%
(FPCore (N) :precision binary64 (atan2 1.0 (+ N (+ 1.0 (* N N)))))
double code(double N) {
return atan2(1.0, (N + (1.0 + (N * N))));
}
real(8) function code(n)
real(8), intent (in) :: n
code = atan2(1.0d0, (n + (1.0d0 + (n * n))))
end function
public static double code(double N) {
return Math.atan2(1.0, (N + (1.0 + (N * N))));
}
def code(N): return math.atan2(1.0, (N + (1.0 + (N * N))))
function code(N) return atan(1.0, Float64(N + Float64(1.0 + Float64(N * N)))) end
function tmp = code(N) tmp = atan2(1.0, (N + (1.0 + (N * N)))); end
code[N_] := N[ArcTan[1.0 / N[(N + N[(1.0 + N[(N * N), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
\\
\tan^{-1}_* \frac{1}{N + \left(1 + N \cdot N\right)}
\end{array}
Initial program 78.1%
diff-atan78.1%
associate--l+78.1%
+-commutative78.1%
*-commutative78.1%
fma-def78.1%
Applied egg-rr78.1%
+-commutative78.1%
associate-+l-99.0%
+-inverses99.0%
metadata-eval99.0%
+-commutative99.0%
Simplified99.0%
Taylor expanded in N around -inf 99.0%
fma-udef99.0%
+-commutative99.0%
*-commutative99.0%
cancel-sign-sub-inv99.0%
metadata-eval99.0%
*-lft-identity99.0%
+-commutative99.0%
distribute-lft1-in99.0%
fma-udef99.0%
Simplified99.0%
+-commutative99.0%
fma-udef99.0%
*-un-lft-identity99.0%
distribute-rgt-in99.0%
+-commutative99.0%
add-sqr-sqrt52.1%
associate-*l*52.1%
fma-def52.1%
Applied egg-rr52.1%
fma-udef52.1%
associate-*r*52.1%
add-sqr-sqrt99.0%
+-commutative99.0%
distribute-rgt-out99.0%
*-un-lft-identity99.0%
associate-+r+99.0%
+-commutative99.0%
associate-+r+99.0%
Applied egg-rr99.0%
Final simplification99.0%
(FPCore (N) :precision binary64 (atan2 1.0 1.0))
double code(double N) {
return atan2(1.0, 1.0);
}
real(8) function code(n)
real(8), intent (in) :: n
code = atan2(1.0d0, 1.0d0)
end function
public static double code(double N) {
return Math.atan2(1.0, 1.0);
}
def code(N): return math.atan2(1.0, 1.0)
function code(N) return atan(1.0, 1.0) end
function tmp = code(N) tmp = atan2(1.0, 1.0); end
code[N_] := N[ArcTan[1.0 / 1.0], $MachinePrecision]
\begin{array}{l}
\\
\tan^{-1}_* \frac{1}{1}
\end{array}
Initial program 78.1%
diff-atan78.1%
associate--l+78.1%
+-commutative78.1%
*-commutative78.1%
fma-def78.1%
Applied egg-rr78.1%
+-commutative78.1%
associate-+l-99.0%
+-inverses99.0%
metadata-eval99.0%
+-commutative99.0%
Simplified99.0%
Taylor expanded in N around -inf 99.0%
fma-udef99.0%
+-commutative99.0%
*-commutative99.0%
cancel-sign-sub-inv99.0%
metadata-eval99.0%
*-lft-identity99.0%
+-commutative99.0%
distribute-lft1-in99.0%
fma-udef99.0%
Simplified99.0%
Taylor expanded in N around 0 56.8%
Final simplification56.8%
(FPCore (N) :precision binary64 (atan (/ 1.0 (+ 1.0 (* N (+ N 1.0))))))
double code(double N) {
return atan((1.0 / (1.0 + (N * (N + 1.0)))));
}
real(8) function code(n)
real(8), intent (in) :: n
code = atan((1.0d0 / (1.0d0 + (n * (n + 1.0d0)))))
end function
public static double code(double N) {
return Math.atan((1.0 / (1.0 + (N * (N + 1.0)))));
}
def code(N): return math.atan((1.0 / (1.0 + (N * (N + 1.0)))))
function code(N) return atan(Float64(1.0 / Float64(1.0 + Float64(N * Float64(N + 1.0))))) end
function tmp = code(N) tmp = atan((1.0 / (1.0 + (N * (N + 1.0))))); end
code[N_] := N[ArcTan[N[(1.0 / N[(1.0 + N[(N * N[(N + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
\\
\tan^{-1} \left(\frac{1}{1 + N \cdot \left(N + 1\right)}\right)
\end{array}
herbie shell --seed 2023277
(FPCore (N)
:name "2atan (example 3.5)"
:precision binary64
:herbie-target
(atan (/ 1.0 (+ 1.0 (* N (+ N 1.0)))))
(- (atan (+ N 1.0)) (atan N)))