In [1]:
import sys
sys.path.append('../')
from numpy import *
from PyGMO import *
import matplotlib.pyplot as plt
from Optimisation import HSS
from Trajectory import Point_Lander_Drag
In [2]:
model = Point_Lander_Drag()
prob = HSS(model, nsegs=10)
In [3]:
zguess = prob.Guess.Ballistic(tf=20)
pop = population(prob)
pop.push_back(zguess)
In [4]:
algo1 = algorithm.scipy_slsqp(max_iter=5000, screen_output=True)
algo2 = algorithm.mbh(algo1, stop=1)
In [5]:
pop = algo2.evolve(pop)
NIT FC OBJFUN GNORM
1 150 -8.482912E+03 1.000000E+00
2 300 -8.482993E+03 1.000000E+00
3 450 -8.483052E+03 1.000000E+00
4 600 -8.483089E+03 1.000000E+00
5 750 -8.482861E+03 1.000000E+00
6 900 -8.482986E+03 1.000000E+00
7 1050 -8.482991E+03 1.000000E+00
8 1200 -8.482915E+03 1.000000E+00
9 1350 -8.482824E+03 1.000000E+00
10 1500 -8.482830E+03 1.000000E+00
11 1650 -8.482828E+03 1.000000E+00
12 1800 -8.482772E+03 1.000000E+00
13 1950 -8.482608E+03 1.000000E+00
14 2100 -8.482396E+03 1.000000E+00
15 2250 -8.481533E+03 1.000000E+00
16 2400 -8.480103E+03 1.000000E+00
17 2550 -8.477296E+03 1.000000E+00
18 2700 -8.467672E+03 1.000000E+00
19 2850 -8.464246E+03 1.000000E+00
20 3000 -8.461687E+03 1.000000E+00
21 3150 -8.455930E+03 1.000000E+00
22 3300 -8.451181E+03 1.000000E+00
23 3450 -8.443552E+03 1.000000E+00
24 3600 -8.437633E+03 1.000000E+00
25 3750 -8.417908E+03 1.000000E+00
26 3900 -8.410769E+03 1.000000E+00
27 4050 -8.397133E+03 1.000000E+00
28 4200 -8.382473E+03 1.000000E+00
29 4350 -8.378681E+03 1.000000E+00
30 4500 -8.370041E+03 1.000000E+00
31 4650 -8.364174E+03 1.000000E+00
32 4800 -8.356865E+03 1.000000E+00
33 4950 -8.338859E+03 1.000000E+00
34 5100 -8.313003E+03 1.000000E+00
35 5250 -8.281588E+03 1.000000E+00
36 5400 -8.274476E+03 1.000000E+00
37 5550 -8.269272E+03 1.000000E+00
38 5700 -8.256115E+03 1.000000E+00
39 5850 -8.192987E+03 1.000000E+00
40 6000 -8.101860E+03 1.000000E+00
41 6150 -8.031202E+03 1.000000E+00
42 6300 -7.991108E+03 1.000000E+00
43 6450 -7.974354E+03 1.000000E+00
44 6600 -7.966255E+03 1.000000E+00
45 6750 -7.958181E+03 1.000000E+00
46 6900 -7.948371E+03 1.000000E+00
47 7050 -7.933764E+03 1.000000E+00
48 7200 -7.918846E+03 1.000000E+00
49 7350 -7.875535E+03 1.000000E+00
50 7500 -7.818820E+03 1.000000E+00
51 7650 -7.781659E+03 1.000000E+00
52 7800 -7.734002E+03 1.000000E+00
53 7950 -7.655873E+03 1.000000E+00
54 8100 -6.396460E+03 1.000000E+00
55 8250 -6.834459E+03 1.000000E+00
56 8401 -6.809197E+03 1.000000E+00
57 8552 -6.721999E+03 1.000000E+00
58 8702 -6.477576E+03 1.000000E+00
59 8852 -6.445635E+03 1.000000E+00
60 9002 -6.439620E+03 1.000000E+00
61 9152 -6.439032E+03 1.000000E+00
62 9302 -6.444588E+03 1.000000E+00
63 9453 -6.456430E+03 1.000000E+00
64 9604 -6.476888E+03 1.000000E+00
65 9755 -6.505774E+03 1.000000E+00
66 9905 -6.548567E+03 1.000000E+00
67 10055 -6.549429E+03 1.000000E+00
68 10205 -6.577393E+03 1.000000E+00
69 10355 -6.579280E+03 1.000000E+00
70 10505 -6.599184E+03 1.000000E+00
71 10655 -6.607468E+03 1.000000E+00
72 10805 -6.612891E+03 1.000000E+00
73 10955 -6.619385E+03 1.000000E+00
74 11105 -6.635718E+03 1.000000E+00
75 11255 -6.640026E+03 1.000000E+00
76 11405 -6.642403E+03 1.000000E+00
77 11555 -6.645822E+03 1.000000E+00
78 11705 -6.660448E+03 1.000000E+00
79 11855 -6.660512E+03 1.000000E+00
80 12005 -6.674186E+03 1.000000E+00
81 12155 -6.680611E+03 1.000000E+00
82 12305 -6.689691E+03 1.000000E+00
83 12455 -6.692324E+03 1.000000E+00
84 12605 -6.701481E+03 1.000000E+00
85 12755 -6.703050E+03 1.000000E+00
86 12905 -6.706335E+03 1.000000E+00
87 13055 -6.709228E+03 1.000000E+00
88 13205 -6.711342E+03 1.000000E+00
89 13355 -6.721172E+03 1.000000E+00
90 13505 -6.722326E+03 1.000000E+00
91 13655 -6.729333E+03 1.000000E+00
92 13805 -6.730348E+03 1.000000E+00
93 13955 -6.734453E+03 1.000000E+00
94 14105 -6.741146E+03 1.000000E+00
95 14255 -6.741466E+03 1.000000E+00
96 14405 -6.743239E+03 1.000000E+00
97 14555 -6.743657E+03 1.000000E+00
98 14705 -6.744543E+03 1.000000E+00
99 14855 -6.744932E+03 1.000000E+00
100 15005 -6.745492E+03 1.000000E+00
101 15155 -6.745858E+03 1.000000E+00
102 15305 -6.747075E+03 1.000000E+00
103 15455 -6.746708E+03 1.000000E+00
104 15605 -6.747046E+03 1.000000E+00
105 15755 -6.747287E+03 1.000000E+00
106 15905 -6.747399E+03 1.000000E+00
107 16055 -6.747762E+03 1.000000E+00
108 16205 -6.748169E+03 1.000000E+00
109 16355 -6.748981E+03 1.000000E+00
110 16505 -6.749689E+03 1.000000E+00
111 16655 -6.750034E+03 1.000000E+00
112 16805 -6.750557E+03 1.000000E+00
113 16955 -6.750691E+03 1.000000E+00
114 17105 -6.751473E+03 1.000000E+00
115 17255 -6.752462E+03 1.000000E+00
116 17405 -6.754044E+03 1.000000E+00
117 17555 -6.755397E+03 1.000000E+00
118 17705 -6.756003E+03 1.000000E+00
119 17855 -6.756616E+03 1.000000E+00
120 18005 -6.757441E+03 1.000000E+00
121 18155 -6.758760E+03 1.000000E+00
122 18305 -6.760164E+03 1.000000E+00
123 18455 -6.760964E+03 1.000000E+00
124 18605 -6.762284E+03 1.000000E+00
125 18755 -6.762487E+03 1.000000E+00
126 18905 -6.763263E+03 1.000000E+00
127 19055 -6.763999E+03 1.000000E+00
128 19205 -6.765010E+03 1.000000E+00
129 19355 -6.766104E+03 1.000000E+00
130 19505 -6.767535E+03 1.000000E+00
131 19655 -6.769373E+03 1.000000E+00
132 19805 -6.772250E+03 1.000000E+00
133 19955 -6.775381E+03 1.000000E+00
134 20105 -6.776771E+03 1.000000E+00
135 20255 -6.778793E+03 1.000000E+00
136 20405 -6.781412E+03 1.000000E+00
137 20555 -6.783912E+03 1.000000E+00
138 20705 -6.785728E+03 1.000000E+00
139 20855 -6.788095E+03 1.000000E+00
140 21005 -6.787426E+03 1.000000E+00
141 21155 -6.788008E+03 1.000000E+00
142 21305 -6.789036E+03 1.000000E+00
143 21455 -6.789050E+03 1.000000E+00
144 21605 -6.790136E+03 1.000000E+00
145 21755 -6.790483E+03 1.000000E+00
146 21905 -6.790989E+03 1.000000E+00
147 22055 -6.791351E+03 1.000000E+00
148 22205 -6.792005E+03 1.000000E+00
149 22355 -6.793261E+03 1.000000E+00
150 22505 -6.794161E+03 1.000000E+00
151 22655 -6.795120E+03 1.000000E+00
152 22805 -6.795542E+03 1.000000E+00
153 22955 -6.796614E+03 1.000000E+00
154 23105 -6.798243E+03 1.000000E+00
155 23255 -6.799823E+03 1.000000E+00
156 23405 -6.800408E+03 1.000000E+00
157 23555 -6.800745E+03 1.000000E+00
158 23705 -6.800975E+03 1.000000E+00
159 23855 -6.801154E+03 1.000000E+00
160 24005 -6.801496E+03 1.000000E+00
161 24155 -6.801686E+03 1.000000E+00
162 24305 -6.802029E+03 1.000000E+00
163 24455 -6.802466E+03 1.000000E+00
164 24605 -6.803207E+03 1.000000E+00
165 24755 -6.804075E+03 1.000000E+00
166 24905 -6.804464E+03 1.000000E+00
167 25055 -6.804464E+03 1.000000E+00
168 25205 -6.804614E+03 1.000000E+00
169 25355 -6.804714E+03 1.000000E+00
170 25505 -6.805183E+03 1.000000E+00
171 25655 -6.805971E+03 1.000000E+00
172 25805 -6.806588E+03 1.000000E+00
173 25955 -6.807893E+03 1.000000E+00
174 26105 -6.807509E+03 1.000000E+00
175 26255 -6.807697E+03 1.000000E+00
176 26405 -6.807898E+03 1.000000E+00
177 26555 -6.807949E+03 1.000000E+00
178 26705 -6.808377E+03 1.000000E+00
179 26855 -6.808763E+03 1.000000E+00
180 27005 -6.809013E+03 1.000000E+00
181 27155 -6.809364E+03 1.000000E+00
182 27305 -6.809239E+03 1.000000E+00
183 27455 -6.809294E+03 1.000000E+00
184 27605 -6.809388E+03 1.000000E+00
185 27755 -6.809559E+03 1.000000E+00
186 27905 -6.809717E+03 1.000000E+00
187 28055 -6.809730E+03 1.000000E+00
188 28205 -6.809792E+03 1.000000E+00
189 28355 -6.809785E+03 1.000000E+00
190 28505 -6.809833E+03 1.000000E+00
191 28655 -6.809890E+03 1.000000E+00
192 28805 -6.809998E+03 1.000000E+00
193 28955 -6.810085E+03 1.000000E+00
194 29105 -6.810165E+03 1.000000E+00
195 29255 -6.810382E+03 1.000000E+00
196 29405 -6.810314E+03 1.000000E+00
197 29555 -6.810495E+03 1.000000E+00
198 29705 -6.810769E+03 1.000000E+00
199 29855 -6.811399E+03 1.000000E+00
200 30005 -6.811989E+03 1.000000E+00
201 30155 -6.813129E+03 1.000000E+00
202 30305 -6.812864E+03 1.000000E+00
203 30455 -6.812836E+03 1.000000E+00
204 30605 -6.812886E+03 1.000000E+00
205 30755 -6.812910E+03 1.000000E+00
206 30905 -6.812962E+03 1.000000E+00
207 31055 -6.813021E+03 1.000000E+00
208 31205 -6.813118E+03 1.000000E+00
209 31355 -6.813166E+03 1.000000E+00
210 31505 -6.813223E+03 1.000000E+00
211 31655 -6.813285E+03 1.000000E+00
212 31805 -6.813282E+03 1.000000E+00
213 31955 -6.813300E+03 1.000000E+00
214 32105 -6.813311E+03 1.000000E+00
215 32255 -6.813375E+03 1.000000E+00
216 32405 -6.813493E+03 1.000000E+00
217 32555 -6.813790E+03 1.000000E+00
218 32705 -6.814325E+03 1.000000E+00
219 32855 -6.814943E+03 1.000000E+00
220 33005 -6.815261E+03 1.000000E+00
221 33155 -6.815691E+03 1.000000E+00
222 33305 -6.815608E+03 1.000000E+00
223 33456 -6.815686E+03 1.000000E+00
224 33606 -6.815802E+03 1.000000E+00
225 33756 -6.815868E+03 1.000000E+00
226 33907 -6.815912E+03 1.000000E+00
227 34058 -6.815976E+03 1.000000E+00
228 34209 -6.816067E+03 1.000000E+00
229 34360 -6.816254E+03 1.000000E+00
230 34511 -6.816770E+03 1.000000E+00
231 34661 -6.817750E+03 1.000000E+00
232 34811 -6.818197E+03 1.000000E+00
233 34961 -6.819067E+03 1.000000E+00
234 35111 -6.819746E+03 1.000000E+00
235 35261 -6.819583E+03 1.000000E+00
236 35411 -6.819739E+03 1.000000E+00
237 35561 -6.819722E+03 1.000000E+00
238 35711 -6.819858E+03 1.000000E+00
239 35861 -6.820100E+03 1.000000E+00
240 36011 -6.820184E+03 1.000000E+00
241 36161 -6.820588E+03 1.000000E+00
242 36311 -6.820485E+03 1.000000E+00
243 36461 -6.820766E+03 1.000000E+00
244 36611 -6.821092E+03 1.000000E+00
245 36761 -6.821173E+03 1.000000E+00
246 36911 -6.821177E+03 1.000000E+00
247 37061 -6.821242E+03 1.000000E+00
248 37211 -6.821274E+03 1.000000E+00
249 37361 -6.821402E+03 1.000000E+00
250 37511 -6.821517E+03 1.000000E+00
251 37661 -6.821595E+03 1.000000E+00
252 37811 -6.821652E+03 1.000000E+00
253 37961 -6.821686E+03 1.000000E+00
254 38111 -6.821715E+03 1.000000E+00
255 38261 -6.821740E+03 1.000000E+00
256 38411 -6.821798E+03 1.000000E+00
257 38561 -6.821895E+03 1.000000E+00
258 38711 -6.822078E+03 1.000000E+00
259 38861 -6.822196E+03 1.000000E+00
260 39011 -6.822381E+03 1.000000E+00
261 39161 -6.822331E+03 1.000000E+00
262 39311 -6.822332E+03 1.000000E+00
263 39461 -6.822336E+03 1.000000E+00
264 39611 -6.822360E+03 1.000000E+00
265 39761 -6.822407E+03 1.000000E+00
266 39911 -6.822513E+03 1.000000E+00
267 40061 -6.822679E+03 1.000000E+00
268 40211 -6.822835E+03 1.000000E+00
269 40361 -6.822891E+03 1.000000E+00
270 40511 -6.823012E+03 1.000000E+00
271 40661 -6.823350E+03 1.000000E+00
272 40811 -6.823827E+03 1.000000E+00
273 40961 -6.824297E+03 1.000000E+00
274 41111 -6.824313E+03 1.000000E+00
275 41261 -6.824342E+03 1.000000E+00
276 41411 -6.824313E+03 1.000000E+00
277 41561 -6.824335E+03 1.000000E+00
278 41711 -6.824355E+03 1.000000E+00
279 41861 -6.824411E+03 1.000000E+00
280 42011 -6.824504E+03 1.000000E+00
281 42161 -6.824676E+03 1.000000E+00
282 42311 -6.824981E+03 1.000000E+00
283 42461 -6.825673E+03 1.000000E+00
284 42611 -6.827455E+03 1.000000E+00
285 42761 -6.828812E+03 1.000000E+00
286 42911 -6.830343E+03 1.000000E+00
287 43061 -6.830492E+03 1.000000E+00
288 43211 -6.830577E+03 1.000000E+00
289 43361 -6.831509E+03 1.000000E+00
290 43511 -6.833406E+03 1.000000E+00
291 43661 -6.833543E+03 1.000000E+00
292 43811 -6.835600E+03 1.000000E+00
293 43961 -6.835874E+03 1.000000E+00
294 44111 -6.836875E+03 1.000000E+00
295 44261 -6.836425E+03 1.000000E+00
296 44411 -6.836495E+03 1.000000E+00
297 44561 -6.836499E+03 1.000000E+00
298 44711 -6.836568E+03 1.000000E+00
299 44861 -6.836731E+03 1.000000E+00
300 45011 -6.837083E+03 1.000000E+00
301 45161 -6.838908E+03 1.000000E+00
302 45311 -6.839566E+03 1.000000E+00
303 45461 -6.840984E+03 1.000000E+00
304 45611 -6.840898E+03 1.000000E+00
305 45761 -6.841545E+03 1.000000E+00
306 45911 -6.841500E+03 1.000000E+00
307 46061 -6.841661E+03 1.000000E+00
308 46211 -6.841700E+03 1.000000E+00
309 46361 -6.841682E+03 1.000000E+00
310 46511 -6.841678E+03 1.000000E+00
311 46661 -6.841682E+03 1.000000E+00
312 46811 -6.841693E+03 1.000000E+00
313 46961 -6.841722E+03 1.000000E+00
314 47111 -6.841797E+03 1.000000E+00
315 47261 -6.841991E+03 1.000000E+00
316 47411 -6.842473E+03 1.000000E+00
317 47561 -6.842692E+03 1.000000E+00
318 47711 -6.843345E+03 1.000000E+00
319 47861 -6.843057E+03 1.000000E+00
320 48011 -6.843049E+03 1.000000E+00
321 48161 -6.843049E+03 1.000000E+00
322 48311 -6.843052E+03 1.000000E+00
323 48461 -6.843059E+03 1.000000E+00
324 48611 -6.843077E+03 1.000000E+00
325 48761 -6.843124E+03 1.000000E+00
326 48911 -6.843247E+03 1.000000E+00
327 49061 -6.843560E+03 1.000000E+00
328 49211 -6.844335E+03 1.000000E+00
329 49361 -6.846075E+03 1.000000E+00
330 49511 -6.847370E+03 1.000000E+00
331 49661 -6.850341E+03 1.000000E+00
332 49811 -6.849566E+03 1.000000E+00
333 49961 -6.849509E+03 1.000000E+00
334 50111 -6.849512E+03 1.000000E+00
335 50261 -6.849524E+03 1.000000E+00
336 50411 -6.849533E+03 1.000000E+00
337 50561 -6.849532E+03 1.000000E+00
338 50711 -6.849530E+03 1.000000E+00
339 50861 -6.849530E+03 1.000000E+00
340 51011 -6.849530E+03 1.000000E+00
341 51161 -6.849530E+03 1.000000E+00
Optimization terminated successfully. (Exit mode 0)
Current function value: -6849.53016555
Iterations: 341
Function evaluations: 51162
Gradient evaluations: 341
NIT FC OBJFUN GNORM
1 150 -6.513791E+03 1.000000E+00
2 300 -6.646865E+03 1.000000E+00
3 450 -6.655354E+03 1.000000E+00
4 600 -6.662407E+03 1.000000E+00
5 750 -6.665427E+03 1.000000E+00
6 900 -6.663142E+03 1.000000E+00
7 1050 -6.655729E+03 1.000000E+00
8 1200 -6.661570E+03 1.000000E+00
9 1350 -6.650444E+03 1.000000E+00
10 1500 -6.632904E+03 1.000000E+00
11 1650 -6.630236E+03 1.000000E+00
12 1800 -6.630822E+03 1.000000E+00
13 1950 -6.631847E+03 1.000000E+00
14 2100 -6.634929E+03 1.000000E+00
15 2250 -6.638360E+03 1.000000E+00
16 2400 -6.641158E+03 1.000000E+00
17 2550 -6.644855E+03 1.000000E+00
18 2700 -6.643832E+03 1.000000E+00
19 2850 -6.644582E+03 1.000000E+00
20 3000 -6.645762E+03 1.000000E+00
21 3150 -6.647120E+03 1.000000E+00
22 3300 -6.650262E+03 1.000000E+00
23 3450 -6.650275E+03 1.000000E+00
24 3600 -6.653782E+03 1.000000E+00
25 3750 -6.653683E+03 1.000000E+00
26 3900 -6.655139E+03 1.000000E+00
27 4050 -6.658597E+03 1.000000E+00
28 4200 -6.664923E+03 1.000000E+00
29 4350 -6.678097E+03 1.000000E+00
30 4500 -6.676005E+03 1.000000E+00
31 4650 -6.677329E+03 1.000000E+00
32 4800 -6.684040E+03 1.000000E+00
33 4950 -6.685596E+03 1.000000E+00
34 5100 -6.696996E+03 1.000000E+00
35 5250 -6.706240E+03 1.000000E+00
36 5400 -6.710042E+03 1.000000E+00
37 5550 -6.713946E+03 1.000000E+00
38 5700 -6.714603E+03 1.000000E+00
39 5850 -6.716178E+03 1.000000E+00
40 6000 -6.717502E+03 1.000000E+00
41 6150 -6.722230E+03 1.000000E+00
42 6300 -6.720588E+03 1.000000E+00
43 6450 -6.722511E+03 1.000000E+00
44 6600 -6.728795E+03 1.000000E+00
45 6750 -6.727254E+03 1.000000E+00
46 6900 -6.728703E+03 1.000000E+00
47 7050 -6.732070E+03 1.000000E+00
48 7200 -6.732116E+03 1.000000E+00
49 7350 -6.735811E+03 1.000000E+00
50 7500 -6.736138E+03 1.000000E+00
51 7650 -6.737898E+03 1.000000E+00
52 7800 -6.740449E+03 1.000000E+00
53 7950 -6.742622E+03 1.000000E+00
54 8100 -6.744615E+03 1.000000E+00
55 8250 -6.746332E+03 1.000000E+00
56 8400 -6.748044E+03 1.000000E+00
57 8550 -6.749506E+03 1.000000E+00
58 8700 -6.749906E+03 1.000000E+00
59 8850 -6.753354E+03 1.000000E+00
60 9000 -6.752916E+03 1.000000E+00
61 9150 -6.753623E+03 1.000000E+00
62 9300 -6.753905E+03 1.000000E+00
63 9450 -6.756034E+03 1.000000E+00
64 9600 -6.756362E+03 1.000000E+00
65 9750 -6.757288E+03 1.000000E+00
66 9900 -6.758082E+03 1.000000E+00
67 10050 -6.759534E+03 1.000000E+00
68 10200 -6.760672E+03 1.000000E+00
69 10350 -6.761394E+03 1.000000E+00
70 10500 -6.761342E+03 1.000000E+00
71 10650 -6.761872E+03 1.000000E+00
72 10800 -6.762053E+03 1.000000E+00
73 10950 -6.762579E+03 1.000000E+00
74 11100 -6.762873E+03 1.000000E+00
75 11250 -6.763093E+03 1.000000E+00
76 11400 -6.763467E+03 1.000000E+00
77 11550 -6.764198E+03 1.000000E+00
78 11700 -6.764490E+03 1.000000E+00
79 11850 -6.765220E+03 1.000000E+00
80 12000 -6.765960E+03 1.000000E+00
81 12150 -6.767597E+03 1.000000E+00
82 12300 -6.768368E+03 1.000000E+00
83 12450 -6.768667E+03 1.000000E+00
84 12600 -6.768929E+03 1.000000E+00
85 12750 -6.769440E+03 1.000000E+00
86 12900 -6.770240E+03 1.000000E+00
87 13050 -6.770645E+03 1.000000E+00
88 13200 -6.771919E+03 1.000000E+00
89 13350 -6.773576E+03 1.000000E+00
90 13500 -6.774261E+03 1.000000E+00
91 13650 -6.778156E+03 1.000000E+00
92 13800 -6.777899E+03 1.000000E+00
93 13950 -6.779291E+03 1.000000E+00
94 14100 -6.780444E+03 1.000000E+00
95 14250 -6.780972E+03 1.000000E+00
96 14400 -6.781314E+03 1.000000E+00
97 14550 -6.782270E+03 1.000000E+00
98 14700 -6.782159E+03 1.000000E+00
99 14850 -6.782516E+03 1.000000E+00
100 15000 -6.782730E+03 1.000000E+00
101 15150 -6.783654E+03 1.000000E+00
102 15300 -6.784961E+03 1.000000E+00
103 15450 -6.788293E+03 1.000000E+00
104 15600 -6.789194E+03 1.000000E+00
105 15750 -6.791212E+03 1.000000E+00
106 15900 -6.790315E+03 1.000000E+00
107 16050 -6.790472E+03 1.000000E+00
108 16200 -6.790650E+03 1.000000E+00
109 16350 -6.790871E+03 1.000000E+00
110 16500 -6.792273E+03 1.000000E+00
111 16650 -6.792213E+03 1.000000E+00
112 16800 -6.792771E+03 1.000000E+00
113 16950 -6.792925E+03 1.000000E+00
114 17100 -6.793025E+03 1.000000E+00
115 17250 -6.793105E+03 1.000000E+00
116 17400 -6.793194E+03 1.000000E+00
117 17550 -6.793391E+03 1.000000E+00
118 17700 -6.793839E+03 1.000000E+00
119 17850 -6.794043E+03 1.000000E+00
120 18000 -6.794344E+03 1.000000E+00
121 18150 -6.794213E+03 1.000000E+00
122 18300 -6.794441E+03 1.000000E+00
123 18450 -6.794696E+03 1.000000E+00
124 18600 -6.794791E+03 1.000000E+00
125 18750 -6.795013E+03 1.000000E+00
126 18900 -6.795085E+03 1.000000E+00
127 19050 -6.795319E+03 1.000000E+00
128 19200 -6.795385E+03 1.000000E+00
129 19350 -6.795514E+03 1.000000E+00
130 19500 -6.795583E+03 1.000000E+00
131 19650 -6.795701E+03 1.000000E+00
132 19800 -6.795759E+03 1.000000E+00
133 19950 -6.796054E+03 1.000000E+00
134 20100 -6.796100E+03 1.000000E+00
135 20250 -6.796332E+03 1.000000E+00
136 20400 -6.796618E+03 1.000000E+00
137 20550 -6.796757E+03 1.000000E+00
138 20700 -6.797034E+03 1.000000E+00
139 20850 -6.797236E+03 1.000000E+00
140 21000 -6.797358E+03 1.000000E+00
141 21150 -6.797504E+03 1.000000E+00
142 21300 -6.797586E+03 1.000000E+00
143 21450 -6.797714E+03 1.000000E+00
144 21600 -6.797726E+03 1.000000E+00
145 21750 -6.797829E+03 1.000000E+00
146 21900 -6.797960E+03 1.000000E+00
147 22050 -6.798216E+03 1.000000E+00
148 22200 -6.798504E+03 1.000000E+00
149 22350 -6.798626E+03 1.000000E+00
150 22500 -6.798720E+03 1.000000E+00
151 22650 -6.799138E+03 1.000000E+00
152 22800 -6.799117E+03 1.000000E+00
153 22950 -6.799214E+03 1.000000E+00
154 23100 -6.799233E+03 1.000000E+00
155 23250 -6.799257E+03 1.000000E+00
156 23400 -6.799329E+03 1.000000E+00
157 23550 -6.799367E+03 1.000000E+00
158 23700 -6.799427E+03 1.000000E+00
159 23850 -6.799449E+03 1.000000E+00
160 24000 -6.799605E+03 1.000000E+00
161 24150 -6.799583E+03 1.000000E+00
162 24300 -6.799684E+03 1.000000E+00
163 24450 -6.799733E+03 1.000000E+00
164 24600 -6.799889E+03 1.000000E+00
165 24750 -6.800089E+03 1.000000E+00
166 24900 -6.800421E+03 1.000000E+00
167 25050 -6.801031E+03 1.000000E+00
168 25200 -6.801343E+03 1.000000E+00
169 25350 -6.801521E+03 1.000000E+00
170 25500 -6.801568E+03 1.000000E+00
171 25650 -6.801629E+03 1.000000E+00
172 25800 -6.801954E+03 1.000000E+00
173 25951 -6.802022E+03 1.000000E+00
174 26101 -6.802342E+03 1.000000E+00
175 26251 -6.802363E+03 1.000000E+00
176 26401 -6.802497E+03 1.000000E+00
177 26551 -6.802536E+03 1.000000E+00
178 26701 -6.802622E+03 1.000000E+00
179 26851 -6.802726E+03 1.000000E+00
180 27002 -6.802816E+03 1.000000E+00
181 27153 -6.802927E+03 1.000000E+00
182 27303 -6.803114E+03 1.000000E+00
183 27453 -6.803156E+03 1.000000E+00
184 27603 -6.803265E+03 1.000000E+00
185 27753 -6.803302E+03 1.000000E+00
186 27903 -6.803345E+03 1.000000E+00
187 28053 -6.803371E+03 1.000000E+00
188 28203 -6.803403E+03 1.000000E+00
189 28353 -6.803506E+03 1.000000E+00
190 28503 -6.803560E+03 1.000000E+00
191 28653 -6.803673E+03 1.000000E+00
192 28803 -6.803709E+03 1.000000E+00
193 28953 -6.803745E+03 1.000000E+00
194 29103 -6.803840E+03 1.000000E+00
195 29253 -6.803956E+03 1.000000E+00
196 29403 -6.804026E+03 1.000000E+00
197 29553 -6.804130E+03 1.000000E+00
198 29703 -6.804140E+03 1.000000E+00
199 29853 -6.804318E+03 1.000000E+00
200 30003 -6.804424E+03 1.000000E+00
201 30154 -6.804494E+03 1.000000E+00
202 30304 -6.804905E+03 1.000000E+00
203 30454 -6.804960E+03 1.000000E+00
204 30604 -6.805146E+03 1.000000E+00
205 30754 -6.805055E+03 1.000000E+00
206 30904 -6.805048E+03 1.000000E+00
207 31054 -6.805134E+03 1.000000E+00
208 31204 -6.805156E+03 1.000000E+00
209 31354 -6.805284E+03 1.000000E+00
210 31504 -6.805398E+03 1.000000E+00
211 31654 -6.805442E+03 1.000000E+00
212 31804 -6.805604E+03 1.000000E+00
213 31954 -6.805763E+03 1.000000E+00
214 32104 -6.805924E+03 1.000000E+00
215 32254 -6.806223E+03 1.000000E+00
216 32405 -6.806405E+03 1.000000E+00
217 32555 -6.806320E+03 1.000000E+00
218 32705 -6.806326E+03 1.000000E+00
219 32855 -6.806511E+03 1.000000E+00
220 33005 -6.806495E+03 1.000000E+00
221 33155 -6.806907E+03 1.000000E+00
222 33305 -6.807346E+03 1.000000E+00
223 33455 -6.807675E+03 1.000000E+00
224 33605 -6.807867E+03 1.000000E+00
225 33755 -6.807933E+03 1.000000E+00
226 33906 -6.807988E+03 1.000000E+00
227 34057 -6.808138E+03 1.000000E+00
228 34208 -6.808384E+03 1.000000E+00
229 34358 -6.808737E+03 1.000000E+00
230 34508 -6.809329E+03 1.000000E+00
231 34658 -6.809436E+03 1.000000E+00
232 34808 -6.809437E+03 1.000000E+00
233 34958 -6.809444E+03 1.000000E+00
234 35108 -6.809497E+03 1.000000E+00
235 35258 -6.809716E+03 1.000000E+00
236 35408 -6.810051E+03 1.000000E+00
237 35558 -6.810384E+03 1.000000E+00
238 35708 -6.812422E+03 1.000000E+00
239 35858 -6.813075E+03 1.000000E+00
240 36008 -6.813812E+03 1.000000E+00
241 36158 -6.813755E+03 1.000000E+00
242 36308 -6.813916E+03 1.000000E+00
243 36458 -6.814231E+03 1.000000E+00
244 36608 -6.814561E+03 1.000000E+00
245 36758 -6.814805E+03 1.000000E+00
246 36908 -6.815267E+03 1.000000E+00
247 37058 -6.815357E+03 1.000000E+00
248 37208 -6.815732E+03 1.000000E+00
249 37358 -6.816694E+03 1.000000E+00
250 37508 -6.817515E+03 1.000000E+00
251 37658 -6.818782E+03 1.000000E+00
252 37808 -6.819052E+03 1.000000E+00
253 37958 -6.821099E+03 1.000000E+00
254 38108 -6.821784E+03 1.000000E+00
255 38258 -6.822210E+03 1.000000E+00
256 38408 -6.822257E+03 1.000000E+00
257 38558 -6.822277E+03 1.000000E+00
258 38708 -6.822346E+03 1.000000E+00
259 38858 -6.822453E+03 1.000000E+00
260 39008 -6.822426E+03 1.000000E+00
261 39158 -6.822661E+03 1.000000E+00
262 39308 -6.823020E+03 1.000000E+00
263 39458 -6.823564E+03 1.000000E+00
264 39608 -6.824933E+03 1.000000E+00
265 39758 -6.826395E+03 1.000000E+00
266 39908 -6.827337E+03 1.000000E+00
267 40058 -6.827360E+03 1.000000E+00
268 40208 -6.828060E+03 1.000000E+00
269 40358 -6.828949E+03 1.000000E+00
270 40508 -6.830300E+03 1.000000E+00
271 40658 -6.830968E+03 1.000000E+00
272 40808 -6.830598E+03 1.000000E+00
273 40958 -6.830615E+03 1.000000E+00
274 41108 -6.830631E+03 1.000000E+00
275 41258 -6.830707E+03 1.000000E+00
276 41408 -6.830801E+03 1.000000E+00
277 41558 -6.831130E+03 1.000000E+00
278 41708 -6.832432E+03 1.000000E+00
279 41858 -6.832549E+03 1.000000E+00
280 42008 -6.833619E+03 1.000000E+00
281 42158 -6.833454E+03 1.000000E+00
282 42308 -6.833406E+03 1.000000E+00
283 42458 -6.833634E+03 1.000000E+00
284 42608 -6.833665E+03 1.000000E+00
285 42758 -6.833908E+03 1.000000E+00
286 42908 -6.834361E+03 1.000000E+00
287 43058 -6.834430E+03 1.000000E+00
288 43208 -6.834951E+03 1.000000E+00
289 43359 -6.835233E+03 1.000000E+00
290 43510 -6.835659E+03 1.000000E+00
291 43660 -6.836993E+03 1.000000E+00
292 43810 -6.837052E+03 1.000000E+00
293 43960 -6.837407E+03 1.000000E+00
294 44110 -6.837356E+03 1.000000E+00
295 44260 -6.837546E+03 1.000000E+00
296 44410 -6.837642E+03 1.000000E+00
297 44560 -6.838096E+03 1.000000E+00
298 44710 -6.838679E+03 1.000000E+00
299 44860 -6.839692E+03 1.000000E+00
300 45010 -6.841190E+03 1.000000E+00
301 45160 -6.842561E+03 1.000000E+00
302 45310 -6.843258E+03 1.000000E+00
303 45460 -6.844295E+03 1.000000E+00
304 45610 -6.845047E+03 1.000000E+00
305 45760 -6.847325E+03 1.000000E+00
306 45910 -6.848522E+03 1.000000E+00
307 46060 -6.849227E+03 1.000000E+00
308 46210 -6.849302E+03 1.000000E+00
309 46360 -6.849386E+03 1.000000E+00
310 46511 -6.849384E+03 1.000000E+00
311 46661 -6.849392E+03 1.000000E+00
312 46811 -6.849415E+03 1.000000E+00
313 46961 -6.849543E+03 1.000000E+00
314 47111 -6.849630E+03 1.000000E+00
315 47261 -6.849839E+03 1.000000E+00
316 47411 -6.850053E+03 1.000000E+00
317 47561 -6.850590E+03 1.000000E+00
318 47711 -6.851944E+03 1.000000E+00
319 47861 -6.854335E+03 1.000000E+00
320 48011 -6.857194E+03 1.000000E+00
321 48161 -6.856897E+03 1.000000E+00
322 48311 -6.856766E+03 1.000000E+00
323 48461 -6.856745E+03 1.000000E+00
324 48611 -6.856842E+03 1.000000E+00
325 48761 -6.856729E+03 1.000000E+00
326 48911 -6.856750E+03 1.000000E+00
327 49061 -6.856770E+03 1.000000E+00
328 49211 -6.856996E+03 1.000000E+00
329 49361 -6.857174E+03 1.000000E+00
330 49512 -6.857574E+03 1.000000E+00
331 49663 -6.858236E+03 1.000000E+00
332 49814 -6.859495E+03 1.000000E+00
333 49964 -6.861289E+03 1.000000E+00
334 50114 -6.861246E+03 1.000000E+00
335 50264 -6.861347E+03 1.000000E+00
336 50415 -6.861313E+03 1.000000E+00
337 50565 -6.861057E+03 1.000000E+00
338 50715 -6.861119E+03 1.000000E+00
339 50865 -6.861190E+03 1.000000E+00
340 51015 -6.861471E+03 1.000000E+00
341 51165 -6.861582E+03 1.000000E+00
342 51315 -6.862227E+03 1.000000E+00
343 51465 -6.862670E+03 1.000000E+00
344 51616 -6.862987E+03 1.000000E+00
345 51766 -6.863913E+03 1.000000E+00
346 51916 -6.864435E+03 1.000000E+00
347 52066 -6.865054E+03 1.000000E+00
348 52216 -6.864785E+03 1.000000E+00
349 52366 -6.864788E+03 1.000000E+00
350 52516 -6.864787E+03 1.000000E+00
351 52666 -6.864796E+03 1.000000E+00
352 52816 -6.864802E+03 1.000000E+00
353 52966 -6.864825E+03 1.000000E+00
354 53116 -6.864852E+03 1.000000E+00
355 53266 -6.864864E+03 1.000000E+00
356 53416 -6.864874E+03 1.000000E+00
357 53566 -6.864870E+03 1.000000E+00
358 53716 -6.864877E+03 1.000000E+00
359 53866 -6.864884E+03 1.000000E+00
360 54016 -6.864910E+03 1.000000E+00
361 54166 -6.864962E+03 1.000000E+00
362 54316 -6.865023E+03 1.000000E+00
363 54466 -6.865095E+03 1.000000E+00
364 54616 -6.865128E+03 1.000000E+00
365 54766 -6.865111E+03 1.000000E+00
366 54916 -6.865120E+03 1.000000E+00
367 55066 -6.865124E+03 1.000000E+00
368 55216 -6.865151E+03 1.000000E+00
369 55366 -6.865191E+03 1.000000E+00
370 55517 -6.865241E+03 1.000000E+00
371 55667 -6.865368E+03 1.000000E+00
372 55817 -6.865519E+03 1.000000E+00
373 55967 -6.865577E+03 1.000000E+00
374 56117 -6.865655E+03 1.000000E+00
375 56268 -6.865706E+03 1.000000E+00
376 56419 -6.865739E+03 1.000000E+00
377 56569 -6.865783E+03 1.000000E+00
378 56719 -6.865754E+03 1.000000E+00
379 56869 -6.865753E+03 1.000000E+00
380 57019 -6.865752E+03 1.000000E+00
381 57169 -6.865752E+03 1.000000E+00
382 57319 -6.865752E+03 1.000000E+00
383 57469 -6.865752E+03 1.000000E+00
384 57619 -6.865752E+03 1.000000E+00
385 57769 -6.865752E+03 1.000000E+00
386 57919 -6.865753E+03 1.000000E+00
387 58069 -6.865754E+03 1.000000E+00
388 58219 -6.865756E+03 1.000000E+00
389 58369 -6.865757E+03 1.000000E+00
390 58519 -6.865759E+03 1.000000E+00
391 58669 -6.865758E+03 1.000000E+00
392 58819 -6.865758E+03 1.000000E+00
393 58969 -6.865758E+03 1.000000E+00
394 59119 -6.865759E+03 1.000000E+00
395 59269 -6.865759E+03 1.000000E+00
396 59419 -6.865760E+03 1.000000E+00
397 59569 -6.865760E+03 1.000000E+00
398 59719 -6.865760E+03 1.000000E+00
399 59869 -6.865760E+03 1.000000E+00
400 60019 -6.865760E+03 1.000000E+00
401 60169 -6.865760E+03 1.000000E+00
402 60319 -6.865761E+03 1.000000E+00
403 60469 -6.865761E+03 1.000000E+00
404 60619 -6.865762E+03 1.000000E+00
405 60769 -6.865762E+03 1.000000E+00
406 60919 -6.865763E+03 1.000000E+00
407 61069 -6.865763E+03 1.000000E+00
408 61219 -6.865763E+03 1.000000E+00
409 61369 -6.865763E+03 1.000000E+00
410 61519 -6.865763E+03 1.000000E+00
411 61669 -6.865764E+03 1.000000E+00
412 61819 -6.865764E+03 1.000000E+00
413 61969 -6.865766E+03 1.000000E+00
414 62119 -6.865770E+03 1.000000E+00
415 62269 -6.865780E+03 1.000000E+00
416 62419 -6.865798E+03 1.000000E+00
417 62569 -6.865817E+03 1.000000E+00
418 62719 -6.865823E+03 1.000000E+00
419 62869 -6.865822E+03 1.000000E+00
420 63019 -6.865821E+03 1.000000E+00
421 63169 -6.865821E+03 1.000000E+00
422 63319 -6.865821E+03 1.000000E+00
423 63469 -6.865821E+03 1.000000E+00
424 63619 -6.865821E+03 1.000000E+00
Optimization terminated successfully. (Exit mode 0)
Current function value: -6865.82063893
Iterations: 424
Function evaluations: 63619
Gradient evaluations: 424
NIT FC OBJFUN GNORM
1 150 -6.679227E+03 1.000000E+00
2 300 -6.658419E+03 1.000000E+00
3 450 -6.637029E+03 1.000000E+00
4 600 -6.608004E+03 1.000000E+00
5 750 -6.605264E+03 1.000000E+00
6 900 -6.611597E+03 1.000000E+00
7 1050 -6.607041E+03 1.000000E+00
8 1200 -6.606563E+03 1.000000E+00
9 1350 -6.606921E+03 1.000000E+00
10 1500 -6.607597E+03 1.000000E+00
11 1650 -6.609162E+03 1.000000E+00
12 1800 -6.613648E+03 1.000000E+00
13 1950 -6.620136E+03 1.000000E+00
14 2100 -6.626873E+03 1.000000E+00
15 2250 -6.625216E+03 1.000000E+00
16 2400 -6.628266E+03 1.000000E+00
17 2550 -6.633064E+03 1.000000E+00
18 2700 -6.640732E+03 1.000000E+00
19 2850 -6.656729E+03 1.000000E+00
20 3000 -6.665625E+03 1.000000E+00
21 3150 -6.672049E+03 1.000000E+00
22 3300 -6.672159E+03 1.000000E+00
23 3450 -6.674680E+03 1.000000E+00
24 3600 -6.676792E+03 1.000000E+00
25 3750 -6.677072E+03 1.000000E+00
26 3900 -6.678027E+03 1.000000E+00
27 4050 -6.679807E+03 1.000000E+00
28 4201 -6.681335E+03 1.000000E+00
29 4351 -6.682735E+03 1.000000E+00
30 4501 -6.686158E+03 1.000000E+00
31 4651 -6.689383E+03 1.000000E+00
32 4801 -6.690814E+03 1.000000E+00
33 4951 -6.694364E+03 1.000000E+00
34 5101 -6.693843E+03 1.000000E+00
35 5251 -6.696002E+03 1.000000E+00
36 5401 -6.704151E+03 1.000000E+00
37 5551 -6.704569E+03 1.000000E+00
38 5701 -6.709076E+03 1.000000E+00
39 5851 -6.716787E+03 1.000000E+00
40 6001 -6.718245E+03 1.000000E+00
41 6151 -6.721828E+03 1.000000E+00
42 6301 -6.727600E+03 1.000000E+00
43 6451 -6.734688E+03 1.000000E+00
44 6601 -6.740024E+03 1.000000E+00
45 6751 -6.743319E+03 1.000000E+00
46 6901 -6.753416E+03 1.000000E+00
47 7051 -6.756452E+03 1.000000E+00
48 7201 -6.765026E+03 1.000000E+00
49 7351 -6.764561E+03 1.000000E+00
50 7501 -6.768556E+03 1.000000E+00
51 7651 -6.771904E+03 1.000000E+00
52 7801 -6.776218E+03 1.000000E+00
53 7951 -6.777861E+03 1.000000E+00
54 8101 -6.780060E+03 1.000000E+00
55 8251 -6.783332E+03 1.000000E+00
56 8401 -6.782723E+03 1.000000E+00
57 8551 -6.783365E+03 1.000000E+00
58 8701 -6.784466E+03 1.000000E+00
59 8851 -6.787267E+03 1.000000E+00
60 9001 -6.787319E+03 1.000000E+00
61 9151 -6.789718E+03 1.000000E+00
62 9302 -6.792585E+03 1.000000E+00
63 9452 -6.794097E+03 1.000000E+00
64 9602 -6.795266E+03 1.000000E+00
65 9752 -6.796309E+03 1.000000E+00
66 9902 -6.798212E+03 1.000000E+00
67 10052 -6.800191E+03 1.000000E+00
68 10202 -6.800754E+03 1.000000E+00
69 10352 -6.801733E+03 1.000000E+00
70 10502 -6.803946E+03 1.000000E+00
71 10652 -6.804230E+03 1.000000E+00
72 10802 -6.807991E+03 1.000000E+00
73 10952 -6.807415E+03 1.000000E+00
74 11102 -6.808127E+03 1.000000E+00
75 11252 -6.811071E+03 1.000000E+00
76 11402 -6.812607E+03 1.000000E+00
77 11552 -6.815197E+03 1.000000E+00
78 11702 -6.820402E+03 1.000000E+00
79 11852 -6.824270E+03 1.000000E+00
80 12002 -6.826451E+03 1.000000E+00
81 12152 -6.827767E+03 1.000000E+00
82 12302 -6.829730E+03 1.000000E+00
83 12452 -6.833825E+03 1.000000E+00
84 12602 -6.834199E+03 1.000000E+00
85 12752 -6.836086E+03 1.000000E+00
86 12902 -6.836360E+03 1.000000E+00
87 13052 -6.837492E+03 1.000000E+00
88 13202 -6.838283E+03 1.000000E+00
89 13352 -6.839383E+03 1.000000E+00
90 13502 -6.840497E+03 1.000000E+00
91 13652 -6.842266E+03 1.000000E+00
92 13802 -6.842982E+03 1.000000E+00
93 13952 -6.844078E+03 1.000000E+00
94 14102 -6.844349E+03 1.000000E+00
95 14252 -6.844489E+03 1.000000E+00
96 14402 -6.844674E+03 1.000000E+00
97 14552 -6.844791E+03 1.000000E+00
98 14702 -6.845123E+03 1.000000E+00
99 14852 -6.845834E+03 1.000000E+00
100 15002 -6.846508E+03 1.000000E+00
101 15152 -6.846811E+03 1.000000E+00
102 15302 -6.847650E+03 1.000000E+00
103 15452 -6.847516E+03 1.000000E+00
104 15602 -6.847655E+03 1.000000E+00
105 15752 -6.847749E+03 1.000000E+00
106 15902 -6.848131E+03 1.000000E+00
107 16052 -6.848372E+03 1.000000E+00
108 16202 -6.848640E+03 1.000000E+00
109 16352 -6.850034E+03 1.000000E+00
110 16502 -6.850858E+03 1.000000E+00
111 16652 -6.851949E+03 1.000000E+00
112 16802 -6.853528E+03 1.000000E+00
113 16952 -6.853195E+03 1.000000E+00
114 17102 -6.853199E+03 1.000000E+00
115 17252 -6.853228E+03 1.000000E+00
116 17402 -6.853269E+03 1.000000E+00
117 17552 -6.853333E+03 1.000000E+00
118 17702 -6.853337E+03 1.000000E+00
119 17852 -6.853364E+03 1.000000E+00
120 18002 -6.853384E+03 1.000000E+00
121 18152 -6.853450E+03 1.000000E+00
122 18302 -6.853504E+03 1.000000E+00
123 18452 -6.853857E+03 1.000000E+00
124 18603 -6.854128E+03 1.000000E+00
125 18753 -6.854401E+03 1.000000E+00
126 18903 -6.854350E+03 1.000000E+00
127 19053 -6.854452E+03 1.000000E+00
128 19203 -6.854575E+03 1.000000E+00
129 19353 -6.854853E+03 1.000000E+00
130 19503 -6.854931E+03 1.000000E+00
131 19653 -6.855079E+03 1.000000E+00
132 19803 -6.855193E+03 1.000000E+00
133 19953 -6.855285E+03 1.000000E+00
134 20103 -6.855524E+03 1.000000E+00
135 20253 -6.855540E+03 1.000000E+00
136 20403 -6.855735E+03 1.000000E+00
137 20553 -6.855832E+03 1.000000E+00
138 20703 -6.856208E+03 1.000000E+00
139 20853 -6.856206E+03 1.000000E+00
140 21003 -6.856325E+03 1.000000E+00
141 21153 -6.856400E+03 1.000000E+00
142 21303 -6.856559E+03 1.000000E+00
143 21453 -6.856802E+03 1.000000E+00
144 21603 -6.857028E+03 1.000000E+00
145 21753 -6.857169E+03 1.000000E+00
146 21903 -6.857454E+03 1.000000E+00
147 22053 -6.857501E+03 1.000000E+00
148 22203 -6.857613E+03 1.000000E+00
149 22353 -6.857594E+03 1.000000E+00
150 22503 -6.857654E+03 1.000000E+00
151 22653 -6.857760E+03 1.000000E+00
152 22803 -6.857822E+03 1.000000E+00
153 22953 -6.857875E+03 1.000000E+00
154 23103 -6.857881E+03 1.000000E+00
155 23253 -6.857919E+03 1.000000E+00
156 23403 -6.857991E+03 1.000000E+00
157 23553 -6.858146E+03 1.000000E+00
158 23703 -6.858387E+03 1.000000E+00
159 23853 -6.858457E+03 1.000000E+00
160 24003 -6.858655E+03 1.000000E+00
161 24153 -6.858708E+03 1.000000E+00
162 24303 -6.858906E+03 1.000000E+00
163 24453 -6.859160E+03 1.000000E+00
164 24603 -6.859237E+03 1.000000E+00
165 24753 -6.859307E+03 1.000000E+00
166 24903 -6.859309E+03 1.000000E+00
167 25053 -6.859333E+03 1.000000E+00
168 25203 -6.859386E+03 1.000000E+00
169 25353 -6.859429E+03 1.000000E+00
170 25503 -6.859617E+03 1.000000E+00
171 25653 -6.859681E+03 1.000000E+00
172 25803 -6.859839E+03 1.000000E+00
173 25953 -6.859867E+03 1.000000E+00
174 26103 -6.859902E+03 1.000000E+00
175 26253 -6.859897E+03 1.000000E+00
176 26403 -6.859907E+03 1.000000E+00
177 26553 -6.859904E+03 1.000000E+00
178 26703 -6.859905E+03 1.000000E+00
179 26853 -6.859909E+03 1.000000E+00
180 27003 -6.859918E+03 1.000000E+00
181 27153 -6.859943E+03 1.000000E+00
182 27303 -6.859957E+03 1.000000E+00
183 27453 -6.860010E+03 1.000000E+00
184 27603 -6.860073E+03 1.000000E+00
185 27753 -6.860071E+03 1.000000E+00
186 27903 -6.860147E+03 1.000000E+00
187 28053 -6.860206E+03 1.000000E+00
188 28203 -6.860222E+03 1.000000E+00
189 28353 -6.860216E+03 1.000000E+00
190 28503 -6.860224E+03 1.000000E+00
191 28653 -6.860234E+03 1.000000E+00
192 28803 -6.860267E+03 1.000000E+00
193 28953 -6.860339E+03 1.000000E+00
194 29103 -6.860461E+03 1.000000E+00
195 29253 -6.860528E+03 1.000000E+00
196 29403 -6.860592E+03 1.000000E+00
197 29553 -6.860565E+03 1.000000E+00
198 29703 -6.860568E+03 1.000000E+00
199 29853 -6.860573E+03 1.000000E+00
200 30003 -6.860590E+03 1.000000E+00
201 30153 -6.860625E+03 1.000000E+00
202 30303 -6.860691E+03 1.000000E+00
203 30453 -6.860760E+03 1.000000E+00
204 30603 -6.860782E+03 1.000000E+00
205 30753 -6.860772E+03 1.000000E+00
206 30903 -6.860769E+03 1.000000E+00
207 31053 -6.860770E+03 1.000000E+00
208 31203 -6.860771E+03 1.000000E+00
209 31353 -6.860772E+03 1.000000E+00
210 31503 -6.860773E+03 1.000000E+00
211 31653 -6.860772E+03 1.000000E+00
212 31803 -6.860774E+03 1.000000E+00
213 31953 -6.860776E+03 1.000000E+00
214 32103 -6.860777E+03 1.000000E+00
215 32253 -6.860783E+03 1.000000E+00
216 32403 -6.860785E+03 1.000000E+00
217 32553 -6.860786E+03 1.000000E+00
218 32703 -6.860795E+03 1.000000E+00
219 32853 -6.860812E+03 1.000000E+00
220 33003 -6.860858E+03 1.000000E+00
221 33153 -6.860954E+03 1.000000E+00
222 33303 -6.860998E+03 1.000000E+00
223 33453 -6.861142E+03 1.000000E+00
224 33603 -6.861077E+03 1.000000E+00
225 33753 -6.861085E+03 1.000000E+00
226 33903 -6.861089E+03 1.000000E+00
227 34053 -6.861091E+03 1.000000E+00
228 34203 -6.861092E+03 1.000000E+00
229 34353 -6.861095E+03 1.000000E+00
230 34503 -6.861100E+03 1.000000E+00
231 34653 -6.861112E+03 1.000000E+00
232 34803 -6.861137E+03 1.000000E+00
233 34953 -6.861193E+03 1.000000E+00
234 35103 -6.861265E+03 1.000000E+00
235 35253 -6.861352E+03 1.000000E+00
236 35403 -6.861383E+03 1.000000E+00
237 35553 -6.861380E+03 1.000000E+00
238 35703 -6.861377E+03 1.000000E+00
239 35853 -6.861377E+03 1.000000E+00
240 36003 -6.861377E+03 1.000000E+00
241 36153 -6.861379E+03 1.000000E+00
242 36303 -6.861380E+03 1.000000E+00
243 36453 -6.861381E+03 1.000000E+00
244 36603 -6.861381E+03 1.000000E+00
245 36753 -6.861380E+03 1.000000E+00
246 36903 -6.861381E+03 1.000000E+00
247 37053 -6.861381E+03 1.000000E+00
248 37203 -6.861382E+03 1.000000E+00
249 37353 -6.861387E+03 1.000000E+00
250 37503 -6.861414E+03 1.000000E+00
251 37654 -6.861473E+03 1.000000E+00
252 37805 -6.861534E+03 1.000000E+00
253 37955 -6.861650E+03 1.000000E+00
254 38105 -6.861722E+03 1.000000E+00
255 38256 -6.861737E+03 1.000000E+00
256 38407 -6.861753E+03 1.000000E+00
257 38557 -6.861664E+03 1.000000E+00
258 38708 -6.861663E+03 1.000000E+00
259 38859 -6.861668E+03 1.000000E+00
260 39010 -6.861662E+03 1.000000E+00
261 39160 -6.861649E+03 1.000000E+00
262 39310 -6.861649E+03 1.000000E+00
263 39460 -6.861648E+03 1.000000E+00
264 39610 -6.861649E+03 1.000000E+00
265 39760 -6.861649E+03 1.000000E+00
266 39910 -6.861649E+03 1.000000E+00
267 40060 -6.861649E+03 1.000000E+00
268 40210 -6.861649E+03 1.000000E+00
269 40360 -6.861649E+03 1.000000E+00
270 40510 -6.861649E+03 1.000000E+00
271 40660 -6.861649E+03 1.000000E+00
272 40810 -6.861649E+03 1.000000E+00
273 40960 -6.861649E+03 1.000000E+00
274 41110 -6.861649E+03 1.000000E+00
275 41260 -6.861649E+03 1.000000E+00
276 41410 -6.861649E+03 1.000000E+00
277 41560 -6.861649E+03 1.000000E+00
278 41710 -6.861649E+03 1.000000E+00
279 41861 -6.861649E+03 1.000000E+00
280 42011 -6.861649E+03 1.000000E+00
281 42161 -6.861650E+03 1.000000E+00
282 42311 -6.861650E+03 1.000000E+00
283 42461 -6.861650E+03 1.000000E+00
Optimization terminated successfully. (Exit mode 0)
Current function value: -6861.64956083
Iterations: 283
Function evaluations: 42461
Gradient evaluations: 283
In [6]:
tf, sb, cb, s, c = prob.Decode(pop.champion.x)
plt.plot(s[:,0], s[:,1], 'k.-')
plt.axes().set_aspect('equal', 'datalim')
plt.show()
plt.plot(c[:,0], 'k.-')
plt.show()
plt.plot(c[:,1], 'k.-')
plt.show()
In [7]:
# AWESOME! Save it before its lost!
z = pop.champion.x
save("../Data/HSS_10_Mars_Base", z)
In [ ]:
print load("../Data/HSS_20_Mars.npy") # Good, it worked!
In [ ]:
Content source: CISprague/Astro.IQ
Similar notebooks: