Admissible Region Inner Boundary (#69)

* Update constants.jl

* AbstractFloat --> Real in admissibleregion.jl

* Add arG(::AdmissibleRegion{T}, ::S)

* Update arenergycoeffs and arenergydis

* Update rangerate

* Update argoldensearch

* Update maxrange

* Update boundary

* Update NEOsRecipesBaseExt.jl

* Solve rangerate ambiguity

* Add tests

* Bump patch version

Author: LuEdRaMo

Project.toml

@@ -1,7 +1,7 @@
 name = "NEOs"
 uuid = "b41c07a2-2abb-11e9-070a-c3c1b239e7df"
 authors = ["Jorge A. Pérez Hernández", "Luis Benet", "Luis Eduardo Ramírez Montoya"]
-version = "0.8.6"
+version = "0.8.7"
 
 [deps]
 Artifacts = "56f22d72-fd6d-98f1-02f0-08ddc0907c33"

ext/NEOsRecipesBaseExt.jl

@@ -2,16 +2,18 @@ module NEOsRecipesBaseExt
 
 using RecipesBase
 using PlanetaryEphemeris: TaylorInterpolant
-using NEOs: OpticalResidual, AdmissibleRegion, NEOSolution, cte, ra, dec, boundary
+using NEOs: OpticalResidual, AdmissibleRegion, NEOSolution, cte, ra, dec, arboundary
 
 @recipe function f(res::AbstractVector{OpticalResidual{T, U}}) where {T <: Real, U <: Number}
     seriestype --> :scatter
     return cte.(ra.(res)), cte.(dec.(res))
 end
 
-@recipe function f(A::AdmissibleRegion{T}; N::Int = 100, ρscale::Symbol = :linear) where {T <: Real}
+@recipe function f(A::AdmissibleRegion{T}; boundary::Symbol = :outer,
+                   N::Int = 100, ρscale::Symbol = :linear) where {T <: Real}
     seriestype --> :path
-    ps = map(t -> boundary(A, t), LinRange(0, 3, N))
+    tmax = boundary == :outer ? 3 : 2
+    ps = map(t -> arboundary(A, t, boundary), LinRange(0, tmax, N))
     xs, ys = first.(ps), last.(ps)
     if ρscale == :log
         xs .= log10.(xs)

src/constants.jl

@@ -203,4 +203,6 @@ const k_gauss = 0.017_202_098_95
 # Earth's sphere of influence radius [AU]
 const R_SI = 0.010044
 # Earth's physical radius [AU]
-const R_EA = 4.2e-5
+const R_EA = 4.24e-5
+# Ratio between the mass of the Earth and the mass of the Sun
+const μ_ES = PE.μ[ea] / PE.μ[su] # 1 / 328_900.5614