spyglass/spyglass_v2_old.ipynb

1.3 MiB

In [1]:
%matplotlib widget

isdark = False

from rayoptics.environment import *
In [2]:
opm = OpticalModel()
sm  = opm['seq_model']
osp = opm['optical_spec']
pm = opm['parax_model']

osp['pupil'] = PupilSpec(osp, key=['object', 'pupil'], value=16)
osp['fov'] = FieldSpec(osp, key=['object', 'angle'], value=0.5, flds=[0., 0.707, 1.], is_relative=True)
osp['wvls'] = WvlSpec([('F', 0.5), (587.5618, 1.0), ('C', 0.5)], ref_wl=1)

opm.radius_mode = True

sm.gaps[0].thi=1e10

def calc_curvature(n, fl):
    return (n-1)*fl

n_bk7 = 1.5168
n_lasf9 = 1.85025
n_f2 =  1.62005

v_bk7 = 64.17
v_lasf9 = 32.16
v_f2 = 36.43
# try for chaining a 3x telescope setup with a second 3x telescope setup
f0 = 150
f0_0 = (v_bk7-v_f2)*f0/v_bk7
f0_1 = -f0_0*v_bk7/v_f2
f0_actual = 1/(1/f0_0+1/f0_1)
f1 = 40
# 1/f = 1/f0 + 1/f1 = 1/f0 - v2/(f0*v1) = (v1-v2)/(v1*f0)
# f0 = (v1-v2)*f/v1
f1_0 = (v_bk7-v_f2)*f1/v_bk7
f1_1 = -f1_0*v_bk7/v_f2
print(f0_actual, f1_0, f1_1)
f2 = 150
f3 = 50

f1_stacked = 2*f1

r0 = calc_curvature(n_lasf9, f0)
r0_0 = calc_curvature(n_bk7, f0_0)
r0_1 = calc_curvature(n_f2, f0_1)
r1 = calc_curvature(n_lasf9, f1)
r1_0 = calc_curvature(n_bk7, f1_0)
r1_1 = calc_curvature(n_f2, f1_1)
r2 = calc_curvature(n_lasf9, f2)
r3 = calc_curvature(n_lasf9, f3)

sm.add_surface([r0_0, 2, 'N-BK7', 'Schott', 42/2])
sm.add_surface([1e9, 2, 'N-F2', 'Schott', 42/2])
#sm.add_surface([-r0_1, 180+2*f1])
sm.add_surface([-r0_1, 180+2*36.])

#sm.add_surface([r1_0, 6, 'N-BK7', 'Schott', 16])
#sm.add_surface([1e9, 2, 'N-F2', 'Schott', 16])
#sm.add_surface([-r1_1, 30])
#sm.add_surface([1e9, 2, 'N-LASF9', 'Schott', 16])
#sm.add_surface([-r1, f2+f3])

#sm.add_surface([r2, 2, 'N-LASF9', 'Schott', 16])
#sm.add_surface([1e9, f2+f3])
#sm.add_surface([1e9, 2, 'N-LASF9', 'Schott', 16])
#sm.add_surface([-r3, f3])
150.0 17.29156926912888 -30.45841339555312
In [3]:
opm.update_model()

layout_plt = plt.figure(FigureClass=InteractiveLayout, opt_model=opm, is_dark=isdark).plot()
Figure
In [4]:
pm.first_order_data()
efl               142.6
ffl              -145.8
pp1              -3.188
bfl                 137
ppk               5.616
f/#               8.914
m            -1.426e-08
red          -7.011e+07
obj_dist          1e+10
obj_ang             0.5
enp_dist             -0
enp_radius            8
na obj            8e-10
n obj                 1
img_dist            137
img_ht            1.245
exp_dist         -117.5
exp_radius        7.825
na img           -0.056
n img                 1
optical invariant      0.06981
In [5]:
# add a focusing mirror
#sm.add_surface([100, 2, 1.85, 16])
#sm.add_surface([1e9, 36])

opm.update_model()

opm.seq_model.gaps[-1].thi = opm.optical_spec.parax_data.fod.bfl - 9
opm.update_model()

layout_plt = plt.figure(FigureClass=InteractiveLayout, opt_model=opm, is_dark=isdark).plot()
Figure
In [6]:
#opm.seq_model.gaps[-1].thi -= 0.1
opm.update_model()
#spot_plt = plt.figure(FigureClass=SpotDiagramFigure, opt_model=opm, scale_type=Fit.User_Scale,
#                      user_scale_value=0.1, is_dark=isdark).plot()
abr_plt = plt.figure(FigureClass=RayFanFigure, opt_model=opm, data_type='Ray',
                     scale_type=Fit.All_Same).plot()
opm.seq_model.gaps[-1].thi
Out[6]:
128.01040075032165
Figure
In [7]:
#opm.seq_model.gaps[-1].thi -= 0.1
opm.update_model()
spot_plt = plt.figure(FigureClass=SpotDiagramFigure, opt_model=opm, is_dark=isdark).plot()

opm.seq_model.gaps[-1].thi
Out[7]:
128.01040075032165
Figure
In [8]:
layout_plt = plt.figure(FigureClass=InteractiveLayout, opt_model=opm, is_dark=isdark).plot()
Figure
In [9]:
def dump_pts(p, wi, ray_pkg, fld, wvl, foc):
    if ray_pkg is not None:
        ray = ray_pkg[0][-1]
        return ray_pkg[0][-1][0]
               
sm.trace_grid(dump_pts, 0, form='list', append_if_none=False)
Out[9]:
([array([[-2.72717191e-01, -1.21207640e-01,  0.00000000e+00],
         [-2.85739536e-01, -9.52465121e-02,  0.00000000e+00],
         [-2.95005176e-01, -6.55567057e-02,  0.00000000e+00],
         [-3.00550228e-01, -3.33944698e-02,  0.00000000e+00],
         [-3.02396198e-01, -4.66287808e-17,  0.00000000e+00],
         [-3.00550228e-01,  3.33944698e-02,  0.00000000e+00],
         [-2.95005176e-01,  6.55567057e-02,  0.00000000e+00],
         [-2.85739536e-01,  9.52465121e-02,  0.00000000e+00],
         [-2.72717191e-01,  1.21207640e-01,  0.00000000e+00],
         [-2.55640053e-01, -1.59775033e-01,  0.00000000e+00],
         [-2.70436427e-01, -1.35218214e-01,  0.00000000e+00],
         [-2.81885611e-01, -1.05707104e-01,  0.00000000e+00],
         [-2.90032188e-01, -7.25080471e-02,  0.00000000e+00],
         [-2.94907642e-01, -3.68634553e-02,  0.00000000e+00],
         [-2.96530718e-01, -5.14398798e-17,  0.00000000e+00],
         [-2.94907642e-01,  3.68634553e-02,  0.00000000e+00],
         [-2.90032188e-01,  7.25080471e-02,  0.00000000e+00],
         [-2.81885611e-01,  1.05707104e-01,  0.00000000e+00],
         [-2.70436427e-01,  1.35218214e-01,  0.00000000e+00],
         [-2.55640053e-01,  1.59775033e-01,  0.00000000e+00],
         [-2.10662688e-01, -2.10662688e-01,  0.00000000e+00],
         [-2.29448470e-01, -1.96670117e-01,  0.00000000e+00],
         [-2.45221517e-01, -1.75158226e-01,  0.00000000e+00],
         [-2.58044187e-01, -1.47453821e-01,  0.00000000e+00],
         [-2.67966537e-01, -1.14842801e-01,  0.00000000e+00],
         [-2.75026911e-01, -7.85791175e-02,  0.00000000e+00],
         [-2.79252390e-01, -3.98931985e-02,  0.00000000e+00],
         [-2.80659098e-01, -5.56418201e-17,  0.00000000e+00],
         [-2.79252390e-01,  3.98931985e-02,  0.00000000e+00],
         [-2.75026911e-01,  7.85791175e-02,  0.00000000e+00],
         [-2.67966537e-01,  1.14842801e-01,  0.00000000e+00],
         [-2.58044187e-01,  1.47453821e-01,  0.00000000e+00],
         [-2.45221517e-01,  1.75158226e-01,  0.00000000e+00],
         [-2.29448470e-01,  1.96670117e-01,  0.00000000e+00],
         [-2.10662688e-01,  2.10662688e-01,  0.00000000e+00],
         [-1.96670117e-01, -2.29448470e-01,  0.00000000e+00],
         [-2.12637761e-01, -2.12637761e-01,  0.00000000e+00],
         [-2.26045304e-01, -1.88371087e-01,  0.00000000e+00],
         [-2.36945398e-01, -1.57963599e-01,  0.00000000e+00],
         [-2.45380313e-01, -1.22690156e-01,  0.00000000e+00],
         [-2.51382432e-01, -8.37941438e-02,  0.00000000e+00],
         [-2.54974626e-01, -4.24957710e-02,  0.00000000e+00],
         [-2.56170517e-01, -5.92513345e-17,  0.00000000e+00],
         [-2.54974626e-01,  4.24957710e-02,  0.00000000e+00],
         [-2.51382432e-01,  8.37941438e-02,  0.00000000e+00],
         [-2.45380313e-01,  1.22690156e-01,  0.00000000e+00],
         [-2.36945398e-01,  1.57963599e-01,  0.00000000e+00],
         [-2.26045304e-01,  1.88371087e-01,  0.00000000e+00],
         [-2.12637761e-01,  2.12637761e-01,  0.00000000e+00],
         [-1.96670117e-01,  2.29448470e-01,  0.00000000e+00],
         [-1.78078709e-01,  2.37438278e-01,  0.00000000e+00],
         [-1.59775033e-01, -2.55640053e-01,  0.00000000e+00],
         [-1.75158226e-01, -2.45221517e-01,  0.00000000e+00],
         [-1.88371087e-01, -2.26045304e-01,  0.00000000e+00],
         [-1.99465993e-01, -1.99465993e-01,  0.00000000e+00],
         [-2.08486269e-01, -1.66789015e-01,  0.00000000e+00],
         [-2.15466701e-01, -1.29280021e-01,  0.00000000e+00],
         [-2.20433941e-01, -8.81735762e-02,  0.00000000e+00],
         [-2.23406813e-01, -4.46813626e-02,  0.00000000e+00],
         [-2.24396530e-01, -6.22825485e-17,  0.00000000e+00],
         [-2.23406813e-01,  4.46813626e-02,  0.00000000e+00],
         [-2.20433941e-01,  8.81735762e-02,  0.00000000e+00],
         [-2.15466701e-01,  1.29280021e-01,  0.00000000e+00],
         [-2.08486269e-01,  1.66789015e-01,  0.00000000e+00],
         [-1.99465993e-01,  1.99465993e-01,  0.00000000e+00],
         [-1.88371087e-01,  2.26045304e-01,  0.00000000e+00],
         [-1.75158226e-01,  2.45221517e-01,  0.00000000e+00],
         [-1.59775033e-01,  2.55640053e-01,  0.00000000e+00],
         [-1.21207640e-01, -2.72717191e-01,  0.00000000e+00],
         [-1.35218214e-01, -2.70436427e-01,  0.00000000e+00],
         [-1.47453821e-01, -2.58044187e-01,  0.00000000e+00],
         [-1.57963599e-01, -2.36945398e-01,  0.00000000e+00],
         [-1.66789015e-01, -2.08486269e-01,  0.00000000e+00],
         [-1.73964375e-01, -1.73964375e-01,  0.00000000e+00],
         [-1.79517224e-01, -1.34637918e-01,  0.00000000e+00],
         [-1.83468669e-01, -9.17343343e-02,  0.00000000e+00],
         [-1.85833618e-01, -4.64584045e-02,  0.00000000e+00],
         [-1.86620952e-01, -6.47471492e-17,  0.00000000e+00],
         [-1.85833618e-01,  4.64584045e-02,  0.00000000e+00],
         [-1.83468669e-01,  9.17343343e-02,  0.00000000e+00],
         [-1.79517224e-01,  1.34637918e-01,  0.00000000e+00],
         [-1.73964375e-01,  1.73964375e-01,  0.00000000e+00],
         [-1.66789015e-01,  2.08486269e-01,  0.00000000e+00],
         [-1.57963599e-01,  2.36945398e-01,  0.00000000e+00],
         [-1.47453821e-01,  2.58044187e-01,  0.00000000e+00],
         [-1.35218214e-01,  2.70436427e-01,  0.00000000e+00],
         [-1.21207640e-01,  2.72717191e-01,  0.00000000e+00],
         [-9.52465121e-02, -2.85739536e-01,  0.00000000e+00],
         [-1.05707104e-01, -2.81885611e-01,  0.00000000e+00],
         [-1.14842801e-01, -2.67966537e-01,  0.00000000e+00],
         [-1.22690156e-01, -2.45380313e-01,  0.00000000e+00],
         [-1.29280021e-01, -2.15466701e-01,  0.00000000e+00],
         [-1.34637918e-01, -1.79517224e-01,  0.00000000e+00],
         [-1.38784345e-01, -1.38784345e-01,  0.00000000e+00],
         [-1.41735008e-01, -9.44900056e-02,  0.00000000e+00],
         [-1.43501003e-01, -4.78336677e-02,  0.00000000e+00],
         [-1.44088936e-01, -6.66545228e-17,  0.00000000e+00],
         [-1.43501003e-01,  4.78336677e-02,  0.00000000e+00],
         [-1.41735008e-01,  9.44900056e-02,  0.00000000e+00],
         [-1.38784345e-01,  1.38784345e-01,  0.00000000e+00],
         [-1.34637918e-01,  1.79517224e-01,  0.00000000e+00],
         [-1.29280021e-01,  2.15466701e-01,  0.00000000e+00],
         [-1.22690156e-01,  2.45380313e-01,  0.00000000e+00],
         [-1.14842801e-01,  2.67966537e-01,  0.00000000e+00],
         [-1.05707104e-01,  2.81885611e-01,  0.00000000e+00],
         [-9.52465121e-02,  2.85739536e-01,  0.00000000e+00],
         [-6.55567057e-02, -2.95005176e-01,  0.00000000e+00],
         [-7.25080471e-02, -2.90032188e-01,  0.00000000e+00],
         [-7.85791175e-02, -2.75026911e-01,  0.00000000e+00],
         [-8.37941438e-02, -2.51382432e-01,  0.00000000e+00],
         [-8.81735762e-02, -2.20433941e-01,  0.00000000e+00],
         [-9.17343343e-02, -1.83468669e-01,  0.00000000e+00],
         [-9.44900056e-02, -1.41735008e-01,  0.00000000e+00],
         [-9.64510030e-02, -9.64510030e-02,  0.00000000e+00],
         [-9.76246820e-02, -4.88123410e-02,  0.00000000e+00],
         [-9.80154233e-02, -6.80118623e-17,  0.00000000e+00],
         [-9.76246820e-02,  4.88123410e-02,  0.00000000e+00],
         [-9.64510030e-02,  9.64510030e-02,  0.00000000e+00],
         [-9.44900056e-02,  1.41735008e-01,  0.00000000e+00],
         [-9.17343343e-02,  1.83468669e-01,  0.00000000e+00],
         [-8.81735762e-02,  2.20433941e-01,  0.00000000e+00],
         [-8.37941438e-02,  2.51382432e-01,  0.00000000e+00],
         [-7.85791175e-02,  2.75026911e-01,  0.00000000e+00],
         [-7.25080471e-02,  2.90032188e-01,  0.00000000e+00],
         [-6.55567057e-02,  2.95005176e-01,  0.00000000e+00],
         [-3.33944698e-02, -3.00550228e-01,  0.00000000e+00],
         [-3.68634553e-02, -2.94907642e-01,  0.00000000e+00],
         [-3.98931985e-02, -2.79252390e-01,  0.00000000e+00],
         [-4.24957710e-02, -2.54974626e-01,  0.00000000e+00],
         [-4.46813626e-02, -2.23406813e-01,  0.00000000e+00],
         [-4.64584045e-02, -1.85833618e-01,  0.00000000e+00],
         [-4.78336677e-02, -1.43501003e-01,  0.00000000e+00],
         [-4.88123410e-02, -9.76246820e-02,  0.00000000e+00],
         [-4.93980901e-02, -4.93980901e-02,  0.00000000e+00],
         [-4.95930981e-02, -6.88242492e-17,  0.00000000e+00],
         [-4.93980901e-02,  4.93980901e-02,  0.00000000e+00],
         [-4.88123410e-02,  9.76246820e-02,  0.00000000e+00],
         [-4.78336677e-02,  1.43501003e-01,  0.00000000e+00],
         [-4.64584045e-02,  1.85833618e-01,  0.00000000e+00],
         [-4.46813626e-02,  2.23406813e-01,  0.00000000e+00],
         [-4.24957710e-02,  2.54974626e-01,  0.00000000e+00],
         [-3.98931985e-02,  2.79252390e-01,  0.00000000e+00],
         [-3.68634553e-02,  2.94907642e-01,  0.00000000e+00],
         [-3.33944698e-02,  3.00550228e-01,  0.00000000e+00],
         [-4.66287808e-17, -3.02396198e-01,  0.00000000e+00],
         [-5.14398798e-17, -2.96530718e-01,  0.00000000e+00],
         [-5.56418201e-17, -2.80659098e-01,  0.00000000e+00],
         [-5.92513345e-17, -2.56170517e-01,  0.00000000e+00],
         [-6.22825485e-17, -2.24396530e-01,  0.00000000e+00],
         [-6.47471492e-17, -1.86620952e-01,  0.00000000e+00],
         [-6.66545228e-17, -1.44088936e-01,  0.00000000e+00],
         [-6.80118623e-17, -9.80154233e-02,  0.00000000e+00],
         [-6.88242492e-17, -4.95930981e-02,  0.00000000e+00],
         [-6.90947097e-17, -6.90947097e-17,  0.00000000e+00],
         [-6.88242492e-17,  4.95930981e-02,  0.00000000e+00],
         [-6.80118623e-17,  9.80154233e-02,  0.00000000e+00],
         [-6.66545228e-17,  1.44088936e-01,  0.00000000e+00],
         [-6.47471492e-17,  1.86620952e-01,  0.00000000e+00],
         [-6.22825485e-17,  2.24396530e-01,  0.00000000e+00],
         [-5.92513345e-17,  2.56170517e-01,  0.00000000e+00],
         [-5.56418201e-17,  2.80659098e-01,  0.00000000e+00],
         [-5.14398798e-17,  2.96530718e-01,  0.00000000e+00],
         [-4.66287808e-17,  3.02396198e-01,  0.00000000e+00],
         [-4.11889756e-17,  2.96797848e-01,  0.00000000e+00],
         [ 3.33944698e-02, -3.00550228e-01,  0.00000000e+00],
         [ 3.68634553e-02, -2.94907642e-01,  0.00000000e+00],
         [ 3.98931985e-02, -2.79252390e-01,  0.00000000e+00],
         [ 4.24957710e-02, -2.54974626e-01,  0.00000000e+00],
         [ 4.46813626e-02, -2.23406813e-01,  0.00000000e+00],
         [ 4.64584045e-02, -1.85833618e-01,  0.00000000e+00],
         [ 4.78336677e-02, -1.43501003e-01,  0.00000000e+00],
         [ 4.88123410e-02, -9.76246820e-02,  0.00000000e+00],
         [ 4.93980901e-02, -4.93980901e-02,  0.00000000e+00],
         [ 4.95930981e-02, -6.88242492e-17,  0.00000000e+00],
         [ 4.93980901e-02,  4.93980901e-02,  0.00000000e+00],
         [ 4.88123410e-02,  9.76246820e-02,  0.00000000e+00],
         [ 4.78336677e-02,  1.43501003e-01,  0.00000000e+00],
         [ 4.64584045e-02,  1.85833618e-01,  0.00000000e+00],
         [ 4.46813626e-02,  2.23406813e-01,  0.00000000e+00],
         [ 4.24957710e-02,  2.54974626e-01,  0.00000000e+00],
         [ 3.98931985e-02,  2.79252390e-01,  0.00000000e+00],
         [ 3.68634553e-02,  2.94907642e-01,  0.00000000e+00],
         [ 3.33944698e-02,  3.00550228e-01,  0.00000000e+00],
         [ 6.55567057e-02, -2.95005176e-01,  0.00000000e+00],
         [ 7.25080471e-02, -2.90032188e-01,  0.00000000e+00],
         [ 7.85791175e-02, -2.75026911e-01,  0.00000000e+00],
         [ 8.37941438e-02, -2.51382432e-01,  0.00000000e+00],
         [ 8.81735762e-02, -2.20433941e-01,  0.00000000e+00],
         [ 9.17343343e-02, -1.83468669e-01,  0.00000000e+00],
         [ 9.44900056e-02, -1.41735008e-01,  0.00000000e+00],
         [ 9.64510030e-02, -9.64510030e-02,  0.00000000e+00],
         [ 9.76246820e-02, -4.88123410e-02,  0.00000000e+00],
         [ 9.80154233e-02, -6.80118623e-17,  0.00000000e+00],
         [ 9.76246820e-02,  4.88123410e-02,  0.00000000e+00],
         [ 9.64510030e-02,  9.64510030e-02,  0.00000000e+00],
         [ 9.44900056e-02,  1.41735008e-01,  0.00000000e+00],
         [ 9.17343343e-02,  1.83468669e-01,  0.00000000e+00],
         [ 8.81735762e-02,  2.20433941e-01,  0.00000000e+00],
         [ 8.37941438e-02,  2.51382432e-01,  0.00000000e+00],
         [ 7.85791175e-02,  2.75026911e-01,  0.00000000e+00],
         [ 7.25080471e-02,  2.90032188e-01,  0.00000000e+00],
         [ 6.55567057e-02,  2.95005176e-01,  0.00000000e+00],
         [ 9.52465121e-02, -2.85739536e-01,  0.00000000e+00],
         [ 1.05707104e-01, -2.81885611e-01,  0.00000000e+00],
         [ 1.14842801e-01, -2.67966537e-01,  0.00000000e+00],
         [ 1.22690156e-01, -2.45380313e-01,  0.00000000e+00],
         [ 1.29280021e-01, -2.15466701e-01,  0.00000000e+00],
         [ 1.34637918e-01, -1.79517224e-01,  0.00000000e+00],
         [ 1.38784345e-01, -1.38784345e-01,  0.00000000e+00],
         [ 1.41735008e-01, -9.44900056e-02,  0.00000000e+00],
         [ 1.43501003e-01, -4.78336677e-02,  0.00000000e+00],
         [ 1.44088936e-01, -6.66545228e-17,  0.00000000e+00],
         [ 1.43501003e-01,  4.78336677e-02,  0.00000000e+00],
         [ 1.41735008e-01,  9.44900056e-02,  0.00000000e+00],
         [ 1.38784345e-01,  1.38784345e-01,  0.00000000e+00],
         [ 1.34637918e-01,  1.79517224e-01,  0.00000000e+00],
         [ 1.29280021e-01,  2.15466701e-01,  0.00000000e+00],
         [ 1.22690156e-01,  2.45380313e-01,  0.00000000e+00],
         [ 1.14842801e-01,  2.67966537e-01,  0.00000000e+00],
         [ 1.05707104e-01,  2.81885611e-01,  0.00000000e+00],
         [ 9.52465121e-02,  2.85739536e-01,  0.00000000e+00],
         [ 1.21207640e-01, -2.72717191e-01,  0.00000000e+00],
         [ 1.35218214e-01, -2.70436427e-01,  0.00000000e+00],
         [ 1.47453821e-01, -2.58044187e-01,  0.00000000e+00],
         [ 1.57963599e-01, -2.36945398e-01,  0.00000000e+00],
         [ 1.66789015e-01, -2.08486269e-01,  0.00000000e+00],
         [ 1.73964375e-01, -1.73964375e-01,  0.00000000e+00],
         [ 1.79517224e-01, -1.34637918e-01,  0.00000000e+00],
         [ 1.83468669e-01, -9.17343343e-02,  0.00000000e+00],
         [ 1.85833618e-01, -4.64584045e-02,  0.00000000e+00],
         [ 1.86620952e-01, -6.47471492e-17,  0.00000000e+00],
         [ 1.85833618e-01,  4.64584045e-02,  0.00000000e+00],
         [ 1.83468669e-01,  9.17343343e-02,  0.00000000e+00],
         [ 1.79517224e-01,  1.34637918e-01,  0.00000000e+00],
         [ 1.73964375e-01,  1.73964375e-01,  0.00000000e+00],
         [ 1.66789015e-01,  2.08486269e-01,  0.00000000e+00],
         [ 1.57963599e-01,  2.36945398e-01,  0.00000000e+00],
         [ 1.47453821e-01,  2.58044187e-01,  0.00000000e+00],
         [ 1.35218214e-01,  2.70436427e-01,  0.00000000e+00],
         [ 1.21207640e-01,  2.72717191e-01,  0.00000000e+00],
         [ 1.59775033e-01, -2.55640053e-01,  0.00000000e+00],
         [ 1.75158226e-01, -2.45221517e-01,  0.00000000e+00],
         [ 1.88371087e-01, -2.26045304e-01,  0.00000000e+00],
         [ 1.99465993e-01, -1.99465993e-01,  0.00000000e+00],
         [ 2.08486269e-01, -1.66789015e-01,  0.00000000e+00],
         [ 2.15466701e-01, -1.29280021e-01,  0.00000000e+00],
         [ 2.20433941e-01, -8.81735762e-02,  0.00000000e+00],
         [ 2.23406813e-01, -4.46813626e-02,  0.00000000e+00],
         [ 2.24396530e-01, -6.22825485e-17,  0.00000000e+00],
         [ 2.23406813e-01,  4.46813626e-02,  0.00000000e+00],
         [ 2.20433941e-01,  8.81735762e-02,  0.00000000e+00],
         [ 2.15466701e-01,  1.29280021e-01,  0.00000000e+00],
         [ 2.08486269e-01,  1.66789015e-01,  0.00000000e+00],
         [ 1.99465993e-01,  1.99465993e-01,  0.00000000e+00],
         [ 1.88371087e-01,  2.26045304e-01,  0.00000000e+00],
         [ 1.75158226e-01,  2.45221517e-01,  0.00000000e+00],
         [ 1.59775033e-01,  2.55640053e-01,  0.00000000e+00],
         [ 1.96670117e-01, -2.29448470e-01,  0.00000000e+00],
         [ 2.12637761e-01, -2.12637761e-01,  0.00000000e+00],
         [ 2.26045304e-01, -1.88371087e-01,  0.00000000e+00],
         [ 2.36945398e-01, -1.57963599e-01,  0.00000000e+00],
         [ 2.45380313e-01, -1.22690156e-01,  0.00000000e+00],
         [ 2.51382432e-01, -8.37941438e-02,  0.00000000e+00],
         [ 2.54974626e-01, -4.24957710e-02,  0.00000000e+00],
         [ 2.56170517e-01, -5.92513345e-17,  0.00000000e+00],
         [ 2.54974626e-01,  4.24957710e-02,  0.00000000e+00],
         [ 2.51382432e-01,  8.37941438e-02,  0.00000000e+00],
         [ 2.45380313e-01,  1.22690156e-01,  0.00000000e+00],
         [ 2.36945398e-01,  1.57963599e-01,  0.00000000e+00],
         [ 2.26045304e-01,  1.88371087e-01,  0.00000000e+00],
         [ 2.12637761e-01,  2.12637761e-01,  0.00000000e+00],
         [ 1.96670117e-01,  2.29448470e-01,  0.00000000e+00],
         [ 1.78078709e-01,  2.37438278e-01,  0.00000000e+00],
         [ 2.10662688e-01, -2.10662688e-01,  0.00000000e+00],
         [ 2.29448470e-01, -1.96670117e-01,  0.00000000e+00],
         [ 2.45221517e-01, -1.75158226e-01,  0.00000000e+00],
         [ 2.58044187e-01, -1.47453821e-01,  0.00000000e+00],
         [ 2.67966537e-01, -1.14842801e-01,  0.00000000e+00],
         [ 2.75026911e-01, -7.85791175e-02,  0.00000000e+00],
         [ 2.79252390e-01, -3.98931985e-02,  0.00000000e+00],
         [ 2.80659098e-01, -5.56418201e-17,  0.00000000e+00],
         [ 2.79252390e-01,  3.98931985e-02,  0.00000000e+00],
         [ 2.75026911e-01,  7.85791175e-02,  0.00000000e+00],
         [ 2.67966537e-01,  1.14842801e-01,  0.00000000e+00],
         [ 2.58044187e-01,  1.47453821e-01,  0.00000000e+00],
         [ 2.45221517e-01,  1.75158226e-01,  0.00000000e+00],
         [ 2.29448470e-01,  1.96670117e-01,  0.00000000e+00],
         [ 2.10662688e-01,  2.10662688e-01,  0.00000000e+00],
         [ 2.37438278e-01, -1.78078709e-01,  0.00000000e+00],
         [ 2.55640053e-01, -1.59775033e-01,  0.00000000e+00],
         [ 2.70436427e-01, -1.35218214e-01,  0.00000000e+00],
         [ 2.81885611e-01, -1.05707104e-01,  0.00000000e+00],
         [ 2.90032188e-01, -7.25080471e-02,  0.00000000e+00],
         [ 2.94907642e-01, -3.68634553e-02,  0.00000000e+00],
         [ 2.96530718e-01, -5.14398798e-17,  0.00000000e+00],
         [ 2.94907642e-01,  3.68634553e-02,  0.00000000e+00],
         [ 2.90032188e-01,  7.25080471e-02,  0.00000000e+00],
         [ 2.81885611e-01,  1.05707104e-01,  0.00000000e+00],
         [ 2.70436427e-01,  1.35218214e-01,  0.00000000e+00],
         [ 2.55640053e-01,  1.59775033e-01,  0.00000000e+00],
         [ 2.37438278e-01,  1.78078709e-01,  0.00000000e+00],
         [ 2.72717191e-01, -1.21207640e-01,  0.00000000e+00],
         [ 2.85739536e-01, -9.52465121e-02,  0.00000000e+00],
         [ 2.95005176e-01, -6.55567057e-02,  0.00000000e+00],
         [ 3.00550228e-01, -3.33944698e-02,  0.00000000e+00],
         [ 3.02396198e-01, -4.66287808e-17,  0.00000000e+00],
         [ 3.00550228e-01,  3.33944698e-02,  0.00000000e+00],
         [ 2.95005176e-01,  6.55567057e-02,  0.00000000e+00],
         [ 2.85739536e-01,  9.52465121e-02,  0.00000000e+00],
         [ 2.72717191e-01,  1.21207640e-01,  0.00000000e+00],
         [ 2.96797848e-01, -4.11889756e-17,  0.00000000e+00]]),
  array([[-2.79188043e-01, -1.24083575e-01,  0.00000000e+00],
         [-2.92195070e-01, -9.73983567e-02,  0.00000000e+00],
         [-3.01449655e-01, -6.69888123e-02,  0.00000000e+00],
         [-3.06988030e-01, -3.41097811e-02,  0.00000000e+00],
         [-3.08831768e-01, -4.76211304e-17,  0.00000000e+00],
         [-3.06988030e-01,  3.41097811e-02,  0.00000000e+00],
         [-3.01449655e-01,  6.69888123e-02,  0.00000000e+00],
         [-2.92195070e-01,  9.73983567e-02,  0.00000000e+00],
         [-2.79188043e-01,  1.24083575e-01,  0.00000000e+00],
         [-2.61376346e-01, -1.63360216e-01,  0.00000000e+00],
         [-2.76154945e-01, -1.38077473e-01,  0.00000000e+00],
         [-2.87590120e-01, -1.07846295e-01,  0.00000000e+00],
         [-2.95726595e-01, -7.39316488e-02,  0.00000000e+00],
         [-3.00595949e-01, -3.75744936e-02,  0.00000000e+00],
         [-3.02216986e-01, -5.24262900e-17,  0.00000000e+00],
         [-3.00595949e-01,  3.75744936e-02,  0.00000000e+00],
         [-2.95726595e-01,  7.39316488e-02,  0.00000000e+00],
         [-2.87590120e-01,  1.07846295e-01,  0.00000000e+00],
         [-2.76154945e-01,  1.38077473e-01,  0.00000000e+00],
         [-2.61376346e-01,  1.63360216e-01,  0.00000000e+00],
         [-2.15697263e-01, -2.15697263e-01,  0.00000000e+00],
         [-2.34460843e-01, -2.00966437e-01,  0.00000000e+00],
         [-2.50214722e-01, -1.78724801e-01,  0.00000000e+00],
         [-2.63021455e-01, -1.50297975e-01,  0.00000000e+00],
         [-2.72931256e-01, -1.16970538e-01,  0.00000000e+00],
         [-2.79982585e-01, -7.99950242e-02,  0.00000000e+00],
         [-2.84202604e-01, -4.06003719e-02,  0.00000000e+00],
         [-2.85607486e-01, -5.66228584e-17,  0.00000000e+00],
         [-2.84202604e-01,  4.06003719e-02,  0.00000000e+00],
         [-2.79982585e-01,  7.99950242e-02,  0.00000000e+00],
         [-2.72931256e-01,  1.16970538e-01,  0.00000000e+00],
         [-2.63021455e-01,  1.50297975e-01,  0.00000000e+00],
         [-2.50214722e-01,  1.78724801e-01,  0.00000000e+00],
         [-2.34460843e-01,  2.00966437e-01,  0.00000000e+00],
         [-2.15697263e-01,  2.15697263e-01,  0.00000000e+00],
         [-2.00966437e-01, -2.34460843e-01,  0.00000000e+00],
         [-2.16914632e-01, -2.16914632e-01,  0.00000000e+00],
         [-2.30305401e-01, -1.91921168e-01,  0.00000000e+00],
         [-2.41191560e-01, -1.60794374e-01,  0.00000000e+00],
         [-2.49615508e-01, -1.24807754e-01,  0.00000000e+00],
         [-2.55609725e-01, -8.52032418e-02,  0.00000000e+00],
         [-2.59197152e-01, -4.31995254e-02,  0.00000000e+00],
         [-2.60391449e-01, -6.02276213e-17,  0.00000000e+00],
         [-2.59197152e-01,  4.31995254e-02,  0.00000000e+00],
         [-2.55609725e-01,  8.52032418e-02,  0.00000000e+00],
         [-2.49615508e-01,  1.24807754e-01,  0.00000000e+00],
         [-2.41191560e-01,  1.60794374e-01,  0.00000000e+00],
         [-2.30305401e-01,  1.91921168e-01,  0.00000000e+00],
         [-2.16914632e-01,  2.16914632e-01,  0.00000000e+00],
         [-2.00966437e-01,  2.34460843e-01,  0.00000000e+00],
         [-1.82396949e-01,  2.43195931e-01,  0.00000000e+00],
         [-1.63360216e-01, -2.61376346e-01,  0.00000000e+00],
         [-1.78724801e-01, -2.50214722e-01,  0.00000000e+00],
         [-1.91921168e-01, -2.30305401e-01,  0.00000000e+00],
         [-2.03001860e-01, -2.03001860e-01,  0.00000000e+00],
         [-2.12010334e-01, -1.69608268e-01,  0.00000000e+00],
         [-2.18981483e-01, -1.31388890e-01,  0.00000000e+00],
         [-2.23942037e-01, -8.95768146e-02,  0.00000000e+00],
         [-2.26910876e-01, -4.53821751e-02,  0.00000000e+00],
         [-2.27899244e-01, -6.32547470e-17,  0.00000000e+00],
         [-2.26910876e-01,  4.53821751e-02,  0.00000000e+00],
         [-2.23942037e-01,  8.95768146e-02,  0.00000000e+00],
         [-2.18981483e-01,  1.31388890e-01,  0.00000000e+00],
         [-2.12010334e-01,  1.69608268e-01,  0.00000000e+00],
         [-2.03001860e-01,  2.03001860e-01,  0.00000000e+00],
         [-1.91921168e-01,  2.30305401e-01,  0.00000000e+00],
         [-1.78724801e-01,  2.50214722e-01,  0.00000000e+00],
         [-1.63360216e-01,  2.61376346e-01,  0.00000000e+00],
         [-1.24083575e-01, -2.79188043e-01,  0.00000000e+00],
         [-1.38077473e-01, -2.76154945e-01,  0.00000000e+00],
         [-1.50297975e-01, -2.63021455e-01,  0.00000000e+00],
         [-1.60794374e-01, -2.41191560e-01,  0.00000000e+00],
         [-1.69608268e-01, -2.12010334e-01,  0.00000000e+00],
         [-1.76774065e-01, -1.76774065e-01,  0.00000000e+00],
         [-1.82319395e-01, -1.36739546e-01,  0.00000000e+00],
         [-1.86265426e-01, -9.31327130e-02,  0.00000000e+00],
         [-1.88627110e-01, -4.71567776e-02,  0.00000000e+00],
         [-1.89413352e-01, -6.57159578e-17,  0.00000000e+00],
         [-1.88627110e-01,  4.71567776e-02,  0.00000000e+00],
         [-1.86265426e-01,  9.31327130e-02,  0.00000000e+00],
         [-1.82319395e-01,  1.36739546e-01,  0.00000000e+00],
         [-1.76774065e-01,  1.76774065e-01,  0.00000000e+00],
         [-1.69608268e-01,  2.12010334e-01,  0.00000000e+00],
         [-1.60794374e-01,  2.41191560e-01,  0.00000000e+00],
         [-1.50297975e-01,  2.63021455e-01,  0.00000000e+00],
         [-1.38077473e-01,  2.76154945e-01,  0.00000000e+00],
         [-1.24083575e-01,  2.79188043e-01,  0.00000000e+00],
         [-9.73983567e-02, -2.92195070e-01,  0.00000000e+00],
         [-1.07846295e-01, -2.87590120e-01,  0.00000000e+00],
         [-1.16970538e-01, -2.72931256e-01,  0.00000000e+00],
         [-1.24807754e-01, -2.49615508e-01,  0.00000000e+00],
         [-1.31388890e-01, -2.18981483e-01,  0.00000000e+00],
         [-1.36739546e-01, -1.82319395e-01,  0.00000000e+00],
         [-1.40880282e-01, -1.40880282e-01,  0.00000000e+00],
         [-1.43826847e-01, -9.58845649e-02,  0.00000000e+00],
         [-1.45590371e-01, -4.85301236e-02,  0.00000000e+00],
         [-1.46177478e-01, -6.76206675e-17,  0.00000000e+00],
         [-1.45590371e-01,  4.85301236e-02,  0.00000000e+00],
         [-1.43826847e-01,  9.58845649e-02,  0.00000000e+00],
         [-1.40880282e-01,  1.40880282e-01,  0.00000000e+00],
         [-1.36739546e-01,  1.82319395e-01,  0.00000000e+00],
         [-1.31388890e-01,  2.18981483e-01,  0.00000000e+00],
         [-1.24807754e-01,  2.49615508e-01,  0.00000000e+00],
         [-1.16970538e-01,  2.72931256e-01,  0.00000000e+00],
         [-1.07846295e-01,  2.87590120e-01,  0.00000000e+00],
         [-9.73983567e-02,  2.92195070e-01,  0.00000000e+00],
         [-6.69888123e-02, -3.01449655e-01,  0.00000000e+00],
         [-7.39316488e-02, -2.95726595e-01,  0.00000000e+00],
         [-7.99950242e-02, -2.79982585e-01,  0.00000000e+00],
         [-8.52032418e-02, -2.55609725e-01,  0.00000000e+00],
         [-8.95768146e-02, -2.23942037e-01,  0.00000000e+00],
         [-9.31327130e-02, -1.86265426e-01,  0.00000000e+00],
         [-9.58845649e-02, -1.43826847e-01,  0.00000000e+00],
         [-9.78428131e-02, -9.78428131e-02,  0.00000000e+00],
         [-9.90148341e-02, -4.95074171e-02,  0.00000000e+00],
         [-9.94050215e-02, -6.89760898e-17,  0.00000000e+00],
         [-9.90148341e-02,  4.95074171e-02,  0.00000000e+00],
         [-9.78428131e-02,  9.78428131e-02,  0.00000000e+00],
         [-9.58845649e-02,  1.43826847e-01,  0.00000000e+00],
         [-9.31327130e-02,  1.86265426e-01,  0.00000000e+00],
         [-8.95768146e-02,  2.23942037e-01,  0.00000000e+00],
         [-8.52032418e-02,  2.55609725e-01,  0.00000000e+00],
         [-7.99950242e-02,  2.79982585e-01,  0.00000000e+00],
         [-7.39316488e-02,  2.95726595e-01,  0.00000000e+00],
         [-6.69888123e-02,  3.01449655e-01,  0.00000000e+00],
         [-3.41097811e-02, -3.06988030e-01,  0.00000000e+00],
         [-3.75744936e-02, -3.00595949e-01,  0.00000000e+00],
         [-4.06003719e-02, -2.84202604e-01,  0.00000000e+00],
         [-4.31995254e-02, -2.59197152e-01,  0.00000000e+00],
         [-4.53821751e-02, -2.26910876e-01,  0.00000000e+00],
         [-4.71567776e-02, -1.88627110e-01,  0.00000000e+00],
         [-4.85301236e-02, -1.45590371e-01,  0.00000000e+00],
         [-4.95074171e-02, -9.90148341e-02,  0.00000000e+00],
         [-5.00923341e-02, -5.00923341e-02,  0.00000000e+00],
         [-5.02870641e-02, -6.97873205e-17,  0.00000000e+00],
         [-5.00923341e-02,  5.00923341e-02,  0.00000000e+00],
         [-4.95074171e-02,  9.90148341e-02,  0.00000000e+00],
         [-4.85301236e-02,  1.45590371e-01,  0.00000000e+00],
         [-4.71567776e-02,  1.88627110e-01,  0.00000000e+00],
         [-4.53821751e-02,  2.26910876e-01,  0.00000000e+00],
         [-4.31995254e-02,  2.59197152e-01,  0.00000000e+00],
         [-4.06003719e-02,  2.84202604e-01,  0.00000000e+00],
         [-3.75744936e-02,  3.00595949e-01,  0.00000000e+00],
         [-3.41097811e-02,  3.06988030e-01,  0.00000000e+00],
         [-4.76211304e-17, -3.08831768e-01,  0.00000000e+00],
         [-5.24262900e-17, -3.02216986e-01,  0.00000000e+00],
         [-5.66228584e-17, -2.85607486e-01,  0.00000000e+00],
         [-6.02276213e-17, -2.60391449e-01,  0.00000000e+00],
         [-6.32547470e-17, -2.27899244e-01,  0.00000000e+00],
         [-6.57159578e-17, -1.89413352e-01,  0.00000000e+00],
         [-6.76206675e-17, -1.46177478e-01,  0.00000000e+00],
         [-6.89760898e-17, -9.94050215e-02,  0.00000000e+00],
         [-6.97873205e-17, -5.02870641e-02,  0.00000000e+00],
         [-7.00573948e-17, -7.00573948e-17,  0.00000000e+00],
         [-6.97873205e-17,  5.02870641e-02,  0.00000000e+00],
         [-6.89760898e-17,  9.94050215e-02,  0.00000000e+00],
         [-6.76206675e-17,  1.46177478e-01,  0.00000000e+00],
         [-6.57159578e-17,  1.89413352e-01,  0.00000000e+00],
         [-6.32547470e-17,  2.27899244e-01,  0.00000000e+00],
         [-6.02276213e-17,  2.60391449e-01,  0.00000000e+00],
         [-5.66228584e-17,  2.85607486e-01,  0.00000000e+00],
         [-5.24262900e-17,  3.02216986e-01,  0.00000000e+00],
         [-4.76211304e-17,  3.08831768e-01,  0.00000000e+00],
         [-4.21877691e-17,  3.03994914e-01,  0.00000000e+00],
         [ 3.41097811e-02, -3.06988030e-01,  0.00000000e+00],
         [ 3.75744936e-02, -3.00595949e-01,  0.00000000e+00],
         [ 4.06003719e-02, -2.84202604e-01,  0.00000000e+00],
         [ 4.31995254e-02, -2.59197152e-01,  0.00000000e+00],
         [ 4.53821751e-02, -2.26910876e-01,  0.00000000e+00],
         [ 4.71567776e-02, -1.88627110e-01,  0.00000000e+00],
         [ 4.85301236e-02, -1.45590371e-01,  0.00000000e+00],
         [ 4.95074171e-02, -9.90148341e-02,  0.00000000e+00],
         [ 5.00923341e-02, -5.00923341e-02,  0.00000000e+00],
         [ 5.02870641e-02, -6.97873205e-17,  0.00000000e+00],
         [ 5.00923341e-02,  5.00923341e-02,  0.00000000e+00],
         [ 4.95074171e-02,  9.90148341e-02,  0.00000000e+00],
         [ 4.85301236e-02,  1.45590371e-01,  0.00000000e+00],
         [ 4.71567776e-02,  1.88627110e-01,  0.00000000e+00],
         [ 4.53821751e-02,  2.26910876e-01,  0.00000000e+00],
         [ 4.31995254e-02,  2.59197152e-01,  0.00000000e+00],
         [ 4.06003719e-02,  2.84202604e-01,  0.00000000e+00],
         [ 3.75744936e-02,  3.00595949e-01,  0.00000000e+00],
         [ 3.41097811e-02,  3.06988030e-01,  0.00000000e+00],
         [ 6.69888123e-02, -3.01449655e-01,  0.00000000e+00],
         [ 7.39316488e-02, -2.95726595e-01,  0.00000000e+00],
         [ 7.99950242e-02, -2.79982585e-01,  0.00000000e+00],
         [ 8.52032418e-02, -2.55609725e-01,  0.00000000e+00],
         [ 8.95768146e-02, -2.23942037e-01,  0.00000000e+00],
         [ 9.31327130e-02, -1.86265426e-01,  0.00000000e+00],
         [ 9.58845649e-02, -1.43826847e-01,  0.00000000e+00],
         [ 9.78428131e-02, -9.78428131e-02,  0.00000000e+00],
         [ 9.90148341e-02, -4.95074171e-02,  0.00000000e+00],
         [ 9.94050215e-02, -6.89760898e-17,  0.00000000e+00],
         [ 9.90148341e-02,  4.95074171e-02,  0.00000000e+00],
         [ 9.78428131e-02,  9.78428131e-02,  0.00000000e+00],
         [ 9.58845649e-02,  1.43826847e-01,  0.00000000e+00],
         [ 9.31327130e-02,  1.86265426e-01,  0.00000000e+00],
         [ 8.95768146e-02,  2.23942037e-01,  0.00000000e+00],
         [ 8.52032418e-02,  2.55609725e-01,  0.00000000e+00],
         [ 7.99950242e-02,  2.79982585e-01,  0.00000000e+00],
         [ 7.39316488e-02,  2.95726595e-01,  0.00000000e+00],
         [ 6.69888123e-02,  3.01449655e-01,  0.00000000e+00],
         [ 9.73983567e-02, -2.92195070e-01,  0.00000000e+00],
         [ 1.07846295e-01, -2.87590120e-01,  0.00000000e+00],
         [ 1.16970538e-01, -2.72931256e-01,  0.00000000e+00],
         [ 1.24807754e-01, -2.49615508e-01,  0.00000000e+00],
         [ 1.31388890e-01, -2.18981483e-01,  0.00000000e+00],
         [ 1.36739546e-01, -1.82319395e-01,  0.00000000e+00],
         [ 1.40880282e-01, -1.40880282e-01,  0.00000000e+00],
         [ 1.43826847e-01, -9.58845649e-02,  0.00000000e+00],
         [ 1.45590371e-01, -4.85301236e-02,  0.00000000e+00],
         [ 1.46177478e-01, -6.76206675e-17,  0.00000000e+00],
         [ 1.45590371e-01,  4.85301236e-02,  0.00000000e+00],
         [ 1.43826847e-01,  9.58845649e-02,  0.00000000e+00],
         [ 1.40880282e-01,  1.40880282e-01,  0.00000000e+00],
         [ 1.36739546e-01,  1.82319395e-01,  0.00000000e+00],
         [ 1.31388890e-01,  2.18981483e-01,  0.00000000e+00],
         [ 1.24807754e-01,  2.49615508e-01,  0.00000000e+00],
         [ 1.16970538e-01,  2.72931256e-01,  0.00000000e+00],
         [ 1.07846295e-01,  2.87590120e-01,  0.00000000e+00],
         [ 9.73983567e-02,  2.92195070e-01,  0.00000000e+00],
         [ 1.24083575e-01, -2.79188043e-01,  0.00000000e+00],
         [ 1.38077473e-01, -2.76154945e-01,  0.00000000e+00],
         [ 1.50297975e-01, -2.63021455e-01,  0.00000000e+00],
         [ 1.60794374e-01, -2.41191560e-01,  0.00000000e+00],
         [ 1.69608268e-01, -2.12010334e-01,  0.00000000e+00],
         [ 1.76774065e-01, -1.76774065e-01,  0.00000000e+00],
         [ 1.82319395e-01, -1.36739546e-01,  0.00000000e+00],
         [ 1.86265426e-01, -9.31327130e-02,  0.00000000e+00],
         [ 1.88627110e-01, -4.71567776e-02,  0.00000000e+00],
         [ 1.89413352e-01, -6.57159578e-17,  0.00000000e+00],
         [ 1.88627110e-01,  4.71567776e-02,  0.00000000e+00],
         [ 1.86265426e-01,  9.31327130e-02,  0.00000000e+00],
         [ 1.82319395e-01,  1.36739546e-01,  0.00000000e+00],
         [ 1.76774065e-01,  1.76774065e-01,  0.00000000e+00],
         [ 1.69608268e-01,  2.12010334e-01,  0.00000000e+00],
         [ 1.60794374e-01,  2.41191560e-01,  0.00000000e+00],
         [ 1.50297975e-01,  2.63021455e-01,  0.00000000e+00],
         [ 1.38077473e-01,  2.76154945e-01,  0.00000000e+00],
         [ 1.24083575e-01,  2.79188043e-01,  0.00000000e+00],
         [ 1.63360216e-01, -2.61376346e-01,  0.00000000e+00],
         [ 1.78724801e-01, -2.50214722e-01,  0.00000000e+00],
         [ 1.91921168e-01, -2.30305401e-01,  0.00000000e+00],
         [ 2.03001860e-01, -2.03001860e-01,  0.00000000e+00],
         [ 2.12010334e-01, -1.69608268e-01,  0.00000000e+00],
         [ 2.18981483e-01, -1.31388890e-01,  0.00000000e+00],
         [ 2.23942037e-01, -8.95768146e-02,  0.00000000e+00],
         [ 2.26910876e-01, -4.53821751e-02,  0.00000000e+00],
         [ 2.27899244e-01, -6.32547470e-17,  0.00000000e+00],
         [ 2.26910876e-01,  4.53821751e-02,  0.00000000e+00],
         [ 2.23942037e-01,  8.95768146e-02,  0.00000000e+00],
         [ 2.18981483e-01,  1.31388890e-01,  0.00000000e+00],
         [ 2.12010334e-01,  1.69608268e-01,  0.00000000e+00],
         [ 2.03001860e-01,  2.03001860e-01,  0.00000000e+00],
         [ 1.91921168e-01,  2.30305401e-01,  0.00000000e+00],
         [ 1.78724801e-01,  2.50214722e-01,  0.00000000e+00],
         [ 1.63360216e-01,  2.61376346e-01,  0.00000000e+00],
         [ 2.00966437e-01, -2.34460843e-01,  0.00000000e+00],
         [ 2.16914632e-01, -2.16914632e-01,  0.00000000e+00],
         [ 2.30305401e-01, -1.91921168e-01,  0.00000000e+00],
         [ 2.41191560e-01, -1.60794374e-01,  0.00000000e+00],
         [ 2.49615508e-01, -1.24807754e-01,  0.00000000e+00],
         [ 2.55609725e-01, -8.52032418e-02,  0.00000000e+00],
         [ 2.59197152e-01, -4.31995254e-02,  0.00000000e+00],
         [ 2.60391449e-01, -6.02276213e-17,  0.00000000e+00],
         [ 2.59197152e-01,  4.31995254e-02,  0.00000000e+00],
         [ 2.55609725e-01,  8.52032418e-02,  0.00000000e+00],
         [ 2.49615508e-01,  1.24807754e-01,  0.00000000e+00],
         [ 2.41191560e-01,  1.60794374e-01,  0.00000000e+00],
         [ 2.30305401e-01,  1.91921168e-01,  0.00000000e+00],
         [ 2.16914632e-01,  2.16914632e-01,  0.00000000e+00],
         [ 2.00966437e-01,  2.34460843e-01,  0.00000000e+00],
         [ 1.82396949e-01,  2.43195931e-01,  0.00000000e+00],
         [ 2.15697263e-01, -2.15697263e-01,  0.00000000e+00],
         [ 2.34460843e-01, -2.00966437e-01,  0.00000000e+00],
         [ 2.50214722e-01, -1.78724801e-01,  0.00000000e+00],
         [ 2.63021455e-01, -1.50297975e-01,  0.00000000e+00],
         [ 2.72931256e-01, -1.16970538e-01,  0.00000000e+00],
         [ 2.79982585e-01, -7.99950242e-02,  0.00000000e+00],
         [ 2.84202604e-01, -4.06003719e-02,  0.00000000e+00],
         [ 2.85607486e-01, -5.66228584e-17,  0.00000000e+00],
         [ 2.84202604e-01,  4.06003719e-02,  0.00000000e+00],
         [ 2.79982585e-01,  7.99950242e-02,  0.00000000e+00],
         [ 2.72931256e-01,  1.16970538e-01,  0.00000000e+00],
         [ 2.63021455e-01,  1.50297975e-01,  0.00000000e+00],
         [ 2.50214722e-01,  1.78724801e-01,  0.00000000e+00],
         [ 2.34460843e-01,  2.00966437e-01,  0.00000000e+00],
         [ 2.15697263e-01,  2.15697263e-01,  0.00000000e+00],
         [ 2.43195931e-01, -1.82396949e-01,  0.00000000e+00],
         [ 2.61376346e-01, -1.63360216e-01,  0.00000000e+00],
         [ 2.76154945e-01, -1.38077473e-01,  0.00000000e+00],
         [ 2.87590120e-01, -1.07846295e-01,  0.00000000e+00],
         [ 2.95726595e-01, -7.39316488e-02,  0.00000000e+00],
         [ 3.00595949e-01, -3.75744936e-02,  0.00000000e+00],
         [ 3.02216986e-01, -5.24262900e-17,  0.00000000e+00],
         [ 3.00595949e-01,  3.75744936e-02,  0.00000000e+00],
         [ 2.95726595e-01,  7.39316488e-02,  0.00000000e+00],
         [ 2.87590120e-01,  1.07846295e-01,  0.00000000e+00],
         [ 2.76154945e-01,  1.38077473e-01,  0.00000000e+00],
         [ 2.61376346e-01,  1.63360216e-01,  0.00000000e+00],
         [ 2.43195931e-01,  1.82396949e-01,  0.00000000e+00],
         [ 2.79188043e-01, -1.24083575e-01,  0.00000000e+00],
         [ 2.92195070e-01, -9.73983567e-02,  0.00000000e+00],
         [ 3.01449655e-01, -6.69888123e-02,  0.00000000e+00],
         [ 3.06988030e-01, -3.41097811e-02,  0.00000000e+00],
         [ 3.08831768e-01, -4.76211304e-17,  0.00000000e+00],
         [ 3.06988030e-01,  3.41097811e-02,  0.00000000e+00],
         [ 3.01449655e-01,  6.69888123e-02,  0.00000000e+00],
         [ 2.92195070e-01,  9.73983567e-02,  0.00000000e+00],
         [ 2.79188043e-01,  1.24083575e-01,  0.00000000e+00],
         [ 3.03994914e-01, -4.21877691e-17,  0.00000000e+00]]),
  array([[-2.86525085e-01, -1.27344482e-01,  0.00000000e+00],
         [-2.99520519e-01, -9.98401731e-02,  0.00000000e+00],
         [-3.08766793e-01, -6.86148429e-02,  0.00000000e+00],
         [-3.14300169e-01, -3.49222410e-02,  0.00000000e+00],
         [-3.16142238e-01, -4.87483877e-17,  0.00000000e+00],
         [-3.14300169e-01,  3.49222410e-02,  0.00000000e+00],
         [-3.08766793e-01,  6.86148429e-02,  0.00000000e+00],
         [-2.99520519e-01,  9.98401731e-02,  0.00000000e+00],
         [-2.86525085e-01,  1.27344482e-01,  0.00000000e+00],
         [-2.67886381e-01, -1.67428988e-01,  0.00000000e+00],
         [-2.82651657e-01, -1.41325829e-01,  0.00000000e+00],
         [-2.94076421e-01, -1.10278658e-01,  0.00000000e+00],
         [-3.02205434e-01, -7.55513585e-02,  0.00000000e+00],
         [-3.07070300e-01, -3.83837875e-02,  0.00000000e+00],
         [-3.08689839e-01, -5.35491511e-17,  0.00000000e+00],
         [-3.07070300e-01,  3.83837875e-02,  0.00000000e+00],
         [-3.02205434e-01,  7.55513585e-02,  0.00000000e+00],
         [-2.94076421e-01,  1.10278658e-01,  0.00000000e+00],
         [-2.82651657e-01,  1.41325829e-01,  0.00000000e+00],
         [-2.67886381e-01,  1.67428988e-01,  0.00000000e+00],
         [-2.21405136e-01, -2.21405136e-01,  0.00000000e+00],
         [-2.40151950e-01, -2.05844529e-01,  0.00000000e+00],
         [-2.55891539e-01, -1.82779671e-01,  0.00000000e+00],
         [-2.68686513e-01, -1.53535150e-01,  0.00000000e+00],
         [-2.78587125e-01, -1.19394482e-01,  0.00000000e+00],
         [-2.85631869e-01, -8.16091055e-02,  0.00000000e+00],
         [-2.89847929e-01, -4.14068470e-02,  0.00000000e+00],
         [-2.91251490e-01, -5.77418055e-17,  0.00000000e+00],
         [-2.89847929e-01,  4.14068470e-02,  0.00000000e+00],
         [-2.85631869e-01,  8.16091055e-02,  0.00000000e+00],
         [-2.78587125e-01,  1.19394482e-01,  0.00000000e+00],
         [-2.68686513e-01,  1.53535150e-01,  0.00000000e+00],
         [-2.55891539e-01,  1.82779671e-01,  0.00000000e+00],
         [-2.40151950e-01,  2.05844529e-01,  0.00000000e+00],
         [-2.21405136e-01,  2.21405136e-01,  0.00000000e+00],
         [-2.05844529e-01, -2.40151950e-01,  0.00000000e+00],
         [-2.21778240e-01, -2.21778240e-01,  0.00000000e+00],
         [-2.35156668e-01, -1.95963890e-01,  0.00000000e+00],
         [-2.46032674e-01, -1.64021783e-01,  0.00000000e+00],
         [-2.54448690e-01, -1.27224345e-01,  0.00000000e+00],
         [-2.60437224e-01, -8.68124081e-02,  0.00000000e+00],
         [-2.64021234e-01, -4.40035390e-02,  0.00000000e+00],
         [-2.65214391e-01, -6.13431507e-17,  0.00000000e+00],
         [-2.64021234e-01,  4.40035390e-02,  0.00000000e+00],
         [-2.60437224e-01,  8.68124081e-02,  0.00000000e+00],
         [-2.54448690e-01,  1.27224345e-01,  0.00000000e+00],
         [-2.46032674e-01,  1.64021783e-01,  0.00000000e+00],
         [-2.35156668e-01,  1.95963890e-01,  0.00000000e+00],
         [-2.21778240e-01,  2.21778240e-01,  0.00000000e+00],
         [-2.05844529e-01,  2.40151950e-01,  0.00000000e+00],
         [-1.87291612e-01,  2.49722150e-01,  0.00000000e+00],
         [-1.67428988e-01, -2.67886381e-01,  0.00000000e+00],
         [-1.82779671e-01, -2.55891539e-01,  0.00000000e+00],
         [-1.95963890e-01, -2.35156668e-01,  0.00000000e+00],
         [-2.07034235e-01, -2.07034235e-01,  0.00000000e+00],
         [-2.16034199e-01, -1.72827359e-01,  0.00000000e+00],
         [-2.22998700e-01, -1.33799220e-01,  0.00000000e+00],
         [-2.27954490e-01, -9.11817962e-02,  0.00000000e+00],
         [-2.30920466e-01, -4.61840933e-02,  0.00000000e+00],
         [-2.31907879e-01, -6.43673668e-17,  0.00000000e+00],
         [-2.30920466e-01,  4.61840933e-02,  0.00000000e+00],
         [-2.27954490e-01,  9.11817962e-02,  0.00000000e+00],
         [-2.22998700e-01,  1.33799220e-01,  0.00000000e+00],
         [-2.16034199e-01,  1.72827359e-01,  0.00000000e+00],
         [-2.07034235e-01,  2.07034235e-01,  0.00000000e+00],
         [-1.95963890e-01,  2.35156668e-01,  0.00000000e+00],
         [-1.82779671e-01,  2.55891539e-01,  0.00000000e+00],
         [-1.67428988e-01,  2.67886381e-01,  0.00000000e+00],
         [-1.27344482e-01, -2.86525085e-01,  0.00000000e+00],
         [-1.41325829e-01, -2.82651657e-01,  0.00000000e+00],
         [-1.53535150e-01, -2.68686513e-01,  0.00000000e+00],
         [-1.64021783e-01, -2.46032674e-01,  0.00000000e+00],
         [-1.72827359e-01, -2.16034199e-01,  0.00000000e+00],
         [-1.79986316e-01, -1.79986316e-01,  0.00000000e+00],
         [-1.85526303e-01, -1.39144728e-01,  0.00000000e+00],
         [-1.89468507e-01, -9.47342535e-02,  0.00000000e+00],
         [-1.91827890e-01, -4.79569726e-02,  0.00000000e+00],
         [-1.92613365e-01, -6.68261852e-17,  0.00000000e+00],
         [-1.91827890e-01,  4.79569726e-02,  0.00000000e+00],
         [-1.89468507e-01,  9.47342535e-02,  0.00000000e+00],
         [-1.85526303e-01,  1.39144728e-01,  0.00000000e+00],
         [-1.79986316e-01,  1.79986316e-01,  0.00000000e+00],
         [-1.72827359e-01,  2.16034199e-01,  0.00000000e+00],
         [-1.64021783e-01,  2.46032674e-01,  0.00000000e+00],
         [-1.53535150e-01,  2.68686513e-01,  0.00000000e+00],
         [-1.41325829e-01,  2.82651657e-01,  0.00000000e+00],
         [-1.27344482e-01,  2.86525085e-01,  0.00000000e+00],
         [-9.98401731e-02, -2.99520519e-01,  0.00000000e+00],
         [-1.10278658e-01, -2.94076421e-01,  0.00000000e+00],
         [-1.19394482e-01, -2.78587125e-01,  0.00000000e+00],
         [-1.27224345e-01, -2.54448690e-01,  0.00000000e+00],
         [-1.33799220e-01, -2.22998700e-01,  0.00000000e+00],
         [-1.39144728e-01, -1.85526303e-01,  0.00000000e+00],
         [-1.43281442e-01, -1.43281442e-01,  0.00000000e+00],
         [-1.46225128e-01, -9.74834186e-02,  0.00000000e+00],
         [-1.47986920e-01, -4.93289733e-02,  0.00000000e+00],
         [-1.48573450e-01, -6.87290270e-17,  0.00000000e+00],
         [-1.47986920e-01,  4.93289733e-02,  0.00000000e+00],
         [-1.46225128e-01,  9.74834186e-02,  0.00000000e+00],
         [-1.43281442e-01,  1.43281442e-01,  0.00000000e+00],
         [-1.39144728e-01,  1.85526303e-01,  0.00000000e+00],
         [-1.33799220e-01,  2.22998700e-01,  0.00000000e+00],
         [-1.27224345e-01,  2.54448690e-01,  0.00000000e+00],
         [-1.19394482e-01,  2.78587125e-01,  0.00000000e+00],
         [-1.10278658e-01,  2.94076421e-01,  0.00000000e+00],
         [-9.98401731e-02,  2.99520519e-01,  0.00000000e+00],
         [-6.86148429e-02, -3.08766793e-01,  0.00000000e+00],
         [-7.55513585e-02, -3.02205434e-01,  0.00000000e+00],
         [-8.16091055e-02, -2.85631869e-01,  0.00000000e+00],
         [-8.68124081e-02, -2.60437224e-01,  0.00000000e+00],
         [-9.11817962e-02, -2.27954490e-01,  0.00000000e+00],
         [-9.47342535e-02, -1.89468507e-01,  0.00000000e+00],
         [-9.74834186e-02, -1.46225128e-01,  0.00000000e+00],
         [-9.94397422e-02, -9.94397422e-02,  0.00000000e+00],
         [-1.00610606e-01, -5.03053032e-02,  0.00000000e+00],
         [-1.01000408e-01, -7.00831113e-17,  0.00000000e+00],
         [-1.00610606e-01,  5.03053032e-02,  0.00000000e+00],
         [-9.94397422e-02,  9.94397422e-02,  0.00000000e+00],
         [-9.74834186e-02,  1.46225128e-01,  0.00000000e+00],
         [-9.47342535e-02,  1.89468507e-01,  0.00000000e+00],
         [-9.11817962e-02,  2.27954490e-01,  0.00000000e+00],
         [-8.68124081e-02,  2.60437224e-01,  0.00000000e+00],
         [-8.16091055e-02,  2.85631869e-01,  0.00000000e+00],
         [-7.55513585e-02,  3.02205434e-01,  0.00000000e+00],
         [-6.86148429e-02,  3.08766793e-01,  0.00000000e+00],
         [-3.49222410e-02, -3.14300169e-01,  0.00000000e+00],
         [-3.83837875e-02, -3.07070300e-01,  0.00000000e+00],
         [-4.14068470e-02, -2.89847929e-01,  0.00000000e+00],
         [-4.40035390e-02, -2.64021234e-01,  0.00000000e+00],
         [-4.61840933e-02, -2.30920466e-01,  0.00000000e+00],
         [-4.79569726e-02, -1.91827890e-01,  0.00000000e+00],
         [-4.93289733e-02, -1.47986920e-01,  0.00000000e+00],
         [-5.03053032e-02, -1.00610606e-01,  0.00000000e+00],
         [-5.08896410e-02, -5.08896410e-02,  0.00000000e+00],
         [-5.10841777e-02, -7.08935378e-17,  0.00000000e+00],
         [-5.08896410e-02,  5.08896410e-02,  0.00000000e+00],
         [-5.03053032e-02,  1.00610606e-01,  0.00000000e+00],
         [-4.93289733e-02,  1.47986920e-01,  0.00000000e+00],
         [-4.79569726e-02,  1.91827890e-01,  0.00000000e+00],
         [-4.61840933e-02,  2.30920466e-01,  0.00000000e+00],
         [-4.40035390e-02,  2.64021234e-01,  0.00000000e+00],
         [-4.14068470e-02,  2.89847929e-01,  0.00000000e+00],
         [-3.83837875e-02,  3.07070300e-01,  0.00000000e+00],
         [-3.49222410e-02,  3.14300169e-01,  0.00000000e+00],
         [-4.87483877e-17, -3.16142238e-01,  0.00000000e+00],
         [-5.35491511e-17, -3.08689839e-01,  0.00000000e+00],
         [-5.77418055e-17, -2.91251490e-01,  0.00000000e+00],
         [-6.13431507e-17, -2.65214391e-01,  0.00000000e+00],
         [-6.43673668e-17, -2.31907879e-01,  0.00000000e+00],
         [-6.68261852e-17, -1.92613365e-01,  0.00000000e+00],
         [-6.87290270e-17, -1.48573450e-01,  0.00000000e+00],
         [-7.00831113e-17, -1.01000408e-01,  0.00000000e+00],
         [-7.08935378e-17, -5.10841777e-02,  0.00000000e+00],
         [-7.11633438e-17, -7.11633438e-17,  0.00000000e+00],
         [-7.08935378e-17,  5.10841777e-02,  0.00000000e+00],
         [-7.00831113e-17,  1.01000408e-01,  0.00000000e+00],
         [-6.87290270e-17,  1.48573450e-01,  0.00000000e+00],
         [-6.68261852e-17,  1.92613365e-01,  0.00000000e+00],
         [-6.43673668e-17,  2.31907879e-01,  0.00000000e+00],
         [-6.13431507e-17,  2.65214391e-01,  0.00000000e+00],
         [-5.77418055e-17,  2.91251490e-01,  0.00000000e+00],
         [-5.35491511e-17,  3.08689839e-01,  0.00000000e+00],
         [-4.87483877e-17,  3.16142238e-01,  0.00000000e+00],
         [-4.33198876e-17,  3.12152687e-01,  0.00000000e+00],
         [ 3.49222410e-02, -3.14300169e-01,  0.00000000e+00],
         [ 3.83837875e-02, -3.07070300e-01,  0.00000000e+00],
         [ 4.14068470e-02, -2.89847929e-01,  0.00000000e+00],
         [ 4.40035390e-02, -2.64021234e-01,  0.00000000e+00],
         [ 4.61840933e-02, -2.30920466e-01,  0.00000000e+00],
         [ 4.79569726e-02, -1.91827890e-01,  0.00000000e+00],
         [ 4.93289733e-02, -1.47986920e-01,  0.00000000e+00],
         [ 5.03053032e-02, -1.00610606e-01,  0.00000000e+00],
         [ 5.08896410e-02, -5.08896410e-02,  0.00000000e+00],
         [ 5.10841777e-02, -7.08935378e-17,  0.00000000e+00],
         [ 5.08896410e-02,  5.08896410e-02,  0.00000000e+00],
         [ 5.03053032e-02,  1.00610606e-01,  0.00000000e+00],
         [ 4.93289733e-02,  1.47986920e-01,  0.00000000e+00],
         [ 4.79569726e-02,  1.91827890e-01,  0.00000000e+00],
         [ 4.61840933e-02,  2.30920466e-01,  0.00000000e+00],
         [ 4.40035390e-02,  2.64021234e-01,  0.00000000e+00],
         [ 4.14068470e-02,  2.89847929e-01,  0.00000000e+00],
         [ 3.83837875e-02,  3.07070300e-01,  0.00000000e+00],
         [ 3.49222410e-02,  3.14300169e-01,  0.00000000e+00],
         [ 6.86148429e-02, -3.08766793e-01,  0.00000000e+00],
         [ 7.55513585e-02, -3.02205434e-01,  0.00000000e+00],
         [ 8.16091055e-02, -2.85631869e-01,  0.00000000e+00],
         [ 8.68124081e-02, -2.60437224e-01,  0.00000000e+00],
         [ 9.11817962e-02, -2.27954490e-01,  0.00000000e+00],
         [ 9.47342535e-02, -1.89468507e-01,  0.00000000e+00],
         [ 9.74834186e-02, -1.46225128e-01,  0.00000000e+00],
         [ 9.94397422e-02, -9.94397422e-02,  0.00000000e+00],
         [ 1.00610606e-01, -5.03053032e-02,  0.00000000e+00],
         [ 1.01000408e-01, -7.00831113e-17,  0.00000000e+00],
         [ 1.00610606e-01,  5.03053032e-02,  0.00000000e+00],
         [ 9.94397422e-02,  9.94397422e-02,  0.00000000e+00],
         [ 9.74834186e-02,  1.46225128e-01,  0.00000000e+00],
         [ 9.47342535e-02,  1.89468507e-01,  0.00000000e+00],
         [ 9.11817962e-02,  2.27954490e-01,  0.00000000e+00],
         [ 8.68124081e-02,  2.60437224e-01,  0.00000000e+00],
         [ 8.16091055e-02,  2.85631869e-01,  0.00000000e+00],
         [ 7.55513585e-02,  3.02205434e-01,  0.00000000e+00],
         [ 6.86148429e-02,  3.08766793e-01,  0.00000000e+00],
         [ 9.98401731e-02, -2.99520519e-01,  0.00000000e+00],
         [ 1.10278658e-01, -2.94076421e-01,  0.00000000e+00],
         [ 1.19394482e-01, -2.78587125e-01,  0.00000000e+00],
         [ 1.27224345e-01, -2.54448690e-01,  0.00000000e+00],
         [ 1.33799220e-01, -2.22998700e-01,  0.00000000e+00],
         [ 1.39144728e-01, -1.85526303e-01,  0.00000000e+00],
         [ 1.43281442e-01, -1.43281442e-01,  0.00000000e+00],
         [ 1.46225128e-01, -9.74834186e-02,  0.00000000e+00],
         [ 1.47986920e-01, -4.93289733e-02,  0.00000000e+00],
         [ 1.48573450e-01, -6.87290270e-17,  0.00000000e+00],
         [ 1.47986920e-01,  4.93289733e-02,  0.00000000e+00],
         [ 1.46225128e-01,  9.74834186e-02,  0.00000000e+00],
         [ 1.43281442e-01,  1.43281442e-01,  0.00000000e+00],
         [ 1.39144728e-01,  1.85526303e-01,  0.00000000e+00],
         [ 1.33799220e-01,  2.22998700e-01,  0.00000000e+00],
         [ 1.27224345e-01,  2.54448690e-01,  0.00000000e+00],
         [ 1.19394482e-01,  2.78587125e-01,  0.00000000e+00],
         [ 1.10278658e-01,  2.94076421e-01,  0.00000000e+00],
         [ 9.98401731e-02,  2.99520519e-01,  0.00000000e+00],
         [ 1.27344482e-01, -2.86525085e-01,  0.00000000e+00],
         [ 1.41325829e-01, -2.82651657e-01,  0.00000000e+00],
         [ 1.53535150e-01, -2.68686513e-01,  0.00000000e+00],
         [ 1.64021783e-01, -2.46032674e-01,  0.00000000e+00],
         [ 1.72827359e-01, -2.16034199e-01,  0.00000000e+00],
         [ 1.79986316e-01, -1.79986316e-01,  0.00000000e+00],
         [ 1.85526303e-01, -1.39144728e-01,  0.00000000e+00],
         [ 1.89468507e-01, -9.47342535e-02,  0.00000000e+00],
         [ 1.91827890e-01, -4.79569726e-02,  0.00000000e+00],
         [ 1.92613365e-01, -6.68261852e-17,  0.00000000e+00],
         [ 1.91827890e-01,  4.79569726e-02,  0.00000000e+00],
         [ 1.89468507e-01,  9.47342535e-02,  0.00000000e+00],
         [ 1.85526303e-01,  1.39144728e-01,  0.00000000e+00],
         [ 1.79986316e-01,  1.79986316e-01,  0.00000000e+00],
         [ 1.72827359e-01,  2.16034199e-01,  0.00000000e+00],
         [ 1.64021783e-01,  2.46032674e-01,  0.00000000e+00],
         [ 1.53535150e-01,  2.68686513e-01,  0.00000000e+00],
         [ 1.41325829e-01,  2.82651657e-01,  0.00000000e+00],
         [ 1.27344482e-01,  2.86525085e-01,  0.00000000e+00],
         [ 1.67428988e-01, -2.67886381e-01,  0.00000000e+00],
         [ 1.82779671e-01, -2.55891539e-01,  0.00000000e+00],
         [ 1.95963890e-01, -2.35156668e-01,  0.00000000e+00],
         [ 2.07034235e-01, -2.07034235e-01,  0.00000000e+00],
         [ 2.16034199e-01, -1.72827359e-01,  0.00000000e+00],
         [ 2.22998700e-01, -1.33799220e-01,  0.00000000e+00],
         [ 2.27954490e-01, -9.11817962e-02,  0.00000000e+00],
         [ 2.30920466e-01, -4.61840933e-02,  0.00000000e+00],
         [ 2.31907879e-01, -6.43673668e-17,  0.00000000e+00],
         [ 2.30920466e-01,  4.61840933e-02,  0.00000000e+00],
         [ 2.27954490e-01,  9.11817962e-02,  0.00000000e+00],
         [ 2.22998700e-01,  1.33799220e-01,  0.00000000e+00],
         [ 2.16034199e-01,  1.72827359e-01,  0.00000000e+00],
         [ 2.07034235e-01,  2.07034235e-01,  0.00000000e+00],
         [ 1.95963890e-01,  2.35156668e-01,  0.00000000e+00],
         [ 1.82779671e-01,  2.55891539e-01,  0.00000000e+00],
         [ 1.67428988e-01,  2.67886381e-01,  0.00000000e+00],
         [ 2.05844529e-01, -2.40151950e-01,  0.00000000e+00],
         [ 2.21778240e-01, -2.21778240e-01,  0.00000000e+00],
         [ 2.35156668e-01, -1.95963890e-01,  0.00000000e+00],
         [ 2.46032674e-01, -1.64021783e-01,  0.00000000e+00],
         [ 2.54448690e-01, -1.27224345e-01,  0.00000000e+00],
         [ 2.60437224e-01, -8.68124081e-02,  0.00000000e+00],
         [ 2.64021234e-01, -4.40035390e-02,  0.00000000e+00],
         [ 2.65214391e-01, -6.13431507e-17,  0.00000000e+00],
         [ 2.64021234e-01,  4.40035390e-02,  0.00000000e+00],
         [ 2.60437224e-01,  8.68124081e-02,  0.00000000e+00],
         [ 2.54448690e-01,  1.27224345e-01,  0.00000000e+00],
         [ 2.46032674e-01,  1.64021783e-01,  0.00000000e+00],
         [ 2.35156668e-01,  1.95963890e-01,  0.00000000e+00],
         [ 2.21778240e-01,  2.21778240e-01,  0.00000000e+00],
         [ 2.05844529e-01,  2.40151950e-01,  0.00000000e+00],
         [ 1.87291612e-01,  2.49722150e-01,  0.00000000e+00],
         [ 2.21405136e-01, -2.21405136e-01,  0.00000000e+00],
         [ 2.40151950e-01, -2.05844529e-01,  0.00000000e+00],
         [ 2.55891539e-01, -1.82779671e-01,  0.00000000e+00],
         [ 2.68686513e-01, -1.53535150e-01,  0.00000000e+00],
         [ 2.78587125e-01, -1.19394482e-01,  0.00000000e+00],
         [ 2.85631869e-01, -8.16091055e-02,  0.00000000e+00],
         [ 2.89847929e-01, -4.14068470e-02,  0.00000000e+00],
         [ 2.91251490e-01, -5.77418055e-17,  0.00000000e+00],
         [ 2.89847929e-01,  4.14068470e-02,  0.00000000e+00],
         [ 2.85631869e-01,  8.16091055e-02,  0.00000000e+00],
         [ 2.78587125e-01,  1.19394482e-01,  0.00000000e+00],
         [ 2.68686513e-01,  1.53535150e-01,  0.00000000e+00],
         [ 2.55891539e-01,  1.82779671e-01,  0.00000000e+00],
         [ 2.40151950e-01,  2.05844529e-01,  0.00000000e+00],
         [ 2.21405136e-01,  2.21405136e-01,  0.00000000e+00],
         [ 2.49722150e-01, -1.87291612e-01,  0.00000000e+00],
         [ 2.67886381e-01, -1.67428988e-01,  0.00000000e+00],
         [ 2.82651657e-01, -1.41325829e-01,  0.00000000e+00],
         [ 2.94076421e-01, -1.10278658e-01,  0.00000000e+00],
         [ 3.02205434e-01, -7.55513585e-02,  0.00000000e+00],
         [ 3.07070300e-01, -3.83837875e-02,  0.00000000e+00],
         [ 3.08689839e-01, -5.35491511e-17,  0.00000000e+00],
         [ 3.07070300e-01,  3.83837875e-02,  0.00000000e+00],
         [ 3.02205434e-01,  7.55513585e-02,  0.00000000e+00],
         [ 2.94076421e-01,  1.10278658e-01,  0.00000000e+00],
         [ 2.82651657e-01,  1.41325829e-01,  0.00000000e+00],
         [ 2.67886381e-01,  1.67428988e-01,  0.00000000e+00],
         [ 2.49722150e-01,  1.87291612e-01,  0.00000000e+00],
         [ 2.86525085e-01, -1.27344482e-01,  0.00000000e+00],
         [ 2.99520519e-01, -9.98401731e-02,  0.00000000e+00],
         [ 3.08766793e-01, -6.86148429e-02,  0.00000000e+00],
         [ 3.14300169e-01, -3.49222410e-02,  0.00000000e+00],
         [ 3.16142238e-01, -4.87483877e-17,  0.00000000e+00],
         [ 3.14300169e-01,  3.49222410e-02,  0.00000000e+00],
         [ 3.08766793e-01,  6.86148429e-02,  0.00000000e+00],
         [ 2.99520519e-01,  9.98401731e-02,  0.00000000e+00],
         [ 2.86525085e-01,  1.27344482e-01,  0.00000000e+00],
         [ 3.12152687e-01, -4.33198876e-17,  0.00000000e+00]])],
 ['#268bd2', '#859900', '#dc322f'])
In [10]:
def dump_dist(p, wi, ray_pkg, fld, wvl, foc):
    if ray_pkg is not None:
        image_pt = fld.ref_sphere[0]
        ray = ray_pkg[0]
        dist = foc / ray[-1][1][2]
        defocused_pt = ray[-1][0] + dist*ray[-1][1]
        t_abr = defocused_pt - image_pt
        return np.sqrt(np.sum(t_abr*t_abr))
        
spot_rms = np.sqrt(np.mean(np.square(sm.trace_grid(dump_dist, 0, form='list', append_if_none=False)[0]), axis=1))
print(spot_rms)
[0.25976753 0.26472785 0.27037767]
In [11]:
def spot_rms(sm):
    return np.sqrt(np.mean(np.square(sm.trace_grid(dump_dist, 0, form='list', append_if_none=False)[0]), axis=1))
'''
def calc_focus(opm, sm):
    old_gap = opm.seq_model.gaps[-1].thi
    offsets = np.linspace(-20, 20, 100)
    spots = np.zeros(offsets.shape)
    for i in range(offsets.shape[0]):
        opm.seq_model.gaps[-1].thi = opm.optical_spec.parax_data.fod.bfl + offsets[i]
        opm.update_model()
        spots[i] = np.sum(spot_rms(sm))
    opm.seq_model.gaps[-1].thi = old_gap
    opm.update_model()
    
    return offsets[np.argmin(spots)]

f = calc_focus(opm, sm)
print(f)
'''
Out[11]:
'\ndef calc_focus(opm, sm):\n    old_gap = opm.seq_model.gaps[-1].thi\n    offsets = np.linspace(-20, 20, 100)\n    spots = np.zeros(offsets.shape)\n    for i in range(offsets.shape[0]):\n        opm.seq_model.gaps[-1].thi = opm.optical_spec.parax_data.fod.bfl + offsets[i]\n        opm.update_model()\n        spots[i] = np.sum(spot_rms(sm))\n    opm.seq_model.gaps[-1].thi = old_gap\n    opm.update_model()\n    \n    return offsets[np.argmin(spots)]\n\nf = calc_focus(opm, sm)\nprint(f)\n'
In [12]:
opm.seq_model.gaps[-1].thi = opm.optical_spec.parax_data.fod.bfl - 6.6666
opm.update_model()

