
(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 27 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
(let* ((t_0 (* (cos phi1) (sin phi2)))
(t_1 (* (cos phi2) (sin phi1)))
(t_2 (* (sin (- lambda1 lambda2)) (cos phi2)))
(t_3 (- t_0 (* t_1 (cos lambda1)))))
(if (<= lambda1 -12500.0)
(atan2 t_2 t_3)
(if (<= lambda1 2.6e-7)
(atan2 t_2 (- t_0 (* t_1 (cos lambda2))))
(atan2 (* (cos phi2) (sin lambda1)) 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(phi1);
double t_2 = sin((lambda1 - lambda2)) * cos(phi2);
double t_3 = t_0 - (t_1 * cos(lambda1));
double tmp;
if (lambda1 <= -12500.0) {
tmp = atan2(t_2, t_3);
} else if (lambda1 <= 2.6e-7) {
tmp = atan2(t_2, (t_0 - (t_1 * cos(lambda2))));
} else {
tmp = atan2((cos(phi2) * sin(lambda1)), t_3);
}
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 = cos(phi1) * sin(phi2)
t_1 = cos(phi2) * sin(phi1)
t_2 = sin((lambda1 - lambda2)) * cos(phi2)
t_3 = t_0 - (t_1 * cos(lambda1))
if (lambda1 <= (-12500.0d0)) then
tmp = atan2(t_2, t_3)
else if (lambda1 <= 2.6d-7) then
tmp = atan2(t_2, (t_0 - (t_1 * cos(lambda2))))
else
tmp = atan2((cos(phi2) * sin(lambda1)), t_3)
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)) * Math.cos(phi2);
double t_3 = t_0 - (t_1 * Math.cos(lambda1));
double tmp;
if (lambda1 <= -12500.0) {
tmp = Math.atan2(t_2, t_3);
} else if (lambda1 <= 2.6e-7) {
tmp = Math.atan2(t_2, (t_0 - (t_1 * Math.cos(lambda2))));
} else {
tmp = Math.atan2((Math.cos(phi2) * Math.sin(lambda1)), t_3);
}
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)) * math.cos(phi2) t_3 = t_0 - (t_1 * math.cos(lambda1)) tmp = 0 if lambda1 <= -12500.0: tmp = math.atan2(t_2, t_3) elif lambda1 <= 2.6e-7: tmp = math.atan2(t_2, (t_0 - (t_1 * math.cos(lambda2)))) else: tmp = math.atan2((math.cos(phi2) * math.sin(lambda1)), t_3) 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(sin(Float64(lambda1 - lambda2)) * cos(phi2)) t_3 = Float64(t_0 - Float64(t_1 * cos(lambda1))) tmp = 0.0 if (lambda1 <= -12500.0) tmp = atan(t_2, t_3); elseif (lambda1 <= 2.6e-7) tmp = atan(t_2, Float64(t_0 - Float64(t_1 * cos(lambda2)))); else tmp = atan(Float64(cos(phi2) * sin(lambda1)), t_3); 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)) * cos(phi2); t_3 = t_0 - (t_1 * cos(lambda1)); tmp = 0.0; if (lambda1 <= -12500.0) tmp = atan2(t_2, t_3); elseif (lambda1 <= 2.6e-7) tmp = atan2(t_2, (t_0 - (t_1 * cos(lambda2)))); else tmp = atan2((cos(phi2) * sin(lambda1)), t_3); 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[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = N[(t$95$0 - N[(t$95$1 * N[Cos[lambda1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[lambda1, -12500.0], N[ArcTan[t$95$2 / t$95$3], $MachinePrecision], If[LessEqual[lambda1, 2.6e-7], N[ArcTan[t$95$2 / N[(t$95$0 - N[(t$95$1 * N[Cos[lambda2], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], N[ArcTan[N[(N[Cos[phi2], $MachinePrecision] * N[Sin[lambda1], $MachinePrecision]), $MachinePrecision] / t$95$3], $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) \cdot \cos \phi_2\\
t_3 := t\_0 - t\_1 \cdot \cos \lambda_1\\
\mathbf{if}\;\lambda_1 \leq -12500:\\
\;\;\;\;\tan^{-1}_* \frac{t\_2}{t\_3}\\
\mathbf{elif}\;\lambda_1 \leq 2.6 \cdot 10^{-7}:\\
\;\;\;\;\tan^{-1}_* \frac{t\_2}{t\_0 - t\_1 \cdot \cos \lambda_2}\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \lambda_1}{t\_3}\\
\end{array}
\end{array}
if lambda1 < -12500Initial program 66.8%
Taylor expanded in lambda2 around 0
lower-cos.f6466.9
Simplified66.9%
if -12500 < lambda1 < 2.59999999999999999e-7Initial program 98.0%
Taylor expanded in lambda1 around 0
cos-negN/A
lower-cos.f6498.0
Simplified98.0%
if 2.59999999999999999e-7 < lambda1 Initial program 62.7%
Taylor expanded in lambda2 around 0
lower-cos.f6462.7
Simplified62.7%
Taylor expanded in lambda2 around 0
lower-sin.f6462.9
Simplified62.9%
Final simplification82.7%
(FPCore (lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (* (cos phi2) (sin phi1)))
(t_1 (* (cos phi2) (sin lambda1)))
(t_2 (* (cos phi1) (sin phi2))))
(if (<= lambda1 -12500.0)
(atan2 t_1 (- t_2 (* t_0 (cos (- lambda1 lambda2)))))
(if (<= lambda1 8.9e-125)
(atan2 (* (cos phi2) (sin (- lambda2))) (- t_2 (* t_0 (cos lambda2))))
(if (<= lambda1 2.6e-7)
(atan2
(* (sin (- lambda1 lambda2)) (cos phi2))
(- t_2 (* (sin phi1) (cos lambda2))))
(atan2 t_1 (- t_2 (* t_0 (cos lambda1)))))))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = cos(phi2) * sin(phi1);
double t_1 = cos(phi2) * sin(lambda1);
double t_2 = cos(phi1) * sin(phi2);
double tmp;
if (lambda1 <= -12500.0) {
tmp = atan2(t_1, (t_2 - (t_0 * cos((lambda1 - lambda2)))));
} else if (lambda1 <= 8.9e-125) {
tmp = atan2((cos(phi2) * sin(-lambda2)), (t_2 - (t_0 * cos(lambda2))));
} else if (lambda1 <= 2.6e-7) {
tmp = atan2((sin((lambda1 - lambda2)) * cos(phi2)), (t_2 - (sin(phi1) * cos(lambda2))));
} else {
tmp = atan2(t_1, (t_2 - (t_0 * cos(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) :: t_0
real(8) :: t_1
real(8) :: t_2
real(8) :: tmp
t_0 = cos(phi2) * sin(phi1)
t_1 = cos(phi2) * sin(lambda1)
t_2 = cos(phi1) * sin(phi2)
if (lambda1 <= (-12500.0d0)) then
tmp = atan2(t_1, (t_2 - (t_0 * cos((lambda1 - lambda2)))))
else if (lambda1 <= 8.9d-125) then
tmp = atan2((cos(phi2) * sin(-lambda2)), (t_2 - (t_0 * cos(lambda2))))
else if (lambda1 <= 2.6d-7) then
tmp = atan2((sin((lambda1 - lambda2)) * cos(phi2)), (t_2 - (sin(phi1) * cos(lambda2))))
else
tmp = atan2(t_1, (t_2 - (t_0 * cos(lambda1))))
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(phi1);
double t_1 = Math.cos(phi2) * Math.sin(lambda1);
double t_2 = Math.cos(phi1) * Math.sin(phi2);
double tmp;
if (lambda1 <= -12500.0) {
tmp = Math.atan2(t_1, (t_2 - (t_0 * Math.cos((lambda1 - lambda2)))));
} else if (lambda1 <= 8.9e-125) {
tmp = Math.atan2((Math.cos(phi2) * Math.sin(-lambda2)), (t_2 - (t_0 * Math.cos(lambda2))));
} else if (lambda1 <= 2.6e-7) {
tmp = Math.atan2((Math.sin((lambda1 - lambda2)) * Math.cos(phi2)), (t_2 - (Math.sin(phi1) * Math.cos(lambda2))));
} else {
tmp = Math.atan2(t_1, (t_2 - (t_0 * Math.cos(lambda1))));
}
return tmp;
}
def code(lambda1, lambda2, phi1, phi2): t_0 = math.cos(phi2) * math.sin(phi1) t_1 = math.cos(phi2) * math.sin(lambda1) t_2 = math.cos(phi1) * math.sin(phi2) tmp = 0 if lambda1 <= -12500.0: tmp = math.atan2(t_1, (t_2 - (t_0 * math.cos((lambda1 - lambda2))))) elif lambda1 <= 8.9e-125: tmp = math.atan2((math.cos(phi2) * math.sin(-lambda2)), (t_2 - (t_0 * math.cos(lambda2)))) elif lambda1 <= 2.6e-7: tmp = math.atan2((math.sin((lambda1 - lambda2)) * math.cos(phi2)), (t_2 - (math.sin(phi1) * math.cos(lambda2)))) else: tmp = math.atan2(t_1, (t_2 - (t_0 * math.cos(lambda1)))) return tmp
function code(lambda1, lambda2, phi1, phi2) t_0 = Float64(cos(phi2) * sin(phi1)) t_1 = Float64(cos(phi2) * sin(lambda1)) t_2 = Float64(cos(phi1) * sin(phi2)) tmp = 0.0 if (lambda1 <= -12500.0) tmp = atan(t_1, Float64(t_2 - Float64(t_0 * cos(Float64(lambda1 - lambda2))))); elseif (lambda1 <= 8.9e-125) tmp = atan(Float64(cos(phi2) * sin(Float64(-lambda2))), Float64(t_2 - Float64(t_0 * cos(lambda2)))); elseif (lambda1 <= 2.6e-7) tmp = atan(Float64(sin(Float64(lambda1 - lambda2)) * cos(phi2)), Float64(t_2 - Float64(sin(phi1) * cos(lambda2)))); else tmp = atan(t_1, Float64(t_2 - Float64(t_0 * cos(lambda1)))); end return tmp end
function tmp_2 = code(lambda1, lambda2, phi1, phi2) t_0 = cos(phi2) * sin(phi1); t_1 = cos(phi2) * sin(lambda1); t_2 = cos(phi1) * sin(phi2); tmp = 0.0; if (lambda1 <= -12500.0) tmp = atan2(t_1, (t_2 - (t_0 * cos((lambda1 - lambda2))))); elseif (lambda1 <= 8.9e-125) tmp = atan2((cos(phi2) * sin(-lambda2)), (t_2 - (t_0 * cos(lambda2)))); elseif (lambda1 <= 2.6e-7) tmp = atan2((sin((lambda1 - lambda2)) * cos(phi2)), (t_2 - (sin(phi1) * cos(lambda2)))); else tmp = atan2(t_1, (t_2 - (t_0 * cos(lambda1)))); end tmp_2 = tmp; end
code[lambda1_, lambda2_, phi1_, phi2_] := Block[{t$95$0 = N[(N[Cos[phi2], $MachinePrecision] * N[Sin[phi1], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[Cos[phi2], $MachinePrecision] * N[Sin[lambda1], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(N[Cos[phi1], $MachinePrecision] * N[Sin[phi2], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[lambda1, -12500.0], N[ArcTan[t$95$1 / N[(t$95$2 - N[(t$95$0 * N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], If[LessEqual[lambda1, 8.9e-125], N[ArcTan[N[(N[Cos[phi2], $MachinePrecision] * N[Sin[(-lambda2)], $MachinePrecision]), $MachinePrecision] / N[(t$95$2 - N[(t$95$0 * N[Cos[lambda2], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], If[LessEqual[lambda1, 2.6e-7], N[ArcTan[N[(N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision] / N[(t$95$2 - N[(N[Sin[phi1], $MachinePrecision] * N[Cos[lambda2], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], N[ArcTan[t$95$1 / N[(t$95$2 - N[(t$95$0 * N[Cos[lambda1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \cos \phi_2 \cdot \sin \phi_1\\
t_1 := \cos \phi_2 \cdot \sin \lambda_1\\
t_2 := \cos \phi_1 \cdot \sin \phi_2\\
\mathbf{if}\;\lambda_1 \leq -12500:\\
\;\;\;\;\tan^{-1}_* \frac{t\_1}{t\_2 - t\_0 \cdot \cos \left(\lambda_1 - \lambda_2\right)}\\
\mathbf{elif}\;\lambda_1 \leq 8.9 \cdot 10^{-125}:\\
\;\;\;\;\tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(-\lambda_2\right)}{t\_2 - t\_0 \cdot \cos \lambda_2}\\
\mathbf{elif}\;\lambda_1 \leq 2.6 \cdot 10^{-7}:\\
\;\;\;\;\tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2}{t\_2 - \sin \phi_1 \cdot \cos \lambda_2}\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1}_* \frac{t\_1}{t\_2 - t\_0 \cdot \cos \lambda_1}\\
\end{array}
\end{array}
if lambda1 < -12500Initial program 66.8%
Taylor expanded in lambda2 around 0
lower-sin.f6464.9
Simplified64.9%
if -12500 < lambda1 < 8.90000000000000001e-125Initial program 97.9%
Taylor expanded in lambda1 around 0
cos-negN/A
lower-cos.f6497.9
Simplified97.9%
Taylor expanded in lambda1 around 0
lower-sin.f64N/A
lower-neg.f6488.1
Simplified88.1%
if 8.90000000000000001e-125 < lambda1 < 2.59999999999999999e-7Initial program 98.7%
Taylor expanded in phi2 around 0
lower-sin.f6495.0
Simplified95.0%
Taylor expanded in lambda1 around 0
cos-negN/A
lower-cos.f6495.0
Simplified95.0%
if 2.59999999999999999e-7 < lambda1 Initial program 62.7%
Taylor expanded in lambda2 around 0
lower-cos.f6462.7
Simplified62.7%
Taylor expanded in lambda2 around 0
lower-sin.f6462.9
Simplified62.9%
Final simplification77.5%
(FPCore (lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (* (cos phi2) (sin phi1)))
(t_1 (* (cos phi1) (sin phi2)))
(t_2
(atan2 (* (cos phi2) (sin lambda1)) (- t_1 (* t_0 (cos lambda1))))))
(if (<= lambda1 -12500.0)
t_2
(if (<= lambda1 8.9e-125)
(atan2 (* (cos phi2) (sin (- lambda2))) (- t_1 (* t_0 (cos lambda2))))
(if (<= lambda1 2.6e-7)
(atan2
(* (sin (- lambda1 lambda2)) (cos phi2))
(- t_1 (* (sin phi1) (cos lambda2))))
t_2)))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = cos(phi2) * sin(phi1);
double t_1 = cos(phi1) * sin(phi2);
double t_2 = atan2((cos(phi2) * sin(lambda1)), (t_1 - (t_0 * cos(lambda1))));
double tmp;
if (lambda1 <= -12500.0) {
tmp = t_2;
} else if (lambda1 <= 8.9e-125) {
tmp = atan2((cos(phi2) * sin(-lambda2)), (t_1 - (t_0 * cos(lambda2))));
} else if (lambda1 <= 2.6e-7) {
tmp = atan2((sin((lambda1 - lambda2)) * cos(phi2)), (t_1 - (sin(phi1) * cos(lambda2))));
} else {
tmp = t_2;
}
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(phi2) * sin(phi1)
t_1 = cos(phi1) * sin(phi2)
t_2 = atan2((cos(phi2) * sin(lambda1)), (t_1 - (t_0 * cos(lambda1))))
if (lambda1 <= (-12500.0d0)) then
tmp = t_2
else if (lambda1 <= 8.9d-125) then
tmp = atan2((cos(phi2) * sin(-lambda2)), (t_1 - (t_0 * cos(lambda2))))
else if (lambda1 <= 2.6d-7) then
tmp = atan2((sin((lambda1 - lambda2)) * cos(phi2)), (t_1 - (sin(phi1) * cos(lambda2))))
else
tmp = t_2
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(phi1);
double t_1 = Math.cos(phi1) * Math.sin(phi2);
double t_2 = Math.atan2((Math.cos(phi2) * Math.sin(lambda1)), (t_1 - (t_0 * Math.cos(lambda1))));
double tmp;
if (lambda1 <= -12500.0) {
tmp = t_2;
} else if (lambda1 <= 8.9e-125) {
tmp = Math.atan2((Math.cos(phi2) * Math.sin(-lambda2)), (t_1 - (t_0 * Math.cos(lambda2))));
} else if (lambda1 <= 2.6e-7) {
tmp = Math.atan2((Math.sin((lambda1 - lambda2)) * Math.cos(phi2)), (t_1 - (Math.sin(phi1) * Math.cos(lambda2))));
} else {
tmp = t_2;
}
return tmp;
}
def code(lambda1, lambda2, phi1, phi2): t_0 = math.cos(phi2) * math.sin(phi1) t_1 = math.cos(phi1) * math.sin(phi2) t_2 = math.atan2((math.cos(phi2) * math.sin(lambda1)), (t_1 - (t_0 * math.cos(lambda1)))) tmp = 0 if lambda1 <= -12500.0: tmp = t_2 elif lambda1 <= 8.9e-125: tmp = math.atan2((math.cos(phi2) * math.sin(-lambda2)), (t_1 - (t_0 * math.cos(lambda2)))) elif lambda1 <= 2.6e-7: tmp = math.atan2((math.sin((lambda1 - lambda2)) * math.cos(phi2)), (t_1 - (math.sin(phi1) * math.cos(lambda2)))) else: tmp = t_2 return tmp
function code(lambda1, lambda2, phi1, phi2) t_0 = Float64(cos(phi2) * sin(phi1)) t_1 = Float64(cos(phi1) * sin(phi2)) t_2 = atan(Float64(cos(phi2) * sin(lambda1)), Float64(t_1 - Float64(t_0 * cos(lambda1)))) tmp = 0.0 if (lambda1 <= -12500.0) tmp = t_2; elseif (lambda1 <= 8.9e-125) tmp = atan(Float64(cos(phi2) * sin(Float64(-lambda2))), Float64(t_1 - Float64(t_0 * cos(lambda2)))); elseif (lambda1 <= 2.6e-7) tmp = atan(Float64(sin(Float64(lambda1 - lambda2)) * cos(phi2)), Float64(t_1 - Float64(sin(phi1) * cos(lambda2)))); else tmp = t_2; end return tmp end
function tmp_2 = code(lambda1, lambda2, phi1, phi2) t_0 = cos(phi2) * sin(phi1); t_1 = cos(phi1) * sin(phi2); t_2 = atan2((cos(phi2) * sin(lambda1)), (t_1 - (t_0 * cos(lambda1)))); tmp = 0.0; if (lambda1 <= -12500.0) tmp = t_2; elseif (lambda1 <= 8.9e-125) tmp = atan2((cos(phi2) * sin(-lambda2)), (t_1 - (t_0 * cos(lambda2)))); elseif (lambda1 <= 2.6e-7) tmp = atan2((sin((lambda1 - lambda2)) * cos(phi2)), (t_1 - (sin(phi1) * cos(lambda2)))); else tmp = t_2; end tmp_2 = tmp; end
code[lambda1_, lambda2_, phi1_, phi2_] := Block[{t$95$0 = N[(N[Cos[phi2], $MachinePrecision] * N[Sin[phi1], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[Cos[phi1], $MachinePrecision] * N[Sin[phi2], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[ArcTan[N[(N[Cos[phi2], $MachinePrecision] * N[Sin[lambda1], $MachinePrecision]), $MachinePrecision] / N[(t$95$1 - N[(t$95$0 * N[Cos[lambda1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]}, If[LessEqual[lambda1, -12500.0], t$95$2, If[LessEqual[lambda1, 8.9e-125], N[ArcTan[N[(N[Cos[phi2], $MachinePrecision] * N[Sin[(-lambda2)], $MachinePrecision]), $MachinePrecision] / N[(t$95$1 - N[(t$95$0 * N[Cos[lambda2], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], If[LessEqual[lambda1, 2.6e-7], N[ArcTan[N[(N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision] / N[(t$95$1 - N[(N[Sin[phi1], $MachinePrecision] * N[Cos[lambda2], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], t$95$2]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \cos \phi_2 \cdot \sin \phi_1\\
t_1 := \cos \phi_1 \cdot \sin \phi_2\\
t_2 := \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \lambda_1}{t\_1 - t\_0 \cdot \cos \lambda_1}\\
\mathbf{if}\;\lambda_1 \leq -12500:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;\lambda_1 \leq 8.9 \cdot 10^{-125}:\\
\;\;\;\;\tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(-\lambda_2\right)}{t\_1 - t\_0 \cdot \cos \lambda_2}\\
\mathbf{elif}\;\lambda_1 \leq 2.6 \cdot 10^{-7}:\\
\;\;\;\;\tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2}{t\_1 - \sin \phi_1 \cdot \cos \lambda_2}\\
\mathbf{else}:\\
\;\;\;\;t\_2\\
\end{array}
\end{array}
if lambda1 < -12500 or 2.59999999999999999e-7 < lambda1 Initial program 64.6%
Taylor expanded in lambda2 around 0
lower-cos.f6464.7
Simplified64.7%
Taylor expanded in lambda2 around 0
lower-sin.f6463.9
Simplified63.9%
if -12500 < lambda1 < 8.90000000000000001e-125Initial program 97.9%
Taylor expanded in lambda1 around 0
cos-negN/A
lower-cos.f6497.9
Simplified97.9%
Taylor expanded in lambda1 around 0
lower-sin.f64N/A
lower-neg.f6488.1
Simplified88.1%
if 8.90000000000000001e-125 < lambda1 < 2.59999999999999999e-7Initial program 98.7%
Taylor expanded in phi2 around 0
lower-sin.f6495.0
Simplified95.0%
Taylor expanded in lambda1 around 0
cos-negN/A
lower-cos.f6495.0
Simplified95.0%
Final simplification77.5%
(FPCore (lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (* (sin phi1) (cos lambda2)))
(t_1 (* (cos phi1) (sin phi2)))
(t_2
(atan2
(* (cos phi2) (sin lambda1))
(- t_1 (* (* (cos phi2) (sin phi1)) (cos lambda1))))))
(if (<= lambda1 -12500.0)
t_2
(if (<= lambda1 8.9e-125)
(atan2 (* (cos phi2) (sin (- lambda2))) (- t_1 (* (cos phi2) t_0)))
(if (<= lambda1 2.6e-7)
(atan2 (* (sin (- lambda1 lambda2)) (cos phi2)) (- t_1 t_0))
t_2)))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = sin(phi1) * cos(lambda2);
double t_1 = cos(phi1) * sin(phi2);
double t_2 = atan2((cos(phi2) * sin(lambda1)), (t_1 - ((cos(phi2) * sin(phi1)) * cos(lambda1))));
double tmp;
if (lambda1 <= -12500.0) {
tmp = t_2;
} else if (lambda1 <= 8.9e-125) {
tmp = atan2((cos(phi2) * sin(-lambda2)), (t_1 - (cos(phi2) * t_0)));
} else if (lambda1 <= 2.6e-7) {
tmp = atan2((sin((lambda1 - lambda2)) * cos(phi2)), (t_1 - t_0));
} else {
tmp = t_2;
}
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 = sin(phi1) * cos(lambda2)
t_1 = cos(phi1) * sin(phi2)
t_2 = atan2((cos(phi2) * sin(lambda1)), (t_1 - ((cos(phi2) * sin(phi1)) * cos(lambda1))))
if (lambda1 <= (-12500.0d0)) then
tmp = t_2
else if (lambda1 <= 8.9d-125) then
tmp = atan2((cos(phi2) * sin(-lambda2)), (t_1 - (cos(phi2) * t_0)))
else if (lambda1 <= 2.6d-7) then
tmp = atan2((sin((lambda1 - lambda2)) * cos(phi2)), (t_1 - t_0))
else
tmp = t_2
end if
code = tmp
end function
public static double code(double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = Math.sin(phi1) * Math.cos(lambda2);
double t_1 = Math.cos(phi1) * Math.sin(phi2);
double t_2 = Math.atan2((Math.cos(phi2) * Math.sin(lambda1)), (t_1 - ((Math.cos(phi2) * Math.sin(phi1)) * Math.cos(lambda1))));
double tmp;
if (lambda1 <= -12500.0) {
tmp = t_2;
} else if (lambda1 <= 8.9e-125) {
tmp = Math.atan2((Math.cos(phi2) * Math.sin(-lambda2)), (t_1 - (Math.cos(phi2) * t_0)));
} else if (lambda1 <= 2.6e-7) {
tmp = Math.atan2((Math.sin((lambda1 - lambda2)) * Math.cos(phi2)), (t_1 - t_0));
} else {
tmp = t_2;
}
return tmp;
}
def code(lambda1, lambda2, phi1, phi2): t_0 = math.sin(phi1) * math.cos(lambda2) t_1 = math.cos(phi1) * math.sin(phi2) t_2 = math.atan2((math.cos(phi2) * math.sin(lambda1)), (t_1 - ((math.cos(phi2) * math.sin(phi1)) * math.cos(lambda1)))) tmp = 0 if lambda1 <= -12500.0: tmp = t_2 elif lambda1 <= 8.9e-125: tmp = math.atan2((math.cos(phi2) * math.sin(-lambda2)), (t_1 - (math.cos(phi2) * t_0))) elif lambda1 <= 2.6e-7: tmp = math.atan2((math.sin((lambda1 - lambda2)) * math.cos(phi2)), (t_1 - t_0)) else: tmp = t_2 return tmp
function code(lambda1, lambda2, phi1, phi2) t_0 = Float64(sin(phi1) * cos(lambda2)) t_1 = Float64(cos(phi1) * sin(phi2)) t_2 = atan(Float64(cos(phi2) * sin(lambda1)), Float64(t_1 - Float64(Float64(cos(phi2) * sin(phi1)) * cos(lambda1)))) tmp = 0.0 if (lambda1 <= -12500.0) tmp = t_2; elseif (lambda1 <= 8.9e-125) tmp = atan(Float64(cos(phi2) * sin(Float64(-lambda2))), Float64(t_1 - Float64(cos(phi2) * t_0))); elseif (lambda1 <= 2.6e-7) tmp = atan(Float64(sin(Float64(lambda1 - lambda2)) * cos(phi2)), Float64(t_1 - t_0)); else tmp = t_2; end return tmp end
function tmp_2 = code(lambda1, lambda2, phi1, phi2) t_0 = sin(phi1) * cos(lambda2); t_1 = cos(phi1) * sin(phi2); t_2 = atan2((cos(phi2) * sin(lambda1)), (t_1 - ((cos(phi2) * sin(phi1)) * cos(lambda1)))); tmp = 0.0; if (lambda1 <= -12500.0) tmp = t_2; elseif (lambda1 <= 8.9e-125) tmp = atan2((cos(phi2) * sin(-lambda2)), (t_1 - (cos(phi2) * t_0))); elseif (lambda1 <= 2.6e-7) tmp = atan2((sin((lambda1 - lambda2)) * cos(phi2)), (t_1 - t_0)); else tmp = t_2; end tmp_2 = tmp; end
code[lambda1_, lambda2_, phi1_, phi2_] := Block[{t$95$0 = N[(N[Sin[phi1], $MachinePrecision] * N[Cos[lambda2], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[Cos[phi1], $MachinePrecision] * N[Sin[phi2], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[ArcTan[N[(N[Cos[phi2], $MachinePrecision] * N[Sin[lambda1], $MachinePrecision]), $MachinePrecision] / N[(t$95$1 - N[(N[(N[Cos[phi2], $MachinePrecision] * N[Sin[phi1], $MachinePrecision]), $MachinePrecision] * N[Cos[lambda1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]}, If[LessEqual[lambda1, -12500.0], t$95$2, If[LessEqual[lambda1, 8.9e-125], N[ArcTan[N[(N[Cos[phi2], $MachinePrecision] * N[Sin[(-lambda2)], $MachinePrecision]), $MachinePrecision] / N[(t$95$1 - N[(N[Cos[phi2], $MachinePrecision] * t$95$0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], If[LessEqual[lambda1, 2.6e-7], N[ArcTan[N[(N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision] / N[(t$95$1 - t$95$0), $MachinePrecision]], $MachinePrecision], t$95$2]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sin \phi_1 \cdot \cos \lambda_2\\
t_1 := \cos \phi_1 \cdot \sin \phi_2\\
t_2 := \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \lambda_1}{t\_1 - \left(\cos \phi_2 \cdot \sin \phi_1\right) \cdot \cos \lambda_1}\\
\mathbf{if}\;\lambda_1 \leq -12500:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;\lambda_1 \leq 8.9 \cdot 10^{-125}:\\
\;\;\;\;\tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(-\lambda_2\right)}{t\_1 - \cos \phi_2 \cdot t\_0}\\
\mathbf{elif}\;\lambda_1 \leq 2.6 \cdot 10^{-7}:\\
\;\;\;\;\tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2}{t\_1 - t\_0}\\
\mathbf{else}:\\
\;\;\;\;t\_2\\
\end{array}
\end{array}
if lambda1 < -12500 or 2.59999999999999999e-7 < lambda1 Initial program 64.6%
Taylor expanded in lambda2 around 0
lower-cos.f6464.7
Simplified64.7%
Taylor expanded in lambda2 around 0
lower-sin.f6463.9
Simplified63.9%
if -12500 < lambda1 < 8.90000000000000001e-125Initial program 97.9%
Taylor expanded in phi1 around inf
lower-*.f64N/A
lower-cos.f64N/A
lower-*.f64N/A
lower-cos.f64N/A
lower--.f64N/A
lower-sin.f6497.9
Simplified97.9%
Taylor expanded in lambda1 around 0
cos-negN/A
lower-cos.f6497.9
Simplified97.9%
Taylor expanded in lambda1 around 0
lower-sin.f64N/A
lower-neg.f6488.1
Simplified88.1%
if 8.90000000000000001e-125 < lambda1 < 2.59999999999999999e-7Initial program 98.7%
Taylor expanded in phi2 around 0
lower-sin.f6495.0
Simplified95.0%
Taylor expanded in lambda1 around 0
cos-negN/A
lower-cos.f6495.0
Simplified95.0%
Final simplification77.5%
(FPCore (lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (* (cos phi1) (sin phi2)))
(t_1 (* (sin (- lambda1 lambda2)) (cos phi2)))
(t_2 (- t_0 (* (* (cos phi2) (sin phi1)) (cos lambda1)))))
(if (<= lambda1 -12500.0)
(atan2 t_1 t_2)
(if (<= lambda1 2.6e-7)
(atan2 t_1 (- t_0 (* (cos phi2) (* (sin phi1) (cos lambda2)))))
(atan2 (* (cos phi2) (sin lambda1)) t_2)))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = cos(phi1) * sin(phi2);
double t_1 = sin((lambda1 - lambda2)) * cos(phi2);
double t_2 = t_0 - ((cos(phi2) * sin(phi1)) * cos(lambda1));
double tmp;
if (lambda1 <= -12500.0) {
tmp = atan2(t_1, t_2);
} else if (lambda1 <= 2.6e-7) {
tmp = atan2(t_1, (t_0 - (cos(phi2) * (sin(phi1) * cos(lambda2)))));
} else {
tmp = atan2((cos(phi2) * sin(lambda1)), t_2);
}
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 = sin((lambda1 - lambda2)) * cos(phi2)
t_2 = t_0 - ((cos(phi2) * sin(phi1)) * cos(lambda1))
if (lambda1 <= (-12500.0d0)) then
tmp = atan2(t_1, t_2)
else if (lambda1 <= 2.6d-7) then
tmp = atan2(t_1, (t_0 - (cos(phi2) * (sin(phi1) * cos(lambda2)))))
else
tmp = atan2((cos(phi2) * sin(lambda1)), t_2)
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.sin((lambda1 - lambda2)) * Math.cos(phi2);
double t_2 = t_0 - ((Math.cos(phi2) * Math.sin(phi1)) * Math.cos(lambda1));
double tmp;
if (lambda1 <= -12500.0) {
tmp = Math.atan2(t_1, t_2);
} else if (lambda1 <= 2.6e-7) {
tmp = Math.atan2(t_1, (t_0 - (Math.cos(phi2) * (Math.sin(phi1) * Math.cos(lambda2)))));
} else {
tmp = Math.atan2((Math.cos(phi2) * Math.sin(lambda1)), t_2);
}
return tmp;
}
def code(lambda1, lambda2, phi1, phi2): t_0 = math.cos(phi1) * math.sin(phi2) t_1 = math.sin((lambda1 - lambda2)) * math.cos(phi2) t_2 = t_0 - ((math.cos(phi2) * math.sin(phi1)) * math.cos(lambda1)) tmp = 0 if lambda1 <= -12500.0: tmp = math.atan2(t_1, t_2) elif lambda1 <= 2.6e-7: tmp = math.atan2(t_1, (t_0 - (math.cos(phi2) * (math.sin(phi1) * math.cos(lambda2))))) else: tmp = math.atan2((math.cos(phi2) * math.sin(lambda1)), t_2) return tmp
function code(lambda1, lambda2, phi1, phi2) t_0 = Float64(cos(phi1) * sin(phi2)) t_1 = Float64(sin(Float64(lambda1 - lambda2)) * cos(phi2)) t_2 = Float64(t_0 - Float64(Float64(cos(phi2) * sin(phi1)) * cos(lambda1))) tmp = 0.0 if (lambda1 <= -12500.0) tmp = atan(t_1, t_2); elseif (lambda1 <= 2.6e-7) tmp = atan(t_1, Float64(t_0 - Float64(cos(phi2) * Float64(sin(phi1) * cos(lambda2))))); else tmp = atan(Float64(cos(phi2) * sin(lambda1)), t_2); end return tmp end
function tmp_2 = code(lambda1, lambda2, phi1, phi2) t_0 = cos(phi1) * sin(phi2); t_1 = sin((lambda1 - lambda2)) * cos(phi2); t_2 = t_0 - ((cos(phi2) * sin(phi1)) * cos(lambda1)); tmp = 0.0; if (lambda1 <= -12500.0) tmp = atan2(t_1, t_2); elseif (lambda1 <= 2.6e-7) tmp = atan2(t_1, (t_0 - (cos(phi2) * (sin(phi1) * cos(lambda2))))); else tmp = atan2((cos(phi2) * sin(lambda1)), t_2); 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[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(t$95$0 - N[(N[(N[Cos[phi2], $MachinePrecision] * N[Sin[phi1], $MachinePrecision]), $MachinePrecision] * N[Cos[lambda1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[lambda1, -12500.0], N[ArcTan[t$95$1 / t$95$2], $MachinePrecision], If[LessEqual[lambda1, 2.6e-7], N[ArcTan[t$95$1 / N[(t$95$0 - N[(N[Cos[phi2], $MachinePrecision] * N[(N[Sin[phi1], $MachinePrecision] * N[Cos[lambda2], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], N[ArcTan[N[(N[Cos[phi2], $MachinePrecision] * N[Sin[lambda1], $MachinePrecision]), $MachinePrecision] / t$95$2], $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \cos \phi_1 \cdot \sin \phi_2\\
t_1 := \sin \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2\\
t_2 := t\_0 - \left(\cos \phi_2 \cdot \sin \phi_1\right) \cdot \cos \lambda_1\\
\mathbf{if}\;\lambda_1 \leq -12500:\\
\;\;\;\;\tan^{-1}_* \frac{t\_1}{t\_2}\\
\mathbf{elif}\;\lambda_1 \leq 2.6 \cdot 10^{-7}:\\
\;\;\;\;\tan^{-1}_* \frac{t\_1}{t\_0 - \cos \phi_2 \cdot \left(\sin \phi_1 \cdot \cos \lambda_2\right)}\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \lambda_1}{t\_2}\\
\end{array}
\end{array}
if lambda1 < -12500Initial program 66.8%
Taylor expanded in lambda2 around 0
lower-cos.f6466.9
Simplified66.9%
if -12500 < lambda1 < 2.59999999999999999e-7Initial program 98.0%
Taylor expanded in phi1 around inf
lower-*.f64N/A
lower-cos.f64N/A
lower-*.f64N/A
lower-cos.f64N/A
lower--.f64N/A
lower-sin.f6498.0
Simplified98.0%
Taylor expanded in lambda1 around 0
cos-negN/A
lower-cos.f6498.0
Simplified98.0%
if 2.59999999999999999e-7 < lambda1 Initial program 62.7%
Taylor expanded in lambda2 around 0
lower-cos.f6462.7
Simplified62.7%
Taylor expanded in lambda2 around 0
lower-sin.f6462.9
Simplified62.9%
Final simplification82.7%
(FPCore (lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (* (cos phi2) (sin phi1)))
(t_1 (* (cos phi2) (sin lambda1)))
(t_2 (* (cos phi1) (sin phi2))))
(if (<= lambda1 -3.2e+93)
(atan2 t_1 (- t_2 (* t_0 (cos (- lambda1 lambda2)))))
(if (<= lambda1 2.6e-7)
(atan2
(* (sin (- lambda1 lambda2)) (cos phi2))
(- t_2 (* (cos phi2) (* (sin phi1) (cos lambda2)))))
(atan2 t_1 (- t_2 (* t_0 (cos lambda1))))))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = cos(phi2) * sin(phi1);
double t_1 = cos(phi2) * sin(lambda1);
double t_2 = cos(phi1) * sin(phi2);
double tmp;
if (lambda1 <= -3.2e+93) {
tmp = atan2(t_1, (t_2 - (t_0 * cos((lambda1 - lambda2)))));
} else if (lambda1 <= 2.6e-7) {
tmp = atan2((sin((lambda1 - lambda2)) * cos(phi2)), (t_2 - (cos(phi2) * (sin(phi1) * cos(lambda2)))));
} else {
tmp = atan2(t_1, (t_2 - (t_0 * cos(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) :: t_0
real(8) :: t_1
real(8) :: t_2
real(8) :: tmp
t_0 = cos(phi2) * sin(phi1)
t_1 = cos(phi2) * sin(lambda1)
t_2 = cos(phi1) * sin(phi2)
if (lambda1 <= (-3.2d+93)) then
tmp = atan2(t_1, (t_2 - (t_0 * cos((lambda1 - lambda2)))))
else if (lambda1 <= 2.6d-7) then
tmp = atan2((sin((lambda1 - lambda2)) * cos(phi2)), (t_2 - (cos(phi2) * (sin(phi1) * cos(lambda2)))))
else
tmp = atan2(t_1, (t_2 - (t_0 * cos(lambda1))))
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(phi1);
double t_1 = Math.cos(phi2) * Math.sin(lambda1);
double t_2 = Math.cos(phi1) * Math.sin(phi2);
double tmp;
if (lambda1 <= -3.2e+93) {
tmp = Math.atan2(t_1, (t_2 - (t_0 * Math.cos((lambda1 - lambda2)))));
} else if (lambda1 <= 2.6e-7) {
tmp = Math.atan2((Math.sin((lambda1 - lambda2)) * Math.cos(phi2)), (t_2 - (Math.cos(phi2) * (Math.sin(phi1) * Math.cos(lambda2)))));
} else {
tmp = Math.atan2(t_1, (t_2 - (t_0 * Math.cos(lambda1))));
}
return tmp;
}
def code(lambda1, lambda2, phi1, phi2): t_0 = math.cos(phi2) * math.sin(phi1) t_1 = math.cos(phi2) * math.sin(lambda1) t_2 = math.cos(phi1) * math.sin(phi2) tmp = 0 if lambda1 <= -3.2e+93: tmp = math.atan2(t_1, (t_2 - (t_0 * math.cos((lambda1 - lambda2))))) elif lambda1 <= 2.6e-7: tmp = math.atan2((math.sin((lambda1 - lambda2)) * math.cos(phi2)), (t_2 - (math.cos(phi2) * (math.sin(phi1) * math.cos(lambda2))))) else: tmp = math.atan2(t_1, (t_2 - (t_0 * math.cos(lambda1)))) return tmp
function code(lambda1, lambda2, phi1, phi2) t_0 = Float64(cos(phi2) * sin(phi1)) t_1 = Float64(cos(phi2) * sin(lambda1)) t_2 = Float64(cos(phi1) * sin(phi2)) tmp = 0.0 if (lambda1 <= -3.2e+93) tmp = atan(t_1, Float64(t_2 - Float64(t_0 * cos(Float64(lambda1 - lambda2))))); elseif (lambda1 <= 2.6e-7) tmp = atan(Float64(sin(Float64(lambda1 - lambda2)) * cos(phi2)), Float64(t_2 - Float64(cos(phi2) * Float64(sin(phi1) * cos(lambda2))))); else tmp = atan(t_1, Float64(t_2 - Float64(t_0 * cos(lambda1)))); end return tmp end
function tmp_2 = code(lambda1, lambda2, phi1, phi2) t_0 = cos(phi2) * sin(phi1); t_1 = cos(phi2) * sin(lambda1); t_2 = cos(phi1) * sin(phi2); tmp = 0.0; if (lambda1 <= -3.2e+93) tmp = atan2(t_1, (t_2 - (t_0 * cos((lambda1 - lambda2))))); elseif (lambda1 <= 2.6e-7) tmp = atan2((sin((lambda1 - lambda2)) * cos(phi2)), (t_2 - (cos(phi2) * (sin(phi1) * cos(lambda2))))); else tmp = atan2(t_1, (t_2 - (t_0 * cos(lambda1)))); end tmp_2 = tmp; end
code[lambda1_, lambda2_, phi1_, phi2_] := Block[{t$95$0 = N[(N[Cos[phi2], $MachinePrecision] * N[Sin[phi1], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[Cos[phi2], $MachinePrecision] * N[Sin[lambda1], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(N[Cos[phi1], $MachinePrecision] * N[Sin[phi2], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[lambda1, -3.2e+93], N[ArcTan[t$95$1 / N[(t$95$2 - N[(t$95$0 * N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], If[LessEqual[lambda1, 2.6e-7], N[ArcTan[N[(N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision] / N[(t$95$2 - N[(N[Cos[phi2], $MachinePrecision] * N[(N[Sin[phi1], $MachinePrecision] * N[Cos[lambda2], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], N[ArcTan[t$95$1 / N[(t$95$2 - N[(t$95$0 * N[Cos[lambda1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \cos \phi_2 \cdot \sin \phi_1\\
t_1 := \cos \phi_2 \cdot \sin \lambda_1\\
t_2 := \cos \phi_1 \cdot \sin \phi_2\\
\mathbf{if}\;\lambda_1 \leq -3.2 \cdot 10^{+93}:\\
\;\;\;\;\tan^{-1}_* \frac{t\_1}{t\_2 - t\_0 \cdot \cos \left(\lambda_1 - \lambda_2\right)}\\
\mathbf{elif}\;\lambda_1 \leq 2.6 \cdot 10^{-7}:\\
\;\;\;\;\tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2}{t\_2 - \cos \phi_2 \cdot \left(\sin \phi_1 \cdot \cos \lambda_2\right)}\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1}_* \frac{t\_1}{t\_2 - t\_0 \cdot \cos \lambda_1}\\
\end{array}
\end{array}
if lambda1 < -3.2000000000000001e93Initial program 72.4%
Taylor expanded in lambda2 around 0
lower-sin.f6471.8
Simplified71.8%
if -3.2000000000000001e93 < lambda1 < 2.59999999999999999e-7Initial program 93.9%
Taylor expanded in phi1 around inf
lower-*.f64N/A
lower-cos.f64N/A
lower-*.f64N/A
lower-cos.f64N/A
lower--.f64N/A
lower-sin.f6493.9
Simplified93.9%
Taylor expanded in lambda1 around 0
cos-negN/A
lower-cos.f6493.9
Simplified93.9%
if 2.59999999999999999e-7 < lambda1 Initial program 62.7%
Taylor expanded in lambda2 around 0
lower-cos.f6462.7
Simplified62.7%
Taylor expanded in lambda2 around 0
lower-sin.f6462.9
Simplified62.9%
Final simplification82.6%
(FPCore (lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (* (cos phi1) (sin phi2)))
(t_1
(atan2
(* (cos phi2) (sin lambda1))
(- t_0 (* (* (cos phi2) (sin phi1)) (cos lambda1))))))
(if (<= lambda1 -7.2e+93)
t_1
(if (<= lambda1 2.6e-7)
(atan2
(* (sin (- lambda1 lambda2)) (cos phi2))
(- t_0 (* (sin phi1) (cos (- lambda1 lambda2)))))
t_1))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = cos(phi1) * sin(phi2);
double t_1 = atan2((cos(phi2) * sin(lambda1)), (t_0 - ((cos(phi2) * sin(phi1)) * cos(lambda1))));
double tmp;
if (lambda1 <= -7.2e+93) {
tmp = t_1;
} else if (lambda1 <= 2.6e-7) {
tmp = atan2((sin((lambda1 - lambda2)) * cos(phi2)), (t_0 - (sin(phi1) * cos((lambda1 - lambda2)))));
} else {
tmp = 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) :: tmp
t_0 = cos(phi1) * sin(phi2)
t_1 = atan2((cos(phi2) * sin(lambda1)), (t_0 - ((cos(phi2) * sin(phi1)) * cos(lambda1))))
if (lambda1 <= (-7.2d+93)) then
tmp = t_1
else if (lambda1 <= 2.6d-7) then
tmp = atan2((sin((lambda1 - lambda2)) * cos(phi2)), (t_0 - (sin(phi1) * cos((lambda1 - lambda2)))))
else
tmp = 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.atan2((Math.cos(phi2) * Math.sin(lambda1)), (t_0 - ((Math.cos(phi2) * Math.sin(phi1)) * Math.cos(lambda1))));
double tmp;
if (lambda1 <= -7.2e+93) {
tmp = t_1;
} else if (lambda1 <= 2.6e-7) {
tmp = Math.atan2((Math.sin((lambda1 - lambda2)) * Math.cos(phi2)), (t_0 - (Math.sin(phi1) * Math.cos((lambda1 - lambda2)))));
} else {
tmp = t_1;
}
return tmp;
}
def code(lambda1, lambda2, phi1, phi2): t_0 = math.cos(phi1) * math.sin(phi2) t_1 = math.atan2((math.cos(phi2) * math.sin(lambda1)), (t_0 - ((math.cos(phi2) * math.sin(phi1)) * math.cos(lambda1)))) tmp = 0 if lambda1 <= -7.2e+93: tmp = t_1 elif lambda1 <= 2.6e-7: tmp = math.atan2((math.sin((lambda1 - lambda2)) * math.cos(phi2)), (t_0 - (math.sin(phi1) * math.cos((lambda1 - lambda2))))) else: tmp = t_1 return tmp
function code(lambda1, lambda2, phi1, phi2) t_0 = Float64(cos(phi1) * sin(phi2)) t_1 = atan(Float64(cos(phi2) * sin(lambda1)), Float64(t_0 - Float64(Float64(cos(phi2) * sin(phi1)) * cos(lambda1)))) tmp = 0.0 if (lambda1 <= -7.2e+93) tmp = t_1; elseif (lambda1 <= 2.6e-7) tmp = atan(Float64(sin(Float64(lambda1 - lambda2)) * cos(phi2)), Float64(t_0 - Float64(sin(phi1) * cos(Float64(lambda1 - lambda2))))); else tmp = t_1; end return tmp end
function tmp_2 = code(lambda1, lambda2, phi1, phi2) t_0 = cos(phi1) * sin(phi2); t_1 = atan2((cos(phi2) * sin(lambda1)), (t_0 - ((cos(phi2) * sin(phi1)) * cos(lambda1)))); tmp = 0.0; if (lambda1 <= -7.2e+93) tmp = t_1; elseif (lambda1 <= 2.6e-7) tmp = atan2((sin((lambda1 - lambda2)) * cos(phi2)), (t_0 - (sin(phi1) * cos((lambda1 - lambda2))))); else tmp = 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[ArcTan[N[(N[Cos[phi2], $MachinePrecision] * N[Sin[lambda1], $MachinePrecision]), $MachinePrecision] / N[(t$95$0 - N[(N[(N[Cos[phi2], $MachinePrecision] * N[Sin[phi1], $MachinePrecision]), $MachinePrecision] * N[Cos[lambda1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]}, If[LessEqual[lambda1, -7.2e+93], t$95$1, If[LessEqual[lambda1, 2.6e-7], N[ArcTan[N[(N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision] / N[(t$95$0 - N[(N[Sin[phi1], $MachinePrecision] * N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \cos \phi_1 \cdot \sin \phi_2\\
t_1 := \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \lambda_1}{t\_0 - \left(\cos \phi_2 \cdot \sin \phi_1\right) \cdot \cos \lambda_1}\\
\mathbf{if}\;\lambda_1 \leq -7.2 \cdot 10^{+93}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;\lambda_1 \leq 2.6 \cdot 10^{-7}:\\
\;\;\;\;\tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2}{t\_0 - \sin \phi_1 \cdot \cos \left(\lambda_1 - \lambda_2\right)}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if lambda1 < -7.1999999999999998e93 or 2.59999999999999999e-7 < lambda1 Initial program 66.7%
Taylor expanded in lambda2 around 0
lower-cos.f6466.7
Simplified66.7%
Taylor expanded in lambda2 around 0
lower-sin.f6466.6
Simplified66.6%
if -7.1999999999999998e93 < lambda1 < 2.59999999999999999e-7Initial program 93.9%
Taylor expanded in phi2 around 0
lower-sin.f6481.5
Simplified81.5%
Final simplification75.3%
(FPCore (lambda1 lambda2 phi1 phi2) :precision binary64 (atan2 (* (sin (- lambda1 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((sin((lambda1 - lambda2)) * cos(phi2)), ((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((sin((lambda1 - lambda2)) * cos(phi2)), ((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.sin((lambda1 - lambda2)) * Math.cos(phi2)), ((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.sin((lambda1 - lambda2)) * math.cos(phi2)), ((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(sin(Float64(lambda1 - lambda2)) * cos(phi2)), 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((sin((lambda1 - lambda2)) * cos(phi2)), ((cos(phi1) * sin(phi2)) - ((cos(phi2) * sin(phi1)) * 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[Cos[phi2], $MachinePrecision] * N[Sin[phi1], $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(\cos \phi_2 \cdot \sin \phi_1\right) \cdot \cos \left(\lambda_1 - \lambda_2\right)}
\end{array}
Initial program 82.6%
Final simplification82.6%
(FPCore (lambda1 lambda2 phi1 phi2) :precision binary64 (atan2 (* (sin (- lambda1 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((sin((lambda1 - lambda2)) * cos(phi2)), ((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((sin((lambda1 - lambda2)) * cos(phi2)), ((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.sin((lambda1 - lambda2)) * Math.cos(phi2)), ((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.sin((lambda1 - lambda2)) * math.cos(phi2)), ((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(sin(Float64(lambda1 - lambda2)) * cos(phi2)), 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((sin((lambda1 - lambda2)) * cos(phi2)), ((cos(phi1) * sin(phi2)) - (cos(phi2) * (sin(phi1) * 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[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{\sin \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2}{\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 82.6%
Taylor expanded in phi1 around inf
lower-*.f64N/A
lower-cos.f64N/A
lower-*.f64N/A
lower-cos.f64N/A
lower--.f64N/A
lower-sin.f6482.6
Simplified82.6%
Final simplification82.6%
(FPCore (lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (* (cos phi1) (sin phi2)))
(t_1 (sin (- lambda1 lambda2)))
(t_2 (atan2 (* t_1 (cos phi2)) (- t_0 (* (cos phi2) (sin phi1))))))
(if (<= phi2 -6.2)
t_2
(if (<= phi2 2.7e-35)
(atan2 t_1 (- t_0 (* (sin phi1) (cos (- lambda1 lambda2)))))
t_2))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = cos(phi1) * sin(phi2);
double t_1 = sin((lambda1 - lambda2));
double t_2 = atan2((t_1 * cos(phi2)), (t_0 - (cos(phi2) * sin(phi1))));
double tmp;
if (phi2 <= -6.2) {
tmp = t_2;
} else if (phi2 <= 2.7e-35) {
tmp = atan2(t_1, (t_0 - (sin(phi1) * cos((lambda1 - lambda2)))));
} else {
tmp = t_2;
}
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 = sin((lambda1 - lambda2))
t_2 = atan2((t_1 * cos(phi2)), (t_0 - (cos(phi2) * sin(phi1))))
if (phi2 <= (-6.2d0)) then
tmp = t_2
else if (phi2 <= 2.7d-35) then
tmp = atan2(t_1, (t_0 - (sin(phi1) * cos((lambda1 - lambda2)))))
else
tmp = t_2
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.sin((lambda1 - lambda2));
double t_2 = Math.atan2((t_1 * Math.cos(phi2)), (t_0 - (Math.cos(phi2) * Math.sin(phi1))));
double tmp;
if (phi2 <= -6.2) {
tmp = t_2;
} else if (phi2 <= 2.7e-35) {
tmp = Math.atan2(t_1, (t_0 - (Math.sin(phi1) * Math.cos((lambda1 - lambda2)))));
} else {
tmp = t_2;
}
return tmp;
}
def code(lambda1, lambda2, phi1, phi2): t_0 = math.cos(phi1) * math.sin(phi2) t_1 = math.sin((lambda1 - lambda2)) t_2 = math.atan2((t_1 * math.cos(phi2)), (t_0 - (math.cos(phi2) * math.sin(phi1)))) tmp = 0 if phi2 <= -6.2: tmp = t_2 elif phi2 <= 2.7e-35: tmp = math.atan2(t_1, (t_0 - (math.sin(phi1) * math.cos((lambda1 - lambda2))))) else: tmp = t_2 return tmp
function code(lambda1, lambda2, phi1, phi2) t_0 = Float64(cos(phi1) * sin(phi2)) t_1 = sin(Float64(lambda1 - lambda2)) t_2 = atan(Float64(t_1 * cos(phi2)), Float64(t_0 - Float64(cos(phi2) * sin(phi1)))) tmp = 0.0 if (phi2 <= -6.2) tmp = t_2; elseif (phi2 <= 2.7e-35) tmp = atan(t_1, Float64(t_0 - Float64(sin(phi1) * cos(Float64(lambda1 - lambda2))))); else tmp = t_2; end return tmp end
function tmp_2 = code(lambda1, lambda2, phi1, phi2) t_0 = cos(phi1) * sin(phi2); t_1 = sin((lambda1 - lambda2)); t_2 = atan2((t_1 * cos(phi2)), (t_0 - (cos(phi2) * sin(phi1)))); tmp = 0.0; if (phi2 <= -6.2) tmp = t_2; elseif (phi2 <= 2.7e-35) tmp = atan2(t_1, (t_0 - (sin(phi1) * cos((lambda1 - lambda2))))); else tmp = t_2; 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[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$2 = N[ArcTan[N[(t$95$1 * N[Cos[phi2], $MachinePrecision]), $MachinePrecision] / N[(t$95$0 - N[(N[Cos[phi2], $MachinePrecision] * N[Sin[phi1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]}, If[LessEqual[phi2, -6.2], t$95$2, If[LessEqual[phi2, 2.7e-35], N[ArcTan[t$95$1 / N[(t$95$0 - N[(N[Sin[phi1], $MachinePrecision] * N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], t$95$2]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \cos \phi_1 \cdot \sin \phi_2\\
t_1 := \sin \left(\lambda_1 - \lambda_2\right)\\
t_2 := \tan^{-1}_* \frac{t\_1 \cdot \cos \phi_2}{t\_0 - \cos \phi_2 \cdot \sin \phi_1}\\
\mathbf{if}\;\phi_2 \leq -6.2:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;\phi_2 \leq 2.7 \cdot 10^{-35}:\\
\;\;\;\;\tan^{-1}_* \frac{t\_1}{t\_0 - \sin \phi_1 \cdot \cos \left(\lambda_1 - \lambda_2\right)}\\
\mathbf{else}:\\
\;\;\;\;t\_2\\
\end{array}
\end{array}
if phi2 < -6.20000000000000018 or 2.6999999999999997e-35 < phi2 Initial program 78.8%
Taylor expanded in lambda1 around 0
cos-negN/A
lower-cos.f6468.8
Simplified68.8%
Taylor expanded in lambda2 around 0
*-commutativeN/A
lower-*.f64N/A
lower-sin.f64N/A
lower-cos.f6457.7
Simplified57.7%
if -6.20000000000000018 < phi2 < 2.6999999999999997e-35Initial program 86.8%
Taylor expanded in phi2 around 0
lower-sin.f6486.3
Simplified86.3%
Taylor expanded in phi2 around 0
lower-sin.f64N/A
lower--.f6486.3
Simplified86.3%
Final simplification71.6%
(FPCore (lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (* (sin (- lambda1 lambda2)) (cos phi2)))
(t_1 (* (cos phi1) (sin phi2))))
(if (<= lambda1 -1e-8)
(atan2 t_0 (- (sin phi2) (* (* (cos phi2) (sin phi1)) (cos lambda1))))
(if (<= lambda1 2.6e-7)
(atan2 t_0 (- t_1 (* (sin phi1) (cos lambda2))))
(atan2
(* (cos phi2) (sin lambda1))
(- t_1 (* (sin phi1) (cos lambda1))))))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = sin((lambda1 - lambda2)) * cos(phi2);
double t_1 = cos(phi1) * sin(phi2);
double tmp;
if (lambda1 <= -1e-8) {
tmp = atan2(t_0, (sin(phi2) - ((cos(phi2) * sin(phi1)) * cos(lambda1))));
} else if (lambda1 <= 2.6e-7) {
tmp = atan2(t_0, (t_1 - (sin(phi1) * cos(lambda2))));
} else {
tmp = atan2((cos(phi2) * sin(lambda1)), (t_1 - (sin(phi1) * cos(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) :: t_0
real(8) :: t_1
real(8) :: tmp
t_0 = sin((lambda1 - lambda2)) * cos(phi2)
t_1 = cos(phi1) * sin(phi2)
if (lambda1 <= (-1d-8)) then
tmp = atan2(t_0, (sin(phi2) - ((cos(phi2) * sin(phi1)) * cos(lambda1))))
else if (lambda1 <= 2.6d-7) then
tmp = atan2(t_0, (t_1 - (sin(phi1) * cos(lambda2))))
else
tmp = atan2((cos(phi2) * sin(lambda1)), (t_1 - (sin(phi1) * cos(lambda1))))
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)) * Math.cos(phi2);
double t_1 = Math.cos(phi1) * Math.sin(phi2);
double tmp;
if (lambda1 <= -1e-8) {
tmp = Math.atan2(t_0, (Math.sin(phi2) - ((Math.cos(phi2) * Math.sin(phi1)) * Math.cos(lambda1))));
} else if (lambda1 <= 2.6e-7) {
tmp = Math.atan2(t_0, (t_1 - (Math.sin(phi1) * Math.cos(lambda2))));
} else {
tmp = Math.atan2((Math.cos(phi2) * Math.sin(lambda1)), (t_1 - (Math.sin(phi1) * Math.cos(lambda1))));
}
return tmp;
}
def code(lambda1, lambda2, phi1, phi2): t_0 = math.sin((lambda1 - lambda2)) * math.cos(phi2) t_1 = math.cos(phi1) * math.sin(phi2) tmp = 0 if lambda1 <= -1e-8: tmp = math.atan2(t_0, (math.sin(phi2) - ((math.cos(phi2) * math.sin(phi1)) * math.cos(lambda1)))) elif lambda1 <= 2.6e-7: tmp = math.atan2(t_0, (t_1 - (math.sin(phi1) * math.cos(lambda2)))) else: tmp = math.atan2((math.cos(phi2) * math.sin(lambda1)), (t_1 - (math.sin(phi1) * math.cos(lambda1)))) return tmp
function code(lambda1, lambda2, phi1, phi2) t_0 = Float64(sin(Float64(lambda1 - lambda2)) * cos(phi2)) t_1 = Float64(cos(phi1) * sin(phi2)) tmp = 0.0 if (lambda1 <= -1e-8) tmp = atan(t_0, Float64(sin(phi2) - Float64(Float64(cos(phi2) * sin(phi1)) * cos(lambda1)))); elseif (lambda1 <= 2.6e-7) tmp = atan(t_0, Float64(t_1 - Float64(sin(phi1) * cos(lambda2)))); else tmp = atan(Float64(cos(phi2) * sin(lambda1)), Float64(t_1 - Float64(sin(phi1) * cos(lambda1)))); end return tmp end
function tmp_2 = code(lambda1, lambda2, phi1, phi2) t_0 = sin((lambda1 - lambda2)) * cos(phi2); t_1 = cos(phi1) * sin(phi2); tmp = 0.0; if (lambda1 <= -1e-8) tmp = atan2(t_0, (sin(phi2) - ((cos(phi2) * sin(phi1)) * cos(lambda1)))); elseif (lambda1 <= 2.6e-7) tmp = atan2(t_0, (t_1 - (sin(phi1) * cos(lambda2)))); else tmp = atan2((cos(phi2) * sin(lambda1)), (t_1 - (sin(phi1) * cos(lambda1)))); end tmp_2 = tmp; end
code[lambda1_, lambda2_, phi1_, phi2_] := Block[{t$95$0 = N[(N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[Cos[phi1], $MachinePrecision] * N[Sin[phi2], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[lambda1, -1e-8], N[ArcTan[t$95$0 / N[(N[Sin[phi2], $MachinePrecision] - N[(N[(N[Cos[phi2], $MachinePrecision] * N[Sin[phi1], $MachinePrecision]), $MachinePrecision] * N[Cos[lambda1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], If[LessEqual[lambda1, 2.6e-7], N[ArcTan[t$95$0 / N[(t$95$1 - N[(N[Sin[phi1], $MachinePrecision] * N[Cos[lambda2], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], N[ArcTan[N[(N[Cos[phi2], $MachinePrecision] * N[Sin[lambda1], $MachinePrecision]), $MachinePrecision] / N[(t$95$1 - N[(N[Sin[phi1], $MachinePrecision] * N[Cos[lambda1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sin \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2\\
t_1 := \cos \phi_1 \cdot \sin \phi_2\\
\mathbf{if}\;\lambda_1 \leq -1 \cdot 10^{-8}:\\
\;\;\;\;\tan^{-1}_* \frac{t\_0}{\sin \phi_2 - \left(\cos \phi_2 \cdot \sin \phi_1\right) \cdot \cos \lambda_1}\\
\mathbf{elif}\;\lambda_1 \leq 2.6 \cdot 10^{-7}:\\
\;\;\;\;\tan^{-1}_* \frac{t\_0}{t\_1 - \sin \phi_1 \cdot \cos \lambda_2}\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \lambda_1}{t\_1 - \sin \phi_1 \cdot \cos \lambda_1}\\
\end{array}
\end{array}
if lambda1 < -1e-8Initial program 66.9%
Taylor expanded in lambda2 around 0
lower-cos.f6466.3
Simplified66.3%
Taylor expanded in phi1 around 0
lower-sin.f6454.8
Simplified54.8%
if -1e-8 < lambda1 < 2.59999999999999999e-7Initial program 99.2%
Taylor expanded in phi2 around 0
lower-sin.f6486.1
Simplified86.1%
Taylor expanded in lambda1 around 0
cos-negN/A
lower-cos.f6486.1
Simplified86.1%
if 2.59999999999999999e-7 < lambda1 Initial program 62.7%
Taylor expanded in phi2 around 0
lower-sin.f6450.6
Simplified50.6%
Taylor expanded in lambda2 around 0
lower-sin.f6450.8
Simplified50.8%
Taylor expanded in lambda2 around 0
lower-cos.f6450.8
Simplified50.8%
Final simplification70.1%
(FPCore (lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (* (sin (- lambda1 lambda2)) (cos phi2)))
(t_1
(atan2
t_0
(- (sin phi2) (* (cos phi2) (* (sin phi1) (cos lambda2)))))))
(if (<= lambda2 -750.0)
t_1
(if (<= lambda2 1.85e+45)
(atan2 t_0 (- (* (cos phi1) (sin phi2)) (* (sin phi1) (cos lambda1))))
t_1))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = sin((lambda1 - lambda2)) * cos(phi2);
double t_1 = atan2(t_0, (sin(phi2) - (cos(phi2) * (sin(phi1) * cos(lambda2)))));
double tmp;
if (lambda2 <= -750.0) {
tmp = t_1;
} else if (lambda2 <= 1.85e+45) {
tmp = atan2(t_0, ((cos(phi1) * sin(phi2)) - (sin(phi1) * cos(lambda1))));
} else {
tmp = 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) :: tmp
t_0 = sin((lambda1 - lambda2)) * cos(phi2)
t_1 = atan2(t_0, (sin(phi2) - (cos(phi2) * (sin(phi1) * cos(lambda2)))))
if (lambda2 <= (-750.0d0)) then
tmp = t_1
else if (lambda2 <= 1.85d+45) then
tmp = atan2(t_0, ((cos(phi1) * sin(phi2)) - (sin(phi1) * cos(lambda1))))
else
tmp = t_1
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)) * Math.cos(phi2);
double t_1 = Math.atan2(t_0, (Math.sin(phi2) - (Math.cos(phi2) * (Math.sin(phi1) * Math.cos(lambda2)))));
double tmp;
if (lambda2 <= -750.0) {
tmp = t_1;
} else if (lambda2 <= 1.85e+45) {
tmp = Math.atan2(t_0, ((Math.cos(phi1) * Math.sin(phi2)) - (Math.sin(phi1) * Math.cos(lambda1))));
} else {
tmp = t_1;
}
return tmp;
}
def code(lambda1, lambda2, phi1, phi2): t_0 = math.sin((lambda1 - lambda2)) * math.cos(phi2) t_1 = math.atan2(t_0, (math.sin(phi2) - (math.cos(phi2) * (math.sin(phi1) * math.cos(lambda2))))) tmp = 0 if lambda2 <= -750.0: tmp = t_1 elif lambda2 <= 1.85e+45: tmp = math.atan2(t_0, ((math.cos(phi1) * math.sin(phi2)) - (math.sin(phi1) * math.cos(lambda1)))) else: tmp = t_1 return tmp
function code(lambda1, lambda2, phi1, phi2) t_0 = Float64(sin(Float64(lambda1 - lambda2)) * cos(phi2)) t_1 = atan(t_0, Float64(sin(phi2) - Float64(cos(phi2) * Float64(sin(phi1) * cos(lambda2))))) tmp = 0.0 if (lambda2 <= -750.0) tmp = t_1; elseif (lambda2 <= 1.85e+45) tmp = atan(t_0, Float64(Float64(cos(phi1) * sin(phi2)) - Float64(sin(phi1) * cos(lambda1)))); else tmp = t_1; end return tmp end
function tmp_2 = code(lambda1, lambda2, phi1, phi2) t_0 = sin((lambda1 - lambda2)) * cos(phi2); t_1 = atan2(t_0, (sin(phi2) - (cos(phi2) * (sin(phi1) * cos(lambda2))))); tmp = 0.0; if (lambda2 <= -750.0) tmp = t_1; elseif (lambda2 <= 1.85e+45) tmp = atan2(t_0, ((cos(phi1) * sin(phi2)) - (sin(phi1) * cos(lambda1)))); else tmp = t_1; end tmp_2 = tmp; end
code[lambda1_, lambda2_, phi1_, phi2_] := Block[{t$95$0 = N[(N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[ArcTan[t$95$0 / N[(N[Sin[phi2], $MachinePrecision] - N[(N[Cos[phi2], $MachinePrecision] * N[(N[Sin[phi1], $MachinePrecision] * N[Cos[lambda2], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]}, If[LessEqual[lambda2, -750.0], t$95$1, If[LessEqual[lambda2, 1.85e+45], N[ArcTan[t$95$0 / N[(N[(N[Cos[phi1], $MachinePrecision] * N[Sin[phi2], $MachinePrecision]), $MachinePrecision] - N[(N[Sin[phi1], $MachinePrecision] * N[Cos[lambda1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sin \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2\\
t_1 := \tan^{-1}_* \frac{t\_0}{\sin \phi_2 - \cos \phi_2 \cdot \left(\sin \phi_1 \cdot \cos \lambda_2\right)}\\
\mathbf{if}\;\lambda_2 \leq -750:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;\lambda_2 \leq 1.85 \cdot 10^{+45}:\\
\;\;\;\;\tan^{-1}_* \frac{t\_0}{\cos \phi_1 \cdot \sin \phi_2 - \sin \phi_1 \cdot \cos \lambda_1}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if lambda2 < -750 or 1.84999999999999989e45 < lambda2 Initial program 68.5%
Taylor expanded in phi1 around inf
lower-*.f64N/A
lower-cos.f64N/A
lower-*.f64N/A
lower-cos.f64N/A
lower--.f64N/A
lower-sin.f6468.5
Simplified68.5%
Taylor expanded in lambda1 around 0
cos-negN/A
lower-cos.f6468.5
Simplified68.5%
Taylor expanded in phi1 around 0
lower-sin.f6458.3
Simplified58.3%
if -750 < lambda2 < 1.84999999999999989e45Initial program 95.8%
Taylor expanded in phi2 around 0
lower-sin.f6480.3
Simplified80.3%
Taylor expanded in lambda2 around 0
lower-cos.f6480.4
Simplified80.4%
Final simplification69.8%
(FPCore (lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (sin (- lambda1 lambda2)))
(t_1 (* t_0 (cos phi2)))
(t_2 (* (cos phi1) (sin phi2))))
(if (<= phi2 -40000000000.0)
(atan2 t_1 (- (sin phi2) (* (cos phi2) (* (sin phi1) (cos lambda2)))))
(if (<= phi2 2.7e-35)
(atan2 t_0 (- t_2 (* (sin phi1) (cos (- lambda1 lambda2)))))
(atan2 t_1 (- t_2 (sin phi1)))))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = sin((lambda1 - lambda2));
double t_1 = t_0 * cos(phi2);
double t_2 = cos(phi1) * sin(phi2);
double tmp;
if (phi2 <= -40000000000.0) {
tmp = atan2(t_1, (sin(phi2) - (cos(phi2) * (sin(phi1) * cos(lambda2)))));
} else if (phi2 <= 2.7e-35) {
tmp = atan2(t_0, (t_2 - (sin(phi1) * cos((lambda1 - lambda2)))));
} else {
tmp = atan2(t_1, (t_2 - sin(phi1)));
}
return tmp;
}
real(8) function code(lambda1, lambda2, phi1, phi2)
real(8), intent (in) :: lambda1
real(8), intent (in) :: lambda2
real(8), intent (in) :: phi1
real(8), intent (in) :: phi2
real(8) :: t_0
real(8) :: t_1
real(8) :: t_2
real(8) :: tmp
t_0 = sin((lambda1 - lambda2))
t_1 = t_0 * cos(phi2)
t_2 = cos(phi1) * sin(phi2)
if (phi2 <= (-40000000000.0d0)) then
tmp = atan2(t_1, (sin(phi2) - (cos(phi2) * (sin(phi1) * cos(lambda2)))))
else if (phi2 <= 2.7d-35) then
tmp = atan2(t_0, (t_2 - (sin(phi1) * cos((lambda1 - lambda2)))))
else
tmp = atan2(t_1, (t_2 - sin(phi1)))
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 = t_0 * Math.cos(phi2);
double t_2 = Math.cos(phi1) * Math.sin(phi2);
double tmp;
if (phi2 <= -40000000000.0) {
tmp = Math.atan2(t_1, (Math.sin(phi2) - (Math.cos(phi2) * (Math.sin(phi1) * Math.cos(lambda2)))));
} else if (phi2 <= 2.7e-35) {
tmp = Math.atan2(t_0, (t_2 - (Math.sin(phi1) * Math.cos((lambda1 - lambda2)))));
} else {
tmp = Math.atan2(t_1, (t_2 - Math.sin(phi1)));
}
return tmp;
}
def code(lambda1, lambda2, phi1, phi2): t_0 = math.sin((lambda1 - lambda2)) t_1 = t_0 * math.cos(phi2) t_2 = math.cos(phi1) * math.sin(phi2) tmp = 0 if phi2 <= -40000000000.0: tmp = math.atan2(t_1, (math.sin(phi2) - (math.cos(phi2) * (math.sin(phi1) * math.cos(lambda2))))) elif phi2 <= 2.7e-35: tmp = math.atan2(t_0, (t_2 - (math.sin(phi1) * math.cos((lambda1 - lambda2))))) else: tmp = math.atan2(t_1, (t_2 - math.sin(phi1))) return tmp
function code(lambda1, lambda2, phi1, phi2) t_0 = sin(Float64(lambda1 - lambda2)) t_1 = Float64(t_0 * cos(phi2)) t_2 = Float64(cos(phi1) * sin(phi2)) tmp = 0.0 if (phi2 <= -40000000000.0) tmp = atan(t_1, Float64(sin(phi2) - Float64(cos(phi2) * Float64(sin(phi1) * cos(lambda2))))); elseif (phi2 <= 2.7e-35) tmp = atan(t_0, Float64(t_2 - Float64(sin(phi1) * cos(Float64(lambda1 - lambda2))))); else tmp = atan(t_1, Float64(t_2 - sin(phi1))); end return tmp end
function tmp_2 = code(lambda1, lambda2, phi1, phi2) t_0 = sin((lambda1 - lambda2)); t_1 = t_0 * cos(phi2); t_2 = cos(phi1) * sin(phi2); tmp = 0.0; if (phi2 <= -40000000000.0) tmp = atan2(t_1, (sin(phi2) - (cos(phi2) * (sin(phi1) * cos(lambda2))))); elseif (phi2 <= 2.7e-35) tmp = atan2(t_0, (t_2 - (sin(phi1) * cos((lambda1 - lambda2))))); else tmp = atan2(t_1, (t_2 - sin(phi1))); 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[(t$95$0 * N[Cos[phi2], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(N[Cos[phi1], $MachinePrecision] * N[Sin[phi2], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[phi2, -40000000000.0], N[ArcTan[t$95$1 / N[(N[Sin[phi2], $MachinePrecision] - N[(N[Cos[phi2], $MachinePrecision] * N[(N[Sin[phi1], $MachinePrecision] * N[Cos[lambda2], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], If[LessEqual[phi2, 2.7e-35], N[ArcTan[t$95$0 / N[(t$95$2 - N[(N[Sin[phi1], $MachinePrecision] * N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], N[ArcTan[t$95$1 / N[(t$95$2 - N[Sin[phi1], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sin \left(\lambda_1 - \lambda_2\right)\\
t_1 := t\_0 \cdot \cos \phi_2\\
t_2 := \cos \phi_1 \cdot \sin \phi_2\\
\mathbf{if}\;\phi_2 \leq -40000000000:\\
\;\;\;\;\tan^{-1}_* \frac{t\_1}{\sin \phi_2 - \cos \phi_2 \cdot \left(\sin \phi_1 \cdot \cos \lambda_2\right)}\\
\mathbf{elif}\;\phi_2 \leq 2.7 \cdot 10^{-35}:\\
\;\;\;\;\tan^{-1}_* \frac{t\_0}{t\_2 - \sin \phi_1 \cdot \cos \left(\lambda_1 - \lambda_2\right)}\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1}_* \frac{t\_1}{t\_2 - \sin \phi_1}\\
\end{array}
\end{array}
if phi2 < -4e10Initial program 83.0%
Taylor expanded in phi1 around inf
lower-*.f64N/A
lower-cos.f64N/A
lower-*.f64N/A
lower-cos.f64N/A
lower--.f64N/A
lower-sin.f6483.0
Simplified83.0%
Taylor expanded in lambda1 around 0
cos-negN/A
lower-cos.f6473.9
Simplified73.9%
Taylor expanded in phi1 around 0
lower-sin.f6461.8
Simplified61.8%
if -4e10 < phi2 < 2.6999999999999997e-35Initial program 87.0%
Taylor expanded in phi2 around 0
lower-sin.f6485.3
Simplified85.3%
Taylor expanded in phi2 around 0
lower-sin.f64N/A
lower--.f6485.3
Simplified85.3%
if 2.6999999999999997e-35 < phi2 Initial program 74.9%
Taylor expanded in phi2 around 0
lower-sin.f6449.6
Simplified49.6%
Taylor expanded in lambda2 around 0
lower-cos.f6449.3
Simplified49.3%
Taylor expanded in lambda1 around 0
lower-sin.f6448.8
Simplified48.8%
Final simplification69.6%
(FPCore (lambda1 lambda2 phi1 phi2) :precision binary64 (atan2 (* (sin (- lambda1 lambda2)) (cos phi2)) (- (* (cos phi1) (sin phi2)) (* (sin phi1) (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((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((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((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((lambda1 - lambda2)))))
function code(lambda1, lambda2, phi1, phi2) return atan(Float64(sin(Float64(lambda1 - lambda2)) * cos(phi2)), Float64(Float64(cos(phi1) * sin(phi2)) - Float64(sin(phi1) * 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((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[Sin[phi1], $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 - \sin \phi_1 \cdot \cos \left(\lambda_1 - \lambda_2\right)}
\end{array}
Initial program 82.6%
Taylor expanded in phi2 around 0
lower-sin.f6469.9
Simplified69.9%
(FPCore (lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (* (cos phi1) (sin phi2)))
(t_1 (sin (- lambda1 lambda2)))
(t_2 (atan2 (* t_1 (cos phi2)) (- t_0 (sin phi1)))))
(if (<= phi2 -6.2)
t_2
(if (<= phi2 2.7e-35)
(atan2 t_1 (- t_0 (* (sin phi1) (cos (- lambda1 lambda2)))))
t_2))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = cos(phi1) * sin(phi2);
double t_1 = sin((lambda1 - lambda2));
double t_2 = atan2((t_1 * cos(phi2)), (t_0 - sin(phi1)));
double tmp;
if (phi2 <= -6.2) {
tmp = t_2;
} else if (phi2 <= 2.7e-35) {
tmp = atan2(t_1, (t_0 - (sin(phi1) * cos((lambda1 - lambda2)))));
} else {
tmp = t_2;
}
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 = sin((lambda1 - lambda2))
t_2 = atan2((t_1 * cos(phi2)), (t_0 - sin(phi1)))
if (phi2 <= (-6.2d0)) then
tmp = t_2
else if (phi2 <= 2.7d-35) then
tmp = atan2(t_1, (t_0 - (sin(phi1) * cos((lambda1 - lambda2)))))
else
tmp = t_2
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.sin((lambda1 - lambda2));
double t_2 = Math.atan2((t_1 * Math.cos(phi2)), (t_0 - Math.sin(phi1)));
double tmp;
if (phi2 <= -6.2) {
tmp = t_2;
} else if (phi2 <= 2.7e-35) {
tmp = Math.atan2(t_1, (t_0 - (Math.sin(phi1) * Math.cos((lambda1 - lambda2)))));
} else {
tmp = t_2;
}
return tmp;
}
def code(lambda1, lambda2, phi1, phi2): t_0 = math.cos(phi1) * math.sin(phi2) t_1 = math.sin((lambda1 - lambda2)) t_2 = math.atan2((t_1 * math.cos(phi2)), (t_0 - math.sin(phi1))) tmp = 0 if phi2 <= -6.2: tmp = t_2 elif phi2 <= 2.7e-35: tmp = math.atan2(t_1, (t_0 - (math.sin(phi1) * math.cos((lambda1 - lambda2))))) else: tmp = t_2 return tmp
function code(lambda1, lambda2, phi1, phi2) t_0 = Float64(cos(phi1) * sin(phi2)) t_1 = sin(Float64(lambda1 - lambda2)) t_2 = atan(Float64(t_1 * cos(phi2)), Float64(t_0 - sin(phi1))) tmp = 0.0 if (phi2 <= -6.2) tmp = t_2; elseif (phi2 <= 2.7e-35) tmp = atan(t_1, Float64(t_0 - Float64(sin(phi1) * cos(Float64(lambda1 - lambda2))))); else tmp = t_2; end return tmp end
function tmp_2 = code(lambda1, lambda2, phi1, phi2) t_0 = cos(phi1) * sin(phi2); t_1 = sin((lambda1 - lambda2)); t_2 = atan2((t_1 * cos(phi2)), (t_0 - sin(phi1))); tmp = 0.0; if (phi2 <= -6.2) tmp = t_2; elseif (phi2 <= 2.7e-35) tmp = atan2(t_1, (t_0 - (sin(phi1) * cos((lambda1 - lambda2))))); else tmp = t_2; 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[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$2 = N[ArcTan[N[(t$95$1 * N[Cos[phi2], $MachinePrecision]), $MachinePrecision] / N[(t$95$0 - N[Sin[phi1], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]}, If[LessEqual[phi2, -6.2], t$95$2, If[LessEqual[phi2, 2.7e-35], N[ArcTan[t$95$1 / N[(t$95$0 - N[(N[Sin[phi1], $MachinePrecision] * N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], t$95$2]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \cos \phi_1 \cdot \sin \phi_2\\
t_1 := \sin \left(\lambda_1 - \lambda_2\right)\\
t_2 := \tan^{-1}_* \frac{t\_1 \cdot \cos \phi_2}{t\_0 - \sin \phi_1}\\
\mathbf{if}\;\phi_2 \leq -6.2:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;\phi_2 \leq 2.7 \cdot 10^{-35}:\\
\;\;\;\;\tan^{-1}_* \frac{t\_1}{t\_0 - \sin \phi_1 \cdot \cos \left(\lambda_1 - \lambda_2\right)}\\
\mathbf{else}:\\
\;\;\;\;t\_2\\
\end{array}
\end{array}
if phi2 < -6.20000000000000018 or 2.6999999999999997e-35 < phi2 Initial program 78.8%
Taylor expanded in phi2 around 0
lower-sin.f6454.5
Simplified54.5%
Taylor expanded in lambda2 around 0
lower-cos.f6454.2
Simplified54.2%
Taylor expanded in lambda1 around 0
lower-sin.f6453.8
Simplified53.8%
if -6.20000000000000018 < phi2 < 2.6999999999999997e-35Initial program 86.8%
Taylor expanded in phi2 around 0
lower-sin.f6486.3
Simplified86.3%
Taylor expanded in phi2 around 0
lower-sin.f64N/A
lower--.f6486.3
Simplified86.3%
(FPCore (lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (sin (- lambda1 lambda2)))
(t_1
(atan2 (* t_0 (cos phi2)) (- (* (cos phi1) (sin phi2)) (sin phi1)))))
(if (<= phi2 -0.0027)
t_1
(if (<= phi2 2.7e-35)
(atan2
t_0
(- (* phi2 (cos phi1)) (* (sin phi1) (cos (- lambda1 lambda2)))))
t_1))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = sin((lambda1 - lambda2));
double t_1 = atan2((t_0 * cos(phi2)), ((cos(phi1) * sin(phi2)) - sin(phi1)));
double tmp;
if (phi2 <= -0.0027) {
tmp = t_1;
} else if (phi2 <= 2.7e-35) {
tmp = atan2(t_0, ((phi2 * cos(phi1)) - (sin(phi1) * cos((lambda1 - lambda2)))));
} else {
tmp = 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) :: tmp
t_0 = sin((lambda1 - lambda2))
t_1 = atan2((t_0 * cos(phi2)), ((cos(phi1) * sin(phi2)) - sin(phi1)))
if (phi2 <= (-0.0027d0)) then
tmp = t_1
else if (phi2 <= 2.7d-35) then
tmp = atan2(t_0, ((phi2 * cos(phi1)) - (sin(phi1) * cos((lambda1 - lambda2)))))
else
tmp = t_1
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.atan2((t_0 * Math.cos(phi2)), ((Math.cos(phi1) * Math.sin(phi2)) - Math.sin(phi1)));
double tmp;
if (phi2 <= -0.0027) {
tmp = t_1;
} else if (phi2 <= 2.7e-35) {
tmp = Math.atan2(t_0, ((phi2 * Math.cos(phi1)) - (Math.sin(phi1) * Math.cos((lambda1 - lambda2)))));
} else {
tmp = t_1;
}
return tmp;
}
def code(lambda1, lambda2, phi1, phi2): t_0 = math.sin((lambda1 - lambda2)) t_1 = math.atan2((t_0 * math.cos(phi2)), ((math.cos(phi1) * math.sin(phi2)) - math.sin(phi1))) tmp = 0 if phi2 <= -0.0027: tmp = t_1 elif phi2 <= 2.7e-35: tmp = math.atan2(t_0, ((phi2 * math.cos(phi1)) - (math.sin(phi1) * math.cos((lambda1 - lambda2))))) else: tmp = t_1 return tmp
function code(lambda1, lambda2, phi1, phi2) t_0 = sin(Float64(lambda1 - lambda2)) t_1 = atan(Float64(t_0 * cos(phi2)), Float64(Float64(cos(phi1) * sin(phi2)) - sin(phi1))) tmp = 0.0 if (phi2 <= -0.0027) tmp = t_1; elseif (phi2 <= 2.7e-35) tmp = atan(t_0, Float64(Float64(phi2 * cos(phi1)) - Float64(sin(phi1) * cos(Float64(lambda1 - lambda2))))); else tmp = t_1; end return tmp end
function tmp_2 = code(lambda1, lambda2, phi1, phi2) t_0 = sin((lambda1 - lambda2)); t_1 = atan2((t_0 * cos(phi2)), ((cos(phi1) * sin(phi2)) - sin(phi1))); tmp = 0.0; if (phi2 <= -0.0027) tmp = t_1; elseif (phi2 <= 2.7e-35) tmp = atan2(t_0, ((phi2 * cos(phi1)) - (sin(phi1) * cos((lambda1 - lambda2))))); else tmp = t_1; 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[ArcTan[N[(t$95$0 * N[Cos[phi2], $MachinePrecision]), $MachinePrecision] / N[(N[(N[Cos[phi1], $MachinePrecision] * N[Sin[phi2], $MachinePrecision]), $MachinePrecision] - N[Sin[phi1], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]}, If[LessEqual[phi2, -0.0027], t$95$1, If[LessEqual[phi2, 2.7e-35], N[ArcTan[t$95$0 / N[(N[(phi2 * N[Cos[phi1], $MachinePrecision]), $MachinePrecision] - N[(N[Sin[phi1], $MachinePrecision] * N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sin \left(\lambda_1 - \lambda_2\right)\\
t_1 := \tan^{-1}_* \frac{t\_0 \cdot \cos \phi_2}{\cos \phi_1 \cdot \sin \phi_2 - \sin \phi_1}\\
\mathbf{if}\;\phi_2 \leq -0.0027:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;\phi_2 \leq 2.7 \cdot 10^{-35}:\\
\;\;\;\;\tan^{-1}_* \frac{t\_0}{\phi_2 \cdot \cos \phi_1 - \sin \phi_1 \cdot \cos \left(\lambda_1 - \lambda_2\right)}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if phi2 < -0.0027000000000000001 or 2.6999999999999997e-35 < phi2 Initial program 78.8%
Taylor expanded in phi2 around 0
lower-sin.f6454.3
Simplified54.3%
Taylor expanded in lambda2 around 0
lower-cos.f6454.0
Simplified54.0%
Taylor expanded in lambda1 around 0
lower-sin.f6453.6
Simplified53.6%
if -0.0027000000000000001 < phi2 < 2.6999999999999997e-35Initial program 86.8%
Taylor expanded in phi2 around 0
lower-sin.f6486.7
Simplified86.7%
Taylor expanded in phi2 around 0
lower-sin.f64N/A
lower--.f6486.8
Simplified86.8%
Taylor expanded in phi2 around 0
lower-*.f64N/A
lower-cos.f6486.8
Simplified86.8%
(FPCore (lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (sin (- lambda1 lambda2)))
(t_1
(atan2
(* t_0 (cos phi2))
(- (sin phi2) (* (sin phi1) (cos lambda1))))))
(if (<= phi2 -0.0027)
t_1
(if (<= phi2 4e-69)
(atan2
t_0
(- (* phi2 (cos phi1)) (* (sin phi1) (cos (- lambda1 lambda2)))))
t_1))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = sin((lambda1 - lambda2));
double t_1 = atan2((t_0 * cos(phi2)), (sin(phi2) - (sin(phi1) * cos(lambda1))));
double tmp;
if (phi2 <= -0.0027) {
tmp = t_1;
} else if (phi2 <= 4e-69) {
tmp = atan2(t_0, ((phi2 * cos(phi1)) - (sin(phi1) * cos((lambda1 - lambda2)))));
} else {
tmp = 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) :: tmp
t_0 = sin((lambda1 - lambda2))
t_1 = atan2((t_0 * cos(phi2)), (sin(phi2) - (sin(phi1) * cos(lambda1))))
if (phi2 <= (-0.0027d0)) then
tmp = t_1
else if (phi2 <= 4d-69) then
tmp = atan2(t_0, ((phi2 * cos(phi1)) - (sin(phi1) * cos((lambda1 - lambda2)))))
else
tmp = t_1
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.atan2((t_0 * Math.cos(phi2)), (Math.sin(phi2) - (Math.sin(phi1) * Math.cos(lambda1))));
double tmp;
if (phi2 <= -0.0027) {
tmp = t_1;
} else if (phi2 <= 4e-69) {
tmp = Math.atan2(t_0, ((phi2 * Math.cos(phi1)) - (Math.sin(phi1) * Math.cos((lambda1 - lambda2)))));
} else {
tmp = t_1;
}
return tmp;
}
def code(lambda1, lambda2, phi1, phi2): t_0 = math.sin((lambda1 - lambda2)) t_1 = math.atan2((t_0 * math.cos(phi2)), (math.sin(phi2) - (math.sin(phi1) * math.cos(lambda1)))) tmp = 0 if phi2 <= -0.0027: tmp = t_1 elif phi2 <= 4e-69: tmp = math.atan2(t_0, ((phi2 * math.cos(phi1)) - (math.sin(phi1) * math.cos((lambda1 - lambda2))))) else: tmp = t_1 return tmp
function code(lambda1, lambda2, phi1, phi2) t_0 = sin(Float64(lambda1 - lambda2)) t_1 = atan(Float64(t_0 * cos(phi2)), Float64(sin(phi2) - Float64(sin(phi1) * cos(lambda1)))) tmp = 0.0 if (phi2 <= -0.0027) tmp = t_1; elseif (phi2 <= 4e-69) tmp = atan(t_0, Float64(Float64(phi2 * cos(phi1)) - Float64(sin(phi1) * cos(Float64(lambda1 - lambda2))))); else tmp = t_1; end return tmp end
function tmp_2 = code(lambda1, lambda2, phi1, phi2) t_0 = sin((lambda1 - lambda2)); t_1 = atan2((t_0 * cos(phi2)), (sin(phi2) - (sin(phi1) * cos(lambda1)))); tmp = 0.0; if (phi2 <= -0.0027) tmp = t_1; elseif (phi2 <= 4e-69) tmp = atan2(t_0, ((phi2 * cos(phi1)) - (sin(phi1) * cos((lambda1 - lambda2))))); else tmp = t_1; 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[ArcTan[N[(t$95$0 * N[Cos[phi2], $MachinePrecision]), $MachinePrecision] / N[(N[Sin[phi2], $MachinePrecision] - N[(N[Sin[phi1], $MachinePrecision] * N[Cos[lambda1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]}, If[LessEqual[phi2, -0.0027], t$95$1, If[LessEqual[phi2, 4e-69], N[ArcTan[t$95$0 / N[(N[(phi2 * N[Cos[phi1], $MachinePrecision]), $MachinePrecision] - N[(N[Sin[phi1], $MachinePrecision] * N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sin \left(\lambda_1 - \lambda_2\right)\\
t_1 := \tan^{-1}_* \frac{t\_0 \cdot \cos \phi_2}{\sin \phi_2 - \sin \phi_1 \cdot \cos \lambda_1}\\
\mathbf{if}\;\phi_2 \leq -0.0027:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;\phi_2 \leq 4 \cdot 10^{-69}:\\
\;\;\;\;\tan^{-1}_* \frac{t\_0}{\phi_2 \cdot \cos \phi_1 - \sin \phi_1 \cdot \cos \left(\lambda_1 - \lambda_2\right)}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if phi2 < -0.0027000000000000001 or 3.9999999999999999e-69 < phi2 Initial program 78.9%
Taylor expanded in phi2 around 0
lower-sin.f6455.7
Simplified55.7%
Taylor expanded in lambda2 around 0
lower-cos.f6455.4
Simplified55.4%
Taylor expanded in phi1 around 0
lower-sin.f6454.2
Simplified54.2%
if -0.0027000000000000001 < phi2 < 3.9999999999999999e-69Initial program 87.3%
Taylor expanded in phi2 around 0
lower-sin.f6487.2
Simplified87.2%
Taylor expanded in phi2 around 0
lower-sin.f64N/A
lower--.f6487.3
Simplified87.3%
Taylor expanded in phi2 around 0
lower-*.f64N/A
lower-cos.f6487.3
Simplified87.3%
(FPCore (lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (cos (- lambda1 lambda2))) (t_1 (sin (- lambda1 lambda2))))
(if (<= phi1 -0.001)
(atan2 t_1 (- (sin phi2) (* (sin phi1) t_0)))
(if (<= phi1 2.55e-14)
(atan2 (* t_1 (cos phi2)) (- (sin phi2) (* phi1 t_0)))
(atan2 t_1 (* t_0 (- (sin phi1))))))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = cos((lambda1 - lambda2));
double t_1 = sin((lambda1 - lambda2));
double tmp;
if (phi1 <= -0.001) {
tmp = atan2(t_1, (sin(phi2) - (sin(phi1) * t_0)));
} else if (phi1 <= 2.55e-14) {
tmp = atan2((t_1 * cos(phi2)), (sin(phi2) - (phi1 * t_0)));
} else {
tmp = atan2(t_1, (t_0 * -sin(phi1)));
}
return tmp;
}
real(8) function code(lambda1, lambda2, phi1, phi2)
real(8), intent (in) :: lambda1
real(8), intent (in) :: lambda2
real(8), intent (in) :: phi1
real(8), intent (in) :: phi2
real(8) :: t_0
real(8) :: t_1
real(8) :: tmp
t_0 = cos((lambda1 - lambda2))
t_1 = sin((lambda1 - lambda2))
if (phi1 <= (-0.001d0)) then
tmp = atan2(t_1, (sin(phi2) - (sin(phi1) * t_0)))
else if (phi1 <= 2.55d-14) then
tmp = atan2((t_1 * cos(phi2)), (sin(phi2) - (phi1 * t_0)))
else
tmp = atan2(t_1, (t_0 * -sin(phi1)))
end if
code = tmp
end function
public static double code(double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = Math.cos((lambda1 - lambda2));
double t_1 = Math.sin((lambda1 - lambda2));
double tmp;
if (phi1 <= -0.001) {
tmp = Math.atan2(t_1, (Math.sin(phi2) - (Math.sin(phi1) * t_0)));
} else if (phi1 <= 2.55e-14) {
tmp = Math.atan2((t_1 * Math.cos(phi2)), (Math.sin(phi2) - (phi1 * t_0)));
} else {
tmp = Math.atan2(t_1, (t_0 * -Math.sin(phi1)));
}
return tmp;
}
def code(lambda1, lambda2, phi1, phi2): t_0 = math.cos((lambda1 - lambda2)) t_1 = math.sin((lambda1 - lambda2)) tmp = 0 if phi1 <= -0.001: tmp = math.atan2(t_1, (math.sin(phi2) - (math.sin(phi1) * t_0))) elif phi1 <= 2.55e-14: tmp = math.atan2((t_1 * math.cos(phi2)), (math.sin(phi2) - (phi1 * t_0))) else: tmp = math.atan2(t_1, (t_0 * -math.sin(phi1))) return tmp
function code(lambda1, lambda2, phi1, phi2) t_0 = cos(Float64(lambda1 - lambda2)) t_1 = sin(Float64(lambda1 - lambda2)) tmp = 0.0 if (phi1 <= -0.001) tmp = atan(t_1, Float64(sin(phi2) - Float64(sin(phi1) * t_0))); elseif (phi1 <= 2.55e-14) tmp = atan(Float64(t_1 * cos(phi2)), Float64(sin(phi2) - Float64(phi1 * t_0))); else tmp = atan(t_1, Float64(t_0 * Float64(-sin(phi1)))); end return tmp end
function tmp_2 = code(lambda1, lambda2, phi1, phi2) t_0 = cos((lambda1 - lambda2)); t_1 = sin((lambda1 - lambda2)); tmp = 0.0; if (phi1 <= -0.001) tmp = atan2(t_1, (sin(phi2) - (sin(phi1) * t_0))); elseif (phi1 <= 2.55e-14) tmp = atan2((t_1 * cos(phi2)), (sin(phi2) - (phi1 * t_0))); else tmp = atan2(t_1, (t_0 * -sin(phi1))); end tmp_2 = tmp; end
code[lambda1_, lambda2_, phi1_, phi2_] := Block[{t$95$0 = N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$1 = N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]}, If[LessEqual[phi1, -0.001], N[ArcTan[t$95$1 / N[(N[Sin[phi2], $MachinePrecision] - N[(N[Sin[phi1], $MachinePrecision] * t$95$0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], If[LessEqual[phi1, 2.55e-14], N[ArcTan[N[(t$95$1 * N[Cos[phi2], $MachinePrecision]), $MachinePrecision] / N[(N[Sin[phi2], $MachinePrecision] - N[(phi1 * t$95$0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], N[ArcTan[t$95$1 / N[(t$95$0 * (-N[Sin[phi1], $MachinePrecision])), $MachinePrecision]], $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \cos \left(\lambda_1 - \lambda_2\right)\\
t_1 := \sin \left(\lambda_1 - \lambda_2\right)\\
\mathbf{if}\;\phi_1 \leq -0.001:\\
\;\;\;\;\tan^{-1}_* \frac{t\_1}{\sin \phi_2 - \sin \phi_1 \cdot t\_0}\\
\mathbf{elif}\;\phi_1 \leq 2.55 \cdot 10^{-14}:\\
\;\;\;\;\tan^{-1}_* \frac{t\_1 \cdot \cos \phi_2}{\sin \phi_2 - \phi_1 \cdot t\_0}\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1}_* \frac{t\_1}{t\_0 \cdot \left(-\sin \phi_1\right)}\\
\end{array}
\end{array}
if phi1 < -1e-3Initial program 82.7%
Taylor expanded in phi2 around 0
lower-sin.f6452.6
Simplified52.6%
Taylor expanded in phi2 around 0
lower-sin.f64N/A
lower--.f6450.2
Simplified50.2%
Taylor expanded in phi1 around 0
lower-sin.f6449.2
Simplified49.2%
if -1e-3 < phi1 < 2.5499999999999999e-14Initial program 82.2%
Taylor expanded in phi2 around 0
lower-sin.f6481.6
Simplified81.6%
Taylor expanded in phi1 around 0
*-commutativeN/A
lower-*.f64N/A
lower-cos.f64N/A
lower--.f6481.6
Simplified81.6%
Taylor expanded in phi1 around 0
lower-sin.f6481.6
Simplified81.6%
if 2.5499999999999999e-14 < phi1 Initial program 83.6%
Taylor expanded in phi2 around 0
lower-sin.f6461.1
Simplified61.1%
Taylor expanded in phi2 around 0
lower-sin.f64N/A
lower--.f6456.2
Simplified56.2%
Taylor expanded in phi1 around 0
lower-sin.f6454.8
Simplified54.8%
Taylor expanded in phi2 around 0
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
lower-*.f64N/A
lower-cos.f64N/A
lower--.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lower-sin.f6454.9
Simplified54.9%
Final simplification67.6%
(FPCore (lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (cos (- lambda1 lambda2))))
(if (<= phi2 4.2e+17)
(atan2
(sin (- lambda1 lambda2))
(fma
t_0
(- (sin phi1))
(fma
(* phi2 phi2)
(* phi2 (fma 0.008333333333333333 (* phi2 phi2) -0.16666666666666666))
phi2)))
(atan2 (sin (- lambda2)) (- (sin phi2) (* (sin phi1) t_0))))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = cos((lambda1 - lambda2));
double tmp;
if (phi2 <= 4.2e+17) {
tmp = atan2(sin((lambda1 - lambda2)), fma(t_0, -sin(phi1), fma((phi2 * phi2), (phi2 * fma(0.008333333333333333, (phi2 * phi2), -0.16666666666666666)), phi2)));
} else {
tmp = atan2(sin(-lambda2), (sin(phi2) - (sin(phi1) * t_0)));
}
return tmp;
}
function code(lambda1, lambda2, phi1, phi2) t_0 = cos(Float64(lambda1 - lambda2)) tmp = 0.0 if (phi2 <= 4.2e+17) tmp = atan(sin(Float64(lambda1 - lambda2)), fma(t_0, Float64(-sin(phi1)), fma(Float64(phi2 * phi2), Float64(phi2 * fma(0.008333333333333333, Float64(phi2 * phi2), -0.16666666666666666)), phi2))); else tmp = atan(sin(Float64(-lambda2)), Float64(sin(phi2) - Float64(sin(phi1) * t_0))); end return tmp end
code[lambda1_, lambda2_, phi1_, phi2_] := Block[{t$95$0 = N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]}, If[LessEqual[phi2, 4.2e+17], N[ArcTan[N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision] / N[(t$95$0 * (-N[Sin[phi1], $MachinePrecision]) + N[(N[(phi2 * phi2), $MachinePrecision] * N[(phi2 * N[(0.008333333333333333 * N[(phi2 * phi2), $MachinePrecision] + -0.16666666666666666), $MachinePrecision]), $MachinePrecision] + phi2), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], N[ArcTan[N[Sin[(-lambda2)], $MachinePrecision] / N[(N[Sin[phi2], $MachinePrecision] - N[(N[Sin[phi1], $MachinePrecision] * t$95$0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \cos \left(\lambda_1 - \lambda_2\right)\\
\mathbf{if}\;\phi_2 \leq 4.2 \cdot 10^{+17}:\\
\;\;\;\;\tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right)}{\mathsf{fma}\left(t\_0, -\sin \phi_1, \mathsf{fma}\left(\phi_2 \cdot \phi_2, \phi_2 \cdot \mathsf{fma}\left(0.008333333333333333, \phi_2 \cdot \phi_2, -0.16666666666666666\right), \phi_2\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1}_* \frac{\sin \left(-\lambda_2\right)}{\sin \phi_2 - \sin \phi_1 \cdot t\_0}\\
\end{array}
\end{array}
if phi2 < 4.2e17Initial program 85.8%
Taylor expanded in phi2 around 0
lower-sin.f6477.9
Simplified77.9%
Taylor expanded in phi2 around 0
lower-sin.f64N/A
lower--.f6464.7
Simplified64.7%
Taylor expanded in phi1 around 0
lower-sin.f6464.3
Simplified64.3%
Taylor expanded in phi2 around 0
sub-negN/A
+-commutativeN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
lower-fma.f64N/A
lower-cos.f64N/A
lower--.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lower-sin.f64N/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
*-lft-identityN/A
lower-fma.f64N/A
Simplified63.7%
if 4.2e17 < phi2 Initial program 73.4%
Taylor expanded in phi2 around 0
lower-sin.f6446.1
Simplified46.1%
Taylor expanded in phi2 around 0
lower-sin.f64N/A
lower--.f6416.5
Simplified16.5%
Taylor expanded in phi1 around 0
lower-sin.f6415.3
Simplified15.3%
Taylor expanded in lambda1 around 0
lower-sin.f64N/A
lower-neg.f6416.3
Simplified16.3%
Final simplification51.7%
(FPCore (lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (sin (- lambda1 lambda2))) (t_1 (- (sin phi1))))
(if (<= phi2 4.2e+17)
(atan2
t_0
(fma
(cos (- lambda1 lambda2))
t_1
(fma
(* phi2 phi2)
(* phi2 (fma 0.008333333333333333 (* phi2 phi2) -0.16666666666666666))
phi2)))
(atan2 t_0 (fma (cos lambda2) t_1 (sin phi2))))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = sin((lambda1 - lambda2));
double t_1 = -sin(phi1);
double tmp;
if (phi2 <= 4.2e+17) {
tmp = atan2(t_0, fma(cos((lambda1 - lambda2)), t_1, fma((phi2 * phi2), (phi2 * fma(0.008333333333333333, (phi2 * phi2), -0.16666666666666666)), phi2)));
} else {
tmp = atan2(t_0, fma(cos(lambda2), t_1, sin(phi2)));
}
return tmp;
}
function code(lambda1, lambda2, phi1, phi2) t_0 = sin(Float64(lambda1 - lambda2)) t_1 = Float64(-sin(phi1)) tmp = 0.0 if (phi2 <= 4.2e+17) tmp = atan(t_0, fma(cos(Float64(lambda1 - lambda2)), t_1, fma(Float64(phi2 * phi2), Float64(phi2 * fma(0.008333333333333333, Float64(phi2 * phi2), -0.16666666666666666)), phi2))); else tmp = atan(t_0, fma(cos(lambda2), t_1, sin(phi2))); end return tmp end
code[lambda1_, lambda2_, phi1_, phi2_] := Block[{t$95$0 = N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$1 = (-N[Sin[phi1], $MachinePrecision])}, If[LessEqual[phi2, 4.2e+17], N[ArcTan[t$95$0 / N[(N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision] * t$95$1 + N[(N[(phi2 * phi2), $MachinePrecision] * N[(phi2 * N[(0.008333333333333333 * N[(phi2 * phi2), $MachinePrecision] + -0.16666666666666666), $MachinePrecision]), $MachinePrecision] + phi2), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], N[ArcTan[t$95$0 / N[(N[Cos[lambda2], $MachinePrecision] * t$95$1 + N[Sin[phi2], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sin \left(\lambda_1 - \lambda_2\right)\\
t_1 := -\sin \phi_1\\
\mathbf{if}\;\phi_2 \leq 4.2 \cdot 10^{+17}:\\
\;\;\;\;\tan^{-1}_* \frac{t\_0}{\mathsf{fma}\left(\cos \left(\lambda_1 - \lambda_2\right), t\_1, \mathsf{fma}\left(\phi_2 \cdot \phi_2, \phi_2 \cdot \mathsf{fma}\left(0.008333333333333333, \phi_2 \cdot \phi_2, -0.16666666666666666\right), \phi_2\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1}_* \frac{t\_0}{\mathsf{fma}\left(\cos \lambda_2, t\_1, \sin \phi_2\right)}\\
\end{array}
\end{array}
if phi2 < 4.2e17Initial program 85.8%
Taylor expanded in phi2 around 0
lower-sin.f6477.9
Simplified77.9%
Taylor expanded in phi2 around 0
lower-sin.f64N/A
lower--.f6464.7
Simplified64.7%
Taylor expanded in phi1 around 0
lower-sin.f6464.3
Simplified64.3%
Taylor expanded in phi2 around 0
sub-negN/A
+-commutativeN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
lower-fma.f64N/A
lower-cos.f64N/A
lower--.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lower-sin.f64N/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
*-lft-identityN/A
lower-fma.f64N/A
Simplified63.7%
if 4.2e17 < phi2 Initial program 73.4%
Taylor expanded in phi2 around 0
lower-sin.f6446.1
Simplified46.1%
Taylor expanded in phi2 around 0
lower-sin.f64N/A
lower--.f6416.5
Simplified16.5%
Taylor expanded in phi1 around 0
lower-sin.f6415.3
Simplified15.3%
Taylor expanded in lambda1 around 0
sub-negN/A
+-commutativeN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
lower-fma.f64N/A
cos-negN/A
lower-cos.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lower-sin.f64N/A
lower-sin.f6415.4
Simplified15.4%
Final simplification51.4%
(FPCore (lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (sin (- lambda1 lambda2))))
(if (<= phi2 4.2e+17)
(atan2
t_0
(fma
(cos (- lambda1 lambda2))
(- (sin phi1))
(fma
(* phi2 phi2)
(* phi2 (fma 0.008333333333333333 (* phi2 phi2) -0.16666666666666666))
phi2)))
(atan2 t_0 (- (sin phi2) (* (sin phi1) (cos lambda1)))))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = sin((lambda1 - lambda2));
double tmp;
if (phi2 <= 4.2e+17) {
tmp = atan2(t_0, fma(cos((lambda1 - lambda2)), -sin(phi1), fma((phi2 * phi2), (phi2 * fma(0.008333333333333333, (phi2 * phi2), -0.16666666666666666)), phi2)));
} else {
tmp = atan2(t_0, (sin(phi2) - (sin(phi1) * cos(lambda1))));
}
return tmp;
}
function code(lambda1, lambda2, phi1, phi2) t_0 = sin(Float64(lambda1 - lambda2)) tmp = 0.0 if (phi2 <= 4.2e+17) tmp = atan(t_0, fma(cos(Float64(lambda1 - lambda2)), Float64(-sin(phi1)), fma(Float64(phi2 * phi2), Float64(phi2 * fma(0.008333333333333333, Float64(phi2 * phi2), -0.16666666666666666)), phi2))); else tmp = atan(t_0, Float64(sin(phi2) - Float64(sin(phi1) * cos(lambda1)))); end return tmp end
code[lambda1_, lambda2_, phi1_, phi2_] := Block[{t$95$0 = N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]}, If[LessEqual[phi2, 4.2e+17], N[ArcTan[t$95$0 / N[(N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision] * (-N[Sin[phi1], $MachinePrecision]) + N[(N[(phi2 * phi2), $MachinePrecision] * N[(phi2 * N[(0.008333333333333333 * N[(phi2 * phi2), $MachinePrecision] + -0.16666666666666666), $MachinePrecision]), $MachinePrecision] + phi2), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], N[ArcTan[t$95$0 / N[(N[Sin[phi2], $MachinePrecision] - N[(N[Sin[phi1], $MachinePrecision] * N[Cos[lambda1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sin \left(\lambda_1 - \lambda_2\right)\\
\mathbf{if}\;\phi_2 \leq 4.2 \cdot 10^{+17}:\\
\;\;\;\;\tan^{-1}_* \frac{t\_0}{\mathsf{fma}\left(\cos \left(\lambda_1 - \lambda_2\right), -\sin \phi_1, \mathsf{fma}\left(\phi_2 \cdot \phi_2, \phi_2 \cdot \mathsf{fma}\left(0.008333333333333333, \phi_2 \cdot \phi_2, -0.16666666666666666\right), \phi_2\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1}_* \frac{t\_0}{\sin \phi_2 - \sin \phi_1 \cdot \cos \lambda_1}\\
\end{array}
\end{array}
if phi2 < 4.2e17Initial program 85.8%
Taylor expanded in phi2 around 0
lower-sin.f6477.9
Simplified77.9%
Taylor expanded in phi2 around 0
lower-sin.f64N/A
lower--.f6464.7
Simplified64.7%
Taylor expanded in phi1 around 0
lower-sin.f6464.3
Simplified64.3%
Taylor expanded in phi2 around 0
sub-negN/A
+-commutativeN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
lower-fma.f64N/A
lower-cos.f64N/A
lower--.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lower-sin.f64N/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
*-lft-identityN/A
lower-fma.f64N/A
Simplified63.7%
if 4.2e17 < phi2 Initial program 73.4%
Taylor expanded in phi2 around 0
lower-sin.f6446.1
Simplified46.1%
Taylor expanded in phi2 around 0
lower-sin.f64N/A
lower--.f6416.5
Simplified16.5%
Taylor expanded in phi1 around 0
lower-sin.f6415.3
Simplified15.3%
Taylor expanded in lambda2 around 0
lower-cos.f6415.2
Simplified15.2%
Final simplification51.4%
(FPCore (lambda1 lambda2 phi1 phi2) :precision binary64 (atan2 (sin (- lambda1 lambda2)) (- (sin phi2) (* (sin phi1) (cos (- lambda1 lambda2))))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
return atan2(sin((lambda1 - lambda2)), (sin(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(sin((lambda1 - lambda2)), (sin(phi2) - (sin(phi1) * 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.sin(phi2) - (Math.sin(phi1) * Math.cos((lambda1 - lambda2)))));
}
def code(lambda1, lambda2, phi1, phi2): return math.atan2(math.sin((lambda1 - lambda2)), (math.sin(phi2) - (math.sin(phi1) * math.cos((lambda1 - lambda2)))))
function code(lambda1, lambda2, phi1, phi2) return atan(sin(Float64(lambda1 - lambda2)), Float64(sin(phi2) - Float64(sin(phi1) * cos(Float64(lambda1 - lambda2))))) end
function tmp = code(lambda1, lambda2, phi1, phi2) tmp = atan2(sin((lambda1 - lambda2)), (sin(phi2) - (sin(phi1) * cos((lambda1 - lambda2))))); end
code[lambda1_, lambda2_, phi1_, phi2_] := N[ArcTan[N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision] / N[(N[Sin[phi2], $MachinePrecision] - N[(N[Sin[phi1], $MachinePrecision] * N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
\\
\tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right)}{\sin \phi_2 - \sin \phi_1 \cdot \cos \left(\lambda_1 - \lambda_2\right)}
\end{array}
Initial program 82.6%
Taylor expanded in phi2 around 0
lower-sin.f6469.9
Simplified69.9%
Taylor expanded in phi2 around 0
lower-sin.f64N/A
lower--.f6452.4
Simplified52.4%
Taylor expanded in phi1 around 0
lower-sin.f6451.9
Simplified51.9%
(FPCore (lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (sin (- lambda1 lambda2)))
(t_1 (cos (- lambda1 lambda2)))
(t_2 (- (sin phi1))))
(if (<= phi2 4.2e+17)
(atan2
t_0
(fma
t_1
t_2
(fma
(* phi2 phi2)
(* phi2 (fma 0.008333333333333333 (* phi2 phi2) -0.16666666666666666))
phi2)))
(atan2 t_0 (* t_1 t_2)))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = sin((lambda1 - lambda2));
double t_1 = cos((lambda1 - lambda2));
double t_2 = -sin(phi1);
double tmp;
if (phi2 <= 4.2e+17) {
tmp = atan2(t_0, fma(t_1, t_2, fma((phi2 * phi2), (phi2 * fma(0.008333333333333333, (phi2 * phi2), -0.16666666666666666)), phi2)));
} else {
tmp = atan2(t_0, (t_1 * t_2));
}
return tmp;
}
function code(lambda1, lambda2, phi1, phi2) t_0 = sin(Float64(lambda1 - lambda2)) t_1 = cos(Float64(lambda1 - lambda2)) t_2 = Float64(-sin(phi1)) tmp = 0.0 if (phi2 <= 4.2e+17) tmp = atan(t_0, fma(t_1, t_2, fma(Float64(phi2 * phi2), Float64(phi2 * fma(0.008333333333333333, Float64(phi2 * phi2), -0.16666666666666666)), phi2))); else tmp = atan(t_0, Float64(t_1 * t_2)); end return tmp end
code[lambda1_, lambda2_, phi1_, phi2_] := Block[{t$95$0 = N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$1 = N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$2 = (-N[Sin[phi1], $MachinePrecision])}, If[LessEqual[phi2, 4.2e+17], N[ArcTan[t$95$0 / N[(t$95$1 * t$95$2 + N[(N[(phi2 * phi2), $MachinePrecision] * N[(phi2 * N[(0.008333333333333333 * N[(phi2 * phi2), $MachinePrecision] + -0.16666666666666666), $MachinePrecision]), $MachinePrecision] + phi2), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], N[ArcTan[t$95$0 / N[(t$95$1 * t$95$2), $MachinePrecision]], $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sin \left(\lambda_1 - \lambda_2\right)\\
t_1 := \cos \left(\lambda_1 - \lambda_2\right)\\
t_2 := -\sin \phi_1\\
\mathbf{if}\;\phi_2 \leq 4.2 \cdot 10^{+17}:\\
\;\;\;\;\tan^{-1}_* \frac{t\_0}{\mathsf{fma}\left(t\_1, t\_2, \mathsf{fma}\left(\phi_2 \cdot \phi_2, \phi_2 \cdot \mathsf{fma}\left(0.008333333333333333, \phi_2 \cdot \phi_2, -0.16666666666666666\right), \phi_2\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1}_* \frac{t\_0}{t\_1 \cdot t\_2}\\
\end{array}
\end{array}
if phi2 < 4.2e17Initial program 85.8%
Taylor expanded in phi2 around 0
lower-sin.f6477.9
Simplified77.9%
Taylor expanded in phi2 around 0
lower-sin.f64N/A
lower--.f6464.7
Simplified64.7%
Taylor expanded in phi1 around 0
lower-sin.f6464.3
Simplified64.3%
Taylor expanded in phi2 around 0
sub-negN/A
+-commutativeN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
lower-fma.f64N/A
lower-cos.f64N/A
lower--.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lower-sin.f64N/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
*-lft-identityN/A
lower-fma.f64N/A
Simplified63.7%
if 4.2e17 < phi2 Initial program 73.4%
Taylor expanded in phi2 around 0
lower-sin.f6446.1
Simplified46.1%
Taylor expanded in phi2 around 0
lower-sin.f64N/A
lower--.f6416.5
Simplified16.5%
Taylor expanded in phi1 around 0
lower-sin.f6415.3
Simplified15.3%
Taylor expanded in phi2 around 0
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
lower-*.f64N/A
lower-cos.f64N/A
lower--.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lower-sin.f6413.4
Simplified13.4%
Final simplification50.9%
(FPCore (lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (sin (- lambda1 lambda2)))
(t_1 (cos (- lambda1 lambda2)))
(t_2 (atan2 t_0 (* t_1 (- (sin phi1))))))
(if (<= phi1 -0.0006)
t_2
(if (<= phi1 2.55e-14) (atan2 t_0 (fma t_1 (- phi1) (sin phi2))) t_2))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = sin((lambda1 - lambda2));
double t_1 = cos((lambda1 - lambda2));
double t_2 = atan2(t_0, (t_1 * -sin(phi1)));
double tmp;
if (phi1 <= -0.0006) {
tmp = t_2;
} else if (phi1 <= 2.55e-14) {
tmp = atan2(t_0, fma(t_1, -phi1, sin(phi2)));
} else {
tmp = t_2;
}
return tmp;
}
function code(lambda1, lambda2, phi1, phi2) t_0 = sin(Float64(lambda1 - lambda2)) t_1 = cos(Float64(lambda1 - lambda2)) t_2 = atan(t_0, Float64(t_1 * Float64(-sin(phi1)))) tmp = 0.0 if (phi1 <= -0.0006) tmp = t_2; elseif (phi1 <= 2.55e-14) tmp = atan(t_0, fma(t_1, Float64(-phi1), sin(phi2))); else tmp = t_2; end return tmp end
code[lambda1_, lambda2_, phi1_, phi2_] := Block[{t$95$0 = N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$1 = N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$2 = N[ArcTan[t$95$0 / N[(t$95$1 * (-N[Sin[phi1], $MachinePrecision])), $MachinePrecision]], $MachinePrecision]}, If[LessEqual[phi1, -0.0006], t$95$2, If[LessEqual[phi1, 2.55e-14], N[ArcTan[t$95$0 / N[(t$95$1 * (-phi1) + N[Sin[phi2], $MachinePrecision]), $MachinePrecision]], $MachinePrecision], t$95$2]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sin \left(\lambda_1 - \lambda_2\right)\\
t_1 := \cos \left(\lambda_1 - \lambda_2\right)\\
t_2 := \tan^{-1}_* \frac{t\_0}{t\_1 \cdot \left(-\sin \phi_1\right)}\\
\mathbf{if}\;\phi_1 \leq -0.0006:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;\phi_1 \leq 2.55 \cdot 10^{-14}:\\
\;\;\;\;\tan^{-1}_* \frac{t\_0}{\mathsf{fma}\left(t\_1, -\phi_1, \sin \phi_2\right)}\\
\mathbf{else}:\\
\;\;\;\;t\_2\\
\end{array}
\end{array}
if phi1 < -5.99999999999999947e-4 or 2.5499999999999999e-14 < phi1 Initial program 83.1%
Taylor expanded in phi2 around 0
lower-sin.f6456.8
Simplified56.8%
Taylor expanded in phi2 around 0
lower-sin.f64N/A
lower--.f6453.1
Simplified53.1%
Taylor expanded in phi1 around 0
lower-sin.f6452.0
Simplified52.0%
Taylor expanded in phi2 around 0
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
lower-*.f64N/A
lower-cos.f64N/A
lower--.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lower-sin.f6451.8
Simplified51.8%
if -5.99999999999999947e-4 < phi1 < 2.5499999999999999e-14Initial program 82.2%
Taylor expanded in phi2 around 0
lower-sin.f6481.6
Simplified81.6%
Taylor expanded in phi2 around 0
lower-sin.f64N/A
lower--.f6451.8
Simplified51.8%
Taylor expanded in phi1 around 0
lower-sin.f6451.8
Simplified51.8%
Taylor expanded in phi1 around 0
+-commutativeN/A
associate-*r*N/A
*-commutativeN/A
lower-fma.f64N/A
lower-cos.f64N/A
lower--.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lower-sin.f6451.8
Simplified51.8%
(FPCore (lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (sin (- lambda1 lambda2)))
(t_1 (cos (- lambda1 lambda2)))
(t_2 (- (sin phi1))))
(if (<= phi2 4.2e+17)
(atan2 t_0 (fma t_1 t_2 phi2))
(atan2 t_0 (* t_1 t_2)))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = sin((lambda1 - lambda2));
double t_1 = cos((lambda1 - lambda2));
double t_2 = -sin(phi1);
double tmp;
if (phi2 <= 4.2e+17) {
tmp = atan2(t_0, fma(t_1, t_2, phi2));
} else {
tmp = atan2(t_0, (t_1 * t_2));
}
return tmp;
}
function code(lambda1, lambda2, phi1, phi2) t_0 = sin(Float64(lambda1 - lambda2)) t_1 = cos(Float64(lambda1 - lambda2)) t_2 = Float64(-sin(phi1)) tmp = 0.0 if (phi2 <= 4.2e+17) tmp = atan(t_0, fma(t_1, t_2, phi2)); else tmp = atan(t_0, Float64(t_1 * t_2)); end return tmp end
code[lambda1_, lambda2_, phi1_, phi2_] := Block[{t$95$0 = N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$1 = N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$2 = (-N[Sin[phi1], $MachinePrecision])}, If[LessEqual[phi2, 4.2e+17], N[ArcTan[t$95$0 / N[(t$95$1 * t$95$2 + phi2), $MachinePrecision]], $MachinePrecision], N[ArcTan[t$95$0 / N[(t$95$1 * t$95$2), $MachinePrecision]], $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sin \left(\lambda_1 - \lambda_2\right)\\
t_1 := \cos \left(\lambda_1 - \lambda_2\right)\\
t_2 := -\sin \phi_1\\
\mathbf{if}\;\phi_2 \leq 4.2 \cdot 10^{+17}:\\
\;\;\;\;\tan^{-1}_* \frac{t\_0}{\mathsf{fma}\left(t\_1, t\_2, \phi_2\right)}\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1}_* \frac{t\_0}{t\_1 \cdot t\_2}\\
\end{array}
\end{array}
if phi2 < 4.2e17Initial program 85.8%
Taylor expanded in phi2 around 0
lower-sin.f6477.9
Simplified77.9%
Taylor expanded in phi2 around 0
lower-sin.f64N/A
lower--.f6464.7
Simplified64.7%
Taylor expanded in phi1 around 0
lower-sin.f6464.3
Simplified64.3%
Taylor expanded in phi2 around 0
sub-negN/A
+-commutativeN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
lower-fma.f64N/A
lower-cos.f64N/A
lower--.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lower-sin.f6463.7
Simplified63.7%
if 4.2e17 < phi2 Initial program 73.4%
Taylor expanded in phi2 around 0
lower-sin.f6446.1
Simplified46.1%
Taylor expanded in phi2 around 0
lower-sin.f64N/A
lower--.f6416.5
Simplified16.5%
Taylor expanded in phi1 around 0
lower-sin.f6415.3
Simplified15.3%
Taylor expanded in phi2 around 0
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
lower-*.f64N/A
lower-cos.f64N/A
lower--.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lower-sin.f6413.4
Simplified13.4%
(FPCore (lambda1 lambda2 phi1 phi2) :precision binary64 (atan2 (sin (- lambda1 lambda2)) (* (cos (- lambda1 lambda2)) (- (sin phi1)))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
return atan2(sin((lambda1 - lambda2)), (cos((lambda1 - lambda2)) * -sin(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(sin((lambda1 - lambda2)), (cos((lambda1 - lambda2)) * -sin(phi1)))
end function
public static double code(double lambda1, double lambda2, double phi1, double phi2) {
return Math.atan2(Math.sin((lambda1 - lambda2)), (Math.cos((lambda1 - lambda2)) * -Math.sin(phi1)));
}
def code(lambda1, lambda2, phi1, phi2): return math.atan2(math.sin((lambda1 - lambda2)), (math.cos((lambda1 - lambda2)) * -math.sin(phi1)))
function code(lambda1, lambda2, phi1, phi2) return atan(sin(Float64(lambda1 - lambda2)), Float64(cos(Float64(lambda1 - lambda2)) * Float64(-sin(phi1)))) end
function tmp = code(lambda1, lambda2, phi1, phi2) tmp = atan2(sin((lambda1 - lambda2)), (cos((lambda1 - lambda2)) * -sin(phi1))); end
code[lambda1_, lambda2_, phi1_, phi2_] := N[ArcTan[N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision] / N[(N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision] * (-N[Sin[phi1], $MachinePrecision])), $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
\\
\tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right)}{\cos \left(\lambda_1 - \lambda_2\right) \cdot \left(-\sin \phi_1\right)}
\end{array}
Initial program 82.6%
Taylor expanded in phi2 around 0
lower-sin.f6469.9
Simplified69.9%
Taylor expanded in phi2 around 0
lower-sin.f64N/A
lower--.f6452.4
Simplified52.4%
Taylor expanded in phi1 around 0
lower-sin.f6451.9
Simplified51.9%
Taylor expanded in phi2 around 0
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
lower-*.f64N/A
lower-cos.f64N/A
lower--.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lower-sin.f6448.1
Simplified48.1%
(FPCore (lambda1 lambda2 phi1 phi2) :precision binary64 (atan2 (sin (- lambda1 lambda2)) (sin phi2)))
double code(double lambda1, double lambda2, double phi1, double phi2) {
return atan2(sin((lambda1 - lambda2)), sin(phi2));
}
real(8) function code(lambda1, lambda2, phi1, phi2)
real(8), intent (in) :: lambda1
real(8), intent (in) :: lambda2
real(8), intent (in) :: phi1
real(8), intent (in) :: phi2
code = atan2(sin((lambda1 - lambda2)), sin(phi2))
end function
public static double code(double lambda1, double lambda2, double phi1, double phi2) {
return Math.atan2(Math.sin((lambda1 - lambda2)), Math.sin(phi2));
}
def code(lambda1, lambda2, phi1, phi2): return math.atan2(math.sin((lambda1 - lambda2)), math.sin(phi2))
function code(lambda1, lambda2, phi1, phi2) return atan(sin(Float64(lambda1 - lambda2)), sin(phi2)) end
function tmp = code(lambda1, lambda2, phi1, phi2) tmp = atan2(sin((lambda1 - lambda2)), sin(phi2)); end
code[lambda1_, lambda2_, phi1_, phi2_] := N[ArcTan[N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision] / N[Sin[phi2], $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
\\
\tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right)}{\sin \phi_2}
\end{array}
Initial program 82.6%
Taylor expanded in phi2 around 0
lower-sin.f6469.9
Simplified69.9%
Taylor expanded in phi2 around 0
lower-sin.f64N/A
lower--.f6452.4
Simplified52.4%
Taylor expanded in phi1 around 0
lower-sin.f6451.9
Simplified51.9%
Taylor expanded in phi1 around 0
lower-sin.f6434.0
Simplified34.0%
herbie shell --seed 2024215
(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))))))