In [2]:
import sys
sys.path.append('../')
from Trajectory import *
from Optimisation import *
import matplotlib.pyplot as plt
from PyGMO import *
%matplotlib inline
%config InlineBackend.figure_format = 'svg'
/usr/lib/python2.7/dist-packages/matplotlib/__init__.py:874: UserWarning: axes.color_cycle is deprecated and replaced with axes.prop_cycle; please use the latter.
warnings.warn(self.msg_depr % (key, alt_key))
In [3]:
model = Point_Lander_Drag()
prob = HSS(model, nsegs=40)
print model
State
Dimensions: 5
Initial: [ 0. 5000. 150. -10. 8165.]
Target: [ 0. 0. 0. 0. 5000.]
Lower Bound: [-20000. 0. -200. -200. 10.]
Upper Bound: [ 20000. 5000. 200. 0. 10000.]
Control
Dimensions: 2
Lower Bound: [ 0. 0.]
Upper Bound: [ 1. 3.14159265]
Time
Lower Bound: 1.0
Upper Bound: 1000.0
In [4]:
zguess = prob.Guess.Ballistic(model.si, tf=30)
tf, s = prob.Guess.Ballistic(model.si, tf=30, nlp=False)
plt.plot(s[:,0], s[:,1], 'k.--')
plt.axes().set_aspect('equal', 'datalim')
plt.show()
In [4]:
pop = population(prob)
pop.push_back(zguess)
In [5]:
algo = algorithm.scipy_slsqp(max_iter=3000, screen_output=True)
alg1 = algorithm.mbh(algo, stop=1)
In [6]:
pop = alg1.evolve(pop)
NIT FC OBJFUN GNORM
1 570 -9.489339E+03 1.000000E+00
2 1140 -9.489232E+03 1.000000E+00
3 1710 -9.488880E+03 1.000000E+00
4 2280 -9.488353E+03 1.000000E+00
5 2850 -9.487473E+03 1.000000E+00
6 3420 -9.486179E+03 1.000000E+00
7 3990 -9.484766E+03 1.000000E+00
8 4560 -9.483547E+03 1.000000E+00
9 5130 -9.481875E+03 1.000000E+00
10 5700 -9.479847E+03 1.000000E+00
11 6270 -9.477928E+03 1.000000E+00
12 6840 -9.474904E+03 1.000000E+00
13 7410 -9.469984E+03 1.000000E+00
14 7980 -9.454509E+03 1.000000E+00
15 8550 -9.441824E+03 1.000000E+00
16 9120 -9.429061E+03 1.000000E+00
17 9690 -9.407144E+03 1.000000E+00
18 10260 -9.163434E+03 1.000000E+00
19 10831 -9.067614E+03 1.000000E+00
20 11402 -8.982269E+03 1.000000E+00
21 11973 -8.916626E+03 1.000000E+00
22 12544 -8.862387E+03 1.000000E+00
23 13115 -8.836760E+03 1.000000E+00
24 13685 -8.825485E+03 1.000000E+00
25 14255 -8.847365E+03 1.000000E+00
26 14825 -8.875432E+03 1.000000E+00
27 15396 -8.890137E+03 1.000000E+00
28 15966 -8.915871E+03 1.000000E+00
29 16536 -8.952934E+03 1.000000E+00
30 17106 -8.987721E+03 1.000000E+00
31 17676 -9.012167E+03 1.000000E+00
32 18246 -9.028978E+03 1.000000E+00
33 18816 -9.061629E+03 1.000000E+00
34 19386 -9.075650E+03 1.000000E+00
35 19956 -9.076712E+03 1.000000E+00
36 20526 -9.082591E+03 1.000000E+00
37 21096 -9.085923E+03 1.000000E+00
38 21666 -9.092415E+03 1.000000E+00
39 22236 -9.094289E+03 1.000000E+00
40 22806 -9.100164E+03 1.000000E+00
41 23376 -9.110147E+03 1.000000E+00
42 23946 -9.113354E+03 1.000000E+00
43 24516 -9.118470E+03 1.000000E+00
44 25086 -9.120773E+03 1.000000E+00
45 25656 -9.125319E+03 1.000000E+00
46 26226 -9.131124E+03 1.000000E+00
47 26796 -9.132326E+03 1.000000E+00
48 27366 -9.133651E+03 1.000000E+00
49 27936 -9.134848E+03 1.000000E+00
50 28506 -9.136200E+03 1.000000E+00
51 29076 -9.139224E+03 1.000000E+00
52 29646 -9.140293E+03 1.000000E+00
53 30216 -9.142998E+03 1.000000E+00
54 30786 -9.146332E+03 1.000000E+00
55 31356 -9.146992E+03 1.000000E+00
56 31926 -9.149663E+03 1.000000E+00
57 32496 -9.152568E+03 1.000000E+00
58 33066 -9.152844E+03 1.000000E+00
59 33636 -9.153947E+03 1.000000E+00
60 34206 -9.154715E+03 1.000000E+00
61 34776 -9.156093E+03 1.000000E+00
62 35346 -9.157734E+03 1.000000E+00
63 35916 -9.160212E+03 1.000000E+00
64 36486 -9.162068E+03 1.000000E+00
65 37056 -9.165104E+03 1.000000E+00
66 37626 -9.168803E+03 1.000000E+00
67 38196 -9.169499E+03 1.000000E+00
68 38766 -9.171564E+03 1.000000E+00
69 39336 -9.172395E+03 1.000000E+00
70 39906 -9.173549E+03 1.000000E+00
71 40476 -9.173982E+03 1.000000E+00
72 41046 -9.176069E+03 1.000000E+00
73 41616 -9.177613E+03 1.000000E+00
74 42186 -9.179154E+03 1.000000E+00
75 42756 -9.180659E+03 1.000000E+00
76 43326 -9.181319E+03 1.000000E+00
77 43896 -9.183059E+03 1.000000E+00
78 44466 -9.183450E+03 1.000000E+00
79 45036 -9.183735E+03 1.000000E+00
80 45606 -9.184393E+03 1.000000E+00
81 46176 -9.184950E+03 1.000000E+00
82 46746 -9.185455E+03 1.000000E+00
83 47316 -9.186018E+03 1.000000E+00
84 47886 -9.186195E+03 1.000000E+00
85 48456 -9.186897E+03 1.000000E+00
86 49026 -9.187464E+03 1.000000E+00
87 49596 -9.187812E+03 1.000000E+00
88 50166 -9.188423E+03 1.000000E+00
89 50736 -9.189001E+03 1.000000E+00
90 51306 -9.189673E+03 1.000000E+00
91 51876 -9.189780E+03 1.000000E+00
92 52446 -9.190323E+03 1.000000E+00
93 53016 -9.190689E+03 1.000000E+00
94 53586 -9.190895E+03 1.000000E+00
95 54156 -9.191552E+03 1.000000E+00
96 54726 -9.191715E+03 1.000000E+00
97 55296 -9.192345E+03 1.000000E+00
98 55866 -9.192741E+03 1.000000E+00
99 56436 -9.192924E+03 1.000000E+00
100 57006 -9.193164E+03 1.000000E+00
101 57576 -9.194169E+03 1.000000E+00
102 58146 -9.194299E+03 1.000000E+00
103 58716 -9.194676E+03 1.000000E+00
104 59286 -9.195105E+03 1.000000E+00
105 59856 -9.195579E+03 1.000000E+00
106 60426 -9.196157E+03 1.000000E+00
107 60996 -9.196946E+03 1.000000E+00
108 61566 -9.197303E+03 1.000000E+00
109 62136 -9.198163E+03 1.000000E+00
110 62706 -9.198335E+03 1.000000E+00
111 63276 -9.198909E+03 1.000000E+00
112 63846 -9.199882E+03 1.000000E+00
113 64416 -9.200540E+03 1.000000E+00
114 64986 -9.200671E+03 1.000000E+00
115 65556 -9.202008E+03 1.000000E+00
116 66126 -9.202598E+03 1.000000E+00
117 66696 -9.203778E+03 1.000000E+00
118 67266 -9.204521E+03 1.000000E+00
119 67836 -9.205213E+03 1.000000E+00
120 68406 -9.205368E+03 1.000000E+00
121 68976 -9.205678E+03 1.000000E+00
122 69546 -9.205973E+03 1.000000E+00
123 70116 -9.206260E+03 1.000000E+00
124 70686 -9.206475E+03 1.000000E+00
125 71256 -9.206797E+03 1.000000E+00
126 71826 -9.207207E+03 1.000000E+00
127 72396 -9.207438E+03 1.000000E+00
128 72966 -9.207657E+03 1.000000E+00
129 73536 -9.208352E+03 1.000000E+00
130 74106 -9.208365E+03 1.000000E+00
131 74676 -9.208878E+03 1.000000E+00
132 75246 -9.209034E+03 1.000000E+00
133 75816 -9.209522E+03 1.000000E+00
134 76386 -9.209654E+03 1.000000E+00
135 76956 -9.209873E+03 1.000000E+00
136 77526 -9.210339E+03 1.000000E+00
137 78096 -9.210946E+03 1.000000E+00
138 78666 -9.211179E+03 1.000000E+00
139 79236 -9.211629E+03 1.000000E+00
140 79806 -9.211786E+03 1.000000E+00
141 80376 -9.212382E+03 1.000000E+00
142 80946 -9.212685E+03 1.000000E+00
143 81516 -9.213006E+03 1.000000E+00
144 82086 -9.213532E+03 1.000000E+00
145 82656 -9.213731E+03 1.000000E+00
146 83226 -9.214227E+03 1.000000E+00
147 83796 -9.214656E+03 1.000000E+00
148 84366 -9.215517E+03 1.000000E+00
149 84936 -9.215774E+03 1.000000E+00
150 85506 -9.216280E+03 1.000000E+00
151 86076 -9.216680E+03 1.000000E+00
152 86646 -9.217180E+03 1.000000E+00
153 87216 -9.217420E+03 1.000000E+00
154 87786 -9.217638E+03 1.000000E+00
155 88356 -9.217752E+03 1.000000E+00
156 88926 -9.217793E+03 1.000000E+00
157 89496 -9.217976E+03 1.000000E+00
158 90066 -9.218067E+03 1.000000E+00
159 90636 -9.218293E+03 1.000000E+00
160 91206 -9.218467E+03 1.000000E+00
161 91776 -9.218648E+03 1.000000E+00
162 92346 -9.218898E+03 1.000000E+00
163 92916 -9.218993E+03 1.000000E+00
164 93486 -9.219126E+03 1.000000E+00
165 94056 -9.219232E+03 1.000000E+00
166 94626 -9.219287E+03 1.000000E+00
167 95196 -9.219412E+03 1.000000E+00
168 95766 -9.219527E+03 1.000000E+00
169 96336 -9.219704E+03 1.000000E+00
170 96906 -9.219829E+03 1.000000E+00
171 97476 -9.219833E+03 1.000000E+00
172 98046 -9.219912E+03 1.000000E+00
173 98616 -9.219937E+03 1.000000E+00
174 99186 -9.220027E+03 1.000000E+00
175 99756 -9.220079E+03 1.000000E+00
176 100326 -9.220238E+03 1.000000E+00
177 100896 -9.220274E+03 1.000000E+00
178 101466 -9.220382E+03 1.000000E+00
179 102036 -9.220470E+03 1.000000E+00
180 102606 -9.220557E+03 1.000000E+00
181 103176 -9.220620E+03 1.000000E+00
182 103746 -9.220774E+03 1.000000E+00
183 104316 -9.220786E+03 1.000000E+00
184 104886 -9.220946E+03 1.000000E+00
185 105456 -9.220977E+03 1.000000E+00
186 106026 -9.221000E+03 1.000000E+00
187 106596 -9.221049E+03 1.000000E+00
188 107166 -9.221206E+03 1.000000E+00
189 107736 -9.221245E+03 1.000000E+00
190 108306 -9.221335E+03 1.000000E+00
191 108876 -9.221393E+03 1.000000E+00
192 109446 -9.221450E+03 1.000000E+00
193 110016 -9.221495E+03 1.000000E+00
194 110586 -9.221591E+03 1.000000E+00
195 111156 -9.221665E+03 1.000000E+00
196 111726 -9.221686E+03 1.000000E+00
197 112296 -9.221702E+03 1.000000E+00
198 112866 -9.221735E+03 1.000000E+00
199 113436 -9.221750E+03 1.000000E+00
200 114006 -9.221762E+03 1.000000E+00
201 114576 -9.221837E+03 1.000000E+00
202 115146 -9.221909E+03 1.000000E+00
203 115716 -9.221990E+03 1.000000E+00
204 116286 -9.222193E+03 1.000000E+00
205 116856 -9.222479E+03 1.000000E+00
206 117426 -9.222575E+03 1.000000E+00
207 117996 -9.222700E+03 1.000000E+00
208 118566 -9.222852E+03 1.000000E+00
209 119136 -9.222900E+03 1.000000E+00
210 119706 -9.223006E+03 1.000000E+00
211 120276 -9.223089E+03 1.000000E+00
212 120846 -9.223202E+03 1.000000E+00
213 121416 -9.223236E+03 1.000000E+00
214 121986 -9.223254E+03 1.000000E+00
215 122556 -9.223307E+03 1.000000E+00
216 123126 -9.223358E+03 1.000000E+00
217 123696 -9.223395E+03 1.000000E+00
218 124266 -9.223497E+03 1.000000E+00
219 124836 -9.223560E+03 1.000000E+00
220 125406 -9.223819E+03 1.000000E+00
221 125976 -9.223948E+03 1.000000E+00
222 126546 -9.224015E+03 1.000000E+00
223 127116 -9.224039E+03 1.000000E+00
224 127686 -9.224074E+03 1.000000E+00
225 128256 -9.224119E+03 1.000000E+00
226 128826 -9.224175E+03 1.000000E+00
227 129396 -9.224221E+03 1.000000E+00
228 129966 -9.224288E+03 1.000000E+00
229 130536 -9.224321E+03 1.000000E+00
230 131106 -9.224359E+03 1.000000E+00
231 131676 -9.224414E+03 1.000000E+00
232 132246 -9.224465E+03 1.000000E+00
233 132816 -9.224547E+03 1.000000E+00
234 133386 -9.224707E+03 1.000000E+00
235 133956 -9.224780E+03 1.000000E+00
236 134526 -9.224838E+03 1.000000E+00
237 135096 -9.224905E+03 1.000000E+00
238 135666 -9.224967E+03 1.000000E+00
239 136236 -9.225041E+03 1.000000E+00
240 136806 -9.225081E+03 1.000000E+00
241 137376 -9.225138E+03 1.000000E+00
242 137946 -9.225160E+03 1.000000E+00
243 138516 -9.225175E+03 1.000000E+00
244 139086 -9.225197E+03 1.000000E+00
245 139656 -9.225243E+03 1.000000E+00
246 140226 -9.225250E+03 1.000000E+00
247 140796 -9.225274E+03 1.000000E+00
248 141366 -9.225319E+03 1.000000E+00
249 141936 -9.225324E+03 1.000000E+00
250 142506 -9.225399E+03 1.000000E+00
251 143076 -9.225414E+03 1.000000E+00
252 143646 -9.225465E+03 1.000000E+00
253 144216 -9.225509E+03 1.000000E+00
254 144786 -9.225559E+03 1.000000E+00
255 145356 -9.225575E+03 1.000000E+00
256 145926 -9.225608E+03 1.000000E+00
257 146496 -9.225619E+03 1.000000E+00
258 147066 -9.225662E+03 1.000000E+00
259 147636 -9.225661E+03 1.000000E+00
260 148206 -9.225674E+03 1.000000E+00
261 148776 -9.225674E+03 1.000000E+00
262 149346 -9.225683E+03 1.000000E+00
263 149916 -9.225689E+03 1.000000E+00
264 150486 -9.225696E+03 1.000000E+00
265 151056 -9.225698E+03 1.000000E+00
266 151626 -9.225698E+03 1.000000E+00
267 152196 -9.225699E+03 1.000000E+00
268 152766 -9.225700E+03 1.000000E+00
269 153336 -9.225701E+03 1.000000E+00
270 153906 -9.225703E+03 1.000000E+00
271 154476 -9.225704E+03 1.000000E+00
272 155046 -9.225706E+03 1.000000E+00
273 155616 -9.225707E+03 1.000000E+00
274 156186 -9.225708E+03 1.000000E+00
275 156756 -9.225710E+03 1.000000E+00
276 157326 -9.225715E+03 1.000000E+00
277 157896 -9.225723E+03 1.000000E+00
278 158466 -9.225733E+03 1.000000E+00
279 159036 -9.225752E+03 1.000000E+00
280 159606 -9.225755E+03 1.000000E+00
281 160176 -9.225768E+03 1.000000E+00
282 160746 -9.225773E+03 1.000000E+00
283 161316 -9.225784E+03 1.000000E+00
284 161886 -9.225791E+03 1.000000E+00
285 162456 -9.225796E+03 1.000000E+00
286 163026 -9.225807E+03 1.000000E+00
287 163596 -9.225812E+03 1.000000E+00
288 164166 -9.225816E+03 1.000000E+00
289 164736 -9.225822E+03 1.000000E+00
290 165306 -9.225824E+03 1.000000E+00
291 165876 -9.225829E+03 1.000000E+00
292 166446 -9.225833E+03 1.000000E+00
293 167016 -9.225842E+03 1.000000E+00
294 167586 -9.225839E+03 1.000000E+00
295 168156 -9.225841E+03 1.000000E+00
296 168726 -9.225840E+03 1.000000E+00
297 169296 -9.225841E+03 1.000000E+00
298 169866 -9.225841E+03 1.000000E+00
299 170436 -9.225841E+03 1.000000E+00
300 171006 -9.225841E+03 1.000000E+00
301 171576 -9.225841E+03 1.000000E+00
302 172146 -9.225841E+03 1.000000E+00
303 172716 -9.225841E+03 1.000000E+00
304 173286 -9.225841E+03 1.000000E+00
305 173856 -9.225841E+03 1.000000E+00
306 174426 -9.225841E+03 1.000000E+00
307 174996 -9.225841E+03 1.000000E+00
Optimization terminated successfully. (Exit mode 0)
Current function value: -9225.84142333
Iterations: 307
Function evaluations: 174997
Gradient evaluations: 307
NIT FC OBJFUN GNORM
1 570 -8.970596E+03 1.000000E+00
2 1140 -9.189939E+03 1.000000E+00
3 1710 -9.194386E+03 1.000000E+00
4 2280 -9.196525E+03 1.000000E+00
5 2850 -9.203566E+03 1.000000E+00
6 3420 -9.195942E+03 1.000000E+00
7 3990 -9.204403E+03 1.000000E+00
8 4560 -9.199990E+03 1.000000E+00
9 5131 -9.196905E+03 1.000000E+00
10 5702 -9.191815E+03 1.000000E+00
11 6272 -9.180099E+03 1.000000E+00
12 6843 -9.178717E+03 1.000000E+00
13 7413 -9.175822E+03 1.000000E+00
14 7983 -9.175617E+03 1.000000E+00
15 8553 -9.175145E+03 1.000000E+00
16 9123 -9.175641E+03 1.000000E+00
17 9693 -9.175320E+03 1.000000E+00
18 10263 -9.175151E+03 1.000000E+00
19 10833 -9.173225E+03 1.000000E+00
20 11403 -9.172777E+03 1.000000E+00
21 11973 -9.172947E+03 1.000000E+00
22 12543 -9.172888E+03 1.000000E+00
23 13113 -9.173298E+03 1.000000E+00
24 13683 -9.174454E+03 1.000000E+00
25 14253 -9.174168E+03 1.000000E+00
26 14823 -9.174737E+03 1.000000E+00
27 15393 -9.176382E+03 1.000000E+00
28 15963 -9.177767E+03 1.000000E+00
29 16533 -9.178196E+03 1.000000E+00
30 17103 -9.179212E+03 1.000000E+00
31 17673 -9.180296E+03 1.000000E+00
32 18243 -9.181468E+03 1.000000E+00
33 18813 -9.182035E+03 1.000000E+00
34 19383 -9.183928E+03 1.000000E+00
35 19953 -9.184190E+03 1.000000E+00
36 20523 -9.184876E+03 1.000000E+00
37 21093 -9.186052E+03 1.000000E+00
38 21663 -9.186167E+03 1.000000E+00
39 22233 -9.186989E+03 1.000000E+00
40 22803 -9.187678E+03 1.000000E+00
41 23373 -9.188516E+03 1.000000E+00
42 23943 -9.189453E+03 1.000000E+00
43 24513 -9.190942E+03 1.000000E+00
44 25083 -9.191441E+03 1.000000E+00
45 25653 -9.192135E+03 1.000000E+00
46 26223 -9.192227E+03 1.000000E+00
47 26793 -9.192743E+03 1.000000E+00
48 27363 -9.193849E+03 1.000000E+00
49 27933 -9.193988E+03 1.000000E+00
50 28503 -9.194574E+03 1.000000E+00
51 29073 -9.195355E+03 1.000000E+00
52 29643 -9.195926E+03 1.000000E+00
53 30213 -9.196567E+03 1.000000E+00
54 30783 -9.197255E+03 1.000000E+00
55 31353 -9.198588E+03 1.000000E+00
56 31923 -9.198686E+03 1.000000E+00
57 32493 -9.199222E+03 1.000000E+00
58 33063 -9.199842E+03 1.000000E+00
59 33634 -9.200526E+03 1.000000E+00
60 34204 -9.201413E+03 1.000000E+00
61 34774 -9.202076E+03 1.000000E+00
62 35344 -9.202844E+03 1.000000E+00
63 35914 -9.203068E+03 1.000000E+00
64 36484 -9.203582E+03 1.000000E+00
65 37054 -9.203886E+03 1.000000E+00
66 37624 -9.204111E+03 1.000000E+00
67 38194 -9.204276E+03 1.000000E+00
68 38764 -9.204592E+03 1.000000E+00
69 39334 -9.205333E+03 1.000000E+00
70 39904 -9.206031E+03 1.000000E+00
71 40474 -9.206518E+03 1.000000E+00
72 41044 -9.207026E+03 1.000000E+00
73 41614 -9.207468E+03 1.000000E+00
74 42184 -9.207522E+03 1.000000E+00
75 42754 -9.207926E+03 1.000000E+00
76 43324 -9.208321E+03 1.000000E+00
77 43894 -9.208424E+03 1.000000E+00
78 44464 -9.208848E+03 1.000000E+00
79 45034 -9.209378E+03 1.000000E+00
80 45604 -9.209562E+03 1.000000E+00
81 46174 -9.210176E+03 1.000000E+00
82 46744 -9.210449E+03 1.000000E+00
83 47314 -9.210851E+03 1.000000E+00
84 47884 -9.211054E+03 1.000000E+00
85 48454 -9.211545E+03 1.000000E+00
86 49024 -9.211843E+03 1.000000E+00
87 49594 -9.212189E+03 1.000000E+00
88 50164 -9.212393E+03 1.000000E+00
89 50734 -9.212478E+03 1.000000E+00
90 51304 -9.212892E+03 1.000000E+00
91 51874 -9.213002E+03 1.000000E+00
92 52444 -9.213262E+03 1.000000E+00
93 53014 -9.213617E+03 1.000000E+00
94 53584 -9.213670E+03 1.000000E+00
95 54154 -9.213837E+03 1.000000E+00
96 54724 -9.214008E+03 1.000000E+00
97 55294 -9.214256E+03 1.000000E+00
98 55864 -9.214506E+03 1.000000E+00
99 56434 -9.214681E+03 1.000000E+00
100 57004 -9.214805E+03 1.000000E+00
101 57574 -9.215059E+03 1.000000E+00
102 58144 -9.215131E+03 1.000000E+00
103 58714 -9.215473E+03 1.000000E+00
104 59284 -9.215506E+03 1.000000E+00
105 59854 -9.215653E+03 1.000000E+00
106 60424 -9.215720E+03 1.000000E+00
107 60994 -9.215895E+03 1.000000E+00
108 61564 -9.215937E+03 1.000000E+00
109 62134 -9.215974E+03 1.000000E+00
110 62704 -9.216088E+03 1.000000E+00
111 63274 -9.216273E+03 1.000000E+00
112 63844 -9.216414E+03 1.000000E+00
113 64414 -9.216477E+03 1.000000E+00
114 64984 -9.216587E+03 1.000000E+00
115 65554 -9.216686E+03 1.000000E+00
116 66124 -9.216741E+03 1.000000E+00
117 66694 -9.216871E+03 1.000000E+00
118 67264 -9.216909E+03 1.000000E+00
119 67834 -9.217112E+03 1.000000E+00
120 68404 -9.217275E+03 1.000000E+00
121 68974 -9.217468E+03 1.000000E+00
122 69544 -9.218094E+03 1.000000E+00
123 70114 -9.218080E+03 1.000000E+00
124 70684 -9.218236E+03 1.000000E+00
125 71254 -9.218217E+03 1.000000E+00
126 71824 -9.218254E+03 1.000000E+00
127 72394 -9.218297E+03 1.000000E+00
128 72964 -9.218368E+03 1.000000E+00
129 73534 -9.218410E+03 1.000000E+00
130 74104 -9.218449E+03 1.000000E+00
131 74674 -9.218553E+03 1.000000E+00
132 75244 -9.218630E+03 1.000000E+00
133 75814 -9.218767E+03 1.000000E+00
134 76384 -9.218978E+03 1.000000E+00
135 76954 -9.219087E+03 1.000000E+00
136 77524 -9.219457E+03 1.000000E+00
137 78094 -9.219672E+03 1.000000E+00
138 78664 -9.219913E+03 1.000000E+00
139 79234 -9.220070E+03 1.000000E+00
140 79804 -9.220224E+03 1.000000E+00
141 80374 -9.220235E+03 1.000000E+00
142 80944 -9.220400E+03 1.000000E+00
143 81514 -9.220453E+03 1.000000E+00
144 82084 -9.220504E+03 1.000000E+00
145 82654 -9.220581E+03 1.000000E+00
146 83224 -9.220621E+03 1.000000E+00
147 83794 -9.220695E+03 1.000000E+00
148 84364 -9.220738E+03 1.000000E+00
149 84934 -9.220801E+03 1.000000E+00
150 85504 -9.220820E+03 1.000000E+00
151 86074 -9.220888E+03 1.000000E+00
152 86644 -9.220915E+03 1.000000E+00
153 87214 -9.220964E+03 1.000000E+00
154 87784 -9.221056E+03 1.000000E+00
155 88354 -9.221101E+03 1.000000E+00
156 88924 -9.221194E+03 1.000000E+00
157 89494 -9.221276E+03 1.000000E+00
158 90064 -9.221353E+03 1.000000E+00
159 90634 -9.221560E+03 1.000000E+00
160 91204 -9.221643E+03 1.000000E+00
161 91774 -9.221765E+03 1.000000E+00
162 92344 -9.221797E+03 1.000000E+00
163 92914 -9.221821E+03 1.000000E+00
164 93484 -9.221881E+03 1.000000E+00
165 94054 -9.221939E+03 1.000000E+00
166 94624 -9.221972E+03 1.000000E+00
167 95194 -9.221992E+03 1.000000E+00
168 95764 -9.222044E+03 1.000000E+00
169 96334 -9.222101E+03 1.000000E+00
170 96904 -9.222222E+03 1.000000E+00
171 97474 -9.222238E+03 1.000000E+00
172 98044 -9.222302E+03 1.000000E+00
173 98614 -9.222287E+03 1.000000E+00
174 99184 -9.222313E+03 1.000000E+00
175 99754 -9.222333E+03 1.000000E+00
176 100324 -9.222417E+03 1.000000E+00
177 100894 -9.222494E+03 1.000000E+00
178 101464 -9.222656E+03 1.000000E+00
179 102034 -9.222797E+03 1.000000E+00
180 102604 -9.223012E+03 1.000000E+00
181 103174 -9.223081E+03 1.000000E+00
182 103744 -9.223143E+03 1.000000E+00
183 104314 -9.223318E+03 1.000000E+00
184 104884 -9.223437E+03 1.000000E+00
185 105454 -9.223619E+03 1.000000E+00
186 106024 -9.223748E+03 1.000000E+00
187 106594 -9.223814E+03 1.000000E+00
188 107164 -9.223811E+03 1.000000E+00
189 107734 -9.223867E+03 1.000000E+00
190 108304 -9.223899E+03 1.000000E+00
191 108874 -9.223956E+03 1.000000E+00
192 109444 -9.223986E+03 1.000000E+00
193 110014 -9.224021E+03 1.000000E+00
194 110584 -9.224088E+03 1.000000E+00
195 111154 -9.224168E+03 1.000000E+00
196 111724 -9.224237E+03 1.000000E+00
197 112294 -9.224340E+03 1.000000E+00
198 112864 -9.224392E+03 1.000000E+00
199 113434 -9.224489E+03 1.000000E+00
200 114004 -9.224588E+03 1.000000E+00
201 114574 -9.224688E+03 1.000000E+00
202 115144 -9.224695E+03 1.000000E+00
203 115714 -9.224737E+03 1.000000E+00
204 116284 -9.224769E+03 1.000000E+00
205 116854 -9.224804E+03 1.000000E+00
206 117424 -9.224865E+03 1.000000E+00
207 117994 -9.224921E+03 1.000000E+00
208 118564 -9.224945E+03 1.000000E+00
209 119134 -9.224977E+03 1.000000E+00
210 119704 -9.224992E+03 1.000000E+00
211 120274 -9.225013E+03 1.000000E+00
212 120844 -9.225102E+03 1.000000E+00
213 121414 -9.225197E+03 1.000000E+00
214 121984 -9.225310E+03 1.000000E+00
215 122554 -9.225368E+03 1.000000E+00
216 123124 -9.225341E+03 1.000000E+00
217 123694 -9.225359E+03 1.000000E+00
218 124264 -9.225370E+03 1.000000E+00
219 124834 -9.225400E+03 1.000000E+00
220 125404 -9.225423E+03 1.000000E+00
221 125974 -9.225446E+03 1.000000E+00
222 126544 -9.225452E+03 1.000000E+00
223 127114 -9.225457E+03 1.000000E+00
224 127684 -9.225460E+03 1.000000E+00
225 128254 -9.225464E+03 1.000000E+00
226 128824 -9.225469E+03 1.000000E+00
227 129394 -9.225469E+03 1.000000E+00
228 129964 -9.225472E+03 1.000000E+00
229 130534 -9.225475E+03 1.000000E+00
230 131104 -9.225483E+03 1.000000E+00
231 131674 -9.225493E+03 1.000000E+00
232 132244 -9.225501E+03 1.000000E+00
233 132814 -9.225515E+03 1.000000E+00
234 133384 -9.225512E+03 1.000000E+00
235 133954 -9.225518E+03 1.000000E+00
236 134524 -9.225519E+03 1.000000E+00
237 135094 -9.225521E+03 1.000000E+00
238 135664 -9.225522E+03 1.000000E+00
239 136234 -9.225523E+03 1.000000E+00
240 136804 -9.225525E+03 1.000000E+00
241 137374 -9.225526E+03 1.000000E+00
242 137944 -9.225527E+03 1.000000E+00
243 138514 -9.225527E+03 1.000000E+00
244 139084 -9.225527E+03 1.000000E+00
245 139654 -9.225528E+03 1.000000E+00
246 140224 -9.225528E+03 1.000000E+00
247 140794 -9.225528E+03 1.000000E+00
248 141364 -9.225528E+03 1.000000E+00
249 141934 -9.225528E+03 1.000000E+00
250 142504 -9.225528E+03 1.000000E+00
251 143074 -9.225528E+03 1.000000E+00
252 143644 -9.225528E+03 1.000000E+00
253 144214 -9.225528E+03 1.000000E+00
254 144784 -9.225528E+03 1.000000E+00
255 145354 -9.225528E+03 1.000000E+00
256 145924 -9.225529E+03 1.000000E+00
257 146494 -9.225529E+03 1.000000E+00
258 147064 -9.225529E+03 1.000000E+00
259 147635 -9.225529E+03 1.000000E+00
260 148205 -9.225529E+03 1.000000E+00
261 148775 -9.225529E+03 1.000000E+00
262 149345 -9.225529E+03 1.000000E+00
263 149915 -9.225529E+03 1.000000E+00
264 150485 -9.225529E+03 1.000000E+00
265 151055 -9.225529E+03 1.000000E+00
266 151625 -9.225529E+03 1.000000E+00
267 152195 -9.225529E+03 1.000000E+00
268 152765 -9.225529E+03 1.000000E+00
269 153335 -9.225529E+03 1.000000E+00
270 153905 -9.225529E+03 1.000000E+00
271 154475 -9.225529E+03 1.000000E+00
272 155045 -9.225529E+03 1.000000E+00
273 155615 -9.225529E+03 1.000000E+00
274 156185 -9.225530E+03 1.000000E+00
275 156755 -9.225530E+03 1.000000E+00
276 157325 -9.225530E+03 1.000000E+00
277 157895 -9.225530E+03 1.000000E+00
278 158465 -9.225530E+03 1.000000E+00
279 159035 -9.225531E+03 1.000000E+00
280 159605 -9.225531E+03 1.000000E+00
281 160175 -9.225531E+03 1.000000E+00
282 160745 -9.225531E+03 1.000000E+00
283 161315 -9.225531E+03 1.000000E+00
284 161885 -9.225531E+03 1.000000E+00
285 162455 -9.225531E+03 1.000000E+00
286 163025 -9.225531E+03 1.000000E+00
287 163595 -9.225531E+03 1.000000E+00
288 164165 -9.225531E+03 1.000000E+00
289 164735 -9.225531E+03 1.000000E+00
290 165305 -9.225531E+03 1.000000E+00
291 165875 -9.225532E+03 1.000000E+00
292 166445 -9.225531E+03 1.000000E+00
293 167015 -9.225531E+03 1.000000E+00
294 167585 -9.225531E+03 1.000000E+00
295 168155 -9.225531E+03 1.000000E+00
296 168725 -9.225531E+03 1.000000E+00
297 169295 -9.225531E+03 1.000000E+00
298 169865 -9.225531E+03 1.000000E+00
299 170435 -9.225531E+03 1.000000E+00
300 171005 -9.225532E+03 1.000000E+00
301 171575 -9.225532E+03 1.000000E+00
302 172145 -9.225532E+03 1.000000E+00
303 172715 -9.225532E+03 1.000000E+00
304 173285 -9.225532E+03 1.000000E+00
305 173855 -9.225532E+03 1.000000E+00
306 174425 -9.225532E+03 1.000000E+00
307 174995 -9.225532E+03 1.000000E+00
Optimization terminated successfully. (Exit mode 0)
Current function value: -9225.53155273
Iterations: 307
Function evaluations: 174995
Gradient evaluations: 307
NIT FC OBJFUN GNORM
1 570 -9.557391E+03 1.000000E+00
2 1140 -9.242828E+03 1.000000E+00
3 1710 -9.239362E+03 1.000000E+00
4 2280 -9.238244E+03 1.000000E+00
5 2850 -9.234085E+03 1.000000E+00
6 3420 -9.232811E+03 1.000000E+00
7 3991 -9.233813E+03 1.000000E+00
8 4562 -9.233784E+03 1.000000E+00
9 5133 -9.230464E+03 1.000000E+00
10 5704 -9.228276E+03 1.000000E+00
11 6274 -9.222656E+03 1.000000E+00
12 6844 -9.222320E+03 1.000000E+00
13 7415 -9.220556E+03 1.000000E+00
14 7986 -9.220782E+03 1.000000E+00
15 8557 -9.219989E+03 1.000000E+00
16 9127 -9.217743E+03 1.000000E+00
17 9697 -9.217360E+03 1.000000E+00
18 10268 -9.217213E+03 1.000000E+00
19 10839 -9.217365E+03 1.000000E+00
20 11410 -9.217328E+03 1.000000E+00
21 11980 -9.216059E+03 1.000000E+00
22 12550 -9.215096E+03 1.000000E+00
23 13120 -9.215044E+03 1.000000E+00
24 13690 -9.215336E+03 1.000000E+00
25 14260 -9.215597E+03 1.000000E+00
26 14830 -9.215822E+03 1.000000E+00
27 15400 -9.216053E+03 1.000000E+00
28 15970 -9.216178E+03 1.000000E+00
29 16540 -9.216627E+03 1.000000E+00
30 17110 -9.216821E+03 1.000000E+00
31 17680 -9.217072E+03 1.000000E+00
32 18250 -9.217424E+03 1.000000E+00
33 18820 -9.217643E+03 1.000000E+00
34 19390 -9.217788E+03 1.000000E+00
35 19960 -9.217892E+03 1.000000E+00
36 20530 -9.218289E+03 1.000000E+00
37 21100 -9.218372E+03 1.000000E+00
38 21670 -9.218517E+03 1.000000E+00
39 22240 -9.218981E+03 1.000000E+00
40 22810 -9.218885E+03 1.000000E+00
41 23380 -9.219110E+03 1.000000E+00
42 23950 -9.219207E+03 1.000000E+00
43 24520 -9.219377E+03 1.000000E+00
44 25090 -9.219457E+03 1.000000E+00
45 25660 -9.219689E+03 1.000000E+00
46 26230 -9.219880E+03 1.000000E+00
47 26800 -9.220103E+03 1.000000E+00
48 27370 -9.220450E+03 1.000000E+00
49 27940 -9.220617E+03 1.000000E+00
50 28510 -9.220712E+03 1.000000E+00
51 29080 -9.220735E+03 1.000000E+00
52 29650 -9.220813E+03 1.000000E+00
53 30220 -9.220904E+03 1.000000E+00
54 30790 -9.220909E+03 1.000000E+00
55 31360 -9.220960E+03 1.000000E+00
56 31930 -9.221009E+03 1.000000E+00
57 32500 -9.221048E+03 1.000000E+00
58 33070 -9.221100E+03 1.000000E+00
59 33640 -9.221187E+03 1.000000E+00
60 34210 -9.221186E+03 1.000000E+00
61 34780 -9.221215E+03 1.000000E+00
62 35350 -9.221215E+03 1.000000E+00
63 35920 -9.221234E+03 1.000000E+00
64 36490 -9.221272E+03 1.000000E+00
65 37060 -9.221274E+03 1.000000E+00
66 37630 -9.221352E+03 1.000000E+00
67 38200 -9.221368E+03 1.000000E+00
68 38770 -9.221395E+03 1.000000E+00
69 39340 -9.221436E+03 1.000000E+00
70 39910 -9.221485E+03 1.000000E+00
71 40480 -9.221535E+03 1.000000E+00
72 41050 -9.221565E+03 1.000000E+00
73 41620 -9.221614E+03 1.000000E+00
74 42190 -9.221700E+03 1.000000E+00
75 42760 -9.221793E+03 1.000000E+00
76 43330 -9.222029E+03 1.000000E+00
77 43900 -9.222050E+03 1.000000E+00
78 44470 -9.222146E+03 1.000000E+00
79 45040 -9.222205E+03 1.000000E+00
80 45610 -9.222259E+03 1.000000E+00
81 46180 -9.222317E+03 1.000000E+00
82 46750 -9.222413E+03 1.000000E+00
83 47320 -9.222577E+03 1.000000E+00
84 47890 -9.222627E+03 1.000000E+00
85 48460 -9.222690E+03 1.000000E+00
86 49030 -9.222827E+03 1.000000E+00
87 49600 -9.222867E+03 1.000000E+00
88 50170 -9.222941E+03 1.000000E+00
89 50740 -9.222997E+03 1.000000E+00
90 51310 -9.223045E+03 1.000000E+00
91 51880 -9.223101E+03 1.000000E+00
92 52450 -9.223096E+03 1.000000E+00
93 53020 -9.223115E+03 1.000000E+00
94 53590 -9.223196E+03 1.000000E+00
95 54160 -9.223230E+03 1.000000E+00
96 54730 -9.223384E+03 1.000000E+00
97 55300 -9.223390E+03 1.000000E+00
98 55870 -9.223497E+03 1.000000E+00
99 56440 -9.223510E+03 1.000000E+00
100 57010 -9.223594E+03 1.000000E+00
101 57580 -9.223717E+03 1.000000E+00
102 58150 -9.223957E+03 1.000000E+00
103 58720 -9.224146E+03 1.000000E+00
104 59290 -9.224217E+03 1.000000E+00
105 59860 -9.224432E+03 1.000000E+00
106 60430 -9.224445E+03 1.000000E+00
107 61000 -9.224540E+03 1.000000E+00
108 61570 -9.224548E+03 1.000000E+00
109 62140 -9.224561E+03 1.000000E+00
110 62710 -9.224565E+03 1.000000E+00
111 63280 -9.224584E+03 1.000000E+00
112 63850 -9.224595E+03 1.000000E+00
113 64420 -9.224623E+03 1.000000E+00
114 64990 -9.224655E+03 1.000000E+00
115 65560 -9.224671E+03 1.000000E+00
116 66130 -9.224701E+03 1.000000E+00
117 66700 -9.224789E+03 1.000000E+00
118 67270 -9.224864E+03 1.000000E+00
119 67840 -9.224998E+03 1.000000E+00
120 68410 -9.225050E+03 1.000000E+00
121 68980 -9.225040E+03 1.000000E+00
122 69550 -9.225047E+03 1.000000E+00
123 70120 -9.225059E+03 1.000000E+00
124 70690 -9.225061E+03 1.000000E+00
125 71260 -9.225068E+03 1.000000E+00
126 71830 -9.225071E+03 1.000000E+00
127 72400 -9.225074E+03 1.000000E+00
128 72970 -9.225080E+03 1.000000E+00
129 73540 -9.225085E+03 1.000000E+00
130 74110 -9.225091E+03 1.000000E+00
131 74680 -9.225096E+03 1.000000E+00
132 75250 -9.225098E+03 1.000000E+00
133 75820 -9.225100E+03 1.000000E+00
134 76390 -9.225105E+03 1.000000E+00
135 76960 -9.225115E+03 1.000000E+00
136 77530 -9.225125E+03 1.000000E+00
137 78100 -9.225135E+03 1.000000E+00
138 78670 -9.225140E+03 1.000000E+00
139 79240 -9.225145E+03 1.000000E+00
140 79810 -9.225151E+03 1.000000E+00
141 80380 -9.225164E+03 1.000000E+00
142 80950 -9.225172E+03 1.000000E+00
143 81520 -9.225179E+03 1.000000E+00
144 82090 -9.225185E+03 1.000000E+00
145 82660 -9.225191E+03 1.000000E+00
146 83230 -9.225213E+03 1.000000E+00
147 83800 -9.225233E+03 1.000000E+00
148 84370 -9.225248E+03 1.000000E+00
149 84940 -9.225293E+03 1.000000E+00
150 85510 -9.225321E+03 1.000000E+00
151 86080 -9.225411E+03 1.000000E+00
152 86650 -9.225500E+03 1.000000E+00
153 87220 -9.225530E+03 1.000000E+00
154 87790 -9.225646E+03 1.000000E+00
155 88360 -9.225632E+03 1.000000E+00
156 88930 -9.225615E+03 1.000000E+00
157 89500 -9.225689E+03 1.000000E+00
158 90070 -9.225783E+03 1.000000E+00
159 90640 -9.225806E+03 1.000000E+00
160 91210 -9.225849E+03 1.000000E+00
161 91780 -9.225851E+03 1.000000E+00
162 92350 -9.225852E+03 1.000000E+00
163 92920 -9.225854E+03 1.000000E+00
164 93490 -9.225856E+03 1.000000E+00
165 94060 -9.225860E+03 1.000000E+00
166 94630 -9.225863E+03 1.000000E+00
167 95200 -9.225866E+03 1.000000E+00
168 95770 -9.225867E+03 1.000000E+00
169 96340 -9.225868E+03 1.000000E+00
170 96910 -9.225868E+03 1.000000E+00
171 97480 -9.225868E+03 1.000000E+00
172 98050 -9.225869E+03 1.000000E+00
173 98620 -9.225870E+03 1.000000E+00
174 99190 -9.225870E+03 1.000000E+00
175 99760 -9.225871E+03 1.000000E+00
176 100330 -9.225872E+03 1.000000E+00
177 100900 -9.225872E+03 1.000000E+00
178 101470 -9.225873E+03 1.000000E+00
179 102040 -9.225874E+03 1.000000E+00
180 102610 -9.225874E+03 1.000000E+00
181 103180 -9.225874E+03 1.000000E+00
182 103750 -9.225874E+03 1.000000E+00
183 104320 -9.225874E+03 1.000000E+00
184 104891 -9.225875E+03 1.000000E+00
185 105462 -9.225877E+03 1.000000E+00
186 106033 -9.225879E+03 1.000000E+00
187 106603 -9.225883E+03 1.000000E+00
188 107173 -9.225885E+03 1.000000E+00
189 107743 -9.225905E+03 1.000000E+00
190 108313 -9.225933E+03 1.000000E+00
191 108883 -9.225945E+03 1.000000E+00
192 109453 -9.225955E+03 1.000000E+00
193 110023 -9.225983E+03 1.000000E+00
194 110593 -9.225998E+03 1.000000E+00
195 111163 -9.226006E+03 1.000000E+00
196 111733 -9.226037E+03 1.000000E+00
197 112303 -9.226039E+03 1.000000E+00
198 112873 -9.226070E+03 1.000000E+00
199 113443 -9.226204E+03 1.000000E+00
200 114013 -9.226263E+03 1.000000E+00
201 114583 -9.226305E+03 1.000000E+00
202 115153 -9.226366E+03 1.000000E+00
203 115723 -9.226387E+03 1.000000E+00
204 116293 -9.226443E+03 1.000000E+00
205 116863 -9.226529E+03 1.000000E+00
206 117433 -9.226643E+03 1.000000E+00
207 118003 -9.226747E+03 1.000000E+00
208 118573 -9.226895E+03 1.000000E+00
209 119143 -9.226922E+03 1.000000E+00
210 119713 -9.226926E+03 1.000000E+00
211 120283 -9.226931E+03 1.000000E+00
212 120853 -9.226937E+03 1.000000E+00
213 121423 -9.226942E+03 1.000000E+00
214 121993 -9.226949E+03 1.000000E+00
215 122563 -9.226953E+03 1.000000E+00
216 123133 -9.226964E+03 1.000000E+00
217 123703 -9.226982E+03 1.000000E+00
218 124273 -9.226987E+03 1.000000E+00
219 124843 -9.227017E+03 1.000000E+00
220 125413 -9.227027E+03 1.000000E+00
221 125983 -9.227041E+03 1.000000E+00
222 126553 -9.227055E+03 1.000000E+00
223 127123 -9.227077E+03 1.000000E+00
224 127693 -9.227146E+03 1.000000E+00
225 128263 -9.227250E+03 1.000000E+00
226 128833 -9.227403E+03 1.000000E+00
227 129403 -9.227519E+03 1.000000E+00
228 129973 -9.227480E+03 1.000000E+00
229 130543 -9.227480E+03 1.000000E+00
230 131113 -9.227476E+03 1.000000E+00
231 131683 -9.227475E+03 1.000000E+00
232 132253 -9.227478E+03 1.000000E+00
233 132823 -9.227480E+03 1.000000E+00
234 133393 -9.227481E+03 1.000000E+00
235 133963 -9.227482E+03 1.000000E+00
236 134533 -9.227482E+03 1.000000E+00
237 135103 -9.227484E+03 1.000000E+00
238 135673 -9.227484E+03 1.000000E+00
239 136243 -9.227485E+03 1.000000E+00
240 136813 -9.227485E+03 1.000000E+00
241 137383 -9.227485E+03 1.000000E+00
242 137953 -9.227485E+03 1.000000E+00
243 138523 -9.227485E+03 1.000000E+00
244 139093 -9.227485E+03 1.000000E+00
245 139663 -9.227485E+03 1.000000E+00
246 140233 -9.227485E+03 1.000000E+00
247 140803 -9.227485E+03 1.000000E+00
248 141373 -9.227486E+03 1.000000E+00
249 141943 -9.227487E+03 1.000000E+00
250 142513 -9.227489E+03 1.000000E+00
251 143083 -9.227491E+03 1.000000E+00
252 143653 -9.227499E+03 1.000000E+00
253 144223 -9.227503E+03 1.000000E+00
254 144793 -9.227503E+03 1.000000E+00
255 145363 -9.227505E+03 1.000000E+00
256 145934 -9.227509E+03 1.000000E+00
257 146505 -9.227512E+03 1.000000E+00
258 147075 -9.227514E+03 1.000000E+00
259 147645 -9.227516E+03 1.000000E+00
260 148215 -9.227517E+03 1.000000E+00
261 148785 -9.227517E+03 1.000000E+00
262 149355 -9.227517E+03 1.000000E+00
263 149925 -9.227517E+03 1.000000E+00
264 150495 -9.227517E+03 1.000000E+00
265 151065 -9.227518E+03 1.000000E+00
266 151635 -9.227518E+03 1.000000E+00
267 152205 -9.227518E+03 1.000000E+00
268 152775 -9.227518E+03 1.000000E+00
269 153345 -9.227518E+03 1.000000E+00
270 153915 -9.227518E+03 1.000000E+00
271 154485 -9.227518E+03 1.000000E+00
272 155055 -9.227518E+03 1.000000E+00
273 155625 -9.227518E+03 1.000000E+00
274 156195 -9.227518E+03 1.000000E+00
Optimization terminated successfully. (Exit mode 0)
Current function value: -9227.51817906
Iterations: 274
Function evaluations: 156196
Gradient evaluations: 274
NIT FC OBJFUN GNORM
1 570 -8.896808E+03 1.000000E+00
2 1140 -9.222950E+03 1.000000E+00
3 1710 -9.226672E+03 1.000000E+00
4 2282 -9.226337E+03 1.000000E+00
5 2852 -9.219706E+03 1.000000E+00
6 3422 -9.221141E+03 1.000000E+00
7 3993 -9.221799E+03 1.000000E+00
8 4564 -9.221455E+03 1.000000E+00
9 5135 -9.221422E+03 1.000000E+00
10 5706 -9.218994E+03 1.000000E+00
11 6277 -9.214148E+03 1.000000E+00
12 6847 -9.208679E+03 1.000000E+00
13 7417 -9.208570E+03 1.000000E+00
14 7987 -9.208859E+03 1.000000E+00
15 8557 -9.208772E+03 1.000000E+00
16 9128 -9.208451E+03 1.000000E+00
17 9699 -9.208430E+03 1.000000E+00
18 10270 -9.208001E+03 1.000000E+00
19 10841 -9.207677E+03 1.000000E+00
20 11411 -9.206264E+03 1.000000E+00
21 11981 -9.205862E+03 1.000000E+00
22 12551 -9.205941E+03 1.000000E+00
23 13121 -9.205929E+03 1.000000E+00
24 13691 -9.205729E+03 1.000000E+00
25 14261 -9.205760E+03 1.000000E+00
26 14831 -9.205953E+03 1.000000E+00
27 15401 -9.206222E+03 1.000000E+00
28 15971 -9.206246E+03 1.000000E+00
29 16541 -9.206543E+03 1.000000E+00
30 17111 -9.206831E+03 1.000000E+00
31 17681 -9.206833E+03 1.000000E+00
32 18251 -9.206996E+03 1.000000E+00
33 18821 -9.207196E+03 1.000000E+00
34 19391 -9.207310E+03 1.000000E+00
35 19961 -9.207496E+03 1.000000E+00
36 20531 -9.207807E+03 1.000000E+00
37 21101 -9.207996E+03 1.000000E+00
38 21671 -9.208296E+03 1.000000E+00
39 22241 -9.208445E+03 1.000000E+00
40 22811 -9.209012E+03 1.000000E+00
41 23381 -9.208988E+03 1.000000E+00
42 23951 -9.209096E+03 1.000000E+00
43 24521 -9.209570E+03 1.000000E+00
44 25091 -9.209553E+03 1.000000E+00
45 25661 -9.209713E+03 1.000000E+00
46 26231 -9.209756E+03 1.000000E+00
47 26801 -9.209845E+03 1.000000E+00
48 27371 -9.209880E+03 1.000000E+00
49 27941 -9.209958E+03 1.000000E+00
50 28511 -9.210063E+03 1.000000E+00
51 29081 -9.210102E+03 1.000000E+00
52 29651 -9.210169E+03 1.000000E+00
53 30221 -9.210251E+03 1.000000E+00
54 30791 -9.210296E+03 1.000000E+00
55 31361 -9.210396E+03 1.000000E+00
56 31931 -9.210519E+03 1.000000E+00
57 32501 -9.210699E+03 1.000000E+00
58 33071 -9.210860E+03 1.000000E+00
59 33641 -9.211056E+03 1.000000E+00
60 34211 -9.211189E+03 1.000000E+00
61 34781 -9.211271E+03 1.000000E+00
62 35351 -9.211356E+03 1.000000E+00
63 35921 -9.211481E+03 1.000000E+00
64 36491 -9.211614E+03 1.000000E+00
65 37061 -9.211808E+03 1.000000E+00
66 37631 -9.211973E+03 1.000000E+00
67 38201 -9.212297E+03 1.000000E+00
68 38771 -9.212706E+03 1.000000E+00
69 39341 -9.213305E+03 1.000000E+00
70 39911 -9.213837E+03 1.000000E+00
71 40481 -9.214226E+03 1.000000E+00
72 41051 -9.215426E+03 1.000000E+00
73 41621 -9.215685E+03 1.000000E+00
74 42191 -9.216141E+03 1.000000E+00
75 42761 -9.216812E+03 1.000000E+00
76 43331 -9.216794E+03 1.000000E+00
77 43901 -9.216945E+03 1.000000E+00
78 44471 -9.217120E+03 1.000000E+00
79 45041 -9.217253E+03 1.000000E+00
80 45611 -9.217484E+03 1.000000E+00
81 46181 -9.217497E+03 1.000000E+00
82 46751 -9.217606E+03 1.000000E+00
83 47321 -9.217620E+03 1.000000E+00
84 47891 -9.217745E+03 1.000000E+00
85 48461 -9.217749E+03 1.000000E+00
86 49031 -9.217840E+03 1.000000E+00
87 49601 -9.217877E+03 1.000000E+00
88 50171 -9.217919E+03 1.000000E+00
89 50741 -9.217965E+03 1.000000E+00
90 51311 -9.217990E+03 1.000000E+00
91 51881 -9.218040E+03 1.000000E+00
92 52451 -9.218071E+03 1.000000E+00
93 53021 -9.218144E+03 1.000000E+00
94 53591 -9.218196E+03 1.000000E+00
95 54161 -9.218310E+03 1.000000E+00
96 54731 -9.218415E+03 1.000000E+00
97 55301 -9.218570E+03 1.000000E+00
98 55871 -9.218619E+03 1.000000E+00
99 56441 -9.218707E+03 1.000000E+00
100 57011 -9.218793E+03 1.000000E+00
101 57581 -9.218809E+03 1.000000E+00
102 58151 -9.218905E+03 1.000000E+00
103 58721 -9.218968E+03 1.000000E+00
104 59291 -9.219114E+03 1.000000E+00
105 59861 -9.219245E+03 1.000000E+00
106 60431 -9.219384E+03 1.000000E+00
107 61001 -9.219409E+03 1.000000E+00
108 61571 -9.219414E+03 1.000000E+00
109 62141 -9.219484E+03 1.000000E+00
110 62711 -9.219539E+03 1.000000E+00
111 63281 -9.219615E+03 1.000000E+00
112 63851 -9.219638E+03 1.000000E+00
113 64421 -9.219685E+03 1.000000E+00
114 64991 -9.219737E+03 1.000000E+00
115 65561 -9.219802E+03 1.000000E+00
116 66131 -9.219850E+03 1.000000E+00
117 66701 -9.219888E+03 1.000000E+00
118 67271 -9.219943E+03 1.000000E+00
119 67841 -9.219989E+03 1.000000E+00
120 68411 -9.220007E+03 1.000000E+00
121 68981 -9.220051E+03 1.000000E+00
122 69551 -9.220110E+03 1.000000E+00
123 70121 -9.220200E+03 1.000000E+00
124 70691 -9.220228E+03 1.000000E+00
125 71261 -9.220281E+03 1.000000E+00
126 71831 -9.220272E+03 1.000000E+00
127 72401 -9.220302E+03 1.000000E+00
128 72971 -9.220336E+03 1.000000E+00
129 73541 -9.220405E+03 1.000000E+00
130 74111 -9.220452E+03 1.000000E+00
131 74681 -9.220536E+03 1.000000E+00
132 75251 -9.220593E+03 1.000000E+00
133 75821 -9.220648E+03 1.000000E+00
134 76391 -9.220719E+03 1.000000E+00
135 76961 -9.220785E+03 1.000000E+00
136 77531 -9.220849E+03 1.000000E+00
137 78101 -9.220904E+03 1.000000E+00
138 78671 -9.221017E+03 1.000000E+00
139 79241 -9.221120E+03 1.000000E+00
140 79811 -9.221251E+03 1.000000E+00
141 80381 -9.221346E+03 1.000000E+00
142 80951 -9.221455E+03 1.000000E+00
143 81521 -9.221606E+03 1.000000E+00
144 82091 -9.221705E+03 1.000000E+00
145 82661 -9.221964E+03 1.000000E+00
146 83231 -9.221929E+03 1.000000E+00
147 83801 -9.222001E+03 1.000000E+00
148 84371 -9.222069E+03 1.000000E+00
149 84941 -9.222120E+03 1.000000E+00
150 85511 -9.222211E+03 1.000000E+00
151 86081 -9.222429E+03 1.000000E+00
152 86651 -9.222577E+03 1.000000E+00
153 87221 -9.222656E+03 1.000000E+00
154 87791 -9.222936E+03 1.000000E+00
155 88361 -9.222995E+03 1.000000E+00
156 88931 -9.223024E+03 1.000000E+00
157 89501 -9.223079E+03 1.000000E+00
158 90071 -9.223105E+03 1.000000E+00
159 90641 -9.223149E+03 1.000000E+00
160 91211 -9.223194E+03 1.000000E+00
161 91781 -9.223231E+03 1.000000E+00
162 92351 -9.223265E+03 1.000000E+00
163 92921 -9.223327E+03 1.000000E+00
164 93491 -9.223404E+03 1.000000E+00
165 94061 -9.223470E+03 1.000000E+00
166 94631 -9.223583E+03 1.000000E+00
167 95201 -9.223816E+03 1.000000E+00
168 95771 -9.223925E+03 1.000000E+00
169 96341 -9.224218E+03 1.000000E+00
170 96911 -9.224315E+03 1.000000E+00
171 97481 -9.224488E+03 1.000000E+00
172 98051 -9.224513E+03 1.000000E+00
173 98621 -9.224585E+03 1.000000E+00
174 99191 -9.224635E+03 1.000000E+00
175 99761 -9.224693E+03 1.000000E+00
176 100331 -9.224775E+03 1.000000E+00
177 100901 -9.224882E+03 1.000000E+00
178 101471 -9.224975E+03 1.000000E+00
179 102041 -9.225101E+03 1.000000E+00
180 102611 -9.225111E+03 1.000000E+00
181 103181 -9.225233E+03 1.000000E+00
182 103751 -9.225291E+03 1.000000E+00
183 104321 -9.225410E+03 1.000000E+00
184 104891 -9.225509E+03 1.000000E+00
185 105461 -9.225604E+03 1.000000E+00
186 106031 -9.225734E+03 1.000000E+00
187 106601 -9.225768E+03 1.000000E+00
188 107171 -9.225796E+03 1.000000E+00
189 107741 -9.225800E+03 1.000000E+00
190 108311 -9.225809E+03 1.000000E+00
191 108881 -9.225838E+03 1.000000E+00
192 109451 -9.225871E+03 1.000000E+00
193 110021 -9.225917E+03 1.000000E+00
194 110591 -9.226004E+03 1.000000E+00
195 111161 -9.226042E+03 1.000000E+00
196 111731 -9.226082E+03 1.000000E+00
197 112301 -9.226136E+03 1.000000E+00
198 112871 -9.226155E+03 1.000000E+00
199 113441 -9.226231E+03 1.000000E+00
200 114011 -9.226275E+03 1.000000E+00
201 114581 -9.226326E+03 1.000000E+00
202 115151 -9.226346E+03 1.000000E+00
203 115721 -9.226374E+03 1.000000E+00
204 116291 -9.226366E+03 1.000000E+00
205 116861 -9.226373E+03 1.000000E+00
206 117431 -9.226373E+03 1.000000E+00
207 118001 -9.226376E+03 1.000000E+00
208 118571 -9.226381E+03 1.000000E+00
209 119141 -9.226391E+03 1.000000E+00
210 119711 -9.226417E+03 1.000000E+00
211 120281 -9.226464E+03 1.000000E+00
212 120852 -9.226529E+03 1.000000E+00
213 121422 -9.226627E+03 1.000000E+00
214 121992 -9.226643E+03 1.000000E+00
215 122562 -9.226643E+03 1.000000E+00
216 123132 -9.226646E+03 1.000000E+00
217 123702 -9.226661E+03 1.000000E+00
218 124272 -9.226684E+03 1.000000E+00
219 124842 -9.226719E+03 1.000000E+00
220 125412 -9.226734E+03 1.000000E+00
221 125982 -9.226741E+03 1.000000E+00
222 126552 -9.226741E+03 1.000000E+00
223 127122 -9.226746E+03 1.000000E+00
224 127692 -9.226751E+03 1.000000E+00
225 128262 -9.226760E+03 1.000000E+00
226 128832 -9.226770E+03 1.000000E+00
227 129402 -9.226781E+03 1.000000E+00
228 129973 -9.226786E+03 1.000000E+00
229 130543 -9.226789E+03 1.000000E+00
230 131113 -9.226800E+03 1.000000E+00
231 131683 -9.226822E+03 1.000000E+00
232 132253 -9.226865E+03 1.000000E+00
233 132823 -9.226900E+03 1.000000E+00
234 133393 -9.226958E+03 1.000000E+00
235 133963 -9.226975E+03 1.000000E+00
236 134533 -9.226980E+03 1.000000E+00
237 135103 -9.226987E+03 1.000000E+00
238 135673 -9.227005E+03 1.000000E+00
239 136243 -9.227123E+03 1.000000E+00
240 136813 -9.227226E+03 1.000000E+00
241 137383 -9.227237E+03 1.000000E+00
242 137953 -9.227247E+03 1.000000E+00
243 138523 -9.227243E+03 1.000000E+00
244 139093 -9.227250E+03 1.000000E+00
245 139663 -9.227253E+03 1.000000E+00
246 140233 -9.227256E+03 1.000000E+00
247 140803 -9.227259E+03 1.000000E+00
248 141373 -9.227267E+03 1.000000E+00
249 141943 -9.227277E+03 1.000000E+00
250 142513 -9.227289E+03 1.000000E+00
251 143083 -9.227296E+03 1.000000E+00
252 143653 -9.227306E+03 1.000000E+00
253 144223 -9.227320E+03 1.000000E+00
254 144793 -9.227347E+03 1.000000E+00
255 145363 -9.227405E+03 1.000000E+00
256 145933 -9.227502E+03 1.000000E+00
257 146504 -9.227548E+03 1.000000E+00
258 147074 -9.227621E+03 1.000000E+00
259 147644 -9.227618E+03 1.000000E+00
260 148214 -9.227633E+03 1.000000E+00
261 148784 -9.227658E+03 1.000000E+00
262 149354 -9.227663E+03 1.000000E+00
263 149924 -9.227662E+03 1.000000E+00
264 150494 -9.227661E+03 1.000000E+00
265 151064 -9.227661E+03 1.000000E+00
266 151634 -9.227661E+03 1.000000E+00
267 152204 -9.227661E+03 1.000000E+00
268 152774 -9.227661E+03 1.000000E+00
269 153344 -9.227662E+03 1.000000E+00
270 153914 -9.227663E+03 1.000000E+00
271 154484 -9.227664E+03 1.000000E+00
272 155054 -9.227666E+03 1.000000E+00
273 155624 -9.227667E+03 1.000000E+00
274 156194 -9.227667E+03 1.000000E+00
275 156764 -9.227668E+03 1.000000E+00
276 157334 -9.227669E+03 1.000000E+00
277 157904 -9.227674E+03 1.000000E+00
278 158474 -9.227682E+03 1.000000E+00
279 159044 -9.227685E+03 1.000000E+00
280 159614 -9.227693E+03 1.000000E+00
281 160184 -9.227690E+03 1.000000E+00
282 160754 -9.227690E+03 1.000000E+00
283 161324 -9.227690E+03 1.000000E+00
284 161894 -9.227690E+03 1.000000E+00
285 162464 -9.227691E+03 1.000000E+00
286 163034 -9.227690E+03 1.000000E+00
287 163604 -9.227691E+03 1.000000E+00
288 164174 -9.227691E+03 1.000000E+00
289 164744 -9.227691E+03 1.000000E+00
290 165314 -9.227692E+03 1.000000E+00
291 165884 -9.227694E+03 1.000000E+00
292 166454 -9.227698E+03 1.000000E+00
293 167024 -9.227701E+03 1.000000E+00
294 167594 -9.227710E+03 1.000000E+00
295 168164 -9.227710E+03 1.000000E+00
296 168734 -9.227712E+03 1.000000E+00
297 169304 -9.227715E+03 1.000000E+00
298 169874 -9.227716E+03 1.000000E+00
299 170444 -9.227716E+03 1.000000E+00
300 171014 -9.227715E+03 1.000000E+00
301 171584 -9.227715E+03 1.000000E+00
302 172154 -9.227715E+03 1.000000E+00
303 172724 -9.227715E+03 1.000000E+00
Optimization terminated successfully. (Exit mode 0)
Current function value: -9227.71534646
Iterations: 303
Function evaluations: 172724
Gradient evaluations: 303
In [18]:
tf, sb, cb, s, c = prob.Decode(pop.champion.x)
In [10]:
z = pop.champion.x
save("../Data/HSS/40Seg/HSS_40_Lunar_Base", z)
In [20]:
plt.close('all')
ax1 = plt.subplot(211)
plt.plot(c[:,0], 'k.-')
plt.ylabel("Throttle")
ax2 = plt.subplot(212, sharex=ax1)
plt.plot(c[:,1], 'k.-')
plt.ylabel('Thrust Angle [rad]')
plt.xlabel("Node Index")
plt.figure()
plt.plot(s[:,0], s[:,1], 'k.-')
plt.axes().set_aspect('equal', 'datalim')
plt.xlabel("Cross-Range [m]")
plt.ylabel("Altitude [m]")
plt.figure()
plt.plot(s[:,2], 'k.-')
plt.plot(s[:,3], 'k.--')
plt.legend(["$v_x$", "$v_y$"], loc="best")
plt.xlabel("Node Index")
plt.ylabel("Velocity [m/s]")
plt.show()
In [ ]:
Content source: CISprague/Astro.IQ
Similar notebooks: