
(FPCore (lambda1 lambda2 phi1 phi2) :precision binary64 (+ lambda1 (atan2 (* (cos phi2) (sin (- lambda1 lambda2))) (+ (cos phi1) (* (cos phi2) (cos (- lambda1 lambda2)))))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
return lambda1 + atan2((cos(phi2) * sin((lambda1 - lambda2))), (cos(phi1) + (cos(phi2) * cos((lambda1 - lambda2)))));
}
real(8) function code(lambda1, lambda2, phi1, phi2)
real(8), intent (in) :: lambda1
real(8), intent (in) :: lambda2
real(8), intent (in) :: phi1
real(8), intent (in) :: phi2
code = lambda1 + atan2((cos(phi2) * sin((lambda1 - lambda2))), (cos(phi1) + (cos(phi2) * cos((lambda1 - lambda2)))))
end function
public static double code(double lambda1, double lambda2, double phi1, double phi2) {
return lambda1 + Math.atan2((Math.cos(phi2) * Math.sin((lambda1 - lambda2))), (Math.cos(phi1) + (Math.cos(phi2) * Math.cos((lambda1 - lambda2)))));
}
def code(lambda1, lambda2, phi1, phi2): return lambda1 + math.atan2((math.cos(phi2) * math.sin((lambda1 - lambda2))), (math.cos(phi1) + (math.cos(phi2) * math.cos((lambda1 - lambda2)))))
function code(lambda1, lambda2, phi1, phi2) return Float64(lambda1 + atan(Float64(cos(phi2) * sin(Float64(lambda1 - lambda2))), Float64(cos(phi1) + Float64(cos(phi2) * cos(Float64(lambda1 - lambda2)))))) end
function tmp = code(lambda1, lambda2, phi1, phi2) tmp = lambda1 + atan2((cos(phi2) * sin((lambda1 - lambda2))), (cos(phi1) + (cos(phi2) * cos((lambda1 - lambda2))))); end
code[lambda1_, lambda2_, phi1_, phi2_] := N[(lambda1 + N[ArcTan[N[(N[Cos[phi2], $MachinePrecision] * N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / N[(N[Cos[phi1], $MachinePrecision] + N[(N[Cos[phi2], $MachinePrecision] * N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 21 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (lambda1 lambda2 phi1 phi2) :precision binary64 (+ lambda1 (atan2 (* (cos phi2) (sin (- lambda1 lambda2))) (+ (cos phi1) (* (cos phi2) (cos (- lambda1 lambda2)))))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
return lambda1 + atan2((cos(phi2) * sin((lambda1 - lambda2))), (cos(phi1) + (cos(phi2) * cos((lambda1 - lambda2)))));
}
real(8) function code(lambda1, lambda2, phi1, phi2)
real(8), intent (in) :: lambda1
real(8), intent (in) :: lambda2
real(8), intent (in) :: phi1
real(8), intent (in) :: phi2
code = lambda1 + atan2((cos(phi2) * sin((lambda1 - lambda2))), (cos(phi1) + (cos(phi2) * cos((lambda1 - lambda2)))))
end function
public static double code(double lambda1, double lambda2, double phi1, double phi2) {
return lambda1 + Math.atan2((Math.cos(phi2) * Math.sin((lambda1 - lambda2))), (Math.cos(phi1) + (Math.cos(phi2) * Math.cos((lambda1 - lambda2)))));
}
def code(lambda1, lambda2, phi1, phi2): return lambda1 + math.atan2((math.cos(phi2) * math.sin((lambda1 - lambda2))), (math.cos(phi1) + (math.cos(phi2) * math.cos((lambda1 - lambda2)))))
function code(lambda1, lambda2, phi1, phi2) return Float64(lambda1 + atan(Float64(cos(phi2) * sin(Float64(lambda1 - lambda2))), Float64(cos(phi1) + Float64(cos(phi2) * cos(Float64(lambda1 - lambda2)))))) end
function tmp = code(lambda1, lambda2, phi1, phi2) tmp = lambda1 + atan2((cos(phi2) * sin((lambda1 - lambda2))), (cos(phi1) + (cos(phi2) * cos((lambda1 - lambda2))))); end
code[lambda1_, lambda2_, phi1_, phi2_] := N[(lambda1 + N[ArcTan[N[(N[Cos[phi2], $MachinePrecision] * N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / N[(N[Cos[phi1], $MachinePrecision] + N[(N[Cos[phi2], $MachinePrecision] * N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)}
\end{array}
(FPCore (lambda1 lambda2 phi1 phi2) :precision binary64 (+ (atan2 (* (sin (- lambda1 lambda2)) (cos phi2)) (/ 1.0 (/ 1.0 (fma (cos (- lambda2 lambda1)) (cos phi2) (cos phi1))))) lambda1))
double code(double lambda1, double lambda2, double phi1, double phi2) {
return atan2((sin((lambda1 - lambda2)) * cos(phi2)), (1.0 / (1.0 / fma(cos((lambda2 - lambda1)), cos(phi2), cos(phi1))))) + lambda1;
}
function code(lambda1, lambda2, phi1, phi2) return Float64(atan(Float64(sin(Float64(lambda1 - lambda2)) * cos(phi2)), Float64(1.0 / Float64(1.0 / fma(cos(Float64(lambda2 - lambda1)), cos(phi2), cos(phi1))))) + lambda1) end
code[lambda1_, lambda2_, phi1_, phi2_] := N[(N[ArcTan[N[(N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision] / N[(1.0 / N[(1.0 / N[(N[Cos[N[(lambda2 - lambda1), $MachinePrecision]], $MachinePrecision] * N[Cos[phi2], $MachinePrecision] + N[Cos[phi1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] + lambda1), $MachinePrecision]
\begin{array}{l}
\\
\tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2}{\frac{1}{\frac{1}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \cos \phi_2, \cos \phi_1\right)}}} + \lambda_1
\end{array}
Initial program 99.1%
lift-+.f64N/A
flip3-+N/A
clear-numN/A
lower-/.f64N/A
clear-numN/A
flip3-+N/A
Applied rewrites99.1%
Final simplification99.1%
(FPCore (lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (* (sin (- lambda1 lambda2)) (cos phi2)))
(t_1
(+
(atan2 t_0 (+ (* (cos (- lambda1 lambda2)) (cos phi2)) (cos phi1)))
lambda1))
(t_2
(atan2 t_0 (fma (cos (- lambda2 lambda1)) (cos phi2) (cos phi1)))))
(if (<= t_1 -40.0)
(+ (atan2 (sin lambda1) (+ (cos lambda1) (cos phi1))) lambda1)
(if (<= t_1 -0.02)
t_2
(if (<= t_1 0.075)
(+ (atan2 t_0 (fma (cos lambda1) (cos phi2) (cos phi1))) lambda1)
(if (<= t_1 3.14)
t_2
(+ (atan2 t_0 (+ (cos phi1) (cos phi2))) lambda1)))))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = sin((lambda1 - lambda2)) * cos(phi2);
double t_1 = atan2(t_0, ((cos((lambda1 - lambda2)) * cos(phi2)) + cos(phi1))) + lambda1;
double t_2 = atan2(t_0, fma(cos((lambda2 - lambda1)), cos(phi2), cos(phi1)));
double tmp;
if (t_1 <= -40.0) {
tmp = atan2(sin(lambda1), (cos(lambda1) + cos(phi1))) + lambda1;
} else if (t_1 <= -0.02) {
tmp = t_2;
} else if (t_1 <= 0.075) {
tmp = atan2(t_0, fma(cos(lambda1), cos(phi2), cos(phi1))) + lambda1;
} else if (t_1 <= 3.14) {
tmp = t_2;
} else {
tmp = atan2(t_0, (cos(phi1) + cos(phi2))) + lambda1;
}
return tmp;
}
function code(lambda1, lambda2, phi1, phi2) t_0 = Float64(sin(Float64(lambda1 - lambda2)) * cos(phi2)) t_1 = Float64(atan(t_0, Float64(Float64(cos(Float64(lambda1 - lambda2)) * cos(phi2)) + cos(phi1))) + lambda1) t_2 = atan(t_0, fma(cos(Float64(lambda2 - lambda1)), cos(phi2), cos(phi1))) tmp = 0.0 if (t_1 <= -40.0) tmp = Float64(atan(sin(lambda1), Float64(cos(lambda1) + cos(phi1))) + lambda1); elseif (t_1 <= -0.02) tmp = t_2; elseif (t_1 <= 0.075) tmp = Float64(atan(t_0, fma(cos(lambda1), cos(phi2), cos(phi1))) + lambda1); elseif (t_1 <= 3.14) tmp = t_2; else tmp = Float64(atan(t_0, Float64(cos(phi1) + cos(phi2))) + lambda1); end return tmp end
code[lambda1_, lambda2_, phi1_, phi2_] := Block[{t$95$0 = N[(N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[ArcTan[t$95$0 / N[(N[(N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision] + N[Cos[phi1], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] + lambda1), $MachinePrecision]}, Block[{t$95$2 = N[ArcTan[t$95$0 / N[(N[Cos[N[(lambda2 - lambda1), $MachinePrecision]], $MachinePrecision] * N[Cos[phi2], $MachinePrecision] + N[Cos[phi1], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]}, If[LessEqual[t$95$1, -40.0], N[(N[ArcTan[N[Sin[lambda1], $MachinePrecision] / N[(N[Cos[lambda1], $MachinePrecision] + N[Cos[phi1], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] + lambda1), $MachinePrecision], If[LessEqual[t$95$1, -0.02], t$95$2, If[LessEqual[t$95$1, 0.075], N[(N[ArcTan[t$95$0 / N[(N[Cos[lambda1], $MachinePrecision] * N[Cos[phi2], $MachinePrecision] + N[Cos[phi1], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] + lambda1), $MachinePrecision], If[LessEqual[t$95$1, 3.14], t$95$2, N[(N[ArcTan[t$95$0 / N[(N[Cos[phi1], $MachinePrecision] + N[Cos[phi2], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] + lambda1), $MachinePrecision]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sin \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2\\
t_1 := \tan^{-1}_* \frac{t\_0}{\cos \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2 + \cos \phi_1} + \lambda_1\\
t_2 := \tan^{-1}_* \frac{t\_0}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \cos \phi_2, \cos \phi_1\right)}\\
\mathbf{if}\;t\_1 \leq -40:\\
\;\;\;\;\tan^{-1}_* \frac{\sin \lambda_1}{\cos \lambda_1 + \cos \phi_1} + \lambda_1\\
\mathbf{elif}\;t\_1 \leq -0.02:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;t\_1 \leq 0.075:\\
\;\;\;\;\tan^{-1}_* \frac{t\_0}{\mathsf{fma}\left(\cos \lambda_1, \cos \phi_2, \cos \phi_1\right)} + \lambda_1\\
\mathbf{elif}\;t\_1 \leq 3.14:\\
\;\;\;\;t\_2\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1}_* \frac{t\_0}{\cos \phi_1 + \cos \phi_2} + \lambda_1\\
\end{array}
\end{array}
if (+.f64 lambda1 (atan2.f64 (*.f64 (cos.f64 phi2) (sin.f64 (-.f64 lambda1 lambda2))) (+.f64 (cos.f64 phi1) (*.f64 (cos.f64 phi2) (cos.f64 (-.f64 lambda1 lambda2)))))) < -40Initial program 100.0%
Taylor expanded in phi2 around 0
lower-sin.f64N/A
lower--.f64100.0
Applied rewrites100.0%
Taylor expanded in phi2 around 0
+-commutativeN/A
sub-negN/A
remove-double-negN/A
distribute-neg-inN/A
+-commutativeN/A
neg-mul-1N/A
cos-negN/A
neg-mul-1N/A
sub-negN/A
lower-+.f64N/A
sub-negN/A
remove-double-negN/A
distribute-neg-inN/A
+-commutativeN/A
neg-mul-1N/A
cos-negN/A
lower-cos.f64N/A
neg-mul-1N/A
sub-negN/A
lower--.f64N/A
lower-cos.f64100.0
Applied rewrites100.0%
Taylor expanded in lambda2 around 0
Applied rewrites100.0%
Taylor expanded in lambda2 around 0
Applied rewrites100.0%
if -40 < (+.f64 lambda1 (atan2.f64 (*.f64 (cos.f64 phi2) (sin.f64 (-.f64 lambda1 lambda2))) (+.f64 (cos.f64 phi1) (*.f64 (cos.f64 phi2) (cos.f64 (-.f64 lambda1 lambda2)))))) < -0.0200000000000000004 or 0.0749999999999999972 < (+.f64 lambda1 (atan2.f64 (*.f64 (cos.f64 phi2) (sin.f64 (-.f64 lambda1 lambda2))) (+.f64 (cos.f64 phi1) (*.f64 (cos.f64 phi2) (cos.f64 (-.f64 lambda1 lambda2)))))) < 3.14000000000000012Initial program 96.8%
Taylor expanded in lambda1 around 0
lower-atan2.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-sin.f64N/A
lower--.f64N/A
lower-cos.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
sub-negN/A
remove-double-negN/A
mul-1-negN/A
distribute-neg-inN/A
+-commutativeN/A
cos-negN/A
lower-cos.f64N/A
mul-1-negN/A
unsub-negN/A
lower--.f64N/A
lower-cos.f64N/A
lower-cos.f6496.6
Applied rewrites96.6%
if -0.0200000000000000004 < (+.f64 lambda1 (atan2.f64 (*.f64 (cos.f64 phi2) (sin.f64 (-.f64 lambda1 lambda2))) (+.f64 (cos.f64 phi1) (*.f64 (cos.f64 phi2) (cos.f64 (-.f64 lambda1 lambda2)))))) < 0.0749999999999999972Initial program 99.3%
Taylor expanded in lambda2 around 0
+-commutativeN/A
lower-fma.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
lower-cos.f6497.9
Applied rewrites97.9%
if 3.14000000000000012 < (+.f64 lambda1 (atan2.f64 (*.f64 (cos.f64 phi2) (sin.f64 (-.f64 lambda1 lambda2))) (+.f64 (cos.f64 phi1) (*.f64 (cos.f64 phi2) (cos.f64 (-.f64 lambda1 lambda2)))))) Initial program 100.0%
Taylor expanded in lambda2 around 0
+-commutativeN/A
lower-fma.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
lower-cos.f64100.0
Applied rewrites100.0%
Taylor expanded in lambda1 around 0
Applied rewrites100.0%
Final simplification98.7%
(FPCore (lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (* (sin (- lambda1 lambda2)) (cos phi2)))
(t_1
(+
(atan2 t_0 (+ (* (cos (- lambda1 lambda2)) (cos phi2)) (cos phi1)))
lambda1))
(t_2
(atan2 t_0 (fma (cos (- lambda2 lambda1)) (cos phi2) (cos phi1)))))
(if (<= t_1 -40.0)
(+ (atan2 (sin lambda1) (+ (cos lambda1) (cos phi1))) lambda1)
(if (<= t_1 -0.02)
t_2
(if (<= t_1 2e-57)
(+
(atan2
t_0
(fma (fma (* -0.5 lambda1) lambda1 1.0) (cos phi2) (cos phi1)))
lambda1)
(if (<= t_1 3.14)
t_2
(+ (atan2 t_0 (+ (cos phi1) (cos phi2))) lambda1)))))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = sin((lambda1 - lambda2)) * cos(phi2);
double t_1 = atan2(t_0, ((cos((lambda1 - lambda2)) * cos(phi2)) + cos(phi1))) + lambda1;
double t_2 = atan2(t_0, fma(cos((lambda2 - lambda1)), cos(phi2), cos(phi1)));
double tmp;
if (t_1 <= -40.0) {
tmp = atan2(sin(lambda1), (cos(lambda1) + cos(phi1))) + lambda1;
} else if (t_1 <= -0.02) {
tmp = t_2;
} else if (t_1 <= 2e-57) {
tmp = atan2(t_0, fma(fma((-0.5 * lambda1), lambda1, 1.0), cos(phi2), cos(phi1))) + lambda1;
} else if (t_1 <= 3.14) {
tmp = t_2;
} else {
tmp = atan2(t_0, (cos(phi1) + cos(phi2))) + lambda1;
}
return tmp;
}
function code(lambda1, lambda2, phi1, phi2) t_0 = Float64(sin(Float64(lambda1 - lambda2)) * cos(phi2)) t_1 = Float64(atan(t_0, Float64(Float64(cos(Float64(lambda1 - lambda2)) * cos(phi2)) + cos(phi1))) + lambda1) t_2 = atan(t_0, fma(cos(Float64(lambda2 - lambda1)), cos(phi2), cos(phi1))) tmp = 0.0 if (t_1 <= -40.0) tmp = Float64(atan(sin(lambda1), Float64(cos(lambda1) + cos(phi1))) + lambda1); elseif (t_1 <= -0.02) tmp = t_2; elseif (t_1 <= 2e-57) tmp = Float64(atan(t_0, fma(fma(Float64(-0.5 * lambda1), lambda1, 1.0), cos(phi2), cos(phi1))) + lambda1); elseif (t_1 <= 3.14) tmp = t_2; else tmp = Float64(atan(t_0, Float64(cos(phi1) + cos(phi2))) + lambda1); end return tmp end
code[lambda1_, lambda2_, phi1_, phi2_] := Block[{t$95$0 = N[(N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[ArcTan[t$95$0 / N[(N[(N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision] + N[Cos[phi1], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] + lambda1), $MachinePrecision]}, Block[{t$95$2 = N[ArcTan[t$95$0 / N[(N[Cos[N[(lambda2 - lambda1), $MachinePrecision]], $MachinePrecision] * N[Cos[phi2], $MachinePrecision] + N[Cos[phi1], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]}, If[LessEqual[t$95$1, -40.0], N[(N[ArcTan[N[Sin[lambda1], $MachinePrecision] / N[(N[Cos[lambda1], $MachinePrecision] + N[Cos[phi1], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] + lambda1), $MachinePrecision], If[LessEqual[t$95$1, -0.02], t$95$2, If[LessEqual[t$95$1, 2e-57], N[(N[ArcTan[t$95$0 / N[(N[(N[(-0.5 * lambda1), $MachinePrecision] * lambda1 + 1.0), $MachinePrecision] * N[Cos[phi2], $MachinePrecision] + N[Cos[phi1], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] + lambda1), $MachinePrecision], If[LessEqual[t$95$1, 3.14], t$95$2, N[(N[ArcTan[t$95$0 / N[(N[Cos[phi1], $MachinePrecision] + N[Cos[phi2], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] + lambda1), $MachinePrecision]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sin \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2\\
t_1 := \tan^{-1}_* \frac{t\_0}{\cos \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2 + \cos \phi_1} + \lambda_1\\
t_2 := \tan^{-1}_* \frac{t\_0}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \cos \phi_2, \cos \phi_1\right)}\\
\mathbf{if}\;t\_1 \leq -40:\\
\;\;\;\;\tan^{-1}_* \frac{\sin \lambda_1}{\cos \lambda_1 + \cos \phi_1} + \lambda_1\\
\mathbf{elif}\;t\_1 \leq -0.02:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;t\_1 \leq 2 \cdot 10^{-57}:\\
\;\;\;\;\tan^{-1}_* \frac{t\_0}{\mathsf{fma}\left(\mathsf{fma}\left(-0.5 \cdot \lambda_1, \lambda_1, 1\right), \cos \phi_2, \cos \phi_1\right)} + \lambda_1\\
\mathbf{elif}\;t\_1 \leq 3.14:\\
\;\;\;\;t\_2\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1}_* \frac{t\_0}{\cos \phi_1 + \cos \phi_2} + \lambda_1\\
\end{array}
\end{array}
if (+.f64 lambda1 (atan2.f64 (*.f64 (cos.f64 phi2) (sin.f64 (-.f64 lambda1 lambda2))) (+.f64 (cos.f64 phi1) (*.f64 (cos.f64 phi2) (cos.f64 (-.f64 lambda1 lambda2)))))) < -40Initial program 100.0%
Taylor expanded in phi2 around 0
lower-sin.f64N/A
lower--.f64100.0
Applied rewrites100.0%
Taylor expanded in phi2 around 0
+-commutativeN/A
sub-negN/A
remove-double-negN/A
distribute-neg-inN/A
+-commutativeN/A
neg-mul-1N/A
cos-negN/A
neg-mul-1N/A
sub-negN/A
lower-+.f64N/A
sub-negN/A
remove-double-negN/A
distribute-neg-inN/A
+-commutativeN/A
neg-mul-1N/A
cos-negN/A
lower-cos.f64N/A
neg-mul-1N/A
sub-negN/A
lower--.f64N/A
lower-cos.f64100.0
Applied rewrites100.0%
Taylor expanded in lambda2 around 0
Applied rewrites100.0%
Taylor expanded in lambda2 around 0
Applied rewrites100.0%
if -40 < (+.f64 lambda1 (atan2.f64 (*.f64 (cos.f64 phi2) (sin.f64 (-.f64 lambda1 lambda2))) (+.f64 (cos.f64 phi1) (*.f64 (cos.f64 phi2) (cos.f64 (-.f64 lambda1 lambda2)))))) < -0.0200000000000000004 or 1.99999999999999991e-57 < (+.f64 lambda1 (atan2.f64 (*.f64 (cos.f64 phi2) (sin.f64 (-.f64 lambda1 lambda2))) (+.f64 (cos.f64 phi1) (*.f64 (cos.f64 phi2) (cos.f64 (-.f64 lambda1 lambda2)))))) < 3.14000000000000012Initial program 97.0%
Taylor expanded in lambda1 around 0
lower-atan2.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-sin.f64N/A
lower--.f64N/A
lower-cos.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
sub-negN/A
remove-double-negN/A
mul-1-negN/A
distribute-neg-inN/A
+-commutativeN/A
cos-negN/A
lower-cos.f64N/A
mul-1-negN/A
unsub-negN/A
lower--.f64N/A
lower-cos.f64N/A
lower-cos.f6495.4
Applied rewrites95.4%
if -0.0200000000000000004 < (+.f64 lambda1 (atan2.f64 (*.f64 (cos.f64 phi2) (sin.f64 (-.f64 lambda1 lambda2))) (+.f64 (cos.f64 phi1) (*.f64 (cos.f64 phi2) (cos.f64 (-.f64 lambda1 lambda2)))))) < 1.99999999999999991e-57Initial program 99.4%
Taylor expanded in lambda2 around 0
+-commutativeN/A
lower-fma.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
lower-cos.f6499.4
Applied rewrites99.4%
Taylor expanded in lambda1 around 0
Applied rewrites98.9%
if 3.14000000000000012 < (+.f64 lambda1 (atan2.f64 (*.f64 (cos.f64 phi2) (sin.f64 (-.f64 lambda1 lambda2))) (+.f64 (cos.f64 phi1) (*.f64 (cos.f64 phi2) (cos.f64 (-.f64 lambda1 lambda2)))))) Initial program 100.0%
Taylor expanded in lambda2 around 0
+-commutativeN/A
lower-fma.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
lower-cos.f64100.0
Applied rewrites100.0%
Taylor expanded in lambda1 around 0
Applied rewrites100.0%
Final simplification98.6%
(FPCore (lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (* (sin (- lambda1 lambda2)) (cos phi2))))
(if (<= (cos phi1) 0.996)
(+
(atan2
t_0
(fma (fma (* -0.5 lambda1) lambda1 1.0) (cos phi2) (cos phi1)))
lambda1)
(+
(atan2
t_0
(fma
(*
(fma
(fma (* phi1 phi1) -0.001388888888888889 0.041666666666666664)
(* phi1 phi1)
-0.5)
phi1)
phi1
(fma (cos (- lambda2 lambda1)) (cos phi2) 1.0)))
lambda1))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = sin((lambda1 - lambda2)) * cos(phi2);
double tmp;
if (cos(phi1) <= 0.996) {
tmp = atan2(t_0, fma(fma((-0.5 * lambda1), lambda1, 1.0), cos(phi2), cos(phi1))) + lambda1;
} else {
tmp = atan2(t_0, fma((fma(fma((phi1 * phi1), -0.001388888888888889, 0.041666666666666664), (phi1 * phi1), -0.5) * phi1), phi1, fma(cos((lambda2 - lambda1)), cos(phi2), 1.0))) + lambda1;
}
return tmp;
}
function code(lambda1, lambda2, phi1, phi2) t_0 = Float64(sin(Float64(lambda1 - lambda2)) * cos(phi2)) tmp = 0.0 if (cos(phi1) <= 0.996) tmp = Float64(atan(t_0, fma(fma(Float64(-0.5 * lambda1), lambda1, 1.0), cos(phi2), cos(phi1))) + lambda1); else tmp = Float64(atan(t_0, fma(Float64(fma(fma(Float64(phi1 * phi1), -0.001388888888888889, 0.041666666666666664), Float64(phi1 * phi1), -0.5) * phi1), phi1, fma(cos(Float64(lambda2 - lambda1)), cos(phi2), 1.0))) + lambda1); end return tmp end
code[lambda1_, lambda2_, phi1_, phi2_] := Block[{t$95$0 = N[(N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[N[Cos[phi1], $MachinePrecision], 0.996], N[(N[ArcTan[t$95$0 / N[(N[(N[(-0.5 * lambda1), $MachinePrecision] * lambda1 + 1.0), $MachinePrecision] * N[Cos[phi2], $MachinePrecision] + N[Cos[phi1], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] + lambda1), $MachinePrecision], N[(N[ArcTan[t$95$0 / N[(N[(N[(N[(N[(phi1 * phi1), $MachinePrecision] * -0.001388888888888889 + 0.041666666666666664), $MachinePrecision] * N[(phi1 * phi1), $MachinePrecision] + -0.5), $MachinePrecision] * phi1), $MachinePrecision] * phi1 + N[(N[Cos[N[(lambda2 - lambda1), $MachinePrecision]], $MachinePrecision] * N[Cos[phi2], $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] + lambda1), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sin \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2\\
\mathbf{if}\;\cos \phi_1 \leq 0.996:\\
\;\;\;\;\tan^{-1}_* \frac{t\_0}{\mathsf{fma}\left(\mathsf{fma}\left(-0.5 \cdot \lambda_1, \lambda_1, 1\right), \cos \phi_2, \cos \phi_1\right)} + \lambda_1\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1}_* \frac{t\_0}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\phi_1 \cdot \phi_1, -0.001388888888888889, 0.041666666666666664\right), \phi_1 \cdot \phi_1, -0.5\right) \cdot \phi_1, \phi_1, \mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \cos \phi_2, 1\right)\right)} + \lambda_1\\
\end{array}
\end{array}
if (cos.f64 phi1) < 0.996Initial program 99.2%
Taylor expanded in lambda2 around 0
+-commutativeN/A
lower-fma.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
lower-cos.f6483.5
Applied rewrites83.5%
Taylor expanded in lambda1 around 0
Applied rewrites82.9%
if 0.996 < (cos.f64 phi1) Initial program 98.9%
Taylor expanded in phi1 around 0
associate-+r+N/A
+-commutativeN/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
lower-fma.f64N/A
Applied rewrites98.9%
Final simplification90.5%
(FPCore (lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (* (sin (- lambda1 lambda2)) (cos phi2))))
(if (<= (cos phi1) 0.996)
(+
(atan2
t_0
(fma (fma (* -0.5 lambda1) lambda1 1.0) (cos phi2) (cos phi1)))
lambda1)
(+
(atan2
t_0
(fma
(* (fma (* phi1 phi1) 0.041666666666666664 -0.5) phi1)
phi1
(fma (cos (- lambda2 lambda1)) (cos phi2) 1.0)))
lambda1))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = sin((lambda1 - lambda2)) * cos(phi2);
double tmp;
if (cos(phi1) <= 0.996) {
tmp = atan2(t_0, fma(fma((-0.5 * lambda1), lambda1, 1.0), cos(phi2), cos(phi1))) + lambda1;
} else {
tmp = atan2(t_0, fma((fma((phi1 * phi1), 0.041666666666666664, -0.5) * phi1), phi1, fma(cos((lambda2 - lambda1)), cos(phi2), 1.0))) + lambda1;
}
return tmp;
}
function code(lambda1, lambda2, phi1, phi2) t_0 = Float64(sin(Float64(lambda1 - lambda2)) * cos(phi2)) tmp = 0.0 if (cos(phi1) <= 0.996) tmp = Float64(atan(t_0, fma(fma(Float64(-0.5 * lambda1), lambda1, 1.0), cos(phi2), cos(phi1))) + lambda1); else tmp = Float64(atan(t_0, fma(Float64(fma(Float64(phi1 * phi1), 0.041666666666666664, -0.5) * phi1), phi1, fma(cos(Float64(lambda2 - lambda1)), cos(phi2), 1.0))) + lambda1); end return tmp end
code[lambda1_, lambda2_, phi1_, phi2_] := Block[{t$95$0 = N[(N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[N[Cos[phi1], $MachinePrecision], 0.996], N[(N[ArcTan[t$95$0 / N[(N[(N[(-0.5 * lambda1), $MachinePrecision] * lambda1 + 1.0), $MachinePrecision] * N[Cos[phi2], $MachinePrecision] + N[Cos[phi1], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] + lambda1), $MachinePrecision], N[(N[ArcTan[t$95$0 / N[(N[(N[(N[(phi1 * phi1), $MachinePrecision] * 0.041666666666666664 + -0.5), $MachinePrecision] * phi1), $MachinePrecision] * phi1 + N[(N[Cos[N[(lambda2 - lambda1), $MachinePrecision]], $MachinePrecision] * N[Cos[phi2], $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] + lambda1), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sin \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2\\
\mathbf{if}\;\cos \phi_1 \leq 0.996:\\
\;\;\;\;\tan^{-1}_* \frac{t\_0}{\mathsf{fma}\left(\mathsf{fma}\left(-0.5 \cdot \lambda_1, \lambda_1, 1\right), \cos \phi_2, \cos \phi_1\right)} + \lambda_1\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1}_* \frac{t\_0}{\mathsf{fma}\left(\mathsf{fma}\left(\phi_1 \cdot \phi_1, 0.041666666666666664, -0.5\right) \cdot \phi_1, \phi_1, \mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \cos \phi_2, 1\right)\right)} + \lambda_1\\
\end{array}
\end{array}
if (cos.f64 phi1) < 0.996Initial program 99.2%
Taylor expanded in lambda2 around 0
+-commutativeN/A
lower-fma.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
lower-cos.f6483.5
Applied rewrites83.5%
Taylor expanded in lambda1 around 0
Applied rewrites82.9%
if 0.996 < (cos.f64 phi1) Initial program 98.9%
Taylor expanded in phi1 around 0
associate-+r+N/A
+-commutativeN/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
lower-fma.f64N/A
lower-*.f64N/A
sub-negN/A
*-commutativeN/A
metadata-evalN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
Applied rewrites98.8%
Final simplification90.5%
(FPCore (lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (* (sin (- lambda1 lambda2)) (cos phi2))))
(if (<= (cos phi1) 0.996)
(+
(atan2
t_0
(fma (fma (* -0.5 lambda1) lambda1 1.0) (cos phi2) (cos phi1)))
lambda1)
(+
(atan2
t_0
(fma (* -0.5 phi1) phi1 (fma (cos (- lambda2 lambda1)) (cos phi2) 1.0)))
lambda1))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = sin((lambda1 - lambda2)) * cos(phi2);
double tmp;
if (cos(phi1) <= 0.996) {
tmp = atan2(t_0, fma(fma((-0.5 * lambda1), lambda1, 1.0), cos(phi2), cos(phi1))) + lambda1;
} else {
tmp = atan2(t_0, fma((-0.5 * phi1), phi1, fma(cos((lambda2 - lambda1)), cos(phi2), 1.0))) + lambda1;
}
return tmp;
}
function code(lambda1, lambda2, phi1, phi2) t_0 = Float64(sin(Float64(lambda1 - lambda2)) * cos(phi2)) tmp = 0.0 if (cos(phi1) <= 0.996) tmp = Float64(atan(t_0, fma(fma(Float64(-0.5 * lambda1), lambda1, 1.0), cos(phi2), cos(phi1))) + lambda1); else tmp = Float64(atan(t_0, fma(Float64(-0.5 * phi1), phi1, fma(cos(Float64(lambda2 - lambda1)), cos(phi2), 1.0))) + lambda1); end return tmp end
code[lambda1_, lambda2_, phi1_, phi2_] := Block[{t$95$0 = N[(N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[N[Cos[phi1], $MachinePrecision], 0.996], N[(N[ArcTan[t$95$0 / N[(N[(N[(-0.5 * lambda1), $MachinePrecision] * lambda1 + 1.0), $MachinePrecision] * N[Cos[phi2], $MachinePrecision] + N[Cos[phi1], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] + lambda1), $MachinePrecision], N[(N[ArcTan[t$95$0 / N[(N[(-0.5 * phi1), $MachinePrecision] * phi1 + N[(N[Cos[N[(lambda2 - lambda1), $MachinePrecision]], $MachinePrecision] * N[Cos[phi2], $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] + lambda1), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sin \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2\\
\mathbf{if}\;\cos \phi_1 \leq 0.996:\\
\;\;\;\;\tan^{-1}_* \frac{t\_0}{\mathsf{fma}\left(\mathsf{fma}\left(-0.5 \cdot \lambda_1, \lambda_1, 1\right), \cos \phi_2, \cos \phi_1\right)} + \lambda_1\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1}_* \frac{t\_0}{\mathsf{fma}\left(-0.5 \cdot \phi_1, \phi_1, \mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \cos \phi_2, 1\right)\right)} + \lambda_1\\
\end{array}
\end{array}
if (cos.f64 phi1) < 0.996Initial program 99.2%
Taylor expanded in lambda2 around 0
+-commutativeN/A
lower-fma.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
lower-cos.f6483.5
Applied rewrites83.5%
Taylor expanded in lambda1 around 0
Applied rewrites82.9%
if 0.996 < (cos.f64 phi1) Initial program 98.9%
Taylor expanded in phi1 around 0
+-commutativeN/A
associate-+l+N/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
+-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
sub-negN/A
remove-double-negN/A
mul-1-negN/A
distribute-neg-inN/A
+-commutativeN/A
cos-negN/A
lower-cos.f64N/A
mul-1-negN/A
unsub-negN/A
lower--.f64N/A
lower-cos.f6498.6
Applied rewrites98.6%
Final simplification90.4%
(FPCore (lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (* (sin (- lambda1 lambda2)) (cos phi2))))
(if (<= (cos phi1) 0.996)
(+
(atan2
t_0
(fma (fma (* -0.5 lambda1) lambda1 1.0) (cos phi2) (cos phi1)))
lambda1)
(+
(atan2 t_0 (+ (* (cos (- lambda2 lambda1)) (cos phi2)) 1.0))
lambda1))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = sin((lambda1 - lambda2)) * cos(phi2);
double tmp;
if (cos(phi1) <= 0.996) {
tmp = atan2(t_0, fma(fma((-0.5 * lambda1), lambda1, 1.0), cos(phi2), cos(phi1))) + lambda1;
} else {
tmp = atan2(t_0, ((cos((lambda2 - lambda1)) * cos(phi2)) + 1.0)) + lambda1;
}
return tmp;
}
function code(lambda1, lambda2, phi1, phi2) t_0 = Float64(sin(Float64(lambda1 - lambda2)) * cos(phi2)) tmp = 0.0 if (cos(phi1) <= 0.996) tmp = Float64(atan(t_0, fma(fma(Float64(-0.5 * lambda1), lambda1, 1.0), cos(phi2), cos(phi1))) + lambda1); else tmp = Float64(atan(t_0, Float64(Float64(cos(Float64(lambda2 - lambda1)) * cos(phi2)) + 1.0)) + lambda1); end return tmp end
code[lambda1_, lambda2_, phi1_, phi2_] := Block[{t$95$0 = N[(N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[N[Cos[phi1], $MachinePrecision], 0.996], N[(N[ArcTan[t$95$0 / N[(N[(N[(-0.5 * lambda1), $MachinePrecision] * lambda1 + 1.0), $MachinePrecision] * N[Cos[phi2], $MachinePrecision] + N[Cos[phi1], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] + lambda1), $MachinePrecision], N[(N[ArcTan[t$95$0 / N[(N[(N[Cos[N[(lambda2 - lambda1), $MachinePrecision]], $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]], $MachinePrecision] + lambda1), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sin \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2\\
\mathbf{if}\;\cos \phi_1 \leq 0.996:\\
\;\;\;\;\tan^{-1}_* \frac{t\_0}{\mathsf{fma}\left(\mathsf{fma}\left(-0.5 \cdot \lambda_1, \lambda_1, 1\right), \cos \phi_2, \cos \phi_1\right)} + \lambda_1\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1}_* \frac{t\_0}{\cos \left(\lambda_2 - \lambda_1\right) \cdot \cos \phi_2 + 1} + \lambda_1\\
\end{array}
\end{array}
if (cos.f64 phi1) < 0.996Initial program 99.2%
Taylor expanded in lambda2 around 0
+-commutativeN/A
lower-fma.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
lower-cos.f6483.5
Applied rewrites83.5%
Taylor expanded in lambda1 around 0
Applied rewrites82.9%
if 0.996 < (cos.f64 phi1) Initial program 98.9%
lift-+.f64N/A
flip3-+N/A
clear-numN/A
lower-/.f64N/A
clear-numN/A
flip3-+N/A
Applied rewrites98.9%
Taylor expanded in phi1 around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
sub-negN/A
remove-double-negN/A
distribute-neg-inN/A
+-commutativeN/A
neg-mul-1N/A
cos-negN/A
lower-cos.f64N/A
neg-mul-1N/A
sub-negN/A
lower--.f64N/A
lower-cos.f6498.4
Applied rewrites98.4%
Applied rewrites98.4%
Final simplification90.3%
(FPCore (lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (* (sin (- lambda1 lambda2)) (cos phi2))))
(if (<= (cos phi1) 0.996)
(+ (atan2 t_0 (+ (cos phi1) (cos phi2))) lambda1)
(+
(atan2 t_0 (+ (* (cos (- lambda2 lambda1)) (cos phi2)) 1.0))
lambda1))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = sin((lambda1 - lambda2)) * cos(phi2);
double tmp;
if (cos(phi1) <= 0.996) {
tmp = atan2(t_0, (cos(phi1) + cos(phi2))) + lambda1;
} else {
tmp = atan2(t_0, ((cos((lambda2 - lambda1)) * cos(phi2)) + 1.0)) + lambda1;
}
return tmp;
}
real(8) function code(lambda1, lambda2, phi1, phi2)
real(8), intent (in) :: lambda1
real(8), intent (in) :: lambda2
real(8), intent (in) :: phi1
real(8), intent (in) :: phi2
real(8) :: t_0
real(8) :: tmp
t_0 = sin((lambda1 - lambda2)) * cos(phi2)
if (cos(phi1) <= 0.996d0) then
tmp = atan2(t_0, (cos(phi1) + cos(phi2))) + lambda1
else
tmp = atan2(t_0, ((cos((lambda2 - lambda1)) * cos(phi2)) + 1.0d0)) + lambda1
end if
code = tmp
end function
public static double code(double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = Math.sin((lambda1 - lambda2)) * Math.cos(phi2);
double tmp;
if (Math.cos(phi1) <= 0.996) {
tmp = Math.atan2(t_0, (Math.cos(phi1) + Math.cos(phi2))) + lambda1;
} else {
tmp = Math.atan2(t_0, ((Math.cos((lambda2 - lambda1)) * Math.cos(phi2)) + 1.0)) + lambda1;
}
return tmp;
}
def code(lambda1, lambda2, phi1, phi2): t_0 = math.sin((lambda1 - lambda2)) * math.cos(phi2) tmp = 0 if math.cos(phi1) <= 0.996: tmp = math.atan2(t_0, (math.cos(phi1) + math.cos(phi2))) + lambda1 else: tmp = math.atan2(t_0, ((math.cos((lambda2 - lambda1)) * math.cos(phi2)) + 1.0)) + lambda1 return tmp
function code(lambda1, lambda2, phi1, phi2) t_0 = Float64(sin(Float64(lambda1 - lambda2)) * cos(phi2)) tmp = 0.0 if (cos(phi1) <= 0.996) tmp = Float64(atan(t_0, Float64(cos(phi1) + cos(phi2))) + lambda1); else tmp = Float64(atan(t_0, Float64(Float64(cos(Float64(lambda2 - lambda1)) * cos(phi2)) + 1.0)) + lambda1); end return tmp end
function tmp_2 = code(lambda1, lambda2, phi1, phi2) t_0 = sin((lambda1 - lambda2)) * cos(phi2); tmp = 0.0; if (cos(phi1) <= 0.996) tmp = atan2(t_0, (cos(phi1) + cos(phi2))) + lambda1; else tmp = atan2(t_0, ((cos((lambda2 - lambda1)) * cos(phi2)) + 1.0)) + lambda1; end tmp_2 = tmp; end
code[lambda1_, lambda2_, phi1_, phi2_] := Block[{t$95$0 = N[(N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[N[Cos[phi1], $MachinePrecision], 0.996], N[(N[ArcTan[t$95$0 / N[(N[Cos[phi1], $MachinePrecision] + N[Cos[phi2], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] + lambda1), $MachinePrecision], N[(N[ArcTan[t$95$0 / N[(N[(N[Cos[N[(lambda2 - lambda1), $MachinePrecision]], $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]], $MachinePrecision] + lambda1), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sin \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2\\
\mathbf{if}\;\cos \phi_1 \leq 0.996:\\
\;\;\;\;\tan^{-1}_* \frac{t\_0}{\cos \phi_1 + \cos \phi_2} + \lambda_1\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1}_* \frac{t\_0}{\cos \left(\lambda_2 - \lambda_1\right) \cdot \cos \phi_2 + 1} + \lambda_1\\
\end{array}
\end{array}
if (cos.f64 phi1) < 0.996Initial program 99.2%
Taylor expanded in lambda2 around 0
+-commutativeN/A
lower-fma.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
lower-cos.f6483.5
Applied rewrites83.5%
Taylor expanded in lambda1 around 0
Applied rewrites82.6%
if 0.996 < (cos.f64 phi1) Initial program 98.9%
lift-+.f64N/A
flip3-+N/A
clear-numN/A
lower-/.f64N/A
clear-numN/A
flip3-+N/A
Applied rewrites98.9%
Taylor expanded in phi1 around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
sub-negN/A
remove-double-negN/A
distribute-neg-inN/A
+-commutativeN/A
neg-mul-1N/A
cos-negN/A
lower-cos.f64N/A
neg-mul-1N/A
sub-negN/A
lower--.f64N/A
lower-cos.f6498.4
Applied rewrites98.4%
Applied rewrites98.4%
Final simplification90.1%
(FPCore (lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (* (sin (- lambda1 lambda2)) (cos phi2))))
(if (<= (cos phi1) 0.996)
(+ (atan2 t_0 (+ (cos phi1) (cos phi2))) lambda1)
(+ (atan2 t_0 (fma (cos (- lambda2 lambda1)) (cos phi2) 1.0)) lambda1))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = sin((lambda1 - lambda2)) * cos(phi2);
double tmp;
if (cos(phi1) <= 0.996) {
tmp = atan2(t_0, (cos(phi1) + cos(phi2))) + lambda1;
} else {
tmp = atan2(t_0, fma(cos((lambda2 - lambda1)), cos(phi2), 1.0)) + lambda1;
}
return tmp;
}
function code(lambda1, lambda2, phi1, phi2) t_0 = Float64(sin(Float64(lambda1 - lambda2)) * cos(phi2)) tmp = 0.0 if (cos(phi1) <= 0.996) tmp = Float64(atan(t_0, Float64(cos(phi1) + cos(phi2))) + lambda1); else tmp = Float64(atan(t_0, fma(cos(Float64(lambda2 - lambda1)), cos(phi2), 1.0)) + lambda1); end return tmp end
code[lambda1_, lambda2_, phi1_, phi2_] := Block[{t$95$0 = N[(N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[N[Cos[phi1], $MachinePrecision], 0.996], N[(N[ArcTan[t$95$0 / N[(N[Cos[phi1], $MachinePrecision] + N[Cos[phi2], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] + lambda1), $MachinePrecision], N[(N[ArcTan[t$95$0 / N[(N[Cos[N[(lambda2 - lambda1), $MachinePrecision]], $MachinePrecision] * N[Cos[phi2], $MachinePrecision] + 1.0), $MachinePrecision]], $MachinePrecision] + lambda1), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sin \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2\\
\mathbf{if}\;\cos \phi_1 \leq 0.996:\\
\;\;\;\;\tan^{-1}_* \frac{t\_0}{\cos \phi_1 + \cos \phi_2} + \lambda_1\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1}_* \frac{t\_0}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \cos \phi_2, 1\right)} + \lambda_1\\
\end{array}
\end{array}
if (cos.f64 phi1) < 0.996Initial program 99.2%
Taylor expanded in lambda2 around 0
+-commutativeN/A
lower-fma.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
lower-cos.f6483.5
Applied rewrites83.5%
Taylor expanded in lambda1 around 0
Applied rewrites82.6%
if 0.996 < (cos.f64 phi1) Initial program 98.9%
Taylor expanded in phi1 around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
sub-negN/A
remove-double-negN/A
mul-1-negN/A
distribute-neg-inN/A
+-commutativeN/A
cos-negN/A
lower-cos.f64N/A
mul-1-negN/A
unsub-negN/A
lower--.f64N/A
lower-cos.f6498.4
Applied rewrites98.4%
Final simplification90.1%
(FPCore (lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (* (sin (- lambda1 lambda2)) (cos phi2))))
(if (<= (cos phi1) 0.996)
(+ (atan2 t_0 (+ (cos phi1) (cos phi2))) lambda1)
(+ (atan2 t_0 (fma (cos lambda2) (cos phi2) 1.0)) lambda1))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = sin((lambda1 - lambda2)) * cos(phi2);
double tmp;
if (cos(phi1) <= 0.996) {
tmp = atan2(t_0, (cos(phi1) + cos(phi2))) + lambda1;
} else {
tmp = atan2(t_0, fma(cos(lambda2), cos(phi2), 1.0)) + lambda1;
}
return tmp;
}
function code(lambda1, lambda2, phi1, phi2) t_0 = Float64(sin(Float64(lambda1 - lambda2)) * cos(phi2)) tmp = 0.0 if (cos(phi1) <= 0.996) tmp = Float64(atan(t_0, Float64(cos(phi1) + cos(phi2))) + lambda1); else tmp = Float64(atan(t_0, fma(cos(lambda2), cos(phi2), 1.0)) + lambda1); end return tmp end
code[lambda1_, lambda2_, phi1_, phi2_] := Block[{t$95$0 = N[(N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[N[Cos[phi1], $MachinePrecision], 0.996], N[(N[ArcTan[t$95$0 / N[(N[Cos[phi1], $MachinePrecision] + N[Cos[phi2], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] + lambda1), $MachinePrecision], N[(N[ArcTan[t$95$0 / N[(N[Cos[lambda2], $MachinePrecision] * N[Cos[phi2], $MachinePrecision] + 1.0), $MachinePrecision]], $MachinePrecision] + lambda1), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sin \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2\\
\mathbf{if}\;\cos \phi_1 \leq 0.996:\\
\;\;\;\;\tan^{-1}_* \frac{t\_0}{\cos \phi_1 + \cos \phi_2} + \lambda_1\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1}_* \frac{t\_0}{\mathsf{fma}\left(\cos \lambda_2, \cos \phi_2, 1\right)} + \lambda_1\\
\end{array}
\end{array}
if (cos.f64 phi1) < 0.996Initial program 99.2%
Taylor expanded in lambda2 around 0
+-commutativeN/A
lower-fma.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
lower-cos.f6483.5
Applied rewrites83.5%
Taylor expanded in lambda1 around 0
Applied rewrites82.6%
if 0.996 < (cos.f64 phi1) Initial program 98.9%
lift-+.f64N/A
flip3-+N/A
clear-numN/A
lower-/.f64N/A
clear-numN/A
flip3-+N/A
Applied rewrites98.9%
Taylor expanded in phi1 around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
sub-negN/A
remove-double-negN/A
distribute-neg-inN/A
+-commutativeN/A
neg-mul-1N/A
cos-negN/A
lower-cos.f64N/A
neg-mul-1N/A
sub-negN/A
lower--.f64N/A
lower-cos.f6498.4
Applied rewrites98.4%
Taylor expanded in lambda1 around 0
Applied rewrites97.9%
Final simplification89.9%
(FPCore (lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (* (sin (- lambda1 lambda2)) (cos phi2))))
(if (<= (cos phi2) 0.995)
(+ (atan2 t_0 (+ (cos phi1) (cos phi2))) lambda1)
(+ (atan2 t_0 (+ (cos phi1) (cos (- lambda2 lambda1)))) lambda1))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = sin((lambda1 - lambda2)) * cos(phi2);
double tmp;
if (cos(phi2) <= 0.995) {
tmp = atan2(t_0, (cos(phi1) + cos(phi2))) + lambda1;
} else {
tmp = atan2(t_0, (cos(phi1) + cos((lambda2 - lambda1)))) + lambda1;
}
return tmp;
}
real(8) function code(lambda1, lambda2, phi1, phi2)
real(8), intent (in) :: lambda1
real(8), intent (in) :: lambda2
real(8), intent (in) :: phi1
real(8), intent (in) :: phi2
real(8) :: t_0
real(8) :: tmp
t_0 = sin((lambda1 - lambda2)) * cos(phi2)
if (cos(phi2) <= 0.995d0) then
tmp = atan2(t_0, (cos(phi1) + cos(phi2))) + lambda1
else
tmp = atan2(t_0, (cos(phi1) + cos((lambda2 - lambda1)))) + lambda1
end if
code = tmp
end function
public static double code(double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = Math.sin((lambda1 - lambda2)) * Math.cos(phi2);
double tmp;
if (Math.cos(phi2) <= 0.995) {
tmp = Math.atan2(t_0, (Math.cos(phi1) + Math.cos(phi2))) + lambda1;
} else {
tmp = Math.atan2(t_0, (Math.cos(phi1) + Math.cos((lambda2 - lambda1)))) + lambda1;
}
return tmp;
}
def code(lambda1, lambda2, phi1, phi2): t_0 = math.sin((lambda1 - lambda2)) * math.cos(phi2) tmp = 0 if math.cos(phi2) <= 0.995: tmp = math.atan2(t_0, (math.cos(phi1) + math.cos(phi2))) + lambda1 else: tmp = math.atan2(t_0, (math.cos(phi1) + math.cos((lambda2 - lambda1)))) + lambda1 return tmp
function code(lambda1, lambda2, phi1, phi2) t_0 = Float64(sin(Float64(lambda1 - lambda2)) * cos(phi2)) tmp = 0.0 if (cos(phi2) <= 0.995) tmp = Float64(atan(t_0, Float64(cos(phi1) + cos(phi2))) + lambda1); else tmp = Float64(atan(t_0, Float64(cos(phi1) + cos(Float64(lambda2 - lambda1)))) + lambda1); end return tmp end
function tmp_2 = code(lambda1, lambda2, phi1, phi2) t_0 = sin((lambda1 - lambda2)) * cos(phi2); tmp = 0.0; if (cos(phi2) <= 0.995) tmp = atan2(t_0, (cos(phi1) + cos(phi2))) + lambda1; else tmp = atan2(t_0, (cos(phi1) + cos((lambda2 - lambda1)))) + lambda1; end tmp_2 = tmp; end
code[lambda1_, lambda2_, phi1_, phi2_] := Block[{t$95$0 = N[(N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[N[Cos[phi2], $MachinePrecision], 0.995], N[(N[ArcTan[t$95$0 / N[(N[Cos[phi1], $MachinePrecision] + N[Cos[phi2], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] + lambda1), $MachinePrecision], N[(N[ArcTan[t$95$0 / N[(N[Cos[phi1], $MachinePrecision] + N[Cos[N[(lambda2 - lambda1), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] + lambda1), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sin \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2\\
\mathbf{if}\;\cos \phi_2 \leq 0.995:\\
\;\;\;\;\tan^{-1}_* \frac{t\_0}{\cos \phi_1 + \cos \phi_2} + \lambda_1\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1}_* \frac{t\_0}{\cos \phi_1 + \cos \left(\lambda_2 - \lambda_1\right)} + \lambda_1\\
\end{array}
\end{array}
if (cos.f64 phi2) < 0.994999999999999996Initial program 99.2%
Taylor expanded in lambda2 around 0
+-commutativeN/A
lower-fma.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
lower-cos.f6481.4
Applied rewrites81.4%
Taylor expanded in lambda1 around 0
Applied rewrites80.4%
if 0.994999999999999996 < (cos.f64 phi2) Initial program 98.9%
Taylor expanded in phi2 around 0
+-commutativeN/A
lower-+.f64N/A
sub-negN/A
remove-double-negN/A
mul-1-negN/A
distribute-neg-inN/A
+-commutativeN/A
cos-negN/A
lower-cos.f64N/A
mul-1-negN/A
unsub-negN/A
lower--.f64N/A
lower-cos.f6498.0
Applied rewrites98.0%
Final simplification89.1%
(FPCore (lambda1 lambda2 phi1 phi2) :precision binary64 (+ (atan2 (* (sin (- lambda1 lambda2)) (cos phi2)) (+ (* (cos (- lambda1 lambda2)) (cos phi2)) (cos phi1))) lambda1))
double code(double lambda1, double lambda2, double phi1, double phi2) {
return atan2((sin((lambda1 - lambda2)) * cos(phi2)), ((cos((lambda1 - lambda2)) * cos(phi2)) + cos(phi1))) + lambda1;
}
real(8) function code(lambda1, lambda2, phi1, phi2)
real(8), intent (in) :: lambda1
real(8), intent (in) :: lambda2
real(8), intent (in) :: phi1
real(8), intent (in) :: phi2
code = atan2((sin((lambda1 - lambda2)) * cos(phi2)), ((cos((lambda1 - lambda2)) * cos(phi2)) + cos(phi1))) + lambda1
end function
public static double code(double lambda1, double lambda2, double phi1, double phi2) {
return Math.atan2((Math.sin((lambda1 - lambda2)) * Math.cos(phi2)), ((Math.cos((lambda1 - lambda2)) * Math.cos(phi2)) + Math.cos(phi1))) + lambda1;
}
def code(lambda1, lambda2, phi1, phi2): return math.atan2((math.sin((lambda1 - lambda2)) * math.cos(phi2)), ((math.cos((lambda1 - lambda2)) * math.cos(phi2)) + math.cos(phi1))) + lambda1
function code(lambda1, lambda2, phi1, phi2) return Float64(atan(Float64(sin(Float64(lambda1 - lambda2)) * cos(phi2)), Float64(Float64(cos(Float64(lambda1 - lambda2)) * cos(phi2)) + cos(phi1))) + lambda1) end
function tmp = code(lambda1, lambda2, phi1, phi2) tmp = atan2((sin((lambda1 - lambda2)) * cos(phi2)), ((cos((lambda1 - lambda2)) * cos(phi2)) + cos(phi1))) + lambda1; end
code[lambda1_, lambda2_, phi1_, phi2_] := N[(N[ArcTan[N[(N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision] / N[(N[(N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision] + N[Cos[phi1], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] + lambda1), $MachinePrecision]
\begin{array}{l}
\\
\tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2}{\cos \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2 + \cos \phi_1} + \lambda_1
\end{array}
Initial program 99.1%
Final simplification99.1%
(FPCore (lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (sin (- lambda1 lambda2))))
(if (<= (cos phi2) 0.995)
(+ (atan2 (* t_0 (cos phi2)) (+ (cos phi1) (cos phi2))) lambda1)
(+ (atan2 t_0 (+ (cos (- lambda1 lambda2)) (cos phi1))) lambda1))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = sin((lambda1 - lambda2));
double tmp;
if (cos(phi2) <= 0.995) {
tmp = atan2((t_0 * cos(phi2)), (cos(phi1) + cos(phi2))) + lambda1;
} else {
tmp = atan2(t_0, (cos((lambda1 - lambda2)) + cos(phi1))) + lambda1;
}
return tmp;
}
real(8) function code(lambda1, lambda2, phi1, phi2)
real(8), intent (in) :: lambda1
real(8), intent (in) :: lambda2
real(8), intent (in) :: phi1
real(8), intent (in) :: phi2
real(8) :: t_0
real(8) :: tmp
t_0 = sin((lambda1 - lambda2))
if (cos(phi2) <= 0.995d0) then
tmp = atan2((t_0 * cos(phi2)), (cos(phi1) + cos(phi2))) + lambda1
else
tmp = atan2(t_0, (cos((lambda1 - lambda2)) + cos(phi1))) + lambda1
end if
code = tmp
end function
public static double code(double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = Math.sin((lambda1 - lambda2));
double tmp;
if (Math.cos(phi2) <= 0.995) {
tmp = Math.atan2((t_0 * Math.cos(phi2)), (Math.cos(phi1) + Math.cos(phi2))) + lambda1;
} else {
tmp = Math.atan2(t_0, (Math.cos((lambda1 - lambda2)) + Math.cos(phi1))) + lambda1;
}
return tmp;
}
def code(lambda1, lambda2, phi1, phi2): t_0 = math.sin((lambda1 - lambda2)) tmp = 0 if math.cos(phi2) <= 0.995: tmp = math.atan2((t_0 * math.cos(phi2)), (math.cos(phi1) + math.cos(phi2))) + lambda1 else: tmp = math.atan2(t_0, (math.cos((lambda1 - lambda2)) + math.cos(phi1))) + lambda1 return tmp
function code(lambda1, lambda2, phi1, phi2) t_0 = sin(Float64(lambda1 - lambda2)) tmp = 0.0 if (cos(phi2) <= 0.995) tmp = Float64(atan(Float64(t_0 * cos(phi2)), Float64(cos(phi1) + cos(phi2))) + lambda1); else tmp = Float64(atan(t_0, Float64(cos(Float64(lambda1 - lambda2)) + cos(phi1))) + lambda1); end return tmp end
function tmp_2 = code(lambda1, lambda2, phi1, phi2) t_0 = sin((lambda1 - lambda2)); tmp = 0.0; if (cos(phi2) <= 0.995) tmp = atan2((t_0 * cos(phi2)), (cos(phi1) + cos(phi2))) + lambda1; else tmp = atan2(t_0, (cos((lambda1 - lambda2)) + cos(phi1))) + lambda1; end tmp_2 = tmp; end
code[lambda1_, lambda2_, phi1_, phi2_] := Block[{t$95$0 = N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]}, If[LessEqual[N[Cos[phi2], $MachinePrecision], 0.995], N[(N[ArcTan[N[(t$95$0 * N[Cos[phi2], $MachinePrecision]), $MachinePrecision] / N[(N[Cos[phi1], $MachinePrecision] + N[Cos[phi2], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] + lambda1), $MachinePrecision], N[(N[ArcTan[t$95$0 / N[(N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision] + N[Cos[phi1], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] + lambda1), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sin \left(\lambda_1 - \lambda_2\right)\\
\mathbf{if}\;\cos \phi_2 \leq 0.995:\\
\;\;\;\;\tan^{-1}_* \frac{t\_0 \cdot \cos \phi_2}{\cos \phi_1 + \cos \phi_2} + \lambda_1\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1}_* \frac{t\_0}{\cos \left(\lambda_1 - \lambda_2\right) + \cos \phi_1} + \lambda_1\\
\end{array}
\end{array}
if (cos.f64 phi2) < 0.994999999999999996Initial program 99.2%
Taylor expanded in lambda2 around 0
+-commutativeN/A
lower-fma.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
lower-cos.f6481.4
Applied rewrites81.4%
Taylor expanded in lambda1 around 0
Applied rewrites80.4%
if 0.994999999999999996 < (cos.f64 phi2) Initial program 98.9%
Taylor expanded in phi2 around 0
lower-sin.f64N/A
lower--.f6497.9
Applied rewrites97.9%
Taylor expanded in phi2 around 0
+-commutativeN/A
sub-negN/A
remove-double-negN/A
distribute-neg-inN/A
+-commutativeN/A
neg-mul-1N/A
cos-negN/A
neg-mul-1N/A
sub-negN/A
lower-+.f64N/A
sub-negN/A
remove-double-negN/A
distribute-neg-inN/A
+-commutativeN/A
neg-mul-1N/A
cos-negN/A
lower-cos.f64N/A
neg-mul-1N/A
sub-negN/A
lower--.f64N/A
lower-cos.f6498.0
Applied rewrites98.0%
Final simplification89.1%
(FPCore (lambda1 lambda2 phi1 phi2) :precision binary64 (+ (atan2 (* (sin (- lambda1 lambda2)) (cos phi2)) (+ (* (cos lambda2) (cos phi2)) (cos phi1))) lambda1))
double code(double lambda1, double lambda2, double phi1, double phi2) {
return atan2((sin((lambda1 - lambda2)) * cos(phi2)), ((cos(lambda2) * cos(phi2)) + cos(phi1))) + lambda1;
}
real(8) function code(lambda1, lambda2, phi1, phi2)
real(8), intent (in) :: lambda1
real(8), intent (in) :: lambda2
real(8), intent (in) :: phi1
real(8), intent (in) :: phi2
code = atan2((sin((lambda1 - lambda2)) * cos(phi2)), ((cos(lambda2) * cos(phi2)) + cos(phi1))) + lambda1
end function
public static double code(double lambda1, double lambda2, double phi1, double phi2) {
return Math.atan2((Math.sin((lambda1 - lambda2)) * Math.cos(phi2)), ((Math.cos(lambda2) * Math.cos(phi2)) + Math.cos(phi1))) + lambda1;
}
def code(lambda1, lambda2, phi1, phi2): return math.atan2((math.sin((lambda1 - lambda2)) * math.cos(phi2)), ((math.cos(lambda2) * math.cos(phi2)) + math.cos(phi1))) + lambda1
function code(lambda1, lambda2, phi1, phi2) return Float64(atan(Float64(sin(Float64(lambda1 - lambda2)) * cos(phi2)), Float64(Float64(cos(lambda2) * cos(phi2)) + cos(phi1))) + lambda1) end
function tmp = code(lambda1, lambda2, phi1, phi2) tmp = atan2((sin((lambda1 - lambda2)) * cos(phi2)), ((cos(lambda2) * cos(phi2)) + cos(phi1))) + lambda1; end
code[lambda1_, lambda2_, phi1_, phi2_] := N[(N[ArcTan[N[(N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision] / N[(N[(N[Cos[lambda2], $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision] + N[Cos[phi1], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] + lambda1), $MachinePrecision]
\begin{array}{l}
\\
\tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2}{\cos \lambda_2 \cdot \cos \phi_2 + \cos \phi_1} + \lambda_1
\end{array}
Initial program 99.1%
Taylor expanded in lambda1 around 0
cos-negN/A
lower-cos.f6498.3
Applied rewrites98.3%
Final simplification98.3%
(FPCore (lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (sin (- lambda1 lambda2))))
(if (<= (cos phi2) 0.988)
(+
(atan2
(* t_0 (cos phi2))
(fma (fma (* lambda1 lambda1) -0.5 1.0) (cos phi2) 1.0))
lambda1)
(+ (atan2 t_0 (+ (cos (- lambda1 lambda2)) (cos phi1))) lambda1))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = sin((lambda1 - lambda2));
double tmp;
if (cos(phi2) <= 0.988) {
tmp = atan2((t_0 * cos(phi2)), fma(fma((lambda1 * lambda1), -0.5, 1.0), cos(phi2), 1.0)) + lambda1;
} else {
tmp = atan2(t_0, (cos((lambda1 - lambda2)) + cos(phi1))) + lambda1;
}
return tmp;
}
function code(lambda1, lambda2, phi1, phi2) t_0 = sin(Float64(lambda1 - lambda2)) tmp = 0.0 if (cos(phi2) <= 0.988) tmp = Float64(atan(Float64(t_0 * cos(phi2)), fma(fma(Float64(lambda1 * lambda1), -0.5, 1.0), cos(phi2), 1.0)) + lambda1); else tmp = Float64(atan(t_0, Float64(cos(Float64(lambda1 - lambda2)) + cos(phi1))) + lambda1); end return tmp end
code[lambda1_, lambda2_, phi1_, phi2_] := Block[{t$95$0 = N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]}, If[LessEqual[N[Cos[phi2], $MachinePrecision], 0.988], N[(N[ArcTan[N[(t$95$0 * N[Cos[phi2], $MachinePrecision]), $MachinePrecision] / N[(N[(N[(lambda1 * lambda1), $MachinePrecision] * -0.5 + 1.0), $MachinePrecision] * N[Cos[phi2], $MachinePrecision] + 1.0), $MachinePrecision]], $MachinePrecision] + lambda1), $MachinePrecision], N[(N[ArcTan[t$95$0 / N[(N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision] + N[Cos[phi1], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] + lambda1), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sin \left(\lambda_1 - \lambda_2\right)\\
\mathbf{if}\;\cos \phi_2 \leq 0.988:\\
\;\;\;\;\tan^{-1}_* \frac{t\_0 \cdot \cos \phi_2}{\mathsf{fma}\left(\mathsf{fma}\left(\lambda_1 \cdot \lambda_1, -0.5, 1\right), \cos \phi_2, 1\right)} + \lambda_1\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1}_* \frac{t\_0}{\cos \left(\lambda_1 - \lambda_2\right) + \cos \phi_1} + \lambda_1\\
\end{array}
\end{array}
if (cos.f64 phi2) < 0.98799999999999999Initial program 99.2%
Taylor expanded in lambda2 around 0
+-commutativeN/A
lower-fma.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
lower-cos.f6480.9
Applied rewrites80.9%
Taylor expanded in lambda1 around 0
Applied rewrites80.3%
Taylor expanded in phi1 around 0
Applied rewrites66.0%
if 0.98799999999999999 < (cos.f64 phi2) Initial program 99.0%
Taylor expanded in phi2 around 0
lower-sin.f64N/A
lower--.f6497.4
Applied rewrites97.4%
Taylor expanded in phi2 around 0
+-commutativeN/A
sub-negN/A
remove-double-negN/A
distribute-neg-inN/A
+-commutativeN/A
neg-mul-1N/A
cos-negN/A
neg-mul-1N/A
sub-negN/A
lower-+.f64N/A
sub-negN/A
remove-double-negN/A
distribute-neg-inN/A
+-commutativeN/A
neg-mul-1N/A
cos-negN/A
lower-cos.f64N/A
neg-mul-1N/A
sub-negN/A
lower--.f64N/A
lower-cos.f6497.5
Applied rewrites97.5%
Final simplification81.9%
(FPCore (lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (cos (- lambda1 lambda2))))
(if (<= (cos phi2) 0.4)
(+
(atan2 (sin lambda1) (fma (fma (* phi2 phi2) -0.5 1.0) t_0 (cos phi1)))
lambda1)
(+ (atan2 (sin (- lambda1 lambda2)) (+ t_0 (cos phi1))) lambda1))))
double code(double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = cos((lambda1 - lambda2));
double tmp;
if (cos(phi2) <= 0.4) {
tmp = atan2(sin(lambda1), fma(fma((phi2 * phi2), -0.5, 1.0), t_0, cos(phi1))) + lambda1;
} else {
tmp = atan2(sin((lambda1 - lambda2)), (t_0 + cos(phi1))) + lambda1;
}
return tmp;
}
function code(lambda1, lambda2, phi1, phi2) t_0 = cos(Float64(lambda1 - lambda2)) tmp = 0.0 if (cos(phi2) <= 0.4) tmp = Float64(atan(sin(lambda1), fma(fma(Float64(phi2 * phi2), -0.5, 1.0), t_0, cos(phi1))) + lambda1); else tmp = Float64(atan(sin(Float64(lambda1 - lambda2)), Float64(t_0 + cos(phi1))) + lambda1); end return tmp end
code[lambda1_, lambda2_, phi1_, phi2_] := Block[{t$95$0 = N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]}, If[LessEqual[N[Cos[phi2], $MachinePrecision], 0.4], N[(N[ArcTan[N[Sin[lambda1], $MachinePrecision] / N[(N[(N[(phi2 * phi2), $MachinePrecision] * -0.5 + 1.0), $MachinePrecision] * t$95$0 + N[Cos[phi1], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] + lambda1), $MachinePrecision], N[(N[ArcTan[N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision] / N[(t$95$0 + N[Cos[phi1], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] + lambda1), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \cos \left(\lambda_1 - \lambda_2\right)\\
\mathbf{if}\;\cos \phi_2 \leq 0.4:\\
\;\;\;\;\tan^{-1}_* \frac{\sin \lambda_1}{\mathsf{fma}\left(\mathsf{fma}\left(\phi_2 \cdot \phi_2, -0.5, 1\right), t\_0, \cos \phi_1\right)} + \lambda_1\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right)}{t\_0 + \cos \phi_1} + \lambda_1\\
\end{array}
\end{array}
if (cos.f64 phi2) < 0.40000000000000002Initial program 99.2%
Taylor expanded in phi2 around 0
lower-sin.f64N/A
lower--.f6456.8
Applied rewrites56.8%
Taylor expanded in phi2 around 0
+-commutativeN/A
sub-negN/A
remove-double-negN/A
distribute-neg-inN/A
+-commutativeN/A
neg-mul-1N/A
cos-negN/A
neg-mul-1N/A
sub-negN/A
lower-+.f64N/A
sub-negN/A
remove-double-negN/A
distribute-neg-inN/A
+-commutativeN/A
neg-mul-1N/A
cos-negN/A
lower-cos.f64N/A
neg-mul-1N/A
sub-negN/A
lower--.f64N/A
lower-cos.f6453.7
Applied rewrites53.7%
Taylor expanded in lambda2 around 0
Applied rewrites53.5%
Taylor expanded in phi2 around 0
+-commutativeN/A
associate-+l+N/A
sub-negN/A
neg-mul-1N/A
cos-negN/A
neg-mul-1N/A
+-commutativeN/A
distribute-neg-inN/A
remove-double-negN/A
sub-negN/A
Applied rewrites59.2%
if 0.40000000000000002 < (cos.f64 phi2) Initial program 99.0%
Taylor expanded in phi2 around 0
lower-sin.f64N/A
lower--.f6489.0
Applied rewrites89.0%
Taylor expanded in phi2 around 0
+-commutativeN/A
sub-negN/A
remove-double-negN/A
distribute-neg-inN/A
+-commutativeN/A
neg-mul-1N/A
cos-negN/A
neg-mul-1N/A
sub-negN/A
lower-+.f64N/A
sub-negN/A
remove-double-negN/A
distribute-neg-inN/A
+-commutativeN/A
neg-mul-1N/A
cos-negN/A
lower-cos.f64N/A
neg-mul-1N/A
sub-negN/A
lower--.f64N/A
lower-cos.f6489.2
Applied rewrites89.2%
Final simplification78.8%
(FPCore (lambda1 lambda2 phi1 phi2) :precision binary64 (+ (atan2 (sin (- lambda1 lambda2)) (+ (cos (- lambda1 lambda2)) (cos phi1))) lambda1))
double code(double lambda1, double lambda2, double phi1, double phi2) {
return atan2(sin((lambda1 - lambda2)), (cos((lambda1 - lambda2)) + cos(phi1))) + lambda1;
}
real(8) function code(lambda1, lambda2, phi1, phi2)
real(8), intent (in) :: lambda1
real(8), intent (in) :: lambda2
real(8), intent (in) :: phi1
real(8), intent (in) :: phi2
code = atan2(sin((lambda1 - lambda2)), (cos((lambda1 - lambda2)) + cos(phi1))) + lambda1
end function
public static double code(double lambda1, double lambda2, double phi1, double phi2) {
return Math.atan2(Math.sin((lambda1 - lambda2)), (Math.cos((lambda1 - lambda2)) + Math.cos(phi1))) + lambda1;
}
def code(lambda1, lambda2, phi1, phi2): return math.atan2(math.sin((lambda1 - lambda2)), (math.cos((lambda1 - lambda2)) + math.cos(phi1))) + lambda1
function code(lambda1, lambda2, phi1, phi2) return Float64(atan(sin(Float64(lambda1 - lambda2)), Float64(cos(Float64(lambda1 - lambda2)) + cos(phi1))) + lambda1) end
function tmp = code(lambda1, lambda2, phi1, phi2) tmp = atan2(sin((lambda1 - lambda2)), (cos((lambda1 - lambda2)) + cos(phi1))) + lambda1; end
code[lambda1_, lambda2_, phi1_, phi2_] := N[(N[ArcTan[N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision] / N[(N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision] + N[Cos[phi1], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] + lambda1), $MachinePrecision]
\begin{array}{l}
\\
\tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right)}{\cos \left(\lambda_1 - \lambda_2\right) + \cos \phi_1} + \lambda_1
\end{array}
Initial program 99.1%
Taylor expanded in phi2 around 0
lower-sin.f64N/A
lower--.f6477.8
Applied rewrites77.8%
Taylor expanded in phi2 around 0
+-commutativeN/A
sub-negN/A
remove-double-negN/A
distribute-neg-inN/A
+-commutativeN/A
neg-mul-1N/A
cos-negN/A
neg-mul-1N/A
sub-negN/A
lower-+.f64N/A
sub-negN/A
remove-double-negN/A
distribute-neg-inN/A
+-commutativeN/A
neg-mul-1N/A
cos-negN/A
lower-cos.f64N/A
neg-mul-1N/A
sub-negN/A
lower--.f64N/A
lower-cos.f6476.9
Applied rewrites76.9%
Final simplification76.9%
(FPCore (lambda1 lambda2 phi1 phi2) :precision binary64 (+ (atan2 (sin (- lambda1 lambda2)) (+ (cos lambda2) (cos phi1))) lambda1))
double code(double lambda1, double lambda2, double phi1, double phi2) {
return atan2(sin((lambda1 - lambda2)), (cos(lambda2) + cos(phi1))) + lambda1;
}
real(8) function code(lambda1, lambda2, phi1, phi2)
real(8), intent (in) :: lambda1
real(8), intent (in) :: lambda2
real(8), intent (in) :: phi1
real(8), intent (in) :: phi2
code = atan2(sin((lambda1 - lambda2)), (cos(lambda2) + cos(phi1))) + lambda1
end function
public static double code(double lambda1, double lambda2, double phi1, double phi2) {
return Math.atan2(Math.sin((lambda1 - lambda2)), (Math.cos(lambda2) + Math.cos(phi1))) + lambda1;
}
def code(lambda1, lambda2, phi1, phi2): return math.atan2(math.sin((lambda1 - lambda2)), (math.cos(lambda2) + math.cos(phi1))) + lambda1
function code(lambda1, lambda2, phi1, phi2) return Float64(atan(sin(Float64(lambda1 - lambda2)), Float64(cos(lambda2) + cos(phi1))) + lambda1) end
function tmp = code(lambda1, lambda2, phi1, phi2) tmp = atan2(sin((lambda1 - lambda2)), (cos(lambda2) + cos(phi1))) + lambda1; end
code[lambda1_, lambda2_, phi1_, phi2_] := N[(N[ArcTan[N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision] / N[(N[Cos[lambda2], $MachinePrecision] + N[Cos[phi1], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] + lambda1), $MachinePrecision]
\begin{array}{l}
\\
\tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right)}{\cos \lambda_2 + \cos \phi_1} + \lambda_1
\end{array}
Initial program 99.1%
Taylor expanded in phi2 around 0
lower-sin.f64N/A
lower--.f6477.8
Applied rewrites77.8%
Taylor expanded in phi2 around 0
+-commutativeN/A
sub-negN/A
remove-double-negN/A
distribute-neg-inN/A
+-commutativeN/A
neg-mul-1N/A
cos-negN/A
neg-mul-1N/A
sub-negN/A
lower-+.f64N/A
sub-negN/A
remove-double-negN/A
distribute-neg-inN/A
+-commutativeN/A
neg-mul-1N/A
cos-negN/A
lower-cos.f64N/A
neg-mul-1N/A
sub-negN/A
lower--.f64N/A
lower-cos.f6476.9
Applied rewrites76.9%
Taylor expanded in lambda1 around 0
Applied rewrites76.6%
Final simplification76.6%
(FPCore (lambda1 lambda2 phi1 phi2) :precision binary64 (+ (atan2 (sin (- lambda1 lambda2)) (+ (cos lambda1) (cos phi1))) lambda1))
double code(double lambda1, double lambda2, double phi1, double phi2) {
return atan2(sin((lambda1 - lambda2)), (cos(lambda1) + cos(phi1))) + lambda1;
}
real(8) function code(lambda1, lambda2, phi1, phi2)
real(8), intent (in) :: lambda1
real(8), intent (in) :: lambda2
real(8), intent (in) :: phi1
real(8), intent (in) :: phi2
code = atan2(sin((lambda1 - lambda2)), (cos(lambda1) + cos(phi1))) + lambda1
end function
public static double code(double lambda1, double lambda2, double phi1, double phi2) {
return Math.atan2(Math.sin((lambda1 - lambda2)), (Math.cos(lambda1) + Math.cos(phi1))) + lambda1;
}
def code(lambda1, lambda2, phi1, phi2): return math.atan2(math.sin((lambda1 - lambda2)), (math.cos(lambda1) + math.cos(phi1))) + lambda1
function code(lambda1, lambda2, phi1, phi2) return Float64(atan(sin(Float64(lambda1 - lambda2)), Float64(cos(lambda1) + cos(phi1))) + lambda1) end
function tmp = code(lambda1, lambda2, phi1, phi2) tmp = atan2(sin((lambda1 - lambda2)), (cos(lambda1) + cos(phi1))) + lambda1; end
code[lambda1_, lambda2_, phi1_, phi2_] := N[(N[ArcTan[N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision] / N[(N[Cos[lambda1], $MachinePrecision] + N[Cos[phi1], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] + lambda1), $MachinePrecision]
\begin{array}{l}
\\
\tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right)}{\cos \lambda_1 + \cos \phi_1} + \lambda_1
\end{array}
Initial program 99.1%
Taylor expanded in phi2 around 0
lower-sin.f64N/A
lower--.f6477.8
Applied rewrites77.8%
Taylor expanded in phi2 around 0
+-commutativeN/A
sub-negN/A
remove-double-negN/A
distribute-neg-inN/A
+-commutativeN/A
neg-mul-1N/A
cos-negN/A
neg-mul-1N/A
sub-negN/A
lower-+.f64N/A
sub-negN/A
remove-double-negN/A
distribute-neg-inN/A
+-commutativeN/A
neg-mul-1N/A
cos-negN/A
lower-cos.f64N/A
neg-mul-1N/A
sub-negN/A
lower--.f64N/A
lower-cos.f6476.9
Applied rewrites76.9%
Taylor expanded in lambda2 around 0
Applied rewrites68.8%
Final simplification68.8%
(FPCore (lambda1 lambda2 phi1 phi2) :precision binary64 (+ (atan2 (sin lambda1) (+ (cos (- lambda1 lambda2)) (cos phi1))) lambda1))
double code(double lambda1, double lambda2, double phi1, double phi2) {
return atan2(sin(lambda1), (cos((lambda1 - lambda2)) + cos(phi1))) + lambda1;
}
real(8) function code(lambda1, lambda2, phi1, phi2)
real(8), intent (in) :: lambda1
real(8), intent (in) :: lambda2
real(8), intent (in) :: phi1
real(8), intent (in) :: phi2
code = atan2(sin(lambda1), (cos((lambda1 - lambda2)) + cos(phi1))) + lambda1
end function
public static double code(double lambda1, double lambda2, double phi1, double phi2) {
return Math.atan2(Math.sin(lambda1), (Math.cos((lambda1 - lambda2)) + Math.cos(phi1))) + lambda1;
}
def code(lambda1, lambda2, phi1, phi2): return math.atan2(math.sin(lambda1), (math.cos((lambda1 - lambda2)) + math.cos(phi1))) + lambda1
function code(lambda1, lambda2, phi1, phi2) return Float64(atan(sin(lambda1), Float64(cos(Float64(lambda1 - lambda2)) + cos(phi1))) + lambda1) end
function tmp = code(lambda1, lambda2, phi1, phi2) tmp = atan2(sin(lambda1), (cos((lambda1 - lambda2)) + cos(phi1))) + lambda1; end
code[lambda1_, lambda2_, phi1_, phi2_] := N[(N[ArcTan[N[Sin[lambda1], $MachinePrecision] / N[(N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision] + N[Cos[phi1], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] + lambda1), $MachinePrecision]
\begin{array}{l}
\\
\tan^{-1}_* \frac{\sin \lambda_1}{\cos \left(\lambda_1 - \lambda_2\right) + \cos \phi_1} + \lambda_1
\end{array}
Initial program 99.1%
Taylor expanded in phi2 around 0
lower-sin.f64N/A
lower--.f6477.8
Applied rewrites77.8%
Taylor expanded in phi2 around 0
+-commutativeN/A
sub-negN/A
remove-double-negN/A
distribute-neg-inN/A
+-commutativeN/A
neg-mul-1N/A
cos-negN/A
neg-mul-1N/A
sub-negN/A
lower-+.f64N/A
sub-negN/A
remove-double-negN/A
distribute-neg-inN/A
+-commutativeN/A
neg-mul-1N/A
cos-negN/A
lower-cos.f64N/A
neg-mul-1N/A
sub-negN/A
lower--.f64N/A
lower-cos.f6476.9
Applied rewrites76.9%
Taylor expanded in lambda2 around 0
Applied rewrites59.3%
Final simplification59.3%
(FPCore (lambda1 lambda2 phi1 phi2) :precision binary64 (+ (atan2 (sin lambda1) (+ (cos lambda1) (cos phi1))) lambda1))
double code(double lambda1, double lambda2, double phi1, double phi2) {
return atan2(sin(lambda1), (cos(lambda1) + cos(phi1))) + lambda1;
}
real(8) function code(lambda1, lambda2, phi1, phi2)
real(8), intent (in) :: lambda1
real(8), intent (in) :: lambda2
real(8), intent (in) :: phi1
real(8), intent (in) :: phi2
code = atan2(sin(lambda1), (cos(lambda1) + cos(phi1))) + lambda1
end function
public static double code(double lambda1, double lambda2, double phi1, double phi2) {
return Math.atan2(Math.sin(lambda1), (Math.cos(lambda1) + Math.cos(phi1))) + lambda1;
}
def code(lambda1, lambda2, phi1, phi2): return math.atan2(math.sin(lambda1), (math.cos(lambda1) + math.cos(phi1))) + lambda1
function code(lambda1, lambda2, phi1, phi2) return Float64(atan(sin(lambda1), Float64(cos(lambda1) + cos(phi1))) + lambda1) end
function tmp = code(lambda1, lambda2, phi1, phi2) tmp = atan2(sin(lambda1), (cos(lambda1) + cos(phi1))) + lambda1; end
code[lambda1_, lambda2_, phi1_, phi2_] := N[(N[ArcTan[N[Sin[lambda1], $MachinePrecision] / N[(N[Cos[lambda1], $MachinePrecision] + N[Cos[phi1], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] + lambda1), $MachinePrecision]
\begin{array}{l}
\\
\tan^{-1}_* \frac{\sin \lambda_1}{\cos \lambda_1 + \cos \phi_1} + \lambda_1
\end{array}
Initial program 99.1%
Taylor expanded in phi2 around 0
lower-sin.f64N/A
lower--.f6477.8
Applied rewrites77.8%
Taylor expanded in phi2 around 0
+-commutativeN/A
sub-negN/A
remove-double-negN/A
distribute-neg-inN/A
+-commutativeN/A
neg-mul-1N/A
cos-negN/A
neg-mul-1N/A
sub-negN/A
lower-+.f64N/A
sub-negN/A
remove-double-negN/A
distribute-neg-inN/A
+-commutativeN/A
neg-mul-1N/A
cos-negN/A
lower-cos.f64N/A
neg-mul-1N/A
sub-negN/A
lower--.f64N/A
lower-cos.f6476.9
Applied rewrites76.9%
Taylor expanded in lambda2 around 0
Applied rewrites59.3%
Taylor expanded in lambda2 around 0
Applied rewrites59.1%
Final simplification59.1%
herbie shell --seed 2024284
(FPCore (lambda1 lambda2 phi1 phi2)
:name "Midpoint on a great circle"
:precision binary64
(+ lambda1 (atan2 (* (cos phi2) (sin (- lambda1 lambda2))) (+ (cos phi1) (* (cos phi2) (cos (- lambda1 lambda2)))))))