diff --git a/julia/src/inner_loop/inner_loop.jl b/julia/src/inner_loop/inner_loop.jl index eca8409..c0ba6f5 100644 --- a/julia/src/inner_loop/inner_loop.jl +++ b/julia/src/inner_loop/inner_loop.jl @@ -44,7 +44,7 @@ function inner_loop(launch_date::DateTime, end time = utc2et(Dates.format(launch_date,"yyyy-mm-ddTHH:MM:SS")) - thrust_profiles = [] + thrust_profiles = Vector{Matrix{Float64}}() for phase in phases planet1_state = [spkssb(ids[phase.from_planet], time, "ECLIPJ2000"); 0.0] diff --git a/julia/test/inner_loop/inner_loop.jl b/julia/test/inner_loop/inner_loop.jl index c612697..26c503d 100644 --- a/julia/test/inner_loop/inner_loop.jl +++ b/julia/test/inner_loop/inner_loop.jl @@ -47,15 +47,14 @@ # The first leg should be valid fresh_sc = copy(sc) - mass, thrusts = inner_loop( launch_date, - sc, - start_mass, - [phase1], - verbose=true, - mbh_specs=(25,50) ) - @test sc.mass > sc.dry_mass - path, final = prop(thrusts[1], earth_state+[zeros(3);v∞s[1]], fresh_sc, μs["Sun"], leg1_tof) - @test final ≈ p1_state + [zeros(3); v∞s[2]] + thrusts = inner_loop( launch_date, + sc, + start_mass, + [phase1], + verbose=true, + mbh_specs=(25,50) ) + path, final = prop(thrusts[1], earth_state+[zeros(3);v∞s[1];0.0], Sc("test"), μs["Sun"], leg1_tof) + @test final[1:6] ≈ p1_state[1:6] + [zeros(3); v∞s[2]] savefig(plot_orbits( [earth, p1_path, path], primary="Sun", labels=["Earth", p1, "Leg 1"],