
(FPCore (lambda1 phi1 phi2 delta theta)
:precision binary64
(+
lambda1
(atan2
(* (* (sin theta) (sin delta)) (cos phi1))
(-
(cos delta)
(*
(sin phi1)
(sin
(asin
(+
(* (sin phi1) (cos delta))
(* (* (cos phi1) (sin delta)) (cos theta))))))))))
double code(double lambda1, double phi1, double phi2, double delta, double theta) {
return lambda1 + atan2(((sin(theta) * sin(delta)) * cos(phi1)), (cos(delta) - (sin(phi1) * sin(asin(((sin(phi1) * cos(delta)) + ((cos(phi1) * sin(delta)) * cos(theta))))))));
}
real(8) function code(lambda1, phi1, phi2, delta, theta)
real(8), intent (in) :: lambda1
real(8), intent (in) :: phi1
real(8), intent (in) :: phi2
real(8), intent (in) :: delta
real(8), intent (in) :: theta
code = lambda1 + atan2(((sin(theta) * sin(delta)) * cos(phi1)), (cos(delta) - (sin(phi1) * sin(asin(((sin(phi1) * cos(delta)) + ((cos(phi1) * sin(delta)) * cos(theta))))))))
end function
public static double code(double lambda1, double phi1, double phi2, double delta, double theta) {
return lambda1 + Math.atan2(((Math.sin(theta) * Math.sin(delta)) * Math.cos(phi1)), (Math.cos(delta) - (Math.sin(phi1) * Math.sin(Math.asin(((Math.sin(phi1) * Math.cos(delta)) + ((Math.cos(phi1) * Math.sin(delta)) * Math.cos(theta))))))));
}
def code(lambda1, phi1, phi2, delta, theta): return lambda1 + math.atan2(((math.sin(theta) * math.sin(delta)) * math.cos(phi1)), (math.cos(delta) - (math.sin(phi1) * math.sin(math.asin(((math.sin(phi1) * math.cos(delta)) + ((math.cos(phi1) * math.sin(delta)) * math.cos(theta))))))))
function code(lambda1, phi1, phi2, delta, theta) return Float64(lambda1 + atan(Float64(Float64(sin(theta) * sin(delta)) * cos(phi1)), Float64(cos(delta) - Float64(sin(phi1) * sin(asin(Float64(Float64(sin(phi1) * cos(delta)) + Float64(Float64(cos(phi1) * sin(delta)) * cos(theta))))))))) end
function tmp = code(lambda1, phi1, phi2, delta, theta) tmp = lambda1 + atan2(((sin(theta) * sin(delta)) * cos(phi1)), (cos(delta) - (sin(phi1) * sin(asin(((sin(phi1) * cos(delta)) + ((cos(phi1) * sin(delta)) * cos(theta)))))))); end
code[lambda1_, phi1_, phi2_, delta_, theta_] := N[(lambda1 + N[ArcTan[N[(N[(N[Sin[theta], $MachinePrecision] * N[Sin[delta], $MachinePrecision]), $MachinePrecision] * N[Cos[phi1], $MachinePrecision]), $MachinePrecision] / N[(N[Cos[delta], $MachinePrecision] - N[(N[Sin[phi1], $MachinePrecision] * N[Sin[N[ArcSin[N[(N[(N[Sin[phi1], $MachinePrecision] * N[Cos[delta], $MachinePrecision]), $MachinePrecision] + N[(N[(N[Cos[phi1], $MachinePrecision] * N[Sin[delta], $MachinePrecision]), $MachinePrecision] * N[Cos[theta], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\lambda_1 + \tan^{-1}_* \frac{\left(\sin theta \cdot \sin delta\right) \cdot \cos \phi_1}{\cos delta - \sin \phi_1 \cdot \sin \sin^{-1} \left(\sin \phi_1 \cdot \cos delta + \left(\cos \phi_1 \cdot \sin delta\right) \cdot \cos theta\right)}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 17 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (lambda1 phi1 phi2 delta theta)
:precision binary64
(+
lambda1
(atan2
(* (* (sin theta) (sin delta)) (cos phi1))
(-
(cos delta)
(*
(sin phi1)
(sin
(asin
(+
(* (sin phi1) (cos delta))
(* (* (cos phi1) (sin delta)) (cos theta))))))))))
double code(double lambda1, double phi1, double phi2, double delta, double theta) {
return lambda1 + atan2(((sin(theta) * sin(delta)) * cos(phi1)), (cos(delta) - (sin(phi1) * sin(asin(((sin(phi1) * cos(delta)) + ((cos(phi1) * sin(delta)) * cos(theta))))))));
}
real(8) function code(lambda1, phi1, phi2, delta, theta)
real(8), intent (in) :: lambda1
real(8), intent (in) :: phi1
real(8), intent (in) :: phi2
real(8), intent (in) :: delta
real(8), intent (in) :: theta
code = lambda1 + atan2(((sin(theta) * sin(delta)) * cos(phi1)), (cos(delta) - (sin(phi1) * sin(asin(((sin(phi1) * cos(delta)) + ((cos(phi1) * sin(delta)) * cos(theta))))))))
end function
public static double code(double lambda1, double phi1, double phi2, double delta, double theta) {
return lambda1 + Math.atan2(((Math.sin(theta) * Math.sin(delta)) * Math.cos(phi1)), (Math.cos(delta) - (Math.sin(phi1) * Math.sin(Math.asin(((Math.sin(phi1) * Math.cos(delta)) + ((Math.cos(phi1) * Math.sin(delta)) * Math.cos(theta))))))));
}
def code(lambda1, phi1, phi2, delta, theta): return lambda1 + math.atan2(((math.sin(theta) * math.sin(delta)) * math.cos(phi1)), (math.cos(delta) - (math.sin(phi1) * math.sin(math.asin(((math.sin(phi1) * math.cos(delta)) + ((math.cos(phi1) * math.sin(delta)) * math.cos(theta))))))))
function code(lambda1, phi1, phi2, delta, theta) return Float64(lambda1 + atan(Float64(Float64(sin(theta) * sin(delta)) * cos(phi1)), Float64(cos(delta) - Float64(sin(phi1) * sin(asin(Float64(Float64(sin(phi1) * cos(delta)) + Float64(Float64(cos(phi1) * sin(delta)) * cos(theta))))))))) end
function tmp = code(lambda1, phi1, phi2, delta, theta) tmp = lambda1 + atan2(((sin(theta) * sin(delta)) * cos(phi1)), (cos(delta) - (sin(phi1) * sin(asin(((sin(phi1) * cos(delta)) + ((cos(phi1) * sin(delta)) * cos(theta)))))))); end
code[lambda1_, phi1_, phi2_, delta_, theta_] := N[(lambda1 + N[ArcTan[N[(N[(N[Sin[theta], $MachinePrecision] * N[Sin[delta], $MachinePrecision]), $MachinePrecision] * N[Cos[phi1], $MachinePrecision]), $MachinePrecision] / N[(N[Cos[delta], $MachinePrecision] - N[(N[Sin[phi1], $MachinePrecision] * N[Sin[N[ArcSin[N[(N[(N[Sin[phi1], $MachinePrecision] * N[Cos[delta], $MachinePrecision]), $MachinePrecision] + N[(N[(N[Cos[phi1], $MachinePrecision] * N[Sin[delta], $MachinePrecision]), $MachinePrecision] * N[Cos[theta], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\lambda_1 + \tan^{-1}_* \frac{\left(\sin theta \cdot \sin delta\right) \cdot \cos \phi_1}{\cos delta - \sin \phi_1 \cdot \sin \sin^{-1} \left(\sin \phi_1 \cdot \cos delta + \left(\cos \phi_1 \cdot \sin delta\right) \cdot \cos theta\right)}
\end{array}
(FPCore (lambda1 phi1 phi2 delta theta)
:precision binary64
(+
(atan2
(* (cos phi1) (* (sin delta) (sin theta)))
(fma
(+ 0.5 (* 0.5 (cos (* phi1 -2.0))))
(cos delta)
(* (* (* (cos phi1) (sin delta)) (cos theta)) (- (sin phi1)))))
lambda1))
double code(double lambda1, double phi1, double phi2, double delta, double theta) {
return atan2((cos(phi1) * (sin(delta) * sin(theta))), fma((0.5 + (0.5 * cos((phi1 * -2.0)))), cos(delta), (((cos(phi1) * sin(delta)) * cos(theta)) * -sin(phi1)))) + lambda1;
}
function code(lambda1, phi1, phi2, delta, theta) return Float64(atan(Float64(cos(phi1) * Float64(sin(delta) * sin(theta))), fma(Float64(0.5 + Float64(0.5 * cos(Float64(phi1 * -2.0)))), cos(delta), Float64(Float64(Float64(cos(phi1) * sin(delta)) * cos(theta)) * Float64(-sin(phi1))))) + lambda1) end
code[lambda1_, phi1_, phi2_, delta_, theta_] := N[(N[ArcTan[N[(N[Cos[phi1], $MachinePrecision] * N[(N[Sin[delta], $MachinePrecision] * N[Sin[theta], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(N[(0.5 + N[(0.5 * N[Cos[N[(phi1 * -2.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[Cos[delta], $MachinePrecision] + N[(N[(N[(N[Cos[phi1], $MachinePrecision] * N[Sin[delta], $MachinePrecision]), $MachinePrecision] * N[Cos[theta], $MachinePrecision]), $MachinePrecision] * (-N[Sin[phi1], $MachinePrecision])), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] + lambda1), $MachinePrecision]
\begin{array}{l}
\\
\tan^{-1}_* \frac{\cos \phi_1 \cdot \left(\sin delta \cdot \sin theta\right)}{\mathsf{fma}\left(0.5 + 0.5 \cdot \cos \left(\phi_1 \cdot -2\right), \cos delta, \left(\left(\cos \phi_1 \cdot \sin delta\right) \cdot \cos theta\right) \cdot \left(-\sin \phi_1\right)\right)} + \lambda_1
\end{array}
Initial program 99.8%
sub-negN/A
+-commutativeN/A
sin-asinN/A
distribute-lft-neg-inN/A
distribute-rgt-inN/A
associate-+l+N/A
Applied egg-rr99.8%
distribute-lft-neg-outN/A
neg-lowering-neg.f64N/A
sqr-sin-aN/A
sub-negN/A
+-commutativeN/A
*-commutativeN/A
distribute-rgt-neg-inN/A
accelerator-lowering-fma.f64N/A
cos-2N/A
cos-sumN/A
cos-lowering-cos.f64N/A
+-lowering-+.f64N/A
metadata-eval99.8
Applied egg-rr99.8%
Taylor expanded in lambda1 around 0
Simplified99.8%
Final simplification99.8%
(FPCore (lambda1 phi1 phi2 delta theta)
:precision binary64
(let* ((t_1 (atan2 (* (sin delta) (* (cos phi1) (sin theta))) (cos delta)))
(t_2 (* (cos phi1) (* (sin delta) (sin theta))))
(t_3 (+ lambda1 (atan2 t_2 (fma delta (* delta -0.5) 1.0))))
(t_4 (* (cos phi1) (sin delta)))
(t_5
(+
lambda1
(atan2
t_2
(-
(cos delta)
(*
(sin phi1)
(sin
(asin (+ (* t_4 (cos theta)) (* (cos delta) (sin phi1)))))))))))
(if (<= t_5 -3.14159265358979)
t_3
(if (<= t_5 -5e-7)
t_1
(if (<= t_5 2e-12)
(* lambda1 (- (/ (atan2 (* (sin theta) t_4) 1.0) lambda1) -1.0))
(if (<= t_5 3.12) t_1 t_3))))))
double code(double lambda1, double phi1, double phi2, double delta, double theta) {
double t_1 = atan2((sin(delta) * (cos(phi1) * sin(theta))), cos(delta));
double t_2 = cos(phi1) * (sin(delta) * sin(theta));
double t_3 = lambda1 + atan2(t_2, fma(delta, (delta * -0.5), 1.0));
double t_4 = cos(phi1) * sin(delta);
double t_5 = lambda1 + atan2(t_2, (cos(delta) - (sin(phi1) * sin(asin(((t_4 * cos(theta)) + (cos(delta) * sin(phi1))))))));
double tmp;
if (t_5 <= -3.14159265358979) {
tmp = t_3;
} else if (t_5 <= -5e-7) {
tmp = t_1;
} else if (t_5 <= 2e-12) {
tmp = lambda1 * ((atan2((sin(theta) * t_4), 1.0) / lambda1) - -1.0);
} else if (t_5 <= 3.12) {
tmp = t_1;
} else {
tmp = t_3;
}
return tmp;
}
function code(lambda1, phi1, phi2, delta, theta) t_1 = atan(Float64(sin(delta) * Float64(cos(phi1) * sin(theta))), cos(delta)) t_2 = Float64(cos(phi1) * Float64(sin(delta) * sin(theta))) t_3 = Float64(lambda1 + atan(t_2, fma(delta, Float64(delta * -0.5), 1.0))) t_4 = Float64(cos(phi1) * sin(delta)) t_5 = Float64(lambda1 + atan(t_2, Float64(cos(delta) - Float64(sin(phi1) * sin(asin(Float64(Float64(t_4 * cos(theta)) + Float64(cos(delta) * sin(phi1))))))))) tmp = 0.0 if (t_5 <= -3.14159265358979) tmp = t_3; elseif (t_5 <= -5e-7) tmp = t_1; elseif (t_5 <= 2e-12) tmp = Float64(lambda1 * Float64(Float64(atan(Float64(sin(theta) * t_4), 1.0) / lambda1) - -1.0)); elseif (t_5 <= 3.12) tmp = t_1; else tmp = t_3; end return tmp end
code[lambda1_, phi1_, phi2_, delta_, theta_] := Block[{t$95$1 = N[ArcTan[N[(N[Sin[delta], $MachinePrecision] * N[(N[Cos[phi1], $MachinePrecision] * N[Sin[theta], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[Cos[delta], $MachinePrecision]], $MachinePrecision]}, Block[{t$95$2 = N[(N[Cos[phi1], $MachinePrecision] * N[(N[Sin[delta], $MachinePrecision] * N[Sin[theta], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = N[(lambda1 + N[ArcTan[t$95$2 / N[(delta * N[(delta * -0.5), $MachinePrecision] + 1.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$4 = N[(N[Cos[phi1], $MachinePrecision] * N[Sin[delta], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$5 = N[(lambda1 + N[ArcTan[t$95$2 / N[(N[Cos[delta], $MachinePrecision] - N[(N[Sin[phi1], $MachinePrecision] * N[Sin[N[ArcSin[N[(N[(t$95$4 * N[Cos[theta], $MachinePrecision]), $MachinePrecision] + N[(N[Cos[delta], $MachinePrecision] * N[Sin[phi1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$5, -3.14159265358979], t$95$3, If[LessEqual[t$95$5, -5e-7], t$95$1, If[LessEqual[t$95$5, 2e-12], N[(lambda1 * N[(N[(N[ArcTan[N[(N[Sin[theta], $MachinePrecision] * t$95$4), $MachinePrecision] / 1.0], $MachinePrecision] / lambda1), $MachinePrecision] - -1.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$5, 3.12], t$95$1, t$95$3]]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \tan^{-1}_* \frac{\sin delta \cdot \left(\cos \phi_1 \cdot \sin theta\right)}{\cos delta}\\
t_2 := \cos \phi_1 \cdot \left(\sin delta \cdot \sin theta\right)\\
t_3 := \lambda_1 + \tan^{-1}_* \frac{t\_2}{\mathsf{fma}\left(delta, delta \cdot -0.5, 1\right)}\\
t_4 := \cos \phi_1 \cdot \sin delta\\
t_5 := \lambda_1 + \tan^{-1}_* \frac{t\_2}{\cos delta - \sin \phi_1 \cdot \sin \sin^{-1} \left(t\_4 \cdot \cos theta + \cos delta \cdot \sin \phi_1\right)}\\
\mathbf{if}\;t\_5 \leq -3.14159265358979:\\
\;\;\;\;t\_3\\
\mathbf{elif}\;t\_5 \leq -5 \cdot 10^{-7}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t\_5 \leq 2 \cdot 10^{-12}:\\
\;\;\;\;\lambda_1 \cdot \left(\frac{\tan^{-1}_* \frac{\sin theta \cdot t\_4}{1}}{\lambda_1} - -1\right)\\
\mathbf{elif}\;t\_5 \leq 3.12:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;t\_3\\
\end{array}
\end{array}
if (+.f64 lambda1 (atan2.f64 (*.f64 (*.f64 (sin.f64 theta) (sin.f64 delta)) (cos.f64 phi1)) (-.f64 (cos.f64 delta) (*.f64 (sin.f64 phi1) (sin.f64 (asin.f64 (+.f64 (*.f64 (sin.f64 phi1) (cos.f64 delta)) (*.f64 (*.f64 (cos.f64 phi1) (sin.f64 delta)) (cos.f64 theta))))))))) < -3.14159265358979001 or 3.1200000000000001 < (+.f64 lambda1 (atan2.f64 (*.f64 (*.f64 (sin.f64 theta) (sin.f64 delta)) (cos.f64 phi1)) (-.f64 (cos.f64 delta) (*.f64 (sin.f64 phi1) (sin.f64 (asin.f64 (+.f64 (*.f64 (sin.f64 phi1) (cos.f64 delta)) (*.f64 (*.f64 (cos.f64 phi1) (sin.f64 delta)) (cos.f64 theta))))))))) Initial program 100.0%
sin-asinN/A
distribute-rgt-inN/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
accelerator-lowering-fma.f64N/A
sqr-sin-aN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
*-lowering-*.f64N/A
Applied egg-rr100.0%
Taylor expanded in delta around 0
associate-+r+N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
Simplified99.4%
Taylor expanded in phi1 around 0
+-commutativeN/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
*-lowering-*.f6499.4
Simplified99.4%
if -3.14159265358979001 < (+.f64 lambda1 (atan2.f64 (*.f64 (*.f64 (sin.f64 theta) (sin.f64 delta)) (cos.f64 phi1)) (-.f64 (cos.f64 delta) (*.f64 (sin.f64 phi1) (sin.f64 (asin.f64 (+.f64 (*.f64 (sin.f64 phi1) (cos.f64 delta)) (*.f64 (*.f64 (cos.f64 phi1) (sin.f64 delta)) (cos.f64 theta))))))))) < -4.99999999999999977e-7 or 1.99999999999999996e-12 < (+.f64 lambda1 (atan2.f64 (*.f64 (*.f64 (sin.f64 theta) (sin.f64 delta)) (cos.f64 phi1)) (-.f64 (cos.f64 delta) (*.f64 (sin.f64 phi1) (sin.f64 (asin.f64 (+.f64 (*.f64 (sin.f64 phi1) (cos.f64 delta)) (*.f64 (*.f64 (cos.f64 phi1) (sin.f64 delta)) (cos.f64 theta))))))))) < 3.1200000000000001Initial program 99.3%
Taylor expanded in lambda1 around 0
atan2-lowering-atan2.f64N/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
sin-lowering-sin.f64N/A
sub-negN/A
+-commutativeN/A
Simplified92.6%
Taylor expanded in phi1 around 0
cos-lowering-cos.f6461.3
Simplified61.3%
if -4.99999999999999977e-7 < (+.f64 lambda1 (atan2.f64 (*.f64 (*.f64 (sin.f64 theta) (sin.f64 delta)) (cos.f64 phi1)) (-.f64 (cos.f64 delta) (*.f64 (sin.f64 phi1) (sin.f64 (asin.f64 (+.f64 (*.f64 (sin.f64 phi1) (cos.f64 delta)) (*.f64 (*.f64 (cos.f64 phi1) (sin.f64 delta)) (cos.f64 theta))))))))) < 1.99999999999999996e-12Initial program 99.5%
Taylor expanded in delta around 0
+-commutativeN/A
associate--l+N/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
1-sub-sinN/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
mul-1-negN/A
distribute-lft-neg-outN/A
distribute-rgt-neg-inN/A
accelerator-lowering-fma.f64N/A
Simplified88.4%
Taylor expanded in phi1 around 0
Simplified80.9%
Taylor expanded in lambda1 around -inf
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
Simplified80.9%
Final simplification88.2%
(FPCore (lambda1 phi1 phi2 delta theta)
:precision binary64
(let* ((t_1 (* (sin delta) (sin theta)))
(t_2 (atan2 t_1 (cos delta)))
(t_3 (* (cos phi1) t_1))
(t_4 (+ lambda1 (atan2 t_3 (fma delta (* delta -0.5) 1.0))))
(t_5 (* (cos phi1) (sin delta)))
(t_6
(+
lambda1
(atan2
t_3
(-
(cos delta)
(*
(sin phi1)
(sin
(asin (+ (* t_5 (cos theta)) (* (cos delta) (sin phi1)))))))))))
(if (<= t_6 -3.14159265358979)
t_4
(if (<= t_6 -5e-7)
t_2
(if (<= t_6 2e-12)
(* lambda1 (- (/ (atan2 (* (sin theta) t_5) 1.0) lambda1) -1.0))
(if (<= t_6 3.0) t_2 t_4))))))
double code(double lambda1, double phi1, double phi2, double delta, double theta) {
double t_1 = sin(delta) * sin(theta);
double t_2 = atan2(t_1, cos(delta));
double t_3 = cos(phi1) * t_1;
double t_4 = lambda1 + atan2(t_3, fma(delta, (delta * -0.5), 1.0));
double t_5 = cos(phi1) * sin(delta);
double t_6 = lambda1 + atan2(t_3, (cos(delta) - (sin(phi1) * sin(asin(((t_5 * cos(theta)) + (cos(delta) * sin(phi1))))))));
double tmp;
if (t_6 <= -3.14159265358979) {
tmp = t_4;
} else if (t_6 <= -5e-7) {
tmp = t_2;
} else if (t_6 <= 2e-12) {
tmp = lambda1 * ((atan2((sin(theta) * t_5), 1.0) / lambda1) - -1.0);
} else if (t_6 <= 3.0) {
tmp = t_2;
} else {
tmp = t_4;
}
return tmp;
}
function code(lambda1, phi1, phi2, delta, theta) t_1 = Float64(sin(delta) * sin(theta)) t_2 = atan(t_1, cos(delta)) t_3 = Float64(cos(phi1) * t_1) t_4 = Float64(lambda1 + atan(t_3, fma(delta, Float64(delta * -0.5), 1.0))) t_5 = Float64(cos(phi1) * sin(delta)) t_6 = Float64(lambda1 + atan(t_3, Float64(cos(delta) - Float64(sin(phi1) * sin(asin(Float64(Float64(t_5 * cos(theta)) + Float64(cos(delta) * sin(phi1))))))))) tmp = 0.0 if (t_6 <= -3.14159265358979) tmp = t_4; elseif (t_6 <= -5e-7) tmp = t_2; elseif (t_6 <= 2e-12) tmp = Float64(lambda1 * Float64(Float64(atan(Float64(sin(theta) * t_5), 1.0) / lambda1) - -1.0)); elseif (t_6 <= 3.0) tmp = t_2; else tmp = t_4; end return tmp end
code[lambda1_, phi1_, phi2_, delta_, theta_] := Block[{t$95$1 = N[(N[Sin[delta], $MachinePrecision] * N[Sin[theta], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[ArcTan[t$95$1 / N[Cos[delta], $MachinePrecision]], $MachinePrecision]}, Block[{t$95$3 = N[(N[Cos[phi1], $MachinePrecision] * t$95$1), $MachinePrecision]}, Block[{t$95$4 = N[(lambda1 + N[ArcTan[t$95$3 / N[(delta * N[(delta * -0.5), $MachinePrecision] + 1.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$5 = N[(N[Cos[phi1], $MachinePrecision] * N[Sin[delta], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$6 = N[(lambda1 + N[ArcTan[t$95$3 / N[(N[Cos[delta], $MachinePrecision] - N[(N[Sin[phi1], $MachinePrecision] * N[Sin[N[ArcSin[N[(N[(t$95$5 * N[Cos[theta], $MachinePrecision]), $MachinePrecision] + N[(N[Cos[delta], $MachinePrecision] * N[Sin[phi1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$6, -3.14159265358979], t$95$4, If[LessEqual[t$95$6, -5e-7], t$95$2, If[LessEqual[t$95$6, 2e-12], N[(lambda1 * N[(N[(N[ArcTan[N[(N[Sin[theta], $MachinePrecision] * t$95$5), $MachinePrecision] / 1.0], $MachinePrecision] / lambda1), $MachinePrecision] - -1.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$6, 3.0], t$95$2, t$95$4]]]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \sin delta \cdot \sin theta\\
t_2 := \tan^{-1}_* \frac{t\_1}{\cos delta}\\
t_3 := \cos \phi_1 \cdot t\_1\\
t_4 := \lambda_1 + \tan^{-1}_* \frac{t\_3}{\mathsf{fma}\left(delta, delta \cdot -0.5, 1\right)}\\
t_5 := \cos \phi_1 \cdot \sin delta\\
t_6 := \lambda_1 + \tan^{-1}_* \frac{t\_3}{\cos delta - \sin \phi_1 \cdot \sin \sin^{-1} \left(t\_5 \cdot \cos theta + \cos delta \cdot \sin \phi_1\right)}\\
\mathbf{if}\;t\_6 \leq -3.14159265358979:\\
\;\;\;\;t\_4\\
\mathbf{elif}\;t\_6 \leq -5 \cdot 10^{-7}:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;t\_6 \leq 2 \cdot 10^{-12}:\\
\;\;\;\;\lambda_1 \cdot \left(\frac{\tan^{-1}_* \frac{\sin theta \cdot t\_5}{1}}{\lambda_1} - -1\right)\\
\mathbf{elif}\;t\_6 \leq 3:\\
\;\;\;\;t\_2\\
\mathbf{else}:\\
\;\;\;\;t\_4\\
\end{array}
\end{array}
if (+.f64 lambda1 (atan2.f64 (*.f64 (*.f64 (sin.f64 theta) (sin.f64 delta)) (cos.f64 phi1)) (-.f64 (cos.f64 delta) (*.f64 (sin.f64 phi1) (sin.f64 (asin.f64 (+.f64 (*.f64 (sin.f64 phi1) (cos.f64 delta)) (*.f64 (*.f64 (cos.f64 phi1) (sin.f64 delta)) (cos.f64 theta))))))))) < -3.14159265358979001 or 3 < (+.f64 lambda1 (atan2.f64 (*.f64 (*.f64 (sin.f64 theta) (sin.f64 delta)) (cos.f64 phi1)) (-.f64 (cos.f64 delta) (*.f64 (sin.f64 phi1) (sin.f64 (asin.f64 (+.f64 (*.f64 (sin.f64 phi1) (cos.f64 delta)) (*.f64 (*.f64 (cos.f64 phi1) (sin.f64 delta)) (cos.f64 theta))))))))) Initial program 100.0%
sin-asinN/A
distribute-rgt-inN/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
accelerator-lowering-fma.f64N/A
sqr-sin-aN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
*-lowering-*.f64N/A
Applied egg-rr100.0%
Taylor expanded in delta around 0
associate-+r+N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
Simplified98.9%
Taylor expanded in phi1 around 0
+-commutativeN/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
*-lowering-*.f6498.9
Simplified98.9%
if -3.14159265358979001 < (+.f64 lambda1 (atan2.f64 (*.f64 (*.f64 (sin.f64 theta) (sin.f64 delta)) (cos.f64 phi1)) (-.f64 (cos.f64 delta) (*.f64 (sin.f64 phi1) (sin.f64 (asin.f64 (+.f64 (*.f64 (sin.f64 phi1) (cos.f64 delta)) (*.f64 (*.f64 (cos.f64 phi1) (sin.f64 delta)) (cos.f64 theta))))))))) < -4.99999999999999977e-7 or 1.99999999999999996e-12 < (+.f64 lambda1 (atan2.f64 (*.f64 (*.f64 (sin.f64 theta) (sin.f64 delta)) (cos.f64 phi1)) (-.f64 (cos.f64 delta) (*.f64 (sin.f64 phi1) (sin.f64 (asin.f64 (+.f64 (*.f64 (sin.f64 phi1) (cos.f64 delta)) (*.f64 (*.f64 (cos.f64 phi1) (sin.f64 delta)) (cos.f64 theta))))))))) < 3Initial program 99.3%
Taylor expanded in lambda1 around 0
atan2-lowering-atan2.f64N/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
sin-lowering-sin.f64N/A
sub-negN/A
+-commutativeN/A
Simplified92.4%
Taylor expanded in phi1 around 0
cos-lowering-cos.f6462.0
Simplified62.0%
Taylor expanded in phi1 around 0
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
sin-lowering-sin.f6458.3
Simplified58.3%
if -4.99999999999999977e-7 < (+.f64 lambda1 (atan2.f64 (*.f64 (*.f64 (sin.f64 theta) (sin.f64 delta)) (cos.f64 phi1)) (-.f64 (cos.f64 delta) (*.f64 (sin.f64 phi1) (sin.f64 (asin.f64 (+.f64 (*.f64 (sin.f64 phi1) (cos.f64 delta)) (*.f64 (*.f64 (cos.f64 phi1) (sin.f64 delta)) (cos.f64 theta))))))))) < 1.99999999999999996e-12Initial program 99.5%
Taylor expanded in delta around 0
+-commutativeN/A
associate--l+N/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
1-sub-sinN/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
mul-1-negN/A
distribute-lft-neg-outN/A
distribute-rgt-neg-inN/A
accelerator-lowering-fma.f64N/A
Simplified88.4%
Taylor expanded in phi1 around 0
Simplified80.9%
Taylor expanded in lambda1 around -inf
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
Simplified80.9%
Final simplification87.6%
(FPCore (lambda1 phi1 phi2 delta theta)
:precision binary64
(let* ((t_1 (* (sin delta) (sin theta)))
(t_2 (* (cos phi1) t_1))
(t_3 (+ lambda1 (atan2 t_2 (fma delta (* delta -0.5) 1.0))))
(t_4
(+
lambda1
(atan2
t_2
(-
(cos delta)
(*
(sin phi1)
(sin
(asin
(+
(* (* (cos phi1) (sin delta)) (cos theta))
(* (cos delta) (sin phi1))))))))))
(t_5 (atan2 t_1 (cos delta))))
(if (<= t_4 -3.14159265358979)
t_3
(if (<= t_4 -5e-7)
t_5
(if (<= t_4 2e-12)
(+ lambda1 (atan2 t_2 1.0))
(if (<= t_4 3.0) t_5 t_3))))))
double code(double lambda1, double phi1, double phi2, double delta, double theta) {
double t_1 = sin(delta) * sin(theta);
double t_2 = cos(phi1) * t_1;
double t_3 = lambda1 + atan2(t_2, fma(delta, (delta * -0.5), 1.0));
double t_4 = lambda1 + atan2(t_2, (cos(delta) - (sin(phi1) * sin(asin((((cos(phi1) * sin(delta)) * cos(theta)) + (cos(delta) * sin(phi1))))))));
double t_5 = atan2(t_1, cos(delta));
double tmp;
if (t_4 <= -3.14159265358979) {
tmp = t_3;
} else if (t_4 <= -5e-7) {
tmp = t_5;
} else if (t_4 <= 2e-12) {
tmp = lambda1 + atan2(t_2, 1.0);
} else if (t_4 <= 3.0) {
tmp = t_5;
} else {
tmp = t_3;
}
return tmp;
}
function code(lambda1, phi1, phi2, delta, theta) t_1 = Float64(sin(delta) * sin(theta)) t_2 = Float64(cos(phi1) * t_1) t_3 = Float64(lambda1 + atan(t_2, fma(delta, Float64(delta * -0.5), 1.0))) t_4 = Float64(lambda1 + atan(t_2, Float64(cos(delta) - Float64(sin(phi1) * sin(asin(Float64(Float64(Float64(cos(phi1) * sin(delta)) * cos(theta)) + Float64(cos(delta) * sin(phi1))))))))) t_5 = atan(t_1, cos(delta)) tmp = 0.0 if (t_4 <= -3.14159265358979) tmp = t_3; elseif (t_4 <= -5e-7) tmp = t_5; elseif (t_4 <= 2e-12) tmp = Float64(lambda1 + atan(t_2, 1.0)); elseif (t_4 <= 3.0) tmp = t_5; else tmp = t_3; end return tmp end
code[lambda1_, phi1_, phi2_, delta_, theta_] := Block[{t$95$1 = N[(N[Sin[delta], $MachinePrecision] * N[Sin[theta], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(N[Cos[phi1], $MachinePrecision] * t$95$1), $MachinePrecision]}, Block[{t$95$3 = N[(lambda1 + N[ArcTan[t$95$2 / N[(delta * N[(delta * -0.5), $MachinePrecision] + 1.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$4 = N[(lambda1 + N[ArcTan[t$95$2 / N[(N[Cos[delta], $MachinePrecision] - N[(N[Sin[phi1], $MachinePrecision] * N[Sin[N[ArcSin[N[(N[(N[(N[Cos[phi1], $MachinePrecision] * N[Sin[delta], $MachinePrecision]), $MachinePrecision] * N[Cos[theta], $MachinePrecision]), $MachinePrecision] + N[(N[Cos[delta], $MachinePrecision] * N[Sin[phi1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$5 = N[ArcTan[t$95$1 / N[Cos[delta], $MachinePrecision]], $MachinePrecision]}, If[LessEqual[t$95$4, -3.14159265358979], t$95$3, If[LessEqual[t$95$4, -5e-7], t$95$5, If[LessEqual[t$95$4, 2e-12], N[(lambda1 + N[ArcTan[t$95$2 / 1.0], $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$4, 3.0], t$95$5, t$95$3]]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \sin delta \cdot \sin theta\\
t_2 := \cos \phi_1 \cdot t\_1\\
t_3 := \lambda_1 + \tan^{-1}_* \frac{t\_2}{\mathsf{fma}\left(delta, delta \cdot -0.5, 1\right)}\\
t_4 := \lambda_1 + \tan^{-1}_* \frac{t\_2}{\cos delta - \sin \phi_1 \cdot \sin \sin^{-1} \left(\left(\cos \phi_1 \cdot \sin delta\right) \cdot \cos theta + \cos delta \cdot \sin \phi_1\right)}\\
t_5 := \tan^{-1}_* \frac{t\_1}{\cos delta}\\
\mathbf{if}\;t\_4 \leq -3.14159265358979:\\
\;\;\;\;t\_3\\
\mathbf{elif}\;t\_4 \leq -5 \cdot 10^{-7}:\\
\;\;\;\;t\_5\\
\mathbf{elif}\;t\_4 \leq 2 \cdot 10^{-12}:\\
\;\;\;\;\lambda_1 + \tan^{-1}_* \frac{t\_2}{1}\\
\mathbf{elif}\;t\_4 \leq 3:\\
\;\;\;\;t\_5\\
\mathbf{else}:\\
\;\;\;\;t\_3\\
\end{array}
\end{array}
if (+.f64 lambda1 (atan2.f64 (*.f64 (*.f64 (sin.f64 theta) (sin.f64 delta)) (cos.f64 phi1)) (-.f64 (cos.f64 delta) (*.f64 (sin.f64 phi1) (sin.f64 (asin.f64 (+.f64 (*.f64 (sin.f64 phi1) (cos.f64 delta)) (*.f64 (*.f64 (cos.f64 phi1) (sin.f64 delta)) (cos.f64 theta))))))))) < -3.14159265358979001 or 3 < (+.f64 lambda1 (atan2.f64 (*.f64 (*.f64 (sin.f64 theta) (sin.f64 delta)) (cos.f64 phi1)) (-.f64 (cos.f64 delta) (*.f64 (sin.f64 phi1) (sin.f64 (asin.f64 (+.f64 (*.f64 (sin.f64 phi1) (cos.f64 delta)) (*.f64 (*.f64 (cos.f64 phi1) (sin.f64 delta)) (cos.f64 theta))))))))) Initial program 100.0%
sin-asinN/A
distribute-rgt-inN/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
accelerator-lowering-fma.f64N/A
sqr-sin-aN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
*-lowering-*.f64N/A
Applied egg-rr100.0%
Taylor expanded in delta around 0
associate-+r+N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
Simplified98.9%
Taylor expanded in phi1 around 0
+-commutativeN/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
accelerator-lowering-fma.f64N/A
*-commutativeN/A
*-lowering-*.f6498.9
Simplified98.9%
if -3.14159265358979001 < (+.f64 lambda1 (atan2.f64 (*.f64 (*.f64 (sin.f64 theta) (sin.f64 delta)) (cos.f64 phi1)) (-.f64 (cos.f64 delta) (*.f64 (sin.f64 phi1) (sin.f64 (asin.f64 (+.f64 (*.f64 (sin.f64 phi1) (cos.f64 delta)) (*.f64 (*.f64 (cos.f64 phi1) (sin.f64 delta)) (cos.f64 theta))))))))) < -4.99999999999999977e-7 or 1.99999999999999996e-12 < (+.f64 lambda1 (atan2.f64 (*.f64 (*.f64 (sin.f64 theta) (sin.f64 delta)) (cos.f64 phi1)) (-.f64 (cos.f64 delta) (*.f64 (sin.f64 phi1) (sin.f64 (asin.f64 (+.f64 (*.f64 (sin.f64 phi1) (cos.f64 delta)) (*.f64 (*.f64 (cos.f64 phi1) (sin.f64 delta)) (cos.f64 theta))))))))) < 3Initial program 99.3%
Taylor expanded in lambda1 around 0
atan2-lowering-atan2.f64N/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
sin-lowering-sin.f64N/A
sub-negN/A
+-commutativeN/A
Simplified92.4%
Taylor expanded in phi1 around 0
cos-lowering-cos.f6462.0
Simplified62.0%
Taylor expanded in phi1 around 0
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
sin-lowering-sin.f6458.3
Simplified58.3%
if -4.99999999999999977e-7 < (+.f64 lambda1 (atan2.f64 (*.f64 (*.f64 (sin.f64 theta) (sin.f64 delta)) (cos.f64 phi1)) (-.f64 (cos.f64 delta) (*.f64 (sin.f64 phi1) (sin.f64 (asin.f64 (+.f64 (*.f64 (sin.f64 phi1) (cos.f64 delta)) (*.f64 (*.f64 (cos.f64 phi1) (sin.f64 delta)) (cos.f64 theta))))))))) < 1.99999999999999996e-12Initial program 99.5%
Taylor expanded in delta around 0
+-commutativeN/A
associate--l+N/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
1-sub-sinN/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
mul-1-negN/A
distribute-lft-neg-outN/A
distribute-rgt-neg-inN/A
accelerator-lowering-fma.f64N/A
Simplified88.4%
Taylor expanded in phi1 around 0
Simplified80.9%
Final simplification87.6%
(FPCore (lambda1 phi1 phi2 delta theta)
:precision binary64
(let* ((t_1 (* (sin delta) (sin theta)))
(t_2 (atan2 t_1 (cos delta)))
(t_3 (* (cos phi1) t_1))
(t_4
(+
lambda1
(atan2
t_3
(-
(cos delta)
(*
(sin phi1)
(sin
(asin
(+
(* (* (cos phi1) (sin delta)) (cos theta))
(* (cos delta) (sin phi1)))))))))))
(if (<= t_4 -500.0)
(+
lambda1
(atan2
(*
(sin theta)
(fma delta (* -0.16666666666666666 (* delta delta)) delta))
1.0))
(if (<= t_4 -5e-7)
t_2
(if (<= t_4 2e-12)
(+ lambda1 (atan2 t_3 1.0))
(if (<= t_4 5.0)
t_2
(+ lambda1 (atan2 (* delta (sin theta)) 1.0))))))))
double code(double lambda1, double phi1, double phi2, double delta, double theta) {
double t_1 = sin(delta) * sin(theta);
double t_2 = atan2(t_1, cos(delta));
double t_3 = cos(phi1) * t_1;
double t_4 = lambda1 + atan2(t_3, (cos(delta) - (sin(phi1) * sin(asin((((cos(phi1) * sin(delta)) * cos(theta)) + (cos(delta) * sin(phi1))))))));
double tmp;
if (t_4 <= -500.0) {
tmp = lambda1 + atan2((sin(theta) * fma(delta, (-0.16666666666666666 * (delta * delta)), delta)), 1.0);
} else if (t_4 <= -5e-7) {
tmp = t_2;
} else if (t_4 <= 2e-12) {
tmp = lambda1 + atan2(t_3, 1.0);
} else if (t_4 <= 5.0) {
tmp = t_2;
} else {
tmp = lambda1 + atan2((delta * sin(theta)), 1.0);
}
return tmp;
}
function code(lambda1, phi1, phi2, delta, theta) t_1 = Float64(sin(delta) * sin(theta)) t_2 = atan(t_1, cos(delta)) t_3 = Float64(cos(phi1) * t_1) t_4 = Float64(lambda1 + atan(t_3, Float64(cos(delta) - Float64(sin(phi1) * sin(asin(Float64(Float64(Float64(cos(phi1) * sin(delta)) * cos(theta)) + Float64(cos(delta) * sin(phi1))))))))) tmp = 0.0 if (t_4 <= -500.0) tmp = Float64(lambda1 + atan(Float64(sin(theta) * fma(delta, Float64(-0.16666666666666666 * Float64(delta * delta)), delta)), 1.0)); elseif (t_4 <= -5e-7) tmp = t_2; elseif (t_4 <= 2e-12) tmp = Float64(lambda1 + atan(t_3, 1.0)); elseif (t_4 <= 5.0) tmp = t_2; else tmp = Float64(lambda1 + atan(Float64(delta * sin(theta)), 1.0)); end return tmp end
code[lambda1_, phi1_, phi2_, delta_, theta_] := Block[{t$95$1 = N[(N[Sin[delta], $MachinePrecision] * N[Sin[theta], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[ArcTan[t$95$1 / N[Cos[delta], $MachinePrecision]], $MachinePrecision]}, Block[{t$95$3 = N[(N[Cos[phi1], $MachinePrecision] * t$95$1), $MachinePrecision]}, Block[{t$95$4 = N[(lambda1 + N[ArcTan[t$95$3 / N[(N[Cos[delta], $MachinePrecision] - N[(N[Sin[phi1], $MachinePrecision] * N[Sin[N[ArcSin[N[(N[(N[(N[Cos[phi1], $MachinePrecision] * N[Sin[delta], $MachinePrecision]), $MachinePrecision] * N[Cos[theta], $MachinePrecision]), $MachinePrecision] + N[(N[Cos[delta], $MachinePrecision] * N[Sin[phi1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$4, -500.0], N[(lambda1 + N[ArcTan[N[(N[Sin[theta], $MachinePrecision] * N[(delta * N[(-0.16666666666666666 * N[(delta * delta), $MachinePrecision]), $MachinePrecision] + delta), $MachinePrecision]), $MachinePrecision] / 1.0], $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$4, -5e-7], t$95$2, If[LessEqual[t$95$4, 2e-12], N[(lambda1 + N[ArcTan[t$95$3 / 1.0], $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$4, 5.0], t$95$2, N[(lambda1 + N[ArcTan[N[(delta * N[Sin[theta], $MachinePrecision]), $MachinePrecision] / 1.0], $MachinePrecision]), $MachinePrecision]]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \sin delta \cdot \sin theta\\
t_2 := \tan^{-1}_* \frac{t\_1}{\cos delta}\\
t_3 := \cos \phi_1 \cdot t\_1\\
t_4 := \lambda_1 + \tan^{-1}_* \frac{t\_3}{\cos delta - \sin \phi_1 \cdot \sin \sin^{-1} \left(\left(\cos \phi_1 \cdot \sin delta\right) \cdot \cos theta + \cos delta \cdot \sin \phi_1\right)}\\
\mathbf{if}\;t\_4 \leq -500:\\
\;\;\;\;\lambda_1 + \tan^{-1}_* \frac{\sin theta \cdot \mathsf{fma}\left(delta, -0.16666666666666666 \cdot \left(delta \cdot delta\right), delta\right)}{1}\\
\mathbf{elif}\;t\_4 \leq -5 \cdot 10^{-7}:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;t\_4 \leq 2 \cdot 10^{-12}:\\
\;\;\;\;\lambda_1 + \tan^{-1}_* \frac{t\_3}{1}\\
\mathbf{elif}\;t\_4 \leq 5:\\
\;\;\;\;t\_2\\
\mathbf{else}:\\
\;\;\;\;\lambda_1 + \tan^{-1}_* \frac{delta \cdot \sin theta}{1}\\
\end{array}
\end{array}
if (+.f64 lambda1 (atan2.f64 (*.f64 (*.f64 (sin.f64 theta) (sin.f64 delta)) (cos.f64 phi1)) (-.f64 (cos.f64 delta) (*.f64 (sin.f64 phi1) (sin.f64 (asin.f64 (+.f64 (*.f64 (sin.f64 phi1) (cos.f64 delta)) (*.f64 (*.f64 (cos.f64 phi1) (sin.f64 delta)) (cos.f64 theta))))))))) < -500Initial program 100.0%
Taylor expanded in delta around 0
+-commutativeN/A
associate--l+N/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
1-sub-sinN/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
mul-1-negN/A
distribute-lft-neg-outN/A
distribute-rgt-neg-inN/A
accelerator-lowering-fma.f64N/A
Simplified100.0%
Taylor expanded in phi1 around 0
Simplified99.4%
Taylor expanded in phi1 around 0
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
sin-lowering-sin.f6499.4
Simplified99.4%
Taylor expanded in delta around 0
+-commutativeN/A
distribute-lft-inN/A
*-rgt-identityN/A
accelerator-lowering-fma.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6499.4
Simplified99.4%
if -500 < (+.f64 lambda1 (atan2.f64 (*.f64 (*.f64 (sin.f64 theta) (sin.f64 delta)) (cos.f64 phi1)) (-.f64 (cos.f64 delta) (*.f64 (sin.f64 phi1) (sin.f64 (asin.f64 (+.f64 (*.f64 (sin.f64 phi1) (cos.f64 delta)) (*.f64 (*.f64 (cos.f64 phi1) (sin.f64 delta)) (cos.f64 theta))))))))) < -4.99999999999999977e-7 or 1.99999999999999996e-12 < (+.f64 lambda1 (atan2.f64 (*.f64 (*.f64 (sin.f64 theta) (sin.f64 delta)) (cos.f64 phi1)) (-.f64 (cos.f64 delta) (*.f64 (sin.f64 phi1) (sin.f64 (asin.f64 (+.f64 (*.f64 (sin.f64 phi1) (cos.f64 delta)) (*.f64 (*.f64 (cos.f64 phi1) (sin.f64 delta)) (cos.f64 theta))))))))) < 5Initial program 99.5%
Taylor expanded in lambda1 around 0
atan2-lowering-atan2.f64N/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
sin-lowering-sin.f64N/A
sub-negN/A
+-commutativeN/A
Simplified94.4%
Taylor expanded in phi1 around 0
cos-lowering-cos.f6467.5
Simplified67.5%
Taylor expanded in phi1 around 0
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
sin-lowering-sin.f6455.6
Simplified55.6%
if -4.99999999999999977e-7 < (+.f64 lambda1 (atan2.f64 (*.f64 (*.f64 (sin.f64 theta) (sin.f64 delta)) (cos.f64 phi1)) (-.f64 (cos.f64 delta) (*.f64 (sin.f64 phi1) (sin.f64 (asin.f64 (+.f64 (*.f64 (sin.f64 phi1) (cos.f64 delta)) (*.f64 (*.f64 (cos.f64 phi1) (sin.f64 delta)) (cos.f64 theta))))))))) < 1.99999999999999996e-12Initial program 99.5%
Taylor expanded in delta around 0
+-commutativeN/A
associate--l+N/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
1-sub-sinN/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
mul-1-negN/A
distribute-lft-neg-outN/A
distribute-rgt-neg-inN/A
accelerator-lowering-fma.f64N/A
Simplified88.4%
Taylor expanded in phi1 around 0
Simplified80.9%
if 5 < (+.f64 lambda1 (atan2.f64 (*.f64 (*.f64 (sin.f64 theta) (sin.f64 delta)) (cos.f64 phi1)) (-.f64 (cos.f64 delta) (*.f64 (sin.f64 phi1) (sin.f64 (asin.f64 (+.f64 (*.f64 (sin.f64 phi1) (cos.f64 delta)) (*.f64 (*.f64 (cos.f64 phi1) (sin.f64 delta)) (cos.f64 theta))))))))) Initial program 100.0%
Taylor expanded in delta around 0
+-commutativeN/A
associate--l+N/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
1-sub-sinN/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
mul-1-negN/A
distribute-lft-neg-outN/A
distribute-rgt-neg-inN/A
accelerator-lowering-fma.f64N/A
Simplified98.7%
Taylor expanded in phi1 around 0
Simplified98.7%
Taylor expanded in phi1 around 0
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
sin-lowering-sin.f6498.7
Simplified98.7%
Taylor expanded in delta around 0
*-lowering-*.f64N/A
sin-lowering-sin.f6498.8
Simplified98.8%
Final simplification84.7%
(FPCore (lambda1 phi1 phi2 delta theta)
:precision binary64
(let* ((t_1 (* (sin delta) (sin theta)))
(t_2
(+
lambda1
(atan2
(* (cos phi1) t_1)
(-
(cos delta)
(*
(sin phi1)
(sin
(asin
(+
(* (* (cos phi1) (sin delta)) (cos theta))
(* (cos delta) (sin phi1))))))))))
(t_3 (atan2 t_1 (cos delta))))
(if (<= t_2 -500.0)
(+
lambda1
(atan2
(*
(sin theta)
(fma delta (* -0.16666666666666666 (* delta delta)) delta))
1.0))
(if (<= t_2 -5e-7)
t_3
(if (<= t_2 2e-12)
(* lambda1 (- (/ (atan2 t_1 1.0) lambda1) -1.0))
(if (<= t_2 5.0)
t_3
(+ lambda1 (atan2 (* delta (sin theta)) 1.0))))))))
double code(double lambda1, double phi1, double phi2, double delta, double theta) {
double t_1 = sin(delta) * sin(theta);
double t_2 = lambda1 + atan2((cos(phi1) * t_1), (cos(delta) - (sin(phi1) * sin(asin((((cos(phi1) * sin(delta)) * cos(theta)) + (cos(delta) * sin(phi1))))))));
double t_3 = atan2(t_1, cos(delta));
double tmp;
if (t_2 <= -500.0) {
tmp = lambda1 + atan2((sin(theta) * fma(delta, (-0.16666666666666666 * (delta * delta)), delta)), 1.0);
} else if (t_2 <= -5e-7) {
tmp = t_3;
} else if (t_2 <= 2e-12) {
tmp = lambda1 * ((atan2(t_1, 1.0) / lambda1) - -1.0);
} else if (t_2 <= 5.0) {
tmp = t_3;
} else {
tmp = lambda1 + atan2((delta * sin(theta)), 1.0);
}
return tmp;
}
function code(lambda1, phi1, phi2, delta, theta) t_1 = Float64(sin(delta) * sin(theta)) t_2 = Float64(lambda1 + atan(Float64(cos(phi1) * t_1), Float64(cos(delta) - Float64(sin(phi1) * sin(asin(Float64(Float64(Float64(cos(phi1) * sin(delta)) * cos(theta)) + Float64(cos(delta) * sin(phi1))))))))) t_3 = atan(t_1, cos(delta)) tmp = 0.0 if (t_2 <= -500.0) tmp = Float64(lambda1 + atan(Float64(sin(theta) * fma(delta, Float64(-0.16666666666666666 * Float64(delta * delta)), delta)), 1.0)); elseif (t_2 <= -5e-7) tmp = t_3; elseif (t_2 <= 2e-12) tmp = Float64(lambda1 * Float64(Float64(atan(t_1, 1.0) / lambda1) - -1.0)); elseif (t_2 <= 5.0) tmp = t_3; else tmp = Float64(lambda1 + atan(Float64(delta * sin(theta)), 1.0)); end return tmp end
code[lambda1_, phi1_, phi2_, delta_, theta_] := Block[{t$95$1 = N[(N[Sin[delta], $MachinePrecision] * N[Sin[theta], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(lambda1 + N[ArcTan[N[(N[Cos[phi1], $MachinePrecision] * t$95$1), $MachinePrecision] / N[(N[Cos[delta], $MachinePrecision] - N[(N[Sin[phi1], $MachinePrecision] * N[Sin[N[ArcSin[N[(N[(N[(N[Cos[phi1], $MachinePrecision] * N[Sin[delta], $MachinePrecision]), $MachinePrecision] * N[Cos[theta], $MachinePrecision]), $MachinePrecision] + N[(N[Cos[delta], $MachinePrecision] * N[Sin[phi1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = N[ArcTan[t$95$1 / N[Cos[delta], $MachinePrecision]], $MachinePrecision]}, If[LessEqual[t$95$2, -500.0], N[(lambda1 + N[ArcTan[N[(N[Sin[theta], $MachinePrecision] * N[(delta * N[(-0.16666666666666666 * N[(delta * delta), $MachinePrecision]), $MachinePrecision] + delta), $MachinePrecision]), $MachinePrecision] / 1.0], $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$2, -5e-7], t$95$3, If[LessEqual[t$95$2, 2e-12], N[(lambda1 * N[(N[(N[ArcTan[t$95$1 / 1.0], $MachinePrecision] / lambda1), $MachinePrecision] - -1.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$2, 5.0], t$95$3, N[(lambda1 + N[ArcTan[N[(delta * N[Sin[theta], $MachinePrecision]), $MachinePrecision] / 1.0], $MachinePrecision]), $MachinePrecision]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \sin delta \cdot \sin theta\\
t_2 := \lambda_1 + \tan^{-1}_* \frac{\cos \phi_1 \cdot t\_1}{\cos delta - \sin \phi_1 \cdot \sin \sin^{-1} \left(\left(\cos \phi_1 \cdot \sin delta\right) \cdot \cos theta + \cos delta \cdot \sin \phi_1\right)}\\
t_3 := \tan^{-1}_* \frac{t\_1}{\cos delta}\\
\mathbf{if}\;t\_2 \leq -500:\\
\;\;\;\;\lambda_1 + \tan^{-1}_* \frac{\sin theta \cdot \mathsf{fma}\left(delta, -0.16666666666666666 \cdot \left(delta \cdot delta\right), delta\right)}{1}\\
\mathbf{elif}\;t\_2 \leq -5 \cdot 10^{-7}:\\
\;\;\;\;t\_3\\
\mathbf{elif}\;t\_2 \leq 2 \cdot 10^{-12}:\\
\;\;\;\;\lambda_1 \cdot \left(\frac{\tan^{-1}_* \frac{t\_1}{1}}{\lambda_1} - -1\right)\\
\mathbf{elif}\;t\_2 \leq 5:\\
\;\;\;\;t\_3\\
\mathbf{else}:\\
\;\;\;\;\lambda_1 + \tan^{-1}_* \frac{delta \cdot \sin theta}{1}\\
\end{array}
\end{array}
if (+.f64 lambda1 (atan2.f64 (*.f64 (*.f64 (sin.f64 theta) (sin.f64 delta)) (cos.f64 phi1)) (-.f64 (cos.f64 delta) (*.f64 (sin.f64 phi1) (sin.f64 (asin.f64 (+.f64 (*.f64 (sin.f64 phi1) (cos.f64 delta)) (*.f64 (*.f64 (cos.f64 phi1) (sin.f64 delta)) (cos.f64 theta))))))))) < -500Initial program 100.0%
Taylor expanded in delta around 0
+-commutativeN/A
associate--l+N/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
1-sub-sinN/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
mul-1-negN/A
distribute-lft-neg-outN/A
distribute-rgt-neg-inN/A
accelerator-lowering-fma.f64N/A
Simplified100.0%
Taylor expanded in phi1 around 0
Simplified99.4%
Taylor expanded in phi1 around 0
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
sin-lowering-sin.f6499.4
Simplified99.4%
Taylor expanded in delta around 0
+-commutativeN/A
distribute-lft-inN/A
*-rgt-identityN/A
accelerator-lowering-fma.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6499.4
Simplified99.4%
if -500 < (+.f64 lambda1 (atan2.f64 (*.f64 (*.f64 (sin.f64 theta) (sin.f64 delta)) (cos.f64 phi1)) (-.f64 (cos.f64 delta) (*.f64 (sin.f64 phi1) (sin.f64 (asin.f64 (+.f64 (*.f64 (sin.f64 phi1) (cos.f64 delta)) (*.f64 (*.f64 (cos.f64 phi1) (sin.f64 delta)) (cos.f64 theta))))))))) < -4.99999999999999977e-7 or 1.99999999999999996e-12 < (+.f64 lambda1 (atan2.f64 (*.f64 (*.f64 (sin.f64 theta) (sin.f64 delta)) (cos.f64 phi1)) (-.f64 (cos.f64 delta) (*.f64 (sin.f64 phi1) (sin.f64 (asin.f64 (+.f64 (*.f64 (sin.f64 phi1) (cos.f64 delta)) (*.f64 (*.f64 (cos.f64 phi1) (sin.f64 delta)) (cos.f64 theta))))))))) < 5Initial program 99.5%
Taylor expanded in lambda1 around 0
atan2-lowering-atan2.f64N/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
sin-lowering-sin.f64N/A
sub-negN/A
+-commutativeN/A
Simplified94.4%
Taylor expanded in phi1 around 0
cos-lowering-cos.f6467.5
Simplified67.5%
Taylor expanded in phi1 around 0
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
sin-lowering-sin.f6455.6
Simplified55.6%
if -4.99999999999999977e-7 < (+.f64 lambda1 (atan2.f64 (*.f64 (*.f64 (sin.f64 theta) (sin.f64 delta)) (cos.f64 phi1)) (-.f64 (cos.f64 delta) (*.f64 (sin.f64 phi1) (sin.f64 (asin.f64 (+.f64 (*.f64 (sin.f64 phi1) (cos.f64 delta)) (*.f64 (*.f64 (cos.f64 phi1) (sin.f64 delta)) (cos.f64 theta))))))))) < 1.99999999999999996e-12Initial program 99.5%
Taylor expanded in delta around 0
+-commutativeN/A
associate--l+N/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
1-sub-sinN/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
mul-1-negN/A
distribute-lft-neg-outN/A
distribute-rgt-neg-inN/A
accelerator-lowering-fma.f64N/A
Simplified88.4%
Taylor expanded in phi1 around 0
Simplified80.9%
Taylor expanded in phi1 around 0
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
sin-lowering-sin.f6479.6
Simplified79.6%
Taylor expanded in lambda1 around -inf
mul-1-negN/A
*-commutativeN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
*-lowering-*.f64N/A
Simplified79.6%
if 5 < (+.f64 lambda1 (atan2.f64 (*.f64 (*.f64 (sin.f64 theta) (sin.f64 delta)) (cos.f64 phi1)) (-.f64 (cos.f64 delta) (*.f64 (sin.f64 phi1) (sin.f64 (asin.f64 (+.f64 (*.f64 (sin.f64 phi1) (cos.f64 delta)) (*.f64 (*.f64 (cos.f64 phi1) (sin.f64 delta)) (cos.f64 theta))))))))) Initial program 100.0%
Taylor expanded in delta around 0
+-commutativeN/A
associate--l+N/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
1-sub-sinN/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
mul-1-negN/A
distribute-lft-neg-outN/A
distribute-rgt-neg-inN/A
accelerator-lowering-fma.f64N/A
Simplified98.7%
Taylor expanded in phi1 around 0
Simplified98.7%
Taylor expanded in phi1 around 0
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
sin-lowering-sin.f6498.7
Simplified98.7%
Taylor expanded in delta around 0
*-lowering-*.f64N/A
sin-lowering-sin.f6498.8
Simplified98.8%
Final simplification84.4%
(FPCore (lambda1 phi1 phi2 delta theta)
:precision binary64
(let* ((t_1 (* (cos phi1) (* (sin delta) (sin theta))))
(t_2
(atan2
t_1
(-
(cos delta)
(*
(sin phi1)
(sin
(asin
(+
(* (* (cos phi1) (sin delta)) (cos theta))
(* (cos delta) (sin phi1)))))))))
(t_3 (+ lambda1 (atan2 t_1 (cos delta)))))
(if (<= t_2 -0.002)
t_3
(if (<= t_2 2e-19) (+ lambda1 (atan2 t_1 (pow (cos phi1) 2.0))) t_3))))
double code(double lambda1, double phi1, double phi2, double delta, double theta) {
double t_1 = cos(phi1) * (sin(delta) * sin(theta));
double t_2 = atan2(t_1, (cos(delta) - (sin(phi1) * sin(asin((((cos(phi1) * sin(delta)) * cos(theta)) + (cos(delta) * sin(phi1))))))));
double t_3 = lambda1 + atan2(t_1, cos(delta));
double tmp;
if (t_2 <= -0.002) {
tmp = t_3;
} else if (t_2 <= 2e-19) {
tmp = lambda1 + atan2(t_1, pow(cos(phi1), 2.0));
} else {
tmp = t_3;
}
return tmp;
}
real(8) function code(lambda1, phi1, phi2, delta, theta)
real(8), intent (in) :: lambda1
real(8), intent (in) :: phi1
real(8), intent (in) :: phi2
real(8), intent (in) :: delta
real(8), intent (in) :: theta
real(8) :: t_1
real(8) :: t_2
real(8) :: t_3
real(8) :: tmp
t_1 = cos(phi1) * (sin(delta) * sin(theta))
t_2 = atan2(t_1, (cos(delta) - (sin(phi1) * sin(asin((((cos(phi1) * sin(delta)) * cos(theta)) + (cos(delta) * sin(phi1))))))))
t_3 = lambda1 + atan2(t_1, cos(delta))
if (t_2 <= (-0.002d0)) then
tmp = t_3
else if (t_2 <= 2d-19) then
tmp = lambda1 + atan2(t_1, (cos(phi1) ** 2.0d0))
else
tmp = t_3
end if
code = tmp
end function
public static double code(double lambda1, double phi1, double phi2, double delta, double theta) {
double t_1 = Math.cos(phi1) * (Math.sin(delta) * Math.sin(theta));
double t_2 = Math.atan2(t_1, (Math.cos(delta) - (Math.sin(phi1) * Math.sin(Math.asin((((Math.cos(phi1) * Math.sin(delta)) * Math.cos(theta)) + (Math.cos(delta) * Math.sin(phi1))))))));
double t_3 = lambda1 + Math.atan2(t_1, Math.cos(delta));
double tmp;
if (t_2 <= -0.002) {
tmp = t_3;
} else if (t_2 <= 2e-19) {
tmp = lambda1 + Math.atan2(t_1, Math.pow(Math.cos(phi1), 2.0));
} else {
tmp = t_3;
}
return tmp;
}
def code(lambda1, phi1, phi2, delta, theta): t_1 = math.cos(phi1) * (math.sin(delta) * math.sin(theta)) t_2 = math.atan2(t_1, (math.cos(delta) - (math.sin(phi1) * math.sin(math.asin((((math.cos(phi1) * math.sin(delta)) * math.cos(theta)) + (math.cos(delta) * math.sin(phi1)))))))) t_3 = lambda1 + math.atan2(t_1, math.cos(delta)) tmp = 0 if t_2 <= -0.002: tmp = t_3 elif t_2 <= 2e-19: tmp = lambda1 + math.atan2(t_1, math.pow(math.cos(phi1), 2.0)) else: tmp = t_3 return tmp
function code(lambda1, phi1, phi2, delta, theta) t_1 = Float64(cos(phi1) * Float64(sin(delta) * sin(theta))) t_2 = atan(t_1, Float64(cos(delta) - Float64(sin(phi1) * sin(asin(Float64(Float64(Float64(cos(phi1) * sin(delta)) * cos(theta)) + Float64(cos(delta) * sin(phi1)))))))) t_3 = Float64(lambda1 + atan(t_1, cos(delta))) tmp = 0.0 if (t_2 <= -0.002) tmp = t_3; elseif (t_2 <= 2e-19) tmp = Float64(lambda1 + atan(t_1, (cos(phi1) ^ 2.0))); else tmp = t_3; end return tmp end
function tmp_2 = code(lambda1, phi1, phi2, delta, theta) t_1 = cos(phi1) * (sin(delta) * sin(theta)); t_2 = atan2(t_1, (cos(delta) - (sin(phi1) * sin(asin((((cos(phi1) * sin(delta)) * cos(theta)) + (cos(delta) * sin(phi1)))))))); t_3 = lambda1 + atan2(t_1, cos(delta)); tmp = 0.0; if (t_2 <= -0.002) tmp = t_3; elseif (t_2 <= 2e-19) tmp = lambda1 + atan2(t_1, (cos(phi1) ^ 2.0)); else tmp = t_3; end tmp_2 = tmp; end
code[lambda1_, phi1_, phi2_, delta_, theta_] := Block[{t$95$1 = N[(N[Cos[phi1], $MachinePrecision] * N[(N[Sin[delta], $MachinePrecision] * N[Sin[theta], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[ArcTan[t$95$1 / N[(N[Cos[delta], $MachinePrecision] - N[(N[Sin[phi1], $MachinePrecision] * N[Sin[N[ArcSin[N[(N[(N[(N[Cos[phi1], $MachinePrecision] * N[Sin[delta], $MachinePrecision]), $MachinePrecision] * N[Cos[theta], $MachinePrecision]), $MachinePrecision] + N[(N[Cos[delta], $MachinePrecision] * N[Sin[phi1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$3 = N[(lambda1 + N[ArcTan[t$95$1 / N[Cos[delta], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$2, -0.002], t$95$3, If[LessEqual[t$95$2, 2e-19], N[(lambda1 + N[ArcTan[t$95$1 / N[Power[N[Cos[phi1], $MachinePrecision], 2.0], $MachinePrecision]], $MachinePrecision]), $MachinePrecision], t$95$3]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \cos \phi_1 \cdot \left(\sin delta \cdot \sin theta\right)\\
t_2 := \tan^{-1}_* \frac{t\_1}{\cos delta - \sin \phi_1 \cdot \sin \sin^{-1} \left(\left(\cos \phi_1 \cdot \sin delta\right) \cdot \cos theta + \cos delta \cdot \sin \phi_1\right)}\\
t_3 := \lambda_1 + \tan^{-1}_* \frac{t\_1}{\cos delta}\\
\mathbf{if}\;t\_2 \leq -0.002:\\
\;\;\;\;t\_3\\
\mathbf{elif}\;t\_2 \leq 2 \cdot 10^{-19}:\\
\;\;\;\;\lambda_1 + \tan^{-1}_* \frac{t\_1}{{\cos \phi_1}^{2}}\\
\mathbf{else}:\\
\;\;\;\;t\_3\\
\end{array}
\end{array}
if (atan2.f64 (*.f64 (*.f64 (sin.f64 theta) (sin.f64 delta)) (cos.f64 phi1)) (-.f64 (cos.f64 delta) (*.f64 (sin.f64 phi1) (sin.f64 (asin.f64 (+.f64 (*.f64 (sin.f64 phi1) (cos.f64 delta)) (*.f64 (*.f64 (cos.f64 phi1) (sin.f64 delta)) (cos.f64 theta)))))))) < -2e-3 or 2e-19 < (atan2.f64 (*.f64 (*.f64 (sin.f64 theta) (sin.f64 delta)) (cos.f64 phi1)) (-.f64 (cos.f64 delta) (*.f64 (sin.f64 phi1) (sin.f64 (asin.f64 (+.f64 (*.f64 (sin.f64 phi1) (cos.f64 delta)) (*.f64 (*.f64 (cos.f64 phi1) (sin.f64 delta)) (cos.f64 theta)))))))) Initial program 99.7%
Taylor expanded in phi1 around 0
cos-lowering-cos.f6485.3
Simplified85.3%
if -2e-3 < (atan2.f64 (*.f64 (*.f64 (sin.f64 theta) (sin.f64 delta)) (cos.f64 phi1)) (-.f64 (cos.f64 delta) (*.f64 (sin.f64 phi1) (sin.f64 (asin.f64 (+.f64 (*.f64 (sin.f64 phi1) (cos.f64 delta)) (*.f64 (*.f64 (cos.f64 phi1) (sin.f64 delta)) (cos.f64 theta)))))))) < 2e-19Initial program 99.8%
sub-negN/A
+-commutativeN/A
sin-asinN/A
distribute-lft-neg-inN/A
distribute-rgt-inN/A
associate-+l+N/A
Applied egg-rr99.8%
Taylor expanded in delta around 0
mul-1-negN/A
sub-negN/A
unpow2N/A
1-sub-sinN/A
unpow2N/A
pow-lowering-pow.f64N/A
cos-lowering-cos.f6496.5
Simplified96.5%
Final simplification91.6%
(FPCore (lambda1 phi1 phi2 delta theta)
:precision binary64
(let* ((t_1 (* (cos phi1) (* (sin delta) (sin theta))))
(t_2
(atan2
t_1
(-
(cos delta)
(*
(sin phi1)
(sin
(asin
(+
(* (* (cos phi1) (sin delta)) (cos theta))
(* (cos delta) (sin phi1)))))))))
(t_3 (+ lambda1 (atan2 t_1 (cos delta)))))
(if (<= t_2 -0.002)
t_3
(if (<= t_2 2e-19)
(+ lambda1 (atan2 t_1 (fma 0.5 (cos (* phi1 2.0)) 0.5)))
t_3))))
double code(double lambda1, double phi1, double phi2, double delta, double theta) {
double t_1 = cos(phi1) * (sin(delta) * sin(theta));
double t_2 = atan2(t_1, (cos(delta) - (sin(phi1) * sin(asin((((cos(phi1) * sin(delta)) * cos(theta)) + (cos(delta) * sin(phi1))))))));
double t_3 = lambda1 + atan2(t_1, cos(delta));
double tmp;
if (t_2 <= -0.002) {
tmp = t_3;
} else if (t_2 <= 2e-19) {
tmp = lambda1 + atan2(t_1, fma(0.5, cos((phi1 * 2.0)), 0.5));
} else {
tmp = t_3;
}
return tmp;
}
function code(lambda1, phi1, phi2, delta, theta) t_1 = Float64(cos(phi1) * Float64(sin(delta) * sin(theta))) t_2 = atan(t_1, Float64(cos(delta) - Float64(sin(phi1) * sin(asin(Float64(Float64(Float64(cos(phi1) * sin(delta)) * cos(theta)) + Float64(cos(delta) * sin(phi1)))))))) t_3 = Float64(lambda1 + atan(t_1, cos(delta))) tmp = 0.0 if (t_2 <= -0.002) tmp = t_3; elseif (t_2 <= 2e-19) tmp = Float64(lambda1 + atan(t_1, fma(0.5, cos(Float64(phi1 * 2.0)), 0.5))); else tmp = t_3; end return tmp end
code[lambda1_, phi1_, phi2_, delta_, theta_] := Block[{t$95$1 = N[(N[Cos[phi1], $MachinePrecision] * N[(N[Sin[delta], $MachinePrecision] * N[Sin[theta], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[ArcTan[t$95$1 / N[(N[Cos[delta], $MachinePrecision] - N[(N[Sin[phi1], $MachinePrecision] * N[Sin[N[ArcSin[N[(N[(N[(N[Cos[phi1], $MachinePrecision] * N[Sin[delta], $MachinePrecision]), $MachinePrecision] * N[Cos[theta], $MachinePrecision]), $MachinePrecision] + N[(N[Cos[delta], $MachinePrecision] * N[Sin[phi1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$3 = N[(lambda1 + N[ArcTan[t$95$1 / N[Cos[delta], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$2, -0.002], t$95$3, If[LessEqual[t$95$2, 2e-19], N[(lambda1 + N[ArcTan[t$95$1 / N[(0.5 * N[Cos[N[(phi1 * 2.0), $MachinePrecision]], $MachinePrecision] + 0.5), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], t$95$3]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \cos \phi_1 \cdot \left(\sin delta \cdot \sin theta\right)\\
t_2 := \tan^{-1}_* \frac{t\_1}{\cos delta - \sin \phi_1 \cdot \sin \sin^{-1} \left(\left(\cos \phi_1 \cdot \sin delta\right) \cdot \cos theta + \cos delta \cdot \sin \phi_1\right)}\\
t_3 := \lambda_1 + \tan^{-1}_* \frac{t\_1}{\cos delta}\\
\mathbf{if}\;t\_2 \leq -0.002:\\
\;\;\;\;t\_3\\
\mathbf{elif}\;t\_2 \leq 2 \cdot 10^{-19}:\\
\;\;\;\;\lambda_1 + \tan^{-1}_* \frac{t\_1}{\mathsf{fma}\left(0.5, \cos \left(\phi_1 \cdot 2\right), 0.5\right)}\\
\mathbf{else}:\\
\;\;\;\;t\_3\\
\end{array}
\end{array}
if (atan2.f64 (*.f64 (*.f64 (sin.f64 theta) (sin.f64 delta)) (cos.f64 phi1)) (-.f64 (cos.f64 delta) (*.f64 (sin.f64 phi1) (sin.f64 (asin.f64 (+.f64 (*.f64 (sin.f64 phi1) (cos.f64 delta)) (*.f64 (*.f64 (cos.f64 phi1) (sin.f64 delta)) (cos.f64 theta)))))))) < -2e-3 or 2e-19 < (atan2.f64 (*.f64 (*.f64 (sin.f64 theta) (sin.f64 delta)) (cos.f64 phi1)) (-.f64 (cos.f64 delta) (*.f64 (sin.f64 phi1) (sin.f64 (asin.f64 (+.f64 (*.f64 (sin.f64 phi1) (cos.f64 delta)) (*.f64 (*.f64 (cos.f64 phi1) (sin.f64 delta)) (cos.f64 theta)))))))) Initial program 99.7%
Taylor expanded in phi1 around 0
cos-lowering-cos.f6485.3
Simplified85.3%
if -2e-3 < (atan2.f64 (*.f64 (*.f64 (sin.f64 theta) (sin.f64 delta)) (cos.f64 phi1)) (-.f64 (cos.f64 delta) (*.f64 (sin.f64 phi1) (sin.f64 (asin.f64 (+.f64 (*.f64 (sin.f64 phi1) (cos.f64 delta)) (*.f64 (*.f64 (cos.f64 phi1) (sin.f64 delta)) (cos.f64 theta)))))))) < 2e-19Initial program 99.8%
sin-asinN/A
distribute-rgt-inN/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
accelerator-lowering-fma.f64N/A
sqr-sin-aN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
*-lowering-*.f64N/A
Applied egg-rr99.8%
Taylor expanded in delta around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
cos-lowering-cos.f64N/A
*-lowering-*.f6496.4
Simplified96.4%
Final simplification91.5%
(FPCore (lambda1 phi1 phi2 delta theta)
:precision binary64
(+
lambda1
(atan2
(* (cos phi1) (* (sin delta) (sin theta)))
(fma
(- (sin phi1))
(fma (cos phi1) (sin delta) (* (cos delta) (sin phi1)))
(cos delta)))))
double code(double lambda1, double phi1, double phi2, double delta, double theta) {
return lambda1 + atan2((cos(phi1) * (sin(delta) * sin(theta))), fma(-sin(phi1), fma(cos(phi1), sin(delta), (cos(delta) * sin(phi1))), cos(delta)));
}
function code(lambda1, phi1, phi2, delta, theta) return Float64(lambda1 + atan(Float64(cos(phi1) * Float64(sin(delta) * sin(theta))), fma(Float64(-sin(phi1)), fma(cos(phi1), sin(delta), Float64(cos(delta) * sin(phi1))), cos(delta)))) end
code[lambda1_, phi1_, phi2_, delta_, theta_] := N[(lambda1 + N[ArcTan[N[(N[Cos[phi1], $MachinePrecision] * N[(N[Sin[delta], $MachinePrecision] * N[Sin[theta], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[((-N[Sin[phi1], $MachinePrecision]) * N[(N[Cos[phi1], $MachinePrecision] * N[Sin[delta], $MachinePrecision] + N[(N[Cos[delta], $MachinePrecision] * N[Sin[phi1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[Cos[delta], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\lambda_1 + \tan^{-1}_* \frac{\cos \phi_1 \cdot \left(\sin delta \cdot \sin theta\right)}{\mathsf{fma}\left(-\sin \phi_1, \mathsf{fma}\left(\cos \phi_1, \sin delta, \cos delta \cdot \sin \phi_1\right), \cos delta\right)}
\end{array}
Initial program 99.8%
sub-negN/A
+-commutativeN/A
sin-asinN/A
distribute-lft-neg-inN/A
distribute-rgt-inN/A
associate-+l+N/A
Applied egg-rr99.8%
Taylor expanded in theta around 0
+-commutativeN/A
+-commutativeN/A
Simplified94.1%
Final simplification94.1%
(FPCore (lambda1 phi1 phi2 delta theta)
:precision binary64
(+
lambda1
(atan2
(* (cos phi1) (* (sin delta) (sin theta)))
(fma
(+ 0.5 (* 0.5 (cos (* phi1 -2.0))))
(cos delta)
(* (sin phi1) (* (cos phi1) (- (sin delta))))))))
double code(double lambda1, double phi1, double phi2, double delta, double theta) {
return lambda1 + atan2((cos(phi1) * (sin(delta) * sin(theta))), fma((0.5 + (0.5 * cos((phi1 * -2.0)))), cos(delta), (sin(phi1) * (cos(phi1) * -sin(delta)))));
}
function code(lambda1, phi1, phi2, delta, theta) return Float64(lambda1 + atan(Float64(cos(phi1) * Float64(sin(delta) * sin(theta))), fma(Float64(0.5 + Float64(0.5 * cos(Float64(phi1 * -2.0)))), cos(delta), Float64(sin(phi1) * Float64(cos(phi1) * Float64(-sin(delta))))))) end
code[lambda1_, phi1_, phi2_, delta_, theta_] := N[(lambda1 + N[ArcTan[N[(N[Cos[phi1], $MachinePrecision] * N[(N[Sin[delta], $MachinePrecision] * N[Sin[theta], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(N[(0.5 + N[(0.5 * N[Cos[N[(phi1 * -2.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[Cos[delta], $MachinePrecision] + N[(N[Sin[phi1], $MachinePrecision] * N[(N[Cos[phi1], $MachinePrecision] * (-N[Sin[delta], $MachinePrecision])), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\lambda_1 + \tan^{-1}_* \frac{\cos \phi_1 \cdot \left(\sin delta \cdot \sin theta\right)}{\mathsf{fma}\left(0.5 + 0.5 \cdot \cos \left(\phi_1 \cdot -2\right), \cos delta, \sin \phi_1 \cdot \left(\cos \phi_1 \cdot \left(-\sin delta\right)\right)\right)}
\end{array}
Initial program 99.8%
sub-negN/A
+-commutativeN/A
sin-asinN/A
distribute-lft-neg-inN/A
distribute-rgt-inN/A
associate-+l+N/A
Applied egg-rr99.8%
distribute-lft-neg-outN/A
neg-lowering-neg.f64N/A
sqr-sin-aN/A
sub-negN/A
+-commutativeN/A
*-commutativeN/A
distribute-rgt-neg-inN/A
accelerator-lowering-fma.f64N/A
cos-2N/A
cos-sumN/A
cos-lowering-cos.f64N/A
+-lowering-+.f64N/A
metadata-eval99.8
Applied egg-rr99.8%
Taylor expanded in theta around 0
associate-+r+N/A
mul-1-negN/A
unsub-negN/A
*-commutativeN/A
cancel-sign-sub-invN/A
mul-1-negN/A
distribute-rgt1-inN/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
Simplified94.0%
Final simplification94.0%
(FPCore (lambda1 phi1 phi2 delta theta) :precision binary64 (+ lambda1 (atan2 (* (cos phi1) (* (sin delta) (sin theta))) (+ (cos delta) (- (* 0.5 (cos (+ phi1 phi1))) 0.5)))))
double code(double lambda1, double phi1, double phi2, double delta, double theta) {
return lambda1 + atan2((cos(phi1) * (sin(delta) * sin(theta))), (cos(delta) + ((0.5 * cos((phi1 + phi1))) - 0.5)));
}
real(8) function code(lambda1, phi1, phi2, delta, theta)
real(8), intent (in) :: lambda1
real(8), intent (in) :: phi1
real(8), intent (in) :: phi2
real(8), intent (in) :: delta
real(8), intent (in) :: theta
code = lambda1 + atan2((cos(phi1) * (sin(delta) * sin(theta))), (cos(delta) + ((0.5d0 * cos((phi1 + phi1))) - 0.5d0)))
end function
public static double code(double lambda1, double phi1, double phi2, double delta, double theta) {
return lambda1 + Math.atan2((Math.cos(phi1) * (Math.sin(delta) * Math.sin(theta))), (Math.cos(delta) + ((0.5 * Math.cos((phi1 + phi1))) - 0.5)));
}
def code(lambda1, phi1, phi2, delta, theta): return lambda1 + math.atan2((math.cos(phi1) * (math.sin(delta) * math.sin(theta))), (math.cos(delta) + ((0.5 * math.cos((phi1 + phi1))) - 0.5)))
function code(lambda1, phi1, phi2, delta, theta) return Float64(lambda1 + atan(Float64(cos(phi1) * Float64(sin(delta) * sin(theta))), Float64(cos(delta) + Float64(Float64(0.5 * cos(Float64(phi1 + phi1))) - 0.5)))) end
function tmp = code(lambda1, phi1, phi2, delta, theta) tmp = lambda1 + atan2((cos(phi1) * (sin(delta) * sin(theta))), (cos(delta) + ((0.5 * cos((phi1 + phi1))) - 0.5))); end
code[lambda1_, phi1_, phi2_, delta_, theta_] := N[(lambda1 + N[ArcTan[N[(N[Cos[phi1], $MachinePrecision] * N[(N[Sin[delta], $MachinePrecision] * N[Sin[theta], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(N[Cos[delta], $MachinePrecision] + N[(N[(0.5 * N[Cos[N[(phi1 + phi1), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] - 0.5), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\lambda_1 + \tan^{-1}_* \frac{\cos \phi_1 \cdot \left(\sin delta \cdot \sin theta\right)}{\cos delta + \left(0.5 \cdot \cos \left(\phi_1 + \phi_1\right) - 0.5\right)}
\end{array}
Initial program 99.8%
Taylor expanded in delta around 0
pow-lowering-pow.f64N/A
sin-lowering-sin.f6490.0
Simplified90.0%
+-commutativeN/A
+-lowering-+.f64N/A
Applied egg-rr90.0%
Final simplification90.0%
(FPCore (lambda1 phi1 phi2 delta theta)
:precision binary64
(let* ((t_1
(+
lambda1
(atan2 (* (cos phi1) (* (sin delta) (sin theta))) (cos delta)))))
(if (<= delta -4.5e-6)
t_1
(if (<= delta 8.6e-18)
(+
lambda1
(atan2 (* (cos phi1) (* delta (sin theta))) (pow (cos phi1) 2.0)))
t_1))))
double code(double lambda1, double phi1, double phi2, double delta, double theta) {
double t_1 = lambda1 + atan2((cos(phi1) * (sin(delta) * sin(theta))), cos(delta));
double tmp;
if (delta <= -4.5e-6) {
tmp = t_1;
} else if (delta <= 8.6e-18) {
tmp = lambda1 + atan2((cos(phi1) * (delta * sin(theta))), pow(cos(phi1), 2.0));
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(lambda1, phi1, phi2, delta, theta)
real(8), intent (in) :: lambda1
real(8), intent (in) :: phi1
real(8), intent (in) :: phi2
real(8), intent (in) :: delta
real(8), intent (in) :: theta
real(8) :: t_1
real(8) :: tmp
t_1 = lambda1 + atan2((cos(phi1) * (sin(delta) * sin(theta))), cos(delta))
if (delta <= (-4.5d-6)) then
tmp = t_1
else if (delta <= 8.6d-18) then
tmp = lambda1 + atan2((cos(phi1) * (delta * sin(theta))), (cos(phi1) ** 2.0d0))
else
tmp = t_1
end if
code = tmp
end function
public static double code(double lambda1, double phi1, double phi2, double delta, double theta) {
double t_1 = lambda1 + Math.atan2((Math.cos(phi1) * (Math.sin(delta) * Math.sin(theta))), Math.cos(delta));
double tmp;
if (delta <= -4.5e-6) {
tmp = t_1;
} else if (delta <= 8.6e-18) {
tmp = lambda1 + Math.atan2((Math.cos(phi1) * (delta * Math.sin(theta))), Math.pow(Math.cos(phi1), 2.0));
} else {
tmp = t_1;
}
return tmp;
}
def code(lambda1, phi1, phi2, delta, theta): t_1 = lambda1 + math.atan2((math.cos(phi1) * (math.sin(delta) * math.sin(theta))), math.cos(delta)) tmp = 0 if delta <= -4.5e-6: tmp = t_1 elif delta <= 8.6e-18: tmp = lambda1 + math.atan2((math.cos(phi1) * (delta * math.sin(theta))), math.pow(math.cos(phi1), 2.0)) else: tmp = t_1 return tmp
function code(lambda1, phi1, phi2, delta, theta) t_1 = Float64(lambda1 + atan(Float64(cos(phi1) * Float64(sin(delta) * sin(theta))), cos(delta))) tmp = 0.0 if (delta <= -4.5e-6) tmp = t_1; elseif (delta <= 8.6e-18) tmp = Float64(lambda1 + atan(Float64(cos(phi1) * Float64(delta * sin(theta))), (cos(phi1) ^ 2.0))); else tmp = t_1; end return tmp end
function tmp_2 = code(lambda1, phi1, phi2, delta, theta) t_1 = lambda1 + atan2((cos(phi1) * (sin(delta) * sin(theta))), cos(delta)); tmp = 0.0; if (delta <= -4.5e-6) tmp = t_1; elseif (delta <= 8.6e-18) tmp = lambda1 + atan2((cos(phi1) * (delta * sin(theta))), (cos(phi1) ^ 2.0)); else tmp = t_1; end tmp_2 = tmp; end
code[lambda1_, phi1_, phi2_, delta_, theta_] := Block[{t$95$1 = N[(lambda1 + N[ArcTan[N[(N[Cos[phi1], $MachinePrecision] * N[(N[Sin[delta], $MachinePrecision] * N[Sin[theta], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[Cos[delta], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[delta, -4.5e-6], t$95$1, If[LessEqual[delta, 8.6e-18], N[(lambda1 + N[ArcTan[N[(N[Cos[phi1], $MachinePrecision] * N[(delta * N[Sin[theta], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[Power[N[Cos[phi1], $MachinePrecision], 2.0], $MachinePrecision]], $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \lambda_1 + \tan^{-1}_* \frac{\cos \phi_1 \cdot \left(\sin delta \cdot \sin theta\right)}{\cos delta}\\
\mathbf{if}\;delta \leq -4.5 \cdot 10^{-6}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;delta \leq 8.6 \cdot 10^{-18}:\\
\;\;\;\;\lambda_1 + \tan^{-1}_* \frac{\cos \phi_1 \cdot \left(delta \cdot \sin theta\right)}{{\cos \phi_1}^{2}}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if delta < -4.50000000000000011e-6 or 8.6000000000000005e-18 < delta Initial program 99.8%
Taylor expanded in phi1 around 0
cos-lowering-cos.f6482.9
Simplified82.9%
if -4.50000000000000011e-6 < delta < 8.6000000000000005e-18Initial program 99.7%
Taylor expanded in delta around 0
+-commutativeN/A
associate--l+N/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
1-sub-sinN/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
mul-1-negN/A
distribute-lft-neg-outN/A
distribute-rgt-neg-inN/A
accelerator-lowering-fma.f64N/A
Simplified99.8%
Taylor expanded in delta around 0
*-lowering-*.f64N/A
sin-lowering-sin.f6499.8
Simplified99.8%
Taylor expanded in delta around 0
pow-lowering-pow.f64N/A
cos-lowering-cos.f6499.5
Simplified99.5%
Final simplification90.4%
(FPCore (lambda1 phi1 phi2 delta theta) :precision binary64 (+ lambda1 (atan2 (* (cos phi1) (* (sin delta) (sin theta))) (cos delta))))
double code(double lambda1, double phi1, double phi2, double delta, double theta) {
return lambda1 + atan2((cos(phi1) * (sin(delta) * sin(theta))), cos(delta));
}
real(8) function code(lambda1, phi1, phi2, delta, theta)
real(8), intent (in) :: lambda1
real(8), intent (in) :: phi1
real(8), intent (in) :: phi2
real(8), intent (in) :: delta
real(8), intent (in) :: theta
code = lambda1 + atan2((cos(phi1) * (sin(delta) * sin(theta))), cos(delta))
end function
public static double code(double lambda1, double phi1, double phi2, double delta, double theta) {
return lambda1 + Math.atan2((Math.cos(phi1) * (Math.sin(delta) * Math.sin(theta))), Math.cos(delta));
}
def code(lambda1, phi1, phi2, delta, theta): return lambda1 + math.atan2((math.cos(phi1) * (math.sin(delta) * math.sin(theta))), math.cos(delta))
function code(lambda1, phi1, phi2, delta, theta) return Float64(lambda1 + atan(Float64(cos(phi1) * Float64(sin(delta) * sin(theta))), cos(delta))) end
function tmp = code(lambda1, phi1, phi2, delta, theta) tmp = lambda1 + atan2((cos(phi1) * (sin(delta) * sin(theta))), cos(delta)); end
code[lambda1_, phi1_, phi2_, delta_, theta_] := N[(lambda1 + N[ArcTan[N[(N[Cos[phi1], $MachinePrecision] * N[(N[Sin[delta], $MachinePrecision] * N[Sin[theta], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[Cos[delta], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\lambda_1 + \tan^{-1}_* \frac{\cos \phi_1 \cdot \left(\sin delta \cdot \sin theta\right)}{\cos delta}
\end{array}
Initial program 99.8%
Taylor expanded in phi1 around 0
cos-lowering-cos.f6487.3
Simplified87.3%
Final simplification87.3%
(FPCore (lambda1 phi1 phi2 delta theta) :precision binary64 (* lambda1 (- (/ (atan2 (* (sin delta) (sin theta)) 1.0) lambda1) -1.0)))
double code(double lambda1, double phi1, double phi2, double delta, double theta) {
return lambda1 * ((atan2((sin(delta) * sin(theta)), 1.0) / lambda1) - -1.0);
}
real(8) function code(lambda1, phi1, phi2, delta, theta)
real(8), intent (in) :: lambda1
real(8), intent (in) :: phi1
real(8), intent (in) :: phi2
real(8), intent (in) :: delta
real(8), intent (in) :: theta
code = lambda1 * ((atan2((sin(delta) * sin(theta)), 1.0d0) / lambda1) - (-1.0d0))
end function
public static double code(double lambda1, double phi1, double phi2, double delta, double theta) {
return lambda1 * ((Math.atan2((Math.sin(delta) * Math.sin(theta)), 1.0) / lambda1) - -1.0);
}
def code(lambda1, phi1, phi2, delta, theta): return lambda1 * ((math.atan2((math.sin(delta) * math.sin(theta)), 1.0) / lambda1) - -1.0)
function code(lambda1, phi1, phi2, delta, theta) return Float64(lambda1 * Float64(Float64(atan(Float64(sin(delta) * sin(theta)), 1.0) / lambda1) - -1.0)) end
function tmp = code(lambda1, phi1, phi2, delta, theta) tmp = lambda1 * ((atan2((sin(delta) * sin(theta)), 1.0) / lambda1) - -1.0); end
code[lambda1_, phi1_, phi2_, delta_, theta_] := N[(lambda1 * N[(N[(N[ArcTan[N[(N[Sin[delta], $MachinePrecision] * N[Sin[theta], $MachinePrecision]), $MachinePrecision] / 1.0], $MachinePrecision] / lambda1), $MachinePrecision] - -1.0), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\lambda_1 \cdot \left(\frac{\tan^{-1}_* \frac{\sin delta \cdot \sin theta}{1}}{\lambda_1} - -1\right)
\end{array}
Initial program 99.8%
Taylor expanded in delta around 0
+-commutativeN/A
associate--l+N/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
1-sub-sinN/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
mul-1-negN/A
distribute-lft-neg-outN/A
distribute-rgt-neg-inN/A
accelerator-lowering-fma.f64N/A
Simplified79.7%
Taylor expanded in phi1 around 0
Simplified76.1%
Taylor expanded in phi1 around 0
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
sin-lowering-sin.f6475.1
Simplified75.1%
Taylor expanded in lambda1 around -inf
mul-1-negN/A
*-commutativeN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
*-lowering-*.f64N/A
Simplified75.1%
Final simplification75.1%
(FPCore (lambda1 phi1 phi2 delta theta) :precision binary64 (+ lambda1 (atan2 (* (sin delta) (sin theta)) 1.0)))
double code(double lambda1, double phi1, double phi2, double delta, double theta) {
return lambda1 + atan2((sin(delta) * sin(theta)), 1.0);
}
real(8) function code(lambda1, phi1, phi2, delta, theta)
real(8), intent (in) :: lambda1
real(8), intent (in) :: phi1
real(8), intent (in) :: phi2
real(8), intent (in) :: delta
real(8), intent (in) :: theta
code = lambda1 + atan2((sin(delta) * sin(theta)), 1.0d0)
end function
public static double code(double lambda1, double phi1, double phi2, double delta, double theta) {
return lambda1 + Math.atan2((Math.sin(delta) * Math.sin(theta)), 1.0);
}
def code(lambda1, phi1, phi2, delta, theta): return lambda1 + math.atan2((math.sin(delta) * math.sin(theta)), 1.0)
function code(lambda1, phi1, phi2, delta, theta) return Float64(lambda1 + atan(Float64(sin(delta) * sin(theta)), 1.0)) end
function tmp = code(lambda1, phi1, phi2, delta, theta) tmp = lambda1 + atan2((sin(delta) * sin(theta)), 1.0); end
code[lambda1_, phi1_, phi2_, delta_, theta_] := N[(lambda1 + N[ArcTan[N[(N[Sin[delta], $MachinePrecision] * N[Sin[theta], $MachinePrecision]), $MachinePrecision] / 1.0], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\lambda_1 + \tan^{-1}_* \frac{\sin delta \cdot \sin theta}{1}
\end{array}
Initial program 99.8%
Taylor expanded in delta around 0
+-commutativeN/A
associate--l+N/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
1-sub-sinN/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
mul-1-negN/A
distribute-lft-neg-outN/A
distribute-rgt-neg-inN/A
accelerator-lowering-fma.f64N/A
Simplified79.7%
Taylor expanded in phi1 around 0
Simplified76.1%
Taylor expanded in phi1 around 0
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
sin-lowering-sin.f6475.1
Simplified75.1%
(FPCore (lambda1 phi1 phi2 delta theta) :precision binary64 (+ lambda1 (atan2 (* delta (sin theta)) 1.0)))
double code(double lambda1, double phi1, double phi2, double delta, double theta) {
return lambda1 + atan2((delta * sin(theta)), 1.0);
}
real(8) function code(lambda1, phi1, phi2, delta, theta)
real(8), intent (in) :: lambda1
real(8), intent (in) :: phi1
real(8), intent (in) :: phi2
real(8), intent (in) :: delta
real(8), intent (in) :: theta
code = lambda1 + atan2((delta * sin(theta)), 1.0d0)
end function
public static double code(double lambda1, double phi1, double phi2, double delta, double theta) {
return lambda1 + Math.atan2((delta * Math.sin(theta)), 1.0);
}
def code(lambda1, phi1, phi2, delta, theta): return lambda1 + math.atan2((delta * math.sin(theta)), 1.0)
function code(lambda1, phi1, phi2, delta, theta) return Float64(lambda1 + atan(Float64(delta * sin(theta)), 1.0)) end
function tmp = code(lambda1, phi1, phi2, delta, theta) tmp = lambda1 + atan2((delta * sin(theta)), 1.0); end
code[lambda1_, phi1_, phi2_, delta_, theta_] := N[(lambda1 + N[ArcTan[N[(delta * N[Sin[theta], $MachinePrecision]), $MachinePrecision] / 1.0], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\lambda_1 + \tan^{-1}_* \frac{delta \cdot \sin theta}{1}
\end{array}
Initial program 99.8%
Taylor expanded in delta around 0
+-commutativeN/A
associate--l+N/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
1-sub-sinN/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
mul-1-negN/A
distribute-lft-neg-outN/A
distribute-rgt-neg-inN/A
accelerator-lowering-fma.f64N/A
Simplified79.7%
Taylor expanded in phi1 around 0
Simplified76.1%
Taylor expanded in phi1 around 0
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
sin-lowering-sin.f6475.1
Simplified75.1%
Taylor expanded in delta around 0
*-lowering-*.f64N/A
sin-lowering-sin.f6472.9
Simplified72.9%
(FPCore (lambda1 phi1 phi2 delta theta) :precision binary64 lambda1)
double code(double lambda1, double phi1, double phi2, double delta, double theta) {
return lambda1;
}
real(8) function code(lambda1, phi1, phi2, delta, theta)
real(8), intent (in) :: lambda1
real(8), intent (in) :: phi1
real(8), intent (in) :: phi2
real(8), intent (in) :: delta
real(8), intent (in) :: theta
code = lambda1
end function
public static double code(double lambda1, double phi1, double phi2, double delta, double theta) {
return lambda1;
}
def code(lambda1, phi1, phi2, delta, theta): return lambda1
function code(lambda1, phi1, phi2, delta, theta) return lambda1 end
function tmp = code(lambda1, phi1, phi2, delta, theta) tmp = lambda1; end
code[lambda1_, phi1_, phi2_, delta_, theta_] := lambda1
\begin{array}{l}
\\
\lambda_1
\end{array}
Initial program 99.8%
Taylor expanded in lambda1 around inf
Simplified69.3%
herbie shell --seed 2024198
(FPCore (lambda1 phi1 phi2 delta theta)
:name "Destination given bearing on a great circle"
:precision binary64
(+ lambda1 (atan2 (* (* (sin theta) (sin delta)) (cos phi1)) (- (cos delta) (* (sin phi1) (sin (asin (+ (* (sin phi1) (cos delta)) (* (* (cos phi1) (sin delta)) (cos theta))))))))))