spot_plt = plt.figure(FigureClass=SpotDiagramFigure, opt_model=opm, is_dark=isdark).plot()
print(spot_rms(sm))
[0.16672045 0.1716864  0.1773614 ]
Figure
In [13]:
import rayoptics.optical.model_constants as mc

old_gap = opm.seq_model.gaps[-1].thi
offsets = np.linspace(-20, 20, 1000)

opm.seq_model.gaps[-1].thi = opm.optical_spec.parax_data.fod.bfl
opm.update_model()

def dump_rays(p, wi, ray_pkg, fld, wvl, foc):
    if ray_pkg is not None:
        image_pt = fld.ref_sphere[0]
        ray = ray_pkg[mc.ray]
        v = ray[-1][mc.d][0:2] / ray[-1][mc.d][2]
        return [ray[-1][mc.p][0:2] - image_pt[0:2], v]

vals, colors = sm.trace_grid(dump_rays, 0, form='list', append_if_none=False)
v1, v2, v3 = vals
v1_p = v1[:,0]
v1_v = v1[:,1]
v2_p = v2[:,0]
v2_v = v2[:,1]
v3_p = v3[:,0]
v3_v = v3[:,1]

v1_rms = np.zeros(offsets.size)
v2_rms = np.zeros(offsets.size)
v3_rms = np.zeros(offsets.size)

for i in range(offsets.size):
    v1_rms[i] = np.sqrt(np.mean(np.square(v1_p+v1_v*offsets[i])))
    v2_rms[i] = np.sqrt(np.mean(np.square(v2_p+v2_v*offsets[i])))
    v3_rms[i] = np.sqrt(np.mean(np.square(v3_p+v3_v*offsets[i])))

opm.seq_model.gaps[-1].thi = old_gap
opm.update_model()

plt.figure()
plt.plot(offsets, v1_rms, color=colors[0])
plt.plot(offsets, v2_rms, color=colors[1])
plt.plot(offsets, v3_rms, color=colors[2])
plt.show()

min_b, min_g, min_r = (offsets[np.argmin(v1_rms)], offsets[np.argmin(v2_rms)], offsets[np.argmin(v3_rms)])
min_b, min_g, min_r
Figure
Out[13]:
(-2.6626626626626617, -2.5025025025025016, -2.3823823823823815)
In [14]:
opm.seq_model.gaps[-1].thi = opm.optical_spec.parax_data.fod.bfl + min_r
opm.update_model()

spot_plt = plt.figure(FigureClass=SpotDiagramFigure, opt_model=opm, is_dark=isdark).plot()
print(spot_rms(sm))
[0.03725505 0.03611947 0.03574278]
Figure
In [15]:
opm.seq_model.gaps[-1].thi = opm.optical_spec.parax_data.fod.bfl + min_g
opm.update_model()

spot_plt = plt.figure(FigureClass=SpotDiagramFigure, opt_model=opm, is_dark=isdark).plot()
print(spot_rms(sm))
[0.03617745 0.03572421 0.03615391]
Figure
In [16]:
opm.seq_model.gaps[-1].thi = opm.optical_spec.parax_data.fod.bfl + min_b
opm.update_model()

spot_plt = plt.figure(FigureClass=SpotDiagramFigure, opt_model=opm, is_dark=isdark).plot()
print(spot_rms(sm))
[0.03573303 0.03621866 0.03768028]
Figure
In [17]:
abr_plt = plt.figure(FigureClass=RayFanFigure, opt_model=opm, data_type='Ray',
                     scale_type=Fit.All_Same).plot()
Figure
In [ ]: