
(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 26 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 lambda1) (cos lambda2) (* (cos lambda1) (- (sin lambda2))))
(cos phi2))
(-
(* (cos phi1) (sin phi2))
(*
(* (cos phi2) (sin phi1))
(fma (cos lambda2) (cos lambda1) (* (sin lambda1) (sin lambda2)))))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
return atan2((fma(sin(lambda1), cos(lambda2), (cos(lambda1) * -sin(lambda2))) * cos(phi2)), ((cos(phi1) * sin(phi2)) - ((cos(phi2) * sin(phi1)) * fma(cos(lambda2), cos(lambda1), (sin(lambda1) * sin(lambda2))))));
}
function code(lambda1, lambda2, phi1, phi2) return atan(Float64(fma(sin(lambda1), cos(lambda2), Float64(cos(lambda1) * Float64(-sin(lambda2)))) * cos(phi2)), Float64(Float64(cos(phi1) * sin(phi2)) - Float64(Float64(cos(phi2) * sin(phi1)) * fma(cos(lambda2), cos(lambda1), Float64(sin(lambda1) * sin(lambda2)))))) end
code[lambda1_, lambda2_, phi1_, phi2_] := N[ArcTan[N[(N[(N[Sin[lambda1], $MachinePrecision] * N[Cos[lambda2], $MachinePrecision] + N[(N[Cos[lambda1], $MachinePrecision] * (-N[Sin[lambda2], $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[lambda1], $MachinePrecision] * N[Sin[lambda2], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
\\
\tan^{-1}_* \frac{\mathsf{fma}\left(\sin \lambda_1, \cos \lambda_2, \cos \lambda_1 \cdot \left(-\sin \lambda_2\right)\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_1 \cdot \sin \lambda_2\right)}
\end{array}
Initial program 81.3%
sin-diff89.7%
fma-neg89.7%
Applied egg-rr89.7%
cos-diff99.7%
*-commutative99.7%
Applied egg-rr99.7%
fma-define99.7%
Simplified99.7%
Final simplification99.7%
(FPCore (lambda1 lambda2 phi1 phi2)
:precision binary64
(atan2
(*
(fma (sin lambda1) (cos lambda2) (* (cos lambda1) (- (sin lambda2))))
(cos phi2))
(-
(* (cos phi1) (sin phi2))
(*
(* (cos phi2) (sin phi1))
(+ (* (sin lambda1) (sin lambda2)) (* (cos lambda2) (cos lambda1)))))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
return atan2((fma(sin(lambda1), cos(lambda2), (cos(lambda1) * -sin(lambda2))) * cos(phi2)), ((cos(phi1) * sin(phi2)) - ((cos(phi2) * sin(phi1)) * ((sin(lambda1) * sin(lambda2)) + (cos(lambda2) * cos(lambda1))))));
}
function code(lambda1, lambda2, phi1, phi2) return atan(Float64(fma(sin(lambda1), cos(lambda2), Float64(cos(lambda1) * Float64(-sin(lambda2)))) * cos(phi2)), Float64(Float64(cos(phi1) * sin(phi2)) - Float64(Float64(cos(phi2) * sin(phi1)) * Float64(Float64(sin(lambda1) * sin(lambda2)) + Float64(cos(lambda2) * cos(lambda1)))))) end
code[lambda1_, lambda2_, phi1_, phi2_] := N[ArcTan[N[(N[(N[Sin[lambda1], $MachinePrecision] * N[Cos[lambda2], $MachinePrecision] + N[(N[Cos[lambda1], $MachinePrecision] * (-N[Sin[lambda2], $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[lambda1], $MachinePrecision] * N[Sin[lambda2], $MachinePrecision]), $MachinePrecision] + N[(N[Cos[lambda2], $MachinePrecision] * N[Cos[lambda1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
\\
\tan^{-1}_* \frac{\mathsf{fma}\left(\sin \lambda_1, \cos \lambda_2, \cos \lambda_1 \cdot \left(-\sin \lambda_2\right)\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_1 \cdot \sin \lambda_2 + \cos \lambda_2 \cdot \cos \lambda_1\right)}
\end{array}
Initial program 81.3%
sin-diff89.7%
fma-neg89.7%
Applied egg-rr89.7%
cos-diff81.7%
+-commutative81.7%
*-commutative81.7%
Applied egg-rr99.7%
Final simplification99.7%
(FPCore (lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (* (cos phi1) (sin phi2)))
(t_1 (expm1 t_0))
(t_2
(*
(fma (sin lambda1) (cos lambda2) (* (cos lambda1) (- (sin lambda2))))
(cos phi2)))
(t_3 (* (cos phi2) (sin phi1)))
(t_4 (* t_3 (cos (- lambda1 lambda2)))))
(if (<= phi2 -1e-5)
(atan2 t_2 (- (log1p t_1) t_4))
(if (<= phi2 4.2e-7)
(atan2
(- (* (sin lambda1) (cos lambda2)) (* (cos lambda1) (sin lambda2)))
(-
t_0
(*
t_3
(fma (cos lambda2) (cos lambda1) (* (sin lambda1) (sin lambda2))))))
(atan2 t_2 (- (log (+ t_1 1.0)) t_4))))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = cos(phi1) * sin(phi2);
double t_1 = expm1(t_0);
double t_2 = fma(sin(lambda1), cos(lambda2), (cos(lambda1) * -sin(lambda2))) * cos(phi2);
double t_3 = cos(phi2) * sin(phi1);
double t_4 = t_3 * cos((lambda1 - lambda2));
double tmp;
if (phi2 <= -1e-5) {
tmp = atan2(t_2, (log1p(t_1) - t_4));
} else if (phi2 <= 4.2e-7) {
tmp = atan2(((sin(lambda1) * cos(lambda2)) - (cos(lambda1) * sin(lambda2))), (t_0 - (t_3 * fma(cos(lambda2), cos(lambda1), (sin(lambda1) * sin(lambda2))))));
} else {
tmp = atan2(t_2, (log((t_1 + 1.0)) - t_4));
}
return tmp;
}
function code(lambda1, lambda2, phi1, phi2) t_0 = Float64(cos(phi1) * sin(phi2)) t_1 = expm1(t_0) t_2 = Float64(fma(sin(lambda1), cos(lambda2), Float64(cos(lambda1) * Float64(-sin(lambda2)))) * cos(phi2)) t_3 = Float64(cos(phi2) * sin(phi1)) t_4 = Float64(t_3 * cos(Float64(lambda1 - lambda2))) tmp = 0.0 if (phi2 <= -1e-5) tmp = atan(t_2, Float64(log1p(t_1) - t_4)); elseif (phi2 <= 4.2e-7) tmp = atan(Float64(Float64(sin(lambda1) * cos(lambda2)) - Float64(cos(lambda1) * sin(lambda2))), Float64(t_0 - Float64(t_3 * fma(cos(lambda2), cos(lambda1), Float64(sin(lambda1) * sin(lambda2)))))); else tmp = atan(t_2, Float64(log(Float64(t_1 + 1.0)) - t_4)); 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[(Exp[t$95$0] - 1), $MachinePrecision]}, Block[{t$95$2 = N[(N[(N[Sin[lambda1], $MachinePrecision] * N[Cos[lambda2], $MachinePrecision] + N[(N[Cos[lambda1], $MachinePrecision] * (-N[Sin[lambda2], $MachinePrecision])), $MachinePrecision]), $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = N[(N[Cos[phi2], $MachinePrecision] * N[Sin[phi1], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$4 = N[(t$95$3 * N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[phi2, -1e-5], N[ArcTan[t$95$2 / N[(N[Log[1 + t$95$1], $MachinePrecision] - t$95$4), $MachinePrecision]], $MachinePrecision], If[LessEqual[phi2, 4.2e-7], N[ArcTan[N[(N[(N[Sin[lambda1], $MachinePrecision] * N[Cos[lambda2], $MachinePrecision]), $MachinePrecision] - N[(N[Cos[lambda1], $MachinePrecision] * N[Sin[lambda2], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(t$95$0 - N[(t$95$3 * N[(N[Cos[lambda2], $MachinePrecision] * N[Cos[lambda1], $MachinePrecision] + N[(N[Sin[lambda1], $MachinePrecision] * N[Sin[lambda2], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], N[ArcTan[t$95$2 / N[(N[Log[N[(t$95$1 + 1.0), $MachinePrecision]], $MachinePrecision] - t$95$4), $MachinePrecision]], $MachinePrecision]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \cos \phi_1 \cdot \sin \phi_2\\
t_1 := \mathsf{expm1}\left(t\_0\right)\\
t_2 := \mathsf{fma}\left(\sin \lambda_1, \cos \lambda_2, \cos \lambda_1 \cdot \left(-\sin \lambda_2\right)\right) \cdot \cos \phi_2\\
t_3 := \cos \phi_2 \cdot \sin \phi_1\\
t_4 := t\_3 \cdot \cos \left(\lambda_1 - \lambda_2\right)\\
\mathbf{if}\;\phi_2 \leq -1 \cdot 10^{-5}:\\
\;\;\;\;\tan^{-1}_* \frac{t\_2}{\mathsf{log1p}\left(t\_1\right) - t\_4}\\
\mathbf{elif}\;\phi_2 \leq 4.2 \cdot 10^{-7}:\\
\;\;\;\;\tan^{-1}_* \frac{\sin \lambda_1 \cdot \cos \lambda_2 - \cos \lambda_1 \cdot \sin \lambda_2}{t\_0 - t\_3 \cdot \mathsf{fma}\left(\cos \lambda_2, \cos \lambda_1, \sin \lambda_1 \cdot \sin \lambda_2\right)}\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1}_* \frac{t\_2}{\log \left(t\_1 + 1\right) - t\_4}\\
\end{array}
\end{array}
if phi2 < -1.00000000000000008e-5Initial program 80.2%
sin-diff93.4%
fma-neg93.4%
Applied egg-rr93.4%
log1p-expm1-u93.4%
Applied egg-rr93.4%
if -1.00000000000000008e-5 < phi2 < 4.2e-7Initial program 83.8%
sin-diff88.8%
fma-neg88.8%
Applied egg-rr88.8%
cos-diff99.8%
*-commutative99.8%
Applied egg-rr99.8%
fma-define99.8%
Simplified99.8%
Taylor expanded in phi2 around 0 99.8%
if 4.2e-7 < phi2 Initial program 77.4%
sin-diff88.2%
fma-neg88.2%
Applied egg-rr88.2%
log1p-expm1-u88.2%
log1p-undefine88.4%
Applied egg-rr88.4%
Final simplification95.3%
(FPCore (lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (* (cos phi1) (sin phi2)))
(t_1 (expm1 t_0))
(t_2
(*
(fma (sin lambda1) (cos lambda2) (* (cos lambda1) (- (sin lambda2))))
(cos phi2)))
(t_3 (* (* (cos phi2) (sin phi1)) (cos (- lambda1 lambda2)))))
(if (<= phi2 -6.2e-5)
(atan2 t_2 (- (log1p t_1) t_3))
(if (<= phi2 0.00065)
(atan2
t_2
(-
t_0
(*
(sin phi1)
(+
(* (sin lambda1) (sin lambda2))
(* (cos lambda2) (cos lambda1))))))
(atan2 t_2 (- (log (+ t_1 1.0)) t_3))))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = cos(phi1) * sin(phi2);
double t_1 = expm1(t_0);
double t_2 = fma(sin(lambda1), cos(lambda2), (cos(lambda1) * -sin(lambda2))) * cos(phi2);
double t_3 = (cos(phi2) * sin(phi1)) * cos((lambda1 - lambda2));
double tmp;
if (phi2 <= -6.2e-5) {
tmp = atan2(t_2, (log1p(t_1) - t_3));
} else if (phi2 <= 0.00065) {
tmp = atan2(t_2, (t_0 - (sin(phi1) * ((sin(lambda1) * sin(lambda2)) + (cos(lambda2) * cos(lambda1))))));
} else {
tmp = atan2(t_2, (log((t_1 + 1.0)) - t_3));
}
return tmp;
}
function code(lambda1, lambda2, phi1, phi2) t_0 = Float64(cos(phi1) * sin(phi2)) t_1 = expm1(t_0) t_2 = Float64(fma(sin(lambda1), cos(lambda2), Float64(cos(lambda1) * Float64(-sin(lambda2)))) * cos(phi2)) t_3 = Float64(Float64(cos(phi2) * sin(phi1)) * cos(Float64(lambda1 - lambda2))) tmp = 0.0 if (phi2 <= -6.2e-5) tmp = atan(t_2, Float64(log1p(t_1) - t_3)); elseif (phi2 <= 0.00065) tmp = atan(t_2, Float64(t_0 - Float64(sin(phi1) * Float64(Float64(sin(lambda1) * sin(lambda2)) + Float64(cos(lambda2) * cos(lambda1)))))); else tmp = atan(t_2, Float64(log(Float64(t_1 + 1.0)) - t_3)); 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[(Exp[t$95$0] - 1), $MachinePrecision]}, Block[{t$95$2 = N[(N[(N[Sin[lambda1], $MachinePrecision] * N[Cos[lambda2], $MachinePrecision] + N[(N[Cos[lambda1], $MachinePrecision] * (-N[Sin[lambda2], $MachinePrecision])), $MachinePrecision]), $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = N[(N[(N[Cos[phi2], $MachinePrecision] * N[Sin[phi1], $MachinePrecision]), $MachinePrecision] * N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[phi2, -6.2e-5], N[ArcTan[t$95$2 / N[(N[Log[1 + t$95$1], $MachinePrecision] - t$95$3), $MachinePrecision]], $MachinePrecision], If[LessEqual[phi2, 0.00065], N[ArcTan[t$95$2 / N[(t$95$0 - N[(N[Sin[phi1], $MachinePrecision] * N[(N[(N[Sin[lambda1], $MachinePrecision] * N[Sin[lambda2], $MachinePrecision]), $MachinePrecision] + N[(N[Cos[lambda2], $MachinePrecision] * N[Cos[lambda1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], N[ArcTan[t$95$2 / N[(N[Log[N[(t$95$1 + 1.0), $MachinePrecision]], $MachinePrecision] - t$95$3), $MachinePrecision]], $MachinePrecision]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \cos \phi_1 \cdot \sin \phi_2\\
t_1 := \mathsf{expm1}\left(t\_0\right)\\
t_2 := \mathsf{fma}\left(\sin \lambda_1, \cos \lambda_2, \cos \lambda_1 \cdot \left(-\sin \lambda_2\right)\right) \cdot \cos \phi_2\\
t_3 := \left(\cos \phi_2 \cdot \sin \phi_1\right) \cdot \cos \left(\lambda_1 - \lambda_2\right)\\
\mathbf{if}\;\phi_2 \leq -6.2 \cdot 10^{-5}:\\
\;\;\;\;\tan^{-1}_* \frac{t\_2}{\mathsf{log1p}\left(t\_1\right) - t\_3}\\
\mathbf{elif}\;\phi_2 \leq 0.00065:\\
\;\;\;\;\tan^{-1}_* \frac{t\_2}{t\_0 - \sin \phi_1 \cdot \left(\sin \lambda_1 \cdot \sin \lambda_2 + \cos \lambda_2 \cdot \cos \lambda_1\right)}\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1}_* \frac{t\_2}{\log \left(t\_1 + 1\right) - t\_3}\\
\end{array}
\end{array}
if phi2 < -6.20000000000000027e-5Initial program 80.2%
sin-diff93.4%
fma-neg93.4%
Applied egg-rr93.4%
log1p-expm1-u93.4%
Applied egg-rr93.4%
if -6.20000000000000027e-5 < phi2 < 6.4999999999999997e-4Initial program 83.8%
sin-diff88.8%
fma-neg88.8%
Applied egg-rr88.8%
cos-diff99.8%
*-commutative99.8%
Applied egg-rr99.8%
fma-define99.8%
Simplified99.8%
Taylor expanded in phi2 around 0 99.8%
if 6.4999999999999997e-4 < phi2 Initial program 77.4%
sin-diff88.2%
fma-neg88.2%
Applied egg-rr88.2%
log1p-expm1-u88.2%
log1p-undefine88.4%
Applied egg-rr88.4%
Final simplification95.3%
(FPCore (lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (* (cos phi1) (sin phi2)))
(t_1
(*
(fma (sin lambda1) (cos lambda2) (* (cos lambda1) (- (sin lambda2))))
(cos phi2)))
(t_2 (* (cos phi2) (sin phi1)))
(t_3 (* t_2 (cos (- lambda1 lambda2)))))
(if (<= phi2 -7e-6)
(atan2 t_1 (- (log1p (expm1 t_0)) t_3))
(if (<= phi2 3e-6)
(atan2
(- (* (sin lambda1) (cos lambda2)) (* (cos lambda1) (sin lambda2)))
(-
t_0
(*
t_2
(fma (cos lambda2) (cos lambda1) (* (sin lambda1) (sin lambda2))))))
(atan2 t_1 (- (log (exp t_0)) t_3))))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = cos(phi1) * sin(phi2);
double t_1 = fma(sin(lambda1), cos(lambda2), (cos(lambda1) * -sin(lambda2))) * cos(phi2);
double t_2 = cos(phi2) * sin(phi1);
double t_3 = t_2 * cos((lambda1 - lambda2));
double tmp;
if (phi2 <= -7e-6) {
tmp = atan2(t_1, (log1p(expm1(t_0)) - t_3));
} else if (phi2 <= 3e-6) {
tmp = atan2(((sin(lambda1) * cos(lambda2)) - (cos(lambda1) * sin(lambda2))), (t_0 - (t_2 * fma(cos(lambda2), cos(lambda1), (sin(lambda1) * sin(lambda2))))));
} else {
tmp = atan2(t_1, (log(exp(t_0)) - t_3));
}
return tmp;
}
function code(lambda1, lambda2, phi1, phi2) t_0 = Float64(cos(phi1) * sin(phi2)) t_1 = Float64(fma(sin(lambda1), cos(lambda2), Float64(cos(lambda1) * Float64(-sin(lambda2)))) * cos(phi2)) t_2 = Float64(cos(phi2) * sin(phi1)) t_3 = Float64(t_2 * cos(Float64(lambda1 - lambda2))) tmp = 0.0 if (phi2 <= -7e-6) tmp = atan(t_1, Float64(log1p(expm1(t_0)) - t_3)); elseif (phi2 <= 3e-6) tmp = atan(Float64(Float64(sin(lambda1) * cos(lambda2)) - Float64(cos(lambda1) * sin(lambda2))), Float64(t_0 - Float64(t_2 * fma(cos(lambda2), cos(lambda1), Float64(sin(lambda1) * sin(lambda2)))))); else tmp = atan(t_1, Float64(log(exp(t_0)) - t_3)); 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[lambda1], $MachinePrecision] * N[Cos[lambda2], $MachinePrecision] + N[(N[Cos[lambda1], $MachinePrecision] * (-N[Sin[lambda2], $MachinePrecision])), $MachinePrecision]), $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(N[Cos[phi2], $MachinePrecision] * N[Sin[phi1], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = N[(t$95$2 * N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[phi2, -7e-6], N[ArcTan[t$95$1 / N[(N[Log[1 + N[(Exp[t$95$0] - 1), $MachinePrecision]], $MachinePrecision] - t$95$3), $MachinePrecision]], $MachinePrecision], If[LessEqual[phi2, 3e-6], N[ArcTan[N[(N[(N[Sin[lambda1], $MachinePrecision] * N[Cos[lambda2], $MachinePrecision]), $MachinePrecision] - N[(N[Cos[lambda1], $MachinePrecision] * N[Sin[lambda2], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(t$95$0 - N[(t$95$2 * N[(N[Cos[lambda2], $MachinePrecision] * N[Cos[lambda1], $MachinePrecision] + N[(N[Sin[lambda1], $MachinePrecision] * N[Sin[lambda2], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], N[ArcTan[t$95$1 / N[(N[Log[N[Exp[t$95$0], $MachinePrecision]], $MachinePrecision] - t$95$3), $MachinePrecision]], $MachinePrecision]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \cos \phi_1 \cdot \sin \phi_2\\
t_1 := \mathsf{fma}\left(\sin \lambda_1, \cos \lambda_2, \cos \lambda_1 \cdot \left(-\sin \lambda_2\right)\right) \cdot \cos \phi_2\\
t_2 := \cos \phi_2 \cdot \sin \phi_1\\
t_3 := t\_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)\\
\mathbf{if}\;\phi_2 \leq -7 \cdot 10^{-6}:\\
\;\;\;\;\tan^{-1}_* \frac{t\_1}{\mathsf{log1p}\left(\mathsf{expm1}\left(t\_0\right)\right) - t\_3}\\
\mathbf{elif}\;\phi_2 \leq 3 \cdot 10^{-6}:\\
\;\;\;\;\tan^{-1}_* \frac{\sin \lambda_1 \cdot \cos \lambda_2 - \cos \lambda_1 \cdot \sin \lambda_2}{t\_0 - t\_2 \cdot \mathsf{fma}\left(\cos \lambda_2, \cos \lambda_1, \sin \lambda_1 \cdot \sin \lambda_2\right)}\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1}_* \frac{t\_1}{\log \left(e^{t\_0}\right) - t\_3}\\
\end{array}
\end{array}
if phi2 < -6.99999999999999989e-6Initial program 80.2%
sin-diff93.4%
fma-neg93.4%
Applied egg-rr93.4%
log1p-expm1-u93.4%
Applied egg-rr93.4%
if -6.99999999999999989e-6 < phi2 < 3.0000000000000001e-6Initial program 83.8%
sin-diff88.8%
fma-neg88.8%
Applied egg-rr88.8%
cos-diff99.8%
*-commutative99.8%
Applied egg-rr99.8%
fma-define99.8%
Simplified99.8%
Taylor expanded in phi2 around 0 99.8%
if 3.0000000000000001e-6 < phi2 Initial program 77.4%
sin-diff88.2%
fma-neg88.2%
Applied egg-rr88.2%
add-log-exp88.3%
Applied egg-rr88.3%
Final simplification95.3%
(FPCore (lambda1 lambda2 phi1 phi2)
:precision binary64
(atan2
(*
(cos phi2)
(- (* (sin lambda1) (cos lambda2)) (* (cos lambda1) (sin lambda2))))
(-
(* (cos phi1) (sin phi2))
(*
(* (cos phi2) (sin phi1))
(+ (* (sin lambda1) (sin lambda2)) (* (cos lambda2) (cos lambda1)))))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
return atan2((cos(phi2) * ((sin(lambda1) * cos(lambda2)) - (cos(lambda1) * sin(lambda2)))), ((cos(phi1) * sin(phi2)) - ((cos(phi2) * sin(phi1)) * ((sin(lambda1) * sin(lambda2)) + (cos(lambda2) * cos(lambda1))))));
}
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) * cos(lambda2)) - (cos(lambda1) * sin(lambda2)))), ((cos(phi1) * sin(phi2)) - ((cos(phi2) * sin(phi1)) * ((sin(lambda1) * sin(lambda2)) + (cos(lambda2) * cos(lambda1))))))
end function
public static double code(double lambda1, double lambda2, double phi1, double phi2) {
return Math.atan2((Math.cos(phi2) * ((Math.sin(lambda1) * Math.cos(lambda2)) - (Math.cos(lambda1) * Math.sin(lambda2)))), ((Math.cos(phi1) * Math.sin(phi2)) - ((Math.cos(phi2) * Math.sin(phi1)) * ((Math.sin(lambda1) * Math.sin(lambda2)) + (Math.cos(lambda2) * Math.cos(lambda1))))));
}
def code(lambda1, lambda2, phi1, phi2): return math.atan2((math.cos(phi2) * ((math.sin(lambda1) * math.cos(lambda2)) - (math.cos(lambda1) * math.sin(lambda2)))), ((math.cos(phi1) * math.sin(phi2)) - ((math.cos(phi2) * math.sin(phi1)) * ((math.sin(lambda1) * math.sin(lambda2)) + (math.cos(lambda2) * math.cos(lambda1))))))
function code(lambda1, lambda2, phi1, phi2) return atan(Float64(cos(phi2) * Float64(Float64(sin(lambda1) * cos(lambda2)) - Float64(cos(lambda1) * sin(lambda2)))), Float64(Float64(cos(phi1) * sin(phi2)) - Float64(Float64(cos(phi2) * sin(phi1)) * Float64(Float64(sin(lambda1) * sin(lambda2)) + Float64(cos(lambda2) * cos(lambda1)))))) end
function tmp = code(lambda1, lambda2, phi1, phi2) tmp = atan2((cos(phi2) * ((sin(lambda1) * cos(lambda2)) - (cos(lambda1) * sin(lambda2)))), ((cos(phi1) * sin(phi2)) - ((cos(phi2) * sin(phi1)) * ((sin(lambda1) * sin(lambda2)) + (cos(lambda2) * cos(lambda1)))))); end
code[lambda1_, lambda2_, phi1_, phi2_] := N[ArcTan[N[(N[Cos[phi2], $MachinePrecision] * N[(N[(N[Sin[lambda1], $MachinePrecision] * N[Cos[lambda2], $MachinePrecision]), $MachinePrecision] - N[(N[Cos[lambda1], $MachinePrecision] * N[Sin[lambda2], $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[(N[Sin[lambda1], $MachinePrecision] * N[Sin[lambda2], $MachinePrecision]), $MachinePrecision] + N[(N[Cos[lambda2], $MachinePrecision] * N[Cos[lambda1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
\\
\tan^{-1}_* \frac{\cos \phi_2 \cdot \left(\sin \lambda_1 \cdot \cos \lambda_2 - \cos \lambda_1 \cdot \sin \lambda_2\right)}{\cos \phi_1 \cdot \sin \phi_2 - \left(\cos \phi_2 \cdot \sin \phi_1\right) \cdot \left(\sin \lambda_1 \cdot \sin \lambda_2 + \cos \lambda_2 \cdot \cos \lambda_1\right)}
\end{array}
Initial program 81.3%
cos-diff81.7%
+-commutative81.7%
*-commutative81.7%
Applied egg-rr81.7%
sin-diff54.4%
Applied egg-rr99.7%
Final simplification99.7%
(FPCore (lambda1 lambda2 phi1 phi2) :precision binary64 (atan2 (* (fma (sin lambda1) (cos lambda2) (* (cos lambda1) (- (sin lambda2)))) (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(lambda1), cos(lambda2), (cos(lambda1) * -sin(lambda2))) * cos(phi2)), ((cos(phi1) * sin(phi2)) - ((cos(phi2) * sin(phi1)) * cos((lambda1 - lambda2)))));
}
function code(lambda1, lambda2, phi1, phi2) return atan(Float64(fma(sin(lambda1), cos(lambda2), Float64(cos(lambda1) * Float64(-sin(lambda2)))) * 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[lambda1], $MachinePrecision] * N[Cos[lambda2], $MachinePrecision] + N[(N[Cos[lambda1], $MachinePrecision] * (-N[Sin[lambda2], $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_1, \cos \lambda_2, \cos \lambda_1 \cdot \left(-\sin \lambda_2\right)\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 81.3%
sin-diff89.7%
fma-neg89.7%
Applied egg-rr89.7%
Final simplification89.7%
(FPCore (lambda1 lambda2 phi1 phi2)
:precision binary64
(if (or (<= phi1 -8.2e-14) (not (<= phi1 4.5e-48)))
(atan2
(* (cos phi2) (sin (- lambda1 lambda2)))
(-
(* (cos phi1) (sin phi2))
(*
(* (cos phi2) (sin phi1))
(+ (* (sin lambda1) (sin lambda2)) (* (cos lambda2) (cos lambda1))))))
(atan2
(*
(fma (sin lambda1) (cos lambda2) (* (cos lambda1) (- (sin lambda2))))
(cos phi2))
(- (sin phi2) (* phi1 (cos (- lambda2 lambda1)))))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
double tmp;
if ((phi1 <= -8.2e-14) || !(phi1 <= 4.5e-48)) {
tmp = atan2((cos(phi2) * sin((lambda1 - lambda2))), ((cos(phi1) * sin(phi2)) - ((cos(phi2) * sin(phi1)) * ((sin(lambda1) * sin(lambda2)) + (cos(lambda2) * cos(lambda1))))));
} else {
tmp = atan2((fma(sin(lambda1), cos(lambda2), (cos(lambda1) * -sin(lambda2))) * cos(phi2)), (sin(phi2) - (phi1 * cos((lambda2 - lambda1)))));
}
return tmp;
}
function code(lambda1, lambda2, phi1, phi2) tmp = 0.0 if ((phi1 <= -8.2e-14) || !(phi1 <= 4.5e-48)) tmp = atan(Float64(cos(phi2) * sin(Float64(lambda1 - lambda2))), Float64(Float64(cos(phi1) * sin(phi2)) - Float64(Float64(cos(phi2) * sin(phi1)) * Float64(Float64(sin(lambda1) * sin(lambda2)) + Float64(cos(lambda2) * cos(lambda1)))))); else tmp = atan(Float64(fma(sin(lambda1), cos(lambda2), Float64(cos(lambda1) * Float64(-sin(lambda2)))) * cos(phi2)), Float64(sin(phi2) - Float64(phi1 * cos(Float64(lambda2 - lambda1))))); end return tmp end
code[lambda1_, lambda2_, phi1_, phi2_] := If[Or[LessEqual[phi1, -8.2e-14], N[Not[LessEqual[phi1, 4.5e-48]], $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[(N[Cos[phi2], $MachinePrecision] * N[Sin[phi1], $MachinePrecision]), $MachinePrecision] * N[(N[(N[Sin[lambda1], $MachinePrecision] * N[Sin[lambda2], $MachinePrecision]), $MachinePrecision] + N[(N[Cos[lambda2], $MachinePrecision] * N[Cos[lambda1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], N[ArcTan[N[(N[(N[Sin[lambda1], $MachinePrecision] * N[Cos[lambda2], $MachinePrecision] + N[(N[Cos[lambda1], $MachinePrecision] * (-N[Sin[lambda2], $MachinePrecision])), $MachinePrecision]), $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision] / N[(N[Sin[phi2], $MachinePrecision] - N[(phi1 * N[Cos[N[(lambda2 - lambda1), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\phi_1 \leq -8.2 \cdot 10^{-14} \lor \neg \left(\phi_1 \leq 4.5 \cdot 10^{-48}\right):\\
\;\;\;\;\tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 \cdot \sin \phi_2 - \left(\cos \phi_2 \cdot \sin \phi_1\right) \cdot \left(\sin \lambda_1 \cdot \sin \lambda_2 + \cos \lambda_2 \cdot \cos \lambda_1\right)}\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1}_* \frac{\mathsf{fma}\left(\sin \lambda_1, \cos \lambda_2, \cos \lambda_1 \cdot \left(-\sin \lambda_2\right)\right) \cdot \cos \phi_2}{\sin \phi_2 - \phi_1 \cdot \cos \left(\lambda_2 - \lambda_1\right)}\\
\end{array}
\end{array}
if phi1 < -8.2000000000000004e-14 or 4.49999999999999988e-48 < phi1 Initial program 80.2%
cos-diff80.9%
+-commutative80.9%
*-commutative80.9%
Applied egg-rr80.9%
if -8.2000000000000004e-14 < phi1 < 4.49999999999999988e-48Initial program 82.8%
*-commutative82.8%
associate-*l*82.8%
Simplified82.8%
Taylor expanded in phi1 around 0 82.8%
sub-neg82.8%
remove-double-neg82.8%
mul-1-neg82.8%
distribute-neg-in82.8%
+-commutative82.8%
*-commutative82.8%
cos-neg82.8%
mul-1-neg82.8%
unsub-neg82.8%
Simplified82.8%
Taylor expanded in phi1 around 0 82.8%
Taylor expanded in phi2 around 0 82.8%
sin-diff99.8%
fma-neg99.8%
Applied egg-rr99.8%
Final simplification88.9%
(FPCore (lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (* (cos phi1) (sin phi2)))
(t_1 (* (cos phi2) (sin (- lambda1 lambda2))))
(t_2 (cos (- lambda1 lambda2)))
(t_3 (* (* (cos phi2) (sin phi1)) t_2)))
(if (<= phi1 -1.85e+18)
(atan2 t_1 (- t_0 (* (cos phi2) (log1p (expm1 (* (sin phi1) t_2))))))
(if (<= phi1 4e-48)
(atan2
(*
(fma (sin lambda1) (cos lambda2) (* (cos lambda1) (- (sin lambda2))))
(cos phi2))
(- (sin phi2) t_3))
(atan2 t_1 (- t_0 t_3))))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = cos(phi1) * sin(phi2);
double t_1 = cos(phi2) * sin((lambda1 - lambda2));
double t_2 = cos((lambda1 - lambda2));
double t_3 = (cos(phi2) * sin(phi1)) * t_2;
double tmp;
if (phi1 <= -1.85e+18) {
tmp = atan2(t_1, (t_0 - (cos(phi2) * log1p(expm1((sin(phi1) * t_2))))));
} else if (phi1 <= 4e-48) {
tmp = atan2((fma(sin(lambda1), cos(lambda2), (cos(lambda1) * -sin(lambda2))) * cos(phi2)), (sin(phi2) - t_3));
} else {
tmp = atan2(t_1, (t_0 - t_3));
}
return tmp;
}
function code(lambda1, lambda2, phi1, phi2) t_0 = Float64(cos(phi1) * sin(phi2)) t_1 = Float64(cos(phi2) * sin(Float64(lambda1 - lambda2))) t_2 = cos(Float64(lambda1 - lambda2)) t_3 = Float64(Float64(cos(phi2) * sin(phi1)) * t_2) tmp = 0.0 if (phi1 <= -1.85e+18) tmp = atan(t_1, Float64(t_0 - Float64(cos(phi2) * log1p(expm1(Float64(sin(phi1) * t_2)))))); elseif (phi1 <= 4e-48) tmp = atan(Float64(fma(sin(lambda1), cos(lambda2), Float64(cos(lambda1) * Float64(-sin(lambda2)))) * cos(phi2)), Float64(sin(phi2) - t_3)); else tmp = atan(t_1, Float64(t_0 - t_3)); 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[phi2], $MachinePrecision] * N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$3 = N[(N[(N[Cos[phi2], $MachinePrecision] * N[Sin[phi1], $MachinePrecision]), $MachinePrecision] * t$95$2), $MachinePrecision]}, If[LessEqual[phi1, -1.85e+18], N[ArcTan[t$95$1 / N[(t$95$0 - N[(N[Cos[phi2], $MachinePrecision] * N[Log[1 + N[(Exp[N[(N[Sin[phi1], $MachinePrecision] * t$95$2), $MachinePrecision]] - 1), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], If[LessEqual[phi1, 4e-48], N[ArcTan[N[(N[(N[Sin[lambda1], $MachinePrecision] * N[Cos[lambda2], $MachinePrecision] + N[(N[Cos[lambda1], $MachinePrecision] * (-N[Sin[lambda2], $MachinePrecision])), $MachinePrecision]), $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision] / N[(N[Sin[phi2], $MachinePrecision] - t$95$3), $MachinePrecision]], $MachinePrecision], N[ArcTan[t$95$1 / N[(t$95$0 - t$95$3), $MachinePrecision]], $MachinePrecision]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \cos \phi_1 \cdot \sin \phi_2\\
t_1 := \cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)\\
t_2 := \cos \left(\lambda_1 - \lambda_2\right)\\
t_3 := \left(\cos \phi_2 \cdot \sin \phi_1\right) \cdot t\_2\\
\mathbf{if}\;\phi_1 \leq -1.85 \cdot 10^{+18}:\\
\;\;\;\;\tan^{-1}_* \frac{t\_1}{t\_0 - \cos \phi_2 \cdot \mathsf{log1p}\left(\mathsf{expm1}\left(\sin \phi_1 \cdot t\_2\right)\right)}\\
\mathbf{elif}\;\phi_1 \leq 4 \cdot 10^{-48}:\\
\;\;\;\;\tan^{-1}_* \frac{\mathsf{fma}\left(\sin \lambda_1, \cos \lambda_2, \cos \lambda_1 \cdot \left(-\sin \lambda_2\right)\right) \cdot \cos \phi_2}{\sin \phi_2 - t\_3}\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1}_* \frac{t\_1}{t\_0 - t\_3}\\
\end{array}
\end{array}
if phi1 < -1.85e18Initial program 80.4%
*-commutative80.4%
associate-*l*80.4%
Simplified80.4%
log1p-expm1-u80.5%
*-commutative80.5%
Applied egg-rr80.5%
if -1.85e18 < phi1 < 3.9999999999999999e-48Initial program 82.6%
sin-diff99.0%
fma-neg99.0%
Applied egg-rr99.0%
Taylor expanded in phi1 around 0 99.0%
if 3.9999999999999999e-48 < phi1 Initial program 80.2%
Final simplification88.5%
(FPCore (lambda1 lambda2 phi1 phi2) :precision binary64 (atan2 (* (cos phi2) (- (* (sin lambda1) (cos lambda2)) (* (cos lambda1) (sin lambda2)))) (- (* (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) * ((sin(lambda1) * cos(lambda2)) - (cos(lambda1) * sin(lambda2)))), ((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) * ((sin(lambda1) * cos(lambda2)) - (cos(lambda1) * sin(lambda2)))), ((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.sin(lambda1) * Math.cos(lambda2)) - (Math.cos(lambda1) * Math.sin(lambda2)))), ((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.sin(lambda1) * math.cos(lambda2)) - (math.cos(lambda1) * math.sin(lambda2)))), ((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(sin(lambda1) * cos(lambda2)) - Float64(cos(lambda1) * sin(lambda2)))), Float64(Float64(cos(phi1) * sin(phi2)) - Float64(Float64(cos(phi2) * sin(phi1)) * cos(Float64(lambda1 - lambda2))))) end
function tmp = code(lambda1, lambda2, phi1, phi2) tmp = atan2((cos(phi2) * ((sin(lambda1) * cos(lambda2)) - (cos(lambda1) * sin(lambda2)))), ((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[Sin[lambda1], $MachinePrecision] * N[Cos[lambda2], $MachinePrecision]), $MachinePrecision] - N[(N[Cos[lambda1], $MachinePrecision] * N[Sin[lambda2], $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[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
\\
\tan^{-1}_* \frac{\cos \phi_2 \cdot \left(\sin \lambda_1 \cdot \cos \lambda_2 - \cos \lambda_1 \cdot \sin \lambda_2\right)}{\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 81.3%
sin-diff54.4%
Applied egg-rr89.7%
Final simplification89.7%
(FPCore (lambda1 lambda2 phi1 phi2)
:precision binary64
(if (or (<= phi1 -1.15e-10) (not (<= phi1 4.5e-48)))
(atan2
(* (cos phi2) (sin (- lambda1 lambda2)))
(-
(* (cos phi1) (sin phi2))
(* (* (cos phi2) (sin phi1)) (cos (- lambda1 lambda2)))))
(atan2
(*
(fma (sin lambda1) (cos lambda2) (* (cos lambda1) (- (sin lambda2))))
(cos phi2))
(- (sin phi2) (* phi1 (cos (- lambda2 lambda1)))))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
double tmp;
if ((phi1 <= -1.15e-10) || !(phi1 <= 4.5e-48)) {
tmp = atan2((cos(phi2) * sin((lambda1 - lambda2))), ((cos(phi1) * sin(phi2)) - ((cos(phi2) * sin(phi1)) * cos((lambda1 - lambda2)))));
} else {
tmp = atan2((fma(sin(lambda1), cos(lambda2), (cos(lambda1) * -sin(lambda2))) * cos(phi2)), (sin(phi2) - (phi1 * cos((lambda2 - lambda1)))));
}
return tmp;
}
function code(lambda1, lambda2, phi1, phi2) tmp = 0.0 if ((phi1 <= -1.15e-10) || !(phi1 <= 4.5e-48)) tmp = atan(Float64(cos(phi2) * sin(Float64(lambda1 - lambda2))), Float64(Float64(cos(phi1) * sin(phi2)) - Float64(Float64(cos(phi2) * sin(phi1)) * cos(Float64(lambda1 - lambda2))))); else tmp = atan(Float64(fma(sin(lambda1), cos(lambda2), Float64(cos(lambda1) * Float64(-sin(lambda2)))) * cos(phi2)), Float64(sin(phi2) - Float64(phi1 * cos(Float64(lambda2 - lambda1))))); end return tmp end
code[lambda1_, lambda2_, phi1_, phi2_] := If[Or[LessEqual[phi1, -1.15e-10], N[Not[LessEqual[phi1, 4.5e-48]], $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[(N[Cos[phi2], $MachinePrecision] * N[Sin[phi1], $MachinePrecision]), $MachinePrecision] * N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], N[ArcTan[N[(N[(N[Sin[lambda1], $MachinePrecision] * N[Cos[lambda2], $MachinePrecision] + N[(N[Cos[lambda1], $MachinePrecision] * (-N[Sin[lambda2], $MachinePrecision])), $MachinePrecision]), $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision] / N[(N[Sin[phi2], $MachinePrecision] - N[(phi1 * N[Cos[N[(lambda2 - lambda1), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\phi_1 \leq -1.15 \cdot 10^{-10} \lor \neg \left(\phi_1 \leq 4.5 \cdot 10^{-48}\right):\\
\;\;\;\;\tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 \cdot \sin \phi_2 - \left(\cos \phi_2 \cdot \sin \phi_1\right) \cdot \cos \left(\lambda_1 - \lambda_2\right)}\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1}_* \frac{\mathsf{fma}\left(\sin \lambda_1, \cos \lambda_2, \cos \lambda_1 \cdot \left(-\sin \lambda_2\right)\right) \cdot \cos \phi_2}{\sin \phi_2 - \phi_1 \cdot \cos \left(\lambda_2 - \lambda_1\right)}\\
\end{array}
\end{array}
if phi1 < -1.15000000000000004e-10 or 4.49999999999999988e-48 < phi1 Initial program 80.1%
if -1.15000000000000004e-10 < phi1 < 4.49999999999999988e-48Initial program 82.8%
*-commutative82.8%
associate-*l*82.8%
Simplified82.8%
Taylor expanded in phi1 around 0 82.8%
sub-neg82.8%
remove-double-neg82.8%
mul-1-neg82.8%
distribute-neg-in82.8%
+-commutative82.8%
*-commutative82.8%
cos-neg82.8%
mul-1-neg82.8%
unsub-neg82.8%
Simplified82.8%
Taylor expanded in phi1 around 0 82.8%
Taylor expanded in phi2 around 0 82.8%
sin-diff99.6%
fma-neg99.7%
Applied egg-rr99.7%
Final simplification88.5%
(FPCore (lambda1 lambda2 phi1 phi2)
:precision binary64
(if (or (<= phi1 -1.85e-10) (not (<= phi1 4.5e-48)))
(atan2
(* (cos phi2) (sin (- lambda1 lambda2)))
(-
(* (cos phi1) (sin phi2))
(* (* (cos phi2) (sin phi1)) (cos (- lambda1 lambda2)))))
(atan2
(*
(cos phi2)
(- (* (sin lambda1) (cos lambda2)) (* (cos lambda1) (sin lambda2))))
(- (sin phi2) (* phi1 (cos (- lambda2 lambda1)))))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
double tmp;
if ((phi1 <= -1.85e-10) || !(phi1 <= 4.5e-48)) {
tmp = atan2((cos(phi2) * sin((lambda1 - lambda2))), ((cos(phi1) * sin(phi2)) - ((cos(phi2) * sin(phi1)) * cos((lambda1 - lambda2)))));
} else {
tmp = atan2((cos(phi2) * ((sin(lambda1) * cos(lambda2)) - (cos(lambda1) * sin(lambda2)))), (sin(phi2) - (phi1 * cos((lambda2 - lambda1)))));
}
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 ((phi1 <= (-1.85d-10)) .or. (.not. (phi1 <= 4.5d-48))) then
tmp = atan2((cos(phi2) * sin((lambda1 - lambda2))), ((cos(phi1) * sin(phi2)) - ((cos(phi2) * sin(phi1)) * cos((lambda1 - lambda2)))))
else
tmp = atan2((cos(phi2) * ((sin(lambda1) * cos(lambda2)) - (cos(lambda1) * sin(lambda2)))), (sin(phi2) - (phi1 * cos((lambda2 - lambda1)))))
end if
code = tmp
end function
public static double code(double lambda1, double lambda2, double phi1, double phi2) {
double tmp;
if ((phi1 <= -1.85e-10) || !(phi1 <= 4.5e-48)) {
tmp = Math.atan2((Math.cos(phi2) * Math.sin((lambda1 - lambda2))), ((Math.cos(phi1) * Math.sin(phi2)) - ((Math.cos(phi2) * Math.sin(phi1)) * Math.cos((lambda1 - lambda2)))));
} else {
tmp = Math.atan2((Math.cos(phi2) * ((Math.sin(lambda1) * Math.cos(lambda2)) - (Math.cos(lambda1) * Math.sin(lambda2)))), (Math.sin(phi2) - (phi1 * Math.cos((lambda2 - lambda1)))));
}
return tmp;
}
def code(lambda1, lambda2, phi1, phi2): tmp = 0 if (phi1 <= -1.85e-10) or not (phi1 <= 4.5e-48): tmp = math.atan2((math.cos(phi2) * math.sin((lambda1 - lambda2))), ((math.cos(phi1) * math.sin(phi2)) - ((math.cos(phi2) * math.sin(phi1)) * math.cos((lambda1 - lambda2))))) else: tmp = math.atan2((math.cos(phi2) * ((math.sin(lambda1) * math.cos(lambda2)) - (math.cos(lambda1) * math.sin(lambda2)))), (math.sin(phi2) - (phi1 * math.cos((lambda2 - lambda1))))) return tmp
function code(lambda1, lambda2, phi1, phi2) tmp = 0.0 if ((phi1 <= -1.85e-10) || !(phi1 <= 4.5e-48)) tmp = atan(Float64(cos(phi2) * sin(Float64(lambda1 - lambda2))), Float64(Float64(cos(phi1) * sin(phi2)) - Float64(Float64(cos(phi2) * sin(phi1)) * cos(Float64(lambda1 - lambda2))))); else tmp = atan(Float64(cos(phi2) * Float64(Float64(sin(lambda1) * cos(lambda2)) - Float64(cos(lambda1) * sin(lambda2)))), Float64(sin(phi2) - Float64(phi1 * cos(Float64(lambda2 - lambda1))))); end return tmp end
function tmp_2 = code(lambda1, lambda2, phi1, phi2) tmp = 0.0; if ((phi1 <= -1.85e-10) || ~((phi1 <= 4.5e-48))) tmp = atan2((cos(phi2) * sin((lambda1 - lambda2))), ((cos(phi1) * sin(phi2)) - ((cos(phi2) * sin(phi1)) * cos((lambda1 - lambda2))))); else tmp = atan2((cos(phi2) * ((sin(lambda1) * cos(lambda2)) - (cos(lambda1) * sin(lambda2)))), (sin(phi2) - (phi1 * cos((lambda2 - lambda1))))); end tmp_2 = tmp; end
code[lambda1_, lambda2_, phi1_, phi2_] := If[Or[LessEqual[phi1, -1.85e-10], N[Not[LessEqual[phi1, 4.5e-48]], $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[(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[(N[(N[Sin[lambda1], $MachinePrecision] * N[Cos[lambda2], $MachinePrecision]), $MachinePrecision] - N[(N[Cos[lambda1], $MachinePrecision] * N[Sin[lambda2], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(N[Sin[phi2], $MachinePrecision] - N[(phi1 * N[Cos[N[(lambda2 - lambda1), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\phi_1 \leq -1.85 \cdot 10^{-10} \lor \neg \left(\phi_1 \leq 4.5 \cdot 10^{-48}\right):\\
\;\;\;\;\tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 \cdot \sin \phi_2 - \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(\sin \lambda_1 \cdot \cos \lambda_2 - \cos \lambda_1 \cdot \sin \lambda_2\right)}{\sin \phi_2 - \phi_1 \cdot \cos \left(\lambda_2 - \lambda_1\right)}\\
\end{array}
\end{array}
if phi1 < -1.85000000000000007e-10 or 4.49999999999999988e-48 < phi1 Initial program 80.1%
if -1.85000000000000007e-10 < phi1 < 4.49999999999999988e-48Initial program 82.8%
*-commutative82.8%
associate-*l*82.8%
Simplified82.8%
Taylor expanded in phi1 around 0 82.8%
sub-neg82.8%
remove-double-neg82.8%
mul-1-neg82.8%
distribute-neg-in82.8%
+-commutative82.8%
*-commutative82.8%
cos-neg82.8%
mul-1-neg82.8%
unsub-neg82.8%
Simplified82.8%
Taylor expanded in phi1 around 0 82.8%
Taylor expanded in phi2 around 0 82.8%
sin-diff99.6%
Applied egg-rr99.6%
Final simplification88.5%
(FPCore (lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (* (cos phi1) (sin phi2)))
(t_1 (* (cos phi2) (sin phi1)))
(t_2 (sin (- lambda1 lambda2))))
(if (or (<= phi2 -1.28e-5) (not (<= phi2 3200000000.0)))
(atan2 (* (cos phi2) t_2) (- t_0 (* (cos lambda1) t_1)))
(atan2 t_2 (- t_0 (* t_1 (cos (- lambda1 lambda2))))))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = cos(phi1) * sin(phi2);
double t_1 = cos(phi2) * sin(phi1);
double t_2 = sin((lambda1 - lambda2));
double tmp;
if ((phi2 <= -1.28e-5) || !(phi2 <= 3200000000.0)) {
tmp = atan2((cos(phi2) * t_2), (t_0 - (cos(lambda1) * t_1)));
} else {
tmp = atan2(t_2, (t_0 - (t_1 * 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) :: t_1
real(8) :: t_2
real(8) :: tmp
t_0 = cos(phi1) * sin(phi2)
t_1 = cos(phi2) * sin(phi1)
t_2 = sin((lambda1 - lambda2))
if ((phi2 <= (-1.28d-5)) .or. (.not. (phi2 <= 3200000000.0d0))) then
tmp = atan2((cos(phi2) * t_2), (t_0 - (cos(lambda1) * t_1)))
else
tmp = atan2(t_2, (t_0 - (t_1 * 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.cos(phi1) * Math.sin(phi2);
double t_1 = Math.cos(phi2) * Math.sin(phi1);
double t_2 = Math.sin((lambda1 - lambda2));
double tmp;
if ((phi2 <= -1.28e-5) || !(phi2 <= 3200000000.0)) {
tmp = Math.atan2((Math.cos(phi2) * t_2), (t_0 - (Math.cos(lambda1) * t_1)));
} else {
tmp = Math.atan2(t_2, (t_0 - (t_1 * Math.cos((lambda1 - lambda2)))));
}
return tmp;
}
def code(lambda1, lambda2, phi1, phi2): t_0 = math.cos(phi1) * math.sin(phi2) t_1 = math.cos(phi2) * math.sin(phi1) t_2 = math.sin((lambda1 - lambda2)) tmp = 0 if (phi2 <= -1.28e-5) or not (phi2 <= 3200000000.0): tmp = math.atan2((math.cos(phi2) * t_2), (t_0 - (math.cos(lambda1) * t_1))) else: tmp = math.atan2(t_2, (t_0 - (t_1 * math.cos((lambda1 - lambda2))))) return tmp
function code(lambda1, lambda2, phi1, phi2) t_0 = Float64(cos(phi1) * sin(phi2)) t_1 = Float64(cos(phi2) * sin(phi1)) t_2 = sin(Float64(lambda1 - lambda2)) tmp = 0.0 if ((phi2 <= -1.28e-5) || !(phi2 <= 3200000000.0)) tmp = atan(Float64(cos(phi2) * t_2), Float64(t_0 - Float64(cos(lambda1) * t_1))); else tmp = atan(t_2, Float64(t_0 - Float64(t_1 * cos(Float64(lambda1 - lambda2))))); end return tmp end
function tmp_2 = code(lambda1, lambda2, phi1, phi2) t_0 = cos(phi1) * sin(phi2); t_1 = cos(phi2) * sin(phi1); t_2 = sin((lambda1 - lambda2)); tmp = 0.0; if ((phi2 <= -1.28e-5) || ~((phi2 <= 3200000000.0))) tmp = atan2((cos(phi2) * t_2), (t_0 - (cos(lambda1) * t_1))); else tmp = atan2(t_2, (t_0 - (t_1 * cos((lambda1 - lambda2))))); end tmp_2 = 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[phi2], $MachinePrecision] * N[Sin[phi1], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]}, If[Or[LessEqual[phi2, -1.28e-5], N[Not[LessEqual[phi2, 3200000000.0]], $MachinePrecision]], N[ArcTan[N[(N[Cos[phi2], $MachinePrecision] * t$95$2), $MachinePrecision] / N[(t$95$0 - N[(N[Cos[lambda1], $MachinePrecision] * t$95$1), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], N[ArcTan[t$95$2 / N[(t$95$0 - N[(t$95$1 * 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 \phi_2 \cdot \sin \phi_1\\
t_2 := \sin \left(\lambda_1 - \lambda_2\right)\\
\mathbf{if}\;\phi_2 \leq -1.28 \cdot 10^{-5} \lor \neg \left(\phi_2 \leq 3200000000\right):\\
\;\;\;\;\tan^{-1}_* \frac{\cos \phi_2 \cdot t\_2}{t\_0 - \cos \lambda_1 \cdot t\_1}\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1}_* \frac{t\_2}{t\_0 - t\_1 \cdot \cos \left(\lambda_1 - \lambda_2\right)}\\
\end{array}
\end{array}
if phi2 < -1.2799999999999999e-5 or 3.2e9 < phi2 Initial program 78.6%
Taylor expanded in lambda2 around 0 67.6%
if -1.2799999999999999e-5 < phi2 < 3.2e9Initial program 83.9%
add-sqr-sqrt43.9%
pow243.9%
Applied egg-rr43.9%
Taylor expanded in phi2 around 0 83.2%
associate-*r*83.2%
distribute-rgt1-in83.2%
Simplified83.2%
Taylor expanded in phi2 around 0 83.4%
Final simplification75.6%
(FPCore (lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (* (cos phi1) (sin phi2)))
(t_1 (* (cos phi2) (sin phi1)))
(t_2 (* (cos phi2) (sin (- lambda1 lambda2)))))
(if (or (<= lambda1 -6e-13) (not (<= lambda1 3.8e-5)))
(atan2 t_2 (- t_0 (* (cos lambda1) t_1)))
(atan2 t_2 (- t_0 (* (cos lambda2) t_1))))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = cos(phi1) * sin(phi2);
double t_1 = cos(phi2) * sin(phi1);
double t_2 = cos(phi2) * sin((lambda1 - lambda2));
double tmp;
if ((lambda1 <= -6e-13) || !(lambda1 <= 3.8e-5)) {
tmp = atan2(t_2, (t_0 - (cos(lambda1) * t_1)));
} else {
tmp = atan2(t_2, (t_0 - (cos(lambda2) * t_1)));
}
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) :: t_2
real(8) :: tmp
t_0 = cos(phi1) * sin(phi2)
t_1 = cos(phi2) * sin(phi1)
t_2 = cos(phi2) * sin((lambda1 - lambda2))
if ((lambda1 <= (-6d-13)) .or. (.not. (lambda1 <= 3.8d-5))) then
tmp = atan2(t_2, (t_0 - (cos(lambda1) * t_1)))
else
tmp = atan2(t_2, (t_0 - (cos(lambda2) * t_1)))
end if
code = tmp
end function
public static double code(double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = Math.cos(phi1) * Math.sin(phi2);
double t_1 = Math.cos(phi2) * Math.sin(phi1);
double t_2 = Math.cos(phi2) * Math.sin((lambda1 - lambda2));
double tmp;
if ((lambda1 <= -6e-13) || !(lambda1 <= 3.8e-5)) {
tmp = Math.atan2(t_2, (t_0 - (Math.cos(lambda1) * t_1)));
} else {
tmp = Math.atan2(t_2, (t_0 - (Math.cos(lambda2) * t_1)));
}
return tmp;
}
def code(lambda1, lambda2, phi1, phi2): t_0 = math.cos(phi1) * math.sin(phi2) t_1 = math.cos(phi2) * math.sin(phi1) t_2 = math.cos(phi2) * math.sin((lambda1 - lambda2)) tmp = 0 if (lambda1 <= -6e-13) or not (lambda1 <= 3.8e-5): tmp = math.atan2(t_2, (t_0 - (math.cos(lambda1) * t_1))) else: tmp = math.atan2(t_2, (t_0 - (math.cos(lambda2) * t_1))) return tmp
function code(lambda1, lambda2, phi1, phi2) t_0 = Float64(cos(phi1) * sin(phi2)) t_1 = Float64(cos(phi2) * sin(phi1)) t_2 = Float64(cos(phi2) * sin(Float64(lambda1 - lambda2))) tmp = 0.0 if ((lambda1 <= -6e-13) || !(lambda1 <= 3.8e-5)) tmp = atan(t_2, Float64(t_0 - Float64(cos(lambda1) * t_1))); else tmp = atan(t_2, Float64(t_0 - Float64(cos(lambda2) * t_1))); end return tmp end
function tmp_2 = code(lambda1, lambda2, phi1, phi2) t_0 = cos(phi1) * sin(phi2); t_1 = cos(phi2) * sin(phi1); t_2 = cos(phi2) * sin((lambda1 - lambda2)); tmp = 0.0; if ((lambda1 <= -6e-13) || ~((lambda1 <= 3.8e-5))) tmp = atan2(t_2, (t_0 - (cos(lambda1) * t_1))); else tmp = atan2(t_2, (t_0 - (cos(lambda2) * t_1))); end tmp_2 = 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[phi2], $MachinePrecision] * N[Sin[phi1], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(N[Cos[phi2], $MachinePrecision] * N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, If[Or[LessEqual[lambda1, -6e-13], N[Not[LessEqual[lambda1, 3.8e-5]], $MachinePrecision]], N[ArcTan[t$95$2 / N[(t$95$0 - N[(N[Cos[lambda1], $MachinePrecision] * t$95$1), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], N[ArcTan[t$95$2 / N[(t$95$0 - N[(N[Cos[lambda2], $MachinePrecision] * t$95$1), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \cos \phi_1 \cdot \sin \phi_2\\
t_1 := \cos \phi_2 \cdot \sin \phi_1\\
t_2 := \cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)\\
\mathbf{if}\;\lambda_1 \leq -6 \cdot 10^{-13} \lor \neg \left(\lambda_1 \leq 3.8 \cdot 10^{-5}\right):\\
\;\;\;\;\tan^{-1}_* \frac{t\_2}{t\_0 - \cos \lambda_1 \cdot t\_1}\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1}_* \frac{t\_2}{t\_0 - \cos \lambda_2 \cdot t\_1}\\
\end{array}
\end{array}
if lambda1 < -5.99999999999999968e-13 or 3.8000000000000002e-5 < lambda1 Initial program 64.6%
Taylor expanded in lambda2 around 0 64.4%
if -5.99999999999999968e-13 < lambda1 < 3.8000000000000002e-5Initial program 99.6%
Taylor expanded in lambda1 around 0 99.6%
cos-neg99.6%
Simplified99.6%
Final simplification81.2%
(FPCore (lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (* phi1 (cos (- lambda2 lambda1))))
(t_1 (sin (- lambda1 lambda2)))
(t_2 (* (cos phi2) t_1))
(t_3
(-
(* (cos phi1) (sin phi2))
(* (* (cos phi2) (sin phi1)) (cos (- lambda1 lambda2))))))
(if (<= phi2 -9.5e+197)
(atan2 t_2 (- (sin phi2) (sqrt (pow t_0 2.0))))
(if (<= phi2 -7e-25)
(atan2 (* (sin lambda1) (cos phi2)) t_3)
(if (<= phi2 -1.46e-47)
(atan2
(*
(cos phi2)
(-
(* lambda1 (+ (cos lambda2) (* lambda1 (* (sin lambda2) 0.5))))
(sin lambda2)))
(- (sin phi2) t_0))
(if (<= phi2 3200000000.0)
(atan2 t_1 t_3)
(atan2
t_2
(-
(sin phi2)
(/
(*
(cos phi2)
(+
(sin (+ lambda2 (- phi1 lambda1)))
(sin (+ phi1 (- lambda1 lambda2)))))
2.0)))))))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = phi1 * cos((lambda2 - lambda1));
double t_1 = sin((lambda1 - lambda2));
double t_2 = cos(phi2) * t_1;
double t_3 = (cos(phi1) * sin(phi2)) - ((cos(phi2) * sin(phi1)) * cos((lambda1 - lambda2)));
double tmp;
if (phi2 <= -9.5e+197) {
tmp = atan2(t_2, (sin(phi2) - sqrt(pow(t_0, 2.0))));
} else if (phi2 <= -7e-25) {
tmp = atan2((sin(lambda1) * cos(phi2)), t_3);
} else if (phi2 <= -1.46e-47) {
tmp = atan2((cos(phi2) * ((lambda1 * (cos(lambda2) + (lambda1 * (sin(lambda2) * 0.5)))) - sin(lambda2))), (sin(phi2) - t_0));
} else if (phi2 <= 3200000000.0) {
tmp = atan2(t_1, t_3);
} else {
tmp = atan2(t_2, (sin(phi2) - ((cos(phi2) * (sin((lambda2 + (phi1 - lambda1))) + sin((phi1 + (lambda1 - lambda2))))) / 2.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) :: t_2
real(8) :: t_3
real(8) :: tmp
t_0 = phi1 * cos((lambda2 - lambda1))
t_1 = sin((lambda1 - lambda2))
t_2 = cos(phi2) * t_1
t_3 = (cos(phi1) * sin(phi2)) - ((cos(phi2) * sin(phi1)) * cos((lambda1 - lambda2)))
if (phi2 <= (-9.5d+197)) then
tmp = atan2(t_2, (sin(phi2) - sqrt((t_0 ** 2.0d0))))
else if (phi2 <= (-7d-25)) then
tmp = atan2((sin(lambda1) * cos(phi2)), t_3)
else if (phi2 <= (-1.46d-47)) then
tmp = atan2((cos(phi2) * ((lambda1 * (cos(lambda2) + (lambda1 * (sin(lambda2) * 0.5d0)))) - sin(lambda2))), (sin(phi2) - t_0))
else if (phi2 <= 3200000000.0d0) then
tmp = atan2(t_1, t_3)
else
tmp = atan2(t_2, (sin(phi2) - ((cos(phi2) * (sin((lambda2 + (phi1 - lambda1))) + sin((phi1 + (lambda1 - lambda2))))) / 2.0d0)))
end if
code = tmp
end function
public static double code(double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = phi1 * Math.cos((lambda2 - lambda1));
double t_1 = Math.sin((lambda1 - lambda2));
double t_2 = Math.cos(phi2) * t_1;
double t_3 = (Math.cos(phi1) * Math.sin(phi2)) - ((Math.cos(phi2) * Math.sin(phi1)) * Math.cos((lambda1 - lambda2)));
double tmp;
if (phi2 <= -9.5e+197) {
tmp = Math.atan2(t_2, (Math.sin(phi2) - Math.sqrt(Math.pow(t_0, 2.0))));
} else if (phi2 <= -7e-25) {
tmp = Math.atan2((Math.sin(lambda1) * Math.cos(phi2)), t_3);
} else if (phi2 <= -1.46e-47) {
tmp = Math.atan2((Math.cos(phi2) * ((lambda1 * (Math.cos(lambda2) + (lambda1 * (Math.sin(lambda2) * 0.5)))) - Math.sin(lambda2))), (Math.sin(phi2) - t_0));
} else if (phi2 <= 3200000000.0) {
tmp = Math.atan2(t_1, t_3);
} else {
tmp = Math.atan2(t_2, (Math.sin(phi2) - ((Math.cos(phi2) * (Math.sin((lambda2 + (phi1 - lambda1))) + Math.sin((phi1 + (lambda1 - lambda2))))) / 2.0)));
}
return tmp;
}
def code(lambda1, lambda2, phi1, phi2): t_0 = phi1 * math.cos((lambda2 - lambda1)) t_1 = math.sin((lambda1 - lambda2)) t_2 = math.cos(phi2) * t_1 t_3 = (math.cos(phi1) * math.sin(phi2)) - ((math.cos(phi2) * math.sin(phi1)) * math.cos((lambda1 - lambda2))) tmp = 0 if phi2 <= -9.5e+197: tmp = math.atan2(t_2, (math.sin(phi2) - math.sqrt(math.pow(t_0, 2.0)))) elif phi2 <= -7e-25: tmp = math.atan2((math.sin(lambda1) * math.cos(phi2)), t_3) elif phi2 <= -1.46e-47: tmp = math.atan2((math.cos(phi2) * ((lambda1 * (math.cos(lambda2) + (lambda1 * (math.sin(lambda2) * 0.5)))) - math.sin(lambda2))), (math.sin(phi2) - t_0)) elif phi2 <= 3200000000.0: tmp = math.atan2(t_1, t_3) else: tmp = math.atan2(t_2, (math.sin(phi2) - ((math.cos(phi2) * (math.sin((lambda2 + (phi1 - lambda1))) + math.sin((phi1 + (lambda1 - lambda2))))) / 2.0))) return tmp
function code(lambda1, lambda2, phi1, phi2) t_0 = Float64(phi1 * cos(Float64(lambda2 - lambda1))) t_1 = sin(Float64(lambda1 - lambda2)) t_2 = Float64(cos(phi2) * t_1) t_3 = Float64(Float64(cos(phi1) * sin(phi2)) - Float64(Float64(cos(phi2) * sin(phi1)) * cos(Float64(lambda1 - lambda2)))) tmp = 0.0 if (phi2 <= -9.5e+197) tmp = atan(t_2, Float64(sin(phi2) - sqrt((t_0 ^ 2.0)))); elseif (phi2 <= -7e-25) tmp = atan(Float64(sin(lambda1) * cos(phi2)), t_3); elseif (phi2 <= -1.46e-47) tmp = atan(Float64(cos(phi2) * Float64(Float64(lambda1 * Float64(cos(lambda2) + Float64(lambda1 * Float64(sin(lambda2) * 0.5)))) - sin(lambda2))), Float64(sin(phi2) - t_0)); elseif (phi2 <= 3200000000.0) tmp = atan(t_1, t_3); else tmp = atan(t_2, Float64(sin(phi2) - Float64(Float64(cos(phi2) * Float64(sin(Float64(lambda2 + Float64(phi1 - lambda1))) + sin(Float64(phi1 + Float64(lambda1 - lambda2))))) / 2.0))); end return tmp end
function tmp_2 = code(lambda1, lambda2, phi1, phi2) t_0 = phi1 * cos((lambda2 - lambda1)); t_1 = sin((lambda1 - lambda2)); t_2 = cos(phi2) * t_1; t_3 = (cos(phi1) * sin(phi2)) - ((cos(phi2) * sin(phi1)) * cos((lambda1 - lambda2))); tmp = 0.0; if (phi2 <= -9.5e+197) tmp = atan2(t_2, (sin(phi2) - sqrt((t_0 ^ 2.0)))); elseif (phi2 <= -7e-25) tmp = atan2((sin(lambda1) * cos(phi2)), t_3); elseif (phi2 <= -1.46e-47) tmp = atan2((cos(phi2) * ((lambda1 * (cos(lambda2) + (lambda1 * (sin(lambda2) * 0.5)))) - sin(lambda2))), (sin(phi2) - t_0)); elseif (phi2 <= 3200000000.0) tmp = atan2(t_1, t_3); else tmp = atan2(t_2, (sin(phi2) - ((cos(phi2) * (sin((lambda2 + (phi1 - lambda1))) + sin((phi1 + (lambda1 - lambda2))))) / 2.0))); end tmp_2 = tmp; end
code[lambda1_, lambda2_, phi1_, phi2_] := Block[{t$95$0 = N[(phi1 * N[Cos[N[(lambda2 - lambda1), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$2 = N[(N[Cos[phi2], $MachinePrecision] * t$95$1), $MachinePrecision]}, Block[{t$95$3 = 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]}, If[LessEqual[phi2, -9.5e+197], N[ArcTan[t$95$2 / N[(N[Sin[phi2], $MachinePrecision] - N[Sqrt[N[Power[t$95$0, 2.0], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]], $MachinePrecision], If[LessEqual[phi2, -7e-25], N[ArcTan[N[(N[Sin[lambda1], $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision] / t$95$3], $MachinePrecision], If[LessEqual[phi2, -1.46e-47], N[ArcTan[N[(N[Cos[phi2], $MachinePrecision] * N[(N[(lambda1 * N[(N[Cos[lambda2], $MachinePrecision] + N[(lambda1 * N[(N[Sin[lambda2], $MachinePrecision] * 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[Sin[lambda2], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(N[Sin[phi2], $MachinePrecision] - t$95$0), $MachinePrecision]], $MachinePrecision], If[LessEqual[phi2, 3200000000.0], N[ArcTan[t$95$1 / t$95$3], $MachinePrecision], N[ArcTan[t$95$2 / N[(N[Sin[phi2], $MachinePrecision] - N[(N[(N[Cos[phi2], $MachinePrecision] * N[(N[Sin[N[(lambda2 + N[(phi1 - lambda1), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] + N[Sin[N[(phi1 + N[(lambda1 - lambda2), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / 2.0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \phi_1 \cdot \cos \left(\lambda_2 - \lambda_1\right)\\
t_1 := \sin \left(\lambda_1 - \lambda_2\right)\\
t_2 := \cos \phi_2 \cdot t\_1\\
t_3 := \cos \phi_1 \cdot \sin \phi_2 - \left(\cos \phi_2 \cdot \sin \phi_1\right) \cdot \cos \left(\lambda_1 - \lambda_2\right)\\
\mathbf{if}\;\phi_2 \leq -9.5 \cdot 10^{+197}:\\
\;\;\;\;\tan^{-1}_* \frac{t\_2}{\sin \phi_2 - \sqrt{{t\_0}^{2}}}\\
\mathbf{elif}\;\phi_2 \leq -7 \cdot 10^{-25}:\\
\;\;\;\;\tan^{-1}_* \frac{\sin \lambda_1 \cdot \cos \phi_2}{t\_3}\\
\mathbf{elif}\;\phi_2 \leq -1.46 \cdot 10^{-47}:\\
\;\;\;\;\tan^{-1}_* \frac{\cos \phi_2 \cdot \left(\lambda_1 \cdot \left(\cos \lambda_2 + \lambda_1 \cdot \left(\sin \lambda_2 \cdot 0.5\right)\right) - \sin \lambda_2\right)}{\sin \phi_2 - t\_0}\\
\mathbf{elif}\;\phi_2 \leq 3200000000:\\
\;\;\;\;\tan^{-1}_* \frac{t\_1}{t\_3}\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1}_* \frac{t\_2}{\sin \phi_2 - \frac{\cos \phi_2 \cdot \left(\sin \left(\lambda_2 + \left(\phi_1 - \lambda_1\right)\right) + \sin \left(\phi_1 + \left(\lambda_1 - \lambda_2\right)\right)\right)}{2}}\\
\end{array}
\end{array}
if phi2 < -9.4999999999999997e197Initial program 84.3%
*-commutative84.3%
associate-*l*84.3%
Simplified84.3%
Taylor expanded in phi1 around 0 59.3%
sub-neg59.3%
remove-double-neg59.3%
mul-1-neg59.3%
distribute-neg-in59.3%
+-commutative59.3%
*-commutative59.3%
cos-neg59.3%
mul-1-neg59.3%
unsub-neg59.3%
Simplified59.3%
Taylor expanded in phi1 around 0 59.3%
Taylor expanded in phi2 around 0 54.0%
add-sqr-sqrt31.3%
sqrt-unprod60.9%
pow260.9%
Applied egg-rr60.9%
if -9.4999999999999997e197 < phi2 < -7.0000000000000004e-25Initial program 78.5%
Taylor expanded in lambda2 around 0 51.7%
if -7.0000000000000004e-25 < phi2 < -1.46e-47Initial program 57.8%
*-commutative57.8%
associate-*l*57.8%
Simplified57.8%
Taylor expanded in phi1 around 0 57.8%
sub-neg57.8%
remove-double-neg57.8%
mul-1-neg57.8%
distribute-neg-in57.8%
+-commutative57.8%
*-commutative57.8%
cos-neg57.8%
mul-1-neg57.8%
unsub-neg57.8%
Simplified57.8%
Taylor expanded in phi1 around 0 57.8%
Taylor expanded in phi2 around 0 57.8%
Taylor expanded in lambda1 around 0 88.6%
+-commutative88.6%
sin-neg88.6%
unsub-neg88.6%
*-commutative88.6%
associate-*r*88.6%
*-commutative88.6%
cos-neg88.6%
*-commutative88.6%
sin-neg88.6%
distribute-lft-neg-out88.6%
distribute-rgt-neg-in88.6%
metadata-eval88.6%
Simplified88.6%
if -1.46e-47 < phi2 < 3.2e9Initial program 85.8%
add-sqr-sqrt44.6%
pow244.6%
Applied egg-rr44.6%
Taylor expanded in phi2 around 0 84.9%
associate-*r*84.9%
distribute-rgt1-in84.9%
Simplified84.9%
Taylor expanded in phi2 around 0 85.1%
if 3.2e9 < phi2 Initial program 77.1%
*-commutative77.1%
associate-*l*77.1%
Simplified77.1%
*-commutative77.1%
sin-cos-mult57.3%
associate-*l/57.3%
associate--r-57.3%
+-commutative57.3%
Applied egg-rr57.3%
Taylor expanded in phi1 around 0 54.9%
Final simplification69.6%
(FPCore (lambda1 lambda2 phi1 phi2) :precision binary64 (atan2 (* (cos phi2) (sin (- lambda1 lambda2))) (- (* (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) * sin((lambda1 - lambda2))), ((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) * sin((lambda1 - lambda2))), ((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.sin((lambda1 - lambda2))), ((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.sin((lambda1 - lambda2))), ((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) * sin(Float64(lambda1 - lambda2))), 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) * sin((lambda1 - lambda2))), ((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[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]
\begin{array}{l}
\\
\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)}
\end{array}
Initial program 81.3%
*-commutative81.3%
associate-*l*81.3%
Simplified81.3%
Final simplification81.3%
(FPCore (lambda1 lambda2 phi1 phi2) :precision binary64 (atan2 (* (cos phi2) (sin (- lambda1 lambda2))) (- (* (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) * sin((lambda1 - lambda2))), ((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) * sin((lambda1 - lambda2))), ((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.sin((lambda1 - lambda2))), ((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.sin((lambda1 - lambda2))), ((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) * sin(Float64(lambda1 - lambda2))), Float64(Float64(cos(phi1) * sin(phi2)) - Float64(Float64(cos(phi2) * sin(phi1)) * cos(Float64(lambda1 - lambda2))))) end
function tmp = code(lambda1, lambda2, phi1, phi2) tmp = atan2((cos(phi2) * sin((lambda1 - lambda2))), ((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[Sin[N[(lambda1 - lambda2), $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[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
\\
\tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\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 81.3%
Final simplification81.3%
(FPCore (lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (sin (- lambda1 lambda2))) (t_1 (* (cos phi2) t_0)))
(if (<= phi2 -8100.0)
(atan2
t_1
(- (sin phi2) (sqrt (pow (* phi1 (cos (- lambda2 lambda1))) 2.0))))
(if (<= phi2 3200000000.0)
(atan2
t_0
(-
(* (cos phi1) (sin phi2))
(* (* (cos phi2) (sin phi1)) (cos (- lambda1 lambda2)))))
(atan2
t_1
(-
(sin phi2)
(/
(*
(cos phi2)
(+
(sin (+ lambda2 (- phi1 lambda1)))
(sin (+ phi1 (- lambda1 lambda2)))))
2.0)))))))
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 <= -8100.0) {
tmp = atan2(t_1, (sin(phi2) - sqrt(pow((phi1 * cos((lambda2 - lambda1))), 2.0))));
} else if (phi2 <= 3200000000.0) {
tmp = atan2(t_0, ((cos(phi1) * sin(phi2)) - ((cos(phi2) * sin(phi1)) * cos((lambda1 - lambda2)))));
} else {
tmp = atan2(t_1, (sin(phi2) - ((cos(phi2) * (sin((lambda2 + (phi1 - lambda1))) + sin((phi1 + (lambda1 - lambda2))))) / 2.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((lambda1 - lambda2))
t_1 = cos(phi2) * t_0
if (phi2 <= (-8100.0d0)) then
tmp = atan2(t_1, (sin(phi2) - sqrt(((phi1 * cos((lambda2 - lambda1))) ** 2.0d0))))
else if (phi2 <= 3200000000.0d0) then
tmp = atan2(t_0, ((cos(phi1) * sin(phi2)) - ((cos(phi2) * sin(phi1)) * cos((lambda1 - lambda2)))))
else
tmp = atan2(t_1, (sin(phi2) - ((cos(phi2) * (sin((lambda2 + (phi1 - lambda1))) + sin((phi1 + (lambda1 - lambda2))))) / 2.0d0)))
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 t_1 = Math.cos(phi2) * t_0;
double tmp;
if (phi2 <= -8100.0) {
tmp = Math.atan2(t_1, (Math.sin(phi2) - Math.sqrt(Math.pow((phi1 * Math.cos((lambda2 - lambda1))), 2.0))));
} else if (phi2 <= 3200000000.0) {
tmp = Math.atan2(t_0, ((Math.cos(phi1) * Math.sin(phi2)) - ((Math.cos(phi2) * Math.sin(phi1)) * Math.cos((lambda1 - lambda2)))));
} else {
tmp = Math.atan2(t_1, (Math.sin(phi2) - ((Math.cos(phi2) * (Math.sin((lambda2 + (phi1 - lambda1))) + Math.sin((phi1 + (lambda1 - lambda2))))) / 2.0)));
}
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 <= -8100.0: tmp = math.atan2(t_1, (math.sin(phi2) - math.sqrt(math.pow((phi1 * math.cos((lambda2 - lambda1))), 2.0)))) elif phi2 <= 3200000000.0: tmp = math.atan2(t_0, ((math.cos(phi1) * math.sin(phi2)) - ((math.cos(phi2) * math.sin(phi1)) * math.cos((lambda1 - lambda2))))) else: tmp = math.atan2(t_1, (math.sin(phi2) - ((math.cos(phi2) * (math.sin((lambda2 + (phi1 - lambda1))) + math.sin((phi1 + (lambda1 - lambda2))))) / 2.0))) 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 <= -8100.0) tmp = atan(t_1, Float64(sin(phi2) - sqrt((Float64(phi1 * cos(Float64(lambda2 - lambda1))) ^ 2.0)))); elseif (phi2 <= 3200000000.0) tmp = atan(t_0, Float64(Float64(cos(phi1) * sin(phi2)) - Float64(Float64(cos(phi2) * sin(phi1)) * cos(Float64(lambda1 - lambda2))))); else tmp = atan(t_1, Float64(sin(phi2) - Float64(Float64(cos(phi2) * Float64(sin(Float64(lambda2 + Float64(phi1 - lambda1))) + sin(Float64(phi1 + Float64(lambda1 - lambda2))))) / 2.0))); end return tmp end
function tmp_2 = code(lambda1, lambda2, phi1, phi2) t_0 = sin((lambda1 - lambda2)); t_1 = cos(phi2) * t_0; tmp = 0.0; if (phi2 <= -8100.0) tmp = atan2(t_1, (sin(phi2) - sqrt(((phi1 * cos((lambda2 - lambda1))) ^ 2.0)))); elseif (phi2 <= 3200000000.0) tmp = atan2(t_0, ((cos(phi1) * sin(phi2)) - ((cos(phi2) * sin(phi1)) * cos((lambda1 - lambda2))))); else tmp = atan2(t_1, (sin(phi2) - ((cos(phi2) * (sin((lambda2 + (phi1 - lambda1))) + sin((phi1 + (lambda1 - lambda2))))) / 2.0))); end tmp_2 = 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, -8100.0], N[ArcTan[t$95$1 / N[(N[Sin[phi2], $MachinePrecision] - N[Sqrt[N[Power[N[(phi1 * N[Cos[N[(lambda2 - lambda1), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]], $MachinePrecision], If[LessEqual[phi2, 3200000000.0], N[ArcTan[t$95$0 / 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], N[ArcTan[t$95$1 / N[(N[Sin[phi2], $MachinePrecision] - N[(N[(N[Cos[phi2], $MachinePrecision] * N[(N[Sin[N[(lambda2 + N[(phi1 - lambda1), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] + N[Sin[N[(phi1 + N[(lambda1 - lambda2), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / 2.0), $MachinePrecision]), $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 -8100:\\
\;\;\;\;\tan^{-1}_* \frac{t\_1}{\sin \phi_2 - \sqrt{{\left(\phi_1 \cdot \cos \left(\lambda_2 - \lambda_1\right)\right)}^{2}}}\\
\mathbf{elif}\;\phi_2 \leq 3200000000:\\
\;\;\;\;\tan^{-1}_* \frac{t\_0}{\cos \phi_1 \cdot \sin \phi_2 - \left(\cos \phi_2 \cdot \sin \phi_1\right) \cdot \cos \left(\lambda_1 - \lambda_2\right)}\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1}_* \frac{t\_1}{\sin \phi_2 - \frac{\cos \phi_2 \cdot \left(\sin \left(\lambda_2 + \left(\phi_1 - \lambda_1\right)\right) + \sin \left(\phi_1 + \left(\lambda_1 - \lambda_2\right)\right)\right)}{2}}\\
\end{array}
\end{array}
if phi2 < -8100Initial program 82.0%
*-commutative82.0%
associate-*l*82.0%
Simplified82.0%
Taylor expanded in phi1 around 0 47.2%
sub-neg47.2%
remove-double-neg47.2%
mul-1-neg47.2%
distribute-neg-in47.2%
+-commutative47.2%
*-commutative47.2%
cos-neg47.2%
mul-1-neg47.2%
unsub-neg47.2%
Simplified47.2%
Taylor expanded in phi1 around 0 47.0%
Taylor expanded in phi2 around 0 41.8%
add-sqr-sqrt25.2%
sqrt-unprod50.3%
pow250.3%
Applied egg-rr50.3%
if -8100 < phi2 < 3.2e9Initial program 83.1%
add-sqr-sqrt43.5%
pow243.5%
Applied egg-rr43.5%
Taylor expanded in phi2 around 0 81.7%
associate-*r*81.7%
distribute-rgt1-in81.7%
Simplified81.7%
Taylor expanded in phi2 around 0 81.6%
if 3.2e9 < phi2 Initial program 77.1%
*-commutative77.1%
associate-*l*77.1%
Simplified77.1%
*-commutative77.1%
sin-cos-mult57.3%
associate-*l/57.3%
associate--r-57.3%
+-commutative57.3%
Applied egg-rr57.3%
Taylor expanded in phi1 around 0 54.9%
Final simplification67.9%
(FPCore (lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (sin (- lambda1 lambda2))))
(if (or (<= phi2 -8100.0) (not (<= phi2 3200000000.0)))
(atan2
(* (cos phi2) t_0)
(- (sin phi2) (sqrt (pow (* phi1 (cos (- lambda2 lambda1))) 2.0))))
(atan2
t_0
(-
(* (cos phi1) (sin phi2))
(* (* (cos 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 <= -8100.0) || !(phi2 <= 3200000000.0)) {
tmp = atan2((cos(phi2) * t_0), (sin(phi2) - sqrt(pow((phi1 * cos((lambda2 - lambda1))), 2.0))));
} else {
tmp = atan2(t_0, ((cos(phi1) * sin(phi2)) - ((cos(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 <= (-8100.0d0)) .or. (.not. (phi2 <= 3200000000.0d0))) then
tmp = atan2((cos(phi2) * t_0), (sin(phi2) - sqrt(((phi1 * cos((lambda2 - lambda1))) ** 2.0d0))))
else
tmp = atan2(t_0, ((cos(phi1) * sin(phi2)) - ((cos(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 <= -8100.0) || !(phi2 <= 3200000000.0)) {
tmp = Math.atan2((Math.cos(phi2) * t_0), (Math.sin(phi2) - Math.sqrt(Math.pow((phi1 * Math.cos((lambda2 - lambda1))), 2.0))));
} else {
tmp = Math.atan2(t_0, ((Math.cos(phi1) * Math.sin(phi2)) - ((Math.cos(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 <= -8100.0) or not (phi2 <= 3200000000.0): tmp = math.atan2((math.cos(phi2) * t_0), (math.sin(phi2) - math.sqrt(math.pow((phi1 * math.cos((lambda2 - lambda1))), 2.0)))) else: tmp = math.atan2(t_0, ((math.cos(phi1) * math.sin(phi2)) - ((math.cos(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 <= -8100.0) || !(phi2 <= 3200000000.0)) tmp = atan(Float64(cos(phi2) * t_0), Float64(sin(phi2) - sqrt((Float64(phi1 * cos(Float64(lambda2 - lambda1))) ^ 2.0)))); else tmp = atan(t_0, Float64(Float64(cos(phi1) * sin(phi2)) - Float64(Float64(cos(phi2) * 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 <= -8100.0) || ~((phi2 <= 3200000000.0))) tmp = atan2((cos(phi2) * t_0), (sin(phi2) - sqrt(((phi1 * cos((lambda2 - lambda1))) ^ 2.0)))); else tmp = atan2(t_0, ((cos(phi1) * sin(phi2)) - ((cos(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, -8100.0], N[Not[LessEqual[phi2, 3200000000.0]], $MachinePrecision]], N[ArcTan[N[(N[Cos[phi2], $MachinePrecision] * t$95$0), $MachinePrecision] / N[(N[Sin[phi2], $MachinePrecision] - N[Sqrt[N[Power[N[(phi1 * N[Cos[N[(lambda2 - lambda1), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]], $MachinePrecision], N[ArcTan[t$95$0 / 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}
\\
\begin{array}{l}
t_0 := \sin \left(\lambda_1 - \lambda_2\right)\\
\mathbf{if}\;\phi_2 \leq -8100 \lor \neg \left(\phi_2 \leq 3200000000\right):\\
\;\;\;\;\tan^{-1}_* \frac{\cos \phi_2 \cdot t\_0}{\sin \phi_2 - \sqrt{{\left(\phi_1 \cdot \cos \left(\lambda_2 - \lambda_1\right)\right)}^{2}}}\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1}_* \frac{t\_0}{\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}
\end{array}
if phi2 < -8100 or 3.2e9 < phi2 Initial program 79.3%
*-commutative79.3%
associate-*l*79.3%
Simplified79.3%
Taylor expanded in phi1 around 0 46.8%
sub-neg46.8%
remove-double-neg46.8%
mul-1-neg46.8%
distribute-neg-in46.8%
+-commutative46.8%
*-commutative46.8%
cos-neg46.8%
mul-1-neg46.8%
unsub-neg46.8%
Simplified46.8%
Taylor expanded in phi1 around 0 46.8%
Taylor expanded in phi2 around 0 45.5%
add-sqr-sqrt28.5%
sqrt-unprod52.6%
pow252.6%
Applied egg-rr52.6%
if -8100 < phi2 < 3.2e9Initial program 83.1%
add-sqr-sqrt43.5%
pow243.5%
Applied egg-rr43.5%
Taylor expanded in phi2 around 0 81.7%
associate-*r*81.7%
distribute-rgt1-in81.7%
Simplified81.7%
Taylor expanded in phi2 around 0 81.6%
Final simplification67.8%
(FPCore (lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (* (cos phi2) (sin (- lambda1 lambda2)))))
(if (or (<= phi1 -1e+24) (not (<= phi1 4.3e+134)))
(atan2 (pow (sqrt t_0) 2.0) (* (sin phi1) (- (cos (- lambda1 lambda2)))))
(atan2 t_0 (- (sin phi2) (* (cos lambda2) phi1))))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = cos(phi2) * sin((lambda1 - lambda2));
double tmp;
if ((phi1 <= -1e+24) || !(phi1 <= 4.3e+134)) {
tmp = atan2(pow(sqrt(t_0), 2.0), (sin(phi1) * -cos((lambda1 - lambda2))));
} else {
tmp = atan2(t_0, (sin(phi2) - (cos(lambda2) * 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) :: tmp
t_0 = cos(phi2) * sin((lambda1 - lambda2))
if ((phi1 <= (-1d+24)) .or. (.not. (phi1 <= 4.3d+134))) then
tmp = atan2((sqrt(t_0) ** 2.0d0), (sin(phi1) * -cos((lambda1 - lambda2))))
else
tmp = atan2(t_0, (sin(phi2) - (cos(lambda2) * 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 tmp;
if ((phi1 <= -1e+24) || !(phi1 <= 4.3e+134)) {
tmp = Math.atan2(Math.pow(Math.sqrt(t_0), 2.0), (Math.sin(phi1) * -Math.cos((lambda1 - lambda2))));
} else {
tmp = Math.atan2(t_0, (Math.sin(phi2) - (Math.cos(lambda2) * phi1)));
}
return tmp;
}
def code(lambda1, lambda2, phi1, phi2): t_0 = math.cos(phi2) * math.sin((lambda1 - lambda2)) tmp = 0 if (phi1 <= -1e+24) or not (phi1 <= 4.3e+134): tmp = math.atan2(math.pow(math.sqrt(t_0), 2.0), (math.sin(phi1) * -math.cos((lambda1 - lambda2)))) else: tmp = math.atan2(t_0, (math.sin(phi2) - (math.cos(lambda2) * phi1))) return tmp
function code(lambda1, lambda2, phi1, phi2) t_0 = Float64(cos(phi2) * sin(Float64(lambda1 - lambda2))) tmp = 0.0 if ((phi1 <= -1e+24) || !(phi1 <= 4.3e+134)) tmp = atan((sqrt(t_0) ^ 2.0), Float64(sin(phi1) * Float64(-cos(Float64(lambda1 - lambda2))))); else tmp = atan(t_0, Float64(sin(phi2) - Float64(cos(lambda2) * phi1))); end return tmp end
function tmp_2 = code(lambda1, lambda2, phi1, phi2) t_0 = cos(phi2) * sin((lambda1 - lambda2)); tmp = 0.0; if ((phi1 <= -1e+24) || ~((phi1 <= 4.3e+134))) tmp = atan2((sqrt(t_0) ^ 2.0), (sin(phi1) * -cos((lambda1 - lambda2)))); else tmp = atan2(t_0, (sin(phi2) - (cos(lambda2) * 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]}, If[Or[LessEqual[phi1, -1e+24], N[Not[LessEqual[phi1, 4.3e+134]], $MachinePrecision]], N[ArcTan[N[Power[N[Sqrt[t$95$0], $MachinePrecision], 2.0], $MachinePrecision] / N[(N[Sin[phi1], $MachinePrecision] * (-N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision])), $MachinePrecision]], $MachinePrecision], N[ArcTan[t$95$0 / N[(N[Sin[phi2], $MachinePrecision] - N[(N[Cos[lambda2], $MachinePrecision] * phi1), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)\\
\mathbf{if}\;\phi_1 \leq -1 \cdot 10^{+24} \lor \neg \left(\phi_1 \leq 4.3 \cdot 10^{+134}\right):\\
\;\;\;\;\tan^{-1}_* \frac{{\left(\sqrt{t\_0}\right)}^{2}}{\sin \phi_1 \cdot \left(-\cos \left(\lambda_1 - \lambda_2\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1}_* \frac{t\_0}{\sin \phi_2 - \cos \lambda_2 \cdot \phi_1}\\
\end{array}
\end{array}
if phi1 < -9.9999999999999998e23 or 4.3e134 < phi1 Initial program 78.6%
add-sqr-sqrt43.2%
pow243.2%
Applied egg-rr43.2%
Taylor expanded in phi2 around 0 28.5%
associate-*r*28.5%
neg-mul-128.5%
Simplified28.5%
if -9.9999999999999998e23 < phi1 < 4.3e134Initial program 82.9%
*-commutative82.9%
associate-*l*82.9%
Simplified82.9%
Taylor expanded in phi1 around 0 67.5%
sub-neg67.5%
remove-double-neg67.5%
mul-1-neg67.5%
distribute-neg-in67.5%
+-commutative67.5%
*-commutative67.5%
cos-neg67.5%
mul-1-neg67.5%
unsub-neg67.5%
Simplified67.5%
Taylor expanded in phi1 around 0 67.5%
Taylor expanded in phi2 around 0 68.0%
Taylor expanded in lambda1 around 0 68.2%
Final simplification53.3%
(FPCore (lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (* (cos phi2) (sin (- lambda1 lambda2)))))
(if (<= phi1 4.6e+134)
(atan2
t_0
(- (sin phi2) (sqrt (pow (* phi1 (cos (- lambda2 lambda1))) 2.0))))
(atan2
(pow (sqrt t_0) 2.0)
(* (sin phi1) (- (cos (- lambda1 lambda2))))))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = cos(phi2) * sin((lambda1 - lambda2));
double tmp;
if (phi1 <= 4.6e+134) {
tmp = atan2(t_0, (sin(phi2) - sqrt(pow((phi1 * cos((lambda2 - lambda1))), 2.0))));
} else {
tmp = atan2(pow(sqrt(t_0), 2.0), (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 = cos(phi2) * sin((lambda1 - lambda2))
if (phi1 <= 4.6d+134) then
tmp = atan2(t_0, (sin(phi2) - sqrt(((phi1 * cos((lambda2 - lambda1))) ** 2.0d0))))
else
tmp = atan2((sqrt(t_0) ** 2.0d0), (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.cos(phi2) * Math.sin((lambda1 - lambda2));
double tmp;
if (phi1 <= 4.6e+134) {
tmp = Math.atan2(t_0, (Math.sin(phi2) - Math.sqrt(Math.pow((phi1 * Math.cos((lambda2 - lambda1))), 2.0))));
} else {
tmp = Math.atan2(Math.pow(Math.sqrt(t_0), 2.0), (Math.sin(phi1) * -Math.cos((lambda1 - lambda2))));
}
return tmp;
}
def code(lambda1, lambda2, phi1, phi2): t_0 = math.cos(phi2) * math.sin((lambda1 - lambda2)) tmp = 0 if phi1 <= 4.6e+134: tmp = math.atan2(t_0, (math.sin(phi2) - math.sqrt(math.pow((phi1 * math.cos((lambda2 - lambda1))), 2.0)))) else: tmp = math.atan2(math.pow(math.sqrt(t_0), 2.0), (math.sin(phi1) * -math.cos((lambda1 - lambda2)))) return tmp
function code(lambda1, lambda2, phi1, phi2) t_0 = Float64(cos(phi2) * sin(Float64(lambda1 - lambda2))) tmp = 0.0 if (phi1 <= 4.6e+134) tmp = atan(t_0, Float64(sin(phi2) - sqrt((Float64(phi1 * cos(Float64(lambda2 - lambda1))) ^ 2.0)))); else tmp = atan((sqrt(t_0) ^ 2.0), Float64(sin(phi1) * Float64(-cos(Float64(lambda1 - lambda2))))); end return tmp end
function tmp_2 = code(lambda1, lambda2, phi1, phi2) t_0 = cos(phi2) * sin((lambda1 - lambda2)); tmp = 0.0; if (phi1 <= 4.6e+134) tmp = atan2(t_0, (sin(phi2) - sqrt(((phi1 * cos((lambda2 - lambda1))) ^ 2.0)))); else tmp = atan2((sqrt(t_0) ^ 2.0), (sin(phi1) * -cos((lambda1 - lambda2)))); 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]}, If[LessEqual[phi1, 4.6e+134], N[ArcTan[t$95$0 / N[(N[Sin[phi2], $MachinePrecision] - N[Sqrt[N[Power[N[(phi1 * N[Cos[N[(lambda2 - lambda1), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]], $MachinePrecision], N[ArcTan[N[Power[N[Sqrt[t$95$0], $MachinePrecision], 2.0], $MachinePrecision] / N[(N[Sin[phi1], $MachinePrecision] * (-N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision])), $MachinePrecision]], $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)\\
\mathbf{if}\;\phi_1 \leq 4.6 \cdot 10^{+134}:\\
\;\;\;\;\tan^{-1}_* \frac{t\_0}{\sin \phi_2 - \sqrt{{\left(\phi_1 \cdot \cos \left(\lambda_2 - \lambda_1\right)\right)}^{2}}}\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1}_* \frac{{\left(\sqrt{t\_0}\right)}^{2}}{\sin \phi_1 \cdot \left(-\cos \left(\lambda_1 - \lambda_2\right)\right)}\\
\end{array}
\end{array}
if phi1 < 4.5999999999999996e134Initial program 82.1%
*-commutative82.1%
associate-*l*82.1%
Simplified82.1%
Taylor expanded in phi1 around 0 54.3%
sub-neg54.3%
remove-double-neg54.3%
mul-1-neg54.3%
distribute-neg-in54.3%
+-commutative54.3%
*-commutative54.3%
cos-neg54.3%
mul-1-neg54.3%
unsub-neg54.3%
Simplified54.3%
Taylor expanded in phi1 around 0 54.3%
Taylor expanded in phi2 around 0 53.4%
add-sqr-sqrt32.5%
sqrt-unprod57.3%
pow257.3%
Applied egg-rr57.3%
if 4.5999999999999996e134 < phi1 Initial program 77.2%
add-sqr-sqrt46.4%
pow246.4%
Applied egg-rr46.4%
Taylor expanded in phi2 around 0 31.4%
associate-*r*31.4%
neg-mul-131.4%
Simplified31.4%
Final simplification53.0%
(FPCore (lambda1 lambda2 phi1 phi2) :precision binary64 (atan2 (* (cos phi2) (sin (- lambda1 lambda2))) (- (sin phi2) (* phi1 (* (cos lambda1) (cos phi2))))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
return atan2((cos(phi2) * sin((lambda1 - lambda2))), (sin(phi2) - (phi1 * (cos(lambda1) * cos(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) - (phi1 * (cos(lambda1) * cos(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) - (phi1 * (Math.cos(lambda1) * Math.cos(phi2)))));
}
def code(lambda1, lambda2, phi1, phi2): return math.atan2((math.cos(phi2) * math.sin((lambda1 - lambda2))), (math.sin(phi2) - (phi1 * (math.cos(lambda1) * math.cos(phi2)))))
function code(lambda1, lambda2, phi1, phi2) return atan(Float64(cos(phi2) * sin(Float64(lambda1 - lambda2))), Float64(sin(phi2) - Float64(phi1 * Float64(cos(lambda1) * cos(phi2))))) end
function tmp = code(lambda1, lambda2, phi1, phi2) tmp = atan2((cos(phi2) * sin((lambda1 - lambda2))), (sin(phi2) - (phi1 * (cos(lambda1) * cos(phi2))))); end
code[lambda1_, lambda2_, phi1_, phi2_] := N[ArcTan[N[(N[Cos[phi2], $MachinePrecision] * N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / N[(N[Sin[phi2], $MachinePrecision] - N[(phi1 * N[(N[Cos[lambda1], $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
\\
\tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\sin \phi_2 - \phi_1 \cdot \left(\cos \lambda_1 \cdot \cos \phi_2\right)}
\end{array}
Initial program 81.3%
*-commutative81.3%
associate-*l*81.3%
Simplified81.3%
Taylor expanded in phi1 around 0 47.1%
sub-neg47.1%
remove-double-neg47.1%
mul-1-neg47.1%
distribute-neg-in47.1%
+-commutative47.1%
*-commutative47.1%
cos-neg47.1%
mul-1-neg47.1%
unsub-neg47.1%
Simplified47.1%
Taylor expanded in phi1 around 0 47.1%
Taylor expanded in lambda2 around 0 46.7%
cos-neg46.7%
Simplified46.7%
Final simplification46.7%
(FPCore (lambda1 lambda2 phi1 phi2) :precision binary64 (atan2 (* (cos phi2) (sin (- lambda1 lambda2))) (- (sin phi2) (* (cos phi2) (* (cos lambda2) phi1)))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
return atan2((cos(phi2) * sin((lambda1 - lambda2))), (sin(phi2) - (cos(phi2) * (cos(lambda2) * phi1))));
}
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) - (cos(phi2) * (cos(lambda2) * phi1))))
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) - (Math.cos(phi2) * (Math.cos(lambda2) * phi1))));
}
def code(lambda1, lambda2, phi1, phi2): return math.atan2((math.cos(phi2) * math.sin((lambda1 - lambda2))), (math.sin(phi2) - (math.cos(phi2) * (math.cos(lambda2) * phi1))))
function code(lambda1, lambda2, phi1, phi2) return atan(Float64(cos(phi2) * sin(Float64(lambda1 - lambda2))), Float64(sin(phi2) - Float64(cos(phi2) * Float64(cos(lambda2) * phi1)))) end
function tmp = code(lambda1, lambda2, phi1, phi2) tmp = atan2((cos(phi2) * sin((lambda1 - lambda2))), (sin(phi2) - (cos(phi2) * (cos(lambda2) * phi1)))); end
code[lambda1_, lambda2_, phi1_, phi2_] := 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[Cos[lambda2], $MachinePrecision] * phi1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
\\
\tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\sin \phi_2 - \cos \phi_2 \cdot \left(\cos \lambda_2 \cdot \phi_1\right)}
\end{array}
Initial program 81.3%
*-commutative81.3%
associate-*l*81.3%
Simplified81.3%
Taylor expanded in phi1 around 0 47.1%
sub-neg47.1%
remove-double-neg47.1%
mul-1-neg47.1%
distribute-neg-in47.1%
+-commutative47.1%
*-commutative47.1%
cos-neg47.1%
mul-1-neg47.1%
unsub-neg47.1%
Simplified47.1%
Taylor expanded in phi1 around 0 47.1%
Taylor expanded in lambda1 around 0 47.8%
associate-*r*47.8%
Simplified47.8%
Final simplification47.8%
(FPCore (lambda1 lambda2 phi1 phi2) :precision binary64 (atan2 (* (cos phi2) (sin (- lambda1 lambda2))) (- (sin phi2) (* phi1 (cos (- lambda2 lambda1))))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
return atan2((cos(phi2) * sin((lambda1 - lambda2))), (sin(phi2) - (phi1 * cos((lambda2 - lambda1)))));
}
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) - (phi1 * cos((lambda2 - lambda1)))))
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) - (phi1 * Math.cos((lambda2 - lambda1)))));
}
def code(lambda1, lambda2, phi1, phi2): return math.atan2((math.cos(phi2) * math.sin((lambda1 - lambda2))), (math.sin(phi2) - (phi1 * math.cos((lambda2 - lambda1)))))
function code(lambda1, lambda2, phi1, phi2) return atan(Float64(cos(phi2) * sin(Float64(lambda1 - lambda2))), Float64(sin(phi2) - Float64(phi1 * cos(Float64(lambda2 - lambda1))))) end
function tmp = code(lambda1, lambda2, phi1, phi2) tmp = atan2((cos(phi2) * sin((lambda1 - lambda2))), (sin(phi2) - (phi1 * cos((lambda2 - lambda1))))); end
code[lambda1_, lambda2_, phi1_, phi2_] := N[ArcTan[N[(N[Cos[phi2], $MachinePrecision] * N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / N[(N[Sin[phi2], $MachinePrecision] - N[(phi1 * N[Cos[N[(lambda2 - lambda1), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
\\
\tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\sin \phi_2 - \phi_1 \cdot \cos \left(\lambda_2 - \lambda_1\right)}
\end{array}
Initial program 81.3%
*-commutative81.3%
associate-*l*81.3%
Simplified81.3%
Taylor expanded in phi1 around 0 47.1%
sub-neg47.1%
remove-double-neg47.1%
mul-1-neg47.1%
distribute-neg-in47.1%
+-commutative47.1%
*-commutative47.1%
cos-neg47.1%
mul-1-neg47.1%
unsub-neg47.1%
Simplified47.1%
Taylor expanded in phi1 around 0 47.1%
Taylor expanded in phi2 around 0 46.5%
Final simplification46.5%
(FPCore (lambda1 lambda2 phi1 phi2) :precision binary64 (atan2 (* (sin lambda1) (cos phi2)) (- (sin phi2) (* phi1 (cos (- lambda2 lambda1))))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
return atan2((sin(lambda1) * cos(phi2)), (sin(phi2) - (phi1 * cos((lambda2 - lambda1)))));
}
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) * cos(phi2)), (sin(phi2) - (phi1 * cos((lambda2 - lambda1)))))
end function
public static double code(double lambda1, double lambda2, double phi1, double phi2) {
return Math.atan2((Math.sin(lambda1) * Math.cos(phi2)), (Math.sin(phi2) - (phi1 * Math.cos((lambda2 - lambda1)))));
}
def code(lambda1, lambda2, phi1, phi2): return math.atan2((math.sin(lambda1) * math.cos(phi2)), (math.sin(phi2) - (phi1 * math.cos((lambda2 - lambda1)))))
function code(lambda1, lambda2, phi1, phi2) return atan(Float64(sin(lambda1) * cos(phi2)), Float64(sin(phi2) - Float64(phi1 * cos(Float64(lambda2 - lambda1))))) end
function tmp = code(lambda1, lambda2, phi1, phi2) tmp = atan2((sin(lambda1) * cos(phi2)), (sin(phi2) - (phi1 * cos((lambda2 - lambda1))))); end
code[lambda1_, lambda2_, phi1_, phi2_] := N[ArcTan[N[(N[Sin[lambda1], $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision] / N[(N[Sin[phi2], $MachinePrecision] - N[(phi1 * N[Cos[N[(lambda2 - lambda1), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
\\
\tan^{-1}_* \frac{\sin \lambda_1 \cdot \cos \phi_2}{\sin \phi_2 - \phi_1 \cdot \cos \left(\lambda_2 - \lambda_1\right)}
\end{array}
Initial program 81.3%
*-commutative81.3%
associate-*l*81.3%
Simplified81.3%
Taylor expanded in phi1 around 0 47.1%
sub-neg47.1%
remove-double-neg47.1%
mul-1-neg47.1%
distribute-neg-in47.1%
+-commutative47.1%
*-commutative47.1%
cos-neg47.1%
mul-1-neg47.1%
unsub-neg47.1%
Simplified47.1%
Taylor expanded in phi1 around 0 47.1%
Taylor expanded in phi2 around 0 46.5%
Taylor expanded in lambda2 around 0 30.3%
Final simplification30.3%
(FPCore (lambda1 lambda2 phi1 phi2) :precision binary64 (atan2 (* (cos phi2) (sin (- lambda1 lambda2))) (- (sin phi2) (* (cos lambda2) phi1))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
return atan2((cos(phi2) * sin((lambda1 - lambda2))), (sin(phi2) - (cos(lambda2) * phi1)));
}
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) - (cos(lambda2) * phi1)))
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) - (Math.cos(lambda2) * phi1)));
}
def code(lambda1, lambda2, phi1, phi2): return math.atan2((math.cos(phi2) * math.sin((lambda1 - lambda2))), (math.sin(phi2) - (math.cos(lambda2) * phi1)))
function code(lambda1, lambda2, phi1, phi2) return atan(Float64(cos(phi2) * sin(Float64(lambda1 - lambda2))), Float64(sin(phi2) - Float64(cos(lambda2) * phi1))) end
function tmp = code(lambda1, lambda2, phi1, phi2) tmp = atan2((cos(phi2) * sin((lambda1 - lambda2))), (sin(phi2) - (cos(lambda2) * phi1))); end
code[lambda1_, lambda2_, phi1_, phi2_] := N[ArcTan[N[(N[Cos[phi2], $MachinePrecision] * N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / N[(N[Sin[phi2], $MachinePrecision] - N[(N[Cos[lambda2], $MachinePrecision] * phi1), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
\\
\tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\sin \phi_2 - \cos \lambda_2 \cdot \phi_1}
\end{array}
Initial program 81.3%
*-commutative81.3%
associate-*l*81.3%
Simplified81.3%
Taylor expanded in phi1 around 0 47.1%
sub-neg47.1%
remove-double-neg47.1%
mul-1-neg47.1%
distribute-neg-in47.1%
+-commutative47.1%
*-commutative47.1%
cos-neg47.1%
mul-1-neg47.1%
unsub-neg47.1%
Simplified47.1%
Taylor expanded in phi1 around 0 47.1%
Taylor expanded in phi2 around 0 46.5%
Taylor expanded in lambda1 around 0 46.5%
Final simplification46.5%
herbie shell --seed 2024054
(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))))))