Matlab code for 2016 AAS paper
Switch branches/tags
Nothing to show
Clone or download
Fetching latest commit…
Cannot retrieve the latest commit at this time.
Permalink
Type Name Latest commit message Commit time
Failed to load latest commit information.
CASTALIA
.gitignore
LICENSE.txt
README.md
ROT1.m
ROT2.m
ROT3.m
angle_contour.m
ast_eoms.m
ast_optimal_eoms.m
ast_stm.m
asteroid_optimal_ode.m
asteroid_shooting.m
asteroid_shooting_min.m
asteroid_shooting_z.m
asteroid_transfer_analysis.m
asteroid_transfer_analysis_vert.m
asteroid_transfer_analysis_z.m
asteroid_transfer_driver.m
asteroid_transfer_driver_vert.m
asteroid_transfer_driver_z.m
bw_ast_eoms.m
colonialone_asteroid.sh
colonialone_asteroid_4convert.sh
colonialone_asteroid_stage1.sh
colonialone_asteroid_z.sh
constraints_nodenom.m
constraints_nodenom_z.m
contour_plane_ex.m
dim2nondim.m
equilibria_behavior.m
equilibria_finder.m
escape_speed.m
events_xcross.m
gradient_test.m
grav_contour.m
grav_potential.m
hpc_test.m
hpc_test_4convert.m
hpc_test_stage1.m
hpc_test_z.m
initial_condition_finder.m
initialize_transfer.m
initialize_transfer_orbit.m
jacobi.m
jacobi_constraint.m
jacobi_continuation.m
linearized_mat.m
load_constants.m
mascon_potential.m
minimum_reach.m
objective.m
objective_z.m
periodic_orbit_diffcorr.m
plot_animation.m
plot_asteroid_animation.m
plot_output.m
poincare_angles.m
polyhedron_potential.m
polyhedron_potential_test.m
polyhedron_shape_input.m
polyhedron_timing.m
pot_contour.m
potential_model_comparison.m
prop_seg.m
radius_contour.m
return_speed.m
rot_asteroid.m
seq_newton.m
seq_newton_z.m
setup_repo.sh
sim_driver.m
surface_slope.m
vertex_plotter.m
zvc_contour.m

README.md

Running Asteroid Shooting README

##ColonialOne login ssh sklumani@login.colonialone.gwu.edu

skulumani and GWU lastpass password

Scheduler commands

sbatch name_of_script.sh squeue -l sinfo

Copy a file from colonialone to the local compute

pscp skulumani@login.colonialone.gwu.edu:/home/skulumani/asteroid/fine_asteroid_1.mat C:\Users\skulumani\Desktop\fine_asteroid_1.mat

Send a file from SEAS computer to colonial one

pscp C:\Users\skulumani\Desktop\fine_asteroid_1.mat skulumani@login.colonialone.gwu.edu:/home/skulumani/asteroid/fine_asteroid_1.mat

Running the Code on local computer

run hpc_test.m

The number of angles is defined in asteroid_shooting.m - this defines the density of the reachability set (discretization)

Run on HPC

Make sure you have the compiled version of polyhedron_potential_mex_1024 Modify file names in colonialone_asteroid.sh Modify file name in asteroid_shooting.m for final mat file

MEX Function

Need to compile polyhedron_potential.m using MEX Use codegen to generate MEX function for polyhedron_potential.m Test function is given by polyhedron_potential_test.m Modify the number of faces using 'false' - full number of faces or 'true' - 1024 faces

Transfer examples

hpc_test.m/colonialone_asteroid.sh - 4 constraint transfer hpc_test_4convert.m/colonialone_asteroid_4convert.sh - Extra vertical weighting hpc_test_z.m/colonialone_asteroid_z.sh - 5 constraint with the added z=0 constraint

asteroid_transfer_analysis.m - Plots the transfers using the saved mat files in ./results

Final transfer

The results in ./results/4constraint_transfer shows the results used in the paper for 2016 AAS

You can generate the the plots by running asteroid_transfer_analysis.m

This is also stored in the git repository by checking out the 2016AAS tag