%pylab inline
import numpy as np
import matplotlib.pyplot as plt
from rayopt import (system_from_text, Analysis, GeometricTrace,
FiniteConjugate, ParaxialTrace, system_from_yaml,
system_to_yaml, material)
np.set_printoptions(precision=5)
Populating the interactive namespace from numpy and matplotlib
# Surface Radius Thickness Material Diameter
lithograph_table = """
0 0.0000 0.1615 AIR 0.145
1 -0.79639 0.0418 Silica 0.189
2 -0.31776 0.0013 AIR 0.197
3 0.31768 0.0247 Silica 0.198
4 -0.90519 0.0014 AIR 0.198
5 -1.53140 0.0138 Silica 0.196
6 0.86996 0.1068 AIR 0.199
7 -0.96937 0.0138 Silica 0.171
8 0.30419 0.0567 AIR 0.169
9 -4.54831 0.0179 Silica 0.176
10 0.36328 0.0559 AIR 0.180
11 -0.13317 0.0179 Silica 0.185
12 0.83902 0.0148 AIR 0.240
13 0.0000 0.0581 Silica 0.273
14 -0.24698 0.0014 AIR 0.273
15 0.97362 0.0653 Silica 0.330
16 -0.34739 0.0176 AIR 0.330
17 1.12092 0.0447 Silica 0.342
18 -0.71734 0.0009 AIR 0.342
19 0.82256 0.0323 Silica 0.332
20 -1.18244 0.0009 AIR 0.332
21 0.31346 0.0345 Silica 0.300
22 1.42431 0.1071 AIR 0.296
23 -2.67307 0.0138 Silica 0.187
24 0.15590 0.0569 AIR 0.162
25 -0.21770 0.0171 Silica 0.154
26 1.01511 0.0839 AIR 0.154
27 -0.11720 0.0176 Silica 0.161
28 3.98510 0.0171 AIR 0.183
Stop
29 -0.40426 0.0339 Silica 0.193
30 -0.21818 0.0019 AIR 0.214
31 -7.96203 0.0650 Silica 0.246
32 -0.23180 0.0014 AIR 0.265
33 1.01215 0.0411 Silica 0.289
34 -0.50372 0.0008 AIR 0.289
35 0.40042 0.0441 Silica 0.287
36 3.75600 0.0014 AIR 0.281
37 0.24527 0.0429 Silica 0.267
38 0.49356 0.0013 AIR 0.267
39 0.17401 0.1101 Silica 0.233
40 0.11185 0.0832 AIR 0.139
41 0.07661 0.0263 Silica 0.082
42 0.18719 0.0227 AIR 0.065
43 0.0000 0.0000 AIR 0.036
"""
columns = "number roc distance material diameter"
s = system_from_text(lithograph_table, columns.split(), scale=25.4e-3*25.4,
shifts={"distance": 1}, wavelengths=[248e-9],
description=".56 NA lithography lens Sasaya, "
"Ushida, and Mercado (1998), (Laikin/Lens Design)")
s.object = FiniteConjugate(radius=s[0].radius,
pupil=dict(type="na", na=.56/4, telecentric=True))
s.stop = 29
s.rescale()
s.update()
#s.reverse()
# https://docs.google.com/viewer?url=patentimages.storage.googleapis.com/pdfs/US5805344.pdf
data = """
0 107.954 vacuum
-617.88 30.375 1.5084
-207.0893 0.934 vacuum
201.9739 68.636 1.5084
-416.6217 0.865 vacuum
460.0439 7.061 1.5084
179.6999 15.608 vacuum
-373.0162 6.952 1.5084
249.496 30.983 vacuum
-2591.2236 11.541 1.5084
229.2357 33.165 vacuum
-82.3025 11.524 1.5084
569.8191 9.159 vacuum
5523.6251 36.703 1.5084
-156.82 0.889 vacuum
610.3354 41.168 1.5084
-221.8862 0.883 vacuum
528.5938 26.903 1.5084
-570.2004 0.883 vacuum
423.5775 21.883 1.5084
-1396.2999 0.883 vacuum
203.9075 22.715 1.5084
835.4548 67.972 vacuum
-735.899 8.386 1.5084
104.6386 23.616 vacuum
-184.6683 11.034 1.5084
288.7053 58.171 vacuum
-74.5663 11.343 1.5084
2318.9745 11.371 vacuum
-283.4504 22.211 1.5084
-142.5176 1.323 vacuum
-5670.469 39.484 1.5084
-146.6908 0.883 vacuum
654.7531 37.168 1.5084
-347.7071 0.883 vacuum
254.9142 31.6 1.5084
2133.198 0.883 vacuum
164.8042 27.885 1.5084
349.3775 0.884 vacuum
108.9816 73.045 1.5084
75.6698 54.069 vacuum
46.2841 16.956 1.5084
99.3161 13.169 vacuum
0 0 vacuum
"""
columns = "roc distance material"
s = system_from_text(data, columns.split(),
shifts={"distance": 1}, wavelengths=[248e-9],
description=".56 NA lithography lens Sasaya, "
"Ushida, and Mercado (1998), (US5805344)",
object=FiniteConjugate(radius=11.7/.25,
pupil=dict(type="na", na=.56/4, telecentric=True)))
s[0].radius = s.object.radius
s.stop = 30
s[s.stop].radius = 74.
s.update()
p = ParaxialTrace(s)
p.update_conjugates()
p.resize()
s.update()
#l.get("glass", "")
s[1].refractive_index(s.wavelengths[0])
s.update()
print(s.pupil((0, 1.), stop=None, l=None))
g = GeometricTrace(s)
#g.rays_point((0, 0), nrays=3, distribution="meridional")
g.rays_clipping((0, 1))
#.rays_paraxial(ParaxialTrace(s))
fig, ax = plt.subplots()
s.plot(ax)
g.plot(ax)
(10053.64184158896, array([[-1421.50959, -1421.50959], [ 1421.50959, 1421.50959]]))
a = Analysis(s, refocus_full=True)
/home/rj/work/nist/pyrayopt/rayopt/elements.py:208: RuntimeWarning: invalid value encountered in less_equal good = np.square(y[:, :2]).sum(1) <= self.radius**2
System: .56 NA lithography lens Sasaya, Ushida, and Mercado (1998), (US5805344) Scale: 1.0 mm Wavelengths: 248 nm Fields: 0, 0.7, 1 Object: Radius: 46.8 Pupil: Pupil Distance: 10053.6 Telecentric: True NA: 0.14 Image: Radius: 11.7 Update Radius: True Pupil: Pupil Distance: -3614.91 Update Radius: True Radius: 2044.48 Stop: 30 Elements: # T Distance Rad Curv Diameter Material n nd Vd 0 S 0 inf 93.6 basic/vacuum 1.000 1.000 inf 1 S 107.95 -617.9 123.12 - 1.508 1.508 inf 2 S 30.375 -207.1 130.67 basic/vacuum 1.000 1.000 inf 3 S 0.934 202 130.72 - 1.508 1.508 inf 4 S 68.636 -416.6 118.2 basic/vacuum 1.000 1.000 inf 5 S 0.865 460 117.84 - 1.508 1.508 inf 6 S 7.061 179.7 115.27 basic/vacuum 1.000 1.000 inf 7 S 15.608 -373 111.78 - 1.508 1.508 inf 8 S 6.952 249.5 111.46 basic/vacuum 1.000 1.000 inf 9 S 30.983 -2591 116.29 - 1.508 1.508 inf 10 S 11.541 229.2 117.66 basic/vacuum 1.000 1.000 inf 11 S 33.165 -82.3 132.25 - 1.508 1.508 inf 12 S 11.524 569.8 141.85 basic/vacuum 1.000 1.000 inf 13 S 9.159 5524 154.52 - 1.508 1.508 inf 14 S 36.703 -156.8 187.83 basic/vacuum 1.000 1.000 inf 15 S 0.889 610.3 188.5 - 1.508 1.508 inf 16 S 41.168 -221.9 204.96 basic/vacuum 1.000 1.000 inf 17 S 0.883 528.6 205.08 - 1.508 1.508 inf 18 S 26.903 -570.2 203.95 basic/vacuum 1.000 1.000 inf 19 S 0.883 423.6 203.73 - 1.508 1.508 inf 20 S 21.883 -1396 196.62 basic/vacuum 1.000 1.000 inf 21 S 0.883 203.9 196.12 - 1.508 1.508 inf 22 S 22.715 835.5 180.3 basic/vacuum 1.000 1.000 inf 23 S 67.972 -735.9 116.34 - 1.508 1.508 inf 24 S 8.386 104.6 111.55 basic/vacuum 1.000 1.000 inf 25 S 23.616 -184.7 104.03 - 1.508 1.508 inf 26 S 11.034 288.7 103.79 basic/vacuum 1.000 1.000 inf 27 S 58.171 -74.57 112.55 - 1.508 1.508 inf 28 S 11.343 2319 119.45 basic/vacuum 1.000 1.000 inf 29 S 11.371 -283.5 130.19 - 1.508 1.508 inf 30 S 22.211 -142.5 147.53 basic/vacuum 1.000 1.000 inf 31 S 1.323 -5670 148.87 - 1.508 1.508 inf 32 S 39.484 -146.7 175.66 basic/vacuum 1.000 1.000 inf 33 S 0.883 654.8 176.02 - 1.508 1.508 inf 34 S 37.168 -347.7 182.88 basic/vacuum 1.000 1.000 inf 35 S 0.883 254.9 182.88 - 1.508 1.508 inf 36 S 31.6 2133 175.47 basic/vacuum 1.000 1.000 inf 37 S 0.883 164.8 175.19 - 1.508 1.508 inf 38 S 27.885 349.4 159.43 basic/vacuum 1.000 1.000 inf 39 S 0.884 109 158.88 - 1.508 1.508 inf 40 S 73.045 75.67 92.929 basic/vacuum 1.000 1.000 inf 41 S 54.069 46.28 53.051 - 1.508 1.508 inf 42 S 16.956 99.32 38.21 basic/vacuum 1.000 1.000 inf 43 S 13.168 inf 23.401 basic/vacuum 1.000 1.000 inf lagrange: -6.6172 track length: 878.88 object, image height: [ 46.8 11.70001] front, back focal length (from PP): [-2141.2186 2141.2186] entry, exit pupil height: [ 1421.50959 2044.48155] entry, exit pupil distance: [ 9945.68784 -3601.73712] front, back focal distance: [ 8456.91785 -522.13643] front, back principal distance: [ 10598.13645 -2663.35503] front, back nodal distance: [ 10598.13645 -2663.35503] front, back numerical aperture: [ 0.14 0.56] front, back f number: [ 0.75315 0.52366] front, back working f number: [ 3.57143 0.89286] front, back airy radius: [ 0.00108 0.00027] transverse, angular magnification: [-0.25 0.69529] # T path n axial y axial nu chief y chief nu 0 S 0 1 0 0.1414 -46.8 0.004655 1 S 108 1.508 15.26 0.154 -46.3 -0.03344 2 S 138.3 1 18.36 0.1089 -46.97 0.08187 3 S 139.3 1.508 18.47 0.06239 -46.89 0.1999 4 S 207.9 1 21.3 0.03639 -37.8 0.246 5 S 208.8 1.508 21.34 0.01281 -37.58 0.2876 6 S 215.8 1 21.4 0.07334 -36.24 0.185 7 S 231.4 1.508 22.54 0.1041 -33.35 0.1396 8 S 238.4 1 23.02 0.151 -32.71 0.07295 9 S 269.4 1.508 27.7 0.1564 -30.45 0.06697 10 S 280.9 1 28.89 0.2205 -29.93 0.0005826 11 S 314.1 1.508 36.21 0.4442 -29.92 -0.1842 12 S 325.6 1 39.6 0.4795 -31.32 -0.2122 13 S 334.8 1.508 43.99 0.4754 -33.27 -0.2091 14 S 371.5 1 55.56 0.2953 -38.35 -0.08476 15 S 372.3 1.508 55.82 0.2488 -38.43 -0.05274 16 S 413.5 1 62.61 0.1053 -39.87 0.0386 17 S 414.4 1.508 62.71 0.04503 -39.83 0.07692 18 S 441.3 1 63.51 -0.01159 -38.46 0.1112 19 S 442.2 1.508 63.5 -0.08781 -38.36 0.1573 20 S 464.1 1 62.23 -0.1105 -36.08 0.1704 21 S 465 1.508 62.13 -0.2654 -35.93 0.26 22 S 487.7 1 58.13 -0.23 -32.02 0.2405 23 S 555.6 1.508 42.5 -0.2006 -15.67 0.2297 24 S 564 1 41.38 0.0004357 -14.39 0.1597 25 S 587.6 1.508 41.39 0.1144 -10.62 0.1305 26 S 598.7 1 42.23 0.1888 -9.666 0.1135 27 S 656.8 1.508 53.21 0.5516 -3.064 0.09259 28 S 668.2 1 57.36 0.5641 -2.368 0.09208 29 S 679.6 1.508 63.77 0.6785 -1.321 0.08971 30 S 701.8 1 73.76 0.4154 6.684e-14 0.08971 31 S 703.1 1.508 74.31 0.422 0.1187 0.08972 32 S 742.6 1 85.36 0.1262 2.467 0.08117 33 S 743.5 1.508 85.47 0.05983 2.539 0.0792 34 S 780.6 1 86.95 -0.0673 4.49 0.07263 35 S 781.5 1.508 86.89 -0.2406 4.554 0.06355 36 S 813.1 1 81.85 -0.2211 5.886 0.06495 37 S 814 1.508 81.65 -0.473 5.943 0.04662 38 S 841.9 1 72.91 -0.3669 6.805 0.05652 39 S 842.8 1.508 72.58 -0.7055 6.855 0.02454 40 S 915.8 1 38.42 -0.4473 8.043 0.07858 41 S 969.9 1.508 14.23 -0.6037 12.29 -0.05644 42 S 986.8 1 7.448 -0.5656 11.66 0.003237 43 S 1000 1-0.0003485 -0.5656 11.7 0.003237 # T SA3 CMA3 AST3 PTZ3 DIS3 TACHC TCHC 0 S 0 0 0 0 0 0 0 1 S -0.01355 -0.009239 -0.006302 0.02112 0.0101 0 0 2 S 0.0002727 0.00417 0.06374 -0.063 0.01134 0 0 3 S -0.05334 0.04003 -0.03004 -0.0646 0.07102 0 -0 4 S 3.674e-05 -0.000839 0.01916 -0.03132 0.2776 -0 0 5 S -0.003975 -0.007892 -0.01567 -0.02836 -0.08742 0 0 6 S 0.04742 -0.004095 0.0003536 0.07261 -0.0063 0 -0 7 S-9.176e-05 -0.00195 -0.04144 0.03498 -0.1373 0 0 8 S 0.1267 -0.03029 0.007241 0.05229 -0.01423 0 -0 9 S -0.03963 -0.02392 -0.01444 0.005035 -0.005681 0 0 10 S 0.4655 -0.1746 0.06551 0.05692 -0.04593 0 -0 11 S -0.03897 0.06465 -0.1073 0.1585 -0.08505 -0 0 12 S 2.999 -1.459 0.7102 0.0229 -0.3567 0 -0 13 S -2.5 1.119 -0.5008 -0.002362 0.2252 0 -0 14 S 0.01476 -0.03998 0.1083 -0.0832 -0.06811 -0 0 15 S -1.373 0.5243 -0.2003 -0.02138 0.08465 0 -0 16 S -0.006939 0.008565 -0.01057 -0.0588 0.08563 -0 0 17 S -0.2379 0.03904 -0.006407 -0.02468 0.005102 0 -0 18 S -0.02665 0.03872 -0.05625 -0.02288 0.115 -0 0 19 S -0.029 -0.004327-0.0006456 -0.0308 -0.004692 0 0 20 S -0.09503 0.1203 -0.1523 -0.009344 0.2046 -0 0 21 S -0.01278 0.0003831-1.149e-05 -0.06399 0.001919 0 -0 22 S -0.1499 0.1889 -0.2381 0.01562 0.2805 -0 0 23 S 0.4412 -0.4014 0.3652 0.01773 -0.3484 -0 0 24 S 0.5082 0.02849 0.001597 0.1247 0.00708 0 0 25 S 0.09129 -0.08865 0.08609 0.07065 -0.1522 -0 0 26 S 0.5804 0.1386 0.0331 0.04519 0.01869 0 0 27 S 0.6953 -0.2048 0.06031 0.175 -0.0693 -0 0 28 S 5.657 0.8748 0.1353 0.005626 0.02178 0 0 29 S -1.724 -0.4919 -0.1403 0.04603 -0.02689 0 0 30 S 0.07981 -0.07005 0.06149 -0.09155 0.02639 -0 0 31 S -2.444 -0.5449 -0.1215 0.002301 -0.02657 0 0 32 S -0.9293 0.1312 -0.01853 -0.08894 0.01518 -0 0 33 S -0.4976 -0.1648 -0.0546 -0.01993 -0.02469 0 0 34 S -0.7246 0.1363 -0.02566 -0.03752 0.01189 -0 0 35 S -0.221 -0.0731 -0.02418 -0.05118 -0.02493 0 0 36 S -0.2786 0.1033 -0.03826 0.006116 0.01191 -0 0 37 S 0.07178 0.02642 0.009728 -0.07917 -0.02556 0 0 38 S -0.2565 0.1232 -0.05919 0.03734 0.01049 -0 0 39 S 0.3262 0.1302 0.05198 -0.1197 -0.02704 0 0 40 S -0.01701 -0.05207 -0.1594 0.1724 0.03996 0 0 41 S 0.04477 -0.1102 0.2713 -0.2819 0.02612 -0 0 42 S -0.4758 0.117 -0.02876 0.1314 -0.02523 -0 0 43 S -0 0 -0 0 0 0 -0 0.0007997-0.0009397-0.0003345-0.0001872-9.308e-05 0 0
from scipy import optimize
import asa
var = []
val = []
valmax = []
valmin = []
for i in range(1, len(s) - 2):
var.extend([(i, "distance"), (i, "curvature")])
val.extend([s.get_path(_) for _ in var[-2:]])
valmin.extend([max(1e-3, val[-2] - .5), val[-1] - 5e-5])
valmax.extend([val[-2] + .5, val[-1] + 5e-5])
p0 = np.array(val)
print p0
valmin = np.array(valmin)
valmax = np.array(valmax)
pscale = valmax - valmin
#print pscale
assert np.all(pscale > 0)
assert np.all(valmax > p0)
assert np.all(valmin < p0)
def update(p):
for v, q in zip(var, p):
s.set_path(v, float(q))
g = []
for h in 0, .7, 1.:
gi = GeometricTrace(s)
gi.rays_point((0, h), distribution="radau", nrays=22,
filter=False, clip=False)
g.append(gi)
path = []
def merit(p):
update(p*pscale)
#if not s.edge_thickness() > 0:
# raise asa.CostParameterError
#s.update()
[gi.propagate() for gi in g]
r = sum(gi.rms() for gi in g)
if np.isnan(r):
#print "fail", p, r
raise asa.CostParameterError
#print p, r
path.append(r)
return r
print merit(p0/pscale)
try:
p = optimize.minimize(merit, p0/pscale, method="SLSQP", tol=1e-7,
bounds=zip(valmin/pscale, valmax/pscale))
#r = asa.asa( merit, p0/pscale, valmin/pscale, valmax/pscale,
# full_output=True, limit_generated=100,
# limit_acceptances=10)
#print r
#print asa.asa_codes[r[2]]
finally:
merit(p0/pscale)
plt.plot(path)
[ 1.07954000e+02 -1.61843724e-03 3.03750000e+01 -4.82883471e-03 9.34000000e-01 4.95113478e-03 6.86360000e+01 -2.40025904e-03 8.65000000e-01 2.17370560e-03 7.06100000e+00 5.56483337e-03 1.56080000e+01 -2.68084871e-03 6.95200000e+00 4.00808029e-03 3.09830000e+01 -3.85918066e-04 1.15410000e+01 4.36232227e-03 3.31650000e+01 -1.21502992e-02 1.15240000e+01 1.75494293e-03 9.15900000e+00 1.81040527e-04 3.67030000e+01 -6.37673766e-03 8.89000000e-01 1.63844339e-03 4.11680000e+01 -4.50681475e-03 8.83000000e-01 1.89181182e-03 2.69030000e+01 -1.75376938e-03 8.83000000e-01 2.36084306e-03 2.18830000e+01 -7.16178523e-04 8.83000000e-01 4.90418450e-03 2.27150000e+01 1.19695285e-03 6.79720000e+01 -1.35888213e-03 8.38600000e+00 9.55670278e-03 2.36160000e+01 -5.41511456e-03 1.10340000e+01 3.46373967e-03 5.81710000e+01 -1.34108840e-02 1.13430000e+01 4.31225095e-04 1.13710000e+01 -3.52795410e-03 2.22110000e+01 -7.01667724e-03 1.32300000e+00 -1.76352256e-04 3.94840000e+01 -6.81706010e-03 8.83000000e-01 1.52729327e-03 3.71680000e+01 -2.87598384e-03 8.83000000e-01 3.92288856e-03 3.16000000e+01 4.68779738e-04 8.83000000e-01 6.06780652e-03 2.78850000e+01 2.86223354e-03 8.84000000e-01 9.17586088e-03 7.30450000e+01 1.32153118e-02 5.40690000e+01 2.16056918e-02] 0.00105331034712
[<matplotlib.lines.Line2D at 0x7f16c5fbec50>]
plt.plot(path)
plt.ylim(0, .001)
(0, 0.001)
plt.plot(path)
plt.ylim(0, .001)
(0, 0.001)
print p.x*pscale - p0
merit(p.x)
print s
#print system_to_yaml(s)
#merit(p0/pscale)
[ 1.87233238e-02 -8.83049577e-06 3.50448931e-02 -5.94900580e-06 -6.60420374e-02 1.10389944e-05 -7.49300118e-02 4.10389023e-06 -1.07305414e-01 -3.35707345e-06 -7.75994532e-02 -3.66721451e-06 -1.42467444e-01 -1.86792307e-06 -8.07214150e-02 -3.72255257e-06 -9.31900831e-02 -3.59861334e-06 -6.48860697e-02 -1.63659185e-05 4.28382062e-02 -3.25872144e-05 -2.68664968e-02 6.25795945e-07 -1.55303733e-01 -1.02464434e-05 7.83799029e-04 1.01298789e-05 -2.93769192e-02 5.32718381e-06 -3.50551700e-02 5.77433365e-06 -6.28560655e-02 1.22920142e-05 -5.23589182e-02 -1.19824120e-05 -7.47142353e-02 8.27748927e-06 -4.27149306e-02 -1.50384355e-05 -4.95661527e-02 2.68787748e-06 7.70724482e-03 -1.15984609e-05 -2.50229349e-02 5.47607010e-06 -5.22277921e-02 -5.51934572e-06 -6.47269422e-02 9.42985137e-07 -1.14023958e-01 -1.83649664e-05 3.26499932e-03 1.31477082e-05 7.40118374e-02 -5.00000000e-05 -3.36279791e-01 -2.23858660e-05 6.40261302e-02 8.96283464e-06 -1.11205628e-01 -3.56116256e-05 3.29322790e-02 1.33780697e-05 1.47676203e-02 -1.84228049e-05 1.17117504e-03 -2.74869892e-05 -4.04608659e-02 -3.96679489e-06 -1.73084125e-02 -1.06503275e-05 -2.84122626e-02 -9.34131155e-07 4.06462385e-02 -1.38459155e-05 1.85631377e-02 -6.29366923e-06 9.56119081e-02 2.65568897e-05 2.89896039e-01 5.76723697e-07] System: .56 NA lithography lens Sasaya, Ushida, and Mercado (1998), (US5805344) Scale: 1.0 mm Wavelengths: 248 nm Object: Index: 1 Entrance: 123 dia at 108 Pupil: 2.85e+03 dia at 1.01e+04 Radius: 46.8 NA: 0.14 Telecentric: True Image: Index: 1 Entrance: 38.3 dia at 13.2 Pupil: 4.1e+03 dia at 3.61e+03 Radius: 11.7 NA: 0.493 Stop: 30 Elements: # T Distance Rad Curv Diameter Material n nd Vd 0 S 0 inf 93.6 basic/vacuum 1.000 nan nan 1 S 107.97 -614.5 123.22 - 1.508 nan nan 2 S 30.41 -206.8 130.79 basic/vacuum 1.000 nan nan 3 S 0.86796 201.5 130.84 - 1.508 nan nan 4 S 68.561 -417.3 118.34 basic/vacuum 1.000 nan nan 5 S 0.75769 460.8 117.98 - 1.508 nan nan 6 S 6.9834 179.8 115.41 basic/vacuum 1.000 nan nan 7 S 15.466 -372.8 111.93 - 1.508 nan nan 8 S 6.8713 249.7 111.6 basic/vacuum 1.000 nan nan 9 S 30.89 -2567 116.47 - 1.508 nan nan 10 S 11.476 230.1 117.84 basic/vacuum 1.000 nan nan 11 S 33.208 -82.08 132.48 - 1.508 nan nan 12 S 11.497 569.6 142.1 basic/vacuum 1.000 nan nan 13 S 9.0037 5855 154.8 - 1.508 nan nan 14 S 36.704 -157.1 188.18 basic/vacuum 1.000 nan nan 15 S 0.85962 608.4 188.86 - 1.508 nan nan 16 S 41.133 -222.2 205.36 basic/vacuum 1.000 nan nan 17 S 0.82014 525.2 205.48 - 1.508 nan nan 18 S 26.851 -566.3 204.35 basic/vacuum 1.000 nan nan 19 S 0.80829 422.1 204.13 - 1.508 nan nan 20 S 21.84 -1368 197.01 basic/vacuum 1.000 nan nan 21 S 0.83343 203.8 196.52 - 1.508 nan nan 22 S 22.723 843.6 180.67 basic/vacuum 1.000 nan nan 23 S 67.947 -738.9 116.61 - 1.508 nan nan 24 S 8.3338 104.7 111.82 basic/vacuum 1.000 nan nan 25 S 23.551 -184.7 104.29 - 1.508 nan nan 26 S 10.92 290.2 104.06 basic/vacuum 1.000 nan nan 27 S 58.174 -74.64 112.89 - 1.508 nan nan 28 S 11.417 2623 119.82 basic/vacuum 1.000 nan nan 29 S 11.035 -281.7 130.6 - 1.508 nan nan 30 S 22.275 -142.7 148 basic/vacuum 1.000 nan nan 31 S 1.2118 -4718 149.34 - 1.508 nan nan 32 S 39.517 -147 176.2 basic/vacuum 1.000 nan nan 33 S 0.89777 662.7 176.57 - 1.508 nan nan 34 S 37.169 -344.4 183.43 basic/vacuum 1.000 nan nan 35 S 0.84254 255.2 183.44 - 1.508 nan nan 36 S 31.583 2183 175.99 basic/vacuum 1.000 nan nan 37 S 0.85459 164.8 175.71 - 1.508 nan nan 38 S 27.926 351.1 159.89 basic/vacuum 1.000 nan nan 39 S 0.90256 109.1 159.34 - 1.508 nan nan 40 S 73.141 75.52 93.174 basic/vacuum 1.000 nan nan 41 S 54.359 46.28 53.142 - 1.508 nan nan 42 S 16.956 99.32 38.258 basic/vacuum 1.000 nan nan 43 S 13.168 inf 23.4 basic/vacuum 1.000 nan nan
p
status: 0 success: True njev: 41 nfev: 3457 fun: 0.00033871264573987424 x: array([ 107.97029232, -16.19353204, 30.39420665, -48.30113857, 0.9339493 , 49.49782862, 68.63085772, -24.01758134, 0.85467672, 21.74704929, 7.05602851, 55.65850218, 15.60038803, -26.78017997, 6.95099837, 40.09248531, 30.99259197, -3.80370601, 11.549283 , 43.63505988, 33.19993943, -121.56762543, 11.54744571, 17.4388284 , 8.8697162 , 1.74989762, 36.73179648, -63.77820973, 0.89300205, 16.41591505, 41.17418074, -45.1213454 , 0.87780647, 19.00046995, 26.90059036, -17.62478518, 0.88248396, 23.71890817, 21.89213953, -7.23401653, 0.89850608, 49.13939694, 22.76223936, 11.94616361, 68.02510219, -13.59857342, 8.42233804, 95.56996598, 23.61724531, -54.18783247, 11.02498266, 34.63284223, 58.15237254, -134.10733868, 11.36795587, 4.184563 , 11.29972627, -35.33906897, 22.22028973, -70.11946219, 1.31706886, -1.77121995, 39.50923295, -68.1731107 , 0.88737822, 15.18988446, 37.17736397, -28.73723745, 0.88568158, 39.16511819, 31.60842 , 4.7320775 , 0.89123658, 60.65426365, 27.90956389, 28.66185107, 0.91277657, 91.76562309, 73.12424419, 132.22370596, 54.20252585, 216.04725809]) message: 'Optimization terminated successfully.' jac: array([ 1.36420487e-04, 5.92054712e-05, 1.10151377e-04, -8.47409756e-05, 9.16996578e-05, 1.11087327e-04, 2.81217617e-05, -1.43347344e-04, 1.68420083e-05, 1.42682467e-04, 1.04551946e-05, -1.33381938e-04, 4.98842128e-05, 1.61978340e-04, 5.23517483e-05, -1.67917155e-04, 1.63135046e-04, 2.48984310e-04, 1.21038356e-04, -2.74534388e-04, 3.47987108e-04, 4.16164978e-04, 8.93496428e-04, -6.33944172e-04, 1.56026365e-03, 7.51568012e-04, 1.03365062e-03, -9.85309787e-04, 6.28474099e-04, 1.18083646e-03, 2.85798524e-04, -1.40605871e-03, 7.51266489e-05, 1.38580635e-03, 1.16220399e-05, -1.42155228e-03, 9.06951027e-06, 1.37669278e-03, 3.93488153e-05, -1.35739053e-03, 1.02116654e-04, 1.28074215e-03, 3.24091794e-04, -1.15082497e-03, 3.75375381e-04, 6.16357349e-04, 1.82680204e-04, -5.59161410e-04, 1.43609941e-05, 5.58804612e-04, 6.69576111e-05, -5.92626016e-04, 2.53504099e-04, 9.62561964e-04, 1.37457105e-03, -1.31649641e-03, 2.37718608e-03, 1.44869718e-03, 2.18245330e-03, -1.71481360e-03, 1.35619034e-03, 2.19251964e-03, 8.83656008e-04, -2.65165756e-03, 1.33420879e-04, 2.70702483e-03, 2.09136779e-05, -2.72298325e-03, 3.06478614e-05, 2.75796666e-03, 2.61076311e-04, -2.39340555e-03, 3.31659769e-04, 2.41484997e-03, 1.02955157e-03, -1.90635628e-03, 9.26193206e-04, 1.90016479e-03, 2.30423639e-03, -5.39767592e-04, 1.34413375e-03, 4.32562301e-05, 0.00000000e+00]) nit: 41