
(FPCore (R lambda1 lambda2 phi1 phi2)
:precision binary64
(*
(acos
(+
(* (sin phi1) (sin phi2))
(* (* (cos phi1) (cos phi2)) (cos (- lambda1 lambda2)))))
R))
double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
return acos(((sin(phi1) * sin(phi2)) + ((cos(phi1) * cos(phi2)) * cos((lambda1 - lambda2))))) * R;
}
real(8) function code(r, lambda1, lambda2, phi1, phi2)
real(8), intent (in) :: r
real(8), intent (in) :: lambda1
real(8), intent (in) :: lambda2
real(8), intent (in) :: phi1
real(8), intent (in) :: phi2
code = acos(((sin(phi1) * sin(phi2)) + ((cos(phi1) * cos(phi2)) * cos((lambda1 - lambda2))))) * r
end function
public static double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
return Math.acos(((Math.sin(phi1) * Math.sin(phi2)) + ((Math.cos(phi1) * Math.cos(phi2)) * Math.cos((lambda1 - lambda2))))) * R;
}
def code(R, lambda1, lambda2, phi1, phi2): return math.acos(((math.sin(phi1) * math.sin(phi2)) + ((math.cos(phi1) * math.cos(phi2)) * math.cos((lambda1 - lambda2))))) * R
function code(R, lambda1, lambda2, phi1, phi2) return Float64(acos(Float64(Float64(sin(phi1) * sin(phi2)) + Float64(Float64(cos(phi1) * cos(phi2)) * cos(Float64(lambda1 - lambda2))))) * R) end
function tmp = code(R, lambda1, lambda2, phi1, phi2) tmp = acos(((sin(phi1) * sin(phi2)) + ((cos(phi1) * cos(phi2)) * cos((lambda1 - lambda2))))) * R; end
code[R_, lambda1_, lambda2_, phi1_, phi2_] := N[(N[ArcCos[N[(N[(N[Sin[phi1], $MachinePrecision] * N[Sin[phi2], $MachinePrecision]), $MachinePrecision] + N[(N[(N[Cos[phi1], $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision] * N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * R), $MachinePrecision]
\begin{array}{l}
\\
\cos^{-1} \left(\sin \phi_1 \cdot \sin \phi_2 + \left(\cos \phi_1 \cdot \cos \phi_2\right) \cdot \cos \left(\lambda_1 - \lambda_2\right)\right) \cdot R
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 30 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (R lambda1 lambda2 phi1 phi2)
:precision binary64
(*
(acos
(+
(* (sin phi1) (sin phi2))
(* (* (cos phi1) (cos phi2)) (cos (- lambda1 lambda2)))))
R))
double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
return acos(((sin(phi1) * sin(phi2)) + ((cos(phi1) * cos(phi2)) * cos((lambda1 - lambda2))))) * R;
}
real(8) function code(r, lambda1, lambda2, phi1, phi2)
real(8), intent (in) :: r
real(8), intent (in) :: lambda1
real(8), intent (in) :: lambda2
real(8), intent (in) :: phi1
real(8), intent (in) :: phi2
code = acos(((sin(phi1) * sin(phi2)) + ((cos(phi1) * cos(phi2)) * cos((lambda1 - lambda2))))) * r
end function
public static double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
return Math.acos(((Math.sin(phi1) * Math.sin(phi2)) + ((Math.cos(phi1) * Math.cos(phi2)) * Math.cos((lambda1 - lambda2))))) * R;
}
def code(R, lambda1, lambda2, phi1, phi2): return math.acos(((math.sin(phi1) * math.sin(phi2)) + ((math.cos(phi1) * math.cos(phi2)) * math.cos((lambda1 - lambda2))))) * R
function code(R, lambda1, lambda2, phi1, phi2) return Float64(acos(Float64(Float64(sin(phi1) * sin(phi2)) + Float64(Float64(cos(phi1) * cos(phi2)) * cos(Float64(lambda1 - lambda2))))) * R) end
function tmp = code(R, lambda1, lambda2, phi1, phi2) tmp = acos(((sin(phi1) * sin(phi2)) + ((cos(phi1) * cos(phi2)) * cos((lambda1 - lambda2))))) * R; end
code[R_, lambda1_, lambda2_, phi1_, phi2_] := N[(N[ArcCos[N[(N[(N[Sin[phi1], $MachinePrecision] * N[Sin[phi2], $MachinePrecision]), $MachinePrecision] + N[(N[(N[Cos[phi1], $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision] * N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * R), $MachinePrecision]
\begin{array}{l}
\\
\cos^{-1} \left(\sin \phi_1 \cdot \sin \phi_2 + \left(\cos \phi_1 \cdot \cos \phi_2\right) \cdot \cos \left(\lambda_1 - \lambda_2\right)\right) \cdot R
\end{array}
(FPCore (R lambda1 lambda2 phi1 phi2)
:precision binary64
(*
(acos
(fma
(* (* (cos phi2) (cos phi1)) (cos lambda2))
(cos lambda1)
(fma
(* (* (sin lambda1) (sin lambda2)) (cos phi2))
(cos phi1)
(* (sin phi2) (sin phi1)))))
R))
double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
return acos(fma(((cos(phi2) * cos(phi1)) * cos(lambda2)), cos(lambda1), fma(((sin(lambda1) * sin(lambda2)) * cos(phi2)), cos(phi1), (sin(phi2) * sin(phi1))))) * R;
}
function code(R, lambda1, lambda2, phi1, phi2) return Float64(acos(fma(Float64(Float64(cos(phi2) * cos(phi1)) * cos(lambda2)), cos(lambda1), fma(Float64(Float64(sin(lambda1) * sin(lambda2)) * cos(phi2)), cos(phi1), Float64(sin(phi2) * sin(phi1))))) * R) end
code[R_, lambda1_, lambda2_, phi1_, phi2_] := N[(N[ArcCos[N[(N[(N[(N[Cos[phi2], $MachinePrecision] * N[Cos[phi1], $MachinePrecision]), $MachinePrecision] * N[Cos[lambda2], $MachinePrecision]), $MachinePrecision] * N[Cos[lambda1], $MachinePrecision] + N[(N[(N[(N[Sin[lambda1], $MachinePrecision] * N[Sin[lambda2], $MachinePrecision]), $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision] * N[Cos[phi1], $MachinePrecision] + N[(N[Sin[phi2], $MachinePrecision] * N[Sin[phi1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * R), $MachinePrecision]
\begin{array}{l}
\\
\cos^{-1} \left(\mathsf{fma}\left(\left(\cos \phi_2 \cdot \cos \phi_1\right) \cdot \cos \lambda_2, \cos \lambda_1, \mathsf{fma}\left(\left(\sin \lambda_1 \cdot \sin \lambda_2\right) \cdot \cos \phi_2, \cos \phi_1, \sin \phi_2 \cdot \sin \phi_1\right)\right)\right) \cdot R
\end{array}
Initial program 73.9%
lift-cos.f64N/A
lift--.f64N/A
cos-diffN/A
*-commutativeN/A
lower-fma.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-sin.f64N/A
lower-sin.f6496.2
Applied rewrites96.2%
lift-+.f64N/A
+-commutativeN/A
lift-*.f64N/A
lift-fma.f64N/A
lift-*.f64N/A
distribute-lft-inN/A
associate-+l+N/A
lift-*.f64N/A
associate-*r*N/A
Applied rewrites96.3%
(FPCore (R lambda1 lambda2 phi1 phi2)
:precision binary64
(*
(acos
(fma
(* (cos lambda1) (* (cos phi2) (cos phi1)))
(cos lambda2)
(fma
(* (* (sin lambda1) (sin lambda2)) (cos phi2))
(cos phi1)
(* (sin phi2) (sin phi1)))))
R))
double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
return acos(fma((cos(lambda1) * (cos(phi2) * cos(phi1))), cos(lambda2), fma(((sin(lambda1) * sin(lambda2)) * cos(phi2)), cos(phi1), (sin(phi2) * sin(phi1))))) * R;
}
function code(R, lambda1, lambda2, phi1, phi2) return Float64(acos(fma(Float64(cos(lambda1) * Float64(cos(phi2) * cos(phi1))), cos(lambda2), fma(Float64(Float64(sin(lambda1) * sin(lambda2)) * cos(phi2)), cos(phi1), Float64(sin(phi2) * sin(phi1))))) * R) end
code[R_, lambda1_, lambda2_, phi1_, phi2_] := N[(N[ArcCos[N[(N[(N[Cos[lambda1], $MachinePrecision] * N[(N[Cos[phi2], $MachinePrecision] * N[Cos[phi1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[Cos[lambda2], $MachinePrecision] + N[(N[(N[(N[Sin[lambda1], $MachinePrecision] * N[Sin[lambda2], $MachinePrecision]), $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision] * N[Cos[phi1], $MachinePrecision] + N[(N[Sin[phi2], $MachinePrecision] * N[Sin[phi1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * R), $MachinePrecision]
\begin{array}{l}
\\
\cos^{-1} \left(\mathsf{fma}\left(\cos \lambda_1 \cdot \left(\cos \phi_2 \cdot \cos \phi_1\right), \cos \lambda_2, \mathsf{fma}\left(\left(\sin \lambda_1 \cdot \sin \lambda_2\right) \cdot \cos \phi_2, \cos \phi_1, \sin \phi_2 \cdot \sin \phi_1\right)\right)\right) \cdot R
\end{array}
Initial program 73.9%
lift-cos.f64N/A
lift--.f64N/A
cos-diffN/A
*-commutativeN/A
lower-fma.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-sin.f64N/A
lower-sin.f6496.2
Applied rewrites96.2%
lift-+.f64N/A
+-commutativeN/A
lift-*.f64N/A
lift-fma.f64N/A
lift-*.f64N/A
distribute-lft-inN/A
associate-+l+N/A
Applied rewrites96.3%
(FPCore (R lambda1 lambda2 phi1 phi2)
:precision binary64
(*
(acos
(+
(* (sin phi1) (sin phi2))
(*
(* (cos phi1) (cos phi2))
(fma (cos lambda2) (cos lambda1) (* (sin lambda2) (sin lambda1))))))
R))
double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
return acos(((sin(phi1) * sin(phi2)) + ((cos(phi1) * cos(phi2)) * fma(cos(lambda2), cos(lambda1), (sin(lambda2) * sin(lambda1)))))) * R;
}
function code(R, lambda1, lambda2, phi1, phi2) return Float64(acos(Float64(Float64(sin(phi1) * sin(phi2)) + Float64(Float64(cos(phi1) * cos(phi2)) * fma(cos(lambda2), cos(lambda1), Float64(sin(lambda2) * sin(lambda1)))))) * R) end
code[R_, lambda1_, lambda2_, phi1_, phi2_] := N[(N[ArcCos[N[(N[(N[Sin[phi1], $MachinePrecision] * N[Sin[phi2], $MachinePrecision]), $MachinePrecision] + N[(N[(N[Cos[phi1], $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision] * N[(N[Cos[lambda2], $MachinePrecision] * N[Cos[lambda1], $MachinePrecision] + N[(N[Sin[lambda2], $MachinePrecision] * N[Sin[lambda1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * R), $MachinePrecision]
\begin{array}{l}
\\
\cos^{-1} \left(\sin \phi_1 \cdot \sin \phi_2 + \left(\cos \phi_1 \cdot \cos \phi_2\right) \cdot \mathsf{fma}\left(\cos \lambda_2, \cos \lambda_1, \sin \lambda_2 \cdot \sin \lambda_1\right)\right) \cdot R
\end{array}
Initial program 73.9%
lift-cos.f64N/A
lift--.f64N/A
cos-diffN/A
*-commutativeN/A
lower-fma.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-sin.f64N/A
lower-sin.f6496.2
Applied rewrites96.2%
(FPCore (R lambda1 lambda2 phi1 phi2)
:precision binary64
(*
(acos
(fma
(*
(fma (sin lambda2) (sin lambda1) (* (cos lambda2) (cos lambda1)))
(cos phi2))
(cos phi1)
(* (sin phi2) (sin phi1))))
R))
double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
return acos(fma((fma(sin(lambda2), sin(lambda1), (cos(lambda2) * cos(lambda1))) * cos(phi2)), cos(phi1), (sin(phi2) * sin(phi1)))) * R;
}
function code(R, lambda1, lambda2, phi1, phi2) return Float64(acos(fma(Float64(fma(sin(lambda2), sin(lambda1), Float64(cos(lambda2) * cos(lambda1))) * cos(phi2)), cos(phi1), Float64(sin(phi2) * sin(phi1)))) * R) end
code[R_, lambda1_, lambda2_, phi1_, phi2_] := N[(N[ArcCos[N[(N[(N[(N[Sin[lambda2], $MachinePrecision] * N[Sin[lambda1], $MachinePrecision] + N[(N[Cos[lambda2], $MachinePrecision] * N[Cos[lambda1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision] * N[Cos[phi1], $MachinePrecision] + N[(N[Sin[phi2], $MachinePrecision] * N[Sin[phi1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * R), $MachinePrecision]
\begin{array}{l}
\\
\cos^{-1} \left(\mathsf{fma}\left(\mathsf{fma}\left(\sin \lambda_2, \sin \lambda_1, \cos \lambda_2 \cdot \cos \lambda_1\right) \cdot \cos \phi_2, \cos \phi_1, \sin \phi_2 \cdot \sin \phi_1\right)\right) \cdot R
\end{array}
Initial program 73.9%
lift-cos.f64N/A
lift--.f64N/A
cos-diffN/A
*-commutativeN/A
lower-fma.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-sin.f64N/A
lower-sin.f6496.2
Applied rewrites96.2%
Taylor expanded in lambda1 around inf
*-commutativeN/A
lower-fma.f64N/A
Applied rewrites96.2%
Final simplification96.2%
(FPCore (R lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (* (cos phi1) (cos phi2))))
(if (<= phi1 -125.0)
(*
(acos (+ (* (sin phi1) (sin phi2)) (* t_0 (cos (- lambda1 lambda2)))))
R)
(if (<= phi1 5.5e-134)
(*
(acos
(+
(* (* (fma -0.16666666666666666 (* phi1 phi1) 1.0) (sin phi2)) phi1)
(*
t_0
(fma (cos lambda2) (cos lambda1) (* (sin lambda2) (sin lambda1))))))
R)
(*
(acos
(fma
(* (cos (- lambda2 lambda1)) (cos phi1))
(cos phi2)
(* (sin phi2) (sin phi1))))
R)))))
double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = cos(phi1) * cos(phi2);
double tmp;
if (phi1 <= -125.0) {
tmp = acos(((sin(phi1) * sin(phi2)) + (t_0 * cos((lambda1 - lambda2))))) * R;
} else if (phi1 <= 5.5e-134) {
tmp = acos((((fma(-0.16666666666666666, (phi1 * phi1), 1.0) * sin(phi2)) * phi1) + (t_0 * fma(cos(lambda2), cos(lambda1), (sin(lambda2) * sin(lambda1)))))) * R;
} else {
tmp = acos(fma((cos((lambda2 - lambda1)) * cos(phi1)), cos(phi2), (sin(phi2) * sin(phi1)))) * R;
}
return tmp;
}
function code(R, lambda1, lambda2, phi1, phi2) t_0 = Float64(cos(phi1) * cos(phi2)) tmp = 0.0 if (phi1 <= -125.0) tmp = Float64(acos(Float64(Float64(sin(phi1) * sin(phi2)) + Float64(t_0 * cos(Float64(lambda1 - lambda2))))) * R); elseif (phi1 <= 5.5e-134) tmp = Float64(acos(Float64(Float64(Float64(fma(-0.16666666666666666, Float64(phi1 * phi1), 1.0) * sin(phi2)) * phi1) + Float64(t_0 * fma(cos(lambda2), cos(lambda1), Float64(sin(lambda2) * sin(lambda1)))))) * R); else tmp = Float64(acos(fma(Float64(cos(Float64(lambda2 - lambda1)) * cos(phi1)), cos(phi2), Float64(sin(phi2) * sin(phi1)))) * R); end return tmp end
code[R_, lambda1_, lambda2_, phi1_, phi2_] := Block[{t$95$0 = N[(N[Cos[phi1], $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[phi1, -125.0], N[(N[ArcCos[N[(N[(N[Sin[phi1], $MachinePrecision] * N[Sin[phi2], $MachinePrecision]), $MachinePrecision] + N[(t$95$0 * N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * R), $MachinePrecision], If[LessEqual[phi1, 5.5e-134], N[(N[ArcCos[N[(N[(N[(N[(-0.16666666666666666 * N[(phi1 * phi1), $MachinePrecision] + 1.0), $MachinePrecision] * N[Sin[phi2], $MachinePrecision]), $MachinePrecision] * phi1), $MachinePrecision] + N[(t$95$0 * N[(N[Cos[lambda2], $MachinePrecision] * N[Cos[lambda1], $MachinePrecision] + N[(N[Sin[lambda2], $MachinePrecision] * N[Sin[lambda1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * R), $MachinePrecision], N[(N[ArcCos[N[(N[(N[Cos[N[(lambda2 - lambda1), $MachinePrecision]], $MachinePrecision] * N[Cos[phi1], $MachinePrecision]), $MachinePrecision] * N[Cos[phi2], $MachinePrecision] + N[(N[Sin[phi2], $MachinePrecision] * N[Sin[phi1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * R), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \cos \phi_1 \cdot \cos \phi_2\\
\mathbf{if}\;\phi_1 \leq -125:\\
\;\;\;\;\cos^{-1} \left(\sin \phi_1 \cdot \sin \phi_2 + t\_0 \cdot \cos \left(\lambda_1 - \lambda_2\right)\right) \cdot R\\
\mathbf{elif}\;\phi_1 \leq 5.5 \cdot 10^{-134}:\\
\;\;\;\;\cos^{-1} \left(\left(\mathsf{fma}\left(-0.16666666666666666, \phi_1 \cdot \phi_1, 1\right) \cdot \sin \phi_2\right) \cdot \phi_1 + t\_0 \cdot \mathsf{fma}\left(\cos \lambda_2, \cos \lambda_1, \sin \lambda_2 \cdot \sin \lambda_1\right)\right) \cdot R\\
\mathbf{else}:\\
\;\;\;\;\cos^{-1} \left(\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right) \cdot \cos \phi_1, \cos \phi_2, \sin \phi_2 \cdot \sin \phi_1\right)\right) \cdot R\\
\end{array}
\end{array}
if phi1 < -125Initial program 84.5%
if -125 < phi1 < 5.5000000000000002e-134Initial program 65.3%
lift-cos.f64N/A
lift--.f64N/A
cos-diffN/A
*-commutativeN/A
lower-fma.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-sin.f64N/A
lower-sin.f6492.9
Applied rewrites92.9%
Taylor expanded in phi1 around 0
*-commutativeN/A
lower-*.f64N/A
associate-*r*N/A
distribute-rgt1-inN/A
lower-*.f64N/A
metadata-evalN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
metadata-evalN/A
lower-sin.f6492.5
Applied rewrites92.5%
if 5.5000000000000002e-134 < phi1 Initial program 75.9%
lift-+.f64N/A
+-commutativeN/A
lift-*.f64N/A
*-commutativeN/A
lift-*.f64N/A
associate-*r*N/A
lower-fma.f64N/A
lower-*.f6475.9
lift-cos.f64N/A
lift--.f64N/A
cos-diffN/A
*-commutativeN/A
*-commutativeN/A
cos-diff-revN/A
lower-cos.f64N/A
lower--.f6475.9
lift-*.f64N/A
*-commutativeN/A
lower-*.f6475.9
Applied rewrites75.9%
(FPCore (R lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (* (cos phi1) (cos phi2))))
(if (<= phi1 -125.0)
(*
(acos (+ (* (sin phi1) (sin phi2)) (* t_0 (cos (- lambda1 lambda2)))))
R)
(if (<= phi1 1.8e-56)
(*
(acos
(+
(* (sin phi2) phi1)
(*
t_0
(fma (cos lambda2) (cos lambda1) (* (sin lambda2) (sin lambda1))))))
R)
(*
(acos
(fma
(* (cos (- lambda2 lambda1)) (cos phi1))
(cos phi2)
(* (sin phi2) (sin phi1))))
R)))))
double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = cos(phi1) * cos(phi2);
double tmp;
if (phi1 <= -125.0) {
tmp = acos(((sin(phi1) * sin(phi2)) + (t_0 * cos((lambda1 - lambda2))))) * R;
} else if (phi1 <= 1.8e-56) {
tmp = acos(((sin(phi2) * phi1) + (t_0 * fma(cos(lambda2), cos(lambda1), (sin(lambda2) * sin(lambda1)))))) * R;
} else {
tmp = acos(fma((cos((lambda2 - lambda1)) * cos(phi1)), cos(phi2), (sin(phi2) * sin(phi1)))) * R;
}
return tmp;
}
function code(R, lambda1, lambda2, phi1, phi2) t_0 = Float64(cos(phi1) * cos(phi2)) tmp = 0.0 if (phi1 <= -125.0) tmp = Float64(acos(Float64(Float64(sin(phi1) * sin(phi2)) + Float64(t_0 * cos(Float64(lambda1 - lambda2))))) * R); elseif (phi1 <= 1.8e-56) tmp = Float64(acos(Float64(Float64(sin(phi2) * phi1) + Float64(t_0 * fma(cos(lambda2), cos(lambda1), Float64(sin(lambda2) * sin(lambda1)))))) * R); else tmp = Float64(acos(fma(Float64(cos(Float64(lambda2 - lambda1)) * cos(phi1)), cos(phi2), Float64(sin(phi2) * sin(phi1)))) * R); end return tmp end
code[R_, lambda1_, lambda2_, phi1_, phi2_] := Block[{t$95$0 = N[(N[Cos[phi1], $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[phi1, -125.0], N[(N[ArcCos[N[(N[(N[Sin[phi1], $MachinePrecision] * N[Sin[phi2], $MachinePrecision]), $MachinePrecision] + N[(t$95$0 * N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * R), $MachinePrecision], If[LessEqual[phi1, 1.8e-56], N[(N[ArcCos[N[(N[(N[Sin[phi2], $MachinePrecision] * phi1), $MachinePrecision] + N[(t$95$0 * N[(N[Cos[lambda2], $MachinePrecision] * N[Cos[lambda1], $MachinePrecision] + N[(N[Sin[lambda2], $MachinePrecision] * N[Sin[lambda1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * R), $MachinePrecision], N[(N[ArcCos[N[(N[(N[Cos[N[(lambda2 - lambda1), $MachinePrecision]], $MachinePrecision] * N[Cos[phi1], $MachinePrecision]), $MachinePrecision] * N[Cos[phi2], $MachinePrecision] + N[(N[Sin[phi2], $MachinePrecision] * N[Sin[phi1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * R), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \cos \phi_1 \cdot \cos \phi_2\\
\mathbf{if}\;\phi_1 \leq -125:\\
\;\;\;\;\cos^{-1} \left(\sin \phi_1 \cdot \sin \phi_2 + t\_0 \cdot \cos \left(\lambda_1 - \lambda_2\right)\right) \cdot R\\
\mathbf{elif}\;\phi_1 \leq 1.8 \cdot 10^{-56}:\\
\;\;\;\;\cos^{-1} \left(\sin \phi_2 \cdot \phi_1 + t\_0 \cdot \mathsf{fma}\left(\cos \lambda_2, \cos \lambda_1, \sin \lambda_2 \cdot \sin \lambda_1\right)\right) \cdot R\\
\mathbf{else}:\\
\;\;\;\;\cos^{-1} \left(\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right) \cdot \cos \phi_1, \cos \phi_2, \sin \phi_2 \cdot \sin \phi_1\right)\right) \cdot R\\
\end{array}
\end{array}
if phi1 < -125Initial program 84.5%
if -125 < phi1 < 1.79999999999999989e-56Initial program 65.3%
lift-cos.f64N/A
lift--.f64N/A
cos-diffN/A
*-commutativeN/A
lower-fma.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-sin.f64N/A
lower-sin.f6493.1
Applied rewrites93.1%
Taylor expanded in phi1 around 0
*-commutativeN/A
lower-*.f64N/A
lower-sin.f6492.6
Applied rewrites92.6%
if 1.79999999999999989e-56 < phi1 Initial program 78.1%
lift-+.f64N/A
+-commutativeN/A
lift-*.f64N/A
*-commutativeN/A
lift-*.f64N/A
associate-*r*N/A
lower-fma.f64N/A
lower-*.f6478.1
lift-cos.f64N/A
lift--.f64N/A
cos-diffN/A
*-commutativeN/A
*-commutativeN/A
cos-diff-revN/A
lower-cos.f64N/A
lower--.f6478.1
lift-*.f64N/A
*-commutativeN/A
lower-*.f6478.1
Applied rewrites78.1%
(FPCore (R lambda1 lambda2 phi1 phi2)
:precision binary64
(if (<= phi1 -0.095)
(*
(acos
(+
(* (sin phi1) (sin phi2))
(* (* (cos phi1) (cos phi2)) (cos (- lambda1 lambda2)))))
R)
(if (<= phi1 1.8e-56)
(*
(acos
(fma
(fma phi1 (* -0.5 phi1) 1.0)
(*
(fma (sin lambda2) (sin lambda1) (* (cos lambda2) (cos lambda1)))
(cos phi2))
(* (sin phi2) phi1)))
R)
(*
(acos
(fma
(* (cos (- lambda2 lambda1)) (cos phi1))
(cos phi2)
(* (sin phi2) (sin phi1))))
R))))
double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
double tmp;
if (phi1 <= -0.095) {
tmp = acos(((sin(phi1) * sin(phi2)) + ((cos(phi1) * cos(phi2)) * cos((lambda1 - lambda2))))) * R;
} else if (phi1 <= 1.8e-56) {
tmp = acos(fma(fma(phi1, (-0.5 * phi1), 1.0), (fma(sin(lambda2), sin(lambda1), (cos(lambda2) * cos(lambda1))) * cos(phi2)), (sin(phi2) * phi1))) * R;
} else {
tmp = acos(fma((cos((lambda2 - lambda1)) * cos(phi1)), cos(phi2), (sin(phi2) * sin(phi1)))) * R;
}
return tmp;
}
function code(R, lambda1, lambda2, phi1, phi2) tmp = 0.0 if (phi1 <= -0.095) tmp = Float64(acos(Float64(Float64(sin(phi1) * sin(phi2)) + Float64(Float64(cos(phi1) * cos(phi2)) * cos(Float64(lambda1 - lambda2))))) * R); elseif (phi1 <= 1.8e-56) tmp = Float64(acos(fma(fma(phi1, Float64(-0.5 * phi1), 1.0), Float64(fma(sin(lambda2), sin(lambda1), Float64(cos(lambda2) * cos(lambda1))) * cos(phi2)), Float64(sin(phi2) * phi1))) * R); else tmp = Float64(acos(fma(Float64(cos(Float64(lambda2 - lambda1)) * cos(phi1)), cos(phi2), Float64(sin(phi2) * sin(phi1)))) * R); end return tmp end
code[R_, lambda1_, lambda2_, phi1_, phi2_] := If[LessEqual[phi1, -0.095], N[(N[ArcCos[N[(N[(N[Sin[phi1], $MachinePrecision] * N[Sin[phi2], $MachinePrecision]), $MachinePrecision] + N[(N[(N[Cos[phi1], $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision] * N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * R), $MachinePrecision], If[LessEqual[phi1, 1.8e-56], N[(N[ArcCos[N[(N[(phi1 * N[(-0.5 * phi1), $MachinePrecision] + 1.0), $MachinePrecision] * N[(N[(N[Sin[lambda2], $MachinePrecision] * N[Sin[lambda1], $MachinePrecision] + N[(N[Cos[lambda2], $MachinePrecision] * N[Cos[lambda1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision] + N[(N[Sin[phi2], $MachinePrecision] * phi1), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * R), $MachinePrecision], N[(N[ArcCos[N[(N[(N[Cos[N[(lambda2 - lambda1), $MachinePrecision]], $MachinePrecision] * N[Cos[phi1], $MachinePrecision]), $MachinePrecision] * N[Cos[phi2], $MachinePrecision] + N[(N[Sin[phi2], $MachinePrecision] * N[Sin[phi1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * R), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\phi_1 \leq -0.095:\\
\;\;\;\;\cos^{-1} \left(\sin \phi_1 \cdot \sin \phi_2 + \left(\cos \phi_1 \cdot \cos \phi_2\right) \cdot \cos \left(\lambda_1 - \lambda_2\right)\right) \cdot R\\
\mathbf{elif}\;\phi_1 \leq 1.8 \cdot 10^{-56}:\\
\;\;\;\;\cos^{-1} \left(\mathsf{fma}\left(\mathsf{fma}\left(\phi_1, -0.5 \cdot \phi_1, 1\right), \mathsf{fma}\left(\sin \lambda_2, \sin \lambda_1, \cos \lambda_2 \cdot \cos \lambda_1\right) \cdot \cos \phi_2, \sin \phi_2 \cdot \phi_1\right)\right) \cdot R\\
\mathbf{else}:\\
\;\;\;\;\cos^{-1} \left(\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right) \cdot \cos \phi_1, \cos \phi_2, \sin \phi_2 \cdot \sin \phi_1\right)\right) \cdot R\\
\end{array}
\end{array}
if phi1 < -0.095000000000000001Initial program 84.5%
if -0.095000000000000001 < phi1 < 1.79999999999999989e-56Initial program 65.3%
Taylor expanded in phi1 around 0
Applied rewrites65.3%
Applied rewrites92.5%
if 1.79999999999999989e-56 < phi1 Initial program 78.1%
lift-+.f64N/A
+-commutativeN/A
lift-*.f64N/A
*-commutativeN/A
lift-*.f64N/A
associate-*r*N/A
lower-fma.f64N/A
lower-*.f6478.1
lift-cos.f64N/A
lift--.f64N/A
cos-diffN/A
*-commutativeN/A
*-commutativeN/A
cos-diff-revN/A
lower-cos.f64N/A
lower--.f6478.1
lift-*.f64N/A
*-commutativeN/A
lower-*.f6478.1
Applied rewrites78.1%
(FPCore (R lambda1 lambda2 phi1 phi2)
:precision binary64
(if (<= phi1 -0.095)
(*
(acos
(+
(* (sin phi1) (sin phi2))
(* (* (cos phi1) (cos phi2)) (cos (- lambda1 lambda2)))))
R)
(if (<= phi1 5.5e-134)
(*
(acos
(fma
(cos phi2)
(fma (cos lambda1) (cos lambda2) (* (sin lambda2) (sin lambda1)))
(* (sin phi2) phi1)))
R)
(*
(acos
(fma
(* (cos (- lambda2 lambda1)) (cos phi1))
(cos phi2)
(* (sin phi2) (sin phi1))))
R))))
double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
double tmp;
if (phi1 <= -0.095) {
tmp = acos(((sin(phi1) * sin(phi2)) + ((cos(phi1) * cos(phi2)) * cos((lambda1 - lambda2))))) * R;
} else if (phi1 <= 5.5e-134) {
tmp = acos(fma(cos(phi2), fma(cos(lambda1), cos(lambda2), (sin(lambda2) * sin(lambda1))), (sin(phi2) * phi1))) * R;
} else {
tmp = acos(fma((cos((lambda2 - lambda1)) * cos(phi1)), cos(phi2), (sin(phi2) * sin(phi1)))) * R;
}
return tmp;
}
function code(R, lambda1, lambda2, phi1, phi2) tmp = 0.0 if (phi1 <= -0.095) tmp = Float64(acos(Float64(Float64(sin(phi1) * sin(phi2)) + Float64(Float64(cos(phi1) * cos(phi2)) * cos(Float64(lambda1 - lambda2))))) * R); elseif (phi1 <= 5.5e-134) tmp = Float64(acos(fma(cos(phi2), fma(cos(lambda1), cos(lambda2), Float64(sin(lambda2) * sin(lambda1))), Float64(sin(phi2) * phi1))) * R); else tmp = Float64(acos(fma(Float64(cos(Float64(lambda2 - lambda1)) * cos(phi1)), cos(phi2), Float64(sin(phi2) * sin(phi1)))) * R); end return tmp end
code[R_, lambda1_, lambda2_, phi1_, phi2_] := If[LessEqual[phi1, -0.095], N[(N[ArcCos[N[(N[(N[Sin[phi1], $MachinePrecision] * N[Sin[phi2], $MachinePrecision]), $MachinePrecision] + N[(N[(N[Cos[phi1], $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision] * N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * R), $MachinePrecision], If[LessEqual[phi1, 5.5e-134], N[(N[ArcCos[N[(N[Cos[phi2], $MachinePrecision] * N[(N[Cos[lambda1], $MachinePrecision] * N[Cos[lambda2], $MachinePrecision] + N[(N[Sin[lambda2], $MachinePrecision] * N[Sin[lambda1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[Sin[phi2], $MachinePrecision] * phi1), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * R), $MachinePrecision], N[(N[ArcCos[N[(N[(N[Cos[N[(lambda2 - lambda1), $MachinePrecision]], $MachinePrecision] * N[Cos[phi1], $MachinePrecision]), $MachinePrecision] * N[Cos[phi2], $MachinePrecision] + N[(N[Sin[phi2], $MachinePrecision] * N[Sin[phi1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * R), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\phi_1 \leq -0.095:\\
\;\;\;\;\cos^{-1} \left(\sin \phi_1 \cdot \sin \phi_2 + \left(\cos \phi_1 \cdot \cos \phi_2\right) \cdot \cos \left(\lambda_1 - \lambda_2\right)\right) \cdot R\\
\mathbf{elif}\;\phi_1 \leq 5.5 \cdot 10^{-134}:\\
\;\;\;\;\cos^{-1} \left(\mathsf{fma}\left(\cos \phi_2, \mathsf{fma}\left(\cos \lambda_1, \cos \lambda_2, \sin \lambda_2 \cdot \sin \lambda_1\right), \sin \phi_2 \cdot \phi_1\right)\right) \cdot R\\
\mathbf{else}:\\
\;\;\;\;\cos^{-1} \left(\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right) \cdot \cos \phi_1, \cos \phi_2, \sin \phi_2 \cdot \sin \phi_1\right)\right) \cdot R\\
\end{array}
\end{array}
if phi1 < -0.095000000000000001Initial program 84.5%
if -0.095000000000000001 < phi1 < 5.5000000000000002e-134Initial program 65.3%
lift-cos.f64N/A
lift--.f64N/A
cos-diffN/A
*-commutativeN/A
lower-fma.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-sin.f64N/A
lower-sin.f6492.9
Applied rewrites92.9%
lift-+.f64N/A
+-commutativeN/A
lift-*.f64N/A
lift-fma.f64N/A
lift-*.f64N/A
distribute-lft-inN/A
associate-+l+N/A
lift-*.f64N/A
associate-*r*N/A
Applied rewrites92.9%
Taylor expanded in phi1 around 0
+-commutativeN/A
associate-*r*N/A
*-commutativeN/A
distribute-rgt-outN/A
lower-fma.f64N/A
lower-cos.f64N/A
lower-fma.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-sin.f64N/A
lower-sin.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-sin.f6492.0
Applied rewrites92.0%
if 5.5000000000000002e-134 < phi1 Initial program 75.9%
lift-+.f64N/A
+-commutativeN/A
lift-*.f64N/A
*-commutativeN/A
lift-*.f64N/A
associate-*r*N/A
lower-fma.f64N/A
lower-*.f6475.9
lift-cos.f64N/A
lift--.f64N/A
cos-diffN/A
*-commutativeN/A
*-commutativeN/A
cos-diff-revN/A
lower-cos.f64N/A
lower--.f6475.9
lift-*.f64N/A
*-commutativeN/A
lower-*.f6475.9
Applied rewrites75.9%
(FPCore (R lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (* (cos phi2) (cos phi1))))
(if (<= phi2 -0.00013)
(* (acos (fma t_0 (cos lambda1) (* (sin phi2) (sin phi1)))) R)
(if (<= phi2 4.2e-16)
(*
(acos
(*
(cos phi1)
(fma (cos lambda1) (cos lambda2) (* (sin lambda2) (sin lambda1)))))
R)
(if (<= phi2 6.6e+82)
(*
(acos
(*
(fma (sin lambda2) (sin lambda1) (* (cos lambda2) (cos lambda1)))
(cos phi2)))
R)
(* (acos (fma (sin phi2) (sin phi1) (* (cos lambda1) t_0))) R))))))
double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = cos(phi2) * cos(phi1);
double tmp;
if (phi2 <= -0.00013) {
tmp = acos(fma(t_0, cos(lambda1), (sin(phi2) * sin(phi1)))) * R;
} else if (phi2 <= 4.2e-16) {
tmp = acos((cos(phi1) * fma(cos(lambda1), cos(lambda2), (sin(lambda2) * sin(lambda1))))) * R;
} else if (phi2 <= 6.6e+82) {
tmp = acos((fma(sin(lambda2), sin(lambda1), (cos(lambda2) * cos(lambda1))) * cos(phi2))) * R;
} else {
tmp = acos(fma(sin(phi2), sin(phi1), (cos(lambda1) * t_0))) * R;
}
return tmp;
}
function code(R, lambda1, lambda2, phi1, phi2) t_0 = Float64(cos(phi2) * cos(phi1)) tmp = 0.0 if (phi2 <= -0.00013) tmp = Float64(acos(fma(t_0, cos(lambda1), Float64(sin(phi2) * sin(phi1)))) * R); elseif (phi2 <= 4.2e-16) tmp = Float64(acos(Float64(cos(phi1) * fma(cos(lambda1), cos(lambda2), Float64(sin(lambda2) * sin(lambda1))))) * R); elseif (phi2 <= 6.6e+82) tmp = Float64(acos(Float64(fma(sin(lambda2), sin(lambda1), Float64(cos(lambda2) * cos(lambda1))) * cos(phi2))) * R); else tmp = Float64(acos(fma(sin(phi2), sin(phi1), Float64(cos(lambda1) * t_0))) * R); end return tmp end
code[R_, lambda1_, lambda2_, phi1_, phi2_] := Block[{t$95$0 = N[(N[Cos[phi2], $MachinePrecision] * N[Cos[phi1], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[phi2, -0.00013], N[(N[ArcCos[N[(t$95$0 * N[Cos[lambda1], $MachinePrecision] + N[(N[Sin[phi2], $MachinePrecision] * N[Sin[phi1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * R), $MachinePrecision], If[LessEqual[phi2, 4.2e-16], N[(N[ArcCos[N[(N[Cos[phi1], $MachinePrecision] * N[(N[Cos[lambda1], $MachinePrecision] * N[Cos[lambda2], $MachinePrecision] + N[(N[Sin[lambda2], $MachinePrecision] * N[Sin[lambda1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * R), $MachinePrecision], If[LessEqual[phi2, 6.6e+82], N[(N[ArcCos[N[(N[(N[Sin[lambda2], $MachinePrecision] * N[Sin[lambda1], $MachinePrecision] + N[(N[Cos[lambda2], $MachinePrecision] * N[Cos[lambda1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * R), $MachinePrecision], N[(N[ArcCos[N[(N[Sin[phi2], $MachinePrecision] * N[Sin[phi1], $MachinePrecision] + N[(N[Cos[lambda1], $MachinePrecision] * t$95$0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * R), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \cos \phi_2 \cdot \cos \phi_1\\
\mathbf{if}\;\phi_2 \leq -0.00013:\\
\;\;\;\;\cos^{-1} \left(\mathsf{fma}\left(t\_0, \cos \lambda_1, \sin \phi_2 \cdot \sin \phi_1\right)\right) \cdot R\\
\mathbf{elif}\;\phi_2 \leq 4.2 \cdot 10^{-16}:\\
\;\;\;\;\cos^{-1} \left(\cos \phi_1 \cdot \mathsf{fma}\left(\cos \lambda_1, \cos \lambda_2, \sin \lambda_2 \cdot \sin \lambda_1\right)\right) \cdot R\\
\mathbf{elif}\;\phi_2 \leq 6.6 \cdot 10^{+82}:\\
\;\;\;\;\cos^{-1} \left(\mathsf{fma}\left(\sin \lambda_2, \sin \lambda_1, \cos \lambda_2 \cdot \cos \lambda_1\right) \cdot \cos \phi_2\right) \cdot R\\
\mathbf{else}:\\
\;\;\;\;\cos^{-1} \left(\mathsf{fma}\left(\sin \phi_2, \sin \phi_1, \cos \lambda_1 \cdot t\_0\right)\right) \cdot R\\
\end{array}
\end{array}
if phi2 < -1.29999999999999989e-4Initial program 81.2%
Taylor expanded in lambda2 around 0
*-commutativeN/A
fp-cancel-sign-sub-invN/A
fp-cancel-sub-sign-invN/A
*-commutativeN/A
sin-neg-revN/A
sin-neg-revN/A
remove-double-negN/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-sin.f64N/A
lower-sin.f6460.5
Applied rewrites60.5%
if -1.29999999999999989e-4 < phi2 < 4.2000000000000002e-16Initial program 67.3%
lift-cos.f64N/A
lift--.f64N/A
cos-diffN/A
*-commutativeN/A
lower-fma.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-sin.f64N/A
lower-sin.f6493.9
Applied rewrites93.9%
lift-+.f64N/A
+-commutativeN/A
lift-*.f64N/A
lift-fma.f64N/A
lift-*.f64N/A
distribute-lft-inN/A
associate-+l+N/A
lift-*.f64N/A
associate-*r*N/A
Applied rewrites93.9%
Taylor expanded in phi2 around 0
associate-*r*N/A
*-commutativeN/A
distribute-rgt-outN/A
lower-*.f64N/A
lower-cos.f64N/A
lower-fma.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-sin.f64N/A
lower-sin.f6493.9
Applied rewrites93.9%
if 4.2000000000000002e-16 < phi2 < 6.5999999999999997e82Initial program 80.6%
lift-cos.f64N/A
lift--.f64N/A
cos-diffN/A
*-commutativeN/A
lower-fma.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-sin.f64N/A
lower-sin.f6498.8
Applied rewrites98.8%
Taylor expanded in phi1 around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-sin.f64N/A
lower-sin.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
lower-cos.f6457.5
Applied rewrites57.5%
if 6.5999999999999997e82 < phi2 Initial program 82.8%
Taylor expanded in lambda2 around 0
*-commutativeN/A
fp-cancel-sign-sub-invN/A
fp-cancel-sub-sign-invN/A
*-commutativeN/A
sin-neg-revN/A
sin-neg-revN/A
remove-double-negN/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-sin.f64N/A
lower-sin.f6471.9
Applied rewrites71.9%
Applied rewrites71.9%
Final simplification80.3%
(FPCore (R lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (* (cos phi2) (cos phi1)))
(t_1
(fma (sin lambda2) (sin lambda1) (* (cos lambda2) (cos lambda1)))))
(if (<= phi2 -0.00013)
(* (acos (fma t_0 (cos lambda1) (* (sin phi2) (sin phi1)))) R)
(if (<= phi2 1.2e-33)
(* (acos (* t_1 (cos phi1))) R)
(if (<= phi2 6.6e+82)
(* (acos (* t_1 (cos phi2))) R)
(* (acos (fma (sin phi2) (sin phi1) (* (cos lambda1) t_0))) R))))))
double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = cos(phi2) * cos(phi1);
double t_1 = fma(sin(lambda2), sin(lambda1), (cos(lambda2) * cos(lambda1)));
double tmp;
if (phi2 <= -0.00013) {
tmp = acos(fma(t_0, cos(lambda1), (sin(phi2) * sin(phi1)))) * R;
} else if (phi2 <= 1.2e-33) {
tmp = acos((t_1 * cos(phi1))) * R;
} else if (phi2 <= 6.6e+82) {
tmp = acos((t_1 * cos(phi2))) * R;
} else {
tmp = acos(fma(sin(phi2), sin(phi1), (cos(lambda1) * t_0))) * R;
}
return tmp;
}
function code(R, lambda1, lambda2, phi1, phi2) t_0 = Float64(cos(phi2) * cos(phi1)) t_1 = fma(sin(lambda2), sin(lambda1), Float64(cos(lambda2) * cos(lambda1))) tmp = 0.0 if (phi2 <= -0.00013) tmp = Float64(acos(fma(t_0, cos(lambda1), Float64(sin(phi2) * sin(phi1)))) * R); elseif (phi2 <= 1.2e-33) tmp = Float64(acos(Float64(t_1 * cos(phi1))) * R); elseif (phi2 <= 6.6e+82) tmp = Float64(acos(Float64(t_1 * cos(phi2))) * R); else tmp = Float64(acos(fma(sin(phi2), sin(phi1), Float64(cos(lambda1) * t_0))) * R); end return tmp end
code[R_, lambda1_, lambda2_, phi1_, phi2_] := Block[{t$95$0 = N[(N[Cos[phi2], $MachinePrecision] * N[Cos[phi1], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[Sin[lambda2], $MachinePrecision] * N[Sin[lambda1], $MachinePrecision] + N[(N[Cos[lambda2], $MachinePrecision] * N[Cos[lambda1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[phi2, -0.00013], N[(N[ArcCos[N[(t$95$0 * N[Cos[lambda1], $MachinePrecision] + N[(N[Sin[phi2], $MachinePrecision] * N[Sin[phi1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * R), $MachinePrecision], If[LessEqual[phi2, 1.2e-33], N[(N[ArcCos[N[(t$95$1 * N[Cos[phi1], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * R), $MachinePrecision], If[LessEqual[phi2, 6.6e+82], N[(N[ArcCos[N[(t$95$1 * N[Cos[phi2], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * R), $MachinePrecision], N[(N[ArcCos[N[(N[Sin[phi2], $MachinePrecision] * N[Sin[phi1], $MachinePrecision] + N[(N[Cos[lambda1], $MachinePrecision] * t$95$0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * R), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \cos \phi_2 \cdot \cos \phi_1\\
t_1 := \mathsf{fma}\left(\sin \lambda_2, \sin \lambda_1, \cos \lambda_2 \cdot \cos \lambda_1\right)\\
\mathbf{if}\;\phi_2 \leq -0.00013:\\
\;\;\;\;\cos^{-1} \left(\mathsf{fma}\left(t\_0, \cos \lambda_1, \sin \phi_2 \cdot \sin \phi_1\right)\right) \cdot R\\
\mathbf{elif}\;\phi_2 \leq 1.2 \cdot 10^{-33}:\\
\;\;\;\;\cos^{-1} \left(t\_1 \cdot \cos \phi_1\right) \cdot R\\
\mathbf{elif}\;\phi_2 \leq 6.6 \cdot 10^{+82}:\\
\;\;\;\;\cos^{-1} \left(t\_1 \cdot \cos \phi_2\right) \cdot R\\
\mathbf{else}:\\
\;\;\;\;\cos^{-1} \left(\mathsf{fma}\left(\sin \phi_2, \sin \phi_1, \cos \lambda_1 \cdot t\_0\right)\right) \cdot R\\
\end{array}
\end{array}
if phi2 < -1.29999999999999989e-4Initial program 81.2%
Taylor expanded in lambda2 around 0
*-commutativeN/A
fp-cancel-sign-sub-invN/A
fp-cancel-sub-sign-invN/A
*-commutativeN/A
sin-neg-revN/A
sin-neg-revN/A
remove-double-negN/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-sin.f64N/A
lower-sin.f6460.5
Applied rewrites60.5%
if -1.29999999999999989e-4 < phi2 < 1.2e-33Initial program 66.9%
lift-cos.f64N/A
lift--.f64N/A
cos-diffN/A
*-commutativeN/A
lower-fma.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-sin.f64N/A
lower-sin.f6493.7
Applied rewrites93.7%
Taylor expanded in phi2 around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-sin.f64N/A
lower-sin.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
lower-cos.f6493.6
Applied rewrites93.6%
if 1.2e-33 < phi2 < 6.5999999999999997e82Initial program 80.4%
lift-cos.f64N/A
lift--.f64N/A
cos-diffN/A
*-commutativeN/A
lower-fma.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-sin.f64N/A
lower-sin.f6499.0
Applied rewrites99.0%
Taylor expanded in phi1 around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-sin.f64N/A
lower-sin.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
lower-cos.f6461.7
Applied rewrites61.7%
if 6.5999999999999997e82 < phi2 Initial program 82.8%
Taylor expanded in lambda2 around 0
*-commutativeN/A
fp-cancel-sign-sub-invN/A
fp-cancel-sub-sign-invN/A
*-commutativeN/A
sin-neg-revN/A
sin-neg-revN/A
remove-double-negN/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-sin.f64N/A
lower-sin.f6471.9
Applied rewrites71.9%
Applied rewrites71.9%
Final simplification79.9%
(FPCore (R lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (* (sin phi2) (sin phi1))) (t_1 (* (cos phi2) (cos phi1))))
(if (<= phi2 -0.00013)
(* (acos (fma t_1 (cos lambda1) t_0)) R)
(if (<= phi2 3.2e-14)
(*
(acos
(*
(fma (sin lambda2) (sin lambda1) (* (cos lambda2) (cos lambda1)))
(cos phi1)))
R)
(if (<= phi2 2.7e+81)
(* (acos (fma (* (cos lambda2) (cos phi2)) (cos phi1) t_0)) R)
(* (acos (fma (sin phi2) (sin phi1) (* (cos lambda1) t_1))) R))))))
double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = sin(phi2) * sin(phi1);
double t_1 = cos(phi2) * cos(phi1);
double tmp;
if (phi2 <= -0.00013) {
tmp = acos(fma(t_1, cos(lambda1), t_0)) * R;
} else if (phi2 <= 3.2e-14) {
tmp = acos((fma(sin(lambda2), sin(lambda1), (cos(lambda2) * cos(lambda1))) * cos(phi1))) * R;
} else if (phi2 <= 2.7e+81) {
tmp = acos(fma((cos(lambda2) * cos(phi2)), cos(phi1), t_0)) * R;
} else {
tmp = acos(fma(sin(phi2), sin(phi1), (cos(lambda1) * t_1))) * R;
}
return tmp;
}
function code(R, lambda1, lambda2, phi1, phi2) t_0 = Float64(sin(phi2) * sin(phi1)) t_1 = Float64(cos(phi2) * cos(phi1)) tmp = 0.0 if (phi2 <= -0.00013) tmp = Float64(acos(fma(t_1, cos(lambda1), t_0)) * R); elseif (phi2 <= 3.2e-14) tmp = Float64(acos(Float64(fma(sin(lambda2), sin(lambda1), Float64(cos(lambda2) * cos(lambda1))) * cos(phi1))) * R); elseif (phi2 <= 2.7e+81) tmp = Float64(acos(fma(Float64(cos(lambda2) * cos(phi2)), cos(phi1), t_0)) * R); else tmp = Float64(acos(fma(sin(phi2), sin(phi1), Float64(cos(lambda1) * t_1))) * R); end return tmp end
code[R_, lambda1_, lambda2_, phi1_, phi2_] := Block[{t$95$0 = N[(N[Sin[phi2], $MachinePrecision] * N[Sin[phi1], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[Cos[phi2], $MachinePrecision] * N[Cos[phi1], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[phi2, -0.00013], N[(N[ArcCos[N[(t$95$1 * N[Cos[lambda1], $MachinePrecision] + t$95$0), $MachinePrecision]], $MachinePrecision] * R), $MachinePrecision], If[LessEqual[phi2, 3.2e-14], N[(N[ArcCos[N[(N[(N[Sin[lambda2], $MachinePrecision] * N[Sin[lambda1], $MachinePrecision] + N[(N[Cos[lambda2], $MachinePrecision] * N[Cos[lambda1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[Cos[phi1], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * R), $MachinePrecision], If[LessEqual[phi2, 2.7e+81], N[(N[ArcCos[N[(N[(N[Cos[lambda2], $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision] * N[Cos[phi1], $MachinePrecision] + t$95$0), $MachinePrecision]], $MachinePrecision] * R), $MachinePrecision], N[(N[ArcCos[N[(N[Sin[phi2], $MachinePrecision] * N[Sin[phi1], $MachinePrecision] + N[(N[Cos[lambda1], $MachinePrecision] * t$95$1), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * R), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sin \phi_2 \cdot \sin \phi_1\\
t_1 := \cos \phi_2 \cdot \cos \phi_1\\
\mathbf{if}\;\phi_2 \leq -0.00013:\\
\;\;\;\;\cos^{-1} \left(\mathsf{fma}\left(t\_1, \cos \lambda_1, t\_0\right)\right) \cdot R\\
\mathbf{elif}\;\phi_2 \leq 3.2 \cdot 10^{-14}:\\
\;\;\;\;\cos^{-1} \left(\mathsf{fma}\left(\sin \lambda_2, \sin \lambda_1, \cos \lambda_2 \cdot \cos \lambda_1\right) \cdot \cos \phi_1\right) \cdot R\\
\mathbf{elif}\;\phi_2 \leq 2.7 \cdot 10^{+81}:\\
\;\;\;\;\cos^{-1} \left(\mathsf{fma}\left(\cos \lambda_2 \cdot \cos \phi_2, \cos \phi_1, t\_0\right)\right) \cdot R\\
\mathbf{else}:\\
\;\;\;\;\cos^{-1} \left(\mathsf{fma}\left(\sin \phi_2, \sin \phi_1, \cos \lambda_1 \cdot t\_1\right)\right) \cdot R\\
\end{array}
\end{array}
if phi2 < -1.29999999999999989e-4Initial program 81.2%
Taylor expanded in lambda2 around 0
*-commutativeN/A
fp-cancel-sign-sub-invN/A
fp-cancel-sub-sign-invN/A
*-commutativeN/A
sin-neg-revN/A
sin-neg-revN/A
remove-double-negN/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-sin.f64N/A
lower-sin.f6460.5
Applied rewrites60.5%
if -1.29999999999999989e-4 < phi2 < 3.2000000000000002e-14Initial program 67.3%
lift-cos.f64N/A
lift--.f64N/A
cos-diffN/A
*-commutativeN/A
lower-fma.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-sin.f64N/A
lower-sin.f6493.9
Applied rewrites93.9%
Taylor expanded in phi2 around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-sin.f64N/A
lower-sin.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
lower-cos.f6493.8
Applied rewrites93.8%
if 3.2000000000000002e-14 < phi2 < 2.6999999999999999e81Initial program 80.6%
Taylor expanded in lambda1 around 0
*-commutativeN/A
fp-cancel-sign-sub-invN/A
fp-cancel-sub-sign-invN/A
*-commutativeN/A
sin-neg-revN/A
sin-neg-revN/A
remove-double-negN/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
cos-negN/A
lower-cos.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-sin.f64N/A
lower-sin.f6471.2
Applied rewrites71.2%
if 2.6999999999999999e81 < phi2 Initial program 82.8%
Taylor expanded in lambda2 around 0
*-commutativeN/A
fp-cancel-sign-sub-invN/A
fp-cancel-sub-sign-invN/A
*-commutativeN/A
sin-neg-revN/A
sin-neg-revN/A
remove-double-negN/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-sin.f64N/A
lower-sin.f6471.9
Applied rewrites71.9%
Applied rewrites71.9%
Final simplification81.2%
(FPCore (R lambda1 lambda2 phi1 phi2)
:precision binary64
(if (or (<= phi1 -0.095) (not (<= phi1 5.5e-134)))
(*
(acos
(fma
(* (cos (- lambda2 lambda1)) (cos phi1))
(cos phi2)
(* (sin phi2) (sin phi1))))
R)
(*
(acos
(*
(cos phi2)
(fma (cos lambda1) (cos lambda2) (* (sin lambda2) (sin lambda1)))))
R)))
double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
double tmp;
if ((phi1 <= -0.095) || !(phi1 <= 5.5e-134)) {
tmp = acos(fma((cos((lambda2 - lambda1)) * cos(phi1)), cos(phi2), (sin(phi2) * sin(phi1)))) * R;
} else {
tmp = acos((cos(phi2) * fma(cos(lambda1), cos(lambda2), (sin(lambda2) * sin(lambda1))))) * R;
}
return tmp;
}
function code(R, lambda1, lambda2, phi1, phi2) tmp = 0.0 if ((phi1 <= -0.095) || !(phi1 <= 5.5e-134)) tmp = Float64(acos(fma(Float64(cos(Float64(lambda2 - lambda1)) * cos(phi1)), cos(phi2), Float64(sin(phi2) * sin(phi1)))) * R); else tmp = Float64(acos(Float64(cos(phi2) * fma(cos(lambda1), cos(lambda2), Float64(sin(lambda2) * sin(lambda1))))) * R); end return tmp end
code[R_, lambda1_, lambda2_, phi1_, phi2_] := If[Or[LessEqual[phi1, -0.095], N[Not[LessEqual[phi1, 5.5e-134]], $MachinePrecision]], N[(N[ArcCos[N[(N[(N[Cos[N[(lambda2 - lambda1), $MachinePrecision]], $MachinePrecision] * N[Cos[phi1], $MachinePrecision]), $MachinePrecision] * N[Cos[phi2], $MachinePrecision] + N[(N[Sin[phi2], $MachinePrecision] * N[Sin[phi1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * R), $MachinePrecision], N[(N[ArcCos[N[(N[Cos[phi2], $MachinePrecision] * N[(N[Cos[lambda1], $MachinePrecision] * N[Cos[lambda2], $MachinePrecision] + N[(N[Sin[lambda2], $MachinePrecision] * N[Sin[lambda1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * R), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\phi_1 \leq -0.095 \lor \neg \left(\phi_1 \leq 5.5 \cdot 10^{-134}\right):\\
\;\;\;\;\cos^{-1} \left(\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right) \cdot \cos \phi_1, \cos \phi_2, \sin \phi_2 \cdot \sin \phi_1\right)\right) \cdot R\\
\mathbf{else}:\\
\;\;\;\;\cos^{-1} \left(\cos \phi_2 \cdot \mathsf{fma}\left(\cos \lambda_1, \cos \lambda_2, \sin \lambda_2 \cdot \sin \lambda_1\right)\right) \cdot R\\
\end{array}
\end{array}
if phi1 < -0.095000000000000001 or 5.5000000000000002e-134 < phi1 Initial program 79.4%
lift-+.f64N/A
+-commutativeN/A
lift-*.f64N/A
*-commutativeN/A
lift-*.f64N/A
associate-*r*N/A
lower-fma.f64N/A
lower-*.f6479.4
lift-cos.f64N/A
lift--.f64N/A
cos-diffN/A
*-commutativeN/A
*-commutativeN/A
cos-diff-revN/A
lower-cos.f64N/A
lower--.f6479.4
lift-*.f64N/A
*-commutativeN/A
lower-*.f6479.4
Applied rewrites79.4%
if -0.095000000000000001 < phi1 < 5.5000000000000002e-134Initial program 65.3%
lift-cos.f64N/A
lift--.f64N/A
cos-diffN/A
*-commutativeN/A
lower-fma.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-sin.f64N/A
lower-sin.f6492.9
Applied rewrites92.9%
lift-+.f64N/A
+-commutativeN/A
lift-*.f64N/A
lift-fma.f64N/A
lift-*.f64N/A
distribute-lft-inN/A
associate-+l+N/A
lift-*.f64N/A
associate-*r*N/A
Applied rewrites92.9%
Taylor expanded in phi1 around 0
associate-*r*N/A
*-commutativeN/A
distribute-rgt-outN/A
lower-*.f64N/A
lower-cos.f64N/A
lower-fma.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-sin.f64N/A
lower-sin.f6491.9
Applied rewrites91.9%
Final simplification84.3%
(FPCore (R lambda1 lambda2 phi1 phi2)
:precision binary64
(if (or (<= phi1 -0.095) (not (<= phi1 5.5e-134)))
(*
(acos
(fma
(sin phi2)
(sin phi1)
(* (* (cos (- lambda2 lambda1)) (cos phi1)) (cos phi2))))
R)
(*
(acos
(*
(cos phi2)
(fma (cos lambda1) (cos lambda2) (* (sin lambda2) (sin lambda1)))))
R)))
double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
double tmp;
if ((phi1 <= -0.095) || !(phi1 <= 5.5e-134)) {
tmp = acos(fma(sin(phi2), sin(phi1), ((cos((lambda2 - lambda1)) * cos(phi1)) * cos(phi2)))) * R;
} else {
tmp = acos((cos(phi2) * fma(cos(lambda1), cos(lambda2), (sin(lambda2) * sin(lambda1))))) * R;
}
return tmp;
}
function code(R, lambda1, lambda2, phi1, phi2) tmp = 0.0 if ((phi1 <= -0.095) || !(phi1 <= 5.5e-134)) tmp = Float64(acos(fma(sin(phi2), sin(phi1), Float64(Float64(cos(Float64(lambda2 - lambda1)) * cos(phi1)) * cos(phi2)))) * R); else tmp = Float64(acos(Float64(cos(phi2) * fma(cos(lambda1), cos(lambda2), Float64(sin(lambda2) * sin(lambda1))))) * R); end return tmp end
code[R_, lambda1_, lambda2_, phi1_, phi2_] := If[Or[LessEqual[phi1, -0.095], N[Not[LessEqual[phi1, 5.5e-134]], $MachinePrecision]], N[(N[ArcCos[N[(N[Sin[phi2], $MachinePrecision] * N[Sin[phi1], $MachinePrecision] + N[(N[(N[Cos[N[(lambda2 - lambda1), $MachinePrecision]], $MachinePrecision] * N[Cos[phi1], $MachinePrecision]), $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * R), $MachinePrecision], N[(N[ArcCos[N[(N[Cos[phi2], $MachinePrecision] * N[(N[Cos[lambda1], $MachinePrecision] * N[Cos[lambda2], $MachinePrecision] + N[(N[Sin[lambda2], $MachinePrecision] * N[Sin[lambda1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * R), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\phi_1 \leq -0.095 \lor \neg \left(\phi_1 \leq 5.5 \cdot 10^{-134}\right):\\
\;\;\;\;\cos^{-1} \left(\mathsf{fma}\left(\sin \phi_2, \sin \phi_1, \left(\cos \left(\lambda_2 - \lambda_1\right) \cdot \cos \phi_1\right) \cdot \cos \phi_2\right)\right) \cdot R\\
\mathbf{else}:\\
\;\;\;\;\cos^{-1} \left(\cos \phi_2 \cdot \mathsf{fma}\left(\cos \lambda_1, \cos \lambda_2, \sin \lambda_2 \cdot \sin \lambda_1\right)\right) \cdot R\\
\end{array}
\end{array}
if phi1 < -0.095000000000000001 or 5.5000000000000002e-134 < phi1 Initial program 79.4%
lift-+.f64N/A
lift-*.f64N/A
*-commutativeN/A
lower-fma.f6479.4
lift-*.f64N/A
*-commutativeN/A
lift-*.f64N/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f6479.4
lift-cos.f64N/A
lift--.f64N/A
cos-diffN/A
*-commutativeN/A
*-commutativeN/A
cos-diff-revN/A
lower-cos.f64N/A
lower--.f6479.4
Applied rewrites79.4%
if -0.095000000000000001 < phi1 < 5.5000000000000002e-134Initial program 65.3%
lift-cos.f64N/A
lift--.f64N/A
cos-diffN/A
*-commutativeN/A
lower-fma.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-sin.f64N/A
lower-sin.f6492.9
Applied rewrites92.9%
lift-+.f64N/A
+-commutativeN/A
lift-*.f64N/A
lift-fma.f64N/A
lift-*.f64N/A
distribute-lft-inN/A
associate-+l+N/A
lift-*.f64N/A
associate-*r*N/A
Applied rewrites92.9%
Taylor expanded in phi1 around 0
associate-*r*N/A
*-commutativeN/A
distribute-rgt-outN/A
lower-*.f64N/A
lower-cos.f64N/A
lower-fma.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-sin.f64N/A
lower-sin.f6491.9
Applied rewrites91.9%
Final simplification84.2%
(FPCore (R lambda1 lambda2 phi1 phi2)
:precision binary64
(if (<= phi1 -0.095)
(*
(acos
(+
(* (sin phi1) (sin phi2))
(* (* (cos phi1) (cos phi2)) (cos (- lambda1 lambda2)))))
R)
(if (<= phi1 5.5e-134)
(*
(acos
(*
(cos phi2)
(fma (cos lambda1) (cos lambda2) (* (sin lambda2) (sin lambda1)))))
R)
(*
(acos
(fma
(* (cos (- lambda2 lambda1)) (cos phi1))
(cos phi2)
(* (sin phi2) (sin phi1))))
R))))
double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
double tmp;
if (phi1 <= -0.095) {
tmp = acos(((sin(phi1) * sin(phi2)) + ((cos(phi1) * cos(phi2)) * cos((lambda1 - lambda2))))) * R;
} else if (phi1 <= 5.5e-134) {
tmp = acos((cos(phi2) * fma(cos(lambda1), cos(lambda2), (sin(lambda2) * sin(lambda1))))) * R;
} else {
tmp = acos(fma((cos((lambda2 - lambda1)) * cos(phi1)), cos(phi2), (sin(phi2) * sin(phi1)))) * R;
}
return tmp;
}
function code(R, lambda1, lambda2, phi1, phi2) tmp = 0.0 if (phi1 <= -0.095) tmp = Float64(acos(Float64(Float64(sin(phi1) * sin(phi2)) + Float64(Float64(cos(phi1) * cos(phi2)) * cos(Float64(lambda1 - lambda2))))) * R); elseif (phi1 <= 5.5e-134) tmp = Float64(acos(Float64(cos(phi2) * fma(cos(lambda1), cos(lambda2), Float64(sin(lambda2) * sin(lambda1))))) * R); else tmp = Float64(acos(fma(Float64(cos(Float64(lambda2 - lambda1)) * cos(phi1)), cos(phi2), Float64(sin(phi2) * sin(phi1)))) * R); end return tmp end
code[R_, lambda1_, lambda2_, phi1_, phi2_] := If[LessEqual[phi1, -0.095], N[(N[ArcCos[N[(N[(N[Sin[phi1], $MachinePrecision] * N[Sin[phi2], $MachinePrecision]), $MachinePrecision] + N[(N[(N[Cos[phi1], $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision] * N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * R), $MachinePrecision], If[LessEqual[phi1, 5.5e-134], N[(N[ArcCos[N[(N[Cos[phi2], $MachinePrecision] * N[(N[Cos[lambda1], $MachinePrecision] * N[Cos[lambda2], $MachinePrecision] + N[(N[Sin[lambda2], $MachinePrecision] * N[Sin[lambda1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * R), $MachinePrecision], N[(N[ArcCos[N[(N[(N[Cos[N[(lambda2 - lambda1), $MachinePrecision]], $MachinePrecision] * N[Cos[phi1], $MachinePrecision]), $MachinePrecision] * N[Cos[phi2], $MachinePrecision] + N[(N[Sin[phi2], $MachinePrecision] * N[Sin[phi1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * R), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\phi_1 \leq -0.095:\\
\;\;\;\;\cos^{-1} \left(\sin \phi_1 \cdot \sin \phi_2 + \left(\cos \phi_1 \cdot \cos \phi_2\right) \cdot \cos \left(\lambda_1 - \lambda_2\right)\right) \cdot R\\
\mathbf{elif}\;\phi_1 \leq 5.5 \cdot 10^{-134}:\\
\;\;\;\;\cos^{-1} \left(\cos \phi_2 \cdot \mathsf{fma}\left(\cos \lambda_1, \cos \lambda_2, \sin \lambda_2 \cdot \sin \lambda_1\right)\right) \cdot R\\
\mathbf{else}:\\
\;\;\;\;\cos^{-1} \left(\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right) \cdot \cos \phi_1, \cos \phi_2, \sin \phi_2 \cdot \sin \phi_1\right)\right) \cdot R\\
\end{array}
\end{array}
if phi1 < -0.095000000000000001Initial program 84.5%
if -0.095000000000000001 < phi1 < 5.5000000000000002e-134Initial program 65.3%
lift-cos.f64N/A
lift--.f64N/A
cos-diffN/A
*-commutativeN/A
lower-fma.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-sin.f64N/A
lower-sin.f6492.9
Applied rewrites92.9%
lift-+.f64N/A
+-commutativeN/A
lift-*.f64N/A
lift-fma.f64N/A
lift-*.f64N/A
distribute-lft-inN/A
associate-+l+N/A
lift-*.f64N/A
associate-*r*N/A
Applied rewrites92.9%
Taylor expanded in phi1 around 0
associate-*r*N/A
*-commutativeN/A
distribute-rgt-outN/A
lower-*.f64N/A
lower-cos.f64N/A
lower-fma.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-sin.f64N/A
lower-sin.f6491.9
Applied rewrites91.9%
if 5.5000000000000002e-134 < phi1 Initial program 75.9%
lift-+.f64N/A
+-commutativeN/A
lift-*.f64N/A
*-commutativeN/A
lift-*.f64N/A
associate-*r*N/A
lower-fma.f64N/A
lower-*.f6475.9
lift-cos.f64N/A
lift--.f64N/A
cos-diffN/A
*-commutativeN/A
*-commutativeN/A
cos-diff-revN/A
lower-cos.f64N/A
lower--.f6475.9
lift-*.f64N/A
*-commutativeN/A
lower-*.f6475.9
Applied rewrites75.9%
(FPCore (R lambda1 lambda2 phi1 phi2)
:precision binary64
(if (or (<= lambda1 -9.6e-8) (not (<= lambda1 3.8e-5)))
(*
(acos
(fma (sin phi2) (sin phi1) (* (cos lambda1) (* (cos phi2) (cos phi1)))))
R)
(*
(acos
(fma (* (cos lambda2) (cos phi2)) (cos phi1) (* (sin phi2) (sin phi1))))
R)))
double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
double tmp;
if ((lambda1 <= -9.6e-8) || !(lambda1 <= 3.8e-5)) {
tmp = acos(fma(sin(phi2), sin(phi1), (cos(lambda1) * (cos(phi2) * cos(phi1))))) * R;
} else {
tmp = acos(fma((cos(lambda2) * cos(phi2)), cos(phi1), (sin(phi2) * sin(phi1)))) * R;
}
return tmp;
}
function code(R, lambda1, lambda2, phi1, phi2) tmp = 0.0 if ((lambda1 <= -9.6e-8) || !(lambda1 <= 3.8e-5)) tmp = Float64(acos(fma(sin(phi2), sin(phi1), Float64(cos(lambda1) * Float64(cos(phi2) * cos(phi1))))) * R); else tmp = Float64(acos(fma(Float64(cos(lambda2) * cos(phi2)), cos(phi1), Float64(sin(phi2) * sin(phi1)))) * R); end return tmp end
code[R_, lambda1_, lambda2_, phi1_, phi2_] := If[Or[LessEqual[lambda1, -9.6e-8], N[Not[LessEqual[lambda1, 3.8e-5]], $MachinePrecision]], N[(N[ArcCos[N[(N[Sin[phi2], $MachinePrecision] * N[Sin[phi1], $MachinePrecision] + N[(N[Cos[lambda1], $MachinePrecision] * N[(N[Cos[phi2], $MachinePrecision] * N[Cos[phi1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * R), $MachinePrecision], N[(N[ArcCos[N[(N[(N[Cos[lambda2], $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision] * N[Cos[phi1], $MachinePrecision] + N[(N[Sin[phi2], $MachinePrecision] * N[Sin[phi1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * R), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\lambda_1 \leq -9.6 \cdot 10^{-8} \lor \neg \left(\lambda_1 \leq 3.8 \cdot 10^{-5}\right):\\
\;\;\;\;\cos^{-1} \left(\mathsf{fma}\left(\sin \phi_2, \sin \phi_1, \cos \lambda_1 \cdot \left(\cos \phi_2 \cdot \cos \phi_1\right)\right)\right) \cdot R\\
\mathbf{else}:\\
\;\;\;\;\cos^{-1} \left(\mathsf{fma}\left(\cos \lambda_2 \cdot \cos \phi_2, \cos \phi_1, \sin \phi_2 \cdot \sin \phi_1\right)\right) \cdot R\\
\end{array}
\end{array}
if lambda1 < -9.59999999999999994e-8 or 3.8000000000000002e-5 < lambda1 Initial program 56.3%
Taylor expanded in lambda2 around 0
*-commutativeN/A
fp-cancel-sign-sub-invN/A
fp-cancel-sub-sign-invN/A
*-commutativeN/A
sin-neg-revN/A
sin-neg-revN/A
remove-double-negN/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-sin.f64N/A
lower-sin.f6456.5
Applied rewrites56.5%
Applied rewrites56.5%
if -9.59999999999999994e-8 < lambda1 < 3.8000000000000002e-5Initial program 93.0%
Taylor expanded in lambda1 around 0
*-commutativeN/A
fp-cancel-sign-sub-invN/A
fp-cancel-sub-sign-invN/A
*-commutativeN/A
sin-neg-revN/A
sin-neg-revN/A
remove-double-negN/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
cos-negN/A
lower-cos.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-sin.f64N/A
lower-sin.f6493.0
Applied rewrites93.0%
Final simplification74.0%
(FPCore (R lambda1 lambda2 phi1 phi2)
:precision binary64
(if (or (<= lambda2 -5.4e-8) (not (<= lambda2 3e-24)))
(* (acos (+ 0.0 (* (* (cos phi1) (cos phi2)) (cos (- lambda1 lambda2))))) R)
(*
(acos
(fma (sin phi2) (sin phi1) (* (cos lambda1) (* (cos phi2) (cos phi1)))))
R)))
double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
double tmp;
if ((lambda2 <= -5.4e-8) || !(lambda2 <= 3e-24)) {
tmp = acos((0.0 + ((cos(phi1) * cos(phi2)) * cos((lambda1 - lambda2))))) * R;
} else {
tmp = acos(fma(sin(phi2), sin(phi1), (cos(lambda1) * (cos(phi2) * cos(phi1))))) * R;
}
return tmp;
}
function code(R, lambda1, lambda2, phi1, phi2) tmp = 0.0 if ((lambda2 <= -5.4e-8) || !(lambda2 <= 3e-24)) tmp = Float64(acos(Float64(0.0 + Float64(Float64(cos(phi1) * cos(phi2)) * cos(Float64(lambda1 - lambda2))))) * R); else tmp = Float64(acos(fma(sin(phi2), sin(phi1), Float64(cos(lambda1) * Float64(cos(phi2) * cos(phi1))))) * R); end return tmp end
code[R_, lambda1_, lambda2_, phi1_, phi2_] := If[Or[LessEqual[lambda2, -5.4e-8], N[Not[LessEqual[lambda2, 3e-24]], $MachinePrecision]], N[(N[ArcCos[N[(0.0 + N[(N[(N[Cos[phi1], $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision] * N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * R), $MachinePrecision], N[(N[ArcCos[N[(N[Sin[phi2], $MachinePrecision] * N[Sin[phi1], $MachinePrecision] + N[(N[Cos[lambda1], $MachinePrecision] * N[(N[Cos[phi2], $MachinePrecision] * N[Cos[phi1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * R), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\lambda_2 \leq -5.4 \cdot 10^{-8} \lor \neg \left(\lambda_2 \leq 3 \cdot 10^{-24}\right):\\
\;\;\;\;\cos^{-1} \left(0 + \left(\cos \phi_1 \cdot \cos \phi_2\right) \cdot \cos \left(\lambda_1 - \lambda_2\right)\right) \cdot R\\
\mathbf{else}:\\
\;\;\;\;\cos^{-1} \left(\mathsf{fma}\left(\sin \phi_2, \sin \phi_1, \cos \lambda_1 \cdot \left(\cos \phi_2 \cdot \cos \phi_1\right)\right)\right) \cdot R\\
\end{array}
\end{array}
if lambda2 < -5.40000000000000005e-8 or 2.99999999999999995e-24 < lambda2 Initial program 59.7%
lift-*.f64N/A
lift-sin.f64N/A
lift-sin.f64N/A
sin-multN/A
lower-/.f64N/A
lower--.f64N/A
cos-diffN/A
lift-cos.f64N/A
lift-cos.f64N/A
*-commutativeN/A
lift-cos.f64N/A
lift-cos.f64N/A
lift-sin.f64N/A
lift-sin.f64N/A
*-commutativeN/A
lift-sin.f64N/A
lift-sin.f64N/A
cos-diff-revN/A
lower-cos.f64N/A
lower--.f64N/A
lower-cos.f64N/A
+-commutativeN/A
lower-+.f6454.0
Applied rewrites54.0%
Taylor expanded in phi2 around 0
cos-neg-revN/A
+-inversesN/A
metadata-eval53.7
Applied rewrites53.7%
if -5.40000000000000005e-8 < lambda2 < 2.99999999999999995e-24Initial program 93.1%
Taylor expanded in lambda2 around 0
*-commutativeN/A
fp-cancel-sign-sub-invN/A
fp-cancel-sub-sign-invN/A
*-commutativeN/A
sin-neg-revN/A
sin-neg-revN/A
remove-double-negN/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-sin.f64N/A
lower-sin.f6493.1
Applied rewrites93.1%
Applied rewrites93.1%
Final simplification70.5%
(FPCore (R lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (* (sin phi2) (sin phi1))) (t_1 (* (cos phi2) (cos phi1))))
(if (<= lambda1 -9.6e-8)
(* (acos (fma t_1 (cos lambda1) t_0)) R)
(if (<= lambda1 3.8e-5)
(* (acos (fma (* (cos lambda2) (cos phi2)) (cos phi1) t_0)) R)
(* (acos (fma (sin phi2) (sin phi1) (* (cos lambda1) t_1))) R)))))
double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = sin(phi2) * sin(phi1);
double t_1 = cos(phi2) * cos(phi1);
double tmp;
if (lambda1 <= -9.6e-8) {
tmp = acos(fma(t_1, cos(lambda1), t_0)) * R;
} else if (lambda1 <= 3.8e-5) {
tmp = acos(fma((cos(lambda2) * cos(phi2)), cos(phi1), t_0)) * R;
} else {
tmp = acos(fma(sin(phi2), sin(phi1), (cos(lambda1) * t_1))) * R;
}
return tmp;
}
function code(R, lambda1, lambda2, phi1, phi2) t_0 = Float64(sin(phi2) * sin(phi1)) t_1 = Float64(cos(phi2) * cos(phi1)) tmp = 0.0 if (lambda1 <= -9.6e-8) tmp = Float64(acos(fma(t_1, cos(lambda1), t_0)) * R); elseif (lambda1 <= 3.8e-5) tmp = Float64(acos(fma(Float64(cos(lambda2) * cos(phi2)), cos(phi1), t_0)) * R); else tmp = Float64(acos(fma(sin(phi2), sin(phi1), Float64(cos(lambda1) * t_1))) * R); end return tmp end
code[R_, lambda1_, lambda2_, phi1_, phi2_] := Block[{t$95$0 = N[(N[Sin[phi2], $MachinePrecision] * N[Sin[phi1], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[Cos[phi2], $MachinePrecision] * N[Cos[phi1], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[lambda1, -9.6e-8], N[(N[ArcCos[N[(t$95$1 * N[Cos[lambda1], $MachinePrecision] + t$95$0), $MachinePrecision]], $MachinePrecision] * R), $MachinePrecision], If[LessEqual[lambda1, 3.8e-5], N[(N[ArcCos[N[(N[(N[Cos[lambda2], $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision] * N[Cos[phi1], $MachinePrecision] + t$95$0), $MachinePrecision]], $MachinePrecision] * R), $MachinePrecision], N[(N[ArcCos[N[(N[Sin[phi2], $MachinePrecision] * N[Sin[phi1], $MachinePrecision] + N[(N[Cos[lambda1], $MachinePrecision] * t$95$1), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * R), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sin \phi_2 \cdot \sin \phi_1\\
t_1 := \cos \phi_2 \cdot \cos \phi_1\\
\mathbf{if}\;\lambda_1 \leq -9.6 \cdot 10^{-8}:\\
\;\;\;\;\cos^{-1} \left(\mathsf{fma}\left(t\_1, \cos \lambda_1, t\_0\right)\right) \cdot R\\
\mathbf{elif}\;\lambda_1 \leq 3.8 \cdot 10^{-5}:\\
\;\;\;\;\cos^{-1} \left(\mathsf{fma}\left(\cos \lambda_2 \cdot \cos \phi_2, \cos \phi_1, t\_0\right)\right) \cdot R\\
\mathbf{else}:\\
\;\;\;\;\cos^{-1} \left(\mathsf{fma}\left(\sin \phi_2, \sin \phi_1, \cos \lambda_1 \cdot t\_1\right)\right) \cdot R\\
\end{array}
\end{array}
if lambda1 < -9.59999999999999994e-8Initial program 58.5%
Taylor expanded in lambda2 around 0
*-commutativeN/A
fp-cancel-sign-sub-invN/A
fp-cancel-sub-sign-invN/A
*-commutativeN/A
sin-neg-revN/A
sin-neg-revN/A
remove-double-negN/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-sin.f64N/A
lower-sin.f6458.5
Applied rewrites58.5%
if -9.59999999999999994e-8 < lambda1 < 3.8000000000000002e-5Initial program 93.0%
Taylor expanded in lambda1 around 0
*-commutativeN/A
fp-cancel-sign-sub-invN/A
fp-cancel-sub-sign-invN/A
*-commutativeN/A
sin-neg-revN/A
sin-neg-revN/A
remove-double-negN/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
cos-negN/A
lower-cos.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-sin.f64N/A
lower-sin.f6493.0
Applied rewrites93.0%
if 3.8000000000000002e-5 < lambda1 Initial program 54.4%
Taylor expanded in lambda2 around 0
*-commutativeN/A
fp-cancel-sign-sub-invN/A
fp-cancel-sub-sign-invN/A
*-commutativeN/A
sin-neg-revN/A
sin-neg-revN/A
remove-double-negN/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-sin.f64N/A
lower-sin.f6454.8
Applied rewrites54.8%
Applied rewrites54.8%
(FPCore (R lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (cos (- lambda1 lambda2))))
(if (<= t_0 0.9999999999999986)
(* (acos (+ 0.0 (* (* (cos phi1) (cos phi2)) t_0))) R)
(* (acos (fma (cos phi2) (cos phi1) (* (sin phi2) (sin phi1)))) R))))
double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = cos((lambda1 - lambda2));
double tmp;
if (t_0 <= 0.9999999999999986) {
tmp = acos((0.0 + ((cos(phi1) * cos(phi2)) * t_0))) * R;
} else {
tmp = acos(fma(cos(phi2), cos(phi1), (sin(phi2) * sin(phi1)))) * R;
}
return tmp;
}
function code(R, lambda1, lambda2, phi1, phi2) t_0 = cos(Float64(lambda1 - lambda2)) tmp = 0.0 if (t_0 <= 0.9999999999999986) tmp = Float64(acos(Float64(0.0 + Float64(Float64(cos(phi1) * cos(phi2)) * t_0))) * R); else tmp = Float64(acos(fma(cos(phi2), cos(phi1), Float64(sin(phi2) * sin(phi1)))) * R); end return tmp end
code[R_, lambda1_, lambda2_, phi1_, phi2_] := Block[{t$95$0 = N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]}, If[LessEqual[t$95$0, 0.9999999999999986], N[(N[ArcCos[N[(0.0 + N[(N[(N[Cos[phi1], $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision] * t$95$0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * R), $MachinePrecision], N[(N[ArcCos[N[(N[Cos[phi2], $MachinePrecision] * N[Cos[phi1], $MachinePrecision] + N[(N[Sin[phi2], $MachinePrecision] * N[Sin[phi1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * R), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \cos \left(\lambda_1 - \lambda_2\right)\\
\mathbf{if}\;t\_0 \leq 0.9999999999999986:\\
\;\;\;\;\cos^{-1} \left(0 + \left(\cos \phi_1 \cdot \cos \phi_2\right) \cdot t\_0\right) \cdot R\\
\mathbf{else}:\\
\;\;\;\;\cos^{-1} \left(\mathsf{fma}\left(\cos \phi_2, \cos \phi_1, \sin \phi_2 \cdot \sin \phi_1\right)\right) \cdot R\\
\end{array}
\end{array}
if (cos.f64 (-.f64 lambda1 lambda2)) < 0.999999999999998557Initial program 70.6%
lift-*.f64N/A
lift-sin.f64N/A
lift-sin.f64N/A
sin-multN/A
lower-/.f64N/A
lower--.f64N/A
cos-diffN/A
lift-cos.f64N/A
lift-cos.f64N/A
*-commutativeN/A
lift-cos.f64N/A
lift-cos.f64N/A
lift-sin.f64N/A
lift-sin.f64N/A
*-commutativeN/A
lift-sin.f64N/A
lift-sin.f64N/A
cos-diff-revN/A
lower-cos.f64N/A
lower--.f64N/A
lower-cos.f64N/A
+-commutativeN/A
lower-+.f6458.1
Applied rewrites58.1%
Taylor expanded in phi2 around 0
cos-neg-revN/A
+-inversesN/A
metadata-eval57.9
Applied rewrites57.9%
if 0.999999999999998557 < (cos.f64 (-.f64 lambda1 lambda2)) Initial program 87.5%
Taylor expanded in lambda2 around 0
*-commutativeN/A
fp-cancel-sign-sub-invN/A
fp-cancel-sub-sign-invN/A
*-commutativeN/A
sin-neg-revN/A
sin-neg-revN/A
remove-double-negN/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-sin.f64N/A
lower-sin.f6487.5
Applied rewrites87.5%
Taylor expanded in lambda1 around 0
Applied rewrites87.5%
(FPCore (R lambda1 lambda2 phi1 phi2) :precision binary64 (* (acos (+ 0.0 (* (* (cos phi1) (cos phi2)) (cos (- lambda1 lambda2))))) R))
double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
return acos((0.0 + ((cos(phi1) * cos(phi2)) * cos((lambda1 - lambda2))))) * R;
}
real(8) function code(r, lambda1, lambda2, phi1, phi2)
real(8), intent (in) :: r
real(8), intent (in) :: lambda1
real(8), intent (in) :: lambda2
real(8), intent (in) :: phi1
real(8), intent (in) :: phi2
code = acos((0.0d0 + ((cos(phi1) * cos(phi2)) * cos((lambda1 - lambda2))))) * r
end function
public static double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
return Math.acos((0.0 + ((Math.cos(phi1) * Math.cos(phi2)) * Math.cos((lambda1 - lambda2))))) * R;
}
def code(R, lambda1, lambda2, phi1, phi2): return math.acos((0.0 + ((math.cos(phi1) * math.cos(phi2)) * math.cos((lambda1 - lambda2))))) * R
function code(R, lambda1, lambda2, phi1, phi2) return Float64(acos(Float64(0.0 + Float64(Float64(cos(phi1) * cos(phi2)) * cos(Float64(lambda1 - lambda2))))) * R) end
function tmp = code(R, lambda1, lambda2, phi1, phi2) tmp = acos((0.0 + ((cos(phi1) * cos(phi2)) * cos((lambda1 - lambda2))))) * R; end
code[R_, lambda1_, lambda2_, phi1_, phi2_] := N[(N[ArcCos[N[(0.0 + N[(N[(N[Cos[phi1], $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision] * N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * R), $MachinePrecision]
\begin{array}{l}
\\
\cos^{-1} \left(0 + \left(\cos \phi_1 \cdot \cos \phi_2\right) \cdot \cos \left(\lambda_1 - \lambda_2\right)\right) \cdot R
\end{array}
Initial program 73.9%
lift-*.f64N/A
lift-sin.f64N/A
lift-sin.f64N/A
sin-multN/A
lower-/.f64N/A
lower--.f64N/A
cos-diffN/A
lift-cos.f64N/A
lift-cos.f64N/A
*-commutativeN/A
lift-cos.f64N/A
lift-cos.f64N/A
lift-sin.f64N/A
lift-sin.f64N/A
*-commutativeN/A
lift-sin.f64N/A
lift-sin.f64N/A
cos-diff-revN/A
lower-cos.f64N/A
lower--.f64N/A
lower-cos.f64N/A
+-commutativeN/A
lower-+.f6457.5
Applied rewrites57.5%
Taylor expanded in phi2 around 0
cos-neg-revN/A
+-inversesN/A
metadata-eval57.4
Applied rewrites57.4%
(FPCore (R lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (* (acos (* (cos lambda2) (cos phi1))) R)))
(if (<= lambda2 -5.4e-8)
t_0
(if (<= lambda2 -5.5e-122)
(* (acos (* (cos lambda1) (cos phi1))) R)
(if (<= lambda2 2.9e-51)
(* (acos (* (cos lambda1) (cos phi2))) R)
t_0)))))
double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = acos((cos(lambda2) * cos(phi1))) * R;
double tmp;
if (lambda2 <= -5.4e-8) {
tmp = t_0;
} else if (lambda2 <= -5.5e-122) {
tmp = acos((cos(lambda1) * cos(phi1))) * R;
} else if (lambda2 <= 2.9e-51) {
tmp = acos((cos(lambda1) * cos(phi2))) * R;
} else {
tmp = t_0;
}
return tmp;
}
real(8) function code(r, lambda1, lambda2, phi1, phi2)
real(8), intent (in) :: r
real(8), intent (in) :: lambda1
real(8), intent (in) :: lambda2
real(8), intent (in) :: phi1
real(8), intent (in) :: phi2
real(8) :: t_0
real(8) :: tmp
t_0 = acos((cos(lambda2) * cos(phi1))) * r
if (lambda2 <= (-5.4d-8)) then
tmp = t_0
else if (lambda2 <= (-5.5d-122)) then
tmp = acos((cos(lambda1) * cos(phi1))) * r
else if (lambda2 <= 2.9d-51) then
tmp = acos((cos(lambda1) * cos(phi2))) * r
else
tmp = t_0
end if
code = tmp
end function
public static double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = Math.acos((Math.cos(lambda2) * Math.cos(phi1))) * R;
double tmp;
if (lambda2 <= -5.4e-8) {
tmp = t_0;
} else if (lambda2 <= -5.5e-122) {
tmp = Math.acos((Math.cos(lambda1) * Math.cos(phi1))) * R;
} else if (lambda2 <= 2.9e-51) {
tmp = Math.acos((Math.cos(lambda1) * Math.cos(phi2))) * R;
} else {
tmp = t_0;
}
return tmp;
}
def code(R, lambda1, lambda2, phi1, phi2): t_0 = math.acos((math.cos(lambda2) * math.cos(phi1))) * R tmp = 0 if lambda2 <= -5.4e-8: tmp = t_0 elif lambda2 <= -5.5e-122: tmp = math.acos((math.cos(lambda1) * math.cos(phi1))) * R elif lambda2 <= 2.9e-51: tmp = math.acos((math.cos(lambda1) * math.cos(phi2))) * R else: tmp = t_0 return tmp
function code(R, lambda1, lambda2, phi1, phi2) t_0 = Float64(acos(Float64(cos(lambda2) * cos(phi1))) * R) tmp = 0.0 if (lambda2 <= -5.4e-8) tmp = t_0; elseif (lambda2 <= -5.5e-122) tmp = Float64(acos(Float64(cos(lambda1) * cos(phi1))) * R); elseif (lambda2 <= 2.9e-51) tmp = Float64(acos(Float64(cos(lambda1) * cos(phi2))) * R); else tmp = t_0; end return tmp end
function tmp_2 = code(R, lambda1, lambda2, phi1, phi2) t_0 = acos((cos(lambda2) * cos(phi1))) * R; tmp = 0.0; if (lambda2 <= -5.4e-8) tmp = t_0; elseif (lambda2 <= -5.5e-122) tmp = acos((cos(lambda1) * cos(phi1))) * R; elseif (lambda2 <= 2.9e-51) tmp = acos((cos(lambda1) * cos(phi2))) * R; else tmp = t_0; end tmp_2 = tmp; end
code[R_, lambda1_, lambda2_, phi1_, phi2_] := Block[{t$95$0 = N[(N[ArcCos[N[(N[Cos[lambda2], $MachinePrecision] * N[Cos[phi1], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * R), $MachinePrecision]}, If[LessEqual[lambda2, -5.4e-8], t$95$0, If[LessEqual[lambda2, -5.5e-122], N[(N[ArcCos[N[(N[Cos[lambda1], $MachinePrecision] * N[Cos[phi1], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * R), $MachinePrecision], If[LessEqual[lambda2, 2.9e-51], N[(N[ArcCos[N[(N[Cos[lambda1], $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * R), $MachinePrecision], t$95$0]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \cos^{-1} \left(\cos \lambda_2 \cdot \cos \phi_1\right) \cdot R\\
\mathbf{if}\;\lambda_2 \leq -5.4 \cdot 10^{-8}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;\lambda_2 \leq -5.5 \cdot 10^{-122}:\\
\;\;\;\;\cos^{-1} \left(\cos \lambda_1 \cdot \cos \phi_1\right) \cdot R\\
\mathbf{elif}\;\lambda_2 \leq 2.9 \cdot 10^{-51}:\\
\;\;\;\;\cos^{-1} \left(\cos \lambda_1 \cdot \cos \phi_2\right) \cdot R\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if lambda2 < -5.40000000000000005e-8 or 2.89999999999999973e-51 < lambda2 Initial program 60.4%
Taylor expanded in phi2 around 0
*-commutativeN/A
lower-*.f64N/A
*-lft-identityN/A
metadata-evalN/A
fp-cancel-sign-sub-invN/A
cos-neg-revN/A
remove-double-negN/A
mul-1-negN/A
distribute-neg-inN/A
+-commutativeN/A
mul-1-negN/A
lower-cos.f64N/A
remove-double-negN/A
fp-cancel-sign-sub-invN/A
metadata-evalN/A
*-lft-identityN/A
lower--.f64N/A
lower-cos.f6443.3
Applied rewrites43.3%
Taylor expanded in lambda1 around 0
Applied rewrites43.0%
if -5.40000000000000005e-8 < lambda2 < -5.50000000000000053e-122Initial program 93.8%
Taylor expanded in phi2 around 0
*-commutativeN/A
lower-*.f64N/A
*-lft-identityN/A
metadata-evalN/A
fp-cancel-sign-sub-invN/A
cos-neg-revN/A
remove-double-negN/A
mul-1-negN/A
distribute-neg-inN/A
+-commutativeN/A
mul-1-negN/A
lower-cos.f64N/A
remove-double-negN/A
fp-cancel-sign-sub-invN/A
metadata-evalN/A
*-lft-identityN/A
lower--.f64N/A
lower-cos.f6451.9
Applied rewrites51.9%
Taylor expanded in lambda2 around 0
Applied rewrites51.9%
if -5.50000000000000053e-122 < lambda2 < 2.89999999999999973e-51Initial program 93.7%
Taylor expanded in lambda2 around 0
*-commutativeN/A
fp-cancel-sign-sub-invN/A
fp-cancel-sub-sign-invN/A
*-commutativeN/A
sin-neg-revN/A
sin-neg-revN/A
remove-double-negN/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-sin.f64N/A
lower-sin.f6493.7
Applied rewrites93.7%
Taylor expanded in phi1 around 0
Applied rewrites45.5%
(FPCore (R lambda1 lambda2 phi1 phi2) :precision binary64 (if (or (<= lambda1 -7.8e-11) (not (<= lambda1 0.000108))) (* (acos (* (cos lambda1) (cos phi2))) R) (* (acos (* (cos lambda2) (cos phi1))) R)))
double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
double tmp;
if ((lambda1 <= -7.8e-11) || !(lambda1 <= 0.000108)) {
tmp = acos((cos(lambda1) * cos(phi2))) * R;
} else {
tmp = acos((cos(lambda2) * cos(phi1))) * R;
}
return tmp;
}
real(8) function code(r, lambda1, lambda2, phi1, phi2)
real(8), intent (in) :: r
real(8), intent (in) :: lambda1
real(8), intent (in) :: lambda2
real(8), intent (in) :: phi1
real(8), intent (in) :: phi2
real(8) :: tmp
if ((lambda1 <= (-7.8d-11)) .or. (.not. (lambda1 <= 0.000108d0))) then
tmp = acos((cos(lambda1) * cos(phi2))) * r
else
tmp = acos((cos(lambda2) * cos(phi1))) * r
end if
code = tmp
end function
public static double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
double tmp;
if ((lambda1 <= -7.8e-11) || !(lambda1 <= 0.000108)) {
tmp = Math.acos((Math.cos(lambda1) * Math.cos(phi2))) * R;
} else {
tmp = Math.acos((Math.cos(lambda2) * Math.cos(phi1))) * R;
}
return tmp;
}
def code(R, lambda1, lambda2, phi1, phi2): tmp = 0 if (lambda1 <= -7.8e-11) or not (lambda1 <= 0.000108): tmp = math.acos((math.cos(lambda1) * math.cos(phi2))) * R else: tmp = math.acos((math.cos(lambda2) * math.cos(phi1))) * R return tmp
function code(R, lambda1, lambda2, phi1, phi2) tmp = 0.0 if ((lambda1 <= -7.8e-11) || !(lambda1 <= 0.000108)) tmp = Float64(acos(Float64(cos(lambda1) * cos(phi2))) * R); else tmp = Float64(acos(Float64(cos(lambda2) * cos(phi1))) * R); end return tmp end
function tmp_2 = code(R, lambda1, lambda2, phi1, phi2) tmp = 0.0; if ((lambda1 <= -7.8e-11) || ~((lambda1 <= 0.000108))) tmp = acos((cos(lambda1) * cos(phi2))) * R; else tmp = acos((cos(lambda2) * cos(phi1))) * R; end tmp_2 = tmp; end
code[R_, lambda1_, lambda2_, phi1_, phi2_] := If[Or[LessEqual[lambda1, -7.8e-11], N[Not[LessEqual[lambda1, 0.000108]], $MachinePrecision]], N[(N[ArcCos[N[(N[Cos[lambda1], $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * R), $MachinePrecision], N[(N[ArcCos[N[(N[Cos[lambda2], $MachinePrecision] * N[Cos[phi1], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * R), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\lambda_1 \leq -7.8 \cdot 10^{-11} \lor \neg \left(\lambda_1 \leq 0.000108\right):\\
\;\;\;\;\cos^{-1} \left(\cos \lambda_1 \cdot \cos \phi_2\right) \cdot R\\
\mathbf{else}:\\
\;\;\;\;\cos^{-1} \left(\cos \lambda_2 \cdot \cos \phi_1\right) \cdot R\\
\end{array}
\end{array}
if lambda1 < -7.80000000000000021e-11 or 1.08e-4 < lambda1 Initial program 56.5%
Taylor expanded in lambda2 around 0
*-commutativeN/A
fp-cancel-sign-sub-invN/A
fp-cancel-sub-sign-invN/A
*-commutativeN/A
sin-neg-revN/A
sin-neg-revN/A
remove-double-negN/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-sin.f64N/A
lower-sin.f6456.2
Applied rewrites56.2%
Taylor expanded in phi1 around 0
Applied rewrites33.8%
if -7.80000000000000021e-11 < lambda1 < 1.08e-4Initial program 93.1%
Taylor expanded in phi2 around 0
*-commutativeN/A
lower-*.f64N/A
*-lft-identityN/A
metadata-evalN/A
fp-cancel-sign-sub-invN/A
cos-neg-revN/A
remove-double-negN/A
mul-1-negN/A
distribute-neg-inN/A
+-commutativeN/A
mul-1-negN/A
lower-cos.f64N/A
remove-double-negN/A
fp-cancel-sign-sub-invN/A
metadata-evalN/A
*-lft-identityN/A
lower--.f64N/A
lower-cos.f6451.7
Applied rewrites51.7%
Taylor expanded in lambda1 around 0
Applied rewrites51.7%
Final simplification42.4%
(FPCore (R lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (cos (- lambda2 lambda1))))
(if (<= phi1 -1.3e-6)
(* (acos (* t_0 (cos phi1))) R)
(* (acos (* t_0 (cos phi2))) R))))
double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = cos((lambda2 - lambda1));
double tmp;
if (phi1 <= -1.3e-6) {
tmp = acos((t_0 * cos(phi1))) * R;
} else {
tmp = acos((t_0 * cos(phi2))) * R;
}
return tmp;
}
real(8) function code(r, lambda1, lambda2, phi1, phi2)
real(8), intent (in) :: r
real(8), intent (in) :: lambda1
real(8), intent (in) :: lambda2
real(8), intent (in) :: phi1
real(8), intent (in) :: phi2
real(8) :: t_0
real(8) :: tmp
t_0 = cos((lambda2 - lambda1))
if (phi1 <= (-1.3d-6)) then
tmp = acos((t_0 * cos(phi1))) * r
else
tmp = acos((t_0 * cos(phi2))) * r
end if
code = tmp
end function
public static double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = Math.cos((lambda2 - lambda1));
double tmp;
if (phi1 <= -1.3e-6) {
tmp = Math.acos((t_0 * Math.cos(phi1))) * R;
} else {
tmp = Math.acos((t_0 * Math.cos(phi2))) * R;
}
return tmp;
}
def code(R, lambda1, lambda2, phi1, phi2): t_0 = math.cos((lambda2 - lambda1)) tmp = 0 if phi1 <= -1.3e-6: tmp = math.acos((t_0 * math.cos(phi1))) * R else: tmp = math.acos((t_0 * math.cos(phi2))) * R return tmp
function code(R, lambda1, lambda2, phi1, phi2) t_0 = cos(Float64(lambda2 - lambda1)) tmp = 0.0 if (phi1 <= -1.3e-6) tmp = Float64(acos(Float64(t_0 * cos(phi1))) * R); else tmp = Float64(acos(Float64(t_0 * cos(phi2))) * R); end return tmp end
function tmp_2 = code(R, lambda1, lambda2, phi1, phi2) t_0 = cos((lambda2 - lambda1)); tmp = 0.0; if (phi1 <= -1.3e-6) tmp = acos((t_0 * cos(phi1))) * R; else tmp = acos((t_0 * cos(phi2))) * R; end tmp_2 = tmp; end
code[R_, lambda1_, lambda2_, phi1_, phi2_] := Block[{t$95$0 = N[Cos[N[(lambda2 - lambda1), $MachinePrecision]], $MachinePrecision]}, If[LessEqual[phi1, -1.3e-6], N[(N[ArcCos[N[(t$95$0 * N[Cos[phi1], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * R), $MachinePrecision], N[(N[ArcCos[N[(t$95$0 * N[Cos[phi2], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * R), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \cos \left(\lambda_2 - \lambda_1\right)\\
\mathbf{if}\;\phi_1 \leq -1.3 \cdot 10^{-6}:\\
\;\;\;\;\cos^{-1} \left(t\_0 \cdot \cos \phi_1\right) \cdot R\\
\mathbf{else}:\\
\;\;\;\;\cos^{-1} \left(t\_0 \cdot \cos \phi_2\right) \cdot R\\
\end{array}
\end{array}
if phi1 < -1.30000000000000005e-6Initial program 82.5%
Taylor expanded in phi2 around 0
*-commutativeN/A
lower-*.f64N/A
*-lft-identityN/A
metadata-evalN/A
fp-cancel-sign-sub-invN/A
cos-neg-revN/A
remove-double-negN/A
mul-1-negN/A
distribute-neg-inN/A
+-commutativeN/A
mul-1-negN/A
lower-cos.f64N/A
remove-double-negN/A
fp-cancel-sign-sub-invN/A
metadata-evalN/A
*-lft-identityN/A
lower--.f64N/A
lower-cos.f6448.8
Applied rewrites48.8%
if -1.30000000000000005e-6 < phi1 Initial program 71.0%
Taylor expanded in phi1 around 0
*-commutativeN/A
lower-*.f64N/A
*-lft-identityN/A
metadata-evalN/A
fp-cancel-sign-sub-invN/A
cos-neg-revN/A
remove-double-negN/A
mul-1-negN/A
distribute-neg-inN/A
+-commutativeN/A
mul-1-negN/A
lower-cos.f64N/A
remove-double-negN/A
fp-cancel-sign-sub-invN/A
metadata-evalN/A
*-lft-identityN/A
lower--.f64N/A
lower-cos.f6450.5
Applied rewrites50.5%
(FPCore (R lambda1 lambda2 phi1 phi2) :precision binary64 (if (<= phi2 5.5e-5) (* (acos (* (cos (- lambda2 lambda1)) (cos phi1))) R) (* (acos (* (cos lambda1) (cos phi2))) R)))
double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
double tmp;
if (phi2 <= 5.5e-5) {
tmp = acos((cos((lambda2 - lambda1)) * cos(phi1))) * R;
} else {
tmp = acos((cos(lambda1) * cos(phi2))) * R;
}
return tmp;
}
real(8) function code(r, lambda1, lambda2, phi1, phi2)
real(8), intent (in) :: r
real(8), intent (in) :: lambda1
real(8), intent (in) :: lambda2
real(8), intent (in) :: phi1
real(8), intent (in) :: phi2
real(8) :: tmp
if (phi2 <= 5.5d-5) then
tmp = acos((cos((lambda2 - lambda1)) * cos(phi1))) * r
else
tmp = acos((cos(lambda1) * cos(phi2))) * r
end if
code = tmp
end function
public static double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
double tmp;
if (phi2 <= 5.5e-5) {
tmp = Math.acos((Math.cos((lambda2 - lambda1)) * Math.cos(phi1))) * R;
} else {
tmp = Math.acos((Math.cos(lambda1) * Math.cos(phi2))) * R;
}
return tmp;
}
def code(R, lambda1, lambda2, phi1, phi2): tmp = 0 if phi2 <= 5.5e-5: tmp = math.acos((math.cos((lambda2 - lambda1)) * math.cos(phi1))) * R else: tmp = math.acos((math.cos(lambda1) * math.cos(phi2))) * R return tmp
function code(R, lambda1, lambda2, phi1, phi2) tmp = 0.0 if (phi2 <= 5.5e-5) tmp = Float64(acos(Float64(cos(Float64(lambda2 - lambda1)) * cos(phi1))) * R); else tmp = Float64(acos(Float64(cos(lambda1) * cos(phi2))) * R); end return tmp end
function tmp_2 = code(R, lambda1, lambda2, phi1, phi2) tmp = 0.0; if (phi2 <= 5.5e-5) tmp = acos((cos((lambda2 - lambda1)) * cos(phi1))) * R; else tmp = acos((cos(lambda1) * cos(phi2))) * R; end tmp_2 = tmp; end
code[R_, lambda1_, lambda2_, phi1_, phi2_] := If[LessEqual[phi2, 5.5e-5], N[(N[ArcCos[N[(N[Cos[N[(lambda2 - lambda1), $MachinePrecision]], $MachinePrecision] * N[Cos[phi1], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * R), $MachinePrecision], N[(N[ArcCos[N[(N[Cos[lambda1], $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * R), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\phi_2 \leq 5.5 \cdot 10^{-5}:\\
\;\;\;\;\cos^{-1} \left(\cos \left(\lambda_2 - \lambda_1\right) \cdot \cos \phi_1\right) \cdot R\\
\mathbf{else}:\\
\;\;\;\;\cos^{-1} \left(\cos \lambda_1 \cdot \cos \phi_2\right) \cdot R\\
\end{array}
\end{array}
if phi2 < 5.5000000000000002e-5Initial program 71.4%
Taylor expanded in phi2 around 0
*-commutativeN/A
lower-*.f64N/A
*-lft-identityN/A
metadata-evalN/A
fp-cancel-sign-sub-invN/A
cos-neg-revN/A
remove-double-negN/A
mul-1-negN/A
distribute-neg-inN/A
+-commutativeN/A
mul-1-negN/A
lower-cos.f64N/A
remove-double-negN/A
fp-cancel-sign-sub-invN/A
metadata-evalN/A
*-lft-identityN/A
lower--.f64N/A
lower-cos.f6453.3
Applied rewrites53.3%
if 5.5000000000000002e-5 < phi2 Initial program 82.0%
Taylor expanded in lambda2 around 0
*-commutativeN/A
fp-cancel-sign-sub-invN/A
fp-cancel-sub-sign-invN/A
*-commutativeN/A
sin-neg-revN/A
sin-neg-revN/A
remove-double-negN/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-sin.f64N/A
lower-sin.f6466.2
Applied rewrites66.2%
Taylor expanded in phi1 around 0
Applied rewrites33.0%
(FPCore (R lambda1 lambda2 phi1 phi2) :precision binary64 (if (<= phi1 -3.6e-48) (* (acos (* (cos lambda2) (cos phi1))) R) (* (acos (* (fma -0.5 (* phi1 phi1) 1.0) (cos (- lambda1 lambda2)))) R)))
double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
double tmp;
if (phi1 <= -3.6e-48) {
tmp = acos((cos(lambda2) * cos(phi1))) * R;
} else {
tmp = acos((fma(-0.5, (phi1 * phi1), 1.0) * cos((lambda1 - lambda2)))) * R;
}
return tmp;
}
function code(R, lambda1, lambda2, phi1, phi2) tmp = 0.0 if (phi1 <= -3.6e-48) tmp = Float64(acos(Float64(cos(lambda2) * cos(phi1))) * R); else tmp = Float64(acos(Float64(fma(-0.5, Float64(phi1 * phi1), 1.0) * cos(Float64(lambda1 - lambda2)))) * R); end return tmp end
code[R_, lambda1_, lambda2_, phi1_, phi2_] := If[LessEqual[phi1, -3.6e-48], N[(N[ArcCos[N[(N[Cos[lambda2], $MachinePrecision] * N[Cos[phi1], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * R), $MachinePrecision], N[(N[ArcCos[N[(N[(-0.5 * N[(phi1 * phi1), $MachinePrecision] + 1.0), $MachinePrecision] * N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * R), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\phi_1 \leq -3.6 \cdot 10^{-48}:\\
\;\;\;\;\cos^{-1} \left(\cos \lambda_2 \cdot \cos \phi_1\right) \cdot R\\
\mathbf{else}:\\
\;\;\;\;\cos^{-1} \left(\mathsf{fma}\left(-0.5, \phi_1 \cdot \phi_1, 1\right) \cdot \cos \left(\lambda_1 - \lambda_2\right)\right) \cdot R\\
\end{array}
\end{array}
if phi1 < -3.6000000000000002e-48Initial program 82.2%
Taylor expanded in phi2 around 0
*-commutativeN/A
lower-*.f64N/A
*-lft-identityN/A
metadata-evalN/A
fp-cancel-sign-sub-invN/A
cos-neg-revN/A
remove-double-negN/A
mul-1-negN/A
distribute-neg-inN/A
+-commutativeN/A
mul-1-negN/A
lower-cos.f64N/A
remove-double-negN/A
fp-cancel-sign-sub-invN/A
metadata-evalN/A
*-lft-identityN/A
lower--.f64N/A
lower-cos.f6448.2
Applied rewrites48.2%
Taylor expanded in lambda1 around 0
Applied rewrites41.5%
if -3.6000000000000002e-48 < phi1 Initial program 70.6%
Taylor expanded in phi2 around 0
*-commutativeN/A
lower-*.f64N/A
*-lft-identityN/A
metadata-evalN/A
fp-cancel-sign-sub-invN/A
cos-neg-revN/A
remove-double-negN/A
mul-1-negN/A
distribute-neg-inN/A
+-commutativeN/A
mul-1-negN/A
lower-cos.f64N/A
remove-double-negN/A
fp-cancel-sign-sub-invN/A
metadata-evalN/A
*-lft-identityN/A
lower--.f64N/A
lower-cos.f6443.7
Applied rewrites43.7%
Taylor expanded in phi1 around 0
Applied rewrites27.1%
(FPCore (R lambda1 lambda2 phi1 phi2)
:precision binary64
(let* ((t_0 (fma -0.5 (* phi1 phi1) 1.0)))
(if (<= phi1 -1.06)
(* (acos (cos phi1)) R)
(if (<= phi1 -1.35e-283)
(* (acos (* t_0 (cos lambda2))) R)
(if (<= phi1 5.4e-141)
(* (acos (* t_0 (cos lambda1))) R)
(* (acos (cos (- phi2 phi1))) R))))))
double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
double t_0 = fma(-0.5, (phi1 * phi1), 1.0);
double tmp;
if (phi1 <= -1.06) {
tmp = acos(cos(phi1)) * R;
} else if (phi1 <= -1.35e-283) {
tmp = acos((t_0 * cos(lambda2))) * R;
} else if (phi1 <= 5.4e-141) {
tmp = acos((t_0 * cos(lambda1))) * R;
} else {
tmp = acos(cos((phi2 - phi1))) * R;
}
return tmp;
}
function code(R, lambda1, lambda2, phi1, phi2) t_0 = fma(-0.5, Float64(phi1 * phi1), 1.0) tmp = 0.0 if (phi1 <= -1.06) tmp = Float64(acos(cos(phi1)) * R); elseif (phi1 <= -1.35e-283) tmp = Float64(acos(Float64(t_0 * cos(lambda2))) * R); elseif (phi1 <= 5.4e-141) tmp = Float64(acos(Float64(t_0 * cos(lambda1))) * R); else tmp = Float64(acos(cos(Float64(phi2 - phi1))) * R); end return tmp end
code[R_, lambda1_, lambda2_, phi1_, phi2_] := Block[{t$95$0 = N[(-0.5 * N[(phi1 * phi1), $MachinePrecision] + 1.0), $MachinePrecision]}, If[LessEqual[phi1, -1.06], N[(N[ArcCos[N[Cos[phi1], $MachinePrecision]], $MachinePrecision] * R), $MachinePrecision], If[LessEqual[phi1, -1.35e-283], N[(N[ArcCos[N[(t$95$0 * N[Cos[lambda2], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * R), $MachinePrecision], If[LessEqual[phi1, 5.4e-141], N[(N[ArcCos[N[(t$95$0 * N[Cos[lambda1], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * R), $MachinePrecision], N[(N[ArcCos[N[Cos[N[(phi2 - phi1), $MachinePrecision]], $MachinePrecision]], $MachinePrecision] * R), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \mathsf{fma}\left(-0.5, \phi_1 \cdot \phi_1, 1\right)\\
\mathbf{if}\;\phi_1 \leq -1.06:\\
\;\;\;\;\cos^{-1} \cos \phi_1 \cdot R\\
\mathbf{elif}\;\phi_1 \leq -1.35 \cdot 10^{-283}:\\
\;\;\;\;\cos^{-1} \left(t\_0 \cdot \cos \lambda_2\right) \cdot R\\
\mathbf{elif}\;\phi_1 \leq 5.4 \cdot 10^{-141}:\\
\;\;\;\;\cos^{-1} \left(t\_0 \cdot \cos \lambda_1\right) \cdot R\\
\mathbf{else}:\\
\;\;\;\;\cos^{-1} \cos \left(\phi_2 - \phi_1\right) \cdot R\\
\end{array}
\end{array}
if phi1 < -1.0600000000000001Initial program 84.5%
Taylor expanded in lambda2 around 0
*-commutativeN/A
fp-cancel-sign-sub-invN/A
fp-cancel-sub-sign-invN/A
*-commutativeN/A
sin-neg-revN/A
sin-neg-revN/A
remove-double-negN/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-sin.f64N/A
lower-sin.f6466.9
Applied rewrites66.9%
Taylor expanded in lambda1 around 0
Applied rewrites45.8%
Applied rewrites27.7%
Taylor expanded in phi2 around 0
Applied rewrites27.8%
if -1.0600000000000001 < phi1 < -1.35e-283Initial program 73.6%
Taylor expanded in phi2 around 0
*-commutativeN/A
lower-*.f64N/A
*-lft-identityN/A
metadata-evalN/A
fp-cancel-sign-sub-invN/A
cos-neg-revN/A
remove-double-negN/A
mul-1-negN/A
distribute-neg-inN/A
+-commutativeN/A
mul-1-negN/A
lower-cos.f64N/A
remove-double-negN/A
fp-cancel-sign-sub-invN/A
metadata-evalN/A
*-lft-identityN/A
lower--.f64N/A
lower-cos.f6443.3
Applied rewrites43.3%
Taylor expanded in phi1 around 0
Applied rewrites43.3%
Taylor expanded in lambda1 around 0
Applied rewrites33.5%
if -1.35e-283 < phi1 < 5.4000000000000005e-141Initial program 55.5%
Taylor expanded in phi2 around 0
*-commutativeN/A
lower-*.f64N/A
*-lft-identityN/A
metadata-evalN/A
fp-cancel-sign-sub-invN/A
cos-neg-revN/A
remove-double-negN/A
mul-1-negN/A
distribute-neg-inN/A
+-commutativeN/A
mul-1-negN/A
lower-cos.f64N/A
remove-double-negN/A
fp-cancel-sign-sub-invN/A
metadata-evalN/A
*-lft-identityN/A
lower--.f64N/A
lower-cos.f6437.9
Applied rewrites37.9%
Taylor expanded in phi1 around 0
Applied rewrites37.9%
Taylor expanded in lambda2 around 0
Applied rewrites22.2%
if 5.4000000000000005e-141 < phi1 Initial program 75.6%
Taylor expanded in lambda2 around 0
*-commutativeN/A
fp-cancel-sign-sub-invN/A
fp-cancel-sub-sign-invN/A
*-commutativeN/A
sin-neg-revN/A
sin-neg-revN/A
remove-double-negN/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-sin.f64N/A
lower-sin.f6453.7
Applied rewrites53.7%
Taylor expanded in lambda1 around 0
Applied rewrites27.6%
Applied rewrites20.7%
(FPCore (R lambda1 lambda2 phi1 phi2) :precision binary64 (if (or (<= phi1 -6.8e-122) (not (<= phi1 5.4e-141))) (* (acos (cos (- phi2 phi1))) R) (* (acos (* (fma -0.5 (* phi1 phi1) 1.0) (cos lambda1))) R)))
double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
double tmp;
if ((phi1 <= -6.8e-122) || !(phi1 <= 5.4e-141)) {
tmp = acos(cos((phi2 - phi1))) * R;
} else {
tmp = acos((fma(-0.5, (phi1 * phi1), 1.0) * cos(lambda1))) * R;
}
return tmp;
}
function code(R, lambda1, lambda2, phi1, phi2) tmp = 0.0 if ((phi1 <= -6.8e-122) || !(phi1 <= 5.4e-141)) tmp = Float64(acos(cos(Float64(phi2 - phi1))) * R); else tmp = Float64(acos(Float64(fma(-0.5, Float64(phi1 * phi1), 1.0) * cos(lambda1))) * R); end return tmp end
code[R_, lambda1_, lambda2_, phi1_, phi2_] := If[Or[LessEqual[phi1, -6.8e-122], N[Not[LessEqual[phi1, 5.4e-141]], $MachinePrecision]], N[(N[ArcCos[N[Cos[N[(phi2 - phi1), $MachinePrecision]], $MachinePrecision]], $MachinePrecision] * R), $MachinePrecision], N[(N[ArcCos[N[(N[(-0.5 * N[(phi1 * phi1), $MachinePrecision] + 1.0), $MachinePrecision] * N[Cos[lambda1], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * R), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\phi_1 \leq -6.8 \cdot 10^{-122} \lor \neg \left(\phi_1 \leq 5.4 \cdot 10^{-141}\right):\\
\;\;\;\;\cos^{-1} \cos \left(\phi_2 - \phi_1\right) \cdot R\\
\mathbf{else}:\\
\;\;\;\;\cos^{-1} \left(\mathsf{fma}\left(-0.5, \phi_1 \cdot \phi_1, 1\right) \cdot \cos \lambda_1\right) \cdot R\\
\end{array}
\end{array}
if phi1 < -6.7999999999999996e-122 or 5.4000000000000005e-141 < phi1 Initial program 77.9%
Taylor expanded in lambda2 around 0
*-commutativeN/A
fp-cancel-sign-sub-invN/A
fp-cancel-sub-sign-invN/A
*-commutativeN/A
sin-neg-revN/A
sin-neg-revN/A
remove-double-negN/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-sin.f64N/A
lower-sin.f6456.3
Applied rewrites56.3%
Taylor expanded in lambda1 around 0
Applied rewrites34.2%
Applied rewrites24.4%
if -6.7999999999999996e-122 < phi1 < 5.4000000000000005e-141Initial program 63.7%
Taylor expanded in phi2 around 0
*-commutativeN/A
lower-*.f64N/A
*-lft-identityN/A
metadata-evalN/A
fp-cancel-sign-sub-invN/A
cos-neg-revN/A
remove-double-negN/A
mul-1-negN/A
distribute-neg-inN/A
+-commutativeN/A
mul-1-negN/A
lower-cos.f64N/A
remove-double-negN/A
fp-cancel-sign-sub-invN/A
metadata-evalN/A
*-lft-identityN/A
lower--.f64N/A
lower-cos.f6444.0
Applied rewrites44.0%
Taylor expanded in phi1 around 0
Applied rewrites44.0%
Taylor expanded in lambda2 around 0
Applied rewrites22.5%
Final simplification23.9%
(FPCore (R lambda1 lambda2 phi1 phi2) :precision binary64 (if (<= phi1 -1.06) (* (acos (cos phi1)) R) (* (acos (* (fma -0.5 (* phi1 phi1) 1.0) (cos (- lambda1 lambda2)))) R)))
double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
double tmp;
if (phi1 <= -1.06) {
tmp = acos(cos(phi1)) * R;
} else {
tmp = acos((fma(-0.5, (phi1 * phi1), 1.0) * cos((lambda1 - lambda2)))) * R;
}
return tmp;
}
function code(R, lambda1, lambda2, phi1, phi2) tmp = 0.0 if (phi1 <= -1.06) tmp = Float64(acos(cos(phi1)) * R); else tmp = Float64(acos(Float64(fma(-0.5, Float64(phi1 * phi1), 1.0) * cos(Float64(lambda1 - lambda2)))) * R); end return tmp end
code[R_, lambda1_, lambda2_, phi1_, phi2_] := If[LessEqual[phi1, -1.06], N[(N[ArcCos[N[Cos[phi1], $MachinePrecision]], $MachinePrecision] * R), $MachinePrecision], N[(N[ArcCos[N[(N[(-0.5 * N[(phi1 * phi1), $MachinePrecision] + 1.0), $MachinePrecision] * N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * R), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\phi_1 \leq -1.06:\\
\;\;\;\;\cos^{-1} \cos \phi_1 \cdot R\\
\mathbf{else}:\\
\;\;\;\;\cos^{-1} \left(\mathsf{fma}\left(-0.5, \phi_1 \cdot \phi_1, 1\right) \cdot \cos \left(\lambda_1 - \lambda_2\right)\right) \cdot R\\
\end{array}
\end{array}
if phi1 < -1.0600000000000001Initial program 84.5%
Taylor expanded in lambda2 around 0
*-commutativeN/A
fp-cancel-sign-sub-invN/A
fp-cancel-sub-sign-invN/A
*-commutativeN/A
sin-neg-revN/A
sin-neg-revN/A
remove-double-negN/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-sin.f64N/A
lower-sin.f6466.9
Applied rewrites66.9%
Taylor expanded in lambda1 around 0
Applied rewrites45.8%
Applied rewrites27.7%
Taylor expanded in phi2 around 0
Applied rewrites27.8%
if -1.0600000000000001 < phi1 Initial program 70.4%
Taylor expanded in phi2 around 0
*-commutativeN/A
lower-*.f64N/A
*-lft-identityN/A
metadata-evalN/A
fp-cancel-sign-sub-invN/A
cos-neg-revN/A
remove-double-negN/A
mul-1-negN/A
distribute-neg-inN/A
+-commutativeN/A
mul-1-negN/A
lower-cos.f64N/A
remove-double-negN/A
fp-cancel-sign-sub-invN/A
metadata-evalN/A
*-lft-identityN/A
lower--.f64N/A
lower-cos.f6443.4
Applied rewrites43.4%
Taylor expanded in phi1 around 0
Applied rewrites27.7%
(FPCore (R lambda1 lambda2 phi1 phi2) :precision binary64 (if (<= phi1 -1.3e-6) (* (acos (cos phi1)) R) (* (acos (cos phi2)) R)))
double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
double tmp;
if (phi1 <= -1.3e-6) {
tmp = acos(cos(phi1)) * R;
} else {
tmp = acos(cos(phi2)) * R;
}
return tmp;
}
real(8) function code(r, lambda1, lambda2, phi1, phi2)
real(8), intent (in) :: r
real(8), intent (in) :: lambda1
real(8), intent (in) :: lambda2
real(8), intent (in) :: phi1
real(8), intent (in) :: phi2
real(8) :: tmp
if (phi1 <= (-1.3d-6)) then
tmp = acos(cos(phi1)) * r
else
tmp = acos(cos(phi2)) * r
end if
code = tmp
end function
public static double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
double tmp;
if (phi1 <= -1.3e-6) {
tmp = Math.acos(Math.cos(phi1)) * R;
} else {
tmp = Math.acos(Math.cos(phi2)) * R;
}
return tmp;
}
def code(R, lambda1, lambda2, phi1, phi2): tmp = 0 if phi1 <= -1.3e-6: tmp = math.acos(math.cos(phi1)) * R else: tmp = math.acos(math.cos(phi2)) * R return tmp
function code(R, lambda1, lambda2, phi1, phi2) tmp = 0.0 if (phi1 <= -1.3e-6) tmp = Float64(acos(cos(phi1)) * R); else tmp = Float64(acos(cos(phi2)) * R); end return tmp end
function tmp_2 = code(R, lambda1, lambda2, phi1, phi2) tmp = 0.0; if (phi1 <= -1.3e-6) tmp = acos(cos(phi1)) * R; else tmp = acos(cos(phi2)) * R; end tmp_2 = tmp; end
code[R_, lambda1_, lambda2_, phi1_, phi2_] := If[LessEqual[phi1, -1.3e-6], N[(N[ArcCos[N[Cos[phi1], $MachinePrecision]], $MachinePrecision] * R), $MachinePrecision], N[(N[ArcCos[N[Cos[phi2], $MachinePrecision]], $MachinePrecision] * R), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\phi_1 \leq -1.3 \cdot 10^{-6}:\\
\;\;\;\;\cos^{-1} \cos \phi_1 \cdot R\\
\mathbf{else}:\\
\;\;\;\;\cos^{-1} \cos \phi_2 \cdot R\\
\end{array}
\end{array}
if phi1 < -1.30000000000000005e-6Initial program 82.5%
Taylor expanded in lambda2 around 0
*-commutativeN/A
fp-cancel-sign-sub-invN/A
fp-cancel-sub-sign-invN/A
*-commutativeN/A
sin-neg-revN/A
sin-neg-revN/A
remove-double-negN/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-sin.f64N/A
lower-sin.f6465.4
Applied rewrites65.4%
Taylor expanded in lambda1 around 0
Applied rewrites44.9%
Applied rewrites27.4%
Taylor expanded in phi2 around 0
Applied rewrites27.4%
if -1.30000000000000005e-6 < phi1 Initial program 71.0%
Taylor expanded in lambda2 around 0
*-commutativeN/A
fp-cancel-sign-sub-invN/A
fp-cancel-sub-sign-invN/A
*-commutativeN/A
sin-neg-revN/A
sin-neg-revN/A
remove-double-negN/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-sin.f64N/A
lower-sin.f6446.7
Applied rewrites46.7%
Taylor expanded in lambda1 around 0
Applied rewrites25.1%
Applied rewrites21.7%
Taylor expanded in phi1 around 0
Applied rewrites17.5%
(FPCore (R lambda1 lambda2 phi1 phi2) :precision binary64 (* (acos (cos (- phi2 phi1))) R))
double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
return acos(cos((phi2 - phi1))) * R;
}
real(8) function code(r, lambda1, lambda2, phi1, phi2)
real(8), intent (in) :: r
real(8), intent (in) :: lambda1
real(8), intent (in) :: lambda2
real(8), intent (in) :: phi1
real(8), intent (in) :: phi2
code = acos(cos((phi2 - phi1))) * r
end function
public static double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
return Math.acos(Math.cos((phi2 - phi1))) * R;
}
def code(R, lambda1, lambda2, phi1, phi2): return math.acos(math.cos((phi2 - phi1))) * R
function code(R, lambda1, lambda2, phi1, phi2) return Float64(acos(cos(Float64(phi2 - phi1))) * R) end
function tmp = code(R, lambda1, lambda2, phi1, phi2) tmp = acos(cos((phi2 - phi1))) * R; end
code[R_, lambda1_, lambda2_, phi1_, phi2_] := N[(N[ArcCos[N[Cos[N[(phi2 - phi1), $MachinePrecision]], $MachinePrecision]], $MachinePrecision] * R), $MachinePrecision]
\begin{array}{l}
\\
\cos^{-1} \cos \left(\phi_2 - \phi_1\right) \cdot R
\end{array}
Initial program 73.9%
Taylor expanded in lambda2 around 0
*-commutativeN/A
fp-cancel-sign-sub-invN/A
fp-cancel-sub-sign-invN/A
*-commutativeN/A
sin-neg-revN/A
sin-neg-revN/A
remove-double-negN/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-sin.f64N/A
lower-sin.f6451.5
Applied rewrites51.5%
Taylor expanded in lambda1 around 0
Applied rewrites30.2%
Applied rewrites23.2%
(FPCore (R lambda1 lambda2 phi1 phi2) :precision binary64 (* (acos (cos phi1)) R))
double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
return acos(cos(phi1)) * R;
}
real(8) function code(r, lambda1, lambda2, phi1, phi2)
real(8), intent (in) :: r
real(8), intent (in) :: lambda1
real(8), intent (in) :: lambda2
real(8), intent (in) :: phi1
real(8), intent (in) :: phi2
code = acos(cos(phi1)) * r
end function
public static double code(double R, double lambda1, double lambda2, double phi1, double phi2) {
return Math.acos(Math.cos(phi1)) * R;
}
def code(R, lambda1, lambda2, phi1, phi2): return math.acos(math.cos(phi1)) * R
function code(R, lambda1, lambda2, phi1, phi2) return Float64(acos(cos(phi1)) * R) end
function tmp = code(R, lambda1, lambda2, phi1, phi2) tmp = acos(cos(phi1)) * R; end
code[R_, lambda1_, lambda2_, phi1_, phi2_] := N[(N[ArcCos[N[Cos[phi1], $MachinePrecision]], $MachinePrecision] * R), $MachinePrecision]
\begin{array}{l}
\\
\cos^{-1} \cos \phi_1 \cdot R
\end{array}
Initial program 73.9%
Taylor expanded in lambda2 around 0
*-commutativeN/A
fp-cancel-sign-sub-invN/A
fp-cancel-sub-sign-invN/A
*-commutativeN/A
sin-neg-revN/A
sin-neg-revN/A
remove-double-negN/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
lower-cos.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-sin.f64N/A
lower-sin.f6451.5
Applied rewrites51.5%
Taylor expanded in lambda1 around 0
Applied rewrites30.2%
Applied rewrites23.2%
Taylor expanded in phi2 around 0
Applied rewrites15.6%
herbie shell --seed 2024343
(FPCore (R lambda1 lambda2 phi1 phi2)
:name "Spherical law of cosines"
:precision binary64
(* (acos (+ (* (sin phi1) (sin phi2)) (* (* (cos phi1) (cos phi2)) (cos (- lambda1 lambda2))))) R))