
(FPCore (lambda1 lambda2 phi1 phi2) :precision binary64 (atan2 (* (sin (- lambda1 lambda2)) (cos phi2)) (- (* (cos phi1) (sin phi2)) (* (* (sin phi1) (cos phi2)) (cos (- lambda1 lambda2))))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
return atan2((sin((lambda1 - lambda2)) * cos(phi2)), ((cos(phi1) * sin(phi2)) - ((sin(phi1) * cos(phi2)) * cos((lambda1 - lambda2)))));
}
real(8) function code(lambda1, lambda2, phi1, phi2)
real(8), intent (in) :: lambda1
real(8), intent (in) :: lambda2
real(8), intent (in) :: phi1
real(8), intent (in) :: phi2
code = atan2((sin((lambda1 - lambda2)) * cos(phi2)), ((cos(phi1) * sin(phi2)) - ((sin(phi1) * cos(phi2)) * cos((lambda1 - lambda2)))))
end function
public static double code(double lambda1, double lambda2, double phi1, double phi2) {
return Math.atan2((Math.sin((lambda1 - lambda2)) * Math.cos(phi2)), ((Math.cos(phi1) * Math.sin(phi2)) - ((Math.sin(phi1) * Math.cos(phi2)) * Math.cos((lambda1 - lambda2)))));
}
def code(lambda1, lambda2, phi1, phi2): return math.atan2((math.sin((lambda1 - lambda2)) * math.cos(phi2)), ((math.cos(phi1) * math.sin(phi2)) - ((math.sin(phi1) * math.cos(phi2)) * math.cos((lambda1 - lambda2)))))
function code(lambda1, lambda2, phi1, phi2) return atan(Float64(sin(Float64(lambda1 - lambda2)) * cos(phi2)), Float64(Float64(cos(phi1) * sin(phi2)) - Float64(Float64(sin(phi1) * cos(phi2)) * cos(Float64(lambda1 - lambda2))))) end
function tmp = code(lambda1, lambda2, phi1, phi2) tmp = atan2((sin((lambda1 - lambda2)) * cos(phi2)), ((cos(phi1) * sin(phi2)) - ((sin(phi1) * cos(phi2)) * cos((lambda1 - lambda2))))); end
code[lambda1_, lambda2_, phi1_, phi2_] := N[ArcTan[N[(N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision] / N[(N[(N[Cos[phi1], $MachinePrecision] * N[Sin[phi2], $MachinePrecision]), $MachinePrecision] - N[(N[(N[Sin[phi1], $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision] * N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
\\
\tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2}{\cos \phi_1 \cdot \sin \phi_2 - \left(\sin \phi_1 \cdot \cos \phi_2\right) \cdot \cos \left(\lambda_1 - \lambda_2\right)}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 36 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (lambda1 lambda2 phi1 phi2) :precision binary64 (atan2 (* (sin (- lambda1 lambda2)) (cos phi2)) (- (* (cos phi1) (sin phi2)) (* (* (sin phi1) (cos phi2)) (cos (- lambda1 lambda2))))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
return atan2((sin((lambda1 - lambda2)) * cos(phi2)), ((cos(phi1) * sin(phi2)) - ((sin(phi1) * cos(phi2)) * cos((lambda1 - lambda2)))));
}
real(8) function code(lambda1, lambda2, phi1, phi2)
real(8), intent (in) :: lambda1
real(8), intent (in) :: lambda2
real(8), intent (in) :: phi1
real(8), intent (in) :: phi2
code = atan2((sin((lambda1 - lambda2)) * cos(phi2)), ((cos(phi1) * sin(phi2)) - ((sin(phi1) * cos(phi2)) * cos((lambda1 - lambda2)))))
end function
public static double code(double lambda1, double lambda2, double phi1, double phi2) {
return Math.atan2((Math.sin((lambda1 - lambda2)) * Math.cos(phi2)), ((Math.cos(phi1) * Math.sin(phi2)) - ((Math.sin(phi1) * Math.cos(phi2)) * Math.cos((lambda1 - lambda2)))));
}
def code(lambda1, lambda2, phi1, phi2): return math.atan2((math.sin((lambda1 - lambda2)) * math.cos(phi2)), ((math.cos(phi1) * math.sin(phi2)) - ((math.sin(phi1) * math.cos(phi2)) * math.cos((lambda1 - lambda2)))))
function code(lambda1, lambda2, phi1, phi2) return atan(Float64(sin(Float64(lambda1 - lambda2)) * cos(phi2)), Float64(Float64(cos(phi1) * sin(phi2)) - Float64(Float64(sin(phi1) * cos(phi2)) * cos(Float64(lambda1 - lambda2))))) end
function tmp = code(lambda1, lambda2, phi1, phi2) tmp = atan2((sin((lambda1 - lambda2)) * cos(phi2)), ((cos(phi1) * sin(phi2)) - ((sin(phi1) * cos(phi2)) * cos((lambda1 - lambda2))))); end
code[lambda1_, lambda2_, phi1_, phi2_] := N[ArcTan[N[(N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision] / N[(N[(N[Cos[phi1], $MachinePrecision] * N[Sin[phi2], $MachinePrecision]), $MachinePrecision] - N[(N[(N[Sin[phi1], $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision] * N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
\\
\tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2}{\cos \phi_1 \cdot \sin \phi_2 - \left(\sin \phi_1 \cdot \cos \phi_2\right) \cdot \cos \left(\lambda_1 - \lambda_2\right)}
\end{array}
(FPCore (lambda1 lambda2 phi1 phi2)
:precision binary64
(atan2
(*
(fma (- (sin lambda2)) (cos lambda1) (* (cos lambda2) (sin lambda1)))
(cos phi2))
(-
(* (cos phi1) (sin phi2))
(*
(* (cos phi2) (sin phi1))
(fma (cos lambda2) (cos lambda1) (* (sin lambda2) (sin lambda1)))))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
return atan2((fma(-sin(lambda2), cos(lambda1), (cos(lambda2) * sin(lambda1))) * cos(phi2)), ((cos(phi1) * sin(phi2)) - ((cos(phi2) * sin(phi1)) * fma(cos(lambda2), cos(lambda1), (sin(lambda2) * sin(lambda1))))));
}
function code(lambda1, lambda2, phi1, phi2) return atan(Float64(fma(Float64(-sin(lambda2)), cos(lambda1), Float64(cos(lambda2) * sin(lambda1))) * cos(phi2)), Float64(Float64(cos(phi1) * sin(phi2)) - Float64(Float64(cos(phi2) * sin(phi1)) * fma(cos(lambda2), cos(lambda1), Float64(sin(lambda2) * sin(lambda1)))))) end
code[lambda1_, lambda2_, phi1_, phi2_] := N[ArcTan[N[(N[((-N[Sin[lambda2], $MachinePrecision]) * N[Cos[lambda1], $MachinePrecision] + N[(N[Cos[lambda2], $MachinePrecision] * N[Sin[lambda1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision] / N[(N[(N[Cos[phi1], $MachinePrecision] * N[Sin[phi2], $MachinePrecision]), $MachinePrecision] - N[(N[(N[Cos[phi2], $MachinePrecision] * N[Sin[phi1], $MachinePrecision]), $MachinePrecision] * N[(N[Cos[lambda2], $MachinePrecision] * N[Cos[lambda1], $MachinePrecision] + N[(N[Sin[lambda2], $MachinePrecision] * N[Sin[lambda1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
\\
\tan^{-1}_* \frac{\mathsf{fma}\left(-\sin \lambda_2, \cos \lambda_1, \cos \lambda_2 \cdot \sin \lambda_1\right) \cdot \cos \phi_2}{\cos \phi_1 \cdot \sin \phi_2 - \left(\cos \phi_2 \cdot \sin \phi_1\right) \cdot \mathsf{fma}\left(\cos \lambda_2, \cos \lambda_1, \sin \lambda_2 \cdot \sin \lambda_1\right)}
\end{array}
Initial program 80.5%
sin-diff89.1%
sub-neg89.1%
Applied egg-rr89.1%
+-commutative89.1%
distribute-rgt-neg-in89.1%
sin-neg89.1%
*-commutative89.1%
fma-define89.2%
sin-neg89.2%
cos-neg89.2%
*-commutative89.2%
cos-neg89.2%
Simplified89.2%
cos-diff99.8%
*-commutative99.8%
Applied egg-rr99.8%
fma-define99.8%
Simplified99.8%
Final simplification99.8%
(FPCore (lambda1 lambda2 phi1 phi2)
:precision binary64
(atan2
(*
(fma (- (sin lambda2)) (cos lambda1) (* (cos lambda2) (sin lambda1)))
(cos phi2))
(-
(* (cos phi1) (sin phi2))
(*
(* (cos phi2) (sin phi1))
(+ (* (sin lambda2) (sin lambda1)) (* (cos lambda1) (cos lambda2)))))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
return atan2((fma(-sin(lambda2), cos(lambda1), (cos(lambda2) * sin(lambda1))) * cos(phi2)), ((cos(phi1) * sin(phi2)) - ((cos(phi2) * sin(phi1)) * ((sin(lambda2) * sin(lambda1)) + (cos(lambda1) * cos(lambda2))))));
}
function code(lambda1, lambda2, phi1, phi2) return atan(Float64(fma(Float64(-sin(lambda2)), cos(lambda1), Float64(cos(lambda2) * sin(lambda1))) * cos(phi2)), Float64(Float64(cos(phi1) * sin(phi2)) - Float64(Float64(cos(phi2) * sin(phi1)) * Float64(Float64(sin(lambda2) * sin(lambda1)) + Float64(cos(lambda1) * cos(lambda2)))))) end
code[lambda1_, lambda2_, phi1_, phi2_] := N[ArcTan[N[(N[((-N[Sin[lambda2], $MachinePrecision]) * N[Cos[lambda1], $MachinePrecision] + N[(N[Cos[lambda2], $MachinePrecision] * N[Sin[lambda1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision] / N[(N[(N[Cos[phi1], $MachinePrecision] * N[Sin[phi2], $MachinePrecision]), $MachinePrecision] - N[(N[(N[Cos[phi2], $MachinePrecision] * N[Sin[phi1], $MachinePrecision]), $MachinePrecision] * N[(N[(N[Sin[lambda2], $MachinePrecision] * N[Sin[lambda1], $MachinePrecision]), $MachinePrecision] + N[(N[Cos[lambda1], $MachinePrecision] * N[Cos[lambda2], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
\\
\tan^{-1}_* \frac{\mathsf{fma}\left(-\sin \lambda_2, \cos \lambda_1, \cos \lambda_2 \cdot \sin \lambda_1\right) \cdot \cos \phi_2}{\cos \phi_1 \cdot \sin \phi_2 - \left(\cos \phi_2 \cdot \sin \phi_1\right) \cdot \left(\sin \lambda_2 \cdot \sin \lambda_1 + \cos \lambda_1 \cdot \cos \lambda_2\right)}
\end{array}
Initial program 80.5%
sin-diff89.1%
sub-neg89.1%
Applied egg-rr89.1%
+-commutative89.1%
distribute-rgt-neg-in89.1%
sin-neg89.1%
*-commutative89.1%
fma-define89.2%
sin-neg89.2%
cos-neg89.2%
*-commutative89.2%
cos-neg89.2%
Simplified89.2%
cos-diff99.8%
+-commutative99.8%
*-commutative99.8%
Applied egg-rr99.8%
Final simplification99.8%
(FPCore (lambda1 lambda2 phi1 phi2)
:precision binary64
(atan2
(*
(cos phi2)
(- (* (cos lambda2) (sin lambda1)) (* (sin lambda2) (cos lambda1))))
(-
(* (cos phi1) (sin phi2))
(*
(* (cos phi2) (sin phi1))
(fma (cos lambda2) (cos lambda1) (* (sin lambda2) (sin lambda1)))))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
return atan2((cos(phi2) * ((cos(lambda2) * sin(lambda1)) - (sin(lambda2) * cos(lambda1)))), ((cos(phi1) * sin(phi2)) - ((cos(phi2) * sin(phi1)) * fma(cos(lambda2), cos(lambda1), (sin(lambda2) * sin(lambda1))))));
}
function code(lambda1, lambda2, phi1, phi2) return atan(Float64(cos(phi2) * Float64(Float64(cos(lambda2) * sin(lambda1)) - Float64(sin(lambda2) * cos(lambda1)))), Float64(Float64(cos(phi1) * sin(phi2)) - Float64(Float64(cos(phi2) * sin(phi1)) * fma(cos(lambda2), cos(lambda1), Float64(sin(lambda2) * sin(lambda1)))))) end
code[lambda1_, lambda2_, phi1_, phi2_] := N[ArcTan[N[(N[Cos[phi2], $MachinePrecision] * N[(N[(N[Cos[lambda2], $MachinePrecision] * N[Sin[lambda1], $MachinePrecision]), $MachinePrecision] - N[(N[Sin[lambda2], $MachinePrecision] * N[Cos[lambda1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(N[(N[Cos[phi1], $MachinePrecision] * N[Sin[phi2], $MachinePrecision]), $MachinePrecision] - N[(N[(N[Cos[phi2], $MachinePrecision] * N[Sin[phi1], $MachinePrecision]), $MachinePrecision] * N[(N[Cos[lambda2], $MachinePrecision] * N[Cos[lambda1], $MachinePrecision] + N[(N[Sin[lambda2], $MachinePrecision] * N[Sin[lambda1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
\\
\tan^{-1}_* \frac{\cos \phi_2 \cdot \left(\cos \lambda_2 \cdot \sin \lambda_1 - \sin \lambda_2 \cdot \cos \lambda_1\right)}{\cos \phi_1 \cdot \sin \phi_2 - \left(\cos \phi_2 \cdot \sin \phi_1\right) \cdot \mathsf{fma}\left(\cos \lambda_2, \cos \lambda_1, \sin \lambda_2 \cdot \sin \lambda_1\right)}
\end{array}
Initial program 80.5%
sin-diff89.1%
sub-neg89.1%
Applied egg-rr89.1%
+-commutative89.1%
distribute-rgt-neg-in89.1%
sin-neg89.1%
*-commutative89.1%
fma-define89.2%
sin-neg89.2%
cos-neg89.2%
*-commutative89.2%
cos-neg89.2%
Simplified89.2%
cos-diff99.8%
*-commutative99.8%
Applied egg-rr99.8%
fma-define99.8%
Simplified99.8%
Taylor expanded in lambda2 around inf 99.8%
neg-mul-199.8%
+-commutative99.8%
sub-neg99.8%
cos-neg99.8%
*-commutative99.8%
cos-neg99.8%
Simplified99.8%
Final simplification99.8%
(FPCore (lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (* (cos lambda2) (sin lambda1))) (t_1 (* (cos phi2) (sin phi1))))
(if (or (<= phi2 -9.6e-5) (not (<= phi2 8e-141)))
(atan2
(* (fma (- (sin lambda2)) (cos lambda1) t_0) (cos phi2))
(- (* (cos phi1) (sin phi2)) (* t_1 (cos (- lambda1 lambda2)))))
(atan2
(* (cos phi2) (- t_0 (* (sin lambda2) (cos lambda1))))
(-
(* phi2 (cos phi1))
(*
t_1
(fma (cos lambda2) (cos lambda1) (* (sin lambda2) (sin lambda1)))))))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = cos(lambda2) * sin(lambda1);
double t_1 = cos(phi2) * sin(phi1);
double tmp;
if ((phi2 <= -9.6e-5) || !(phi2 <= 8e-141)) {
tmp = atan2((fma(-sin(lambda2), cos(lambda1), t_0) * cos(phi2)), ((cos(phi1) * sin(phi2)) - (t_1 * cos((lambda1 - lambda2)))));
} else {
tmp = atan2((cos(phi2) * (t_0 - (sin(lambda2) * cos(lambda1)))), ((phi2 * cos(phi1)) - (t_1 * fma(cos(lambda2), cos(lambda1), (sin(lambda2) * sin(lambda1))))));
}
return tmp;
}
function code(lambda1, lambda2, phi1, phi2) t_0 = Float64(cos(lambda2) * sin(lambda1)) t_1 = Float64(cos(phi2) * sin(phi1)) tmp = 0.0 if ((phi2 <= -9.6e-5) || !(phi2 <= 8e-141)) tmp = atan(Float64(fma(Float64(-sin(lambda2)), cos(lambda1), t_0) * cos(phi2)), Float64(Float64(cos(phi1) * sin(phi2)) - Float64(t_1 * cos(Float64(lambda1 - lambda2))))); else tmp = atan(Float64(cos(phi2) * Float64(t_0 - Float64(sin(lambda2) * cos(lambda1)))), Float64(Float64(phi2 * cos(phi1)) - Float64(t_1 * fma(cos(lambda2), cos(lambda1), Float64(sin(lambda2) * sin(lambda1)))))); end return tmp end
code[lambda1_, lambda2_, phi1_, phi2_] := Block[{t$95$0 = N[(N[Cos[lambda2], $MachinePrecision] * N[Sin[lambda1], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[Cos[phi2], $MachinePrecision] * N[Sin[phi1], $MachinePrecision]), $MachinePrecision]}, If[Or[LessEqual[phi2, -9.6e-5], N[Not[LessEqual[phi2, 8e-141]], $MachinePrecision]], N[ArcTan[N[(N[((-N[Sin[lambda2], $MachinePrecision]) * N[Cos[lambda1], $MachinePrecision] + t$95$0), $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision] / N[(N[(N[Cos[phi1], $MachinePrecision] * N[Sin[phi2], $MachinePrecision]), $MachinePrecision] - N[(t$95$1 * N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], N[ArcTan[N[(N[Cos[phi2], $MachinePrecision] * N[(t$95$0 - N[(N[Sin[lambda2], $MachinePrecision] * N[Cos[lambda1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(N[(phi2 * N[Cos[phi1], $MachinePrecision]), $MachinePrecision] - N[(t$95$1 * N[(N[Cos[lambda2], $MachinePrecision] * N[Cos[lambda1], $MachinePrecision] + N[(N[Sin[lambda2], $MachinePrecision] * N[Sin[lambda1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \cos \lambda_2 \cdot \sin \lambda_1\\
t_1 := \cos \phi_2 \cdot \sin \phi_1\\
\mathbf{if}\;\phi_2 \leq -9.6 \cdot 10^{-5} \lor \neg \left(\phi_2 \leq 8 \cdot 10^{-141}\right):\\
\;\;\;\;\tan^{-1}_* \frac{\mathsf{fma}\left(-\sin \lambda_2, \cos \lambda_1, t\_0\right) \cdot \cos \phi_2}{\cos \phi_1 \cdot \sin \phi_2 - t\_1 \cdot \cos \left(\lambda_1 - \lambda_2\right)}\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1}_* \frac{\cos \phi_2 \cdot \left(t\_0 - \sin \lambda_2 \cdot \cos \lambda_1\right)}{\phi_2 \cdot \cos \phi_1 - t\_1 \cdot \mathsf{fma}\left(\cos \lambda_2, \cos \lambda_1, \sin \lambda_2 \cdot \sin \lambda_1\right)}\\
\end{array}
\end{array}
if phi2 < -9.6000000000000002e-5 or 8.0000000000000003e-141 < phi2 Initial program 76.6%
sin-diff87.7%
sub-neg87.7%
Applied egg-rr87.7%
+-commutative87.7%
distribute-rgt-neg-in87.7%
sin-neg87.7%
*-commutative87.7%
fma-define87.8%
sin-neg87.8%
cos-neg87.8%
*-commutative87.8%
cos-neg87.8%
Simplified87.8%
if -9.6000000000000002e-5 < phi2 < 8.0000000000000003e-141Initial program 86.1%
sin-diff91.1%
sub-neg91.1%
Applied egg-rr91.1%
+-commutative91.1%
distribute-rgt-neg-in91.1%
sin-neg91.1%
*-commutative91.1%
fma-define91.1%
sin-neg91.1%
cos-neg91.1%
*-commutative91.1%
cos-neg91.1%
Simplified91.1%
cos-diff99.9%
*-commutative99.9%
Applied egg-rr99.9%
fma-define99.9%
Simplified99.9%
Taylor expanded in lambda2 around inf 99.9%
neg-mul-199.9%
+-commutative99.9%
sub-neg99.9%
cos-neg99.9%
*-commutative99.9%
cos-neg99.9%
Simplified99.9%
Taylor expanded in phi2 around 0 99.9%
Final simplification92.8%
(FPCore (lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (* (cos phi1) (sin phi2)))
(t_1
(*
(fma (- (sin lambda2)) (cos lambda1) (* (cos lambda2) (sin lambda1)))
(cos phi2))))
(if (or (<= phi2 -7.0) (not (<= phi2 8e-141)))
(atan2
t_1
(- t_0 (* (* (cos phi2) (sin phi1)) (cos (- lambda1 lambda2)))))
(atan2
t_1
(-
t_0
(*
(sin phi1)
(+
(* (sin lambda2) (sin lambda1))
(* (cos lambda1) (cos lambda2)))))))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = cos(phi1) * sin(phi2);
double t_1 = fma(-sin(lambda2), cos(lambda1), (cos(lambda2) * sin(lambda1))) * cos(phi2);
double tmp;
if ((phi2 <= -7.0) || !(phi2 <= 8e-141)) {
tmp = atan2(t_1, (t_0 - ((cos(phi2) * sin(phi1)) * cos((lambda1 - lambda2)))));
} else {
tmp = atan2(t_1, (t_0 - (sin(phi1) * ((sin(lambda2) * sin(lambda1)) + (cos(lambda1) * cos(lambda2))))));
}
return tmp;
}
function code(lambda1, lambda2, phi1, phi2) t_0 = Float64(cos(phi1) * sin(phi2)) t_1 = Float64(fma(Float64(-sin(lambda2)), cos(lambda1), Float64(cos(lambda2) * sin(lambda1))) * cos(phi2)) tmp = 0.0 if ((phi2 <= -7.0) || !(phi2 <= 8e-141)) tmp = atan(t_1, Float64(t_0 - Float64(Float64(cos(phi2) * sin(phi1)) * cos(Float64(lambda1 - lambda2))))); else tmp = atan(t_1, Float64(t_0 - Float64(sin(phi1) * Float64(Float64(sin(lambda2) * sin(lambda1)) + Float64(cos(lambda1) * cos(lambda2)))))); end return tmp end
code[lambda1_, lambda2_, phi1_, phi2_] := Block[{t$95$0 = N[(N[Cos[phi1], $MachinePrecision] * N[Sin[phi2], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[((-N[Sin[lambda2], $MachinePrecision]) * N[Cos[lambda1], $MachinePrecision] + N[(N[Cos[lambda2], $MachinePrecision] * N[Sin[lambda1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision]}, If[Or[LessEqual[phi2, -7.0], N[Not[LessEqual[phi2, 8e-141]], $MachinePrecision]], N[ArcTan[t$95$1 / N[(t$95$0 - N[(N[(N[Cos[phi2], $MachinePrecision] * N[Sin[phi1], $MachinePrecision]), $MachinePrecision] * N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], N[ArcTan[t$95$1 / N[(t$95$0 - N[(N[Sin[phi1], $MachinePrecision] * N[(N[(N[Sin[lambda2], $MachinePrecision] * N[Sin[lambda1], $MachinePrecision]), $MachinePrecision] + N[(N[Cos[lambda1], $MachinePrecision] * N[Cos[lambda2], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \cos \phi_1 \cdot \sin \phi_2\\
t_1 := \mathsf{fma}\left(-\sin \lambda_2, \cos \lambda_1, \cos \lambda_2 \cdot \sin \lambda_1\right) \cdot \cos \phi_2\\
\mathbf{if}\;\phi_2 \leq -7 \lor \neg \left(\phi_2 \leq 8 \cdot 10^{-141}\right):\\
\;\;\;\;\tan^{-1}_* \frac{t\_1}{t\_0 - \left(\cos \phi_2 \cdot \sin \phi_1\right) \cdot \cos \left(\lambda_1 - \lambda_2\right)}\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1}_* \frac{t\_1}{t\_0 - \sin \phi_1 \cdot \left(\sin \lambda_2 \cdot \sin \lambda_1 + \cos \lambda_1 \cdot \cos \lambda_2\right)}\\
\end{array}
\end{array}
if phi2 < -7 or 8.0000000000000003e-141 < phi2 Initial program 76.3%
sin-diff87.6%
sub-neg87.6%
Applied egg-rr87.6%
+-commutative87.6%
distribute-rgt-neg-in87.6%
sin-neg87.6%
*-commutative87.6%
fma-define87.6%
sin-neg87.6%
cos-neg87.6%
*-commutative87.6%
cos-neg87.6%
Simplified87.6%
if -7 < phi2 < 8.0000000000000003e-141Initial program 86.3%
sin-diff91.3%
sub-neg91.3%
Applied egg-rr91.3%
+-commutative91.3%
distribute-rgt-neg-in91.3%
sin-neg91.3%
*-commutative91.3%
fma-define91.3%
sin-neg91.3%
cos-neg91.3%
*-commutative91.3%
cos-neg91.3%
Simplified91.3%
cos-diff99.9%
*-commutative99.9%
Applied egg-rr99.9%
fma-define99.9%
Simplified99.9%
Taylor expanded in phi2 around 0 99.8%
Final simplification92.8%
(FPCore (lambda1 lambda2 phi1 phi2)
:precision binary64
(atan2
(*
(cos phi2)
(- (* (cos lambda2) (sin lambda1)) (* (sin lambda2) (cos lambda1))))
(-
(* (cos phi1) (sin phi2))
(*
(cos phi2)
(*
(sin phi1)
(+ (* (sin lambda2) (sin lambda1)) (* (cos lambda1) (cos lambda2))))))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
return atan2((cos(phi2) * ((cos(lambda2) * sin(lambda1)) - (sin(lambda2) * cos(lambda1)))), ((cos(phi1) * sin(phi2)) - (cos(phi2) * (sin(phi1) * ((sin(lambda2) * sin(lambda1)) + (cos(lambda1) * cos(lambda2)))))));
}
real(8) function code(lambda1, lambda2, phi1, phi2)
real(8), intent (in) :: lambda1
real(8), intent (in) :: lambda2
real(8), intent (in) :: phi1
real(8), intent (in) :: phi2
code = atan2((cos(phi2) * ((cos(lambda2) * sin(lambda1)) - (sin(lambda2) * cos(lambda1)))), ((cos(phi1) * sin(phi2)) - (cos(phi2) * (sin(phi1) * ((sin(lambda2) * sin(lambda1)) + (cos(lambda1) * cos(lambda2)))))))
end function
public static double code(double lambda1, double lambda2, double phi1, double phi2) {
return Math.atan2((Math.cos(phi2) * ((Math.cos(lambda2) * Math.sin(lambda1)) - (Math.sin(lambda2) * Math.cos(lambda1)))), ((Math.cos(phi1) * Math.sin(phi2)) - (Math.cos(phi2) * (Math.sin(phi1) * ((Math.sin(lambda2) * Math.sin(lambda1)) + (Math.cos(lambda1) * Math.cos(lambda2)))))));
}
def code(lambda1, lambda2, phi1, phi2): return math.atan2((math.cos(phi2) * ((math.cos(lambda2) * math.sin(lambda1)) - (math.sin(lambda2) * math.cos(lambda1)))), ((math.cos(phi1) * math.sin(phi2)) - (math.cos(phi2) * (math.sin(phi1) * ((math.sin(lambda2) * math.sin(lambda1)) + (math.cos(lambda1) * math.cos(lambda2)))))))
function code(lambda1, lambda2, phi1, phi2) return atan(Float64(cos(phi2) * Float64(Float64(cos(lambda2) * sin(lambda1)) - Float64(sin(lambda2) * cos(lambda1)))), Float64(Float64(cos(phi1) * sin(phi2)) - Float64(cos(phi2) * Float64(sin(phi1) * Float64(Float64(sin(lambda2) * sin(lambda1)) + Float64(cos(lambda1) * cos(lambda2))))))) end
function tmp = code(lambda1, lambda2, phi1, phi2) tmp = atan2((cos(phi2) * ((cos(lambda2) * sin(lambda1)) - (sin(lambda2) * cos(lambda1)))), ((cos(phi1) * sin(phi2)) - (cos(phi2) * (sin(phi1) * ((sin(lambda2) * sin(lambda1)) + (cos(lambda1) * cos(lambda2))))))); end
code[lambda1_, lambda2_, phi1_, phi2_] := N[ArcTan[N[(N[Cos[phi2], $MachinePrecision] * N[(N[(N[Cos[lambda2], $MachinePrecision] * N[Sin[lambda1], $MachinePrecision]), $MachinePrecision] - N[(N[Sin[lambda2], $MachinePrecision] * N[Cos[lambda1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(N[(N[Cos[phi1], $MachinePrecision] * N[Sin[phi2], $MachinePrecision]), $MachinePrecision] - N[(N[Cos[phi2], $MachinePrecision] * N[(N[Sin[phi1], $MachinePrecision] * N[(N[(N[Sin[lambda2], $MachinePrecision] * N[Sin[lambda1], $MachinePrecision]), $MachinePrecision] + N[(N[Cos[lambda1], $MachinePrecision] * N[Cos[lambda2], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
\\
\tan^{-1}_* \frac{\cos \phi_2 \cdot \left(\cos \lambda_2 \cdot \sin \lambda_1 - \sin \lambda_2 \cdot \cos \lambda_1\right)}{\cos \phi_1 \cdot \sin \phi_2 - \cos \phi_2 \cdot \left(\sin \phi_1 \cdot \left(\sin \lambda_2 \cdot \sin \lambda_1 + \cos \lambda_1 \cdot \cos \lambda_2\right)\right)}
\end{array}
Initial program 80.5%
sin-diff89.1%
sub-neg89.1%
Applied egg-rr89.1%
+-commutative89.1%
distribute-rgt-neg-in89.1%
sin-neg89.1%
*-commutative89.1%
fma-define89.2%
sin-neg89.2%
cos-neg89.2%
*-commutative89.2%
cos-neg89.2%
Simplified89.2%
cos-diff99.8%
*-commutative99.8%
Applied egg-rr99.8%
fma-define99.8%
Simplified99.8%
Taylor expanded in lambda2 around inf 99.8%
neg-mul-199.8%
+-commutative99.8%
sub-neg99.8%
cos-neg99.8%
*-commutative99.8%
cos-neg99.8%
Simplified99.8%
Taylor expanded in phi1 around inf 99.7%
Final simplification99.7%
(FPCore (lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (* (cos lambda2) (sin lambda1))) (t_1 (* (cos phi1) (sin phi2))))
(if (or (<= phi2 -7.0) (not (<= phi2 8e-141)))
(atan2
(* (fma (- (sin lambda2)) (cos lambda1) t_0) (cos phi2))
(- t_1 (* (* (cos phi2) (sin phi1)) (cos (- lambda1 lambda2)))))
(atan2
(* (cos phi2) (- t_0 (* (sin lambda2) (cos lambda1))))
(-
t_1
(*
(sin phi1)
(+
(* (sin lambda2) (sin lambda1))
(* (cos lambda1) (cos lambda2)))))))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = cos(lambda2) * sin(lambda1);
double t_1 = cos(phi1) * sin(phi2);
double tmp;
if ((phi2 <= -7.0) || !(phi2 <= 8e-141)) {
tmp = atan2((fma(-sin(lambda2), cos(lambda1), t_0) * cos(phi2)), (t_1 - ((cos(phi2) * sin(phi1)) * cos((lambda1 - lambda2)))));
} else {
tmp = atan2((cos(phi2) * (t_0 - (sin(lambda2) * cos(lambda1)))), (t_1 - (sin(phi1) * ((sin(lambda2) * sin(lambda1)) + (cos(lambda1) * cos(lambda2))))));
}
return tmp;
}
function code(lambda1, lambda2, phi1, phi2) t_0 = Float64(cos(lambda2) * sin(lambda1)) t_1 = Float64(cos(phi1) * sin(phi2)) tmp = 0.0 if ((phi2 <= -7.0) || !(phi2 <= 8e-141)) tmp = atan(Float64(fma(Float64(-sin(lambda2)), cos(lambda1), t_0) * cos(phi2)), Float64(t_1 - Float64(Float64(cos(phi2) * sin(phi1)) * cos(Float64(lambda1 - lambda2))))); else tmp = atan(Float64(cos(phi2) * Float64(t_0 - Float64(sin(lambda2) * cos(lambda1)))), Float64(t_1 - Float64(sin(phi1) * Float64(Float64(sin(lambda2) * sin(lambda1)) + Float64(cos(lambda1) * cos(lambda2)))))); end return tmp end
code[lambda1_, lambda2_, phi1_, phi2_] := Block[{t$95$0 = N[(N[Cos[lambda2], $MachinePrecision] * N[Sin[lambda1], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[Cos[phi1], $MachinePrecision] * N[Sin[phi2], $MachinePrecision]), $MachinePrecision]}, If[Or[LessEqual[phi2, -7.0], N[Not[LessEqual[phi2, 8e-141]], $MachinePrecision]], N[ArcTan[N[(N[((-N[Sin[lambda2], $MachinePrecision]) * N[Cos[lambda1], $MachinePrecision] + t$95$0), $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision] / N[(t$95$1 - N[(N[(N[Cos[phi2], $MachinePrecision] * N[Sin[phi1], $MachinePrecision]), $MachinePrecision] * N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], N[ArcTan[N[(N[Cos[phi2], $MachinePrecision] * N[(t$95$0 - N[(N[Sin[lambda2], $MachinePrecision] * N[Cos[lambda1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(t$95$1 - N[(N[Sin[phi1], $MachinePrecision] * N[(N[(N[Sin[lambda2], $MachinePrecision] * N[Sin[lambda1], $MachinePrecision]), $MachinePrecision] + N[(N[Cos[lambda1], $MachinePrecision] * N[Cos[lambda2], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \cos \lambda_2 \cdot \sin \lambda_1\\
t_1 := \cos \phi_1 \cdot \sin \phi_2\\
\mathbf{if}\;\phi_2 \leq -7 \lor \neg \left(\phi_2 \leq 8 \cdot 10^{-141}\right):\\
\;\;\;\;\tan^{-1}_* \frac{\mathsf{fma}\left(-\sin \lambda_2, \cos \lambda_1, t\_0\right) \cdot \cos \phi_2}{t\_1 - \left(\cos \phi_2 \cdot \sin \phi_1\right) \cdot \cos \left(\lambda_1 - \lambda_2\right)}\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1}_* \frac{\cos \phi_2 \cdot \left(t\_0 - \sin \lambda_2 \cdot \cos \lambda_1\right)}{t\_1 - \sin \phi_1 \cdot \left(\sin \lambda_2 \cdot \sin \lambda_1 + \cos \lambda_1 \cdot \cos \lambda_2\right)}\\
\end{array}
\end{array}
if phi2 < -7 or 8.0000000000000003e-141 < phi2 Initial program 76.3%
sin-diff87.6%
sub-neg87.6%
Applied egg-rr87.6%
+-commutative87.6%
distribute-rgt-neg-in87.6%
sin-neg87.6%
*-commutative87.6%
fma-define87.6%
sin-neg87.6%
cos-neg87.6%
*-commutative87.6%
cos-neg87.6%
Simplified87.6%
if -7 < phi2 < 8.0000000000000003e-141Initial program 86.3%
sin-diff91.3%
sub-neg91.3%
Applied egg-rr91.3%
+-commutative91.3%
distribute-rgt-neg-in91.3%
sin-neg91.3%
*-commutative91.3%
fma-define91.3%
sin-neg91.3%
cos-neg91.3%
*-commutative91.3%
cos-neg91.3%
Simplified91.3%
cos-diff99.9%
*-commutative99.9%
Applied egg-rr99.9%
fma-define99.9%
Simplified99.9%
Taylor expanded in lambda2 around inf 99.9%
neg-mul-199.9%
+-commutative99.9%
sub-neg99.9%
cos-neg99.9%
*-commutative99.9%
cos-neg99.9%
Simplified99.9%
Taylor expanded in phi2 around 0 99.8%
Final simplification92.8%
(FPCore (lambda1 lambda2 phi1 phi2) :precision binary64 (atan2 (* (fma (- (sin lambda2)) (cos lambda1) (* (cos lambda2) (sin lambda1))) (cos phi2)) (- (* (cos phi1) (sin phi2)) (* (* (cos phi2) (sin phi1)) (cos (- lambda1 lambda2))))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
return atan2((fma(-sin(lambda2), cos(lambda1), (cos(lambda2) * sin(lambda1))) * cos(phi2)), ((cos(phi1) * sin(phi2)) - ((cos(phi2) * sin(phi1)) * cos((lambda1 - lambda2)))));
}
function code(lambda1, lambda2, phi1, phi2) return atan(Float64(fma(Float64(-sin(lambda2)), cos(lambda1), Float64(cos(lambda2) * sin(lambda1))) * cos(phi2)), Float64(Float64(cos(phi1) * sin(phi2)) - Float64(Float64(cos(phi2) * sin(phi1)) * cos(Float64(lambda1 - lambda2))))) end
code[lambda1_, lambda2_, phi1_, phi2_] := N[ArcTan[N[(N[((-N[Sin[lambda2], $MachinePrecision]) * N[Cos[lambda1], $MachinePrecision] + N[(N[Cos[lambda2], $MachinePrecision] * N[Sin[lambda1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision] / N[(N[(N[Cos[phi1], $MachinePrecision] * N[Sin[phi2], $MachinePrecision]), $MachinePrecision] - N[(N[(N[Cos[phi2], $MachinePrecision] * N[Sin[phi1], $MachinePrecision]), $MachinePrecision] * N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
\\
\tan^{-1}_* \frac{\mathsf{fma}\left(-\sin \lambda_2, \cos \lambda_1, \cos \lambda_2 \cdot \sin \lambda_1\right) \cdot \cos \phi_2}{\cos \phi_1 \cdot \sin \phi_2 - \left(\cos \phi_2 \cdot \sin \phi_1\right) \cdot \cos \left(\lambda_1 - \lambda_2\right)}
\end{array}
Initial program 80.5%
sin-diff89.1%
sub-neg89.1%
Applied egg-rr89.1%
+-commutative89.1%
distribute-rgt-neg-in89.1%
sin-neg89.1%
*-commutative89.1%
fma-define89.2%
sin-neg89.2%
cos-neg89.2%
*-commutative89.2%
cos-neg89.2%
Simplified89.2%
Final simplification89.2%
(FPCore (lambda1 lambda2 phi1 phi2) :precision binary64 (atan2 (* (fma (- (sin lambda2)) (cos lambda1) (* (cos lambda2) (sin lambda1))) (cos phi2)) (- (* (cos phi1) (sin phi2)) (* (sin phi1) (* (cos phi2) (cos (- lambda1 lambda2)))))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
return atan2((fma(-sin(lambda2), cos(lambda1), (cos(lambda2) * sin(lambda1))) * cos(phi2)), ((cos(phi1) * sin(phi2)) - (sin(phi1) * (cos(phi2) * cos((lambda1 - lambda2))))));
}
function code(lambda1, lambda2, phi1, phi2) return atan(Float64(fma(Float64(-sin(lambda2)), cos(lambda1), Float64(cos(lambda2) * sin(lambda1))) * cos(phi2)), Float64(Float64(cos(phi1) * sin(phi2)) - Float64(sin(phi1) * Float64(cos(phi2) * cos(Float64(lambda1 - lambda2)))))) end
code[lambda1_, lambda2_, phi1_, phi2_] := N[ArcTan[N[(N[((-N[Sin[lambda2], $MachinePrecision]) * N[Cos[lambda1], $MachinePrecision] + N[(N[Cos[lambda2], $MachinePrecision] * N[Sin[lambda1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision] / N[(N[(N[Cos[phi1], $MachinePrecision] * N[Sin[phi2], $MachinePrecision]), $MachinePrecision] - N[(N[Sin[phi1], $MachinePrecision] * N[(N[Cos[phi2], $MachinePrecision] * N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
\\
\tan^{-1}_* \frac{\mathsf{fma}\left(-\sin \lambda_2, \cos \lambda_1, \cos \lambda_2 \cdot \sin \lambda_1\right) \cdot \cos \phi_2}{\cos \phi_1 \cdot \sin \phi_2 - \sin \phi_1 \cdot \left(\cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)\right)}
\end{array}
Initial program 80.5%
sin-diff89.1%
sub-neg89.1%
Applied egg-rr89.1%
+-commutative89.1%
distribute-rgt-neg-in89.1%
sin-neg89.1%
*-commutative89.1%
fma-define89.2%
sin-neg89.2%
cos-neg89.2%
*-commutative89.2%
cos-neg89.2%
Simplified89.2%
Taylor expanded in phi1 around inf 89.1%
associate-*r*89.2%
*-commutative89.2%
Simplified89.2%
(FPCore (lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (* (cos phi2) (sin phi1))))
(if (or (<= lambda1 -4.3e-5) (not (<= lambda1 1.75e-31)))
(atan2
(*
(cos phi2)
(- (* (cos lambda2) (sin lambda1)) (* (sin lambda2) (cos lambda1))))
(- (* (cos phi1) (sin phi2)) (* (cos lambda1) t_0)))
(atan2
(* (cos phi2) (sin (- lambda1 lambda2)))
(-
(* (sin phi2) (cbrt (pow (cos phi1) 3.0)))
(* t_0 (cos (- lambda1 lambda2))))))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = cos(phi2) * sin(phi1);
double tmp;
if ((lambda1 <= -4.3e-5) || !(lambda1 <= 1.75e-31)) {
tmp = atan2((cos(phi2) * ((cos(lambda2) * sin(lambda1)) - (sin(lambda2) * cos(lambda1)))), ((cos(phi1) * sin(phi2)) - (cos(lambda1) * t_0)));
} else {
tmp = atan2((cos(phi2) * sin((lambda1 - lambda2))), ((sin(phi2) * cbrt(pow(cos(phi1), 3.0))) - (t_0 * cos((lambda1 - lambda2)))));
}
return tmp;
}
public static double code(double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = Math.cos(phi2) * Math.sin(phi1);
double tmp;
if ((lambda1 <= -4.3e-5) || !(lambda1 <= 1.75e-31)) {
tmp = Math.atan2((Math.cos(phi2) * ((Math.cos(lambda2) * Math.sin(lambda1)) - (Math.sin(lambda2) * Math.cos(lambda1)))), ((Math.cos(phi1) * Math.sin(phi2)) - (Math.cos(lambda1) * t_0)));
} else {
tmp = Math.atan2((Math.cos(phi2) * Math.sin((lambda1 - lambda2))), ((Math.sin(phi2) * Math.cbrt(Math.pow(Math.cos(phi1), 3.0))) - (t_0 * Math.cos((lambda1 - lambda2)))));
}
return tmp;
}
function code(lambda1, lambda2, phi1, phi2) t_0 = Float64(cos(phi2) * sin(phi1)) tmp = 0.0 if ((lambda1 <= -4.3e-5) || !(lambda1 <= 1.75e-31)) tmp = atan(Float64(cos(phi2) * Float64(Float64(cos(lambda2) * sin(lambda1)) - Float64(sin(lambda2) * cos(lambda1)))), Float64(Float64(cos(phi1) * sin(phi2)) - Float64(cos(lambda1) * t_0))); else tmp = atan(Float64(cos(phi2) * sin(Float64(lambda1 - lambda2))), Float64(Float64(sin(phi2) * cbrt((cos(phi1) ^ 3.0))) - Float64(t_0 * cos(Float64(lambda1 - lambda2))))); end return tmp end
code[lambda1_, lambda2_, phi1_, phi2_] := Block[{t$95$0 = N[(N[Cos[phi2], $MachinePrecision] * N[Sin[phi1], $MachinePrecision]), $MachinePrecision]}, If[Or[LessEqual[lambda1, -4.3e-5], N[Not[LessEqual[lambda1, 1.75e-31]], $MachinePrecision]], N[ArcTan[N[(N[Cos[phi2], $MachinePrecision] * N[(N[(N[Cos[lambda2], $MachinePrecision] * N[Sin[lambda1], $MachinePrecision]), $MachinePrecision] - N[(N[Sin[lambda2], $MachinePrecision] * N[Cos[lambda1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(N[(N[Cos[phi1], $MachinePrecision] * N[Sin[phi2], $MachinePrecision]), $MachinePrecision] - N[(N[Cos[lambda1], $MachinePrecision] * t$95$0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], N[ArcTan[N[(N[Cos[phi2], $MachinePrecision] * N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / N[(N[(N[Sin[phi2], $MachinePrecision] * N[Power[N[Power[N[Cos[phi1], $MachinePrecision], 3.0], $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision] - N[(t$95$0 * N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \cos \phi_2 \cdot \sin \phi_1\\
\mathbf{if}\;\lambda_1 \leq -4.3 \cdot 10^{-5} \lor \neg \left(\lambda_1 \leq 1.75 \cdot 10^{-31}\right):\\
\;\;\;\;\tan^{-1}_* \frac{\cos \phi_2 \cdot \left(\cos \lambda_2 \cdot \sin \lambda_1 - \sin \lambda_2 \cdot \cos \lambda_1\right)}{\cos \phi_1 \cdot \sin \phi_2 - \cos \lambda_1 \cdot t\_0}\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\sin \phi_2 \cdot \sqrt[3]{{\cos \phi_1}^{3}} - t\_0 \cdot \cos \left(\lambda_1 - \lambda_2\right)}\\
\end{array}
\end{array}
if lambda1 < -4.3000000000000002e-5 or 1.74999999999999993e-31 < lambda1 Initial program 57.0%
sin-diff76.1%
sub-neg76.1%
Applied egg-rr76.1%
+-commutative76.1%
distribute-rgt-neg-in76.1%
sin-neg76.1%
*-commutative76.1%
fma-define76.1%
sin-neg76.1%
cos-neg76.1%
*-commutative76.1%
cos-neg76.1%
Simplified76.1%
cos-diff99.8%
*-commutative99.8%
Applied egg-rr99.8%
fma-define99.8%
Simplified99.8%
Taylor expanded in lambda2 around inf 99.8%
neg-mul-199.8%
+-commutative99.8%
sub-neg99.8%
cos-neg99.8%
*-commutative99.8%
cos-neg99.8%
Simplified99.8%
Taylor expanded in lambda2 around 0 76.3%
if -4.3000000000000002e-5 < lambda1 < 1.74999999999999993e-31Initial program 99.8%
Taylor expanded in lambda1 around inf 71.1%
mul-1-neg71.1%
unsub-neg71.1%
Simplified71.1%
add-cbrt-cube71.2%
pow371.2%
Applied egg-rr71.2%
Taylor expanded in lambda1 around 0 99.8%
neg-mul-199.8%
sub-neg99.8%
Simplified99.8%
Final simplification89.3%
(FPCore (lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (* (* (cos phi2) (sin phi1)) (cos (- lambda1 lambda2))))
(t_1 (sin (- lambda1 lambda2))))
(if (<= phi1 -0.055)
(atan2
(* (cos phi2) (expm1 (log1p t_1)))
(- (* (cos phi1) (sin phi2)) t_0))
(if (<= phi1 7e-6)
(atan2
(*
(fma (- (sin lambda2)) (cos lambda1) (* (cos lambda2) (sin lambda1)))
(cos phi2))
(- (sin phi2) t_0))
(atan2
(* (cos phi2) t_1)
(- (* (sin phi2) (cbrt (pow (cos phi1) 3.0))) t_0))))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = (cos(phi2) * sin(phi1)) * cos((lambda1 - lambda2));
double t_1 = sin((lambda1 - lambda2));
double tmp;
if (phi1 <= -0.055) {
tmp = atan2((cos(phi2) * expm1(log1p(t_1))), ((cos(phi1) * sin(phi2)) - t_0));
} else if (phi1 <= 7e-6) {
tmp = atan2((fma(-sin(lambda2), cos(lambda1), (cos(lambda2) * sin(lambda1))) * cos(phi2)), (sin(phi2) - t_0));
} else {
tmp = atan2((cos(phi2) * t_1), ((sin(phi2) * cbrt(pow(cos(phi1), 3.0))) - t_0));
}
return tmp;
}
function code(lambda1, lambda2, phi1, phi2) t_0 = Float64(Float64(cos(phi2) * sin(phi1)) * cos(Float64(lambda1 - lambda2))) t_1 = sin(Float64(lambda1 - lambda2)) tmp = 0.0 if (phi1 <= -0.055) tmp = atan(Float64(cos(phi2) * expm1(log1p(t_1))), Float64(Float64(cos(phi1) * sin(phi2)) - t_0)); elseif (phi1 <= 7e-6) tmp = atan(Float64(fma(Float64(-sin(lambda2)), cos(lambda1), Float64(cos(lambda2) * sin(lambda1))) * cos(phi2)), Float64(sin(phi2) - t_0)); else tmp = atan(Float64(cos(phi2) * t_1), Float64(Float64(sin(phi2) * cbrt((cos(phi1) ^ 3.0))) - t_0)); end return tmp end
code[lambda1_, lambda2_, phi1_, phi2_] := Block[{t$95$0 = N[(N[(N[Cos[phi2], $MachinePrecision] * N[Sin[phi1], $MachinePrecision]), $MachinePrecision] * N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]}, If[LessEqual[phi1, -0.055], N[ArcTan[N[(N[Cos[phi2], $MachinePrecision] * N[(Exp[N[Log[1 + t$95$1], $MachinePrecision]] - 1), $MachinePrecision]), $MachinePrecision] / N[(N[(N[Cos[phi1], $MachinePrecision] * N[Sin[phi2], $MachinePrecision]), $MachinePrecision] - t$95$0), $MachinePrecision]], $MachinePrecision], If[LessEqual[phi1, 7e-6], N[ArcTan[N[(N[((-N[Sin[lambda2], $MachinePrecision]) * N[Cos[lambda1], $MachinePrecision] + N[(N[Cos[lambda2], $MachinePrecision] * N[Sin[lambda1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision] / N[(N[Sin[phi2], $MachinePrecision] - t$95$0), $MachinePrecision]], $MachinePrecision], N[ArcTan[N[(N[Cos[phi2], $MachinePrecision] * t$95$1), $MachinePrecision] / N[(N[(N[Sin[phi2], $MachinePrecision] * N[Power[N[Power[N[Cos[phi1], $MachinePrecision], 3.0], $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision] - t$95$0), $MachinePrecision]], $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(\cos \phi_2 \cdot \sin \phi_1\right) \cdot \cos \left(\lambda_1 - \lambda_2\right)\\
t_1 := \sin \left(\lambda_1 - \lambda_2\right)\\
\mathbf{if}\;\phi_1 \leq -0.055:\\
\;\;\;\;\tan^{-1}_* \frac{\cos \phi_2 \cdot \mathsf{expm1}\left(\mathsf{log1p}\left(t\_1\right)\right)}{\cos \phi_1 \cdot \sin \phi_2 - t\_0}\\
\mathbf{elif}\;\phi_1 \leq 7 \cdot 10^{-6}:\\
\;\;\;\;\tan^{-1}_* \frac{\mathsf{fma}\left(-\sin \lambda_2, \cos \lambda_1, \cos \lambda_2 \cdot \sin \lambda_1\right) \cdot \cos \phi_2}{\sin \phi_2 - t\_0}\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1}_* \frac{\cos \phi_2 \cdot t\_1}{\sin \phi_2 \cdot \sqrt[3]{{\cos \phi_1}^{3}} - t\_0}\\
\end{array}
\end{array}
if phi1 < -0.0550000000000000003Initial program 75.2%
expm1-log1p-u75.3%
expm1-undefine63.3%
Applied egg-rr63.3%
expm1-define75.3%
Simplified75.3%
if -0.0550000000000000003 < phi1 < 6.99999999999999989e-6Initial program 84.4%
sin-diff99.9%
sub-neg99.9%
Applied egg-rr99.9%
+-commutative99.9%
distribute-rgt-neg-in99.9%
sin-neg99.9%
*-commutative99.9%
fma-define99.9%
sin-neg99.9%
cos-neg99.9%
*-commutative99.9%
cos-neg99.9%
Simplified99.9%
Taylor expanded in phi1 around 0 99.9%
if 6.99999999999999989e-6 < phi1 Initial program 80.1%
Taylor expanded in lambda1 around inf 60.6%
mul-1-neg60.6%
unsub-neg60.6%
Simplified60.6%
add-cbrt-cube60.6%
pow360.6%
Applied egg-rr60.6%
Taylor expanded in lambda1 around 0 80.1%
neg-mul-180.1%
sub-neg80.1%
Simplified80.1%
Final simplification87.9%
(FPCore (lambda1 lambda2 phi1 phi2) :precision binary64 (atan2 (* (cos phi2) (- (* (cos lambda2) (sin lambda1)) (* (sin lambda2) (cos lambda1)))) (- (* (cos phi1) (sin phi2)) (* (cos phi2) (* (sin phi1) (cos (- lambda1 lambda2)))))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
return atan2((cos(phi2) * ((cos(lambda2) * sin(lambda1)) - (sin(lambda2) * cos(lambda1)))), ((cos(phi1) * sin(phi2)) - (cos(phi2) * (sin(phi1) * cos((lambda1 - lambda2))))));
}
real(8) function code(lambda1, lambda2, phi1, phi2)
real(8), intent (in) :: lambda1
real(8), intent (in) :: lambda2
real(8), intent (in) :: phi1
real(8), intent (in) :: phi2
code = atan2((cos(phi2) * ((cos(lambda2) * sin(lambda1)) - (sin(lambda2) * cos(lambda1)))), ((cos(phi1) * sin(phi2)) - (cos(phi2) * (sin(phi1) * cos((lambda1 - lambda2))))))
end function
public static double code(double lambda1, double lambda2, double phi1, double phi2) {
return Math.atan2((Math.cos(phi2) * ((Math.cos(lambda2) * Math.sin(lambda1)) - (Math.sin(lambda2) * Math.cos(lambda1)))), ((Math.cos(phi1) * Math.sin(phi2)) - (Math.cos(phi2) * (Math.sin(phi1) * Math.cos((lambda1 - lambda2))))));
}
def code(lambda1, lambda2, phi1, phi2): return math.atan2((math.cos(phi2) * ((math.cos(lambda2) * math.sin(lambda1)) - (math.sin(lambda2) * math.cos(lambda1)))), ((math.cos(phi1) * math.sin(phi2)) - (math.cos(phi2) * (math.sin(phi1) * math.cos((lambda1 - lambda2))))))
function code(lambda1, lambda2, phi1, phi2) return atan(Float64(cos(phi2) * Float64(Float64(cos(lambda2) * sin(lambda1)) - Float64(sin(lambda2) * cos(lambda1)))), Float64(Float64(cos(phi1) * sin(phi2)) - Float64(cos(phi2) * Float64(sin(phi1) * cos(Float64(lambda1 - lambda2)))))) end
function tmp = code(lambda1, lambda2, phi1, phi2) tmp = atan2((cos(phi2) * ((cos(lambda2) * sin(lambda1)) - (sin(lambda2) * cos(lambda1)))), ((cos(phi1) * sin(phi2)) - (cos(phi2) * (sin(phi1) * cos((lambda1 - lambda2)))))); end
code[lambda1_, lambda2_, phi1_, phi2_] := N[ArcTan[N[(N[Cos[phi2], $MachinePrecision] * N[(N[(N[Cos[lambda2], $MachinePrecision] * N[Sin[lambda1], $MachinePrecision]), $MachinePrecision] - N[(N[Sin[lambda2], $MachinePrecision] * N[Cos[lambda1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(N[(N[Cos[phi1], $MachinePrecision] * N[Sin[phi2], $MachinePrecision]), $MachinePrecision] - N[(N[Cos[phi2], $MachinePrecision] * N[(N[Sin[phi1], $MachinePrecision] * N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
\\
\tan^{-1}_* \frac{\cos \phi_2 \cdot \left(\cos \lambda_2 \cdot \sin \lambda_1 - \sin \lambda_2 \cdot \cos \lambda_1\right)}{\cos \phi_1 \cdot \sin \phi_2 - \cos \phi_2 \cdot \left(\sin \phi_1 \cdot \cos \left(\lambda_1 - \lambda_2\right)\right)}
\end{array}
Initial program 80.5%
*-commutative80.5%
associate-*l*80.5%
Simplified80.5%
sin-diff89.1%
Applied egg-rr89.1%
Final simplification89.1%
(FPCore (lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (* (* (cos phi2) (sin phi1)) (cos (- lambda1 lambda2))))
(t_1 (sin (- lambda1 lambda2))))
(if (<= phi1 -3.1e-128)
(atan2
(* (cos phi2) (expm1 (log1p t_1)))
(- (* (cos phi1) (sin phi2)) t_0))
(if (<= phi1 3.2e-15)
(atan2
(*
(fma (- (sin lambda2)) (cos lambda1) (* (cos lambda2) (sin lambda1)))
(cos phi2))
(sin phi2))
(atan2
(* (cos phi2) t_1)
(- (* (sin phi2) (cbrt (pow (cos phi1) 3.0))) t_0))))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = (cos(phi2) * sin(phi1)) * cos((lambda1 - lambda2));
double t_1 = sin((lambda1 - lambda2));
double tmp;
if (phi1 <= -3.1e-128) {
tmp = atan2((cos(phi2) * expm1(log1p(t_1))), ((cos(phi1) * sin(phi2)) - t_0));
} else if (phi1 <= 3.2e-15) {
tmp = atan2((fma(-sin(lambda2), cos(lambda1), (cos(lambda2) * sin(lambda1))) * cos(phi2)), sin(phi2));
} else {
tmp = atan2((cos(phi2) * t_1), ((sin(phi2) * cbrt(pow(cos(phi1), 3.0))) - t_0));
}
return tmp;
}
function code(lambda1, lambda2, phi1, phi2) t_0 = Float64(Float64(cos(phi2) * sin(phi1)) * cos(Float64(lambda1 - lambda2))) t_1 = sin(Float64(lambda1 - lambda2)) tmp = 0.0 if (phi1 <= -3.1e-128) tmp = atan(Float64(cos(phi2) * expm1(log1p(t_1))), Float64(Float64(cos(phi1) * sin(phi2)) - t_0)); elseif (phi1 <= 3.2e-15) tmp = atan(Float64(fma(Float64(-sin(lambda2)), cos(lambda1), Float64(cos(lambda2) * sin(lambda1))) * cos(phi2)), sin(phi2)); else tmp = atan(Float64(cos(phi2) * t_1), Float64(Float64(sin(phi2) * cbrt((cos(phi1) ^ 3.0))) - t_0)); end return tmp end
code[lambda1_, lambda2_, phi1_, phi2_] := Block[{t$95$0 = N[(N[(N[Cos[phi2], $MachinePrecision] * N[Sin[phi1], $MachinePrecision]), $MachinePrecision] * N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]}, If[LessEqual[phi1, -3.1e-128], N[ArcTan[N[(N[Cos[phi2], $MachinePrecision] * N[(Exp[N[Log[1 + t$95$1], $MachinePrecision]] - 1), $MachinePrecision]), $MachinePrecision] / N[(N[(N[Cos[phi1], $MachinePrecision] * N[Sin[phi2], $MachinePrecision]), $MachinePrecision] - t$95$0), $MachinePrecision]], $MachinePrecision], If[LessEqual[phi1, 3.2e-15], N[ArcTan[N[(N[((-N[Sin[lambda2], $MachinePrecision]) * N[Cos[lambda1], $MachinePrecision] + N[(N[Cos[lambda2], $MachinePrecision] * N[Sin[lambda1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision] / N[Sin[phi2], $MachinePrecision]], $MachinePrecision], N[ArcTan[N[(N[Cos[phi2], $MachinePrecision] * t$95$1), $MachinePrecision] / N[(N[(N[Sin[phi2], $MachinePrecision] * N[Power[N[Power[N[Cos[phi1], $MachinePrecision], 3.0], $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision] - t$95$0), $MachinePrecision]], $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(\cos \phi_2 \cdot \sin \phi_1\right) \cdot \cos \left(\lambda_1 - \lambda_2\right)\\
t_1 := \sin \left(\lambda_1 - \lambda_2\right)\\
\mathbf{if}\;\phi_1 \leq -3.1 \cdot 10^{-128}:\\
\;\;\;\;\tan^{-1}_* \frac{\cos \phi_2 \cdot \mathsf{expm1}\left(\mathsf{log1p}\left(t\_1\right)\right)}{\cos \phi_1 \cdot \sin \phi_2 - t\_0}\\
\mathbf{elif}\;\phi_1 \leq 3.2 \cdot 10^{-15}:\\
\;\;\;\;\tan^{-1}_* \frac{\mathsf{fma}\left(-\sin \lambda_2, \cos \lambda_1, \cos \lambda_2 \cdot \sin \lambda_1\right) \cdot \cos \phi_2}{\sin \phi_2}\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1}_* \frac{\cos \phi_2 \cdot t\_1}{\sin \phi_2 \cdot \sqrt[3]{{\cos \phi_1}^{3}} - t\_0}\\
\end{array}
\end{array}
if phi1 < -3.10000000000000003e-128Initial program 79.4%
expm1-log1p-u79.5%
expm1-undefine64.7%
Applied egg-rr64.7%
expm1-define79.5%
Simplified79.5%
if -3.10000000000000003e-128 < phi1 < 3.1999999999999999e-15Initial program 82.0%
add-sqr-sqrt59.0%
sqrt-unprod82.0%
pow282.0%
Applied egg-rr82.0%
Taylor expanded in phi1 around 0 82.0%
sin-diff99.9%
sub-neg99.9%
Applied egg-rr99.9%
+-commutative99.9%
distribute-rgt-neg-in99.9%
sin-neg99.9%
*-commutative99.9%
fma-define99.9%
sin-neg99.9%
cos-neg99.9%
*-commutative99.9%
cos-neg99.9%
Simplified99.9%
if 3.1999999999999999e-15 < phi1 Initial program 80.1%
Taylor expanded in lambda1 around inf 60.6%
mul-1-neg60.6%
unsub-neg60.6%
Simplified60.6%
add-cbrt-cube60.6%
pow360.6%
Applied egg-rr60.6%
Taylor expanded in lambda1 around 0 80.1%
neg-mul-180.1%
sub-neg80.1%
Simplified80.1%
Final simplification87.3%
(FPCore (lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (cos (- lambda1 lambda2))) (t_1 (sin (- lambda1 lambda2))))
(if (<= phi1 -3.1e-128)
(atan2
(* (cos phi2) (expm1 (log1p t_1)))
(- (* (cos phi1) (sin phi2)) (* (* (cos phi2) (sin phi1)) t_0)))
(if (<= phi1 7e-12)
(atan2
(*
(fma (- (sin lambda2)) (cos lambda1) (* (cos lambda2) (sin lambda1)))
(cos phi2))
(sin phi2))
(atan2
(* (cos phi2) t_1)
(fma (sin phi2) (cos phi1) (* (sin phi1) (* t_0 (- (cos phi2))))))))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = cos((lambda1 - lambda2));
double t_1 = sin((lambda1 - lambda2));
double tmp;
if (phi1 <= -3.1e-128) {
tmp = atan2((cos(phi2) * expm1(log1p(t_1))), ((cos(phi1) * sin(phi2)) - ((cos(phi2) * sin(phi1)) * t_0)));
} else if (phi1 <= 7e-12) {
tmp = atan2((fma(-sin(lambda2), cos(lambda1), (cos(lambda2) * sin(lambda1))) * cos(phi2)), sin(phi2));
} else {
tmp = atan2((cos(phi2) * t_1), fma(sin(phi2), cos(phi1), (sin(phi1) * (t_0 * -cos(phi2)))));
}
return tmp;
}
function code(lambda1, lambda2, phi1, phi2) t_0 = cos(Float64(lambda1 - lambda2)) t_1 = sin(Float64(lambda1 - lambda2)) tmp = 0.0 if (phi1 <= -3.1e-128) tmp = atan(Float64(cos(phi2) * expm1(log1p(t_1))), Float64(Float64(cos(phi1) * sin(phi2)) - Float64(Float64(cos(phi2) * sin(phi1)) * t_0))); elseif (phi1 <= 7e-12) tmp = atan(Float64(fma(Float64(-sin(lambda2)), cos(lambda1), Float64(cos(lambda2) * sin(lambda1))) * cos(phi2)), sin(phi2)); else tmp = atan(Float64(cos(phi2) * t_1), fma(sin(phi2), cos(phi1), Float64(sin(phi1) * Float64(t_0 * Float64(-cos(phi2)))))); end return tmp end
code[lambda1_, lambda2_, phi1_, phi2_] := Block[{t$95$0 = N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$1 = N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]}, If[LessEqual[phi1, -3.1e-128], N[ArcTan[N[(N[Cos[phi2], $MachinePrecision] * N[(Exp[N[Log[1 + t$95$1], $MachinePrecision]] - 1), $MachinePrecision]), $MachinePrecision] / N[(N[(N[Cos[phi1], $MachinePrecision] * N[Sin[phi2], $MachinePrecision]), $MachinePrecision] - N[(N[(N[Cos[phi2], $MachinePrecision] * N[Sin[phi1], $MachinePrecision]), $MachinePrecision] * t$95$0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], If[LessEqual[phi1, 7e-12], N[ArcTan[N[(N[((-N[Sin[lambda2], $MachinePrecision]) * N[Cos[lambda1], $MachinePrecision] + N[(N[Cos[lambda2], $MachinePrecision] * N[Sin[lambda1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision] / N[Sin[phi2], $MachinePrecision]], $MachinePrecision], N[ArcTan[N[(N[Cos[phi2], $MachinePrecision] * t$95$1), $MachinePrecision] / N[(N[Sin[phi2], $MachinePrecision] * N[Cos[phi1], $MachinePrecision] + N[(N[Sin[phi1], $MachinePrecision] * N[(t$95$0 * (-N[Cos[phi2], $MachinePrecision])), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \cos \left(\lambda_1 - \lambda_2\right)\\
t_1 := \sin \left(\lambda_1 - \lambda_2\right)\\
\mathbf{if}\;\phi_1 \leq -3.1 \cdot 10^{-128}:\\
\;\;\;\;\tan^{-1}_* \frac{\cos \phi_2 \cdot \mathsf{expm1}\left(\mathsf{log1p}\left(t\_1\right)\right)}{\cos \phi_1 \cdot \sin \phi_2 - \left(\cos \phi_2 \cdot \sin \phi_1\right) \cdot t\_0}\\
\mathbf{elif}\;\phi_1 \leq 7 \cdot 10^{-12}:\\
\;\;\;\;\tan^{-1}_* \frac{\mathsf{fma}\left(-\sin \lambda_2, \cos \lambda_1, \cos \lambda_2 \cdot \sin \lambda_1\right) \cdot \cos \phi_2}{\sin \phi_2}\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1}_* \frac{\cos \phi_2 \cdot t\_1}{\mathsf{fma}\left(\sin \phi_2, \cos \phi_1, \sin \phi_1 \cdot \left(t\_0 \cdot \left(-\cos \phi_2\right)\right)\right)}\\
\end{array}
\end{array}
if phi1 < -3.10000000000000003e-128Initial program 79.4%
expm1-log1p-u79.5%
expm1-undefine64.7%
Applied egg-rr64.7%
expm1-define79.5%
Simplified79.5%
if -3.10000000000000003e-128 < phi1 < 7.0000000000000001e-12Initial program 82.0%
add-sqr-sqrt59.0%
sqrt-unprod82.0%
pow282.0%
Applied egg-rr82.0%
Taylor expanded in phi1 around 0 82.0%
sin-diff99.9%
sub-neg99.9%
Applied egg-rr99.9%
+-commutative99.9%
distribute-rgt-neg-in99.9%
sin-neg99.9%
*-commutative99.9%
fma-define99.9%
sin-neg99.9%
cos-neg99.9%
*-commutative99.9%
cos-neg99.9%
Simplified99.9%
if 7.0000000000000001e-12 < phi1 Initial program 80.1%
add-sqr-sqrt61.1%
sqrt-unprod65.9%
pow265.9%
Applied egg-rr65.9%
*-commutative65.9%
sqrt-pow180.1%
metadata-eval80.1%
pow180.1%
fma-neg80.1%
*-commutative80.1%
*-commutative80.1%
associate-*l*80.1%
Applied egg-rr80.1%
Final simplification87.3%
(FPCore (lambda1 lambda2 phi1 phi2)
:precision binary64
(if (or (<= phi1 -3.1e-128) (not (<= phi1 3e-12)))
(atan2
(* (cos phi2) (sin (- lambda1 lambda2)))
(fma
(sin phi2)
(cos phi1)
(* (sin phi1) (* (cos (- lambda1 lambda2)) (- (cos phi2))))))
(atan2
(*
(fma (- (sin lambda2)) (cos lambda1) (* (cos lambda2) (sin lambda1)))
(cos phi2))
(sin phi2))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
double tmp;
if ((phi1 <= -3.1e-128) || !(phi1 <= 3e-12)) {
tmp = atan2((cos(phi2) * sin((lambda1 - lambda2))), fma(sin(phi2), cos(phi1), (sin(phi1) * (cos((lambda1 - lambda2)) * -cos(phi2)))));
} else {
tmp = atan2((fma(-sin(lambda2), cos(lambda1), (cos(lambda2) * sin(lambda1))) * cos(phi2)), sin(phi2));
}
return tmp;
}
function code(lambda1, lambda2, phi1, phi2) tmp = 0.0 if ((phi1 <= -3.1e-128) || !(phi1 <= 3e-12)) tmp = atan(Float64(cos(phi2) * sin(Float64(lambda1 - lambda2))), fma(sin(phi2), cos(phi1), Float64(sin(phi1) * Float64(cos(Float64(lambda1 - lambda2)) * Float64(-cos(phi2)))))); else tmp = atan(Float64(fma(Float64(-sin(lambda2)), cos(lambda1), Float64(cos(lambda2) * sin(lambda1))) * cos(phi2)), sin(phi2)); end return tmp end
code[lambda1_, lambda2_, phi1_, phi2_] := If[Or[LessEqual[phi1, -3.1e-128], N[Not[LessEqual[phi1, 3e-12]], $MachinePrecision]], N[ArcTan[N[(N[Cos[phi2], $MachinePrecision] * N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / N[(N[Sin[phi2], $MachinePrecision] * N[Cos[phi1], $MachinePrecision] + N[(N[Sin[phi1], $MachinePrecision] * N[(N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision] * (-N[Cos[phi2], $MachinePrecision])), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], N[ArcTan[N[(N[((-N[Sin[lambda2], $MachinePrecision]) * N[Cos[lambda1], $MachinePrecision] + N[(N[Cos[lambda2], $MachinePrecision] * N[Sin[lambda1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision] / N[Sin[phi2], $MachinePrecision]], $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\phi_1 \leq -3.1 \cdot 10^{-128} \lor \neg \left(\phi_1 \leq 3 \cdot 10^{-12}\right):\\
\;\;\;\;\tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\mathsf{fma}\left(\sin \phi_2, \cos \phi_1, \sin \phi_1 \cdot \left(\cos \left(\lambda_1 - \lambda_2\right) \cdot \left(-\cos \phi_2\right)\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1}_* \frac{\mathsf{fma}\left(-\sin \lambda_2, \cos \lambda_1, \cos \lambda_2 \cdot \sin \lambda_1\right) \cdot \cos \phi_2}{\sin \phi_2}\\
\end{array}
\end{array}
if phi1 < -3.10000000000000003e-128 or 3.0000000000000001e-12 < phi1 Initial program 79.7%
add-sqr-sqrt54.8%
sqrt-unprod64.9%
pow264.9%
Applied egg-rr64.9%
*-commutative64.9%
sqrt-pow179.7%
metadata-eval79.7%
pow179.7%
fma-neg79.7%
*-commutative79.7%
*-commutative79.7%
associate-*l*79.7%
Applied egg-rr79.7%
if -3.10000000000000003e-128 < phi1 < 3.0000000000000001e-12Initial program 82.0%
add-sqr-sqrt59.0%
sqrt-unprod82.0%
pow282.0%
Applied egg-rr82.0%
Taylor expanded in phi1 around 0 82.0%
sin-diff99.9%
sub-neg99.9%
Applied egg-rr99.9%
+-commutative99.9%
distribute-rgt-neg-in99.9%
sin-neg99.9%
*-commutative99.9%
fma-define99.9%
sin-neg99.9%
cos-neg99.9%
*-commutative99.9%
cos-neg99.9%
Simplified99.9%
Final simplification87.3%
(FPCore (lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (* (cos phi2) (sin phi1)))
(t_1 (* (cos phi1) (sin phi2)))
(t_2 (* (cos lambda2) (sin lambda1))))
(if (<= lambda1 -120000000.0)
(atan2
(* (fma (- (sin lambda2)) (cos lambda1) t_2) (cos phi2))
(sin phi2))
(if (<= lambda1 5.6e-8)
(atan2
(* (cos phi2) (sin (- lambda1 lambda2)))
(- t_1 (* (cos lambda2) t_0)))
(if (<= lambda1 1.22e+119)
(atan2
(* (cos phi2) (- t_2 (* (sin lambda2) (cos lambda1))))
(sin phi2))
(atan2
(* (sin lambda1) (cos phi2))
(- t_1 (* t_0 (cos (- lambda1 lambda2))))))))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = cos(phi2) * sin(phi1);
double t_1 = cos(phi1) * sin(phi2);
double t_2 = cos(lambda2) * sin(lambda1);
double tmp;
if (lambda1 <= -120000000.0) {
tmp = atan2((fma(-sin(lambda2), cos(lambda1), t_2) * cos(phi2)), sin(phi2));
} else if (lambda1 <= 5.6e-8) {
tmp = atan2((cos(phi2) * sin((lambda1 - lambda2))), (t_1 - (cos(lambda2) * t_0)));
} else if (lambda1 <= 1.22e+119) {
tmp = atan2((cos(phi2) * (t_2 - (sin(lambda2) * cos(lambda1)))), sin(phi2));
} else {
tmp = atan2((sin(lambda1) * cos(phi2)), (t_1 - (t_0 * cos((lambda1 - lambda2)))));
}
return tmp;
}
function code(lambda1, lambda2, phi1, phi2) t_0 = Float64(cos(phi2) * sin(phi1)) t_1 = Float64(cos(phi1) * sin(phi2)) t_2 = Float64(cos(lambda2) * sin(lambda1)) tmp = 0.0 if (lambda1 <= -120000000.0) tmp = atan(Float64(fma(Float64(-sin(lambda2)), cos(lambda1), t_2) * cos(phi2)), sin(phi2)); elseif (lambda1 <= 5.6e-8) tmp = atan(Float64(cos(phi2) * sin(Float64(lambda1 - lambda2))), Float64(t_1 - Float64(cos(lambda2) * t_0))); elseif (lambda1 <= 1.22e+119) tmp = atan(Float64(cos(phi2) * Float64(t_2 - Float64(sin(lambda2) * cos(lambda1)))), sin(phi2)); else tmp = atan(Float64(sin(lambda1) * cos(phi2)), Float64(t_1 - Float64(t_0 * cos(Float64(lambda1 - lambda2))))); end return tmp end
code[lambda1_, lambda2_, phi1_, phi2_] := Block[{t$95$0 = N[(N[Cos[phi2], $MachinePrecision] * N[Sin[phi1], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[Cos[phi1], $MachinePrecision] * N[Sin[phi2], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(N[Cos[lambda2], $MachinePrecision] * N[Sin[lambda1], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[lambda1, -120000000.0], N[ArcTan[N[(N[((-N[Sin[lambda2], $MachinePrecision]) * N[Cos[lambda1], $MachinePrecision] + t$95$2), $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision] / N[Sin[phi2], $MachinePrecision]], $MachinePrecision], If[LessEqual[lambda1, 5.6e-8], N[ArcTan[N[(N[Cos[phi2], $MachinePrecision] * N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / N[(t$95$1 - N[(N[Cos[lambda2], $MachinePrecision] * t$95$0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], If[LessEqual[lambda1, 1.22e+119], N[ArcTan[N[(N[Cos[phi2], $MachinePrecision] * N[(t$95$2 - N[(N[Sin[lambda2], $MachinePrecision] * N[Cos[lambda1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[Sin[phi2], $MachinePrecision]], $MachinePrecision], N[ArcTan[N[(N[Sin[lambda1], $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision] / N[(t$95$1 - N[(t$95$0 * N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \cos \phi_2 \cdot \sin \phi_1\\
t_1 := \cos \phi_1 \cdot \sin \phi_2\\
t_2 := \cos \lambda_2 \cdot \sin \lambda_1\\
\mathbf{if}\;\lambda_1 \leq -120000000:\\
\;\;\;\;\tan^{-1}_* \frac{\mathsf{fma}\left(-\sin \lambda_2, \cos \lambda_1, t\_2\right) \cdot \cos \phi_2}{\sin \phi_2}\\
\mathbf{elif}\;\lambda_1 \leq 5.6 \cdot 10^{-8}:\\
\;\;\;\;\tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{t\_1 - \cos \lambda_2 \cdot t\_0}\\
\mathbf{elif}\;\lambda_1 \leq 1.22 \cdot 10^{+119}:\\
\;\;\;\;\tan^{-1}_* \frac{\cos \phi_2 \cdot \left(t\_2 - \sin \lambda_2 \cdot \cos \lambda_1\right)}{\sin \phi_2}\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1}_* \frac{\sin \lambda_1 \cdot \cos \phi_2}{t\_1 - t\_0 \cdot \cos \left(\lambda_1 - \lambda_2\right)}\\
\end{array}
\end{array}
if lambda1 < -1.2e8Initial program 54.9%
add-sqr-sqrt29.4%
sqrt-unprod51.0%
pow251.0%
Applied egg-rr51.0%
Taylor expanded in phi1 around 0 40.8%
sin-diff76.4%
sub-neg76.4%
Applied egg-rr61.4%
+-commutative76.4%
distribute-rgt-neg-in76.4%
sin-neg76.4%
*-commutative76.4%
fma-define76.5%
sin-neg76.5%
cos-neg76.5%
*-commutative76.5%
cos-neg76.5%
Simplified61.5%
if -1.2e8 < lambda1 < 5.5999999999999999e-8Initial program 98.7%
Taylor expanded in lambda1 around 0 98.7%
cos-neg98.7%
Simplified98.7%
if 5.5999999999999999e-8 < lambda1 < 1.2200000000000001e119Initial program 49.5%
add-sqr-sqrt39.5%
sqrt-unprod46.1%
pow246.1%
Applied egg-rr46.1%
Taylor expanded in phi1 around 0 42.2%
sin-diff76.1%
Applied egg-rr68.7%
if 1.2200000000000001e119 < lambda1 Initial program 62.9%
Taylor expanded in lambda2 around 0 67.6%
Final simplification84.3%
(FPCore (lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (* (cos phi1) (sin phi2))) (t_1 (* (cos lambda2) (sin lambda1))))
(if (<= lambda1 -120000000.0)
(atan2
(* (fma (- (sin lambda2)) (cos lambda1) t_1) (cos phi2))
(sin phi2))
(if (<= lambda1 5.6e-8)
(atan2
(* (cos phi2) (sin (- lambda1 lambda2)))
(- t_0 (* (cos phi2) (* (cos lambda2) (sin phi1)))))
(if (<= lambda1 7.4e+117)
(atan2
(* (cos phi2) (- t_1 (* (sin lambda2) (cos lambda1))))
(sin phi2))
(atan2
(* (sin lambda1) (cos phi2))
(- t_0 (* (* (cos phi2) (sin phi1)) (cos (- lambda1 lambda2))))))))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = cos(phi1) * sin(phi2);
double t_1 = cos(lambda2) * sin(lambda1);
double tmp;
if (lambda1 <= -120000000.0) {
tmp = atan2((fma(-sin(lambda2), cos(lambda1), t_1) * cos(phi2)), sin(phi2));
} else if (lambda1 <= 5.6e-8) {
tmp = atan2((cos(phi2) * sin((lambda1 - lambda2))), (t_0 - (cos(phi2) * (cos(lambda2) * sin(phi1)))));
} else if (lambda1 <= 7.4e+117) {
tmp = atan2((cos(phi2) * (t_1 - (sin(lambda2) * cos(lambda1)))), sin(phi2));
} else {
tmp = atan2((sin(lambda1) * cos(phi2)), (t_0 - ((cos(phi2) * sin(phi1)) * cos((lambda1 - lambda2)))));
}
return tmp;
}
function code(lambda1, lambda2, phi1, phi2) t_0 = Float64(cos(phi1) * sin(phi2)) t_1 = Float64(cos(lambda2) * sin(lambda1)) tmp = 0.0 if (lambda1 <= -120000000.0) tmp = atan(Float64(fma(Float64(-sin(lambda2)), cos(lambda1), t_1) * cos(phi2)), sin(phi2)); elseif (lambda1 <= 5.6e-8) tmp = atan(Float64(cos(phi2) * sin(Float64(lambda1 - lambda2))), Float64(t_0 - Float64(cos(phi2) * Float64(cos(lambda2) * sin(phi1))))); elseif (lambda1 <= 7.4e+117) tmp = atan(Float64(cos(phi2) * Float64(t_1 - Float64(sin(lambda2) * cos(lambda1)))), sin(phi2)); else tmp = atan(Float64(sin(lambda1) * cos(phi2)), Float64(t_0 - Float64(Float64(cos(phi2) * sin(phi1)) * cos(Float64(lambda1 - lambda2))))); end return tmp end
code[lambda1_, lambda2_, phi1_, phi2_] := Block[{t$95$0 = N[(N[Cos[phi1], $MachinePrecision] * N[Sin[phi2], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[Cos[lambda2], $MachinePrecision] * N[Sin[lambda1], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[lambda1, -120000000.0], N[ArcTan[N[(N[((-N[Sin[lambda2], $MachinePrecision]) * N[Cos[lambda1], $MachinePrecision] + t$95$1), $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision] / N[Sin[phi2], $MachinePrecision]], $MachinePrecision], If[LessEqual[lambda1, 5.6e-8], N[ArcTan[N[(N[Cos[phi2], $MachinePrecision] * N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / N[(t$95$0 - N[(N[Cos[phi2], $MachinePrecision] * N[(N[Cos[lambda2], $MachinePrecision] * N[Sin[phi1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], If[LessEqual[lambda1, 7.4e+117], N[ArcTan[N[(N[Cos[phi2], $MachinePrecision] * N[(t$95$1 - N[(N[Sin[lambda2], $MachinePrecision] * N[Cos[lambda1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[Sin[phi2], $MachinePrecision]], $MachinePrecision], N[ArcTan[N[(N[Sin[lambda1], $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision] / N[(t$95$0 - N[(N[(N[Cos[phi2], $MachinePrecision] * N[Sin[phi1], $MachinePrecision]), $MachinePrecision] * N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \cos \phi_1 \cdot \sin \phi_2\\
t_1 := \cos \lambda_2 \cdot \sin \lambda_1\\
\mathbf{if}\;\lambda_1 \leq -120000000:\\
\;\;\;\;\tan^{-1}_* \frac{\mathsf{fma}\left(-\sin \lambda_2, \cos \lambda_1, t\_1\right) \cdot \cos \phi_2}{\sin \phi_2}\\
\mathbf{elif}\;\lambda_1 \leq 5.6 \cdot 10^{-8}:\\
\;\;\;\;\tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{t\_0 - \cos \phi_2 \cdot \left(\cos \lambda_2 \cdot \sin \phi_1\right)}\\
\mathbf{elif}\;\lambda_1 \leq 7.4 \cdot 10^{+117}:\\
\;\;\;\;\tan^{-1}_* \frac{\cos \phi_2 \cdot \left(t\_1 - \sin \lambda_2 \cdot \cos \lambda_1\right)}{\sin \phi_2}\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1}_* \frac{\sin \lambda_1 \cdot \cos \phi_2}{t\_0 - \left(\cos \phi_2 \cdot \sin \phi_1\right) \cdot \cos \left(\lambda_1 - \lambda_2\right)}\\
\end{array}
\end{array}
if lambda1 < -1.2e8Initial program 54.9%
add-sqr-sqrt29.4%
sqrt-unprod51.0%
pow251.0%
Applied egg-rr51.0%
Taylor expanded in phi1 around 0 40.8%
sin-diff76.4%
sub-neg76.4%
Applied egg-rr61.4%
+-commutative76.4%
distribute-rgt-neg-in76.4%
sin-neg76.4%
*-commutative76.4%
fma-define76.5%
sin-neg76.5%
cos-neg76.5%
*-commutative76.5%
cos-neg76.5%
Simplified61.5%
if -1.2e8 < lambda1 < 5.5999999999999999e-8Initial program 98.7%
*-commutative98.7%
associate-*l*98.7%
Simplified98.7%
Taylor expanded in lambda1 around 0 98.7%
cos-neg98.7%
Simplified98.7%
if 5.5999999999999999e-8 < lambda1 < 7.3999999999999997e117Initial program 49.5%
add-sqr-sqrt39.5%
sqrt-unprod46.1%
pow246.1%
Applied egg-rr46.1%
Taylor expanded in phi1 around 0 42.2%
sin-diff76.1%
Applied egg-rr68.7%
if 7.3999999999999997e117 < lambda1 Initial program 62.9%
Taylor expanded in lambda2 around 0 67.6%
Final simplification84.3%
(FPCore (lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0
(atan2
(* (cos phi2) (sin (- lambda2)))
(-
(* (cos phi1) (sin phi2))
(* (cos phi2) (* (cos lambda2) (sin phi1)))))))
(if (<= lambda2 -3.6e+28)
t_0
(if (<= lambda2 1.26e+36)
(atan2
(* (cos phi2) (sin (- lambda1 lambda2)))
(- (sin phi2) (* (cos phi2) (* (sin phi1) (cos (- lambda1 lambda2))))))
(if (<= lambda2 1.8e+106)
(atan2
(*
(fma
(- (sin lambda2))
(cos lambda1)
(* (cos lambda2) (sin lambda1)))
(cos phi2))
(sin phi2))
t_0)))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = atan2((cos(phi2) * sin(-lambda2)), ((cos(phi1) * sin(phi2)) - (cos(phi2) * (cos(lambda2) * sin(phi1)))));
double tmp;
if (lambda2 <= -3.6e+28) {
tmp = t_0;
} else if (lambda2 <= 1.26e+36) {
tmp = atan2((cos(phi2) * sin((lambda1 - lambda2))), (sin(phi2) - (cos(phi2) * (sin(phi1) * cos((lambda1 - lambda2))))));
} else if (lambda2 <= 1.8e+106) {
tmp = atan2((fma(-sin(lambda2), cos(lambda1), (cos(lambda2) * sin(lambda1))) * cos(phi2)), sin(phi2));
} else {
tmp = t_0;
}
return tmp;
}
function code(lambda1, lambda2, phi1, phi2) t_0 = atan(Float64(cos(phi2) * sin(Float64(-lambda2))), Float64(Float64(cos(phi1) * sin(phi2)) - Float64(cos(phi2) * Float64(cos(lambda2) * sin(phi1))))) tmp = 0.0 if (lambda2 <= -3.6e+28) tmp = t_0; elseif (lambda2 <= 1.26e+36) tmp = atan(Float64(cos(phi2) * sin(Float64(lambda1 - lambda2))), Float64(sin(phi2) - Float64(cos(phi2) * Float64(sin(phi1) * cos(Float64(lambda1 - lambda2)))))); elseif (lambda2 <= 1.8e+106) tmp = atan(Float64(fma(Float64(-sin(lambda2)), cos(lambda1), Float64(cos(lambda2) * sin(lambda1))) * cos(phi2)), sin(phi2)); else tmp = t_0; end return tmp end
code[lambda1_, lambda2_, phi1_, phi2_] := Block[{t$95$0 = N[ArcTan[N[(N[Cos[phi2], $MachinePrecision] * N[Sin[(-lambda2)], $MachinePrecision]), $MachinePrecision] / N[(N[(N[Cos[phi1], $MachinePrecision] * N[Sin[phi2], $MachinePrecision]), $MachinePrecision] - N[(N[Cos[phi2], $MachinePrecision] * N[(N[Cos[lambda2], $MachinePrecision] * N[Sin[phi1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]}, If[LessEqual[lambda2, -3.6e+28], t$95$0, If[LessEqual[lambda2, 1.26e+36], N[ArcTan[N[(N[Cos[phi2], $MachinePrecision] * N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / N[(N[Sin[phi2], $MachinePrecision] - N[(N[Cos[phi2], $MachinePrecision] * N[(N[Sin[phi1], $MachinePrecision] * N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], If[LessEqual[lambda2, 1.8e+106], N[ArcTan[N[(N[((-N[Sin[lambda2], $MachinePrecision]) * N[Cos[lambda1], $MachinePrecision] + N[(N[Cos[lambda2], $MachinePrecision] * N[Sin[lambda1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision] / N[Sin[phi2], $MachinePrecision]], $MachinePrecision], t$95$0]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(-\lambda_2\right)}{\cos \phi_1 \cdot \sin \phi_2 - \cos \phi_2 \cdot \left(\cos \lambda_2 \cdot \sin \phi_1\right)}\\
\mathbf{if}\;\lambda_2 \leq -3.6 \cdot 10^{+28}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;\lambda_2 \leq 1.26 \cdot 10^{+36}:\\
\;\;\;\;\tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\sin \phi_2 - \cos \phi_2 \cdot \left(\sin \phi_1 \cdot \cos \left(\lambda_1 - \lambda_2\right)\right)}\\
\mathbf{elif}\;\lambda_2 \leq 1.8 \cdot 10^{+106}:\\
\;\;\;\;\tan^{-1}_* \frac{\mathsf{fma}\left(-\sin \lambda_2, \cos \lambda_1, \cos \lambda_2 \cdot \sin \lambda_1\right) \cdot \cos \phi_2}{\sin \phi_2}\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if lambda2 < -3.5999999999999999e28 or 1.8e106 < lambda2 Initial program 64.2%
*-commutative64.2%
associate-*l*64.2%
Simplified64.2%
Taylor expanded in lambda1 around 0 64.3%
cos-neg64.3%
Simplified64.3%
Taylor expanded in lambda1 around 0 66.4%
if -3.5999999999999999e28 < lambda2 < 1.25999999999999994e36Initial program 97.9%
*-commutative97.9%
associate-*l*97.9%
Simplified97.9%
Taylor expanded in phi1 around 0 82.5%
if 1.25999999999999994e36 < lambda2 < 1.8e106Initial program 38.4%
add-sqr-sqrt14.8%
sqrt-unprod34.3%
pow234.3%
Applied egg-rr34.3%
Taylor expanded in phi1 around 0 33.5%
sin-diff74.2%
sub-neg74.2%
Applied egg-rr68.6%
+-commutative74.2%
distribute-rgt-neg-in74.2%
sin-neg74.2%
*-commutative74.2%
fma-define74.3%
sin-neg74.3%
cos-neg74.3%
*-commutative74.3%
cos-neg74.3%
Simplified68.7%
Final simplification75.4%
(FPCore (lambda1 lambda2 phi1 phi2)
:precision binary64
(if (or (<= phi1 -3.1e-128) (not (<= phi1 6e-13)))
(atan2
(* (cos phi2) (sin (- lambda1 lambda2)))
(-
(* (cos phi1) (sin phi2))
(* (sin phi1) (* (cos phi2) (cos (- lambda1 lambda2))))))
(atan2
(*
(fma (- (sin lambda2)) (cos lambda1) (* (cos lambda2) (sin lambda1)))
(cos phi2))
(sin phi2))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
double tmp;
if ((phi1 <= -3.1e-128) || !(phi1 <= 6e-13)) {
tmp = atan2((cos(phi2) * sin((lambda1 - lambda2))), ((cos(phi1) * sin(phi2)) - (sin(phi1) * (cos(phi2) * cos((lambda1 - lambda2))))));
} else {
tmp = atan2((fma(-sin(lambda2), cos(lambda1), (cos(lambda2) * sin(lambda1))) * cos(phi2)), sin(phi2));
}
return tmp;
}
function code(lambda1, lambda2, phi1, phi2) tmp = 0.0 if ((phi1 <= -3.1e-128) || !(phi1 <= 6e-13)) tmp = atan(Float64(cos(phi2) * sin(Float64(lambda1 - lambda2))), Float64(Float64(cos(phi1) * sin(phi2)) - Float64(sin(phi1) * Float64(cos(phi2) * cos(Float64(lambda1 - lambda2)))))); else tmp = atan(Float64(fma(Float64(-sin(lambda2)), cos(lambda1), Float64(cos(lambda2) * sin(lambda1))) * cos(phi2)), sin(phi2)); end return tmp end
code[lambda1_, lambda2_, phi1_, phi2_] := If[Or[LessEqual[phi1, -3.1e-128], N[Not[LessEqual[phi1, 6e-13]], $MachinePrecision]], N[ArcTan[N[(N[Cos[phi2], $MachinePrecision] * N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / N[(N[(N[Cos[phi1], $MachinePrecision] * N[Sin[phi2], $MachinePrecision]), $MachinePrecision] - N[(N[Sin[phi1], $MachinePrecision] * N[(N[Cos[phi2], $MachinePrecision] * N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], N[ArcTan[N[(N[((-N[Sin[lambda2], $MachinePrecision]) * N[Cos[lambda1], $MachinePrecision] + N[(N[Cos[lambda2], $MachinePrecision] * N[Sin[lambda1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision] / N[Sin[phi2], $MachinePrecision]], $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\phi_1 \leq -3.1 \cdot 10^{-128} \lor \neg \left(\phi_1 \leq 6 \cdot 10^{-13}\right):\\
\;\;\;\;\tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 \cdot \sin \phi_2 - \sin \phi_1 \cdot \left(\cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1}_* \frac{\mathsf{fma}\left(-\sin \lambda_2, \cos \lambda_1, \cos \lambda_2 \cdot \sin \lambda_1\right) \cdot \cos \phi_2}{\sin \phi_2}\\
\end{array}
\end{array}
if phi1 < -3.10000000000000003e-128 or 5.99999999999999968e-13 < phi1 Initial program 79.7%
Taylor expanded in phi1 around inf 79.6%
associate-*r*82.7%
*-commutative82.7%
Simplified79.6%
if -3.10000000000000003e-128 < phi1 < 5.99999999999999968e-13Initial program 82.0%
add-sqr-sqrt59.0%
sqrt-unprod82.0%
pow282.0%
Applied egg-rr82.0%
Taylor expanded in phi1 around 0 82.0%
sin-diff99.9%
sub-neg99.9%
Applied egg-rr99.9%
+-commutative99.9%
distribute-rgt-neg-in99.9%
sin-neg99.9%
*-commutative99.9%
fma-define99.9%
sin-neg99.9%
cos-neg99.9%
*-commutative99.9%
cos-neg99.9%
Simplified99.9%
Final simplification87.2%
(FPCore (lambda1 lambda2 phi1 phi2)
:precision binary64
(if (or (<= phi1 -3.1e-128) (not (<= phi1 4.8e-14)))
(atan2
(* (cos phi2) (sin (- lambda1 lambda2)))
(-
(* (cos phi1) (sin phi2))
(* (cos phi2) (* (sin phi1) (cos (- lambda1 lambda2))))))
(atan2
(*
(fma (- (sin lambda2)) (cos lambda1) (* (cos lambda2) (sin lambda1)))
(cos phi2))
(sin phi2))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
double tmp;
if ((phi1 <= -3.1e-128) || !(phi1 <= 4.8e-14)) {
tmp = atan2((cos(phi2) * sin((lambda1 - lambda2))), ((cos(phi1) * sin(phi2)) - (cos(phi2) * (sin(phi1) * cos((lambda1 - lambda2))))));
} else {
tmp = atan2((fma(-sin(lambda2), cos(lambda1), (cos(lambda2) * sin(lambda1))) * cos(phi2)), sin(phi2));
}
return tmp;
}
function code(lambda1, lambda2, phi1, phi2) tmp = 0.0 if ((phi1 <= -3.1e-128) || !(phi1 <= 4.8e-14)) tmp = atan(Float64(cos(phi2) * sin(Float64(lambda1 - lambda2))), Float64(Float64(cos(phi1) * sin(phi2)) - Float64(cos(phi2) * Float64(sin(phi1) * cos(Float64(lambda1 - lambda2)))))); else tmp = atan(Float64(fma(Float64(-sin(lambda2)), cos(lambda1), Float64(cos(lambda2) * sin(lambda1))) * cos(phi2)), sin(phi2)); end return tmp end
code[lambda1_, lambda2_, phi1_, phi2_] := If[Or[LessEqual[phi1, -3.1e-128], N[Not[LessEqual[phi1, 4.8e-14]], $MachinePrecision]], N[ArcTan[N[(N[Cos[phi2], $MachinePrecision] * N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / N[(N[(N[Cos[phi1], $MachinePrecision] * N[Sin[phi2], $MachinePrecision]), $MachinePrecision] - N[(N[Cos[phi2], $MachinePrecision] * N[(N[Sin[phi1], $MachinePrecision] * N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], N[ArcTan[N[(N[((-N[Sin[lambda2], $MachinePrecision]) * N[Cos[lambda1], $MachinePrecision] + N[(N[Cos[lambda2], $MachinePrecision] * N[Sin[lambda1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision] / N[Sin[phi2], $MachinePrecision]], $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\phi_1 \leq -3.1 \cdot 10^{-128} \lor \neg \left(\phi_1 \leq 4.8 \cdot 10^{-14}\right):\\
\;\;\;\;\tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 \cdot \sin \phi_2 - \cos \phi_2 \cdot \left(\sin \phi_1 \cdot \cos \left(\lambda_1 - \lambda_2\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1}_* \frac{\mathsf{fma}\left(-\sin \lambda_2, \cos \lambda_1, \cos \lambda_2 \cdot \sin \lambda_1\right) \cdot \cos \phi_2}{\sin \phi_2}\\
\end{array}
\end{array}
if phi1 < -3.10000000000000003e-128 or 4.8e-14 < phi1 Initial program 79.7%
*-commutative79.7%
associate-*l*79.6%
Simplified79.6%
if -3.10000000000000003e-128 < phi1 < 4.8e-14Initial program 82.0%
add-sqr-sqrt59.0%
sqrt-unprod82.0%
pow282.0%
Applied egg-rr82.0%
Taylor expanded in phi1 around 0 82.0%
sin-diff99.9%
sub-neg99.9%
Applied egg-rr99.9%
+-commutative99.9%
distribute-rgt-neg-in99.9%
sin-neg99.9%
*-commutative99.9%
fma-define99.9%
sin-neg99.9%
cos-neg99.9%
*-commutative99.9%
cos-neg99.9%
Simplified99.9%
Final simplification87.2%
(FPCore (lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (* (cos phi1) (sin phi2)))
(t_1 (cos (- lambda1 lambda2)))
(t_2 (* (cos phi2) (sin (- lambda1 lambda2)))))
(if (<= phi1 -2.35e-129)
(atan2 t_2 (- t_0 (* (* (cos phi2) (sin phi1)) t_1)))
(if (<= phi1 1.25e-17)
(atan2
(*
(fma (- (sin lambda2)) (cos lambda1) (* (cos lambda2) (sin lambda1)))
(cos phi2))
(sin phi2))
(atan2 t_2 (- t_0 (* (sin phi1) (* (cos phi2) t_1))))))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = cos(phi1) * sin(phi2);
double t_1 = cos((lambda1 - lambda2));
double t_2 = cos(phi2) * sin((lambda1 - lambda2));
double tmp;
if (phi1 <= -2.35e-129) {
tmp = atan2(t_2, (t_0 - ((cos(phi2) * sin(phi1)) * t_1)));
} else if (phi1 <= 1.25e-17) {
tmp = atan2((fma(-sin(lambda2), cos(lambda1), (cos(lambda2) * sin(lambda1))) * cos(phi2)), sin(phi2));
} else {
tmp = atan2(t_2, (t_0 - (sin(phi1) * (cos(phi2) * t_1))));
}
return tmp;
}
function code(lambda1, lambda2, phi1, phi2) t_0 = Float64(cos(phi1) * sin(phi2)) t_1 = cos(Float64(lambda1 - lambda2)) t_2 = Float64(cos(phi2) * sin(Float64(lambda1 - lambda2))) tmp = 0.0 if (phi1 <= -2.35e-129) tmp = atan(t_2, Float64(t_0 - Float64(Float64(cos(phi2) * sin(phi1)) * t_1))); elseif (phi1 <= 1.25e-17) tmp = atan(Float64(fma(Float64(-sin(lambda2)), cos(lambda1), Float64(cos(lambda2) * sin(lambda1))) * cos(phi2)), sin(phi2)); else tmp = atan(t_2, Float64(t_0 - Float64(sin(phi1) * Float64(cos(phi2) * t_1)))); end return tmp end
code[lambda1_, lambda2_, phi1_, phi2_] := Block[{t$95$0 = N[(N[Cos[phi1], $MachinePrecision] * N[Sin[phi2], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$2 = N[(N[Cos[phi2], $MachinePrecision] * N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[phi1, -2.35e-129], N[ArcTan[t$95$2 / N[(t$95$0 - N[(N[(N[Cos[phi2], $MachinePrecision] * N[Sin[phi1], $MachinePrecision]), $MachinePrecision] * t$95$1), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], If[LessEqual[phi1, 1.25e-17], N[ArcTan[N[(N[((-N[Sin[lambda2], $MachinePrecision]) * N[Cos[lambda1], $MachinePrecision] + N[(N[Cos[lambda2], $MachinePrecision] * N[Sin[lambda1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision] / N[Sin[phi2], $MachinePrecision]], $MachinePrecision], N[ArcTan[t$95$2 / N[(t$95$0 - N[(N[Sin[phi1], $MachinePrecision] * N[(N[Cos[phi2], $MachinePrecision] * t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \cos \phi_1 \cdot \sin \phi_2\\
t_1 := \cos \left(\lambda_1 - \lambda_2\right)\\
t_2 := \cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)\\
\mathbf{if}\;\phi_1 \leq -2.35 \cdot 10^{-129}:\\
\;\;\;\;\tan^{-1}_* \frac{t\_2}{t\_0 - \left(\cos \phi_2 \cdot \sin \phi_1\right) \cdot t\_1}\\
\mathbf{elif}\;\phi_1 \leq 1.25 \cdot 10^{-17}:\\
\;\;\;\;\tan^{-1}_* \frac{\mathsf{fma}\left(-\sin \lambda_2, \cos \lambda_1, \cos \lambda_2 \cdot \sin \lambda_1\right) \cdot \cos \phi_2}{\sin \phi_2}\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1}_* \frac{t\_2}{t\_0 - \sin \phi_1 \cdot \left(\cos \phi_2 \cdot t\_1\right)}\\
\end{array}
\end{array}
if phi1 < -2.3500000000000001e-129Initial program 79.4%
if -2.3500000000000001e-129 < phi1 < 1.25e-17Initial program 82.0%
add-sqr-sqrt59.0%
sqrt-unprod82.0%
pow282.0%
Applied egg-rr82.0%
Taylor expanded in phi1 around 0 82.0%
sin-diff99.9%
sub-neg99.9%
Applied egg-rr99.9%
+-commutative99.9%
distribute-rgt-neg-in99.9%
sin-neg99.9%
*-commutative99.9%
fma-define99.9%
sin-neg99.9%
cos-neg99.9%
*-commutative99.9%
cos-neg99.9%
Simplified99.9%
if 1.25e-17 < phi1 Initial program 80.1%
Taylor expanded in phi1 around inf 80.1%
associate-*r*83.0%
*-commutative83.0%
Simplified80.1%
Final simplification87.3%
(FPCore (lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (* (sin phi1) (cos (- lambda1 lambda2))))
(t_1 (* (cos phi2) (sin (- lambda1 lambda2)))))
(if (<= phi1 -8.6e-129)
(atan2 t_1 (- (sin phi2) (* (cos phi2) t_0)))
(if (<= phi1 3.2e-9)
(atan2
(*
(fma (- (sin lambda2)) (cos lambda1) (* (cos lambda2) (sin lambda1)))
(cos phi2))
(sin phi2))
(atan2 t_1 (- (* (cos phi1) (sin phi2)) t_0))))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = sin(phi1) * cos((lambda1 - lambda2));
double t_1 = cos(phi2) * sin((lambda1 - lambda2));
double tmp;
if (phi1 <= -8.6e-129) {
tmp = atan2(t_1, (sin(phi2) - (cos(phi2) * t_0)));
} else if (phi1 <= 3.2e-9) {
tmp = atan2((fma(-sin(lambda2), cos(lambda1), (cos(lambda2) * sin(lambda1))) * cos(phi2)), sin(phi2));
} else {
tmp = atan2(t_1, ((cos(phi1) * sin(phi2)) - t_0));
}
return tmp;
}
function code(lambda1, lambda2, phi1, phi2) t_0 = Float64(sin(phi1) * cos(Float64(lambda1 - lambda2))) t_1 = Float64(cos(phi2) * sin(Float64(lambda1 - lambda2))) tmp = 0.0 if (phi1 <= -8.6e-129) tmp = atan(t_1, Float64(sin(phi2) - Float64(cos(phi2) * t_0))); elseif (phi1 <= 3.2e-9) tmp = atan(Float64(fma(Float64(-sin(lambda2)), cos(lambda1), Float64(cos(lambda2) * sin(lambda1))) * cos(phi2)), sin(phi2)); else tmp = atan(t_1, Float64(Float64(cos(phi1) * sin(phi2)) - t_0)); end return tmp end
code[lambda1_, lambda2_, phi1_, phi2_] := Block[{t$95$0 = N[(N[Sin[phi1], $MachinePrecision] * N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[Cos[phi2], $MachinePrecision] * N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[phi1, -8.6e-129], N[ArcTan[t$95$1 / N[(N[Sin[phi2], $MachinePrecision] - N[(N[Cos[phi2], $MachinePrecision] * t$95$0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], If[LessEqual[phi1, 3.2e-9], N[ArcTan[N[(N[((-N[Sin[lambda2], $MachinePrecision]) * N[Cos[lambda1], $MachinePrecision] + N[(N[Cos[lambda2], $MachinePrecision] * N[Sin[lambda1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision] / N[Sin[phi2], $MachinePrecision]], $MachinePrecision], N[ArcTan[t$95$1 / N[(N[(N[Cos[phi1], $MachinePrecision] * N[Sin[phi2], $MachinePrecision]), $MachinePrecision] - t$95$0), $MachinePrecision]], $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sin \phi_1 \cdot \cos \left(\lambda_1 - \lambda_2\right)\\
t_1 := \cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)\\
\mathbf{if}\;\phi_1 \leq -8.6 \cdot 10^{-129}:\\
\;\;\;\;\tan^{-1}_* \frac{t\_1}{\sin \phi_2 - \cos \phi_2 \cdot t\_0}\\
\mathbf{elif}\;\phi_1 \leq 3.2 \cdot 10^{-9}:\\
\;\;\;\;\tan^{-1}_* \frac{\mathsf{fma}\left(-\sin \lambda_2, \cos \lambda_1, \cos \lambda_2 \cdot \sin \lambda_1\right) \cdot \cos \phi_2}{\sin \phi_2}\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1}_* \frac{t\_1}{\cos \phi_1 \cdot \sin \phi_2 - t\_0}\\
\end{array}
\end{array}
if phi1 < -8.59999999999999962e-129Initial program 79.4%
*-commutative79.4%
associate-*l*79.4%
Simplified79.4%
Taylor expanded in phi1 around 0 59.3%
if -8.59999999999999962e-129 < phi1 < 3.20000000000000012e-9Initial program 82.0%
add-sqr-sqrt59.0%
sqrt-unprod82.0%
pow282.0%
Applied egg-rr82.0%
Taylor expanded in phi1 around 0 82.0%
sin-diff99.9%
sub-neg99.9%
Applied egg-rr99.9%
+-commutative99.9%
distribute-rgt-neg-in99.9%
sin-neg99.9%
*-commutative99.9%
fma-define99.9%
sin-neg99.9%
cos-neg99.9%
*-commutative99.9%
cos-neg99.9%
Simplified99.9%
if 3.20000000000000012e-9 < phi1 Initial program 80.1%
add-sqr-sqrt61.1%
sqrt-unprod65.9%
pow265.9%
Applied egg-rr65.9%
Taylor expanded in phi2 around 0 54.2%
Final simplification73.4%
(FPCore (lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (* (cos phi2) (sin (- lambda1 lambda2))))
(t_1
(atan2
t_0
(-
(sin phi2)
(* (cos phi2) (* (sin phi1) (cos (- lambda1 lambda2))))))))
(if (<= phi1 -3.1e-128)
t_1
(if (<= phi1 1.5e-16)
(atan2
(*
(cos phi2)
(- (* (cos lambda2) (sin lambda1)) (* (sin lambda2) (cos lambda1))))
(sin phi2))
(if (<= phi1 1.56e+224)
t_1
(atan2
t_0
(- (* (cos phi1) (sin phi2)) (* (cos phi2) (sin phi1)))))))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = cos(phi2) * sin((lambda1 - lambda2));
double t_1 = atan2(t_0, (sin(phi2) - (cos(phi2) * (sin(phi1) * cos((lambda1 - lambda2))))));
double tmp;
if (phi1 <= -3.1e-128) {
tmp = t_1;
} else if (phi1 <= 1.5e-16) {
tmp = atan2((cos(phi2) * ((cos(lambda2) * sin(lambda1)) - (sin(lambda2) * cos(lambda1)))), sin(phi2));
} else if (phi1 <= 1.56e+224) {
tmp = t_1;
} else {
tmp = atan2(t_0, ((cos(phi1) * sin(phi2)) - (cos(phi2) * sin(phi1))));
}
return tmp;
}
real(8) function code(lambda1, lambda2, phi1, phi2)
real(8), intent (in) :: lambda1
real(8), intent (in) :: lambda2
real(8), intent (in) :: phi1
real(8), intent (in) :: phi2
real(8) :: t_0
real(8) :: t_1
real(8) :: tmp
t_0 = cos(phi2) * sin((lambda1 - lambda2))
t_1 = atan2(t_0, (sin(phi2) - (cos(phi2) * (sin(phi1) * cos((lambda1 - lambda2))))))
if (phi1 <= (-3.1d-128)) then
tmp = t_1
else if (phi1 <= 1.5d-16) then
tmp = atan2((cos(phi2) * ((cos(lambda2) * sin(lambda1)) - (sin(lambda2) * cos(lambda1)))), sin(phi2))
else if (phi1 <= 1.56d+224) then
tmp = t_1
else
tmp = atan2(t_0, ((cos(phi1) * sin(phi2)) - (cos(phi2) * sin(phi1))))
end if
code = tmp
end function
public static double code(double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = Math.cos(phi2) * Math.sin((lambda1 - lambda2));
double t_1 = Math.atan2(t_0, (Math.sin(phi2) - (Math.cos(phi2) * (Math.sin(phi1) * Math.cos((lambda1 - lambda2))))));
double tmp;
if (phi1 <= -3.1e-128) {
tmp = t_1;
} else if (phi1 <= 1.5e-16) {
tmp = Math.atan2((Math.cos(phi2) * ((Math.cos(lambda2) * Math.sin(lambda1)) - (Math.sin(lambda2) * Math.cos(lambda1)))), Math.sin(phi2));
} else if (phi1 <= 1.56e+224) {
tmp = t_1;
} else {
tmp = Math.atan2(t_0, ((Math.cos(phi1) * Math.sin(phi2)) - (Math.cos(phi2) * Math.sin(phi1))));
}
return tmp;
}
def code(lambda1, lambda2, phi1, phi2): t_0 = math.cos(phi2) * math.sin((lambda1 - lambda2)) t_1 = math.atan2(t_0, (math.sin(phi2) - (math.cos(phi2) * (math.sin(phi1) * math.cos((lambda1 - lambda2)))))) tmp = 0 if phi1 <= -3.1e-128: tmp = t_1 elif phi1 <= 1.5e-16: tmp = math.atan2((math.cos(phi2) * ((math.cos(lambda2) * math.sin(lambda1)) - (math.sin(lambda2) * math.cos(lambda1)))), math.sin(phi2)) elif phi1 <= 1.56e+224: tmp = t_1 else: tmp = math.atan2(t_0, ((math.cos(phi1) * math.sin(phi2)) - (math.cos(phi2) * math.sin(phi1)))) return tmp
function code(lambda1, lambda2, phi1, phi2) t_0 = Float64(cos(phi2) * sin(Float64(lambda1 - lambda2))) t_1 = atan(t_0, Float64(sin(phi2) - Float64(cos(phi2) * Float64(sin(phi1) * cos(Float64(lambda1 - lambda2)))))) tmp = 0.0 if (phi1 <= -3.1e-128) tmp = t_1; elseif (phi1 <= 1.5e-16) tmp = atan(Float64(cos(phi2) * Float64(Float64(cos(lambda2) * sin(lambda1)) - Float64(sin(lambda2) * cos(lambda1)))), sin(phi2)); elseif (phi1 <= 1.56e+224) tmp = t_1; else tmp = atan(t_0, Float64(Float64(cos(phi1) * sin(phi2)) - Float64(cos(phi2) * sin(phi1)))); end return tmp end
function tmp_2 = code(lambda1, lambda2, phi1, phi2) t_0 = cos(phi2) * sin((lambda1 - lambda2)); t_1 = atan2(t_0, (sin(phi2) - (cos(phi2) * (sin(phi1) * cos((lambda1 - lambda2)))))); tmp = 0.0; if (phi1 <= -3.1e-128) tmp = t_1; elseif (phi1 <= 1.5e-16) tmp = atan2((cos(phi2) * ((cos(lambda2) * sin(lambda1)) - (sin(lambda2) * cos(lambda1)))), sin(phi2)); elseif (phi1 <= 1.56e+224) tmp = t_1; else tmp = atan2(t_0, ((cos(phi1) * sin(phi2)) - (cos(phi2) * sin(phi1)))); end tmp_2 = tmp; end
code[lambda1_, lambda2_, phi1_, phi2_] := Block[{t$95$0 = N[(N[Cos[phi2], $MachinePrecision] * N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[ArcTan[t$95$0 / N[(N[Sin[phi2], $MachinePrecision] - N[(N[Cos[phi2], $MachinePrecision] * N[(N[Sin[phi1], $MachinePrecision] * N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]}, If[LessEqual[phi1, -3.1e-128], t$95$1, If[LessEqual[phi1, 1.5e-16], N[ArcTan[N[(N[Cos[phi2], $MachinePrecision] * N[(N[(N[Cos[lambda2], $MachinePrecision] * N[Sin[lambda1], $MachinePrecision]), $MachinePrecision] - N[(N[Sin[lambda2], $MachinePrecision] * N[Cos[lambda1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[Sin[phi2], $MachinePrecision]], $MachinePrecision], If[LessEqual[phi1, 1.56e+224], t$95$1, N[ArcTan[t$95$0 / N[(N[(N[Cos[phi1], $MachinePrecision] * N[Sin[phi2], $MachinePrecision]), $MachinePrecision] - N[(N[Cos[phi2], $MachinePrecision] * N[Sin[phi1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)\\
t_1 := \tan^{-1}_* \frac{t\_0}{\sin \phi_2 - \cos \phi_2 \cdot \left(\sin \phi_1 \cdot \cos \left(\lambda_1 - \lambda_2\right)\right)}\\
\mathbf{if}\;\phi_1 \leq -3.1 \cdot 10^{-128}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;\phi_1 \leq 1.5 \cdot 10^{-16}:\\
\;\;\;\;\tan^{-1}_* \frac{\cos \phi_2 \cdot \left(\cos \lambda_2 \cdot \sin \lambda_1 - \sin \lambda_2 \cdot \cos \lambda_1\right)}{\sin \phi_2}\\
\mathbf{elif}\;\phi_1 \leq 1.56 \cdot 10^{+224}:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1}_* \frac{t\_0}{\cos \phi_1 \cdot \sin \phi_2 - \cos \phi_2 \cdot \sin \phi_1}\\
\end{array}
\end{array}
if phi1 < -3.10000000000000003e-128 or 1.49999999999999997e-16 < phi1 < 1.5599999999999999e224Initial program 80.3%
*-commutative80.3%
associate-*l*80.3%
Simplified80.3%
Taylor expanded in phi1 around 0 59.6%
if -3.10000000000000003e-128 < phi1 < 1.49999999999999997e-16Initial program 82.0%
add-sqr-sqrt59.0%
sqrt-unprod82.0%
pow282.0%
Applied egg-rr82.0%
Taylor expanded in phi1 around 0 82.0%
sin-diff99.9%
Applied egg-rr99.9%
if 1.5599999999999999e224 < phi1 Initial program 74.4%
*-commutative74.4%
associate-*l*74.4%
Simplified74.4%
Taylor expanded in lambda1 around 0 65.2%
cos-neg65.2%
Simplified65.2%
Taylor expanded in lambda2 around 0 61.5%
Final simplification74.9%
(FPCore (lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (* (sin phi1) (cos (- lambda1 lambda2))))
(t_1 (* (cos phi2) (sin (- lambda1 lambda2)))))
(if (<= phi1 -3e-128)
(atan2 t_1 (- (sin phi2) (* (cos phi2) t_0)))
(if (<= phi1 5e-13)
(atan2
(*
(cos phi2)
(- (* (cos lambda2) (sin lambda1)) (* (sin lambda2) (cos lambda1))))
(sin phi2))
(atan2 t_1 (- (* (cos phi1) (sin phi2)) t_0))))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = sin(phi1) * cos((lambda1 - lambda2));
double t_1 = cos(phi2) * sin((lambda1 - lambda2));
double tmp;
if (phi1 <= -3e-128) {
tmp = atan2(t_1, (sin(phi2) - (cos(phi2) * t_0)));
} else if (phi1 <= 5e-13) {
tmp = atan2((cos(phi2) * ((cos(lambda2) * sin(lambda1)) - (sin(lambda2) * cos(lambda1)))), sin(phi2));
} else {
tmp = atan2(t_1, ((cos(phi1) * sin(phi2)) - t_0));
}
return tmp;
}
real(8) function code(lambda1, lambda2, phi1, phi2)
real(8), intent (in) :: lambda1
real(8), intent (in) :: lambda2
real(8), intent (in) :: phi1
real(8), intent (in) :: phi2
real(8) :: t_0
real(8) :: t_1
real(8) :: tmp
t_0 = sin(phi1) * cos((lambda1 - lambda2))
t_1 = cos(phi2) * sin((lambda1 - lambda2))
if (phi1 <= (-3d-128)) then
tmp = atan2(t_1, (sin(phi2) - (cos(phi2) * t_0)))
else if (phi1 <= 5d-13) then
tmp = atan2((cos(phi2) * ((cos(lambda2) * sin(lambda1)) - (sin(lambda2) * cos(lambda1)))), sin(phi2))
else
tmp = atan2(t_1, ((cos(phi1) * sin(phi2)) - t_0))
end if
code = tmp
end function
public static double code(double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = Math.sin(phi1) * Math.cos((lambda1 - lambda2));
double t_1 = Math.cos(phi2) * Math.sin((lambda1 - lambda2));
double tmp;
if (phi1 <= -3e-128) {
tmp = Math.atan2(t_1, (Math.sin(phi2) - (Math.cos(phi2) * t_0)));
} else if (phi1 <= 5e-13) {
tmp = Math.atan2((Math.cos(phi2) * ((Math.cos(lambda2) * Math.sin(lambda1)) - (Math.sin(lambda2) * Math.cos(lambda1)))), Math.sin(phi2));
} else {
tmp = Math.atan2(t_1, ((Math.cos(phi1) * Math.sin(phi2)) - t_0));
}
return tmp;
}
def code(lambda1, lambda2, phi1, phi2): t_0 = math.sin(phi1) * math.cos((lambda1 - lambda2)) t_1 = math.cos(phi2) * math.sin((lambda1 - lambda2)) tmp = 0 if phi1 <= -3e-128: tmp = math.atan2(t_1, (math.sin(phi2) - (math.cos(phi2) * t_0))) elif phi1 <= 5e-13: tmp = math.atan2((math.cos(phi2) * ((math.cos(lambda2) * math.sin(lambda1)) - (math.sin(lambda2) * math.cos(lambda1)))), math.sin(phi2)) else: tmp = math.atan2(t_1, ((math.cos(phi1) * math.sin(phi2)) - t_0)) return tmp
function code(lambda1, lambda2, phi1, phi2) t_0 = Float64(sin(phi1) * cos(Float64(lambda1 - lambda2))) t_1 = Float64(cos(phi2) * sin(Float64(lambda1 - lambda2))) tmp = 0.0 if (phi1 <= -3e-128) tmp = atan(t_1, Float64(sin(phi2) - Float64(cos(phi2) * t_0))); elseif (phi1 <= 5e-13) tmp = atan(Float64(cos(phi2) * Float64(Float64(cos(lambda2) * sin(lambda1)) - Float64(sin(lambda2) * cos(lambda1)))), sin(phi2)); else tmp = atan(t_1, Float64(Float64(cos(phi1) * sin(phi2)) - t_0)); end return tmp end
function tmp_2 = code(lambda1, lambda2, phi1, phi2) t_0 = sin(phi1) * cos((lambda1 - lambda2)); t_1 = cos(phi2) * sin((lambda1 - lambda2)); tmp = 0.0; if (phi1 <= -3e-128) tmp = atan2(t_1, (sin(phi2) - (cos(phi2) * t_0))); elseif (phi1 <= 5e-13) tmp = atan2((cos(phi2) * ((cos(lambda2) * sin(lambda1)) - (sin(lambda2) * cos(lambda1)))), sin(phi2)); else tmp = atan2(t_1, ((cos(phi1) * sin(phi2)) - t_0)); end tmp_2 = tmp; end
code[lambda1_, lambda2_, phi1_, phi2_] := Block[{t$95$0 = N[(N[Sin[phi1], $MachinePrecision] * N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[Cos[phi2], $MachinePrecision] * N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[phi1, -3e-128], N[ArcTan[t$95$1 / N[(N[Sin[phi2], $MachinePrecision] - N[(N[Cos[phi2], $MachinePrecision] * t$95$0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], If[LessEqual[phi1, 5e-13], N[ArcTan[N[(N[Cos[phi2], $MachinePrecision] * N[(N[(N[Cos[lambda2], $MachinePrecision] * N[Sin[lambda1], $MachinePrecision]), $MachinePrecision] - N[(N[Sin[lambda2], $MachinePrecision] * N[Cos[lambda1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[Sin[phi2], $MachinePrecision]], $MachinePrecision], N[ArcTan[t$95$1 / N[(N[(N[Cos[phi1], $MachinePrecision] * N[Sin[phi2], $MachinePrecision]), $MachinePrecision] - t$95$0), $MachinePrecision]], $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sin \phi_1 \cdot \cos \left(\lambda_1 - \lambda_2\right)\\
t_1 := \cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)\\
\mathbf{if}\;\phi_1 \leq -3 \cdot 10^{-128}:\\
\;\;\;\;\tan^{-1}_* \frac{t\_1}{\sin \phi_2 - \cos \phi_2 \cdot t\_0}\\
\mathbf{elif}\;\phi_1 \leq 5 \cdot 10^{-13}:\\
\;\;\;\;\tan^{-1}_* \frac{\cos \phi_2 \cdot \left(\cos \lambda_2 \cdot \sin \lambda_1 - \sin \lambda_2 \cdot \cos \lambda_1\right)}{\sin \phi_2}\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1}_* \frac{t\_1}{\cos \phi_1 \cdot \sin \phi_2 - t\_0}\\
\end{array}
\end{array}
if phi1 < -2.99999999999999978e-128Initial program 79.4%
*-commutative79.4%
associate-*l*79.4%
Simplified79.4%
Taylor expanded in phi1 around 0 59.3%
if -2.99999999999999978e-128 < phi1 < 4.9999999999999999e-13Initial program 82.0%
add-sqr-sqrt59.0%
sqrt-unprod82.0%
pow282.0%
Applied egg-rr82.0%
Taylor expanded in phi1 around 0 82.0%
sin-diff99.9%
Applied egg-rr99.9%
if 4.9999999999999999e-13 < phi1 Initial program 80.1%
add-sqr-sqrt61.1%
sqrt-unprod65.9%
pow265.9%
Applied egg-rr65.9%
Taylor expanded in phi2 around 0 54.2%
Final simplification73.4%
(FPCore (lambda1 lambda2 phi1 phi2)
:precision binary64
(if (or (<= phi2 -3.6e-7) (not (<= phi2 800000.0)))
(atan2
(*
(cos phi2)
(- (* (cos lambda2) (sin lambda1)) (* (sin lambda2) (cos lambda1))))
(sin phi2))
(atan2
(sin (- lambda1 lambda2))
(- (* (cos phi1) (sin phi2)) (* (sin phi1) (cos (- lambda1 lambda2)))))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
double tmp;
if ((phi2 <= -3.6e-7) || !(phi2 <= 800000.0)) {
tmp = atan2((cos(phi2) * ((cos(lambda2) * sin(lambda1)) - (sin(lambda2) * cos(lambda1)))), sin(phi2));
} else {
tmp = atan2(sin((lambda1 - lambda2)), ((cos(phi1) * sin(phi2)) - (sin(phi1) * cos((lambda1 - lambda2)))));
}
return tmp;
}
real(8) function code(lambda1, lambda2, phi1, phi2)
real(8), intent (in) :: lambda1
real(8), intent (in) :: lambda2
real(8), intent (in) :: phi1
real(8), intent (in) :: phi2
real(8) :: tmp
if ((phi2 <= (-3.6d-7)) .or. (.not. (phi2 <= 800000.0d0))) then
tmp = atan2((cos(phi2) * ((cos(lambda2) * sin(lambda1)) - (sin(lambda2) * cos(lambda1)))), sin(phi2))
else
tmp = atan2(sin((lambda1 - lambda2)), ((cos(phi1) * sin(phi2)) - (sin(phi1) * cos((lambda1 - lambda2)))))
end if
code = tmp
end function
public static double code(double lambda1, double lambda2, double phi1, double phi2) {
double tmp;
if ((phi2 <= -3.6e-7) || !(phi2 <= 800000.0)) {
tmp = Math.atan2((Math.cos(phi2) * ((Math.cos(lambda2) * Math.sin(lambda1)) - (Math.sin(lambda2) * Math.cos(lambda1)))), Math.sin(phi2));
} else {
tmp = Math.atan2(Math.sin((lambda1 - lambda2)), ((Math.cos(phi1) * Math.sin(phi2)) - (Math.sin(phi1) * Math.cos((lambda1 - lambda2)))));
}
return tmp;
}
def code(lambda1, lambda2, phi1, phi2): tmp = 0 if (phi2 <= -3.6e-7) or not (phi2 <= 800000.0): tmp = math.atan2((math.cos(phi2) * ((math.cos(lambda2) * math.sin(lambda1)) - (math.sin(lambda2) * math.cos(lambda1)))), math.sin(phi2)) else: tmp = math.atan2(math.sin((lambda1 - lambda2)), ((math.cos(phi1) * math.sin(phi2)) - (math.sin(phi1) * math.cos((lambda1 - lambda2))))) return tmp
function code(lambda1, lambda2, phi1, phi2) tmp = 0.0 if ((phi2 <= -3.6e-7) || !(phi2 <= 800000.0)) tmp = atan(Float64(cos(phi2) * Float64(Float64(cos(lambda2) * sin(lambda1)) - Float64(sin(lambda2) * cos(lambda1)))), sin(phi2)); else tmp = atan(sin(Float64(lambda1 - lambda2)), Float64(Float64(cos(phi1) * sin(phi2)) - Float64(sin(phi1) * cos(Float64(lambda1 - lambda2))))); end return tmp end
function tmp_2 = code(lambda1, lambda2, phi1, phi2) tmp = 0.0; if ((phi2 <= -3.6e-7) || ~((phi2 <= 800000.0))) tmp = atan2((cos(phi2) * ((cos(lambda2) * sin(lambda1)) - (sin(lambda2) * cos(lambda1)))), sin(phi2)); else tmp = atan2(sin((lambda1 - lambda2)), ((cos(phi1) * sin(phi2)) - (sin(phi1) * cos((lambda1 - lambda2))))); end tmp_2 = tmp; end
code[lambda1_, lambda2_, phi1_, phi2_] := If[Or[LessEqual[phi2, -3.6e-7], N[Not[LessEqual[phi2, 800000.0]], $MachinePrecision]], N[ArcTan[N[(N[Cos[phi2], $MachinePrecision] * N[(N[(N[Cos[lambda2], $MachinePrecision] * N[Sin[lambda1], $MachinePrecision]), $MachinePrecision] - N[(N[Sin[lambda2], $MachinePrecision] * N[Cos[lambda1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[Sin[phi2], $MachinePrecision]], $MachinePrecision], N[ArcTan[N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision] / N[(N[(N[Cos[phi1], $MachinePrecision] * N[Sin[phi2], $MachinePrecision]), $MachinePrecision] - N[(N[Sin[phi1], $MachinePrecision] * N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\phi_2 \leq -3.6 \cdot 10^{-7} \lor \neg \left(\phi_2 \leq 800000\right):\\
\;\;\;\;\tan^{-1}_* \frac{\cos \phi_2 \cdot \left(\cos \lambda_2 \cdot \sin \lambda_1 - \sin \lambda_2 \cdot \cos \lambda_1\right)}{\sin \phi_2}\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 \cdot \sin \phi_2 - \sin \phi_1 \cdot \cos \left(\lambda_1 - \lambda_2\right)}\\
\end{array}
\end{array}
if phi2 < -3.59999999999999994e-7 or 8e5 < phi2 Initial program 74.3%
add-sqr-sqrt53.0%
sqrt-unprod65.1%
pow265.1%
Applied egg-rr65.1%
Taylor expanded in phi1 around 0 45.9%
sin-diff85.8%
Applied egg-rr57.1%
if -3.59999999999999994e-7 < phi2 < 8e5Initial program 87.5%
Taylor expanded in phi2 around 0 87.0%
Taylor expanded in phi2 around 0 87.1%
Final simplification71.3%
(FPCore (lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (sin (- lambda1 lambda2))) (t_1 (* (cos phi2) t_0)))
(if (<= phi2 -3.8e+32)
(atan2 (* t_0 (expm1 (log1p (cos phi2)))) (sin phi2))
(if (<= phi2 800000.0)
(atan2
t_1
(- (* phi2 (cos phi1)) (* (sin phi1) (cos (- lambda1 lambda2)))))
(atan2 t_1 (sin phi2))))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = sin((lambda1 - lambda2));
double t_1 = cos(phi2) * t_0;
double tmp;
if (phi2 <= -3.8e+32) {
tmp = atan2((t_0 * expm1(log1p(cos(phi2)))), sin(phi2));
} else if (phi2 <= 800000.0) {
tmp = atan2(t_1, ((phi2 * cos(phi1)) - (sin(phi1) * cos((lambda1 - lambda2)))));
} else {
tmp = atan2(t_1, sin(phi2));
}
return tmp;
}
public static double code(double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = Math.sin((lambda1 - lambda2));
double t_1 = Math.cos(phi2) * t_0;
double tmp;
if (phi2 <= -3.8e+32) {
tmp = Math.atan2((t_0 * Math.expm1(Math.log1p(Math.cos(phi2)))), Math.sin(phi2));
} else if (phi2 <= 800000.0) {
tmp = Math.atan2(t_1, ((phi2 * Math.cos(phi1)) - (Math.sin(phi1) * Math.cos((lambda1 - lambda2)))));
} else {
tmp = Math.atan2(t_1, Math.sin(phi2));
}
return tmp;
}
def code(lambda1, lambda2, phi1, phi2): t_0 = math.sin((lambda1 - lambda2)) t_1 = math.cos(phi2) * t_0 tmp = 0 if phi2 <= -3.8e+32: tmp = math.atan2((t_0 * math.expm1(math.log1p(math.cos(phi2)))), math.sin(phi2)) elif phi2 <= 800000.0: tmp = math.atan2(t_1, ((phi2 * math.cos(phi1)) - (math.sin(phi1) * math.cos((lambda1 - lambda2))))) else: tmp = math.atan2(t_1, math.sin(phi2)) return tmp
function code(lambda1, lambda2, phi1, phi2) t_0 = sin(Float64(lambda1 - lambda2)) t_1 = Float64(cos(phi2) * t_0) tmp = 0.0 if (phi2 <= -3.8e+32) tmp = atan(Float64(t_0 * expm1(log1p(cos(phi2)))), sin(phi2)); elseif (phi2 <= 800000.0) tmp = atan(t_1, Float64(Float64(phi2 * cos(phi1)) - Float64(sin(phi1) * cos(Float64(lambda1 - lambda2))))); else tmp = atan(t_1, sin(phi2)); end return tmp end
code[lambda1_, lambda2_, phi1_, phi2_] := Block[{t$95$0 = N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$1 = N[(N[Cos[phi2], $MachinePrecision] * t$95$0), $MachinePrecision]}, If[LessEqual[phi2, -3.8e+32], N[ArcTan[N[(t$95$0 * N[(Exp[N[Log[1 + N[Cos[phi2], $MachinePrecision]], $MachinePrecision]] - 1), $MachinePrecision]), $MachinePrecision] / N[Sin[phi2], $MachinePrecision]], $MachinePrecision], If[LessEqual[phi2, 800000.0], N[ArcTan[t$95$1 / N[(N[(phi2 * N[Cos[phi1], $MachinePrecision]), $MachinePrecision] - N[(N[Sin[phi1], $MachinePrecision] * N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], N[ArcTan[t$95$1 / N[Sin[phi2], $MachinePrecision]], $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sin \left(\lambda_1 - \lambda_2\right)\\
t_1 := \cos \phi_2 \cdot t\_0\\
\mathbf{if}\;\phi_2 \leq -3.8 \cdot 10^{+32}:\\
\;\;\;\;\tan^{-1}_* \frac{t\_0 \cdot \mathsf{expm1}\left(\mathsf{log1p}\left(\cos \phi_2\right)\right)}{\sin \phi_2}\\
\mathbf{elif}\;\phi_2 \leq 800000:\\
\;\;\;\;\tan^{-1}_* \frac{t\_1}{\phi_2 \cdot \cos \phi_1 - \sin \phi_1 \cdot \cos \left(\lambda_1 - \lambda_2\right)}\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1}_* \frac{t\_1}{\sin \phi_2}\\
\end{array}
\end{array}
if phi2 < -3.8000000000000003e32Initial program 73.9%
add-sqr-sqrt53.1%
sqrt-unprod63.6%
pow263.6%
Applied egg-rr63.6%
Taylor expanded in phi1 around 0 46.4%
expm1-log1p-u46.4%
expm1-undefine46.3%
Applied egg-rr46.3%
expm1-define46.4%
Simplified46.4%
if -3.8000000000000003e32 < phi2 < 8e5Initial program 86.2%
add-sqr-sqrt60.4%
sqrt-unprod77.6%
pow277.6%
Applied egg-rr77.6%
Taylor expanded in phi2 around 0 83.9%
if 8e5 < phi2 Initial program 75.5%
add-sqr-sqrt51.6%
sqrt-unprod66.2%
pow266.2%
Applied egg-rr66.2%
Taylor expanded in phi1 around 0 46.7%
Final simplification65.4%
(FPCore (lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (sin (- lambda1 lambda2))))
(if (<= phi2 -6e-7)
(atan2 (* t_0 (expm1 (log1p (cos phi2)))) (sin phi2))
(if (<= phi2 920000.0)
(atan2
t_0
(- (* (cos phi1) (sin phi2)) (* (sin phi1) (cos (- lambda1 lambda2)))))
(atan2 (* (cos phi2) t_0) (sin phi2))))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = sin((lambda1 - lambda2));
double tmp;
if (phi2 <= -6e-7) {
tmp = atan2((t_0 * expm1(log1p(cos(phi2)))), sin(phi2));
} else if (phi2 <= 920000.0) {
tmp = atan2(t_0, ((cos(phi1) * sin(phi2)) - (sin(phi1) * cos((lambda1 - lambda2)))));
} else {
tmp = atan2((cos(phi2) * t_0), sin(phi2));
}
return tmp;
}
public static double code(double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = Math.sin((lambda1 - lambda2));
double tmp;
if (phi2 <= -6e-7) {
tmp = Math.atan2((t_0 * Math.expm1(Math.log1p(Math.cos(phi2)))), Math.sin(phi2));
} else if (phi2 <= 920000.0) {
tmp = Math.atan2(t_0, ((Math.cos(phi1) * Math.sin(phi2)) - (Math.sin(phi1) * Math.cos((lambda1 - lambda2)))));
} else {
tmp = Math.atan2((Math.cos(phi2) * t_0), Math.sin(phi2));
}
return tmp;
}
def code(lambda1, lambda2, phi1, phi2): t_0 = math.sin((lambda1 - lambda2)) tmp = 0 if phi2 <= -6e-7: tmp = math.atan2((t_0 * math.expm1(math.log1p(math.cos(phi2)))), math.sin(phi2)) elif phi2 <= 920000.0: tmp = math.atan2(t_0, ((math.cos(phi1) * math.sin(phi2)) - (math.sin(phi1) * math.cos((lambda1 - lambda2))))) else: tmp = math.atan2((math.cos(phi2) * t_0), math.sin(phi2)) return tmp
function code(lambda1, lambda2, phi1, phi2) t_0 = sin(Float64(lambda1 - lambda2)) tmp = 0.0 if (phi2 <= -6e-7) tmp = atan(Float64(t_0 * expm1(log1p(cos(phi2)))), sin(phi2)); elseif (phi2 <= 920000.0) tmp = atan(t_0, Float64(Float64(cos(phi1) * sin(phi2)) - Float64(sin(phi1) * cos(Float64(lambda1 - lambda2))))); else tmp = atan(Float64(cos(phi2) * t_0), sin(phi2)); end return tmp end
code[lambda1_, lambda2_, phi1_, phi2_] := Block[{t$95$0 = N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]}, If[LessEqual[phi2, -6e-7], N[ArcTan[N[(t$95$0 * N[(Exp[N[Log[1 + N[Cos[phi2], $MachinePrecision]], $MachinePrecision]] - 1), $MachinePrecision]), $MachinePrecision] / N[Sin[phi2], $MachinePrecision]], $MachinePrecision], If[LessEqual[phi2, 920000.0], N[ArcTan[t$95$0 / N[(N[(N[Cos[phi1], $MachinePrecision] * N[Sin[phi2], $MachinePrecision]), $MachinePrecision] - N[(N[Sin[phi1], $MachinePrecision] * N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], N[ArcTan[N[(N[Cos[phi2], $MachinePrecision] * t$95$0), $MachinePrecision] / N[Sin[phi2], $MachinePrecision]], $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sin \left(\lambda_1 - \lambda_2\right)\\
\mathbf{if}\;\phi_2 \leq -6 \cdot 10^{-7}:\\
\;\;\;\;\tan^{-1}_* \frac{t\_0 \cdot \mathsf{expm1}\left(\mathsf{log1p}\left(\cos \phi_2\right)\right)}{\sin \phi_2}\\
\mathbf{elif}\;\phi_2 \leq 920000:\\
\;\;\;\;\tan^{-1}_* \frac{t\_0}{\cos \phi_1 \cdot \sin \phi_2 - \sin \phi_1 \cdot \cos \left(\lambda_1 - \lambda_2\right)}\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1}_* \frac{\cos \phi_2 \cdot t\_0}{\sin \phi_2}\\
\end{array}
\end{array}
if phi2 < -5.9999999999999997e-7Initial program 73.1%
add-sqr-sqrt54.4%
sqrt-unprod64.0%
pow264.0%
Applied egg-rr64.0%
Taylor expanded in phi1 around 0 45.0%
expm1-log1p-u45.0%
expm1-undefine44.9%
Applied egg-rr44.9%
expm1-define45.0%
Simplified45.0%
if -5.9999999999999997e-7 < phi2 < 9.2e5Initial program 87.5%
Taylor expanded in phi2 around 0 87.0%
Taylor expanded in phi2 around 0 87.1%
if 9.2e5 < phi2 Initial program 75.5%
add-sqr-sqrt51.6%
sqrt-unprod66.2%
pow266.2%
Applied egg-rr66.2%
Taylor expanded in phi1 around 0 46.7%
Final simplification65.4%
(FPCore (lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (sin (- lambda1 lambda2))))
(if (<= phi2 -6e-7)
(atan2 (* t_0 (expm1 (log1p (cos phi2)))) (sin phi2))
(if (<= phi2 800000.0)
(atan2 t_0 (- (sin phi2) (* (sin phi1) (cos (- lambda1 lambda2)))))
(atan2 (* (cos phi2) t_0) (sin phi2))))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = sin((lambda1 - lambda2));
double tmp;
if (phi2 <= -6e-7) {
tmp = atan2((t_0 * expm1(log1p(cos(phi2)))), sin(phi2));
} else if (phi2 <= 800000.0) {
tmp = atan2(t_0, (sin(phi2) - (sin(phi1) * cos((lambda1 - lambda2)))));
} else {
tmp = atan2((cos(phi2) * t_0), sin(phi2));
}
return tmp;
}
public static double code(double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = Math.sin((lambda1 - lambda2));
double tmp;
if (phi2 <= -6e-7) {
tmp = Math.atan2((t_0 * Math.expm1(Math.log1p(Math.cos(phi2)))), Math.sin(phi2));
} else if (phi2 <= 800000.0) {
tmp = Math.atan2(t_0, (Math.sin(phi2) - (Math.sin(phi1) * Math.cos((lambda1 - lambda2)))));
} else {
tmp = Math.atan2((Math.cos(phi2) * t_0), Math.sin(phi2));
}
return tmp;
}
def code(lambda1, lambda2, phi1, phi2): t_0 = math.sin((lambda1 - lambda2)) tmp = 0 if phi2 <= -6e-7: tmp = math.atan2((t_0 * math.expm1(math.log1p(math.cos(phi2)))), math.sin(phi2)) elif phi2 <= 800000.0: tmp = math.atan2(t_0, (math.sin(phi2) - (math.sin(phi1) * math.cos((lambda1 - lambda2))))) else: tmp = math.atan2((math.cos(phi2) * t_0), math.sin(phi2)) return tmp
function code(lambda1, lambda2, phi1, phi2) t_0 = sin(Float64(lambda1 - lambda2)) tmp = 0.0 if (phi2 <= -6e-7) tmp = atan(Float64(t_0 * expm1(log1p(cos(phi2)))), sin(phi2)); elseif (phi2 <= 800000.0) tmp = atan(t_0, Float64(sin(phi2) - Float64(sin(phi1) * cos(Float64(lambda1 - lambda2))))); else tmp = atan(Float64(cos(phi2) * t_0), sin(phi2)); end return tmp end
code[lambda1_, lambda2_, phi1_, phi2_] := Block[{t$95$0 = N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]}, If[LessEqual[phi2, -6e-7], N[ArcTan[N[(t$95$0 * N[(Exp[N[Log[1 + N[Cos[phi2], $MachinePrecision]], $MachinePrecision]] - 1), $MachinePrecision]), $MachinePrecision] / N[Sin[phi2], $MachinePrecision]], $MachinePrecision], If[LessEqual[phi2, 800000.0], N[ArcTan[t$95$0 / N[(N[Sin[phi2], $MachinePrecision] - N[(N[Sin[phi1], $MachinePrecision] * N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], N[ArcTan[N[(N[Cos[phi2], $MachinePrecision] * t$95$0), $MachinePrecision] / N[Sin[phi2], $MachinePrecision]], $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sin \left(\lambda_1 - \lambda_2\right)\\
\mathbf{if}\;\phi_2 \leq -6 \cdot 10^{-7}:\\
\;\;\;\;\tan^{-1}_* \frac{t\_0 \cdot \mathsf{expm1}\left(\mathsf{log1p}\left(\cos \phi_2\right)\right)}{\sin \phi_2}\\
\mathbf{elif}\;\phi_2 \leq 800000:\\
\;\;\;\;\tan^{-1}_* \frac{t\_0}{\sin \phi_2 - \sin \phi_1 \cdot \cos \left(\lambda_1 - \lambda_2\right)}\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1}_* \frac{\cos \phi_2 \cdot t\_0}{\sin \phi_2}\\
\end{array}
\end{array}
if phi2 < -5.9999999999999997e-7Initial program 73.1%
add-sqr-sqrt54.4%
sqrt-unprod64.0%
pow264.0%
Applied egg-rr64.0%
Taylor expanded in phi1 around 0 45.0%
expm1-log1p-u45.0%
expm1-undefine44.9%
Applied egg-rr44.9%
expm1-define45.0%
Simplified45.0%
if -5.9999999999999997e-7 < phi2 < 8e5Initial program 87.5%
Taylor expanded in phi2 around 0 87.0%
Taylor expanded in phi1 around 0 85.9%
Taylor expanded in phi2 around 0 86.7%
if 8e5 < phi2 Initial program 75.5%
add-sqr-sqrt51.6%
sqrt-unprod66.2%
pow266.2%
Applied egg-rr66.2%
Taylor expanded in phi1 around 0 46.7%
Final simplification65.1%
(FPCore (lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (sin (- lambda1 lambda2))))
(if (or (<= phi2 -6e-7) (not (<= phi2 800000.0)))
(atan2 (* (cos phi2) t_0) (sin phi2))
(atan2 t_0 (- (sin phi2) (* (sin phi1) (cos (- lambda1 lambda2))))))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = sin((lambda1 - lambda2));
double tmp;
if ((phi2 <= -6e-7) || !(phi2 <= 800000.0)) {
tmp = atan2((cos(phi2) * t_0), sin(phi2));
} else {
tmp = atan2(t_0, (sin(phi2) - (sin(phi1) * cos((lambda1 - lambda2)))));
}
return tmp;
}
real(8) function code(lambda1, lambda2, phi1, phi2)
real(8), intent (in) :: lambda1
real(8), intent (in) :: lambda2
real(8), intent (in) :: phi1
real(8), intent (in) :: phi2
real(8) :: t_0
real(8) :: tmp
t_0 = sin((lambda1 - lambda2))
if ((phi2 <= (-6d-7)) .or. (.not. (phi2 <= 800000.0d0))) then
tmp = atan2((cos(phi2) * t_0), sin(phi2))
else
tmp = atan2(t_0, (sin(phi2) - (sin(phi1) * cos((lambda1 - lambda2)))))
end if
code = tmp
end function
public static double code(double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = Math.sin((lambda1 - lambda2));
double tmp;
if ((phi2 <= -6e-7) || !(phi2 <= 800000.0)) {
tmp = Math.atan2((Math.cos(phi2) * t_0), Math.sin(phi2));
} else {
tmp = Math.atan2(t_0, (Math.sin(phi2) - (Math.sin(phi1) * Math.cos((lambda1 - lambda2)))));
}
return tmp;
}
def code(lambda1, lambda2, phi1, phi2): t_0 = math.sin((lambda1 - lambda2)) tmp = 0 if (phi2 <= -6e-7) or not (phi2 <= 800000.0): tmp = math.atan2((math.cos(phi2) * t_0), math.sin(phi2)) else: tmp = math.atan2(t_0, (math.sin(phi2) - (math.sin(phi1) * math.cos((lambda1 - lambda2))))) return tmp
function code(lambda1, lambda2, phi1, phi2) t_0 = sin(Float64(lambda1 - lambda2)) tmp = 0.0 if ((phi2 <= -6e-7) || !(phi2 <= 800000.0)) tmp = atan(Float64(cos(phi2) * t_0), sin(phi2)); else tmp = atan(t_0, Float64(sin(phi2) - Float64(sin(phi1) * cos(Float64(lambda1 - lambda2))))); end return tmp end
function tmp_2 = code(lambda1, lambda2, phi1, phi2) t_0 = sin((lambda1 - lambda2)); tmp = 0.0; if ((phi2 <= -6e-7) || ~((phi2 <= 800000.0))) tmp = atan2((cos(phi2) * t_0), sin(phi2)); else tmp = atan2(t_0, (sin(phi2) - (sin(phi1) * cos((lambda1 - lambda2))))); end tmp_2 = tmp; end
code[lambda1_, lambda2_, phi1_, phi2_] := Block[{t$95$0 = N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]}, If[Or[LessEqual[phi2, -6e-7], N[Not[LessEqual[phi2, 800000.0]], $MachinePrecision]], N[ArcTan[N[(N[Cos[phi2], $MachinePrecision] * t$95$0), $MachinePrecision] / N[Sin[phi2], $MachinePrecision]], $MachinePrecision], N[ArcTan[t$95$0 / N[(N[Sin[phi2], $MachinePrecision] - N[(N[Sin[phi1], $MachinePrecision] * N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sin \left(\lambda_1 - \lambda_2\right)\\
\mathbf{if}\;\phi_2 \leq -6 \cdot 10^{-7} \lor \neg \left(\phi_2 \leq 800000\right):\\
\;\;\;\;\tan^{-1}_* \frac{\cos \phi_2 \cdot t\_0}{\sin \phi_2}\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1}_* \frac{t\_0}{\sin \phi_2 - \sin \phi_1 \cdot \cos \left(\lambda_1 - \lambda_2\right)}\\
\end{array}
\end{array}
if phi2 < -5.9999999999999997e-7 or 8e5 < phi2 Initial program 74.3%
add-sqr-sqrt53.0%
sqrt-unprod65.1%
pow265.1%
Applied egg-rr65.1%
Taylor expanded in phi1 around 0 45.9%
if -5.9999999999999997e-7 < phi2 < 8e5Initial program 87.5%
Taylor expanded in phi2 around 0 87.0%
Taylor expanded in phi1 around 0 85.9%
Taylor expanded in phi2 around 0 86.7%
Final simplification65.1%
(FPCore (lambda1 lambda2 phi1 phi2) :precision binary64 (if (or (<= lambda2 -4.3e-27) (not (<= lambda2 4e-88))) (atan2 (* (cos phi2) (sin (- lambda2))) (sin phi2)) (atan2 (* (sin lambda1) (cos phi2)) (sin phi2))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
double tmp;
if ((lambda2 <= -4.3e-27) || !(lambda2 <= 4e-88)) {
tmp = atan2((cos(phi2) * sin(-lambda2)), sin(phi2));
} else {
tmp = atan2((sin(lambda1) * cos(phi2)), sin(phi2));
}
return tmp;
}
real(8) function code(lambda1, lambda2, phi1, phi2)
real(8), intent (in) :: lambda1
real(8), intent (in) :: lambda2
real(8), intent (in) :: phi1
real(8), intent (in) :: phi2
real(8) :: tmp
if ((lambda2 <= (-4.3d-27)) .or. (.not. (lambda2 <= 4d-88))) then
tmp = atan2((cos(phi2) * sin(-lambda2)), sin(phi2))
else
tmp = atan2((sin(lambda1) * cos(phi2)), sin(phi2))
end if
code = tmp
end function
public static double code(double lambda1, double lambda2, double phi1, double phi2) {
double tmp;
if ((lambda2 <= -4.3e-27) || !(lambda2 <= 4e-88)) {
tmp = Math.atan2((Math.cos(phi2) * Math.sin(-lambda2)), Math.sin(phi2));
} else {
tmp = Math.atan2((Math.sin(lambda1) * Math.cos(phi2)), Math.sin(phi2));
}
return tmp;
}
def code(lambda1, lambda2, phi1, phi2): tmp = 0 if (lambda2 <= -4.3e-27) or not (lambda2 <= 4e-88): tmp = math.atan2((math.cos(phi2) * math.sin(-lambda2)), math.sin(phi2)) else: tmp = math.atan2((math.sin(lambda1) * math.cos(phi2)), math.sin(phi2)) return tmp
function code(lambda1, lambda2, phi1, phi2) tmp = 0.0 if ((lambda2 <= -4.3e-27) || !(lambda2 <= 4e-88)) tmp = atan(Float64(cos(phi2) * sin(Float64(-lambda2))), sin(phi2)); else tmp = atan(Float64(sin(lambda1) * cos(phi2)), sin(phi2)); end return tmp end
function tmp_2 = code(lambda1, lambda2, phi1, phi2) tmp = 0.0; if ((lambda2 <= -4.3e-27) || ~((lambda2 <= 4e-88))) tmp = atan2((cos(phi2) * sin(-lambda2)), sin(phi2)); else tmp = atan2((sin(lambda1) * cos(phi2)), sin(phi2)); end tmp_2 = tmp; end
code[lambda1_, lambda2_, phi1_, phi2_] := If[Or[LessEqual[lambda2, -4.3e-27], N[Not[LessEqual[lambda2, 4e-88]], $MachinePrecision]], N[ArcTan[N[(N[Cos[phi2], $MachinePrecision] * N[Sin[(-lambda2)], $MachinePrecision]), $MachinePrecision] / N[Sin[phi2], $MachinePrecision]], $MachinePrecision], N[ArcTan[N[(N[Sin[lambda1], $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision] / N[Sin[phi2], $MachinePrecision]], $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\lambda_2 \leq -4.3 \cdot 10^{-27} \lor \neg \left(\lambda_2 \leq 4 \cdot 10^{-88}\right):\\
\;\;\;\;\tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(-\lambda_2\right)}{\sin \phi_2}\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1}_* \frac{\sin \lambda_1 \cdot \cos \phi_2}{\sin \phi_2}\\
\end{array}
\end{array}
if lambda2 < -4.30000000000000002e-27 or 3.99999999999999974e-88 < lambda2 Initial program 67.3%
add-sqr-sqrt37.9%
sqrt-unprod55.3%
pow255.3%
Applied egg-rr55.3%
Taylor expanded in phi1 around 0 43.3%
Taylor expanded in lambda1 around 0 40.9%
if -4.30000000000000002e-27 < lambda2 < 3.99999999999999974e-88Initial program 99.8%
add-sqr-sqrt83.5%
sqrt-unprod94.7%
pow294.7%
Applied egg-rr94.7%
Taylor expanded in phi1 around 0 59.3%
Taylor expanded in lambda2 around 0 51.5%
Final simplification45.2%
(FPCore (lambda1 lambda2 phi1 phi2)
:precision binary64
(if (or (<= phi2 -1.1) (not (<= phi2 920000.0)))
(atan2 (* (sin lambda1) (cos phi2)) (sin phi2))
(atan2
(sin (- lambda1 lambda2))
(+ phi2 (* -0.16666666666666666 (pow phi2 3.0))))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
double tmp;
if ((phi2 <= -1.1) || !(phi2 <= 920000.0)) {
tmp = atan2((sin(lambda1) * cos(phi2)), sin(phi2));
} else {
tmp = atan2(sin((lambda1 - lambda2)), (phi2 + (-0.16666666666666666 * pow(phi2, 3.0))));
}
return tmp;
}
real(8) function code(lambda1, lambda2, phi1, phi2)
real(8), intent (in) :: lambda1
real(8), intent (in) :: lambda2
real(8), intent (in) :: phi1
real(8), intent (in) :: phi2
real(8) :: tmp
if ((phi2 <= (-1.1d0)) .or. (.not. (phi2 <= 920000.0d0))) then
tmp = atan2((sin(lambda1) * cos(phi2)), sin(phi2))
else
tmp = atan2(sin((lambda1 - lambda2)), (phi2 + ((-0.16666666666666666d0) * (phi2 ** 3.0d0))))
end if
code = tmp
end function
public static double code(double lambda1, double lambda2, double phi1, double phi2) {
double tmp;
if ((phi2 <= -1.1) || !(phi2 <= 920000.0)) {
tmp = Math.atan2((Math.sin(lambda1) * Math.cos(phi2)), Math.sin(phi2));
} else {
tmp = Math.atan2(Math.sin((lambda1 - lambda2)), (phi2 + (-0.16666666666666666 * Math.pow(phi2, 3.0))));
}
return tmp;
}
def code(lambda1, lambda2, phi1, phi2): tmp = 0 if (phi2 <= -1.1) or not (phi2 <= 920000.0): tmp = math.atan2((math.sin(lambda1) * math.cos(phi2)), math.sin(phi2)) else: tmp = math.atan2(math.sin((lambda1 - lambda2)), (phi2 + (-0.16666666666666666 * math.pow(phi2, 3.0)))) return tmp
function code(lambda1, lambda2, phi1, phi2) tmp = 0.0 if ((phi2 <= -1.1) || !(phi2 <= 920000.0)) tmp = atan(Float64(sin(lambda1) * cos(phi2)), sin(phi2)); else tmp = atan(sin(Float64(lambda1 - lambda2)), Float64(phi2 + Float64(-0.16666666666666666 * (phi2 ^ 3.0)))); end return tmp end
function tmp_2 = code(lambda1, lambda2, phi1, phi2) tmp = 0.0; if ((phi2 <= -1.1) || ~((phi2 <= 920000.0))) tmp = atan2((sin(lambda1) * cos(phi2)), sin(phi2)); else tmp = atan2(sin((lambda1 - lambda2)), (phi2 + (-0.16666666666666666 * (phi2 ^ 3.0)))); end tmp_2 = tmp; end
code[lambda1_, lambda2_, phi1_, phi2_] := If[Or[LessEqual[phi2, -1.1], N[Not[LessEqual[phi2, 920000.0]], $MachinePrecision]], N[ArcTan[N[(N[Sin[lambda1], $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision] / N[Sin[phi2], $MachinePrecision]], $MachinePrecision], N[ArcTan[N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision] / N[(phi2 + N[(-0.16666666666666666 * N[Power[phi2, 3.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\phi_2 \leq -1.1 \lor \neg \left(\phi_2 \leq 920000\right):\\
\;\;\;\;\tan^{-1}_* \frac{\sin \lambda_1 \cdot \cos \phi_2}{\sin \phi_2}\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right)}{\phi_2 + -0.16666666666666666 \cdot {\phi_2}^{3}}\\
\end{array}
\end{array}
if phi2 < -1.1000000000000001 or 9.2e5 < phi2 Initial program 74.5%
add-sqr-sqrt52.9%
sqrt-unprod65.2%
pow265.2%
Applied egg-rr65.2%
Taylor expanded in phi1 around 0 45.6%
Taylor expanded in lambda2 around 0 27.4%
if -1.1000000000000001 < phi2 < 9.2e5Initial program 87.1%
add-sqr-sqrt60.1%
sqrt-unprod78.0%
pow278.0%
Applied egg-rr78.0%
Taylor expanded in phi1 around 0 54.3%
Taylor expanded in phi2 around 0 53.7%
Taylor expanded in phi2 around 0 54.5%
distribute-rgt-in54.5%
*-lft-identity54.5%
associate-*l*54.5%
pow-plus54.5%
metadata-eval54.5%
Simplified54.5%
Final simplification40.4%
(FPCore (lambda1 lambda2 phi1 phi2) :precision binary64 (atan2 (* (cos phi2) (sin (- lambda1 lambda2))) (sin phi2)))
double code(double lambda1, double lambda2, double phi1, double phi2) {
return atan2((cos(phi2) * sin((lambda1 - lambda2))), sin(phi2));
}
real(8) function code(lambda1, lambda2, phi1, phi2)
real(8), intent (in) :: lambda1
real(8), intent (in) :: lambda2
real(8), intent (in) :: phi1
real(8), intent (in) :: phi2
code = atan2((cos(phi2) * sin((lambda1 - lambda2))), sin(phi2))
end function
public static double code(double lambda1, double lambda2, double phi1, double phi2) {
return Math.atan2((Math.cos(phi2) * Math.sin((lambda1 - lambda2))), Math.sin(phi2));
}
def code(lambda1, lambda2, phi1, phi2): return math.atan2((math.cos(phi2) * math.sin((lambda1 - lambda2))), math.sin(phi2))
function code(lambda1, lambda2, phi1, phi2) return atan(Float64(cos(phi2) * sin(Float64(lambda1 - lambda2))), sin(phi2)) end
function tmp = code(lambda1, lambda2, phi1, phi2) tmp = atan2((cos(phi2) * sin((lambda1 - lambda2))), sin(phi2)); end
code[lambda1_, lambda2_, phi1_, phi2_] := N[ArcTan[N[(N[Cos[phi2], $MachinePrecision] * N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / N[Sin[phi2], $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
\\
\tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\sin \phi_2}
\end{array}
Initial program 80.5%
add-sqr-sqrt56.4%
sqrt-unprod71.3%
pow271.3%
Applied egg-rr71.3%
Taylor expanded in phi1 around 0 49.8%
Final simplification49.8%
(FPCore (lambda1 lambda2 phi1 phi2) :precision binary64 (if (<= phi2 2.9e+72) (atan2 (sin (- lambda1 lambda2)) phi2) (atan2 (sin (- lambda2)) (sin phi2))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
double tmp;
if (phi2 <= 2.9e+72) {
tmp = atan2(sin((lambda1 - lambda2)), phi2);
} else {
tmp = atan2(sin(-lambda2), sin(phi2));
}
return tmp;
}
real(8) function code(lambda1, lambda2, phi1, phi2)
real(8), intent (in) :: lambda1
real(8), intent (in) :: lambda2
real(8), intent (in) :: phi1
real(8), intent (in) :: phi2
real(8) :: tmp
if (phi2 <= 2.9d+72) then
tmp = atan2(sin((lambda1 - lambda2)), phi2)
else
tmp = atan2(sin(-lambda2), sin(phi2))
end if
code = tmp
end function
public static double code(double lambda1, double lambda2, double phi1, double phi2) {
double tmp;
if (phi2 <= 2.9e+72) {
tmp = Math.atan2(Math.sin((lambda1 - lambda2)), phi2);
} else {
tmp = Math.atan2(Math.sin(-lambda2), Math.sin(phi2));
}
return tmp;
}
def code(lambda1, lambda2, phi1, phi2): tmp = 0 if phi2 <= 2.9e+72: tmp = math.atan2(math.sin((lambda1 - lambda2)), phi2) else: tmp = math.atan2(math.sin(-lambda2), math.sin(phi2)) return tmp
function code(lambda1, lambda2, phi1, phi2) tmp = 0.0 if (phi2 <= 2.9e+72) tmp = atan(sin(Float64(lambda1 - lambda2)), phi2); else tmp = atan(sin(Float64(-lambda2)), sin(phi2)); end return tmp end
function tmp_2 = code(lambda1, lambda2, phi1, phi2) tmp = 0.0; if (phi2 <= 2.9e+72) tmp = atan2(sin((lambda1 - lambda2)), phi2); else tmp = atan2(sin(-lambda2), sin(phi2)); end tmp_2 = tmp; end
code[lambda1_, lambda2_, phi1_, phi2_] := If[LessEqual[phi2, 2.9e+72], N[ArcTan[N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision] / phi2], $MachinePrecision], N[ArcTan[N[Sin[(-lambda2)], $MachinePrecision] / N[Sin[phi2], $MachinePrecision]], $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\phi_2 \leq 2.9 \cdot 10^{+72}:\\
\;\;\;\;\tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right)}{\phi_2}\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1}_* \frac{\sin \left(-\lambda_2\right)}{\sin \phi_2}\\
\end{array}
\end{array}
if phi2 < 2.90000000000000017e72Initial program 83.2%
add-sqr-sqrt59.5%
sqrt-unprod74.3%
pow274.3%
Applied egg-rr74.3%
Taylor expanded in phi1 around 0 52.0%
Taylor expanded in phi2 around 0 38.3%
Taylor expanded in phi2 around 0 38.3%
if 2.90000000000000017e72 < phi2 Initial program 70.6%
add-sqr-sqrt44.7%
sqrt-unprod60.3%
pow260.3%
Applied egg-rr60.3%
Taylor expanded in phi1 around 0 41.6%
Taylor expanded in phi2 around 0 16.2%
Taylor expanded in lambda1 around 0 16.6%
(FPCore (lambda1 lambda2 phi1 phi2) :precision binary64 (if (<= phi2 3000000.0) (atan2 (sin (- lambda1 lambda2)) phi2) (atan2 (sin lambda1) (sin phi2))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
double tmp;
if (phi2 <= 3000000.0) {
tmp = atan2(sin((lambda1 - lambda2)), phi2);
} else {
tmp = atan2(sin(lambda1), sin(phi2));
}
return tmp;
}
real(8) function code(lambda1, lambda2, phi1, phi2)
real(8), intent (in) :: lambda1
real(8), intent (in) :: lambda2
real(8), intent (in) :: phi1
real(8), intent (in) :: phi2
real(8) :: tmp
if (phi2 <= 3000000.0d0) then
tmp = atan2(sin((lambda1 - lambda2)), phi2)
else
tmp = atan2(sin(lambda1), sin(phi2))
end if
code = tmp
end function
public static double code(double lambda1, double lambda2, double phi1, double phi2) {
double tmp;
if (phi2 <= 3000000.0) {
tmp = Math.atan2(Math.sin((lambda1 - lambda2)), phi2);
} else {
tmp = Math.atan2(Math.sin(lambda1), Math.sin(phi2));
}
return tmp;
}
def code(lambda1, lambda2, phi1, phi2): tmp = 0 if phi2 <= 3000000.0: tmp = math.atan2(math.sin((lambda1 - lambda2)), phi2) else: tmp = math.atan2(math.sin(lambda1), math.sin(phi2)) return tmp
function code(lambda1, lambda2, phi1, phi2) tmp = 0.0 if (phi2 <= 3000000.0) tmp = atan(sin(Float64(lambda1 - lambda2)), phi2); else tmp = atan(sin(lambda1), sin(phi2)); end return tmp end
function tmp_2 = code(lambda1, lambda2, phi1, phi2) tmp = 0.0; if (phi2 <= 3000000.0) tmp = atan2(sin((lambda1 - lambda2)), phi2); else tmp = atan2(sin(lambda1), sin(phi2)); end tmp_2 = tmp; end
code[lambda1_, lambda2_, phi1_, phi2_] := If[LessEqual[phi2, 3000000.0], N[ArcTan[N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision] / phi2], $MachinePrecision], N[ArcTan[N[Sin[lambda1], $MachinePrecision] / N[Sin[phi2], $MachinePrecision]], $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\phi_2 \leq 3000000:\\
\;\;\;\;\tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right)}{\phi_2}\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1}_* \frac{\sin \lambda_1}{\sin \phi_2}\\
\end{array}
\end{array}
if phi2 < 3e6Initial program 82.5%
add-sqr-sqrt58.3%
sqrt-unprod73.3%
pow273.3%
Applied egg-rr73.3%
Taylor expanded in phi1 around 0 51.2%
Taylor expanded in phi2 around 0 40.6%
Taylor expanded in phi2 around 0 40.8%
if 3e6 < phi2 Initial program 75.2%
add-sqr-sqrt50.9%
sqrt-unprod65.6%
pow265.6%
Applied egg-rr65.6%
Taylor expanded in phi1 around 0 45.9%
Taylor expanded in phi2 around 0 14.0%
Taylor expanded in lambda2 around 0 12.5%
(FPCore (lambda1 lambda2 phi1 phi2) :precision binary64 (atan2 (sin (- lambda1 lambda2)) (sin phi2)))
double code(double lambda1, double lambda2, double phi1, double phi2) {
return atan2(sin((lambda1 - lambda2)), sin(phi2));
}
real(8) function code(lambda1, lambda2, phi1, phi2)
real(8), intent (in) :: lambda1
real(8), intent (in) :: lambda2
real(8), intent (in) :: phi1
real(8), intent (in) :: phi2
code = atan2(sin((lambda1 - lambda2)), sin(phi2))
end function
public static double code(double lambda1, double lambda2, double phi1, double phi2) {
return Math.atan2(Math.sin((lambda1 - lambda2)), Math.sin(phi2));
}
def code(lambda1, lambda2, phi1, phi2): return math.atan2(math.sin((lambda1 - lambda2)), math.sin(phi2))
function code(lambda1, lambda2, phi1, phi2) return atan(sin(Float64(lambda1 - lambda2)), sin(phi2)) end
function tmp = code(lambda1, lambda2, phi1, phi2) tmp = atan2(sin((lambda1 - lambda2)), sin(phi2)); end
code[lambda1_, lambda2_, phi1_, phi2_] := N[ArcTan[N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision] / N[Sin[phi2], $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
\\
\tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right)}{\sin \phi_2}
\end{array}
Initial program 80.5%
add-sqr-sqrt56.4%
sqrt-unprod71.3%
pow271.3%
Applied egg-rr71.3%
Taylor expanded in phi1 around 0 49.8%
Taylor expanded in phi2 around 0 33.7%
(FPCore (lambda1 lambda2 phi1 phi2) :precision binary64 (atan2 (sin (- lambda1 lambda2)) phi2))
double code(double lambda1, double lambda2, double phi1, double phi2) {
return atan2(sin((lambda1 - lambda2)), phi2);
}
real(8) function code(lambda1, lambda2, phi1, phi2)
real(8), intent (in) :: lambda1
real(8), intent (in) :: lambda2
real(8), intent (in) :: phi1
real(8), intent (in) :: phi2
code = atan2(sin((lambda1 - lambda2)), phi2)
end function
public static double code(double lambda1, double lambda2, double phi1, double phi2) {
return Math.atan2(Math.sin((lambda1 - lambda2)), phi2);
}
def code(lambda1, lambda2, phi1, phi2): return math.atan2(math.sin((lambda1 - lambda2)), phi2)
function code(lambda1, lambda2, phi1, phi2) return atan(sin(Float64(lambda1 - lambda2)), phi2) end
function tmp = code(lambda1, lambda2, phi1, phi2) tmp = atan2(sin((lambda1 - lambda2)), phi2); end
code[lambda1_, lambda2_, phi1_, phi2_] := N[ArcTan[N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision] / phi2], $MachinePrecision]
\begin{array}{l}
\\
\tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right)}{\phi_2}
\end{array}
Initial program 80.5%
add-sqr-sqrt56.4%
sqrt-unprod71.3%
pow271.3%
Applied egg-rr71.3%
Taylor expanded in phi1 around 0 49.8%
Taylor expanded in phi2 around 0 33.7%
Taylor expanded in phi2 around 0 31.0%
herbie shell --seed 2024143
(FPCore (lambda1 lambda2 phi1 phi2)
:name "Bearing on a great circle"
:precision binary64
(atan2 (* (sin (- lambda1 lambda2)) (cos phi2)) (- (* (cos phi1) (sin phi2)) (* (* (sin phi1) (cos phi2)) (cos (- lambda1 lambda2))))))