
(FPCore (R lambda1 lambda2 phi1 phi2)
:precision binary64
(*
(acos
(+
(* (sin phi1) (sin phi2))
(* (* (cos phi1) (cos phi2)) (cos (- lambda1 lambda2)))))
R))
double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
return acos(((sin(phi1) * sin(phi2)) + ((cos(phi1) * cos(phi2)) * cos((lambda1 - lambda2))))) * R;
}
real(8) function code(r, lambda1, lambda2, phi1, phi2)
real(8), intent (in) :: r
real(8), intent (in) :: lambda1
real(8), intent (in) :: lambda2
real(8), intent (in) :: phi1
real(8), intent (in) :: phi2
code = acos(((sin(phi1) * sin(phi2)) + ((cos(phi1) * cos(phi2)) * cos((lambda1 - lambda2))))) * r
end function
public static double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
return Math.acos(((Math.sin(phi1) * Math.sin(phi2)) + ((Math.cos(phi1) * Math.cos(phi2)) * Math.cos((lambda1 - lambda2))))) * R;
}
def code(R, lambda1, lambda2, phi1, phi2): return math.acos(((math.sin(phi1) * math.sin(phi2)) + ((math.cos(phi1) * math.cos(phi2)) * math.cos((lambda1 - lambda2))))) * R
function code(R, lambda1, lambda2, phi1, phi2) return Float64(acos(Float64(Float64(sin(phi1) * sin(phi2)) + Float64(Float64(cos(phi1) * cos(phi2)) * cos(Float64(lambda1 - lambda2))))) * R) end
function tmp = code(R, lambda1, lambda2, phi1, phi2) tmp = acos(((sin(phi1) * sin(phi2)) + ((cos(phi1) * cos(phi2)) * cos((lambda1 - lambda2))))) * R; end
code[R_, lambda1_, lambda2_, phi1_, phi2_] := N[(N[ArcCos[N[(N[(N[Sin[phi1], $MachinePrecision] * N[Sin[phi2], $MachinePrecision]), $MachinePrecision] + N[(N[(N[Cos[phi1], $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision] * N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * R), $MachinePrecision]
\begin{array}{l}
\\
\cos^{-1} \left(\sin \phi_1 \cdot \sin \phi_2 + \left(\cos \phi_1 \cdot \cos \phi_2\right) \cdot \cos \left(\lambda_1 - \lambda_2\right)\right) \cdot R
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 28 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (R lambda1 lambda2 phi1 phi2)
:precision binary64
(*
(acos
(+
(* (sin phi1) (sin phi2))
(* (* (cos phi1) (cos phi2)) (cos (- lambda1 lambda2)))))
R))
double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
return acos(((sin(phi1) * sin(phi2)) + ((cos(phi1) * cos(phi2)) * cos((lambda1 - lambda2))))) * R;
}
real(8) function code(r, lambda1, lambda2, phi1, phi2)
real(8), intent (in) :: r
real(8), intent (in) :: lambda1
real(8), intent (in) :: lambda2
real(8), intent (in) :: phi1
real(8), intent (in) :: phi2
code = acos(((sin(phi1) * sin(phi2)) + ((cos(phi1) * cos(phi2)) * cos((lambda1 - lambda2))))) * r
end function
public static double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
return Math.acos(((Math.sin(phi1) * Math.sin(phi2)) + ((Math.cos(phi1) * Math.cos(phi2)) * Math.cos((lambda1 - lambda2))))) * R;
}
def code(R, lambda1, lambda2, phi1, phi2): return math.acos(((math.sin(phi1) * math.sin(phi2)) + ((math.cos(phi1) * math.cos(phi2)) * math.cos((lambda1 - lambda2))))) * R
function code(R, lambda1, lambda2, phi1, phi2) return Float64(acos(Float64(Float64(sin(phi1) * sin(phi2)) + Float64(Float64(cos(phi1) * cos(phi2)) * cos(Float64(lambda1 - lambda2))))) * R) end
function tmp = code(R, lambda1, lambda2, phi1, phi2) tmp = acos(((sin(phi1) * sin(phi2)) + ((cos(phi1) * cos(phi2)) * cos((lambda1 - lambda2))))) * R; end
code[R_, lambda1_, lambda2_, phi1_, phi2_] := N[(N[ArcCos[N[(N[(N[Sin[phi1], $MachinePrecision] * N[Sin[phi2], $MachinePrecision]), $MachinePrecision] + N[(N[(N[Cos[phi1], $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision] * N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * R), $MachinePrecision]
\begin{array}{l}
\\
\cos^{-1} \left(\sin \phi_1 \cdot \sin \phi_2 + \left(\cos \phi_1 \cdot \cos \phi_2\right) \cdot \cos \left(\lambda_1 - \lambda_2\right)\right) \cdot R
\end{array}
(FPCore (R lambda1 lambda2 phi1 phi2)
:precision binary64
(*
(acos
(fma
(cos phi2)
(*
(cos phi1)
(fma (sin lambda2) (sin lambda1) (* (cos lambda1) (cos lambda2))))
(* (sin phi1) (sin phi2))))
R))
double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
return acos(fma(cos(phi2), (cos(phi1) * fma(sin(lambda2), sin(lambda1), (cos(lambda1) * cos(lambda2)))), (sin(phi1) * sin(phi2)))) * R;
}
function code(R, lambda1, lambda2, phi1, phi2) return Float64(acos(fma(cos(phi2), Float64(cos(phi1) * fma(sin(lambda2), sin(lambda1), Float64(cos(lambda1) * cos(lambda2)))), Float64(sin(phi1) * sin(phi2)))) * R) end
code[R_, lambda1_, lambda2_, phi1_, phi2_] := N[(N[ArcCos[N[(N[Cos[phi2], $MachinePrecision] * N[(N[Cos[phi1], $MachinePrecision] * N[(N[Sin[lambda2], $MachinePrecision] * N[Sin[lambda1], $MachinePrecision] + N[(N[Cos[lambda1], $MachinePrecision] * N[Cos[lambda2], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[Sin[phi1], $MachinePrecision] * N[Sin[phi2], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * R), $MachinePrecision]
\begin{array}{l}
\\
\cos^{-1} \left(\mathsf{fma}\left(\cos \phi_2, \cos \phi_1 \cdot \mathsf{fma}\left(\sin \lambda_2, \sin \lambda_1, \cos \lambda_1 \cdot \cos \lambda_2\right), \sin \phi_1 \cdot \sin \phi_2\right)\right) \cdot R
\end{array}
Initial program 74.0%
cos-diff95.1%
distribute-lft-in95.1%
Applied egg-rr95.1%
distribute-lft-out95.1%
*-commutative95.1%
fma-def95.1%
Simplified95.1%
Taylor expanded in phi1 around 0 95.1%
fma-def95.2%
fma-udef95.1%
*-commutative95.1%
+-commutative95.1%
fma-def95.1%
*-commutative95.1%
Simplified95.1%
Final simplification95.1%
(FPCore (R lambda1 lambda2 phi1 phi2)
:precision binary64
(*
R
(acos
(+
(* (sin phi1) (sin phi2))
(*
(* (cos phi2) (cos phi1))
(fma (cos lambda1) (cos lambda2) (* (sin lambda2) (sin lambda1))))))))
double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
return R * acos(((sin(phi1) * sin(phi2)) + ((cos(phi2) * cos(phi1)) * fma(cos(lambda1), cos(lambda2), (sin(lambda2) * sin(lambda1))))));
}
function code(R, lambda1, lambda2, phi1, phi2) return Float64(R * acos(Float64(Float64(sin(phi1) * sin(phi2)) + Float64(Float64(cos(phi2) * cos(phi1)) * fma(cos(lambda1), cos(lambda2), Float64(sin(lambda2) * sin(lambda1))))))) end
code[R_, lambda1_, lambda2_, phi1_, phi2_] := N[(R * N[ArcCos[N[(N[(N[Sin[phi1], $MachinePrecision] * N[Sin[phi2], $MachinePrecision]), $MachinePrecision] + N[(N[(N[Cos[phi2], $MachinePrecision] * N[Cos[phi1], $MachinePrecision]), $MachinePrecision] * N[(N[Cos[lambda1], $MachinePrecision] * N[Cos[lambda2], $MachinePrecision] + N[(N[Sin[lambda2], $MachinePrecision] * N[Sin[lambda1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
R \cdot \cos^{-1} \left(\sin \phi_1 \cdot \sin \phi_2 + \left(\cos \phi_2 \cdot \cos \phi_1\right) \cdot \mathsf{fma}\left(\cos \lambda_1, \cos \lambda_2, \sin \lambda_2 \cdot \sin \lambda_1\right)\right)
\end{array}
Initial program 74.0%
cos-diff95.1%
distribute-lft-in95.1%
Applied egg-rr95.1%
distribute-lft-out95.1%
*-commutative95.1%
fma-def95.1%
Simplified95.1%
Final simplification95.1%
(FPCore (R lambda1 lambda2 phi1 phi2)
:precision binary64
(*
R
(acos
(fma
(sin phi1)
(sin phi2)
(*
(cos phi1)
(*
(cos phi2)
(+ (* (cos lambda1) (cos lambda2)) (* (sin lambda2) (sin lambda1)))))))))
double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
return R * acos(fma(sin(phi1), sin(phi2), (cos(phi1) * (cos(phi2) * ((cos(lambda1) * cos(lambda2)) + (sin(lambda2) * sin(lambda1)))))));
}
function code(R, lambda1, lambda2, phi1, phi2) return Float64(R * acos(fma(sin(phi1), sin(phi2), Float64(cos(phi1) * Float64(cos(phi2) * Float64(Float64(cos(lambda1) * cos(lambda2)) + Float64(sin(lambda2) * sin(lambda1)))))))) end
code[R_, lambda1_, lambda2_, phi1_, phi2_] := N[(R * N[ArcCos[N[(N[Sin[phi1], $MachinePrecision] * N[Sin[phi2], $MachinePrecision] + N[(N[Cos[phi1], $MachinePrecision] * N[(N[Cos[phi2], $MachinePrecision] * N[(N[(N[Cos[lambda1], $MachinePrecision] * N[Cos[lambda2], $MachinePrecision]), $MachinePrecision] + N[(N[Sin[lambda2], $MachinePrecision] * N[Sin[lambda1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
R \cdot \cos^{-1} \left(\mathsf{fma}\left(\sin \phi_1, \sin \phi_2, \cos \phi_1 \cdot \left(\cos \phi_2 \cdot \left(\cos \lambda_1 \cdot \cos \lambda_2 + \sin \lambda_2 \cdot \sin \lambda_1\right)\right)\right)\right)
\end{array}
Initial program 74.0%
fma-def74.0%
associate-*l*74.0%
Simplified74.0%
cos-diff95.1%
+-commutative95.1%
Applied egg-rr95.1%
Final simplification95.1%
(FPCore (R lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (cos (- lambda1 lambda2))) (t_1 (* (cos phi2) (cos phi1))))
(if (<= phi1 -5e+17)
(* R (acos (fma (sin phi1) (sin phi2) (* (cos phi1) (* (cos phi2) t_0)))))
(if (<= phi1 0.2)
(*
R
(acos
(+
(*
t_1
(fma (cos lambda1) (cos lambda2) (* (sin lambda2) (sin lambda1))))
(* phi1 (sin phi2)))))
(* R (acos (+ (* (sin phi1) (sin phi2)) (* t_1 t_0))))))))
double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = cos((lambda1 - lambda2));
double t_1 = cos(phi2) * cos(phi1);
double tmp;
if (phi1 <= -5e+17) {
tmp = R * acos(fma(sin(phi1), sin(phi2), (cos(phi1) * (cos(phi2) * t_0))));
} else if (phi1 <= 0.2) {
tmp = R * acos(((t_1 * fma(cos(lambda1), cos(lambda2), (sin(lambda2) * sin(lambda1)))) + (phi1 * sin(phi2))));
} else {
tmp = R * acos(((sin(phi1) * sin(phi2)) + (t_1 * t_0)));
}
return tmp;
}
function code(R, lambda1, lambda2, phi1, phi2) t_0 = cos(Float64(lambda1 - lambda2)) t_1 = Float64(cos(phi2) * cos(phi1)) tmp = 0.0 if (phi1 <= -5e+17) tmp = Float64(R * acos(fma(sin(phi1), sin(phi2), Float64(cos(phi1) * Float64(cos(phi2) * t_0))))); elseif (phi1 <= 0.2) tmp = Float64(R * acos(Float64(Float64(t_1 * fma(cos(lambda1), cos(lambda2), Float64(sin(lambda2) * sin(lambda1)))) + Float64(phi1 * sin(phi2))))); else tmp = Float64(R * acos(Float64(Float64(sin(phi1) * sin(phi2)) + Float64(t_1 * t_0)))); end return tmp end
code[R_, lambda1_, lambda2_, phi1_, phi2_] := Block[{t$95$0 = N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$1 = N[(N[Cos[phi2], $MachinePrecision] * N[Cos[phi1], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[phi1, -5e+17], N[(R * N[ArcCos[N[(N[Sin[phi1], $MachinePrecision] * N[Sin[phi2], $MachinePrecision] + N[(N[Cos[phi1], $MachinePrecision] * N[(N[Cos[phi2], $MachinePrecision] * t$95$0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], If[LessEqual[phi1, 0.2], N[(R * N[ArcCos[N[(N[(t$95$1 * N[(N[Cos[lambda1], $MachinePrecision] * N[Cos[lambda2], $MachinePrecision] + N[(N[Sin[lambda2], $MachinePrecision] * N[Sin[lambda1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(phi1 * N[Sin[phi2], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], N[(R * N[ArcCos[N[(N[(N[Sin[phi1], $MachinePrecision] * N[Sin[phi2], $MachinePrecision]), $MachinePrecision] + N[(t$95$1 * t$95$0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \cos \left(\lambda_1 - \lambda_2\right)\\
t_1 := \cos \phi_2 \cdot \cos \phi_1\\
\mathbf{if}\;\phi_1 \leq -5 \cdot 10^{+17}:\\
\;\;\;\;R \cdot \cos^{-1} \left(\mathsf{fma}\left(\sin \phi_1, \sin \phi_2, \cos \phi_1 \cdot \left(\cos \phi_2 \cdot t_0\right)\right)\right)\\
\mathbf{elif}\;\phi_1 \leq 0.2:\\
\;\;\;\;R \cdot \cos^{-1} \left(t_1 \cdot \mathsf{fma}\left(\cos \lambda_1, \cos \lambda_2, \sin \lambda_2 \cdot \sin \lambda_1\right) + \phi_1 \cdot \sin \phi_2\right)\\
\mathbf{else}:\\
\;\;\;\;R \cdot \cos^{-1} \left(\sin \phi_1 \cdot \sin \phi_2 + t_1 \cdot t_0\right)\\
\end{array}
\end{array}
if phi1 < -5e17Initial program 83.9%
fma-def83.9%
associate-*l*83.9%
Simplified83.9%
if -5e17 < phi1 < 0.20000000000000001Initial program 67.5%
cos-diff91.2%
distribute-lft-in91.2%
Applied egg-rr91.2%
distribute-lft-out91.2%
*-commutative91.2%
fma-def91.3%
Simplified91.3%
Taylor expanded in phi1 around 0 89.6%
if 0.20000000000000001 < phi1 Initial program 78.3%
Final simplification85.5%
(FPCore (R lambda1 lambda2 phi1 phi2)
:precision binary64
(*
R
(acos
(+
(* (sin phi1) (sin phi2))
(*
(+ (* (cos lambda1) (cos lambda2)) (* (sin lambda2) (sin lambda1)))
(* (cos phi2) (cos phi1)))))))
double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
return R * acos(((sin(phi1) * sin(phi2)) + (((cos(lambda1) * cos(lambda2)) + (sin(lambda2) * sin(lambda1))) * (cos(phi2) * cos(phi1)))));
}
real(8) function code(r, lambda1, lambda2, phi1, phi2)
real(8), intent (in) :: r
real(8), intent (in) :: lambda1
real(8), intent (in) :: lambda2
real(8), intent (in) :: phi1
real(8), intent (in) :: phi2
code = r * acos(((sin(phi1) * sin(phi2)) + (((cos(lambda1) * cos(lambda2)) + (sin(lambda2) * sin(lambda1))) * (cos(phi2) * cos(phi1)))))
end function
public static double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
return R * Math.acos(((Math.sin(phi1) * Math.sin(phi2)) + (((Math.cos(lambda1) * Math.cos(lambda2)) + (Math.sin(lambda2) * Math.sin(lambda1))) * (Math.cos(phi2) * Math.cos(phi1)))));
}
def code(R, lambda1, lambda2, phi1, phi2): return R * math.acos(((math.sin(phi1) * math.sin(phi2)) + (((math.cos(lambda1) * math.cos(lambda2)) + (math.sin(lambda2) * math.sin(lambda1))) * (math.cos(phi2) * math.cos(phi1)))))
function code(R, lambda1, lambda2, phi1, phi2) return Float64(R * acos(Float64(Float64(sin(phi1) * sin(phi2)) + Float64(Float64(Float64(cos(lambda1) * cos(lambda2)) + Float64(sin(lambda2) * sin(lambda1))) * Float64(cos(phi2) * cos(phi1)))))) end
function tmp = code(R, lambda1, lambda2, phi1, phi2) tmp = R * acos(((sin(phi1) * sin(phi2)) + (((cos(lambda1) * cos(lambda2)) + (sin(lambda2) * sin(lambda1))) * (cos(phi2) * cos(phi1))))); end
code[R_, lambda1_, lambda2_, phi1_, phi2_] := N[(R * N[ArcCos[N[(N[(N[Sin[phi1], $MachinePrecision] * N[Sin[phi2], $MachinePrecision]), $MachinePrecision] + N[(N[(N[(N[Cos[lambda1], $MachinePrecision] * N[Cos[lambda2], $MachinePrecision]), $MachinePrecision] + N[(N[Sin[lambda2], $MachinePrecision] * N[Sin[lambda1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[Cos[phi2], $MachinePrecision] * N[Cos[phi1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
R \cdot \cos^{-1} \left(\sin \phi_1 \cdot \sin \phi_2 + \left(\cos \lambda_1 \cdot \cos \lambda_2 + \sin \lambda_2 \cdot \sin \lambda_1\right) \cdot \left(\cos \phi_2 \cdot \cos \phi_1\right)\right)
\end{array}
Initial program 74.0%
cos-diff95.1%
+-commutative95.1%
Applied egg-rr95.1%
Final simplification95.1%
(FPCore (R lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (cos (- lambda1 lambda2))))
(if (<= phi1 -0.00265)
(* R (acos (fma (sin phi1) (sin phi2) (* (cos phi1) (* (cos phi2) t_0)))))
(if (<= phi1 0.025)
(*
R
(acos
(+
(* phi1 (sin phi2))
(*
(+ (* -0.5 (* phi1 phi1)) 1.0)
(*
(cos phi2)
(fma
(sin lambda2)
(sin lambda1)
(* (cos lambda1) (cos lambda2))))))))
(*
R
(acos
(+ (* (sin phi1) (sin phi2)) (* (* (cos phi2) (cos phi1)) t_0))))))))
double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = cos((lambda1 - lambda2));
double tmp;
if (phi1 <= -0.00265) {
tmp = R * acos(fma(sin(phi1), sin(phi2), (cos(phi1) * (cos(phi2) * t_0))));
} else if (phi1 <= 0.025) {
tmp = R * acos(((phi1 * sin(phi2)) + (((-0.5 * (phi1 * phi1)) + 1.0) * (cos(phi2) * fma(sin(lambda2), sin(lambda1), (cos(lambda1) * cos(lambda2)))))));
} else {
tmp = R * acos(((sin(phi1) * sin(phi2)) + ((cos(phi2) * cos(phi1)) * t_0)));
}
return tmp;
}
function code(R, lambda1, lambda2, phi1, phi2) t_0 = cos(Float64(lambda1 - lambda2)) tmp = 0.0 if (phi1 <= -0.00265) tmp = Float64(R * acos(fma(sin(phi1), sin(phi2), Float64(cos(phi1) * Float64(cos(phi2) * t_0))))); elseif (phi1 <= 0.025) tmp = Float64(R * acos(Float64(Float64(phi1 * sin(phi2)) + Float64(Float64(Float64(-0.5 * Float64(phi1 * phi1)) + 1.0) * Float64(cos(phi2) * fma(sin(lambda2), sin(lambda1), Float64(cos(lambda1) * cos(lambda2)))))))); else tmp = Float64(R * acos(Float64(Float64(sin(phi1) * sin(phi2)) + Float64(Float64(cos(phi2) * cos(phi1)) * t_0)))); end return tmp end
code[R_, lambda1_, lambda2_, phi1_, phi2_] := Block[{t$95$0 = N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]}, If[LessEqual[phi1, -0.00265], N[(R * N[ArcCos[N[(N[Sin[phi1], $MachinePrecision] * N[Sin[phi2], $MachinePrecision] + N[(N[Cos[phi1], $MachinePrecision] * N[(N[Cos[phi2], $MachinePrecision] * t$95$0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], If[LessEqual[phi1, 0.025], N[(R * N[ArcCos[N[(N[(phi1 * N[Sin[phi2], $MachinePrecision]), $MachinePrecision] + N[(N[(N[(-0.5 * N[(phi1 * phi1), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision] * N[(N[Cos[phi2], $MachinePrecision] * N[(N[Sin[lambda2], $MachinePrecision] * N[Sin[lambda1], $MachinePrecision] + N[(N[Cos[lambda1], $MachinePrecision] * N[Cos[lambda2], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], N[(R * N[ArcCos[N[(N[(N[Sin[phi1], $MachinePrecision] * N[Sin[phi2], $MachinePrecision]), $MachinePrecision] + N[(N[(N[Cos[phi2], $MachinePrecision] * N[Cos[phi1], $MachinePrecision]), $MachinePrecision] * t$95$0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \cos \left(\lambda_1 - \lambda_2\right)\\
\mathbf{if}\;\phi_1 \leq -0.00265:\\
\;\;\;\;R \cdot \cos^{-1} \left(\mathsf{fma}\left(\sin \phi_1, \sin \phi_2, \cos \phi_1 \cdot \left(\cos \phi_2 \cdot t_0\right)\right)\right)\\
\mathbf{elif}\;\phi_1 \leq 0.025:\\
\;\;\;\;R \cdot \cos^{-1} \left(\phi_1 \cdot \sin \phi_2 + \left(-0.5 \cdot \left(\phi_1 \cdot \phi_1\right) + 1\right) \cdot \left(\cos \phi_2 \cdot \mathsf{fma}\left(\sin \lambda_2, \sin \lambda_1, \cos \lambda_1 \cdot \cos \lambda_2\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;R \cdot \cos^{-1} \left(\sin \phi_1 \cdot \sin \phi_2 + \left(\cos \phi_2 \cdot \cos \phi_1\right) \cdot t_0\right)\\
\end{array}
\end{array}
if phi1 < -0.00265000000000000001Initial program 84.5%
fma-def84.5%
associate-*l*84.5%
Simplified84.5%
if -0.00265000000000000001 < phi1 < 0.025000000000000001Initial program 66.9%
cos-diff91.1%
distribute-lft-in91.1%
Applied egg-rr91.1%
distribute-lft-out91.1%
*-commutative91.1%
fma-def91.2%
Simplified91.2%
Taylor expanded in phi1 around 0 91.2%
fma-def91.2%
fma-udef91.1%
*-commutative91.1%
+-commutative91.1%
fma-def91.2%
*-commutative91.2%
Simplified91.2%
Taylor expanded in phi1 around 0 90.0%
associate-+r+90.0%
+-commutative90.0%
associate-*r*90.0%
distribute-lft1-in90.0%
unpow290.0%
+-commutative90.0%
Simplified90.0%
if 0.025000000000000001 < phi1 Initial program 78.3%
Final simplification85.7%
(FPCore (R lambda1 lambda2 phi1 phi2)
:precision binary64
(if (or (<= phi2 -1.22e-6) (not (<= phi2 5e-7)))
(*
R
(acos
(fma
(sin phi1)
(sin phi2)
(* (cos phi1) (* (cos phi2) (cos (- lambda1 lambda2)))))))
(*
R
(acos
(+
(* phi2 (sin phi1))
(*
(cos phi1)
(fma (sin lambda1) (sin lambda2) (* (cos lambda1) (cos lambda2)))))))))
double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
double tmp;
if ((phi2 <= -1.22e-6) || !(phi2 <= 5e-7)) {
tmp = R * acos(fma(sin(phi1), sin(phi2), (cos(phi1) * (cos(phi2) * cos((lambda1 - lambda2))))));
} else {
tmp = R * acos(((phi2 * sin(phi1)) + (cos(phi1) * fma(sin(lambda1), sin(lambda2), (cos(lambda1) * cos(lambda2))))));
}
return tmp;
}
function code(R, lambda1, lambda2, phi1, phi2) tmp = 0.0 if ((phi2 <= -1.22e-6) || !(phi2 <= 5e-7)) tmp = Float64(R * acos(fma(sin(phi1), sin(phi2), Float64(cos(phi1) * Float64(cos(phi2) * cos(Float64(lambda1 - lambda2))))))); else tmp = Float64(R * acos(Float64(Float64(phi2 * sin(phi1)) + Float64(cos(phi1) * fma(sin(lambda1), sin(lambda2), Float64(cos(lambda1) * cos(lambda2))))))); end return tmp end
code[R_, lambda1_, lambda2_, phi1_, phi2_] := If[Or[LessEqual[phi2, -1.22e-6], N[Not[LessEqual[phi2, 5e-7]], $MachinePrecision]], N[(R * N[ArcCos[N[(N[Sin[phi1], $MachinePrecision] * N[Sin[phi2], $MachinePrecision] + N[(N[Cos[phi1], $MachinePrecision] * N[(N[Cos[phi2], $MachinePrecision] * N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], N[(R * N[ArcCos[N[(N[(phi2 * N[Sin[phi1], $MachinePrecision]), $MachinePrecision] + N[(N[Cos[phi1], $MachinePrecision] * N[(N[Sin[lambda1], $MachinePrecision] * N[Sin[lambda2], $MachinePrecision] + N[(N[Cos[lambda1], $MachinePrecision] * N[Cos[lambda2], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\phi_2 \leq -1.22 \cdot 10^{-6} \lor \neg \left(\phi_2 \leq 5 \cdot 10^{-7}\right):\\
\;\;\;\;R \cdot \cos^{-1} \left(\mathsf{fma}\left(\sin \phi_1, \sin \phi_2, \cos \phi_1 \cdot \left(\cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;R \cdot \cos^{-1} \left(\phi_2 \cdot \sin \phi_1 + \cos \phi_1 \cdot \mathsf{fma}\left(\sin \lambda_1, \sin \lambda_2, \cos \lambda_1 \cdot \cos \lambda_2\right)\right)\\
\end{array}
\end{array}
if phi2 < -1.21999999999999997e-6 or 4.99999999999999977e-7 < phi2 Initial program 79.1%
fma-def79.1%
associate-*l*79.1%
Simplified79.1%
if -1.21999999999999997e-6 < phi2 < 4.99999999999999977e-7Initial program 69.0%
Taylor expanded in phi2 around 0 69.0%
sub-neg69.0%
+-commutative69.0%
neg-mul-169.0%
neg-mul-169.0%
remove-double-neg69.0%
mul-1-neg69.0%
distribute-neg-in69.0%
+-commutative69.0%
cos-neg69.0%
+-commutative69.0%
mul-1-neg69.0%
unsub-neg69.0%
Simplified69.0%
cos-diff91.0%
*-commutative91.0%
*-commutative91.0%
+-commutative91.0%
fma-def91.0%
Applied egg-rr91.0%
Taylor expanded in phi2 around 0 91.0%
Final simplification85.1%
(FPCore (R lambda1 lambda2 phi1 phi2)
:precision binary64
(if (<= phi2 -1.2e-6)
(*
R
(-
(* PI 0.5)
(asin
(fma
(cos (- lambda2 lambda1))
(* (cos phi2) (cos phi1))
(* (sin phi1) (sin phi2))))))
(if (<= phi2 5.5e-7)
(*
R
(acos
(+
(* phi2 (sin phi1))
(*
(cos phi1)
(fma (sin lambda1) (sin lambda2) (* (cos lambda1) (cos lambda2)))))))
(*
R
(acos
(fma
(sin phi1)
(sin phi2)
(* (cos phi1) (* (cos phi2) (cos (- lambda1 lambda2))))))))))
double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
double tmp;
if (phi2 <= -1.2e-6) {
tmp = R * ((((double) M_PI) * 0.5) - asin(fma(cos((lambda2 - lambda1)), (cos(phi2) * cos(phi1)), (sin(phi1) * sin(phi2)))));
} else if (phi2 <= 5.5e-7) {
tmp = R * acos(((phi2 * sin(phi1)) + (cos(phi1) * fma(sin(lambda1), sin(lambda2), (cos(lambda1) * cos(lambda2))))));
} else {
tmp = R * acos(fma(sin(phi1), sin(phi2), (cos(phi1) * (cos(phi2) * cos((lambda1 - lambda2))))));
}
return tmp;
}
function code(R, lambda1, lambda2, phi1, phi2) tmp = 0.0 if (phi2 <= -1.2e-6) tmp = Float64(R * Float64(Float64(pi * 0.5) - asin(fma(cos(Float64(lambda2 - lambda1)), Float64(cos(phi2) * cos(phi1)), Float64(sin(phi1) * sin(phi2)))))); elseif (phi2 <= 5.5e-7) tmp = Float64(R * acos(Float64(Float64(phi2 * sin(phi1)) + Float64(cos(phi1) * fma(sin(lambda1), sin(lambda2), Float64(cos(lambda1) * cos(lambda2))))))); else tmp = Float64(R * acos(fma(sin(phi1), sin(phi2), Float64(cos(phi1) * Float64(cos(phi2) * cos(Float64(lambda1 - lambda2))))))); end return tmp end
code[R_, lambda1_, lambda2_, phi1_, phi2_] := If[LessEqual[phi2, -1.2e-6], N[(R * N[(N[(Pi * 0.5), $MachinePrecision] - N[ArcSin[N[(N[Cos[N[(lambda2 - lambda1), $MachinePrecision]], $MachinePrecision] * N[(N[Cos[phi2], $MachinePrecision] * N[Cos[phi1], $MachinePrecision]), $MachinePrecision] + N[(N[Sin[phi1], $MachinePrecision] * N[Sin[phi2], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[phi2, 5.5e-7], N[(R * N[ArcCos[N[(N[(phi2 * N[Sin[phi1], $MachinePrecision]), $MachinePrecision] + N[(N[Cos[phi1], $MachinePrecision] * N[(N[Sin[lambda1], $MachinePrecision] * N[Sin[lambda2], $MachinePrecision] + N[(N[Cos[lambda1], $MachinePrecision] * N[Cos[lambda2], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], N[(R * N[ArcCos[N[(N[Sin[phi1], $MachinePrecision] * N[Sin[phi2], $MachinePrecision] + N[(N[Cos[phi1], $MachinePrecision] * N[(N[Cos[phi2], $MachinePrecision] * N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\phi_2 \leq -1.2 \cdot 10^{-6}:\\
\;\;\;\;R \cdot \left(\pi \cdot 0.5 - \sin^{-1} \left(\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \cos \phi_2 \cdot \cos \phi_1, \sin \phi_1 \cdot \sin \phi_2\right)\right)\right)\\
\mathbf{elif}\;\phi_2 \leq 5.5 \cdot 10^{-7}:\\
\;\;\;\;R \cdot \cos^{-1} \left(\phi_2 \cdot \sin \phi_1 + \cos \phi_1 \cdot \mathsf{fma}\left(\sin \lambda_1, \sin \lambda_2, \cos \lambda_1 \cdot \cos \lambda_2\right)\right)\\
\mathbf{else}:\\
\;\;\;\;R \cdot \cos^{-1} \left(\mathsf{fma}\left(\sin \phi_1, \sin \phi_2, \cos \phi_1 \cdot \left(\cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)\right)\right)\right)\\
\end{array}
\end{array}
if phi2 < -1.1999999999999999e-6Initial program 80.6%
acos-asin80.6%
sub-neg80.6%
div-inv80.6%
metadata-eval80.6%
+-commutative80.6%
*-commutative80.6%
fma-def80.7%
Applied egg-rr80.7%
sub-neg80.7%
*-commutative80.7%
sub-neg80.7%
+-commutative80.7%
neg-mul-180.7%
neg-mul-180.7%
remove-double-neg80.7%
mul-1-neg80.7%
distribute-neg-in80.7%
+-commutative80.7%
cos-neg80.7%
+-commutative80.7%
mul-1-neg80.7%
unsub-neg80.7%
Simplified80.7%
if -1.1999999999999999e-6 < phi2 < 5.5000000000000003e-7Initial program 69.0%
Taylor expanded in phi2 around 0 69.0%
sub-neg69.0%
+-commutative69.0%
neg-mul-169.0%
neg-mul-169.0%
remove-double-neg69.0%
mul-1-neg69.0%
distribute-neg-in69.0%
+-commutative69.0%
cos-neg69.0%
+-commutative69.0%
mul-1-neg69.0%
unsub-neg69.0%
Simplified69.0%
cos-diff91.0%
*-commutative91.0%
*-commutative91.0%
+-commutative91.0%
fma-def91.0%
Applied egg-rr91.0%
Taylor expanded in phi2 around 0 91.0%
if 5.5000000000000003e-7 < phi2 Initial program 77.3%
fma-def77.3%
associate-*l*77.3%
Simplified77.3%
Final simplification85.1%
(FPCore (R lambda1 lambda2 phi1 phi2)
:precision binary64
(if (or (<= phi2 -7.2e-7) (not (<= phi2 3.6e-6)))
(*
R
(acos
(fma
(sin phi1)
(sin phi2)
(* (cos phi1) (* (cos phi2) (cos (- lambda1 lambda2)))))))
(*
R
(acos
(+
(* phi2 (sin phi1))
(*
(cos phi1)
(+
(* (cos lambda1) (cos lambda2))
(* (sin lambda2) (sin lambda1)))))))))
double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
double tmp;
if ((phi2 <= -7.2e-7) || !(phi2 <= 3.6e-6)) {
tmp = R * acos(fma(sin(phi1), sin(phi2), (cos(phi1) * (cos(phi2) * cos((lambda1 - lambda2))))));
} else {
tmp = R * acos(((phi2 * sin(phi1)) + (cos(phi1) * ((cos(lambda1) * cos(lambda2)) + (sin(lambda2) * sin(lambda1))))));
}
return tmp;
}
function code(R, lambda1, lambda2, phi1, phi2) tmp = 0.0 if ((phi2 <= -7.2e-7) || !(phi2 <= 3.6e-6)) tmp = Float64(R * acos(fma(sin(phi1), sin(phi2), Float64(cos(phi1) * Float64(cos(phi2) * cos(Float64(lambda1 - lambda2))))))); else tmp = Float64(R * acos(Float64(Float64(phi2 * sin(phi1)) + Float64(cos(phi1) * Float64(Float64(cos(lambda1) * cos(lambda2)) + Float64(sin(lambda2) * sin(lambda1))))))); end return tmp end
code[R_, lambda1_, lambda2_, phi1_, phi2_] := If[Or[LessEqual[phi2, -7.2e-7], N[Not[LessEqual[phi2, 3.6e-6]], $MachinePrecision]], N[(R * N[ArcCos[N[(N[Sin[phi1], $MachinePrecision] * N[Sin[phi2], $MachinePrecision] + N[(N[Cos[phi1], $MachinePrecision] * N[(N[Cos[phi2], $MachinePrecision] * N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], N[(R * N[ArcCos[N[(N[(phi2 * N[Sin[phi1], $MachinePrecision]), $MachinePrecision] + N[(N[Cos[phi1], $MachinePrecision] * N[(N[(N[Cos[lambda1], $MachinePrecision] * N[Cos[lambda2], $MachinePrecision]), $MachinePrecision] + N[(N[Sin[lambda2], $MachinePrecision] * N[Sin[lambda1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\phi_2 \leq -7.2 \cdot 10^{-7} \lor \neg \left(\phi_2 \leq 3.6 \cdot 10^{-6}\right):\\
\;\;\;\;R \cdot \cos^{-1} \left(\mathsf{fma}\left(\sin \phi_1, \sin \phi_2, \cos \phi_1 \cdot \left(\cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;R \cdot \cos^{-1} \left(\phi_2 \cdot \sin \phi_1 + \cos \phi_1 \cdot \left(\cos \lambda_1 \cdot \cos \lambda_2 + \sin \lambda_2 \cdot \sin \lambda_1\right)\right)\\
\end{array}
\end{array}
if phi2 < -7.19999999999999989e-7 or 3.59999999999999984e-6 < phi2 Initial program 79.1%
fma-def79.1%
associate-*l*79.1%
Simplified79.1%
if -7.19999999999999989e-7 < phi2 < 3.59999999999999984e-6Initial program 69.0%
cos-diff91.1%
distribute-lft-in91.1%
Applied egg-rr91.1%
distribute-lft-out91.1%
*-commutative91.1%
fma-def91.1%
Simplified91.1%
Taylor expanded in phi1 around 0 91.1%
fma-def91.1%
fma-udef91.1%
*-commutative91.1%
+-commutative91.1%
fma-def91.1%
*-commutative91.1%
Simplified91.1%
Taylor expanded in phi2 around 0 91.0%
Final simplification85.1%
(FPCore (R lambda1 lambda2 phi1 phi2)
:precision binary64
(if (or (<= phi2 -8.2e-8) (not (<= phi2 6.2e-8)))
(*
R
(acos
(fma
(sin phi1)
(sin phi2)
(* (cos phi1) (* (cos phi2) (cos (- lambda1 lambda2)))))))
(*
R
(acos
(*
(cos phi1)
(fma (sin lambda2) (sin lambda1) (* (cos lambda1) (cos lambda2))))))))
double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
double tmp;
if ((phi2 <= -8.2e-8) || !(phi2 <= 6.2e-8)) {
tmp = R * acos(fma(sin(phi1), sin(phi2), (cos(phi1) * (cos(phi2) * cos((lambda1 - lambda2))))));
} else {
tmp = R * acos((cos(phi1) * fma(sin(lambda2), sin(lambda1), (cos(lambda1) * cos(lambda2)))));
}
return tmp;
}
function code(R, lambda1, lambda2, phi1, phi2) tmp = 0.0 if ((phi2 <= -8.2e-8) || !(phi2 <= 6.2e-8)) tmp = Float64(R * acos(fma(sin(phi1), sin(phi2), Float64(cos(phi1) * Float64(cos(phi2) * cos(Float64(lambda1 - lambda2))))))); else tmp = Float64(R * acos(Float64(cos(phi1) * fma(sin(lambda2), sin(lambda1), Float64(cos(lambda1) * cos(lambda2)))))); end return tmp end
code[R_, lambda1_, lambda2_, phi1_, phi2_] := If[Or[LessEqual[phi2, -8.2e-8], N[Not[LessEqual[phi2, 6.2e-8]], $MachinePrecision]], N[(R * N[ArcCos[N[(N[Sin[phi1], $MachinePrecision] * N[Sin[phi2], $MachinePrecision] + N[(N[Cos[phi1], $MachinePrecision] * N[(N[Cos[phi2], $MachinePrecision] * N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], N[(R * N[ArcCos[N[(N[Cos[phi1], $MachinePrecision] * N[(N[Sin[lambda2], $MachinePrecision] * N[Sin[lambda1], $MachinePrecision] + N[(N[Cos[lambda1], $MachinePrecision] * N[Cos[lambda2], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\phi_2 \leq -8.2 \cdot 10^{-8} \lor \neg \left(\phi_2 \leq 6.2 \cdot 10^{-8}\right):\\
\;\;\;\;R \cdot \cos^{-1} \left(\mathsf{fma}\left(\sin \phi_1, \sin \phi_2, \cos \phi_1 \cdot \left(\cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;R \cdot \cos^{-1} \left(\cos \phi_1 \cdot \mathsf{fma}\left(\sin \lambda_2, \sin \lambda_1, \cos \lambda_1 \cdot \cos \lambda_2\right)\right)\\
\end{array}
\end{array}
if phi2 < -8.20000000000000063e-8 or 6.2e-8 < phi2 Initial program 79.1%
fma-def79.1%
associate-*l*79.1%
Simplified79.1%
if -8.20000000000000063e-8 < phi2 < 6.2e-8Initial program 69.0%
cos-diff91.1%
distribute-lft-in91.1%
Applied egg-rr91.1%
distribute-lft-out91.1%
*-commutative91.1%
fma-def91.1%
Simplified91.1%
Taylor expanded in phi1 around 0 91.1%
fma-def91.1%
fma-udef91.1%
*-commutative91.1%
+-commutative91.1%
fma-def91.1%
*-commutative91.1%
Simplified91.1%
Taylor expanded in phi2 around 0 90.2%
+-commutative90.2%
fma-def90.2%
Simplified90.2%
Final simplification84.7%
(FPCore (R lambda1 lambda2 phi1 phi2)
:precision binary64
(if (or (<= phi2 -2.65e-8) (not (<= phi2 9e-8)))
(*
R
(acos
(+
(* (sin phi1) (sin phi2))
(* (* (cos phi2) (cos phi1)) (cos (- lambda1 lambda2))))))
(*
R
(acos
(*
(cos phi1)
(fma (sin lambda2) (sin lambda1) (* (cos lambda1) (cos lambda2))))))))
double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
double tmp;
if ((phi2 <= -2.65e-8) || !(phi2 <= 9e-8)) {
tmp = R * acos(((sin(phi1) * sin(phi2)) + ((cos(phi2) * cos(phi1)) * cos((lambda1 - lambda2)))));
} else {
tmp = R * acos((cos(phi1) * fma(sin(lambda2), sin(lambda1), (cos(lambda1) * cos(lambda2)))));
}
return tmp;
}
function code(R, lambda1, lambda2, phi1, phi2) tmp = 0.0 if ((phi2 <= -2.65e-8) || !(phi2 <= 9e-8)) tmp = Float64(R * acos(Float64(Float64(sin(phi1) * sin(phi2)) + Float64(Float64(cos(phi2) * cos(phi1)) * cos(Float64(lambda1 - lambda2)))))); else tmp = Float64(R * acos(Float64(cos(phi1) * fma(sin(lambda2), sin(lambda1), Float64(cos(lambda1) * cos(lambda2)))))); end return tmp end
code[R_, lambda1_, lambda2_, phi1_, phi2_] := If[Or[LessEqual[phi2, -2.65e-8], N[Not[LessEqual[phi2, 9e-8]], $MachinePrecision]], N[(R * N[ArcCos[N[(N[(N[Sin[phi1], $MachinePrecision] * N[Sin[phi2], $MachinePrecision]), $MachinePrecision] + N[(N[(N[Cos[phi2], $MachinePrecision] * N[Cos[phi1], $MachinePrecision]), $MachinePrecision] * N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], N[(R * N[ArcCos[N[(N[Cos[phi1], $MachinePrecision] * N[(N[Sin[lambda2], $MachinePrecision] * N[Sin[lambda1], $MachinePrecision] + N[(N[Cos[lambda1], $MachinePrecision] * N[Cos[lambda2], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\phi_2 \leq -2.65 \cdot 10^{-8} \lor \neg \left(\phi_2 \leq 9 \cdot 10^{-8}\right):\\
\;\;\;\;R \cdot \cos^{-1} \left(\sin \phi_1 \cdot \sin \phi_2 + \left(\cos \phi_2 \cdot \cos \phi_1\right) \cdot \cos \left(\lambda_1 - \lambda_2\right)\right)\\
\mathbf{else}:\\
\;\;\;\;R \cdot \cos^{-1} \left(\cos \phi_1 \cdot \mathsf{fma}\left(\sin \lambda_2, \sin \lambda_1, \cos \lambda_1 \cdot \cos \lambda_2\right)\right)\\
\end{array}
\end{array}
if phi2 < -2.6499999999999999e-8 or 8.99999999999999986e-8 < phi2 Initial program 78.7%
if -2.6499999999999999e-8 < phi2 < 8.99999999999999986e-8Initial program 69.2%
cos-diff91.0%
distribute-lft-in91.0%
Applied egg-rr91.0%
distribute-lft-out91.0%
*-commutative91.0%
fma-def91.0%
Simplified91.0%
Taylor expanded in phi1 around 0 91.0%
fma-def91.0%
fma-udef91.0%
*-commutative91.0%
+-commutative91.0%
fma-def91.0%
*-commutative91.0%
Simplified91.0%
Taylor expanded in phi2 around 0 90.4%
+-commutative90.4%
fma-def90.4%
Simplified90.4%
Final simplification84.5%
(FPCore (R lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (* (sin phi1) (sin phi2)))
(t_1
(+ (* (cos lambda1) (cos lambda2)) (* (sin lambda2) (sin lambda1)))))
(if (<= lambda1 -6.3e+88)
(* R (acos (* (cos phi1) t_1)))
(if (<= lambda1 -5e-5)
(* R (acos (+ t_0 (* (cos lambda1) (* (cos phi2) (cos phi1))))))
(if (<= lambda1 6.6e-18)
(* R (acos (+ t_0 (* (cos phi1) (* (cos phi2) (cos lambda2))))))
(* R (acos (* (cos phi2) t_1))))))))
double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = sin(phi1) * sin(phi2);
double t_1 = (cos(lambda1) * cos(lambda2)) + (sin(lambda2) * sin(lambda1));
double tmp;
if (lambda1 <= -6.3e+88) {
tmp = R * acos((cos(phi1) * t_1));
} else if (lambda1 <= -5e-5) {
tmp = R * acos((t_0 + (cos(lambda1) * (cos(phi2) * cos(phi1)))));
} else if (lambda1 <= 6.6e-18) {
tmp = R * acos((t_0 + (cos(phi1) * (cos(phi2) * cos(lambda2)))));
} else {
tmp = R * acos((cos(phi2) * t_1));
}
return tmp;
}
real(8) function code(r, lambda1, lambda2, phi1, phi2)
real(8), intent (in) :: r
real(8), intent (in) :: lambda1
real(8), intent (in) :: lambda2
real(8), intent (in) :: phi1
real(8), intent (in) :: phi2
real(8) :: t_0
real(8) :: t_1
real(8) :: tmp
t_0 = sin(phi1) * sin(phi2)
t_1 = (cos(lambda1) * cos(lambda2)) + (sin(lambda2) * sin(lambda1))
if (lambda1 <= (-6.3d+88)) then
tmp = r * acos((cos(phi1) * t_1))
else if (lambda1 <= (-5d-5)) then
tmp = r * acos((t_0 + (cos(lambda1) * (cos(phi2) * cos(phi1)))))
else if (lambda1 <= 6.6d-18) then
tmp = r * acos((t_0 + (cos(phi1) * (cos(phi2) * cos(lambda2)))))
else
tmp = r * acos((cos(phi2) * t_1))
end if
code = tmp
end function
public static double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = Math.sin(phi1) * Math.sin(phi2);
double t_1 = (Math.cos(lambda1) * Math.cos(lambda2)) + (Math.sin(lambda2) * Math.sin(lambda1));
double tmp;
if (lambda1 <= -6.3e+88) {
tmp = R * Math.acos((Math.cos(phi1) * t_1));
} else if (lambda1 <= -5e-5) {
tmp = R * Math.acos((t_0 + (Math.cos(lambda1) * (Math.cos(phi2) * Math.cos(phi1)))));
} else if (lambda1 <= 6.6e-18) {
tmp = R * Math.acos((t_0 + (Math.cos(phi1) * (Math.cos(phi2) * Math.cos(lambda2)))));
} else {
tmp = R * Math.acos((Math.cos(phi2) * t_1));
}
return tmp;
}
def code(R, lambda1, lambda2, phi1, phi2): t_0 = math.sin(phi1) * math.sin(phi2) t_1 = (math.cos(lambda1) * math.cos(lambda2)) + (math.sin(lambda2) * math.sin(lambda1)) tmp = 0 if lambda1 <= -6.3e+88: tmp = R * math.acos((math.cos(phi1) * t_1)) elif lambda1 <= -5e-5: tmp = R * math.acos((t_0 + (math.cos(lambda1) * (math.cos(phi2) * math.cos(phi1))))) elif lambda1 <= 6.6e-18: tmp = R * math.acos((t_0 + (math.cos(phi1) * (math.cos(phi2) * math.cos(lambda2))))) else: tmp = R * math.acos((math.cos(phi2) * t_1)) return tmp
function code(R, lambda1, lambda2, phi1, phi2) t_0 = Float64(sin(phi1) * sin(phi2)) t_1 = Float64(Float64(cos(lambda1) * cos(lambda2)) + Float64(sin(lambda2) * sin(lambda1))) tmp = 0.0 if (lambda1 <= -6.3e+88) tmp = Float64(R * acos(Float64(cos(phi1) * t_1))); elseif (lambda1 <= -5e-5) tmp = Float64(R * acos(Float64(t_0 + Float64(cos(lambda1) * Float64(cos(phi2) * cos(phi1)))))); elseif (lambda1 <= 6.6e-18) tmp = Float64(R * acos(Float64(t_0 + Float64(cos(phi1) * Float64(cos(phi2) * cos(lambda2)))))); else tmp = Float64(R * acos(Float64(cos(phi2) * t_1))); end return tmp end
function tmp_2 = code(R, lambda1, lambda2, phi1, phi2) t_0 = sin(phi1) * sin(phi2); t_1 = (cos(lambda1) * cos(lambda2)) + (sin(lambda2) * sin(lambda1)); tmp = 0.0; if (lambda1 <= -6.3e+88) tmp = R * acos((cos(phi1) * t_1)); elseif (lambda1 <= -5e-5) tmp = R * acos((t_0 + (cos(lambda1) * (cos(phi2) * cos(phi1))))); elseif (lambda1 <= 6.6e-18) tmp = R * acos((t_0 + (cos(phi1) * (cos(phi2) * cos(lambda2))))); else tmp = R * acos((cos(phi2) * t_1)); end tmp_2 = tmp; end
code[R_, lambda1_, lambda2_, phi1_, phi2_] := Block[{t$95$0 = N[(N[Sin[phi1], $MachinePrecision] * N[Sin[phi2], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[(N[Cos[lambda1], $MachinePrecision] * N[Cos[lambda2], $MachinePrecision]), $MachinePrecision] + N[(N[Sin[lambda2], $MachinePrecision] * N[Sin[lambda1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[lambda1, -6.3e+88], N[(R * N[ArcCos[N[(N[Cos[phi1], $MachinePrecision] * t$95$1), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], If[LessEqual[lambda1, -5e-5], N[(R * N[ArcCos[N[(t$95$0 + N[(N[Cos[lambda1], $MachinePrecision] * N[(N[Cos[phi2], $MachinePrecision] * N[Cos[phi1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], If[LessEqual[lambda1, 6.6e-18], N[(R * N[ArcCos[N[(t$95$0 + N[(N[Cos[phi1], $MachinePrecision] * N[(N[Cos[phi2], $MachinePrecision] * N[Cos[lambda2], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], N[(R * N[ArcCos[N[(N[Cos[phi2], $MachinePrecision] * t$95$1), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sin \phi_1 \cdot \sin \phi_2\\
t_1 := \cos \lambda_1 \cdot \cos \lambda_2 + \sin \lambda_2 \cdot \sin \lambda_1\\
\mathbf{if}\;\lambda_1 \leq -6.3 \cdot 10^{+88}:\\
\;\;\;\;R \cdot \cos^{-1} \left(\cos \phi_1 \cdot t_1\right)\\
\mathbf{elif}\;\lambda_1 \leq -5 \cdot 10^{-5}:\\
\;\;\;\;R \cdot \cos^{-1} \left(t_0 + \cos \lambda_1 \cdot \left(\cos \phi_2 \cdot \cos \phi_1\right)\right)\\
\mathbf{elif}\;\lambda_1 \leq 6.6 \cdot 10^{-18}:\\
\;\;\;\;R \cdot \cos^{-1} \left(t_0 + \cos \phi_1 \cdot \left(\cos \phi_2 \cdot \cos \lambda_2\right)\right)\\
\mathbf{else}:\\
\;\;\;\;R \cdot \cos^{-1} \left(\cos \phi_2 \cdot t_1\right)\\
\end{array}
\end{array}
if lambda1 < -6.30000000000000001e88Initial program 49.3%
cos-diff99.0%
distribute-lft-in99.0%
Applied egg-rr99.0%
distribute-lft-out99.0%
*-commutative99.0%
fma-def99.0%
Simplified99.0%
Taylor expanded in phi1 around 0 99.1%
fma-def99.1%
fma-udef99.1%
*-commutative99.1%
+-commutative99.1%
fma-def99.1%
*-commutative99.1%
Simplified99.1%
Taylor expanded in phi2 around 0 61.5%
if -6.30000000000000001e88 < lambda1 < -5.00000000000000024e-5Initial program 65.7%
Taylor expanded in lambda2 around 0 65.7%
if -5.00000000000000024e-5 < lambda1 < 6.6000000000000003e-18Initial program 90.6%
Taylor expanded in lambda1 around 0 90.6%
associate-*r*90.6%
*-commutative90.6%
associate-*l*90.6%
cos-neg90.6%
Simplified90.6%
if 6.6000000000000003e-18 < lambda1 Initial program 59.3%
cos-diff99.2%
distribute-lft-in99.2%
Applied egg-rr99.2%
distribute-lft-out99.2%
*-commutative99.2%
fma-def99.3%
Simplified99.3%
Taylor expanded in phi1 around 0 99.3%
fma-def99.3%
fma-udef99.3%
*-commutative99.3%
+-commutative99.3%
fma-def99.3%
*-commutative99.3%
Simplified99.3%
Taylor expanded in phi1 around 0 62.0%
Final simplification76.7%
(FPCore (R lambda1 lambda2 phi1 phi2)
:precision binary64
(if (or (<= phi2 -3.5e-10) (not (<= phi2 7e-8)))
(*
R
(acos
(+
(* (sin phi1) (sin phi2))
(* (* (cos phi2) (cos phi1)) (cos (- lambda1 lambda2))))))
(*
R
(acos
(*
(cos phi1)
(+ (* (cos lambda1) (cos lambda2)) (* (sin lambda2) (sin lambda1))))))))
double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
double tmp;
if ((phi2 <= -3.5e-10) || !(phi2 <= 7e-8)) {
tmp = R * acos(((sin(phi1) * sin(phi2)) + ((cos(phi2) * cos(phi1)) * cos((lambda1 - lambda2)))));
} else {
tmp = R * acos((cos(phi1) * ((cos(lambda1) * cos(lambda2)) + (sin(lambda2) * sin(lambda1)))));
}
return tmp;
}
real(8) function code(r, lambda1, lambda2, phi1, phi2)
real(8), intent (in) :: r
real(8), intent (in) :: lambda1
real(8), intent (in) :: lambda2
real(8), intent (in) :: phi1
real(8), intent (in) :: phi2
real(8) :: tmp
if ((phi2 <= (-3.5d-10)) .or. (.not. (phi2 <= 7d-8))) then
tmp = r * acos(((sin(phi1) * sin(phi2)) + ((cos(phi2) * cos(phi1)) * cos((lambda1 - lambda2)))))
else
tmp = r * acos((cos(phi1) * ((cos(lambda1) * cos(lambda2)) + (sin(lambda2) * sin(lambda1)))))
end if
code = tmp
end function
public static double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
double tmp;
if ((phi2 <= -3.5e-10) || !(phi2 <= 7e-8)) {
tmp = R * Math.acos(((Math.sin(phi1) * Math.sin(phi2)) + ((Math.cos(phi2) * Math.cos(phi1)) * Math.cos((lambda1 - lambda2)))));
} else {
tmp = R * Math.acos((Math.cos(phi1) * ((Math.cos(lambda1) * Math.cos(lambda2)) + (Math.sin(lambda2) * Math.sin(lambda1)))));
}
return tmp;
}
def code(R, lambda1, lambda2, phi1, phi2): tmp = 0 if (phi2 <= -3.5e-10) or not (phi2 <= 7e-8): tmp = R * math.acos(((math.sin(phi1) * math.sin(phi2)) + ((math.cos(phi2) * math.cos(phi1)) * math.cos((lambda1 - lambda2))))) else: tmp = R * math.acos((math.cos(phi1) * ((math.cos(lambda1) * math.cos(lambda2)) + (math.sin(lambda2) * math.sin(lambda1))))) return tmp
function code(R, lambda1, lambda2, phi1, phi2) tmp = 0.0 if ((phi2 <= -3.5e-10) || !(phi2 <= 7e-8)) tmp = Float64(R * acos(Float64(Float64(sin(phi1) * sin(phi2)) + Float64(Float64(cos(phi2) * cos(phi1)) * cos(Float64(lambda1 - lambda2)))))); else tmp = Float64(R * acos(Float64(cos(phi1) * Float64(Float64(cos(lambda1) * cos(lambda2)) + Float64(sin(lambda2) * sin(lambda1)))))); end return tmp end
function tmp_2 = code(R, lambda1, lambda2, phi1, phi2) tmp = 0.0; if ((phi2 <= -3.5e-10) || ~((phi2 <= 7e-8))) tmp = R * acos(((sin(phi1) * sin(phi2)) + ((cos(phi2) * cos(phi1)) * cos((lambda1 - lambda2))))); else tmp = R * acos((cos(phi1) * ((cos(lambda1) * cos(lambda2)) + (sin(lambda2) * sin(lambda1))))); end tmp_2 = tmp; end
code[R_, lambda1_, lambda2_, phi1_, phi2_] := If[Or[LessEqual[phi2, -3.5e-10], N[Not[LessEqual[phi2, 7e-8]], $MachinePrecision]], N[(R * N[ArcCos[N[(N[(N[Sin[phi1], $MachinePrecision] * N[Sin[phi2], $MachinePrecision]), $MachinePrecision] + N[(N[(N[Cos[phi2], $MachinePrecision] * N[Cos[phi1], $MachinePrecision]), $MachinePrecision] * N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], N[(R * N[ArcCos[N[(N[Cos[phi1], $MachinePrecision] * N[(N[(N[Cos[lambda1], $MachinePrecision] * N[Cos[lambda2], $MachinePrecision]), $MachinePrecision] + N[(N[Sin[lambda2], $MachinePrecision] * N[Sin[lambda1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\phi_2 \leq -3.5 \cdot 10^{-10} \lor \neg \left(\phi_2 \leq 7 \cdot 10^{-8}\right):\\
\;\;\;\;R \cdot \cos^{-1} \left(\sin \phi_1 \cdot \sin \phi_2 + \left(\cos \phi_2 \cdot \cos \phi_1\right) \cdot \cos \left(\lambda_1 - \lambda_2\right)\right)\\
\mathbf{else}:\\
\;\;\;\;R \cdot \cos^{-1} \left(\cos \phi_1 \cdot \left(\cos \lambda_1 \cdot \cos \lambda_2 + \sin \lambda_2 \cdot \sin \lambda_1\right)\right)\\
\end{array}
\end{array}
if phi2 < -3.4999999999999998e-10 or 7.00000000000000048e-8 < phi2 Initial program 78.7%
if -3.4999999999999998e-10 < phi2 < 7.00000000000000048e-8Initial program 69.2%
cos-diff91.0%
distribute-lft-in91.0%
Applied egg-rr91.0%
distribute-lft-out91.0%
*-commutative91.0%
fma-def91.0%
Simplified91.0%
Taylor expanded in phi1 around 0 91.0%
fma-def91.0%
fma-udef91.0%
*-commutative91.0%
+-commutative91.0%
fma-def91.0%
*-commutative91.0%
Simplified91.0%
Taylor expanded in phi2 around 0 90.4%
Final simplification84.5%
(FPCore (R lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0
(+ (* (cos lambda1) (cos lambda2)) (* (sin lambda2) (sin lambda1)))))
(if (<= phi2 -7.3e-6)
(*
R
(acos
(+
(* (sin phi1) (sin phi2))
(* (cos phi1) (* (cos phi2) (cos lambda2))))))
(if (<= phi2 0.022)
(* R (acos (* (cos phi1) t_0)))
(* R (acos (* (cos phi2) t_0)))))))
double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = (cos(lambda1) * cos(lambda2)) + (sin(lambda2) * sin(lambda1));
double tmp;
if (phi2 <= -7.3e-6) {
tmp = R * acos(((sin(phi1) * sin(phi2)) + (cos(phi1) * (cos(phi2) * cos(lambda2)))));
} else if (phi2 <= 0.022) {
tmp = R * acos((cos(phi1) * t_0));
} else {
tmp = R * acos((cos(phi2) * t_0));
}
return tmp;
}
real(8) function code(r, lambda1, lambda2, phi1, phi2)
real(8), intent (in) :: r
real(8), intent (in) :: lambda1
real(8), intent (in) :: lambda2
real(8), intent (in) :: phi1
real(8), intent (in) :: phi2
real(8) :: t_0
real(8) :: tmp
t_0 = (cos(lambda1) * cos(lambda2)) + (sin(lambda2) * sin(lambda1))
if (phi2 <= (-7.3d-6)) then
tmp = r * acos(((sin(phi1) * sin(phi2)) + (cos(phi1) * (cos(phi2) * cos(lambda2)))))
else if (phi2 <= 0.022d0) then
tmp = r * acos((cos(phi1) * t_0))
else
tmp = r * acos((cos(phi2) * t_0))
end if
code = tmp
end function
public static double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = (Math.cos(lambda1) * Math.cos(lambda2)) + (Math.sin(lambda2) * Math.sin(lambda1));
double tmp;
if (phi2 <= -7.3e-6) {
tmp = R * Math.acos(((Math.sin(phi1) * Math.sin(phi2)) + (Math.cos(phi1) * (Math.cos(phi2) * Math.cos(lambda2)))));
} else if (phi2 <= 0.022) {
tmp = R * Math.acos((Math.cos(phi1) * t_0));
} else {
tmp = R * Math.acos((Math.cos(phi2) * t_0));
}
return tmp;
}
def code(R, lambda1, lambda2, phi1, phi2): t_0 = (math.cos(lambda1) * math.cos(lambda2)) + (math.sin(lambda2) * math.sin(lambda1)) tmp = 0 if phi2 <= -7.3e-6: tmp = R * math.acos(((math.sin(phi1) * math.sin(phi2)) + (math.cos(phi1) * (math.cos(phi2) * math.cos(lambda2))))) elif phi2 <= 0.022: tmp = R * math.acos((math.cos(phi1) * t_0)) else: tmp = R * math.acos((math.cos(phi2) * t_0)) return tmp
function code(R, lambda1, lambda2, phi1, phi2) t_0 = Float64(Float64(cos(lambda1) * cos(lambda2)) + Float64(sin(lambda2) * sin(lambda1))) tmp = 0.0 if (phi2 <= -7.3e-6) tmp = Float64(R * acos(Float64(Float64(sin(phi1) * sin(phi2)) + Float64(cos(phi1) * Float64(cos(phi2) * cos(lambda2)))))); elseif (phi2 <= 0.022) tmp = Float64(R * acos(Float64(cos(phi1) * t_0))); else tmp = Float64(R * acos(Float64(cos(phi2) * t_0))); end return tmp end
function tmp_2 = code(R, lambda1, lambda2, phi1, phi2) t_0 = (cos(lambda1) * cos(lambda2)) + (sin(lambda2) * sin(lambda1)); tmp = 0.0; if (phi2 <= -7.3e-6) tmp = R * acos(((sin(phi1) * sin(phi2)) + (cos(phi1) * (cos(phi2) * cos(lambda2))))); elseif (phi2 <= 0.022) tmp = R * acos((cos(phi1) * t_0)); else tmp = R * acos((cos(phi2) * t_0)); end tmp_2 = tmp; end
code[R_, lambda1_, lambda2_, phi1_, phi2_] := Block[{t$95$0 = N[(N[(N[Cos[lambda1], $MachinePrecision] * N[Cos[lambda2], $MachinePrecision]), $MachinePrecision] + N[(N[Sin[lambda2], $MachinePrecision] * N[Sin[lambda1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[phi2, -7.3e-6], N[(R * N[ArcCos[N[(N[(N[Sin[phi1], $MachinePrecision] * N[Sin[phi2], $MachinePrecision]), $MachinePrecision] + N[(N[Cos[phi1], $MachinePrecision] * N[(N[Cos[phi2], $MachinePrecision] * N[Cos[lambda2], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], If[LessEqual[phi2, 0.022], N[(R * N[ArcCos[N[(N[Cos[phi1], $MachinePrecision] * t$95$0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], N[(R * N[ArcCos[N[(N[Cos[phi2], $MachinePrecision] * t$95$0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \cos \lambda_1 \cdot \cos \lambda_2 + \sin \lambda_2 \cdot \sin \lambda_1\\
\mathbf{if}\;\phi_2 \leq -7.3 \cdot 10^{-6}:\\
\;\;\;\;R \cdot \cos^{-1} \left(\sin \phi_1 \cdot \sin \phi_2 + \cos \phi_1 \cdot \left(\cos \phi_2 \cdot \cos \lambda_2\right)\right)\\
\mathbf{elif}\;\phi_2 \leq 0.022:\\
\;\;\;\;R \cdot \cos^{-1} \left(\cos \phi_1 \cdot t_0\right)\\
\mathbf{else}:\\
\;\;\;\;R \cdot \cos^{-1} \left(\cos \phi_2 \cdot t_0\right)\\
\end{array}
\end{array}
if phi2 < -7.30000000000000041e-6Initial program 80.6%
Taylor expanded in lambda1 around 0 60.2%
associate-*r*60.2%
*-commutative60.2%
associate-*l*60.2%
cos-neg60.2%
Simplified60.2%
if -7.30000000000000041e-6 < phi2 < 0.021999999999999999Initial program 69.3%
cos-diff91.2%
distribute-lft-in91.2%
Applied egg-rr91.2%
distribute-lft-out91.2%
*-commutative91.2%
fma-def91.2%
Simplified91.2%
Taylor expanded in phi1 around 0 91.2%
fma-def91.2%
fma-udef91.2%
*-commutative91.2%
+-commutative91.2%
fma-def91.2%
*-commutative91.2%
Simplified91.2%
Taylor expanded in phi2 around 0 89.8%
if 0.021999999999999999 < phi2 Initial program 76.9%
cos-diff99.2%
distribute-lft-in99.2%
Applied egg-rr99.2%
distribute-lft-out99.2%
*-commutative99.2%
fma-def99.2%
Simplified99.2%
Taylor expanded in phi1 around 0 99.2%
fma-def99.2%
fma-udef99.2%
*-commutative99.2%
+-commutative99.2%
fma-def99.2%
*-commutative99.2%
Simplified99.2%
Taylor expanded in phi1 around 0 62.2%
Final simplification75.7%
(FPCore (R lambda1 lambda2 phi1 phi2)
:precision binary64
(if (<= phi2 0.022)
(*
R
(acos
(*
(cos phi1)
(+ (* (cos lambda1) (cos lambda2)) (* (sin lambda2) (sin lambda1))))))
(*
R
(acos
(+ (* (sin phi1) (sin phi2)) (* (cos phi2) (cos (- lambda2 lambda1))))))))
double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
double tmp;
if (phi2 <= 0.022) {
tmp = R * acos((cos(phi1) * ((cos(lambda1) * cos(lambda2)) + (sin(lambda2) * sin(lambda1)))));
} else {
tmp = R * acos(((sin(phi1) * sin(phi2)) + (cos(phi2) * cos((lambda2 - lambda1)))));
}
return tmp;
}
real(8) function code(r, lambda1, lambda2, phi1, phi2)
real(8), intent (in) :: r
real(8), intent (in) :: lambda1
real(8), intent (in) :: lambda2
real(8), intent (in) :: phi1
real(8), intent (in) :: phi2
real(8) :: tmp
if (phi2 <= 0.022d0) then
tmp = r * acos((cos(phi1) * ((cos(lambda1) * cos(lambda2)) + (sin(lambda2) * sin(lambda1)))))
else
tmp = r * acos(((sin(phi1) * sin(phi2)) + (cos(phi2) * cos((lambda2 - lambda1)))))
end if
code = tmp
end function
public static double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
double tmp;
if (phi2 <= 0.022) {
tmp = R * Math.acos((Math.cos(phi1) * ((Math.cos(lambda1) * Math.cos(lambda2)) + (Math.sin(lambda2) * Math.sin(lambda1)))));
} else {
tmp = R * Math.acos(((Math.sin(phi1) * Math.sin(phi2)) + (Math.cos(phi2) * Math.cos((lambda2 - lambda1)))));
}
return tmp;
}
def code(R, lambda1, lambda2, phi1, phi2): tmp = 0 if phi2 <= 0.022: tmp = R * math.acos((math.cos(phi1) * ((math.cos(lambda1) * math.cos(lambda2)) + (math.sin(lambda2) * math.sin(lambda1))))) else: tmp = R * math.acos(((math.sin(phi1) * math.sin(phi2)) + (math.cos(phi2) * math.cos((lambda2 - lambda1))))) return tmp
function code(R, lambda1, lambda2, phi1, phi2) tmp = 0.0 if (phi2 <= 0.022) tmp = Float64(R * acos(Float64(cos(phi1) * Float64(Float64(cos(lambda1) * cos(lambda2)) + Float64(sin(lambda2) * sin(lambda1)))))); else tmp = Float64(R * acos(Float64(Float64(sin(phi1) * sin(phi2)) + Float64(cos(phi2) * cos(Float64(lambda2 - lambda1)))))); end return tmp end
function tmp_2 = code(R, lambda1, lambda2, phi1, phi2) tmp = 0.0; if (phi2 <= 0.022) tmp = R * acos((cos(phi1) * ((cos(lambda1) * cos(lambda2)) + (sin(lambda2) * sin(lambda1))))); else tmp = R * acos(((sin(phi1) * sin(phi2)) + (cos(phi2) * cos((lambda2 - lambda1))))); end tmp_2 = tmp; end
code[R_, lambda1_, lambda2_, phi1_, phi2_] := If[LessEqual[phi2, 0.022], N[(R * N[ArcCos[N[(N[Cos[phi1], $MachinePrecision] * N[(N[(N[Cos[lambda1], $MachinePrecision] * N[Cos[lambda2], $MachinePrecision]), $MachinePrecision] + N[(N[Sin[lambda2], $MachinePrecision] * N[Sin[lambda1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], N[(R * N[ArcCos[N[(N[(N[Sin[phi1], $MachinePrecision] * N[Sin[phi2], $MachinePrecision]), $MachinePrecision] + N[(N[Cos[phi2], $MachinePrecision] * N[Cos[N[(lambda2 - lambda1), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\phi_2 \leq 0.022:\\
\;\;\;\;R \cdot \cos^{-1} \left(\cos \phi_1 \cdot \left(\cos \lambda_1 \cdot \cos \lambda_2 + \sin \lambda_2 \cdot \sin \lambda_1\right)\right)\\
\mathbf{else}:\\
\;\;\;\;R \cdot \cos^{-1} \left(\sin \phi_1 \cdot \sin \phi_2 + \cos \phi_2 \cdot \cos \left(\lambda_2 - \lambda_1\right)\right)\\
\end{array}
\end{array}
if phi2 < 0.021999999999999999Initial program 73.1%
cos-diff93.9%
distribute-lft-in93.8%
Applied egg-rr93.8%
distribute-lft-out93.9%
*-commutative93.9%
fma-def93.9%
Simplified93.9%
Taylor expanded in phi1 around 0 93.9%
fma-def93.9%
fma-udef93.9%
*-commutative93.9%
+-commutative93.9%
fma-def93.9%
*-commutative93.9%
Simplified93.9%
Taylor expanded in phi2 around 0 65.7%
if 0.021999999999999999 < phi2 Initial program 76.9%
Taylor expanded in phi1 around 0 50.0%
sub-neg50.0%
+-commutative50.0%
neg-mul-150.0%
neg-mul-150.0%
remove-double-neg50.0%
mul-1-neg50.0%
distribute-neg-in50.0%
+-commutative50.0%
cos-neg50.0%
+-commutative50.0%
mul-1-neg50.0%
unsub-neg50.0%
Simplified50.0%
Final simplification62.1%
(FPCore (R lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0
(+ (* (cos lambda1) (cos lambda2)) (* (sin lambda2) (sin lambda1)))))
(if (<= phi2 0.022)
(* R (acos (* (cos phi1) t_0)))
(* R (acos (* (cos phi2) t_0))))))
double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = (cos(lambda1) * cos(lambda2)) + (sin(lambda2) * sin(lambda1));
double tmp;
if (phi2 <= 0.022) {
tmp = R * acos((cos(phi1) * t_0));
} else {
tmp = R * acos((cos(phi2) * t_0));
}
return tmp;
}
real(8) function code(r, lambda1, lambda2, phi1, phi2)
real(8), intent (in) :: r
real(8), intent (in) :: lambda1
real(8), intent (in) :: lambda2
real(8), intent (in) :: phi1
real(8), intent (in) :: phi2
real(8) :: t_0
real(8) :: tmp
t_0 = (cos(lambda1) * cos(lambda2)) + (sin(lambda2) * sin(lambda1))
if (phi2 <= 0.022d0) then
tmp = r * acos((cos(phi1) * t_0))
else
tmp = r * acos((cos(phi2) * t_0))
end if
code = tmp
end function
public static double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = (Math.cos(lambda1) * Math.cos(lambda2)) + (Math.sin(lambda2) * Math.sin(lambda1));
double tmp;
if (phi2 <= 0.022) {
tmp = R * Math.acos((Math.cos(phi1) * t_0));
} else {
tmp = R * Math.acos((Math.cos(phi2) * t_0));
}
return tmp;
}
def code(R, lambda1, lambda2, phi1, phi2): t_0 = (math.cos(lambda1) * math.cos(lambda2)) + (math.sin(lambda2) * math.sin(lambda1)) tmp = 0 if phi2 <= 0.022: tmp = R * math.acos((math.cos(phi1) * t_0)) else: tmp = R * math.acos((math.cos(phi2) * t_0)) return tmp
function code(R, lambda1, lambda2, phi1, phi2) t_0 = Float64(Float64(cos(lambda1) * cos(lambda2)) + Float64(sin(lambda2) * sin(lambda1))) tmp = 0.0 if (phi2 <= 0.022) tmp = Float64(R * acos(Float64(cos(phi1) * t_0))); else tmp = Float64(R * acos(Float64(cos(phi2) * t_0))); end return tmp end
function tmp_2 = code(R, lambda1, lambda2, phi1, phi2) t_0 = (cos(lambda1) * cos(lambda2)) + (sin(lambda2) * sin(lambda1)); tmp = 0.0; if (phi2 <= 0.022) tmp = R * acos((cos(phi1) * t_0)); else tmp = R * acos((cos(phi2) * t_0)); end tmp_2 = tmp; end
code[R_, lambda1_, lambda2_, phi1_, phi2_] := Block[{t$95$0 = N[(N[(N[Cos[lambda1], $MachinePrecision] * N[Cos[lambda2], $MachinePrecision]), $MachinePrecision] + N[(N[Sin[lambda2], $MachinePrecision] * N[Sin[lambda1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[phi2, 0.022], N[(R * N[ArcCos[N[(N[Cos[phi1], $MachinePrecision] * t$95$0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], N[(R * N[ArcCos[N[(N[Cos[phi2], $MachinePrecision] * t$95$0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \cos \lambda_1 \cdot \cos \lambda_2 + \sin \lambda_2 \cdot \sin \lambda_1\\
\mathbf{if}\;\phi_2 \leq 0.022:\\
\;\;\;\;R \cdot \cos^{-1} \left(\cos \phi_1 \cdot t_0\right)\\
\mathbf{else}:\\
\;\;\;\;R \cdot \cos^{-1} \left(\cos \phi_2 \cdot t_0\right)\\
\end{array}
\end{array}
if phi2 < 0.021999999999999999Initial program 73.1%
cos-diff93.9%
distribute-lft-in93.8%
Applied egg-rr93.8%
distribute-lft-out93.9%
*-commutative93.9%
fma-def93.9%
Simplified93.9%
Taylor expanded in phi1 around 0 93.9%
fma-def93.9%
fma-udef93.9%
*-commutative93.9%
+-commutative93.9%
fma-def93.9%
*-commutative93.9%
Simplified93.9%
Taylor expanded in phi2 around 0 65.7%
if 0.021999999999999999 < phi2 Initial program 76.9%
cos-diff99.2%
distribute-lft-in99.2%
Applied egg-rr99.2%
distribute-lft-out99.2%
*-commutative99.2%
fma-def99.2%
Simplified99.2%
Taylor expanded in phi1 around 0 99.2%
fma-def99.2%
fma-udef99.2%
*-commutative99.2%
+-commutative99.2%
fma-def99.2%
*-commutative99.2%
Simplified99.2%
Taylor expanded in phi1 around 0 62.2%
Final simplification64.9%
(FPCore (R lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (cos (- lambda2 lambda1))))
(if (<= phi1 -0.0002)
(*
R
(acos
(+
(/ (- (cos (- phi1 phi2)) (cos (+ phi2 phi1))) 2.0)
(* (cos phi1) t_0))))
(* R (acos (+ (* (sin phi1) (sin phi2)) (* (cos phi2) t_0)))))))
double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = cos((lambda2 - lambda1));
double tmp;
if (phi1 <= -0.0002) {
tmp = R * acos((((cos((phi1 - phi2)) - cos((phi2 + phi1))) / 2.0) + (cos(phi1) * t_0)));
} else {
tmp = R * acos(((sin(phi1) * sin(phi2)) + (cos(phi2) * t_0)));
}
return tmp;
}
real(8) function code(r, lambda1, lambda2, phi1, phi2)
real(8), intent (in) :: r
real(8), intent (in) :: lambda1
real(8), intent (in) :: lambda2
real(8), intent (in) :: phi1
real(8), intent (in) :: phi2
real(8) :: t_0
real(8) :: tmp
t_0 = cos((lambda2 - lambda1))
if (phi1 <= (-0.0002d0)) then
tmp = r * acos((((cos((phi1 - phi2)) - cos((phi2 + phi1))) / 2.0d0) + (cos(phi1) * t_0)))
else
tmp = r * acos(((sin(phi1) * sin(phi2)) + (cos(phi2) * t_0)))
end if
code = tmp
end function
public static double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = Math.cos((lambda2 - lambda1));
double tmp;
if (phi1 <= -0.0002) {
tmp = R * Math.acos((((Math.cos((phi1 - phi2)) - Math.cos((phi2 + phi1))) / 2.0) + (Math.cos(phi1) * t_0)));
} else {
tmp = R * Math.acos(((Math.sin(phi1) * Math.sin(phi2)) + (Math.cos(phi2) * t_0)));
}
return tmp;
}
def code(R, lambda1, lambda2, phi1, phi2): t_0 = math.cos((lambda2 - lambda1)) tmp = 0 if phi1 <= -0.0002: tmp = R * math.acos((((math.cos((phi1 - phi2)) - math.cos((phi2 + phi1))) / 2.0) + (math.cos(phi1) * t_0))) else: tmp = R * math.acos(((math.sin(phi1) * math.sin(phi2)) + (math.cos(phi2) * t_0))) return tmp
function code(R, lambda1, lambda2, phi1, phi2) t_0 = cos(Float64(lambda2 - lambda1)) tmp = 0.0 if (phi1 <= -0.0002) tmp = Float64(R * acos(Float64(Float64(Float64(cos(Float64(phi1 - phi2)) - cos(Float64(phi2 + phi1))) / 2.0) + Float64(cos(phi1) * t_0)))); else tmp = Float64(R * acos(Float64(Float64(sin(phi1) * sin(phi2)) + Float64(cos(phi2) * t_0)))); end return tmp end
function tmp_2 = code(R, lambda1, lambda2, phi1, phi2) t_0 = cos((lambda2 - lambda1)); tmp = 0.0; if (phi1 <= -0.0002) tmp = R * acos((((cos((phi1 - phi2)) - cos((phi2 + phi1))) / 2.0) + (cos(phi1) * t_0))); else tmp = R * acos(((sin(phi1) * sin(phi2)) + (cos(phi2) * t_0))); end tmp_2 = tmp; end
code[R_, lambda1_, lambda2_, phi1_, phi2_] := Block[{t$95$0 = N[Cos[N[(lambda2 - lambda1), $MachinePrecision]], $MachinePrecision]}, If[LessEqual[phi1, -0.0002], N[(R * N[ArcCos[N[(N[(N[(N[Cos[N[(phi1 - phi2), $MachinePrecision]], $MachinePrecision] - N[Cos[N[(phi2 + phi1), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / 2.0), $MachinePrecision] + N[(N[Cos[phi1], $MachinePrecision] * t$95$0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], N[(R * N[ArcCos[N[(N[(N[Sin[phi1], $MachinePrecision] * N[Sin[phi2], $MachinePrecision]), $MachinePrecision] + N[(N[Cos[phi2], $MachinePrecision] * t$95$0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \cos \left(\lambda_2 - \lambda_1\right)\\
\mathbf{if}\;\phi_1 \leq -0.0002:\\
\;\;\;\;R \cdot \cos^{-1} \left(\frac{\cos \left(\phi_1 - \phi_2\right) - \cos \left(\phi_2 + \phi_1\right)}{2} + \cos \phi_1 \cdot t_0\right)\\
\mathbf{else}:\\
\;\;\;\;R \cdot \cos^{-1} \left(\sin \phi_1 \cdot \sin \phi_2 + \cos \phi_2 \cdot t_0\right)\\
\end{array}
\end{array}
if phi1 < -2.0000000000000001e-4Initial program 83.5%
Taylor expanded in phi2 around 0 51.1%
sub-neg51.1%
+-commutative51.1%
neg-mul-151.1%
neg-mul-151.1%
remove-double-neg51.1%
mul-1-neg51.1%
distribute-neg-in51.1%
+-commutative51.1%
cos-neg51.1%
+-commutative51.1%
mul-1-neg51.1%
unsub-neg51.1%
Simplified51.1%
sin-mult52.6%
Applied egg-rr52.6%
if -2.0000000000000001e-4 < phi1 Initial program 70.9%
Taylor expanded in phi1 around 0 50.6%
sub-neg50.6%
+-commutative50.6%
neg-mul-150.6%
neg-mul-150.6%
remove-double-neg50.6%
mul-1-neg50.6%
distribute-neg-in50.6%
+-commutative50.6%
cos-neg50.6%
+-commutative50.6%
mul-1-neg50.6%
unsub-neg50.6%
Simplified50.6%
Final simplification51.1%
(FPCore (R lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (* (sin phi1) (sin phi2))))
(if (<= phi2 -0.0092)
(* R (acos (+ t_0 (* (cos phi2) (cos phi1)))))
(if (<= phi2 9600000.0)
(*
R
(acos
(+ (* phi2 (sin phi1)) (* (cos phi1) (cos (- lambda2 lambda1))))))
(if (or (<= phi2 3.7e+220) (not (<= phi2 9.5e+269)))
(*
R
(acos (* (cos phi2) (+ (cos lambda2) (* (sin lambda2) lambda1)))))
(* R (acos (+ t_0 (* (cos phi2) (cos lambda1))))))))))
double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = sin(phi1) * sin(phi2);
double tmp;
if (phi2 <= -0.0092) {
tmp = R * acos((t_0 + (cos(phi2) * cos(phi1))));
} else if (phi2 <= 9600000.0) {
tmp = R * acos(((phi2 * sin(phi1)) + (cos(phi1) * cos((lambda2 - lambda1)))));
} else if ((phi2 <= 3.7e+220) || !(phi2 <= 9.5e+269)) {
tmp = R * acos((cos(phi2) * (cos(lambda2) + (sin(lambda2) * lambda1))));
} else {
tmp = R * acos((t_0 + (cos(phi2) * cos(lambda1))));
}
return tmp;
}
real(8) function code(r, lambda1, lambda2, phi1, phi2)
real(8), intent (in) :: r
real(8), intent (in) :: lambda1
real(8), intent (in) :: lambda2
real(8), intent (in) :: phi1
real(8), intent (in) :: phi2
real(8) :: t_0
real(8) :: tmp
t_0 = sin(phi1) * sin(phi2)
if (phi2 <= (-0.0092d0)) then
tmp = r * acos((t_0 + (cos(phi2) * cos(phi1))))
else if (phi2 <= 9600000.0d0) then
tmp = r * acos(((phi2 * sin(phi1)) + (cos(phi1) * cos((lambda2 - lambda1)))))
else if ((phi2 <= 3.7d+220) .or. (.not. (phi2 <= 9.5d+269))) then
tmp = r * acos((cos(phi2) * (cos(lambda2) + (sin(lambda2) * lambda1))))
else
tmp = r * acos((t_0 + (cos(phi2) * cos(lambda1))))
end if
code = tmp
end function
public static double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = Math.sin(phi1) * Math.sin(phi2);
double tmp;
if (phi2 <= -0.0092) {
tmp = R * Math.acos((t_0 + (Math.cos(phi2) * Math.cos(phi1))));
} else if (phi2 <= 9600000.0) {
tmp = R * Math.acos(((phi2 * Math.sin(phi1)) + (Math.cos(phi1) * Math.cos((lambda2 - lambda1)))));
} else if ((phi2 <= 3.7e+220) || !(phi2 <= 9.5e+269)) {
tmp = R * Math.acos((Math.cos(phi2) * (Math.cos(lambda2) + (Math.sin(lambda2) * lambda1))));
} else {
tmp = R * Math.acos((t_0 + (Math.cos(phi2) * Math.cos(lambda1))));
}
return tmp;
}
def code(R, lambda1, lambda2, phi1, phi2): t_0 = math.sin(phi1) * math.sin(phi2) tmp = 0 if phi2 <= -0.0092: tmp = R * math.acos((t_0 + (math.cos(phi2) * math.cos(phi1)))) elif phi2 <= 9600000.0: tmp = R * math.acos(((phi2 * math.sin(phi1)) + (math.cos(phi1) * math.cos((lambda2 - lambda1))))) elif (phi2 <= 3.7e+220) or not (phi2 <= 9.5e+269): tmp = R * math.acos((math.cos(phi2) * (math.cos(lambda2) + (math.sin(lambda2) * lambda1)))) else: tmp = R * math.acos((t_0 + (math.cos(phi2) * math.cos(lambda1)))) return tmp
function code(R, lambda1, lambda2, phi1, phi2) t_0 = Float64(sin(phi1) * sin(phi2)) tmp = 0.0 if (phi2 <= -0.0092) tmp = Float64(R * acos(Float64(t_0 + Float64(cos(phi2) * cos(phi1))))); elseif (phi2 <= 9600000.0) tmp = Float64(R * acos(Float64(Float64(phi2 * sin(phi1)) + Float64(cos(phi1) * cos(Float64(lambda2 - lambda1)))))); elseif ((phi2 <= 3.7e+220) || !(phi2 <= 9.5e+269)) tmp = Float64(R * acos(Float64(cos(phi2) * Float64(cos(lambda2) + Float64(sin(lambda2) * lambda1))))); else tmp = Float64(R * acos(Float64(t_0 + Float64(cos(phi2) * cos(lambda1))))); end return tmp end
function tmp_2 = code(R, lambda1, lambda2, phi1, phi2) t_0 = sin(phi1) * sin(phi2); tmp = 0.0; if (phi2 <= -0.0092) tmp = R * acos((t_0 + (cos(phi2) * cos(phi1)))); elseif (phi2 <= 9600000.0) tmp = R * acos(((phi2 * sin(phi1)) + (cos(phi1) * cos((lambda2 - lambda1))))); elseif ((phi2 <= 3.7e+220) || ~((phi2 <= 9.5e+269))) tmp = R * acos((cos(phi2) * (cos(lambda2) + (sin(lambda2) * lambda1)))); else tmp = R * acos((t_0 + (cos(phi2) * cos(lambda1)))); end tmp_2 = tmp; end
code[R_, lambda1_, lambda2_, phi1_, phi2_] := Block[{t$95$0 = N[(N[Sin[phi1], $MachinePrecision] * N[Sin[phi2], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[phi2, -0.0092], N[(R * N[ArcCos[N[(t$95$0 + N[(N[Cos[phi2], $MachinePrecision] * N[Cos[phi1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], If[LessEqual[phi2, 9600000.0], N[(R * N[ArcCos[N[(N[(phi2 * N[Sin[phi1], $MachinePrecision]), $MachinePrecision] + N[(N[Cos[phi1], $MachinePrecision] * N[Cos[N[(lambda2 - lambda1), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], If[Or[LessEqual[phi2, 3.7e+220], N[Not[LessEqual[phi2, 9.5e+269]], $MachinePrecision]], N[(R * N[ArcCos[N[(N[Cos[phi2], $MachinePrecision] * N[(N[Cos[lambda2], $MachinePrecision] + N[(N[Sin[lambda2], $MachinePrecision] * lambda1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], N[(R * N[ArcCos[N[(t$95$0 + N[(N[Cos[phi2], $MachinePrecision] * N[Cos[lambda1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sin \phi_1 \cdot \sin \phi_2\\
\mathbf{if}\;\phi_2 \leq -0.0092:\\
\;\;\;\;R \cdot \cos^{-1} \left(t_0 + \cos \phi_2 \cdot \cos \phi_1\right)\\
\mathbf{elif}\;\phi_2 \leq 9600000:\\
\;\;\;\;R \cdot \cos^{-1} \left(\phi_2 \cdot \sin \phi_1 + \cos \phi_1 \cdot \cos \left(\lambda_2 - \lambda_1\right)\right)\\
\mathbf{elif}\;\phi_2 \leq 3.7 \cdot 10^{+220} \lor \neg \left(\phi_2 \leq 9.5 \cdot 10^{+269}\right):\\
\;\;\;\;R \cdot \cos^{-1} \left(\cos \phi_2 \cdot \left(\cos \lambda_2 + \sin \lambda_2 \cdot \lambda_1\right)\right)\\
\mathbf{else}:\\
\;\;\;\;R \cdot \cos^{-1} \left(t_0 + \cos \phi_2 \cdot \cos \lambda_1\right)\\
\end{array}
\end{array}
if phi2 < -0.0091999999999999998Initial program 80.6%
Taylor expanded in lambda1 around 0 51.2%
mul-1-neg51.2%
unsub-neg51.2%
cos-neg51.2%
*-commutative51.2%
sin-neg51.2%
distribute-rgt-neg-out51.2%
Simplified51.2%
Taylor expanded in lambda2 around 0 40.3%
if -0.0091999999999999998 < phi2 < 9.6e6Initial program 69.1%
Taylor expanded in phi2 around 0 68.2%
sub-neg68.2%
+-commutative68.2%
neg-mul-168.2%
neg-mul-168.2%
remove-double-neg68.2%
mul-1-neg68.2%
distribute-neg-in68.2%
+-commutative68.2%
cos-neg68.2%
+-commutative68.2%
mul-1-neg68.2%
unsub-neg68.2%
Simplified68.2%
Taylor expanded in phi2 around 0 68.2%
if 9.6e6 < phi2 < 3.7e220 or 9.50000000000000044e269 < phi2 Initial program 77.3%
Taylor expanded in lambda1 around 0 64.1%
mul-1-neg64.1%
unsub-neg64.1%
cos-neg64.1%
*-commutative64.1%
sin-neg64.1%
distribute-rgt-neg-out64.1%
Simplified64.1%
Taylor expanded in phi1 around 0 64.1%
fma-def64.0%
*-commutative64.0%
Simplified64.0%
Taylor expanded in phi1 around 0 40.2%
if 3.7e220 < phi2 < 9.50000000000000044e269Initial program 78.2%
Taylor expanded in lambda2 around 0 49.8%
Taylor expanded in phi1 around 0 41.0%
Final simplification54.7%
(FPCore (R lambda1 lambda2 phi1 phi2)
:precision binary64
(if (<= phi2 9600000.0)
(*
R
(acos (+ (* phi2 (sin phi1)) (* (cos phi1) (cos (- lambda2 lambda1))))))
(if (or (<= phi2 1.4e+219) (not (<= phi2 1.7e+271)))
(* R (acos (* (cos phi2) (+ (cos lambda2) (* (sin lambda2) lambda1)))))
(* R (acos (+ (* (sin phi1) (sin phi2)) (* (cos phi2) (cos lambda1))))))))
double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
double tmp;
if (phi2 <= 9600000.0) {
tmp = R * acos(((phi2 * sin(phi1)) + (cos(phi1) * cos((lambda2 - lambda1)))));
} else if ((phi2 <= 1.4e+219) || !(phi2 <= 1.7e+271)) {
tmp = R * acos((cos(phi2) * (cos(lambda2) + (sin(lambda2) * lambda1))));
} else {
tmp = R * acos(((sin(phi1) * sin(phi2)) + (cos(phi2) * cos(lambda1))));
}
return tmp;
}
real(8) function code(r, lambda1, lambda2, phi1, phi2)
real(8), intent (in) :: r
real(8), intent (in) :: lambda1
real(8), intent (in) :: lambda2
real(8), intent (in) :: phi1
real(8), intent (in) :: phi2
real(8) :: tmp
if (phi2 <= 9600000.0d0) then
tmp = r * acos(((phi2 * sin(phi1)) + (cos(phi1) * cos((lambda2 - lambda1)))))
else if ((phi2 <= 1.4d+219) .or. (.not. (phi2 <= 1.7d+271))) then
tmp = r * acos((cos(phi2) * (cos(lambda2) + (sin(lambda2) * lambda1))))
else
tmp = r * acos(((sin(phi1) * sin(phi2)) + (cos(phi2) * cos(lambda1))))
end if
code = tmp
end function
public static double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
double tmp;
if (phi2 <= 9600000.0) {
tmp = R * Math.acos(((phi2 * Math.sin(phi1)) + (Math.cos(phi1) * Math.cos((lambda2 - lambda1)))));
} else if ((phi2 <= 1.4e+219) || !(phi2 <= 1.7e+271)) {
tmp = R * Math.acos((Math.cos(phi2) * (Math.cos(lambda2) + (Math.sin(lambda2) * lambda1))));
} else {
tmp = R * Math.acos(((Math.sin(phi1) * Math.sin(phi2)) + (Math.cos(phi2) * Math.cos(lambda1))));
}
return tmp;
}
def code(R, lambda1, lambda2, phi1, phi2): tmp = 0 if phi2 <= 9600000.0: tmp = R * math.acos(((phi2 * math.sin(phi1)) + (math.cos(phi1) * math.cos((lambda2 - lambda1))))) elif (phi2 <= 1.4e+219) or not (phi2 <= 1.7e+271): tmp = R * math.acos((math.cos(phi2) * (math.cos(lambda2) + (math.sin(lambda2) * lambda1)))) else: tmp = R * math.acos(((math.sin(phi1) * math.sin(phi2)) + (math.cos(phi2) * math.cos(lambda1)))) return tmp
function code(R, lambda1, lambda2, phi1, phi2) tmp = 0.0 if (phi2 <= 9600000.0) tmp = Float64(R * acos(Float64(Float64(phi2 * sin(phi1)) + Float64(cos(phi1) * cos(Float64(lambda2 - lambda1)))))); elseif ((phi2 <= 1.4e+219) || !(phi2 <= 1.7e+271)) tmp = Float64(R * acos(Float64(cos(phi2) * Float64(cos(lambda2) + Float64(sin(lambda2) * lambda1))))); else tmp = Float64(R * acos(Float64(Float64(sin(phi1) * sin(phi2)) + Float64(cos(phi2) * cos(lambda1))))); end return tmp end
function tmp_2 = code(R, lambda1, lambda2, phi1, phi2) tmp = 0.0; if (phi2 <= 9600000.0) tmp = R * acos(((phi2 * sin(phi1)) + (cos(phi1) * cos((lambda2 - lambda1))))); elseif ((phi2 <= 1.4e+219) || ~((phi2 <= 1.7e+271))) tmp = R * acos((cos(phi2) * (cos(lambda2) + (sin(lambda2) * lambda1)))); else tmp = R * acos(((sin(phi1) * sin(phi2)) + (cos(phi2) * cos(lambda1)))); end tmp_2 = tmp; end
code[R_, lambda1_, lambda2_, phi1_, phi2_] := If[LessEqual[phi2, 9600000.0], N[(R * N[ArcCos[N[(N[(phi2 * N[Sin[phi1], $MachinePrecision]), $MachinePrecision] + N[(N[Cos[phi1], $MachinePrecision] * N[Cos[N[(lambda2 - lambda1), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], If[Or[LessEqual[phi2, 1.4e+219], N[Not[LessEqual[phi2, 1.7e+271]], $MachinePrecision]], N[(R * N[ArcCos[N[(N[Cos[phi2], $MachinePrecision] * N[(N[Cos[lambda2], $MachinePrecision] + N[(N[Sin[lambda2], $MachinePrecision] * lambda1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], N[(R * N[ArcCos[N[(N[(N[Sin[phi1], $MachinePrecision] * N[Sin[phi2], $MachinePrecision]), $MachinePrecision] + N[(N[Cos[phi2], $MachinePrecision] * N[Cos[lambda1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\phi_2 \leq 9600000:\\
\;\;\;\;R \cdot \cos^{-1} \left(\phi_2 \cdot \sin \phi_1 + \cos \phi_1 \cdot \cos \left(\lambda_2 - \lambda_1\right)\right)\\
\mathbf{elif}\;\phi_2 \leq 1.4 \cdot 10^{+219} \lor \neg \left(\phi_2 \leq 1.7 \cdot 10^{+271}\right):\\
\;\;\;\;R \cdot \cos^{-1} \left(\cos \phi_2 \cdot \left(\cos \lambda_2 + \sin \lambda_2 \cdot \lambda_1\right)\right)\\
\mathbf{else}:\\
\;\;\;\;R \cdot \cos^{-1} \left(\sin \phi_1 \cdot \sin \phi_2 + \cos \phi_2 \cdot \cos \lambda_1\right)\\
\end{array}
\end{array}
if phi2 < 9.6e6Initial program 73.0%
Taylor expanded in phi2 around 0 50.9%
sub-neg50.9%
+-commutative50.9%
neg-mul-150.9%
neg-mul-150.9%
remove-double-neg50.9%
mul-1-neg50.9%
distribute-neg-in50.9%
+-commutative50.9%
cos-neg50.9%
+-commutative50.9%
mul-1-neg50.9%
unsub-neg50.9%
Simplified50.9%
Taylor expanded in phi2 around 0 47.1%
if 9.6e6 < phi2 < 1.40000000000000008e219 or 1.70000000000000013e271 < phi2 Initial program 77.3%
Taylor expanded in lambda1 around 0 64.1%
mul-1-neg64.1%
unsub-neg64.1%
cos-neg64.1%
*-commutative64.1%
sin-neg64.1%
distribute-rgt-neg-out64.1%
Simplified64.1%
Taylor expanded in phi1 around 0 64.1%
fma-def64.0%
*-commutative64.0%
Simplified64.0%
Taylor expanded in phi1 around 0 40.2%
if 1.40000000000000008e219 < phi2 < 1.70000000000000013e271Initial program 78.2%
Taylor expanded in lambda2 around 0 49.8%
Taylor expanded in phi1 around 0 41.0%
Final simplification45.6%
(FPCore (R lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (* (sin phi1) (sin phi2))))
(if (<= phi2 9600000.0)
(* R (acos (+ t_0 (* (cos phi1) (cos (- lambda2 lambda1))))))
(if (or (<= phi2 1.4e+221) (not (<= phi2 1.9e+270)))
(* R (acos (* (cos phi2) (+ (cos lambda2) (* (sin lambda2) lambda1)))))
(* R (acos (+ t_0 (* (cos phi2) (cos lambda1)))))))))
double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = sin(phi1) * sin(phi2);
double tmp;
if (phi2 <= 9600000.0) {
tmp = R * acos((t_0 + (cos(phi1) * cos((lambda2 - lambda1)))));
} else if ((phi2 <= 1.4e+221) || !(phi2 <= 1.9e+270)) {
tmp = R * acos((cos(phi2) * (cos(lambda2) + (sin(lambda2) * lambda1))));
} else {
tmp = R * acos((t_0 + (cos(phi2) * cos(lambda1))));
}
return tmp;
}
real(8) function code(r, lambda1, lambda2, phi1, phi2)
real(8), intent (in) :: r
real(8), intent (in) :: lambda1
real(8), intent (in) :: lambda2
real(8), intent (in) :: phi1
real(8), intent (in) :: phi2
real(8) :: t_0
real(8) :: tmp
t_0 = sin(phi1) * sin(phi2)
if (phi2 <= 9600000.0d0) then
tmp = r * acos((t_0 + (cos(phi1) * cos((lambda2 - lambda1)))))
else if ((phi2 <= 1.4d+221) .or. (.not. (phi2 <= 1.9d+270))) then
tmp = r * acos((cos(phi2) * (cos(lambda2) + (sin(lambda2) * lambda1))))
else
tmp = r * acos((t_0 + (cos(phi2) * cos(lambda1))))
end if
code = tmp
end function
public static double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = Math.sin(phi1) * Math.sin(phi2);
double tmp;
if (phi2 <= 9600000.0) {
tmp = R * Math.acos((t_0 + (Math.cos(phi1) * Math.cos((lambda2 - lambda1)))));
} else if ((phi2 <= 1.4e+221) || !(phi2 <= 1.9e+270)) {
tmp = R * Math.acos((Math.cos(phi2) * (Math.cos(lambda2) + (Math.sin(lambda2) * lambda1))));
} else {
tmp = R * Math.acos((t_0 + (Math.cos(phi2) * Math.cos(lambda1))));
}
return tmp;
}
def code(R, lambda1, lambda2, phi1, phi2): t_0 = math.sin(phi1) * math.sin(phi2) tmp = 0 if phi2 <= 9600000.0: tmp = R * math.acos((t_0 + (math.cos(phi1) * math.cos((lambda2 - lambda1))))) elif (phi2 <= 1.4e+221) or not (phi2 <= 1.9e+270): tmp = R * math.acos((math.cos(phi2) * (math.cos(lambda2) + (math.sin(lambda2) * lambda1)))) else: tmp = R * math.acos((t_0 + (math.cos(phi2) * math.cos(lambda1)))) return tmp
function code(R, lambda1, lambda2, phi1, phi2) t_0 = Float64(sin(phi1) * sin(phi2)) tmp = 0.0 if (phi2 <= 9600000.0) tmp = Float64(R * acos(Float64(t_0 + Float64(cos(phi1) * cos(Float64(lambda2 - lambda1)))))); elseif ((phi2 <= 1.4e+221) || !(phi2 <= 1.9e+270)) tmp = Float64(R * acos(Float64(cos(phi2) * Float64(cos(lambda2) + Float64(sin(lambda2) * lambda1))))); else tmp = Float64(R * acos(Float64(t_0 + Float64(cos(phi2) * cos(lambda1))))); end return tmp end
function tmp_2 = code(R, lambda1, lambda2, phi1, phi2) t_0 = sin(phi1) * sin(phi2); tmp = 0.0; if (phi2 <= 9600000.0) tmp = R * acos((t_0 + (cos(phi1) * cos((lambda2 - lambda1))))); elseif ((phi2 <= 1.4e+221) || ~((phi2 <= 1.9e+270))) tmp = R * acos((cos(phi2) * (cos(lambda2) + (sin(lambda2) * lambda1)))); else tmp = R * acos((t_0 + (cos(phi2) * cos(lambda1)))); end tmp_2 = tmp; end
code[R_, lambda1_, lambda2_, phi1_, phi2_] := Block[{t$95$0 = N[(N[Sin[phi1], $MachinePrecision] * N[Sin[phi2], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[phi2, 9600000.0], N[(R * N[ArcCos[N[(t$95$0 + N[(N[Cos[phi1], $MachinePrecision] * N[Cos[N[(lambda2 - lambda1), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], If[Or[LessEqual[phi2, 1.4e+221], N[Not[LessEqual[phi2, 1.9e+270]], $MachinePrecision]], N[(R * N[ArcCos[N[(N[Cos[phi2], $MachinePrecision] * N[(N[Cos[lambda2], $MachinePrecision] + N[(N[Sin[lambda2], $MachinePrecision] * lambda1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], N[(R * N[ArcCos[N[(t$95$0 + N[(N[Cos[phi2], $MachinePrecision] * N[Cos[lambda1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sin \phi_1 \cdot \sin \phi_2\\
\mathbf{if}\;\phi_2 \leq 9600000:\\
\;\;\;\;R \cdot \cos^{-1} \left(t_0 + \cos \phi_1 \cdot \cos \left(\lambda_2 - \lambda_1\right)\right)\\
\mathbf{elif}\;\phi_2 \leq 1.4 \cdot 10^{+221} \lor \neg \left(\phi_2 \leq 1.9 \cdot 10^{+270}\right):\\
\;\;\;\;R \cdot \cos^{-1} \left(\cos \phi_2 \cdot \left(\cos \lambda_2 + \sin \lambda_2 \cdot \lambda_1\right)\right)\\
\mathbf{else}:\\
\;\;\;\;R \cdot \cos^{-1} \left(t_0 + \cos \phi_2 \cdot \cos \lambda_1\right)\\
\end{array}
\end{array}
if phi2 < 9.6e6Initial program 73.0%
Taylor expanded in phi2 around 0 50.9%
sub-neg50.9%
+-commutative50.9%
neg-mul-150.9%
neg-mul-150.9%
remove-double-neg50.9%
mul-1-neg50.9%
distribute-neg-in50.9%
+-commutative50.9%
cos-neg50.9%
+-commutative50.9%
mul-1-neg50.9%
unsub-neg50.9%
Simplified50.9%
if 9.6e6 < phi2 < 1.39999999999999994e221 or 1.90000000000000009e270 < phi2 Initial program 77.3%
Taylor expanded in lambda1 around 0 64.1%
mul-1-neg64.1%
unsub-neg64.1%
cos-neg64.1%
*-commutative64.1%
sin-neg64.1%
distribute-rgt-neg-out64.1%
Simplified64.1%
Taylor expanded in phi1 around 0 64.1%
fma-def64.0%
*-commutative64.0%
Simplified64.0%
Taylor expanded in phi1 around 0 40.2%
if 1.39999999999999994e221 < phi2 < 1.90000000000000009e270Initial program 78.2%
Taylor expanded in lambda2 around 0 49.8%
Taylor expanded in phi1 around 0 41.0%
Final simplification48.5%
(FPCore (R lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (cos (- lambda2 lambda1))) (t_1 (* (sin phi1) (sin phi2))))
(if (<= phi1 -9e-5)
(* R (acos (+ t_1 (* (cos phi1) t_0))))
(* R (acos (+ t_1 (* (cos phi2) t_0)))))))
double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = cos((lambda2 - lambda1));
double t_1 = sin(phi1) * sin(phi2);
double tmp;
if (phi1 <= -9e-5) {
tmp = R * acos((t_1 + (cos(phi1) * t_0)));
} else {
tmp = R * acos((t_1 + (cos(phi2) * t_0)));
}
return tmp;
}
real(8) function code(r, lambda1, lambda2, phi1, phi2)
real(8), intent (in) :: r
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((lambda2 - lambda1))
t_1 = sin(phi1) * sin(phi2)
if (phi1 <= (-9d-5)) then
tmp = r * acos((t_1 + (cos(phi1) * t_0)))
else
tmp = r * acos((t_1 + (cos(phi2) * t_0)))
end if
code = tmp
end function
public static double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = Math.cos((lambda2 - lambda1));
double t_1 = Math.sin(phi1) * Math.sin(phi2);
double tmp;
if (phi1 <= -9e-5) {
tmp = R * Math.acos((t_1 + (Math.cos(phi1) * t_0)));
} else {
tmp = R * Math.acos((t_1 + (Math.cos(phi2) * t_0)));
}
return tmp;
}
def code(R, lambda1, lambda2, phi1, phi2): t_0 = math.cos((lambda2 - lambda1)) t_1 = math.sin(phi1) * math.sin(phi2) tmp = 0 if phi1 <= -9e-5: tmp = R * math.acos((t_1 + (math.cos(phi1) * t_0))) else: tmp = R * math.acos((t_1 + (math.cos(phi2) * t_0))) return tmp
function code(R, lambda1, lambda2, phi1, phi2) t_0 = cos(Float64(lambda2 - lambda1)) t_1 = Float64(sin(phi1) * sin(phi2)) tmp = 0.0 if (phi1 <= -9e-5) tmp = Float64(R * acos(Float64(t_1 + Float64(cos(phi1) * t_0)))); else tmp = Float64(R * acos(Float64(t_1 + Float64(cos(phi2) * t_0)))); end return tmp end
function tmp_2 = code(R, lambda1, lambda2, phi1, phi2) t_0 = cos((lambda2 - lambda1)); t_1 = sin(phi1) * sin(phi2); tmp = 0.0; if (phi1 <= -9e-5) tmp = R * acos((t_1 + (cos(phi1) * t_0))); else tmp = R * acos((t_1 + (cos(phi2) * t_0))); end tmp_2 = tmp; end
code[R_, lambda1_, lambda2_, phi1_, phi2_] := Block[{t$95$0 = N[Cos[N[(lambda2 - lambda1), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$1 = N[(N[Sin[phi1], $MachinePrecision] * N[Sin[phi2], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[phi1, -9e-5], N[(R * N[ArcCos[N[(t$95$1 + N[(N[Cos[phi1], $MachinePrecision] * t$95$0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], N[(R * N[ArcCos[N[(t$95$1 + N[(N[Cos[phi2], $MachinePrecision] * t$95$0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \cos \left(\lambda_2 - \lambda_1\right)\\
t_1 := \sin \phi_1 \cdot \sin \phi_2\\
\mathbf{if}\;\phi_1 \leq -9 \cdot 10^{-5}:\\
\;\;\;\;R \cdot \cos^{-1} \left(t_1 + \cos \phi_1 \cdot t_0\right)\\
\mathbf{else}:\\
\;\;\;\;R \cdot \cos^{-1} \left(t_1 + \cos \phi_2 \cdot t_0\right)\\
\end{array}
\end{array}
if phi1 < -9.00000000000000057e-5Initial program 83.5%
Taylor expanded in phi2 around 0 51.1%
sub-neg51.1%
+-commutative51.1%
neg-mul-151.1%
neg-mul-151.1%
remove-double-neg51.1%
mul-1-neg51.1%
distribute-neg-in51.1%
+-commutative51.1%
cos-neg51.1%
+-commutative51.1%
mul-1-neg51.1%
unsub-neg51.1%
Simplified51.1%
if -9.00000000000000057e-5 < phi1 Initial program 70.9%
Taylor expanded in phi1 around 0 50.6%
sub-neg50.6%
+-commutative50.6%
neg-mul-150.6%
neg-mul-150.6%
remove-double-neg50.6%
mul-1-neg50.6%
distribute-neg-in50.6%
+-commutative50.6%
cos-neg50.6%
+-commutative50.6%
mul-1-neg50.6%
unsub-neg50.6%
Simplified50.6%
Final simplification50.7%
(FPCore (R lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (cos (- lambda2 lambda1)))
(t_1 (+ (cos lambda2) (* (sin lambda2) lambda1)))
(t_2 (* R (acos (* (cos phi1) t_1)))))
(if (<= phi2 -1.62e-133)
t_2
(if (<= phi2 2.9e-84)
(* R (acos (+ (* (cos phi1) t_0) (* phi2 phi1))))
(if (<= phi2 3.9e-15)
t_2
(if (<= phi2 9600000.0)
(* R (acos (+ t_0 (* phi2 phi1))))
(* R (acos (* (cos phi2) t_1)))))))))
double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = cos((lambda2 - lambda1));
double t_1 = cos(lambda2) + (sin(lambda2) * lambda1);
double t_2 = R * acos((cos(phi1) * t_1));
double tmp;
if (phi2 <= -1.62e-133) {
tmp = t_2;
} else if (phi2 <= 2.9e-84) {
tmp = R * acos(((cos(phi1) * t_0) + (phi2 * phi1)));
} else if (phi2 <= 3.9e-15) {
tmp = t_2;
} else if (phi2 <= 9600000.0) {
tmp = R * acos((t_0 + (phi2 * phi1)));
} else {
tmp = R * acos((cos(phi2) * t_1));
}
return tmp;
}
real(8) function code(r, lambda1, lambda2, phi1, phi2)
real(8), intent (in) :: r
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((lambda2 - lambda1))
t_1 = cos(lambda2) + (sin(lambda2) * lambda1)
t_2 = r * acos((cos(phi1) * t_1))
if (phi2 <= (-1.62d-133)) then
tmp = t_2
else if (phi2 <= 2.9d-84) then
tmp = r * acos(((cos(phi1) * t_0) + (phi2 * phi1)))
else if (phi2 <= 3.9d-15) then
tmp = t_2
else if (phi2 <= 9600000.0d0) then
tmp = r * acos((t_0 + (phi2 * phi1)))
else
tmp = r * acos((cos(phi2) * t_1))
end if
code = tmp
end function
public static double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = Math.cos((lambda2 - lambda1));
double t_1 = Math.cos(lambda2) + (Math.sin(lambda2) * lambda1);
double t_2 = R * Math.acos((Math.cos(phi1) * t_1));
double tmp;
if (phi2 <= -1.62e-133) {
tmp = t_2;
} else if (phi2 <= 2.9e-84) {
tmp = R * Math.acos(((Math.cos(phi1) * t_0) + (phi2 * phi1)));
} else if (phi2 <= 3.9e-15) {
tmp = t_2;
} else if (phi2 <= 9600000.0) {
tmp = R * Math.acos((t_0 + (phi2 * phi1)));
} else {
tmp = R * Math.acos((Math.cos(phi2) * t_1));
}
return tmp;
}
def code(R, lambda1, lambda2, phi1, phi2): t_0 = math.cos((lambda2 - lambda1)) t_1 = math.cos(lambda2) + (math.sin(lambda2) * lambda1) t_2 = R * math.acos((math.cos(phi1) * t_1)) tmp = 0 if phi2 <= -1.62e-133: tmp = t_2 elif phi2 <= 2.9e-84: tmp = R * math.acos(((math.cos(phi1) * t_0) + (phi2 * phi1))) elif phi2 <= 3.9e-15: tmp = t_2 elif phi2 <= 9600000.0: tmp = R * math.acos((t_0 + (phi2 * phi1))) else: tmp = R * math.acos((math.cos(phi2) * t_1)) return tmp
function code(R, lambda1, lambda2, phi1, phi2) t_0 = cos(Float64(lambda2 - lambda1)) t_1 = Float64(cos(lambda2) + Float64(sin(lambda2) * lambda1)) t_2 = Float64(R * acos(Float64(cos(phi1) * t_1))) tmp = 0.0 if (phi2 <= -1.62e-133) tmp = t_2; elseif (phi2 <= 2.9e-84) tmp = Float64(R * acos(Float64(Float64(cos(phi1) * t_0) + Float64(phi2 * phi1)))); elseif (phi2 <= 3.9e-15) tmp = t_2; elseif (phi2 <= 9600000.0) tmp = Float64(R * acos(Float64(t_0 + Float64(phi2 * phi1)))); else tmp = Float64(R * acos(Float64(cos(phi2) * t_1))); end return tmp end
function tmp_2 = code(R, lambda1, lambda2, phi1, phi2) t_0 = cos((lambda2 - lambda1)); t_1 = cos(lambda2) + (sin(lambda2) * lambda1); t_2 = R * acos((cos(phi1) * t_1)); tmp = 0.0; if (phi2 <= -1.62e-133) tmp = t_2; elseif (phi2 <= 2.9e-84) tmp = R * acos(((cos(phi1) * t_0) + (phi2 * phi1))); elseif (phi2 <= 3.9e-15) tmp = t_2; elseif (phi2 <= 9600000.0) tmp = R * acos((t_0 + (phi2 * phi1))); else tmp = R * acos((cos(phi2) * t_1)); end tmp_2 = tmp; end
code[R_, lambda1_, lambda2_, phi1_, phi2_] := Block[{t$95$0 = N[Cos[N[(lambda2 - lambda1), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$1 = N[(N[Cos[lambda2], $MachinePrecision] + N[(N[Sin[lambda2], $MachinePrecision] * lambda1), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(R * N[ArcCos[N[(N[Cos[phi1], $MachinePrecision] * t$95$1), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[phi2, -1.62e-133], t$95$2, If[LessEqual[phi2, 2.9e-84], N[(R * N[ArcCos[N[(N[(N[Cos[phi1], $MachinePrecision] * t$95$0), $MachinePrecision] + N[(phi2 * phi1), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], If[LessEqual[phi2, 3.9e-15], t$95$2, If[LessEqual[phi2, 9600000.0], N[(R * N[ArcCos[N[(t$95$0 + N[(phi2 * phi1), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], N[(R * N[ArcCos[N[(N[Cos[phi2], $MachinePrecision] * t$95$1), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \cos \left(\lambda_2 - \lambda_1\right)\\
t_1 := \cos \lambda_2 + \sin \lambda_2 \cdot \lambda_1\\
t_2 := R \cdot \cos^{-1} \left(\cos \phi_1 \cdot t_1\right)\\
\mathbf{if}\;\phi_2 \leq -1.62 \cdot 10^{-133}:\\
\;\;\;\;t_2\\
\mathbf{elif}\;\phi_2 \leq 2.9 \cdot 10^{-84}:\\
\;\;\;\;R \cdot \cos^{-1} \left(\cos \phi_1 \cdot t_0 + \phi_2 \cdot \phi_1\right)\\
\mathbf{elif}\;\phi_2 \leq 3.9 \cdot 10^{-15}:\\
\;\;\;\;t_2\\
\mathbf{elif}\;\phi_2 \leq 9600000:\\
\;\;\;\;R \cdot \cos^{-1} \left(t_0 + \phi_2 \cdot \phi_1\right)\\
\mathbf{else}:\\
\;\;\;\;R \cdot \cos^{-1} \left(\cos \phi_2 \cdot t_1\right)\\
\end{array}
\end{array}
if phi2 < -1.61999999999999998e-133 or 2.90000000000000019e-84 < phi2 < 3.90000000000000026e-15Initial program 74.7%
Taylor expanded in lambda1 around 0 45.7%
mul-1-neg45.7%
unsub-neg45.7%
cos-neg45.7%
*-commutative45.7%
sin-neg45.7%
distribute-rgt-neg-out45.7%
Simplified45.7%
Taylor expanded in phi1 around 0 45.7%
fma-def45.7%
*-commutative45.7%
Simplified45.7%
Taylor expanded in phi2 around 0 19.8%
if -1.61999999999999998e-133 < phi2 < 2.90000000000000019e-84Initial program 71.9%
Taylor expanded in phi2 around 0 71.9%
sub-neg71.9%
+-commutative71.9%
neg-mul-171.9%
neg-mul-171.9%
remove-double-neg71.9%
mul-1-neg71.9%
distribute-neg-in71.9%
+-commutative71.9%
cos-neg71.9%
+-commutative71.9%
mul-1-neg71.9%
unsub-neg71.9%
Simplified71.9%
Taylor expanded in phi1 around 0 61.2%
Taylor expanded in phi2 around 0 61.2%
*-commutative61.2%
Simplified61.2%
if 3.90000000000000026e-15 < phi2 < 9.6e6Initial program 54.1%
Taylor expanded in phi2 around 0 30.4%
sub-neg30.4%
+-commutative30.4%
neg-mul-130.4%
neg-mul-130.4%
remove-double-neg30.4%
mul-1-neg30.4%
distribute-neg-in30.4%
+-commutative30.4%
cos-neg30.4%
+-commutative30.4%
mul-1-neg30.4%
unsub-neg30.4%
Simplified30.4%
Taylor expanded in phi1 around 0 4.9%
Taylor expanded in phi2 around 0 4.9%
*-commutative4.9%
Simplified4.9%
Taylor expanded in phi1 around 0 4.9%
if 9.6e6 < phi2 Initial program 77.4%
Taylor expanded in lambda1 around 0 62.5%
mul-1-neg62.5%
unsub-neg62.5%
cos-neg62.5%
*-commutative62.5%
sin-neg62.5%
distribute-rgt-neg-out62.5%
Simplified62.5%
Taylor expanded in phi1 around 0 62.5%
fma-def62.5%
*-commutative62.5%
Simplified62.5%
Taylor expanded in phi1 around 0 40.4%
Final simplification37.8%
(FPCore (R lambda1 lambda2 phi1 phi2)
:precision binary64
(if (<= phi2 9600000.0)
(*
R
(acos (+ (* phi2 (sin phi1)) (* (cos phi1) (cos (- lambda2 lambda1))))))
(* R (acos (* (cos phi2) (+ (cos lambda2) (* (sin lambda2) lambda1)))))))
double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
double tmp;
if (phi2 <= 9600000.0) {
tmp = R * acos(((phi2 * sin(phi1)) + (cos(phi1) * cos((lambda2 - lambda1)))));
} else {
tmp = R * acos((cos(phi2) * (cos(lambda2) + (sin(lambda2) * lambda1))));
}
return tmp;
}
real(8) function code(r, lambda1, lambda2, phi1, phi2)
real(8), intent (in) :: r
real(8), intent (in) :: lambda1
real(8), intent (in) :: lambda2
real(8), intent (in) :: phi1
real(8), intent (in) :: phi2
real(8) :: tmp
if (phi2 <= 9600000.0d0) then
tmp = r * acos(((phi2 * sin(phi1)) + (cos(phi1) * cos((lambda2 - lambda1)))))
else
tmp = r * acos((cos(phi2) * (cos(lambda2) + (sin(lambda2) * lambda1))))
end if
code = tmp
end function
public static double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
double tmp;
if (phi2 <= 9600000.0) {
tmp = R * Math.acos(((phi2 * Math.sin(phi1)) + (Math.cos(phi1) * Math.cos((lambda2 - lambda1)))));
} else {
tmp = R * Math.acos((Math.cos(phi2) * (Math.cos(lambda2) + (Math.sin(lambda2) * lambda1))));
}
return tmp;
}
def code(R, lambda1, lambda2, phi1, phi2): tmp = 0 if phi2 <= 9600000.0: tmp = R * math.acos(((phi2 * math.sin(phi1)) + (math.cos(phi1) * math.cos((lambda2 - lambda1))))) else: tmp = R * math.acos((math.cos(phi2) * (math.cos(lambda2) + (math.sin(lambda2) * lambda1)))) return tmp
function code(R, lambda1, lambda2, phi1, phi2) tmp = 0.0 if (phi2 <= 9600000.0) tmp = Float64(R * acos(Float64(Float64(phi2 * sin(phi1)) + Float64(cos(phi1) * cos(Float64(lambda2 - lambda1)))))); else tmp = Float64(R * acos(Float64(cos(phi2) * Float64(cos(lambda2) + Float64(sin(lambda2) * lambda1))))); end return tmp end
function tmp_2 = code(R, lambda1, lambda2, phi1, phi2) tmp = 0.0; if (phi2 <= 9600000.0) tmp = R * acos(((phi2 * sin(phi1)) + (cos(phi1) * cos((lambda2 - lambda1))))); else tmp = R * acos((cos(phi2) * (cos(lambda2) + (sin(lambda2) * lambda1)))); end tmp_2 = tmp; end
code[R_, lambda1_, lambda2_, phi1_, phi2_] := If[LessEqual[phi2, 9600000.0], N[(R * N[ArcCos[N[(N[(phi2 * N[Sin[phi1], $MachinePrecision]), $MachinePrecision] + N[(N[Cos[phi1], $MachinePrecision] * N[Cos[N[(lambda2 - lambda1), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], N[(R * N[ArcCos[N[(N[Cos[phi2], $MachinePrecision] * N[(N[Cos[lambda2], $MachinePrecision] + N[(N[Sin[lambda2], $MachinePrecision] * lambda1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\phi_2 \leq 9600000:\\
\;\;\;\;R \cdot \cos^{-1} \left(\phi_2 \cdot \sin \phi_1 + \cos \phi_1 \cdot \cos \left(\lambda_2 - \lambda_1\right)\right)\\
\mathbf{else}:\\
\;\;\;\;R \cdot \cos^{-1} \left(\cos \phi_2 \cdot \left(\cos \lambda_2 + \sin \lambda_2 \cdot \lambda_1\right)\right)\\
\end{array}
\end{array}
if phi2 < 9.6e6Initial program 73.0%
Taylor expanded in phi2 around 0 50.9%
sub-neg50.9%
+-commutative50.9%
neg-mul-150.9%
neg-mul-150.9%
remove-double-neg50.9%
mul-1-neg50.9%
distribute-neg-in50.9%
+-commutative50.9%
cos-neg50.9%
+-commutative50.9%
mul-1-neg50.9%
unsub-neg50.9%
Simplified50.9%
Taylor expanded in phi2 around 0 47.1%
if 9.6e6 < phi2 Initial program 77.4%
Taylor expanded in lambda1 around 0 62.5%
mul-1-neg62.5%
unsub-neg62.5%
cos-neg62.5%
*-commutative62.5%
sin-neg62.5%
distribute-rgt-neg-out62.5%
Simplified62.5%
Taylor expanded in phi1 around 0 62.5%
fma-def62.5%
*-commutative62.5%
Simplified62.5%
Taylor expanded in phi1 around 0 40.4%
Final simplification45.6%
(FPCore (R lambda1 lambda2 phi1 phi2) :precision binary64 (if (<= lambda1 -0.00018) (* R (acos (+ (* phi2 phi1) (* (cos phi1) (cos lambda1))))) (* R (acos (* (cos phi1) (+ (cos lambda2) (* (sin lambda2) lambda1)))))))
double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
double tmp;
if (lambda1 <= -0.00018) {
tmp = R * acos(((phi2 * phi1) + (cos(phi1) * cos(lambda1))));
} else {
tmp = R * acos((cos(phi1) * (cos(lambda2) + (sin(lambda2) * lambda1))));
}
return tmp;
}
real(8) function code(r, lambda1, lambda2, phi1, phi2)
real(8), intent (in) :: r
real(8), intent (in) :: lambda1
real(8), intent (in) :: lambda2
real(8), intent (in) :: phi1
real(8), intent (in) :: phi2
real(8) :: tmp
if (lambda1 <= (-0.00018d0)) then
tmp = r * acos(((phi2 * phi1) + (cos(phi1) * cos(lambda1))))
else
tmp = r * acos((cos(phi1) * (cos(lambda2) + (sin(lambda2) * lambda1))))
end if
code = tmp
end function
public static double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
double tmp;
if (lambda1 <= -0.00018) {
tmp = R * Math.acos(((phi2 * phi1) + (Math.cos(phi1) * Math.cos(lambda1))));
} else {
tmp = R * Math.acos((Math.cos(phi1) * (Math.cos(lambda2) + (Math.sin(lambda2) * lambda1))));
}
return tmp;
}
def code(R, lambda1, lambda2, phi1, phi2): tmp = 0 if lambda1 <= -0.00018: tmp = R * math.acos(((phi2 * phi1) + (math.cos(phi1) * math.cos(lambda1)))) else: tmp = R * math.acos((math.cos(phi1) * (math.cos(lambda2) + (math.sin(lambda2) * lambda1)))) return tmp
function code(R, lambda1, lambda2, phi1, phi2) tmp = 0.0 if (lambda1 <= -0.00018) tmp = Float64(R * acos(Float64(Float64(phi2 * phi1) + Float64(cos(phi1) * cos(lambda1))))); else tmp = Float64(R * acos(Float64(cos(phi1) * Float64(cos(lambda2) + Float64(sin(lambda2) * lambda1))))); end return tmp end
function tmp_2 = code(R, lambda1, lambda2, phi1, phi2) tmp = 0.0; if (lambda1 <= -0.00018) tmp = R * acos(((phi2 * phi1) + (cos(phi1) * cos(lambda1)))); else tmp = R * acos((cos(phi1) * (cos(lambda2) + (sin(lambda2) * lambda1)))); end tmp_2 = tmp; end
code[R_, lambda1_, lambda2_, phi1_, phi2_] := If[LessEqual[lambda1, -0.00018], N[(R * N[ArcCos[N[(N[(phi2 * phi1), $MachinePrecision] + N[(N[Cos[phi1], $MachinePrecision] * N[Cos[lambda1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], N[(R * N[ArcCos[N[(N[Cos[phi1], $MachinePrecision] * N[(N[Cos[lambda2], $MachinePrecision] + N[(N[Sin[lambda2], $MachinePrecision] * lambda1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\lambda_1 \leq -0.00018:\\
\;\;\;\;R \cdot \cos^{-1} \left(\phi_2 \cdot \phi_1 + \cos \phi_1 \cdot \cos \lambda_1\right)\\
\mathbf{else}:\\
\;\;\;\;R \cdot \cos^{-1} \left(\cos \phi_1 \cdot \left(\cos \lambda_2 + \sin \lambda_2 \cdot \lambda_1\right)\right)\\
\end{array}
\end{array}
if lambda1 < -1.80000000000000011e-4Initial program 53.9%
Taylor expanded in phi2 around 0 37.4%
sub-neg37.4%
+-commutative37.4%
neg-mul-137.4%
neg-mul-137.4%
remove-double-neg37.4%
mul-1-neg37.4%
distribute-neg-in37.4%
+-commutative37.4%
cos-neg37.4%
+-commutative37.4%
mul-1-neg37.4%
unsub-neg37.4%
Simplified37.4%
Taylor expanded in phi1 around 0 21.0%
Taylor expanded in phi2 around 0 19.7%
*-commutative19.7%
Simplified19.7%
Taylor expanded in lambda2 around 0 19.9%
cos-neg19.9%
*-commutative19.9%
Simplified19.9%
if -1.80000000000000011e-4 < lambda1 Initial program 79.8%
Taylor expanded in lambda1 around 0 61.7%
mul-1-neg61.7%
unsub-neg61.7%
cos-neg61.7%
*-commutative61.7%
sin-neg61.7%
distribute-rgt-neg-out61.7%
Simplified61.7%
Taylor expanded in phi1 around 0 61.7%
fma-def61.7%
*-commutative61.7%
Simplified61.7%
Taylor expanded in phi2 around 0 33.0%
Final simplification30.1%
(FPCore (R lambda1 lambda2 phi1 phi2)
:precision binary64
(if (<= phi1 -0.036)
(* R (acos (+ (* phi2 phi1) (* (cos phi1) (cos lambda1)))))
(*
R
(acos
(+
(* phi2 phi1)
(* (+ (* -0.5 (* phi1 phi1)) 1.0) (cos (- lambda2 lambda1))))))))
double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
double tmp;
if (phi1 <= -0.036) {
tmp = R * acos(((phi2 * phi1) + (cos(phi1) * cos(lambda1))));
} else {
tmp = R * acos(((phi2 * phi1) + (((-0.5 * (phi1 * phi1)) + 1.0) * cos((lambda2 - lambda1)))));
}
return tmp;
}
real(8) function code(r, lambda1, lambda2, phi1, phi2)
real(8), intent (in) :: r
real(8), intent (in) :: lambda1
real(8), intent (in) :: lambda2
real(8), intent (in) :: phi1
real(8), intent (in) :: phi2
real(8) :: tmp
if (phi1 <= (-0.036d0)) then
tmp = r * acos(((phi2 * phi1) + (cos(phi1) * cos(lambda1))))
else
tmp = r * acos(((phi2 * phi1) + ((((-0.5d0) * (phi1 * phi1)) + 1.0d0) * cos((lambda2 - lambda1)))))
end if
code = tmp
end function
public static double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
double tmp;
if (phi1 <= -0.036) {
tmp = R * Math.acos(((phi2 * phi1) + (Math.cos(phi1) * Math.cos(lambda1))));
} else {
tmp = R * Math.acos(((phi2 * phi1) + (((-0.5 * (phi1 * phi1)) + 1.0) * Math.cos((lambda2 - lambda1)))));
}
return tmp;
}
def code(R, lambda1, lambda2, phi1, phi2): tmp = 0 if phi1 <= -0.036: tmp = R * math.acos(((phi2 * phi1) + (math.cos(phi1) * math.cos(lambda1)))) else: tmp = R * math.acos(((phi2 * phi1) + (((-0.5 * (phi1 * phi1)) + 1.0) * math.cos((lambda2 - lambda1))))) return tmp
function code(R, lambda1, lambda2, phi1, phi2) tmp = 0.0 if (phi1 <= -0.036) tmp = Float64(R * acos(Float64(Float64(phi2 * phi1) + Float64(cos(phi1) * cos(lambda1))))); else tmp = Float64(R * acos(Float64(Float64(phi2 * phi1) + Float64(Float64(Float64(-0.5 * Float64(phi1 * phi1)) + 1.0) * cos(Float64(lambda2 - lambda1)))))); end return tmp end
function tmp_2 = code(R, lambda1, lambda2, phi1, phi2) tmp = 0.0; if (phi1 <= -0.036) tmp = R * acos(((phi2 * phi1) + (cos(phi1) * cos(lambda1)))); else tmp = R * acos(((phi2 * phi1) + (((-0.5 * (phi1 * phi1)) + 1.0) * cos((lambda2 - lambda1))))); end tmp_2 = tmp; end
code[R_, lambda1_, lambda2_, phi1_, phi2_] := If[LessEqual[phi1, -0.036], N[(R * N[ArcCos[N[(N[(phi2 * phi1), $MachinePrecision] + N[(N[Cos[phi1], $MachinePrecision] * N[Cos[lambda1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], N[(R * N[ArcCos[N[(N[(phi2 * phi1), $MachinePrecision] + N[(N[(N[(-0.5 * N[(phi1 * phi1), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision] * N[Cos[N[(lambda2 - lambda1), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\phi_1 \leq -0.036:\\
\;\;\;\;R \cdot \cos^{-1} \left(\phi_2 \cdot \phi_1 + \cos \phi_1 \cdot \cos \lambda_1\right)\\
\mathbf{else}:\\
\;\;\;\;R \cdot \cos^{-1} \left(\phi_2 \cdot \phi_1 + \left(-0.5 \cdot \left(\phi_1 \cdot \phi_1\right) + 1\right) \cdot \cos \left(\lambda_2 - \lambda_1\right)\right)\\
\end{array}
\end{array}
if phi1 < -0.0359999999999999973Initial program 84.5%
Taylor expanded in phi2 around 0 51.5%
sub-neg51.5%
+-commutative51.5%
neg-mul-151.5%
neg-mul-151.5%
remove-double-neg51.5%
mul-1-neg51.5%
distribute-neg-in51.5%
+-commutative51.5%
cos-neg51.5%
+-commutative51.5%
mul-1-neg51.5%
unsub-neg51.5%
Simplified51.5%
Taylor expanded in phi1 around 0 21.6%
Taylor expanded in phi2 around 0 21.6%
*-commutative21.6%
Simplified21.6%
Taylor expanded in lambda2 around 0 16.4%
cos-neg16.4%
*-commutative16.4%
Simplified16.4%
if -0.0359999999999999973 < phi1 Initial program 70.6%
Taylor expanded in phi2 around 0 40.7%
sub-neg40.7%
+-commutative40.7%
neg-mul-140.7%
neg-mul-140.7%
remove-double-neg40.7%
mul-1-neg40.7%
distribute-neg-in40.7%
+-commutative40.7%
cos-neg40.7%
+-commutative40.7%
mul-1-neg40.7%
unsub-neg40.7%
Simplified40.7%
Taylor expanded in phi1 around 0 31.3%
Taylor expanded in phi2 around 0 29.0%
*-commutative29.0%
Simplified29.0%
Taylor expanded in phi1 around 0 21.6%
associate-*r*21.6%
distribute-lft1-in21.6%
unpow221.6%
Simplified21.6%
Final simplification20.4%
(FPCore (R lambda1 lambda2 phi1 phi2) :precision binary64 (if (<= lambda1 -0.000145) (* R (acos (+ (* phi2 phi1) (* (cos phi1) (cos lambda1))))) (* R (acos (+ (* phi2 phi1) (* (cos phi1) (cos lambda2)))))))
double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
double tmp;
if (lambda1 <= -0.000145) {
tmp = R * acos(((phi2 * phi1) + (cos(phi1) * cos(lambda1))));
} else {
tmp = R * acos(((phi2 * phi1) + (cos(phi1) * cos(lambda2))));
}
return tmp;
}
real(8) function code(r, lambda1, lambda2, phi1, phi2)
real(8), intent (in) :: r
real(8), intent (in) :: lambda1
real(8), intent (in) :: lambda2
real(8), intent (in) :: phi1
real(8), intent (in) :: phi2
real(8) :: tmp
if (lambda1 <= (-0.000145d0)) then
tmp = r * acos(((phi2 * phi1) + (cos(phi1) * cos(lambda1))))
else
tmp = r * acos(((phi2 * phi1) + (cos(phi1) * cos(lambda2))))
end if
code = tmp
end function
public static double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
double tmp;
if (lambda1 <= -0.000145) {
tmp = R * Math.acos(((phi2 * phi1) + (Math.cos(phi1) * Math.cos(lambda1))));
} else {
tmp = R * Math.acos(((phi2 * phi1) + (Math.cos(phi1) * Math.cos(lambda2))));
}
return tmp;
}
def code(R, lambda1, lambda2, phi1, phi2): tmp = 0 if lambda1 <= -0.000145: tmp = R * math.acos(((phi2 * phi1) + (math.cos(phi1) * math.cos(lambda1)))) else: tmp = R * math.acos(((phi2 * phi1) + (math.cos(phi1) * math.cos(lambda2)))) return tmp
function code(R, lambda1, lambda2, phi1, phi2) tmp = 0.0 if (lambda1 <= -0.000145) tmp = Float64(R * acos(Float64(Float64(phi2 * phi1) + Float64(cos(phi1) * cos(lambda1))))); else tmp = Float64(R * acos(Float64(Float64(phi2 * phi1) + Float64(cos(phi1) * cos(lambda2))))); end return tmp end
function tmp_2 = code(R, lambda1, lambda2, phi1, phi2) tmp = 0.0; if (lambda1 <= -0.000145) tmp = R * acos(((phi2 * phi1) + (cos(phi1) * cos(lambda1)))); else tmp = R * acos(((phi2 * phi1) + (cos(phi1) * cos(lambda2)))); end tmp_2 = tmp; end
code[R_, lambda1_, lambda2_, phi1_, phi2_] := If[LessEqual[lambda1, -0.000145], N[(R * N[ArcCos[N[(N[(phi2 * phi1), $MachinePrecision] + N[(N[Cos[phi1], $MachinePrecision] * N[Cos[lambda1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], N[(R * N[ArcCos[N[(N[(phi2 * phi1), $MachinePrecision] + N[(N[Cos[phi1], $MachinePrecision] * N[Cos[lambda2], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\lambda_1 \leq -0.000145:\\
\;\;\;\;R \cdot \cos^{-1} \left(\phi_2 \cdot \phi_1 + \cos \phi_1 \cdot \cos \lambda_1\right)\\
\mathbf{else}:\\
\;\;\;\;R \cdot \cos^{-1} \left(\phi_2 \cdot \phi_1 + \cos \phi_1 \cdot \cos \lambda_2\right)\\
\end{array}
\end{array}
if lambda1 < -1.45e-4Initial program 53.9%
Taylor expanded in phi2 around 0 37.4%
sub-neg37.4%
+-commutative37.4%
neg-mul-137.4%
neg-mul-137.4%
remove-double-neg37.4%
mul-1-neg37.4%
distribute-neg-in37.4%
+-commutative37.4%
cos-neg37.4%
+-commutative37.4%
mul-1-neg37.4%
unsub-neg37.4%
Simplified37.4%
Taylor expanded in phi1 around 0 21.0%
Taylor expanded in phi2 around 0 19.7%
*-commutative19.7%
Simplified19.7%
Taylor expanded in lambda2 around 0 19.9%
cos-neg19.9%
*-commutative19.9%
Simplified19.9%
if -1.45e-4 < lambda1 Initial program 79.8%
Taylor expanded in phi2 around 0 45.0%
sub-neg45.0%
+-commutative45.0%
neg-mul-145.0%
neg-mul-145.0%
remove-double-neg45.0%
mul-1-neg45.0%
distribute-neg-in45.0%
+-commutative45.0%
cos-neg45.0%
+-commutative45.0%
mul-1-neg45.0%
unsub-neg45.0%
Simplified45.0%
Taylor expanded in phi1 around 0 31.2%
Taylor expanded in phi2 around 0 29.4%
*-commutative29.4%
Simplified29.4%
Taylor expanded in lambda1 around 0 23.3%
Final simplification22.5%
(FPCore (R lambda1 lambda2 phi1 phi2) :precision binary64 (* R (acos (+ (* (cos phi1) (cos (- lambda2 lambda1))) (* phi2 phi1)))))
double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
return R * acos(((cos(phi1) * cos((lambda2 - lambda1))) + (phi2 * phi1)));
}
real(8) function code(r, lambda1, lambda2, phi1, phi2)
real(8), intent (in) :: r
real(8), intent (in) :: lambda1
real(8), intent (in) :: lambda2
real(8), intent (in) :: phi1
real(8), intent (in) :: phi2
code = r * acos(((cos(phi1) * cos((lambda2 - lambda1))) + (phi2 * phi1)))
end function
public static double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
return R * Math.acos(((Math.cos(phi1) * Math.cos((lambda2 - lambda1))) + (phi2 * phi1)));
}
def code(R, lambda1, lambda2, phi1, phi2): return R * math.acos(((math.cos(phi1) * math.cos((lambda2 - lambda1))) + (phi2 * phi1)))
function code(R, lambda1, lambda2, phi1, phi2) return Float64(R * acos(Float64(Float64(cos(phi1) * cos(Float64(lambda2 - lambda1))) + Float64(phi2 * phi1)))) end
function tmp = code(R, lambda1, lambda2, phi1, phi2) tmp = R * acos(((cos(phi1) * cos((lambda2 - lambda1))) + (phi2 * phi1))); end
code[R_, lambda1_, lambda2_, phi1_, phi2_] := N[(R * N[ArcCos[N[(N[(N[Cos[phi1], $MachinePrecision] * N[Cos[N[(lambda2 - lambda1), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] + N[(phi2 * phi1), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
R \cdot \cos^{-1} \left(\cos \phi_1 \cdot \cos \left(\lambda_2 - \lambda_1\right) + \phi_2 \cdot \phi_1\right)
\end{array}
Initial program 74.0%
Taylor expanded in phi2 around 0 43.3%
sub-neg43.3%
+-commutative43.3%
neg-mul-143.3%
neg-mul-143.3%
remove-double-neg43.3%
mul-1-neg43.3%
distribute-neg-in43.3%
+-commutative43.3%
cos-neg43.3%
+-commutative43.3%
mul-1-neg43.3%
unsub-neg43.3%
Simplified43.3%
Taylor expanded in phi1 around 0 28.9%
Taylor expanded in phi2 around 0 27.2%
*-commutative27.2%
Simplified27.2%
Final simplification27.2%
(FPCore (R lambda1 lambda2 phi1 phi2) :precision binary64 (* R (acos (+ (cos (- lambda2 lambda1)) (* phi2 phi1)))))
double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
return R * acos((cos((lambda2 - lambda1)) + (phi2 * phi1)));
}
real(8) function code(r, lambda1, lambda2, phi1, phi2)
real(8), intent (in) :: r
real(8), intent (in) :: lambda1
real(8), intent (in) :: lambda2
real(8), intent (in) :: phi1
real(8), intent (in) :: phi2
code = r * acos((cos((lambda2 - lambda1)) + (phi2 * phi1)))
end function
public static double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
return R * Math.acos((Math.cos((lambda2 - lambda1)) + (phi2 * phi1)));
}
def code(R, lambda1, lambda2, phi1, phi2): return R * math.acos((math.cos((lambda2 - lambda1)) + (phi2 * phi1)))
function code(R, lambda1, lambda2, phi1, phi2) return Float64(R * acos(Float64(cos(Float64(lambda2 - lambda1)) + Float64(phi2 * phi1)))) end
function tmp = code(R, lambda1, lambda2, phi1, phi2) tmp = R * acos((cos((lambda2 - lambda1)) + (phi2 * phi1))); end
code[R_, lambda1_, lambda2_, phi1_, phi2_] := N[(R * N[ArcCos[N[(N[Cos[N[(lambda2 - lambda1), $MachinePrecision]], $MachinePrecision] + N[(phi2 * phi1), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
R \cdot \cos^{-1} \left(\cos \left(\lambda_2 - \lambda_1\right) + \phi_2 \cdot \phi_1\right)
\end{array}
Initial program 74.0%
Taylor expanded in phi2 around 0 43.3%
sub-neg43.3%
+-commutative43.3%
neg-mul-143.3%
neg-mul-143.3%
remove-double-neg43.3%
mul-1-neg43.3%
distribute-neg-in43.3%
+-commutative43.3%
cos-neg43.3%
+-commutative43.3%
mul-1-neg43.3%
unsub-neg43.3%
Simplified43.3%
Taylor expanded in phi1 around 0 28.9%
Taylor expanded in phi2 around 0 27.2%
*-commutative27.2%
Simplified27.2%
Taylor expanded in phi1 around 0 18.5%
Final simplification18.5%
herbie shell --seed 2023207
(FPCore (R lambda1 lambda2 phi1 phi2)
:name "Spherical law of cosines"
:precision binary64
(* (acos (+ (* (sin phi1) (sin phi2)) (* (* (cos phi1) (cos phi2)) (cos (- lambda1 lambda2))))) R))