Profile code

[gomezzz]

Please post some cprofile outputs to see where our runtime is going, it is starting to be quite slow :pray:

[rasmusmarak]

Test 1:

(From here on after referred to as the "Basic problem").

Search space: Cartesian coordinates [r_x,r_y,r_z,v_x,v_y,v_z]

• lower_bound = [-6000,-6000,-6000,0,0,0]
• upper_bound = [6000,6000,6000,2.5,2.5,2.5]

Integration time:

• t0 = 0
• tf = 20*3600 (i.e 20h)

Problem parameters:

• body density = 533
• body mu = 665.666
• target altitude = 5000 [m]

Algorithm parameters:

• initial step = 600 [seconds]
• integration algorithm = DP8713M (using Desolver)
• population size = 10
• number of evolutions = 32
• optimization algorithm = SADE (self-adaptive differential evolution, using pygmo)

Results:

• Algorithm run time = 38min 9.2s
• champion integration time = 4.5s
• champion fitness value = 1599897.43729276
• champion chromosome = [ 3.54273804e+03, -2.05937976e+03, -2.90055903e+03, 3.45286575e-01, 4.45044180e-01, 1.43148382e-01]
• F = 0.259343
• CR = 0.626794
• dx = 8507.09
• df = 1.06143e+08

Trajectory with mesh:

Trajectory (x and y axis):

C-profile: Sorted by number of calls: output_calls.txt Sorted by time: output_time.txt

[rasmusmarak]

Test 2:

Basic Problem structure and the following changes:

• Changed to osculating orbital elements.

Setup: Search space: Osculating orbital elements [a, e, o, w, i, ea]

• lower_bound = [5000, 0, 0, 0, 0, 0]
• upper_bound = [5000, 1, 2pi, 2pi, pi, 2*pi]

Integration time:

• t0 = 0
• tf = 20*3600 (i.e 20h)

Problem parameters:

• body density = 533
• body mu = 665.666
• target altitude = 5000 [m]

Algorithm parameters:

• initial step = 600 [seconds]
• integration algorithm = DP8713M (using Desolver)
• population size = 10
• number of evolutions = 32
• optimization algorithm = SADE (self-adaptive differential evolution, using pygmo)

Results:

• Algorithm run time = 57min 40.5s
• champion integration time = 7.2s
• champion fitness value = 11824891.30506848
• champion chromosome = [1.11372603e+04, 4.68773238e-01, 2.20821729e+00, 4.51629780e-01, 3.33853297e+00, 3.48571918e-01]
• F = 0.774538
• CR = 0.26791
• dx = 2724.48
• df = 3.53187e+06

Trajectory with mesh:

Trajectory (x and y axis):

C-profile: Sorted by number of calls: output_calls.txt Sorted by time: output_time.txt

[rasmusmarak]

Test 3:

Basic problem structure and the following changes:

• changed to osculating orbital elements
• added collision detection

Setup: Search space: Osculating orbital elements [a, e, o, w, i, ea]

• lower_bound = [5000, 0, 0, 0, 0, 0]
• upper_bound = [5000, 1, 2pi, 2pi, pi, 2*pi]

Integration time:

• t0 = 0
• tf = 20*3600 (i.e 20h)

Problem parameters:

• body density = 533
• body mu = 665.666
• target altitude = 5000 [m]
• radius bounding sphere = 3500 [m]
• collision penalty = 1e30

Algorithm parameters:

• initial step = 600 [seconds]
• integration algorithm = DP8713M (using Desolver)
• population size = 10
• number of evolutions = 32
• optimization algorithm = SADE (self-adaptive differential evolution, using pygmo)

Results:

• Algorithm run time = 35min 53.4s
• champion integration time = 4.4s
• champion fitness value = 10030080.69184285
• champion chromosome = [1.02346115e+04, 5.04262474e-01, 1.40521347e+00, 3.03148072e+00, 2.68878957e-01, 5.74690265e+00]
• F = 0.405411
• CR = 0.752927
• dx = 3445.18
• df = 2.10126e+07

Trajectory with mesh:

Trajectory (x and y axis):

C-profile: Sorted by number of calls: output_calls.txt Sorted by time: output_time.txt

[rasmusmarak]

Test 4:

Basic problem structure and the following changes:

• changed to osculating orbital elements
• added collision detection
• added rotation of satellite position using quaternions

Setup: Search space: Osculating orbital elements [a, e, o, w, i, ea]

• lower_bound = [5000, 0, 0, 0, 0, 0]
• upper_bound = [5000, 1, 2pi, 2pi, pi, 2*pi]

Integration time:

• t0 = 0
• tf = 20*3600 (i.e 20h)

Problem parameters:

• body density = 533
• body mu = 665.666
• target altitude = 5000 [m]
• radius bounding sphere = 3500 [m]
• collision penalty = 1e30
• spin axis declination = 64 [degrees]
• spin axis right ascension = 69 [degrees]
• spin period = 12.06*3600 [seconds] (i.e 12.06h)

Algorithm parameters:

• initial step = 600 [seconds]
• integration algorithm = DP8713M (using Desolver)
• population size = 10
• number of evolutions = 32
• optimization algorithm = SADE (self-adaptive differential evolution, using pygmo)

Results:

• Algorithm run time = 15min 35.9s
• champion integration time = 3.5s
• champion fitness value = 20938771.53471591
• champion chromosome = [9.72752753e+03, 3.64128116e-02, 1.28319173e+00, 1.69622584e+00, 2.18503270e+00, 2.66184179e+00]
• F = 0.763792
• CR = 0.686998
• dx = 1606.6
• df = 6.64294e+07

Trajectory with mesh:

Trajectory (x and y axis):

C-profile: Sorted by number of calls: output_calls.txt Sorted by time: output_time.txt

[gomezzz]

In future plots ideally add fitness over generations

[rasmusmarak]

Profile code:

Number of spacecrafts = 1 Number of maneuvers = 1 Mission time = 2419200 [s] (ie 4 weeks, or 28 days)

Population size = 120 Number of generations = 10 Number of threads = 40

Penalty scaling factor = 0.1 Max velocity scaling factor = 40

Radius bounding spheres = [4000, 15000]. (i.e inner and outer) Rotation of body = True

output_calls.txt output_time.txt