#%matplotlib inline
%matplotlib notebook
#%matplotlib widget
import sys
import os
try:
hostname = os.environ['HOSTNAME']
except KeyError:
hostname = 'juno'
print('HOSTNAME:', hostname)
#os.path.exists('/home/chust/tch')
if hostname == 'juno' and not '/home/chust/tch' in sys.path:
sys.path.append('/home/chust/tch')
print('sys.path:', sys.path, '\n')
import matplotlib.pyplot as plt
import numpy as np
import glob
from lib_lfr.global_params import *
from lib_lfr.load_routines import get_lfr_decom_l1_asm
from lib_lfr.calib_routines import calibrate_sm
from lib_lfr.time_routines import index_from_date
from lib_io.store_routines import json2data
dir_root = '/WIN/Users/chust/DD CHUST' if hostname == 'PC-CHUST2.lpp.polytechnnique.fr' else \
'/home/chust/DD pc-p-chust' if hostname == 'pc-p-chust' else \
'/home/thomas/Labo' if hostname == 'pc-tango2' else 'ERROR!'
dir_DATA1 = '/DATA1' if hostname == 'PC-CHUST2.lpp.polytechnnique.fr' else \
'/home/chust/DATA' if hostname in ['pc-p-chust', 'juno'] else \
'/home/thomas/Labo/DATA' if hostname == 'pc-tango2' else 'ERROR!'
dir_DATA3 = '/run/media/chust/Samsung_T5'
dir_plots = '/home/thomas/Labo/Plots' if hostname == 'pc-tango2' else '/home/chust/Plots'
dir_storage = '/home/thomas/Labo/Storage' if hostname == 'pc-tango2' else '/home/chust/Storage'
dir_DATA = dir_DATA1 + '/SO/Tests/2022_03_04_Bruno_FSW_3.3'
dir_DATA = dir_DATA1 + '/SO/Tests/2022_03_10_Bruno_FSW_3.3'
print('dir_root:', dir_root)
print('dir_DATA:', dir_DATA)
HOSTNAME: pc-p-chust sys.path: ['/home/chust/DD pc-p-chust/Missions/Solar Orbiter/LFR/Python/Alexis', '/home/chust/DD pc-p-chust/Missions/Solar Orbiter/LFR/Python/Alexis/lfrcompliance_AJE', '/home/chust/DD pc-p-chust/Missions/Solar Orbiter/LFR/Python/tch', '/usr/lib64/python310.zip', '/usr/lib64/python3.10', '/usr/lib64/python3.10/lib-dynload', '', '/home/chust/venv_py3/lib64/python3.10/site-packages', '/home/chust/venv_py3/lib/python3.10/site-packages', '/usr/lib64/python3.10/site-packages', '/usr/lib/python3.10/site-packages'] dir_root: /home/chust/DD pc-p-chust dir_DATA: /home/chust/DATA/SO/Tests/2022_03_10_Bruno_FSW_3.3
/home/chust/venv_py3/lib64/python3.10/site-packages/spacepy/time.py:2365: UserWarning: Leapseconds may be out of date. Use spacepy.toolbox.update(leapsecs=True)
warnings.warn('Leapseconds may be out of date.'
Mapping Discospace :
Analog Discovery 1 : E1/E2
Analog Discovery 2 : B1/B2
Analog Discovery 3 : B3
files = glob.glob(dir_root + '/Missions/Solar Orbiter/LFR/Ground Segment/ROB kcoefficients' +
'/kcoefficients update #3/ASM_calibration_default-full-tables*')
for file in files:
print(os.path.basename(file))
ASM_calibration_default-full-tables_FSW_R3.3_AC_DIFF_G5_or_DC_DIFF_G1_R=0-0-0_kappa=1e-04_2021-09-08.json
[freq_cal, B_cal, E_cal] = json2data(file)
index_freq = 3*[None]
print(" freq_cal | B_cal | E_cal")
for F in [0, 1, 2]:
index_freq[F] = freq_cal[F].nonzero()[0]
print("F%s: %s %3d"%(F, freq_cal[F].shape, len(index_freq[F])), B_cal[F].shape, E_cal[F].shape)
index_freq
freq_cal | B_cal | E_cal F0: (128,) 88 (128, 3, 3) (128, 2, 2) F1: (128,) 104 (128, 3, 3) (128, 2, 2) F2: (128,) 96 (128, 3, 3) (128, 2, 2)
[array([ 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28,
29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41,
42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54,
55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67,
68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80,
81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93,
94, 95, 96, 97, 98, 99, 100, 101, 102, 103]),
array([ 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17,
18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30,
31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43,
44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56,
57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69,
70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82,
83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95,
96, 97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 107, 108]),
array([ 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18,
19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31,
32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44,
45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57,
58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70,
71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83,
84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96,
97, 98, 99, 100, 101])]
F=0
freq_cal[F][index_freq[F]]
array([1632., 1728., 1824., 1920., 2016., 2112., 2208., 2304., 2400.,
2496., 2592., 2688., 2784., 2880., 2976., 3072., 3168., 3264.,
3360., 3456., 3552., 3648., 3744., 3840., 3936., 4032., 4128.,
4224., 4320., 4416., 4512., 4608., 4704., 4800., 4896., 4992.,
5088., 5184., 5280., 5376., 5472., 5568., 5664., 5760., 5856.,
5952., 6048., 6144., 6240., 6336., 6432., 6528., 6624., 6720.,
6816., 6912., 7008., 7104., 7200., 7296., 7392., 7488., 7584.,
7680., 7776., 7872., 7968., 8064., 8160., 8256., 8352., 8448.,
8544., 8640., 8736., 8832., 8928., 9024., 9120., 9216., 9312.,
9408., 9504., 9600., 9696., 9792., 9888., 9984.])
F=1
freq_cal[F][index_freq[F]]
array([ 96., 112., 128., 144., 160., 176., 192., 208., 224.,
240., 256., 272., 288., 304., 320., 336., 352., 368.,
384., 400., 416., 432., 448., 464., 480., 496., 512.,
528., 544., 560., 576., 592., 608., 624., 640., 656.,
672., 688., 704., 720., 736., 752., 768., 784., 800.,
816., 832., 848., 864., 880., 896., 912., 928., 944.,
960., 976., 992., 1008., 1024., 1040., 1056., 1072., 1088.,
1104., 1120., 1136., 1152., 1168., 1184., 1200., 1216., 1232.,
1248., 1264., 1280., 1296., 1312., 1328., 1344., 1360., 1376.,
1392., 1408., 1424., 1440., 1456., 1472., 1488., 1504., 1520.,
1536., 1552., 1568., 1584., 1600., 1616., 1632., 1648., 1664.,
1680., 1696., 1712., 1728., 1744.])
F=2
freq_cal[F][index_freq[F]]
array([ 7., 8., 9., 10., 11., 12., 13., 14., 15., 16., 17.,
18., 19., 20., 21., 22., 23., 24., 25., 26., 27., 28.,
29., 30., 31., 32., 33., 34., 35., 36., 37., 38., 39.,
40., 41., 42., 43., 44., 45., 46., 47., 48., 49., 50.,
51., 52., 53., 54., 55., 56., 57., 58., 59., 60., 61.,
62., 63., 64., 65., 66., 67., 68., 69., 70., 71., 72.,
73., 74., 75., 76., 77., 78., 79., 80., 81., 82., 83.,
84., 85., 86., 87., 88., 89., 90., 91., 92., 93., 94.,
95., 96., 97., 98., 99., 100., 101., 102.])
B_cal[2][index_freq[F], ...]
array([[[-7.01193715e+01-3.87430971e+01j,
7.08770553e+01+4.08102044e+01j,
-7.30251013e+01-4.15923997e+01j],
[ 8.02652159e+01+4.42486630e+01j,
-1.17983937e+00-8.28774020e-01j,
-9.78426607e+01-5.59562190e+01j],
[ 4.60580066e+01+2.54956193e+01j,
8.91185498e+01+5.14343038e+01j,
5.59840682e+01+3.19725413e+01j]],
[[-5.17677864e+01-3.38375868e+01j,
5.29531867e+01+3.55105416e+01j,
-5.45868337e+01-3.60581463e+01j],
[ 6.01655366e+01+3.87369166e+01j,
-1.13310427e+00-6.12070330e-01j,
-7.30494768e+01-4.81175809e+01j],
[ 3.41707691e+01+2.22528201e+01j,
6.62215161e+01+4.47048178e+01j,
4.18446978e+01+2.75807385e+01j]],
[[-4.01541684e+01-2.99032809e+01j,
4.05368471e+01+3.16169126e+01j,
-4.15076464e+01-3.20855147e+01j],
[ 4.57244400e+01+3.40893532e+01j,
-6.65442672e-01-4.86793669e-01j,
-5.60343478e+01-4.28007312e+01j],
[ 2.62352246e+01+1.98248686e+01j,
5.07571760e+01+3.96245149e+01j,
3.20261362e+01+2.45369717e+01j]],
[[-3.11878217e+01-2.65611401e+01j,
3.15442700e+01+2.81223606e+01j,
-3.27321606e+01-2.84551063e+01j],
[ 3.59673466e+01+3.03922054e+01j,
-6.39379127e-01-3.77433628e-01j,
-4.37630725e+01-3.81101339e+01j],
[ 2.05245433e+01+1.74799656e+01j,
3.97795333e+01+3.53731994e+01j,
2.51527783e+01+2.19625481e+01j]],
[[-2.47253318e+01-2.40712004e+01j,
2.48846851e+01+2.52000780e+01j,
-2.59230810e+01-2.56040369e+01j],
[ 2.84993534e+01+2.77100394e+01j,
-4.46531622e-01-5.37392520e-01j,
-3.46252929e+01-3.44275109e+01j],
[ 1.61811475e+01+1.58663733e+01j,
3.13807266e+01+3.19348607e+01j,
1.98128851e+01+1.97064351e+01j]],
[[-1.98985764e+01-2.18552852e+01j,
2.00242331e+01+2.30150725e+01j,
-2.09273248e+01-2.35004764e+01j],
[ 2.27833278e+01+2.49962770e+01j,
-1.90027323e-01-3.81002097e-01j,
-2.78845634e+01-3.14019626e+01j],
[ 1.30655312e+01+1.44001169e+01j,
2.51305258e+01+2.90450396e+01j,
1.59453034e+01+1.80346421e+01j]],
[[-1.61675072e+01-2.00169601e+01j,
1.60762142e+01+2.10595014e+01j,
-1.68251253e+01-2.15250105e+01j],
[ 1.85277630e+01+2.29162107e+01j,
-2.34372198e-01-2.95546338e-01j,
-2.25546774e+01-2.89752172e+01j],
[ 1.04971145e+01+1.31998703e+01j,
2.02471409e+01+2.66061277e+01j,
1.28887359e+01+1.65691833e+01j]],
[[-1.30882223e+01-1.84083544e+01j,
1.30499329e+01+1.93850560e+01j,
-1.36810079e+01-1.98759736e+01j],
[ 1.49649880e+01+2.11947893e+01j,
-2.18208146e-01-3.37967577e-01j,
-1.83050196e+01-2.66151514e+01j],
[ 8.59871600e+00+1.21404317e+01j,
1.63621966e+01+2.44413473e+01j,
1.04447760e+01+1.52365267e+01j]],
[[-1.06789697e+01-1.70495890e+01j,
1.05436594e+01+1.79506945e+01j,
-1.11140244e+01-1.84412537e+01j],
[ 1.21601154e+01+1.95404003e+01j,
-1.47094595e-01-3.11576039e-01j,
-1.48902433e+01-2.46172563e+01j],
[ 7.00091261e+00+1.12105956e+01j,
1.32545052e+01+2.26507728e+01j,
8.50564329e+00+1.41759183e+01j]],
[[-8.65513100e+00-1.58303344e+01j,
8.57580430e+00+1.67545760e+01j,
-8.94822227e+00-1.70621525e+01j],
[ 9.94040199e+00+1.81091052e+01j,
-1.53602197e-01-1.83628917e-01j,
-1.21054757e+01-2.28382434e+01j],
[ 5.65986217e+00+1.04073958e+01j,
1.07232919e+01+2.11049190e+01j,
6.87910133e+00+1.31043733e+01j]],
[[-7.01546665e+00-1.47786260e+01j,
6.92456824e+00+1.55977948e+01j,
-7.27489881e+00-1.59745216e+01j],
[ 8.03630036e+00+1.69600436e+01j,
-1.76159702e-01-3.13231807e-01j,
-9.84708886e+00-2.13763645e+01j],
[ 4.59895830e+00+9.75969869e+00j,
8.62495263e+00+1.95622720e+01j,
5.56437530e+00+1.22649616e+01j]],
[[-5.67265482e+00-1.38487978e+01j,
5.48368673e+00+1.45815924e+01j,
-5.83578791e+00-1.47919008e+01j],
[ 6.49393975e+00+1.58121749e+01j,
-4.87851717e-02-3.03209119e-01j,
-7.91112631e+00-1.99349080e+01j],
[ 3.69727546e+00+9.08075565e+00j,
6.88560540e+00+1.82502466e+01j,
4.51132034e+00+1.14277483e+01j]],
[[-4.47129531e+00-1.29319520e+01j,
4.32466061e+00+1.36203019e+01j,
-4.60016265e+00-1.39551229e+01j],
[ 5.17568388e+00+1.48434293e+01j,
-7.85633704e-02-2.44938494e-01j,
-6.24190963e+00-1.87417357e+01j],
[ 2.92622550e+00+8.53955400e+00j,
5.40364005e+00+1.71447536e+01j,
3.52770760e+00+1.07453795e+01j]],
[[-3.51950103e+00-1.21575732e+01j,
3.32764523e+00+1.27960193e+01j,
-3.59265366e+00-1.31166807e+01j],
[ 4.00863024e+00+1.39166893e+01j,
-4.31501128e-02-1.94721306e-01j,
-4.83581645e+00-1.75801119e+01j],
[ 2.29501571e+00+7.98497347e+00j,
4.13536189e+00+1.60967409e+01j,
2.73211605e+00+1.00437608e+01j]],
[[-2.66810027e+00-1.14534616e+01j,
2.44354029e+00+1.20174606e+01j,
-2.71680403e+00-1.23628967e+01j],
[ 3.06562233e+00+1.30830685e+01j,
-2.48190919e-02-1.80194939e-01j,
-3.64608239e+00-1.65472515e+01j],
[ 1.70828723e+00+7.52854195e+00j,
3.02716752e+00+1.51753226e+01j,
2.09714686e+00+9.45488331e+00j]],
[[-1.92123676e+00-1.07941312e+01j,
1.74382324e+00+1.13854567e+01j,
-1.95008848e+00-1.16313523e+01j],
[ 2.22473160e+00+1.23302946e+01j,
1.53510097e-02-1.66650585e-01j,
-2.61279074e+00-1.56206457e+01j],
[ 1.23936262e+00+7.08236188e+00j,
2.11660138e+00+1.42576314e+01j,
1.49505328e+00+8.93884894e+00j]],
[[-1.30633894e+00-1.01758333e+01j,
1.09061711e+00+1.07099794e+01j,
-1.29300433e+00-1.10032372e+01j],
[ 1.50204953e+00+1.16350973e+01j,
-1.56688793e-02-1.75084644e-01j,
-1.72079555e+00-1.47252407e+01j],
[ 8.29083212e-01+6.67716300e+00j,
1.30750920e+00+1.34835280e+01j,
9.95384794e-01+8.43623078e+00j]],
[[-7.42579535e-01-9.60572864e+00j,
5.49620549e-01+1.01044079e+01j,
-7.17175056e-01-1.03751184e+01j],
[ 8.59511058e-01+1.09820451e+01j,
-4.15557111e-03-1.59825784e-01j,
-9.71748940e-01-1.39347345e+01j],
[ 4.84323154e-01+6.28820447e+00j,
6.11200876e-01+1.26990659e+01j,
5.46909837e-01+7.95778332e+00j]],
[[-2.66134156e-01-9.08602898e+00j,
3.88095442e-02+9.56137941e+00j,
-1.86082786e-01-9.83732688e+00j],
[ 3.19898759e-01+1.03894843e+01j,
-1.20544897e-02-1.44893012e-01j,
-2.60871653e-01-1.31759976e+01j],
[ 1.57866493e-01+5.95380356e+00j,
-1.35489928e-02+1.20027904e+01j,
1.23532900e-01+7.53814201e+00j]],
[[ 1.65027570e-01-8.60411537e+00j,
-3.74679966e-01+9.00954252e+00j,
2.72427637e-01-9.29494610e+00j],
[-1.96591601e-01+9.82632487e+00j,
1.30703969e-02-1.53549860e-01j,
3.16285696e-01-1.24535874e+01j],
[-1.21390489e-01+5.62860525e+00j,
-5.47645853e-01+1.13653154e+01j,
-1.94206225e-01+7.17317725e+00j]],
[[ 5.52074395e-01-8.13233164e+00j,
-7.71644774e-01+8.54146115e+00j,
6.31489238e-01-8.80978945e+00j],
[-6.08869760e-01+9.32415416e+00j,
1.60377038e-02-1.38247281e-01j,
8.39805177e-01-1.18275705e+01j],
[-3.76683156e-01+5.32992003e+00j,
-1.00892737e+00+1.07446599e+01j,
-4.93734559e-01+6.77260474e+00j]],
[[ 8.57082220e-01-7.69715209e+00j,
-1.11188474e+00+8.09866260e+00j,
9.83396318e-01-8.34693613e+00j],
[-9.86503493e-01+8.81538372e+00j,
2.15485144e-02-1.37517108e-01j,
1.31476874e+00-1.11999705e+01j],
[-5.77985861e-01+5.05131678e+00j,
-1.43352507e+00+1.01782602e+01j,
-7.59207577e-01+6.40681168e+00j]],
[[ 1.14325933e+00-7.29387311e+00j,
-1.39177597e+00+7.66231031e+00j,
1.29185287e+00-7.90736402e+00j],
[-1.31880113e+00+8.34862812e+00j,
1.86610240e-02-1.24697755e-01j,
1.72822469e+00-1.06163002e+01j],
[-7.57357970e-01+4.77121974e+00j,
-1.81642007e+00+9.62074669e+00j,
-9.96352111e-01+6.07722349e+00j]],
[[ 1.42076215e+00-6.88601151e+00j,
-1.66682921e+00+7.24592755e+00j,
1.56535711e+00-7.52267272e+00j],
[-1.61315416e+00+7.92015349e+00j,
3.53488783e-02-1.21990999e-01j,
2.09964546e+00-1.00742979e+01j],
[-9.48880795e-01+4.53170609e+00j,
-2.12436266e+00+9.11144096e+00j,
-1.22076976e+00+5.76222865e+00j]],
[[ 1.63910389e+00-6.54374965e+00j,
-1.90237155e+00+6.88373610e+00j,
1.80929179e+00-7.15792009e+00j],
[-1.89887073e+00+7.48923462e+00j,
2.44157543e-02-1.12462059e-01j,
2.40411085e+00-9.55214324e+00j],
[-1.08479374e+00+4.31247293e+00j,
-2.44968007e+00+8.63685207e+00j,
-1.40232387e+00+5.47492857e+00j]],
[[ 1.85288693e+00-6.20807719e+00j,
-2.11020636e+00+6.50506240e+00j,
2.03805990e+00-6.75831818e+00j],
[-2.11617455e+00+7.09385325e+00j,
2.17587415e-02-1.21647519e-01j,
2.72191906e+00-9.06811229e+00j],
[-1.23974605e+00+4.07513170e+00j,
-2.70474812e+00+8.20142112e+00j,
-1.59136380e+00+5.17985920e+00j]],
[[ 2.02481954e+00-5.88415049e+00j,
-2.28794090e+00+6.16724225e+00j,
2.21741107e+00-6.40402051e+00j],
[-2.32790805e+00+6.72874398e+00j,
2.40656556e-02-1.07686079e-01j,
2.98676886e+00-8.61078947e+00j],
[-1.35540376e+00+3.84774633e+00j,
-2.92102603e+00+7.75453805e+00j,
-1.73773834e+00+4.91868810e+00j]],
[[ 2.20814583e+00-5.56445466e+00j,
-2.47367869e+00+5.83571630e+00j,
2.38824010e+00-6.07685359e+00j],
[-2.50801485e+00+6.37021338e+00j,
4.00364806e-02-1.15462299e-01j,
3.22598464e+00-8.14409754e+00j],
[-1.46034196e+00+3.64524823e+00j,
-3.13709479e+00+7.34182555e+00j,
-1.87380117e+00+4.66499099e+00j]],
[[ 2.34175134e+00-5.26628932e+00j,
-2.60784579e+00+5.53649529e+00j,
2.55013396e+00-5.76591627e+00j],
[-2.67959232e+00+6.02974562e+00j,
3.57252391e-02-8.67884828e-02j,
3.42839176e+00-7.73210199e+00j],
[-1.55163373e+00+3.44991209e+00j,
-3.31488332e+00+6.93450355e+00j,
-1.97237265e+00+4.42011620e+00j]],
[[ 2.46571453e+00-4.98532868e+00j,
-2.74025187e+00+5.22239662e+00j,
2.69666264e+00-5.45272947e+00j],
[-2.82470074e+00+5.70173082e+00j,
3.62120567e-02-7.70872897e-02j,
3.61129761e+00-7.32761761e+00j],
[-1.64132245e+00+3.25413563e+00j,
-3.48042154e+00+6.54633052e+00j,
-2.08046916e+00+4.20490696e+00j]],
[[ 2.58490329e+00-4.70496561e+00j,
-2.85240678e+00+4.94352873e+00j,
2.82841070e+00-5.17970741e+00j],
[-2.97316545e+00+5.37930375e+00j,
3.65060256e-02-8.44552201e-02j,
3.79090371e+00-6.93542475e+00j],
[-1.71317337e+00+3.07231251e+00j,
-3.62174998e+00+6.19636081e+00j,
-2.18149622e+00+3.96943828e+00j]],
[[ 2.68333962e+00-4.44480919e+00j,
-2.96164164e+00+4.66184543e+00j,
2.94174466e+00-4.87179462e+00j],
[-3.08126525e+00+5.07586432e+00j,
3.15955773e-02-8.75379390e-02j,
3.92834640e+00-6.56732068e+00j],
[-1.78312647e+00+2.91079836e+00j,
-3.75283239e+00+5.84472899e+00j,
-2.26503701e+00+3.72610643e+00j]],
[[ 2.77914986e+00-4.20753900e+00j,
-3.04444030e+00+4.38795197e+00j,
3.02075636e+00-4.60436768e+00j],
[-3.17986519e+00+4.77864092e+00j,
4.14327630e-02-7.47209294e-02j,
4.05882695e+00-6.18738699e+00j],
[-1.84348082e+00+2.75053808e+00j,
-3.87228709e+00+5.51153704e+00j,
-2.33393663e+00+3.52347864e+00j]],
[[ 2.85891671e+00-3.95318718e+00j,
-3.12627210e+00+4.13455138e+00j,
3.10953418e+00-4.36568888e+00j],
[-3.26505198e+00+4.52544503e+00j,
4.90677652e-02-7.47239695e-02j,
4.17379429e+00-5.84966792e+00j],
[-1.89277562e+00+2.59025811e+00j,
-3.97749587e+00+5.19051957e+00j,
-2.38601266e+00+3.33814620e+00j]],
[[ 2.92488068e+00-3.72563055e+00j,
-3.19562620e+00+3.88873146e+00j,
3.18287471e+00-4.10569933e+00j],
[-3.35935400e+00+4.25634303e+00j,
5.07120311e-02-5.81836349e-02j,
4.27233059e+00-5.51633911e+00j],
[-1.94113471e+00+2.44000373e+00j,
-4.06908118e+00+4.88696449e+00j,
-2.45106185e+00+3.15429782e+00j]],
[[ 2.98590256e+00-3.50194085e+00j,
-3.25844436e+00+3.65568125e+00j,
3.25332253e+00-3.86395651e+00j],
[-3.43588965e+00+4.00211017e+00j,
4.07042351e-02-5.41888595e-02j,
4.36525743e+00-5.19314094e+00j],
[-1.97774181e+00+2.29211481e+00j,
-4.13212630e+00+4.59206607e+00j,
-2.50561227e+00+2.96504677e+00j]],
[[ 3.03948934e+00-3.28040314e+00j,
-3.30956338e+00+3.43193628e+00j,
3.31761488e+00-3.63636780e+00j],
[-3.49257930e+00+3.75678071e+00j,
3.97380815e-02-5.77983445e-02j,
4.44934989e+00-4.87726572e+00j],
[-2.01295702e+00+2.15191590e+00j,
-4.18842508e+00+4.30204653e+00j,
-2.55196557e+00+2.77950572e+00j]],
[[ 3.08453701e+00-3.06943400e+00j,
-3.34716545e+00+3.21436338e+00j,
3.36754156e+00-3.41070609e+00j],
[-3.53614768e+00+3.51349041e+00j,
5.28766537e-02-5.74197746e-02j,
4.51736845e+00-4.57186199e+00j],
[-2.05034185e+00+2.01872022e+00j,
-4.24870075e+00+4.01786413e+00j,
-2.59385045e+00+2.60527766e+00j]],
[[ 3.13128952e+00-2.87879342e+00j,
-3.38705306e+00+3.00060222e+00j,
3.40207457e+00-3.17930450e+00j],
[-3.58714407e+00+3.28346809e+00j,
5.14585057e-02-4.89020000e-02j,
4.56989945e+00-4.29230114e+00j],
[-2.06978813e+00+1.87303090e+00j,
-4.29461924e+00+3.74489688e+00j,
-2.62369276e+00+2.44110067e+00j]],
[[ 3.16524833e+00-2.67179236e+00j,
-3.41241135e+00+2.79644678e+00j,
3.44002205e+00-2.98297382e+00j],
[-3.62499216e+00+3.05880575e+00j,
6.01007635e-02-5.23147149e-02j,
4.62081518e+00-3.99687100e+00j],
[-2.10112646e+00+1.74539720e+00j,
-4.33455290e+00+3.48737341e+00j,
-2.64813480e+00+2.26933125e+00j]],
[[ 3.17958926e+00-2.49027426e+00j,
-3.44388169e+00+2.58957906e+00j,
3.47908108e+00-2.77773926e+00j],
[-3.64984212e+00+2.84903131e+00j,
7.32424099e-02-3.98286096e-02j,
4.64218303e+00-3.70956721e+00j],
[-2.11641764e+00+1.61135974e+00j,
-4.37283980e+00+3.23154269e+00j,
-2.66134132e+00+2.11374635e+00j]],
[[ 3.18688789e+00-2.31781447e+00j,
-3.47136953e+00+2.38792744e+00j,
3.50681595e+00-2.57531415e+00j],
[-3.66793604e+00+2.64201412e+00j,
8.05354636e-02-3.09578037e-02j,
4.66095209e+00-3.43952430e+00j],
[-2.12731611e+00+1.47959426e+00j,
-4.39873693e+00+2.98456087e+00j,
-2.67490048e+00+1.96486235e+00j]],
[[ 3.20944738e+00-2.10985565e+00j,
-3.47467432e+00+2.21793431e+00j,
3.50102303e+00-2.41486849e+00j],
[-3.68632276e+00+2.41052972e+00j,
6.82201721e-02-7.75530024e-02j,
4.74450324e+00-3.19877273e+00j],
[-2.16462878e+00+1.37266058e+00j,
-4.38888177e+00+2.76836991e+00j,
-2.71066979e+00+1.79691854e+00j]],
[[ 3.22370551e+00-1.92854963e+00j,
-3.47433252e+00+2.02461586e+00j,
3.50845013e+00-2.22484247e+00j],
[-3.69698509e+00+2.20640527e+00j,
5.90980052e-02-8.17938787e-02j,
4.76718778e+00-2.94737772e+00j],
[-2.17261272e+00+1.26078036e+00j,
-4.38441370e+00+2.51990625e+00j,
-2.71590958e+00+1.65874995e+00j]],
[[ 3.22716024e+00-1.76495040e+00j,
-3.46956866e+00+1.83957806e+00j,
3.51833533e+00-2.02970760e+00j],
[-3.69610732e+00+2.01965293e+00j,
6.53365569e-02-6.20702030e-02j,
4.75622813e+00-2.69087105e+00j],
[-2.16114334e+00+1.14324997e+00j,
-4.38237695e+00+2.28667374e+00j,
-2.71153066e+00+1.52688099e+00j]],
[[ 3.21918931e+00-1.62108219e+00j,
-3.46308646e+00+1.67546156e+00j,
3.53111386e+00-1.83168612e+00j],
[-3.68293449e+00+1.85111923e+00j,
8.74794709e-02-1.54595430e-02j,
4.71372045e+00-2.43356251e+00j],
[-2.13022290e+00+1.01764757e+00j,
-4.39029390e+00+2.08348321e+00j,
-2.70361222e+00+1.39352886e+00j]],
[[ 3.22797389e+00-1.46020556e+00j,
-3.46868248e+00+1.49984557e+00j,
3.51625395e+00-1.66359045e+00j],
[-3.68874452e+00+1.65996811e+00j,
6.21897959e-02-2.72277681e-02j,
4.72303644e+00-2.22642221e+00j],
[-2.12143187e+00+9.24838575e-01j,
-4.38064458e+00+1.85230298e+00j,
-2.70550112e+00+1.26163598e+00j]],
[[ 3.22551723e+00-1.29907219e+00j,
-3.46310167e+00+1.32802218e+00j,
3.50071197e+00-1.49508753e+00j],
[-3.68403108e+00+1.47320785e+00j,
4.72742927e-02-3.43283160e-02j,
4.72166531e+00-2.01589647e+00j],
[-2.11604129e+00+8.31928095e-01j,
-4.36256964e+00+1.63303489e+00j,
-2.70287142e+00+1.13345437e+00j]],
[[ 3.20624324e+00-1.14087162e+00j,
-3.44313236e+00+1.16236207e+00j,
3.49113183e+00-1.31753882e+00j],
[-3.66562619e+00+1.29620213e+00j,
5.49750101e-02-2.60870433e-02j,
4.69601648e+00-1.79385802e+00j],
[-2.11168741e+00+7.31789685e-01j,
-4.34124903e+00+1.43546257e+00j,
-2.69501966e+00+1.00838459e+00j]],
[[ 3.18235754e+00-9.93992947e-01j,
-3.41792910e+00+1.01165368e+00j,
3.47431520e+00-1.16639528e+00j],
[-3.64389106e+00+1.12994476e+00j,
6.17830763e-02-1.71201448e-02j,
4.66067726e+00-1.57724314e+00j],
[-2.10317028e+00+6.31158995e-01j,
-4.31881310e+00+1.23989260e+00j,
-2.67076486e+00+8.83428243e-01j]],
[[ 3.15508946e+00-8.51605683e-01j,
-3.38862607e+00+8.64065322e-01j,
3.45160477e+00-1.01816318e+00j],
[-3.61692815e+00+9.67308156e-01j,
6.44247415e-02-1.11646689e-02j,
4.62226494e+00-1.36414527e+00j],
[-2.08855798e+00+5.34377712e-01j,
-4.28666641e+00+1.04601646e+00j,
-2.64434625e+00+7.61364646e-01j]],
[[ 3.13105611e+00-7.11544730e-01j,
-3.35821630e+00+7.08870192e-01j,
3.42189004e+00-8.51684812e-01j],
[-3.58590278e+00+8.03742465e-01j,
5.36087861e-02-2.03866921e-02j,
4.59940254e+00-1.15120611e+00j],
[-2.06340956e+00+4.50141076e-01j,
-4.23407017e+00+8.52922507e-01j,
-2.63675794e+00+6.43493211e-01j]],
[[ 3.10195774e+00-5.78024168e-01j,
-3.32326876e+00+5.63821009e-01j,
3.39307419e+00-6.97687343e-01j],
[-3.55174527e+00+6.45403453e-01j,
4.71128044e-02-2.17452847e-02j,
4.55983056e+00-9.52925681e-01j],
[-2.03844420e+00+3.59490662e-01j,
-4.18646328e+00+6.65478900e-01j,
-2.61803822e+00+5.30179731e-01j]],
[[ 3.06806770e+00-4.50865683e-01j,
-3.28453023e+00+4.27311269e-01j,
3.36541994e+00-5.53566497e-01j],
[-3.51448930e+00+4.91838688e-01j,
4.36988785e-02-1.66287484e-02j,
4.50609860e+00-7.68930799e-01j],
[-2.01245285e+00+2.63246318e-01j,
-4.14258552e+00+4.82898899e-01j,
-2.58929175e+00+4.21681497e-01j]],
[[ 3.03353246e+00-3.16751141e-01j,
-3.24415350e+00+2.85582683e-01j,
3.32589868e+00-4.07136697e-01j],
[-3.46967459e+00+3.42456859e-01j,
4.49983547e-02-1.42632558e-02j,
4.45850368e+00-5.75196344e-01j],
[-1.98991144e+00+1.81122807e-01j,
-4.08907850e+00+3.14052015e-01j,
-2.56275229e+00+3.10827291e-01j]],
[[ 2.99380695e+00-1.84003240e-01j,
-3.19856211e+00+1.46174233e-01j,
3.27924589e+00-2.64413664e-01j],
[-3.41899542e+00+1.97239281e-01j,
4.70971098e-02-1.20595678e-02j,
4.40518968e+00-3.82981578e-01j],
[-1.96588590e+00+1.02582335e-01j,
-4.02906708e+00+1.50910869e-01j,
-2.53239581e+00+2.01577742e-01j]],
[[ 2.94596130e+00-5.75474139e-02j,
-3.14958170e+00+1.99945571e-02j,
3.23422661e+00-1.39977539e-01j],
[-3.37210469e+00+5.68381761e-02j,
4.92082311e-02-5.88342717e-03j,
4.33930639e+00-1.99438715e-01j],
[-1.93543513e+00+1.57158519e-02j,
-3.97084843e+00-2.12093987e-02j,
-2.48859769e+00+9.80808109e-02j]],
[[ 2.89466003e+00+6.63907166e-02j,
-3.09630205e+00-1.04909156e-01j,
3.18632357e+00-1.55953968e-02j],
[-3.31879007e+00-8.13793873e-02j,
5.04930859e-02-7.23547299e-04j,
4.26997608e+00-2.04889071e-02j],
[-1.90278628e+00-6.70385911e-02j,
-3.90593225e+00-1.86270314e-01j,
-2.44430272e+00-2.68996335e-03j]],
[[ 2.84367502e+00+1.88498886e-01j,
-3.03950671e+00-2.31713892e-01j,
3.13968284e+00+1.13918408e-01j],
[-3.25783220e+00-2.18544610e-01j,
5.11085041e-02+5.37596450e-04j,
4.20524530e+00+1.53389475e-01j],
[-1.87168145e+00-1.41583823e-01j,
-3.83690801e+00-3.36709973e-01j,
-2.40759529e+00-1.00773221e-01j]],
[[ 2.78894077e+00+3.05565642e-01j,
-2.97892321e+00-3.56645237e-01j,
3.08507230e+00+2.42572334e-01j],
[-3.19438155e+00-3.50010350e-01j,
5.11691085e-02+2.34494871e-03j,
4.13384826e+00+3.22687434e-01j],
[-1.83635502e+00-2.16262822e-01j,
-3.76181061e+00-4.86158914e-01j,
-2.36615172e+00-1.97423150e-01j]],
[[ 2.73202144e+00+4.15880885e-01j,
-2.91553096e+00-4.79782240e-01j,
3.02068086e+00+3.69800435e-01j],
[-3.13106339e+00-4.74032125e-01j,
5.03761910e-02+5.02765067e-03j,
4.05653280e+00+4.85703009e-01j],
[-1.79578586e+00-2.91198287e-01j,
-3.68048324e+00-6.34949246e-01j,
-2.31952703e+00-2.92017951e-01j]],
[[ 2.67299727e+00+5.24289707e-01j,
-2.84978149e+00-5.95081529e-01j,
2.95414121e+00+4.88514588e-01j],
[-3.06190500e+00-5.96082231e-01j,
4.87953392e-02+7.20189802e-03j,
3.97385613e+00+6.45960090e-01j],
[-1.75491732e+00-3.63664454e-01j,
-3.59623456e+00-7.80160858e-01j,
-2.27078546e+00-3.83647330e-01j]],
[[ 2.61204819e+00+6.31965599e-01j,
-2.78263473e+00-7.03151750e-01j,
2.88652800e+00+5.99550271e-01j],
[-2.98646541e+00-7.17127664e-01j,
4.63720846e-02+8.73223093e-03j,
3.88577615e+00+8.04652344e-01j],
[-1.71420934e+00-4.34421364e-01j,
-3.50937546e+00-9.23300241e-01j,
-2.22021617e+00-4.73202929e-01j]],
[[ 2.54726416e+00+7.34861556e-01j,
-2.71276707e+00-8.10398649e-01j,
2.81788509e+00+7.09877511e-01j],
[-2.90977422e+00-8.35539832e-01j,
4.39875271e-02+1.10391628e-02j,
3.79283970e+00+9.54507928e-01j],
[-1.67074445e+00-5.01844396e-01j,
-3.41812697e+00-1.06077658e+00j,
-2.16628191e+00-5.59733191e-01j]],
[[ 2.47878145e+00+8.33541540e-01j,
-2.63986030e+00-9.17693485e-01j,
2.74806049e+00+8.20654670e-01j],
[-2.83180392e+00-9.52477082e-01j,
4.15889880e-02+1.41725677e-02j,
3.69551908e+00+1.09636230e+00j],
[-1.62461370e+00-5.66268077e-01j,
-3.32255985e+00-1.19347484e+00j,
-2.10893717e+00-6.43745743e-01j]],
[[ 2.40683385e+00+9.28715687e-01j,
-2.56361662e+00-1.01971312e+00j,
2.67381596e+00+9.27390980e-01j],
[-2.74993650e+00-1.06398092e+00j,
3.99982172e-02+1.60464426e-02j,
3.59460013e+00+1.23602255e+00j],
[-1.57620619e+00-6.28584766e-01j,
-3.22315368e+00-1.32078641e+00j,
-2.04940639e+00-7.24002268e-01j]],
[[ 2.33144724e+00+1.02081865e+00j,
-2.48443262e+00-1.11665476e+00j,
2.59525668e+00+1.03033817e+00j],
[-2.66455847e+00-1.17030626e+00j,
3.96030545e-02+1.66038249e-02j,
3.49001018e+00+1.37493333e+00j],
[-1.52549290e+00-6.89130822e-01j,
-3.12022979e+00-1.44333894e+00j,
-1.98810983e+00-8.00791322e-01j]],
[[ 2.25333758e+00+1.10932881e+00j,
-2.40138061e+00-1.21018661e+00j,
2.51291502e+00+1.12898386e+00j],
[-2.57558058e+00-1.27242351e+00j,
3.87472107e-02+1.77000363e-02j,
3.38015186e+00+1.50940671e+00j],
[-1.47348055e+00-7.47023028e-01j,
-3.01297288e+00-1.56121834e+00j,
-1.92440201e+00-8.76289129e-01j]],
[[ 2.17369712e+00+1.19431132e+00j,
-2.31481636e+00-1.30037438e+00j,
2.42756488e+00+1.22277160e+00j],
[-2.48405218e+00-1.37018887e+00j,
3.69625146e-02+1.96678177e-02j,
3.26574295e+00+1.63970111e+00j],
[-1.42141372e+00-8.02236429e-01j,
-2.90235534e+00-1.67447710e+00j,
-1.85885568e+00-9.51771488e-01j]],
[[ 2.09028871e+00+1.27618757e+00j,
-2.22451838e+00-1.38756791e+00j,
2.33871309e+00+1.31354358e+00j],
[-2.38814825e+00-1.46429298e+00j,
3.49736348e-02+2.11930859e-02j,
3.14660625e+00+1.76396373e+00j],
[-1.36657073e+00-8.55305925e-01j,
-2.78648796e+00-1.78329006e+00j,
-1.79046817e+00-1.02418797e+00j]],
[[ 2.00344873e+00+1.35611860e+00j,
-2.13156504e+00-1.47438151e+00j,
2.24847518e+00+1.40372385e+00j],
[-2.28821752e+00-1.55590051e+00j,
3.25552075e-02+2.12367063e-02j,
3.02695387e+00+1.87970275e+00j],
[-1.30819912e+00-9.06354324e-01j,
-2.66519961e+00-1.88898699e+00j,
-1.72126618e+00-1.09340826e+00j]],
[[ 1.91291003e+00+1.43184497e+00j,
-2.03451448e+00-1.55670583e+00j,
2.15397567e+00+1.48964964e+00j],
[-2.18403569e+00-1.64263825e+00j,
3.01887295e-02+2.11275609e-02j,
2.90223858e+00+1.98996306e+00j],
[-1.24750865e+00-9.54556469e-01j,
-2.53903073e+00-1.98871125e+00j,
-1.64877760e+00-1.15939790e+00j]],
[[ 1.82199728e+00+1.50324995e+00j,
-1.93443388e+00-1.63108517e+00j,
2.05856315e+00+1.57101181e+00j],
[-2.07970836e+00-1.72641102e+00j,
2.90720219e-02+2.40135477e-02j,
2.76848956e+00+2.09908207e+00j],
[-1.18509563e+00-1.00246568e+00j,
-2.41670390e+00-2.08208353e+00j,
-1.57490472e+00-1.21831874e+00j]],
[[ 1.72752922e+00+1.57057171e+00j,
-1.83064550e+00-1.70078035e+00j,
1.95910862e+00+1.64810260e+00j],
[-1.97113986e+00-1.80553095e+00j,
2.77773813e-02+2.68912001e-02j,
2.62920294e+00+2.20222261e+00j],
[-1.12023239e+00-1.04750824e+00j,
-2.28980979e+00-2.17015311e+00j,
-1.49823191e+00-1.27388999e+00j]],
[[ 1.62885974e+00+1.63459295e+00j,
-1.72626324e+00-1.76660735e+00j,
1.85620699e+00+1.72032615e+00j],
[-1.86039892e+00-1.87653532e+00j,
2.77125707e-02+2.79871472e-02j,
2.49328824e+00+2.29935896e+00j],
[-1.05414798e+00-1.08818813e+00j,
-2.15666092e+00-2.25303551e+00j,
-1.41960810e+00-1.32865716e+00j]],
[[ 1.52616410e+00+1.69466649e+00j,
-1.61937987e+00-1.82837761e+00j,
1.74947116e+00+1.78798447e+00j],
[-1.74649166e+00-1.94139376e+00j,
2.82383139e-02+2.84253445e-02j,
2.35562784e+00+2.39132148e+00j],
[-9.86274657e-01-1.12536751e+00j,
-2.01793695e+00-2.33025863e+00j,
-1.33819102e+00-1.38127375e+00j]],
[[ 1.42190147e+00+1.74956183e+00j,
-1.50914213e+00-1.88556957e+00j,
1.63923716e+00+1.85094181e+00j],
[-1.62837669e+00-2.00347411e+00j,
2.74759094e-02+2.92538455e-02j,
2.21126958e+00+2.47739664e+00j],
[-9.17329213e-01-1.16038738e+00j,
-1.87617835e+00-2.40219601e+00j,
-1.25395449e+00-1.42991361e+00j]],
[[ 1.31791360e+00+1.79919268e+00j,
-1.39574612e+00-1.93832359e+00j,
1.52577785e+00+1.90922635e+00j],
[-1.50523782e+00-2.06465007e+00j,
2.42036861e-02+3.04687040e-02j,
2.05900909e+00+2.55687215e+00j],
[-8.48120711e-01-1.19427941e+00j,
-1.73258453e+00-2.46994848e+00j,
-1.16712727e+00-1.47418933e+00j]],
[[ 1.21116300e+00+1.84378466e+00j,
-1.27942231e+00-1.98543083e+00j,
1.40916194e+00+1.96186951e+00j],
[-1.37856429e+00-2.11945406e+00j,
2.08663459e-02+3.14250537e-02j,
1.90243435e+00+2.62876774e+00j],
[-7.76996818e-01-1.22476373e+00j,
-1.58509932e+00-2.53053589e+00j,
-1.07789426e+00-1.51419372e+00j]],
[[ 1.10031730e+00+1.88637444e+00j,
-1.16145534e+00-2.02975569e+00j,
1.29288800e+00+2.01055989e+00j],
[-1.25256602e+00-2.16650194e+00j,
1.91008462e-02+3.14390167e-02j,
1.74389732e+00+2.69718882e+00j],
[-7.04511457e-01-1.25111561e+00j,
-1.43592686e+00-2.58427221e+00j,
-9.86298635e-01-1.55317243e+00j]],
[[ 9.86493217e-01+1.92361024e+00j,
-1.04058908e+00-2.06843648e+00j,
1.17373921e+00+2.05376508e+00j],
[-1.12377532e+00-2.20716132e+00j,
1.74901181e-02+3.13229227e-02j,
1.58110902e+00+2.75811977e+00j],
[-6.30303505e-01-1.27380369e+00j,
-1.28326813e+00-2.63062016e+00j,
-8.92153437e-01-1.58785471e+00j]],
[[ 8.69890761e-01+1.95476288e+00j,
-9.17056722e-01-2.10112855e+00j,
1.05041274e+00+2.09184158e+00j],
[-9.91814162e-01-2.24222933e+00j,
1.59999654e-02+3.21317851e-02j,
1.41461637e+00+2.81077095e+00j],
[-5.54235786e-01-1.29360898e+00j,
-1.12640775e+00-2.66992732e+00j,
-7.95718429e-01-1.61757463e+00j]],
[[ 7.51145871e-01+1.97932989e+00j,
-7.91553750e-01-2.12746349e+00j,
9.23100789e-01+2.12432866e+00j],
[-8.57175079e-01-2.27148666e+00j,
1.45037291e-02+3.41995881e-02j,
1.24545925e+00+2.85456739e+00j],
[-4.76505303e-01-1.31049507e+00j,
-9.65842173e-01-2.70171808e+00j,
-6.97568594e-01-1.64191873e+00j]],
[[ 6.30186079e-01+1.99749708e+00j,
-6.63610194e-01-2.14714287e+00j,
7.93071988e-01+2.15006846e+00j],
[-7.19845401e-01-2.29357215e+00j,
1.27888237e-02+3.62118083e-02j,
1.07259612e+00+2.88955237e+00j],
[-3.97181164e-01-1.32321428e+00j,
-8.02250518e-01-2.72488276e+00j,
-5.97353552e-01-1.66106110e+00j]],
[[ 5.08452754e-01+2.01101382e+00j,
-5.34054212e-01-2.16119328e+00j,
6.62835615e-01+2.17220073e+00j],
[-5.80175520e-01-2.30815418e+00j,
1.14316092e-02+3.64912365e-02j,
8.96437240e-01+2.91750751e+00j],
[-3.16662076e-01-1.33172927e+00j,
-6.36908904e-01-2.74074972e+00j,
-4.97559304e-01-1.67531510e+00j]],
[[ 3.85140399e-01+2.01798300e+00j,
-4.02688191e-01-2.16820309e+00j,
5.30495138e-01+2.18770873e+00j],
[-4.38604857e-01-2.31485644e+00j,
1.01278172e-02+3.64940815e-02j,
7.17249597e-01+2.93606254e+00j],
[-2.35016595e-01-1.33570744e+00j,
-4.69427079e-01-2.74753852e+00j,
-3.96518151e-01-1.68422061e+00j]],
[[ 2.60235183e-01+2.01843745e+00j,
-2.69196322e-01-2.16807132e+00j,
3.95435246e-01+2.19558947e+00j],
[-2.95448610e-01-2.31389069e+00j,
8.59430732e-03+3.63077483e-02j,
5.35516768e-01+2.94553610e+00j],
[-1.52310695e-01-1.33497266e+00j,
-3.00203860e-01-2.74483047e+00j,
-2.93302092e-01-1.68823966e+00j]],
[[ 1.34113055e-01+2.01348041e+00j,
-1.33113609e-01-2.16122689e+00j,
2.57373937e-01+2.19486846e+00j],
[-1.51448072e-01-2.30613193e+00j,
6.41312398e-03+3.54461096e-02j,
3.52330061e-01+2.94802412e+00j],
[-6.88630653e-02-1.32935864e+00j,
-1.30648072e-01-2.73275715e+00j,
-1.86500688e-01-1.68861644e+00j]],
[[ 6.57726081e-03+2.00080196e+00j,
4.31844832e-03-2.14599375e+00j,
1.17628396e-01+2.18579413e+00j],
[-6.04962356e-03-2.28958971e+00j,
4.30839370e-03+3.44482867e-02j,
1.66356108e-01+2.93960540e+00j],
[ 1.52988642e-02-1.31860469e+00j,
4.04128603e-02-2.71034981e+00j,
-7.82790294e-02-1.68247364e+00j]],
[[-1.21818146e-01+1.97997502e+00j,
1.41985305e-01-2.12240033e+00j,
-2.21700372e-02+2.16936352e+00j],
[ 1.40735629e-01-2.26604113e+00j,
1.62102968e-03+3.40350901e-02j,
-2.21736436e-02+2.91825646e+00j],
[ 1.00214388e-01-1.30284832e+00j,
2.13642415e-01-2.67856295e+00j,
3.02176968e-02-1.67008378e+00j]],
[[-2.50590400e-01+1.95066231e+00j,
2.79483611e-01-2.09018274e+00j,
-1.61865206e-01+2.14534775e+00j],
[ 2.88735624e-01-2.23490035e+00j,
-1.73151910e-03+3.40517327e-02j,
-2.12297852e-01+2.88339730e+00j],
[ 1.85603657e-01-1.28178904e+00j,
3.88565268e-01-2.63663075e+00j,
1.38768025e-01-1.65125796e+00j]],
[[-3.79507164e-01+1.91255380e+00j,
4.17081237e-01-2.04859712e+00j,
-3.02341668e-01+2.11207039e+00j],
[ 4.37225585e-01-2.19360764e+00j,
-5.13094482e-03+3.37611867e-02j,
-4.03126777e-01+2.83581947e+00j],
[ 2.70922467e-01-1.25482246e+00j,
5.63577557e-01-2.58261242e+00j,
2.47903361e-01-1.62515030e+00j]],
[[-5.09146985e-01+1.86715072e+00j,
5.55289545e-01-1.99889461e+00j,
-4.44860898e-01+2.07006858e+00j],
[ 5.86218069e-01-2.14092945e+00j,
-7.62194819e-03+3.32036341e-02j,
-5.94678809e-01+2.77891484e+00j],
[ 3.55486315e-01-1.22332427e+00j,
7.37932200e-01-2.51762673e+00j,
3.56945579e-01-1.59148249e+00j]],
[[-6.38671578e-01+1.81274253e+00j,
6.93221727e-01-1.93947971e+00j,
-5.88041067e-01+2.01800060e+00j],
[ 7.34636144e-01-2.07674349e+00j,
-9.76403184e-03+3.24950398e-02j,
-7.86471894e-01+2.70957169e+00j],
[ 4.39446491e-01-1.18612619e+00j,
9.11242524e-01-2.44016036e+00j,
4.65721645e-01-1.54990437e+00j]],
[[-7.67356710e-01+1.74800456e+00j,
8.30139519e-01-1.86909075e+00j,
-7.30888555e-01+1.95482060e+00j],
[ 8.81813447e-01-2.00072912e+00j,
-1.19255556e-02+3.16283763e-02j,
-9.77895215e-01+2.62537991e+00j],
[ 5.22758219e-01-1.14224216e+00j,
1.08309453e+00-2.34883990e+00j,
5.74107294e-01-1.49994129e+00j]],
[[-8.94238741e-01+1.67182793e+00j,
9.65149719e-01-1.78730045e+00j,
-8.73570191e-01+1.87954229e+00j],
[ 1.02653988e+00-1.91391223e+00j,
-1.48557830e-02+3.07178518e-02j,
-1.16854306e+00+2.52556288e+00j],
[ 6.06034158e-01-1.09159839e+00j,
1.25417883e+00-2.24402911e+00j,
6.82777307e-01-1.44170764e+00j]]])
case = 'F0-3168Hz-2Vpp_on_E1E2B1B2B3-UnitKCOEFFMatrices'
dir_Packets = dir_DATA + '/' + case + '/decom'
files = glob.glob(dir_Packets + '/*.af0')
files.sort()
print(len(files))
for file in files:
print(os.path.basename(file))
asm_unit = get_lfr_decom_l1_asm(file)
1 2022_03_10_15_12_03_packet_record_NORMAL.af0
case = 'F0-3168Hz-2Vpp_on_E1E2B1B2B3-DefaultKCOEFFMatrices'
dir_Packets = dir_DATA + '/' + case + '/decom'
files = glob.glob(dir_Packets + '/*.af0')
files.sort()
print(len(files))
for file in files:
print(os.path.basename(file))
asm_default = get_lfr_decom_l1_asm(file)
1 2022_03_10_15_08_20_packet_record_NORMAL.af0
print('asm mode :', asm_unit['mode'], '@F%d'%(LFR_Fs2F(asm_unit['fe'])))
print('asm time :', asm_unit['relative_time'])
print('asm shape:', asm_unit['asm_idl'].shape)
nspec, nfreq, ndim, ndim = asm_unit['asm_idl'].shape
asm mode : NORMAL @F0 asm time : [ 0. 3.99992371 7.99984741 11.99977112 15.99969482 19.99961853 23.99954224 27.99946594 31.99938965 35.99931335 39.99923706 43.99916077 47.99908447 51.99899292 55.99891663 59.99884033] asm shape: (16, 88, 5, 5)
print('asm mode :', asm_default['mode'], '@F%d'%(LFR_Fs2F(asm_default['fe'])))
print('asm time :', asm_default['relative_time'])
print('asm shape:', asm_default['asm_idl'].shape)
asm mode : NORMAL @F0 asm time : [ 0. 3.99993896 7.99984741 11.99977112 15.99969482 19.99961853 23.99954224 27.9994812 31.99940491 35.99932861 39.99925232 43.99919128 47.99911499 51.9990387 55.9989624 59.99890137] asm shape: (16, 88, 5, 5)
F = 0
calib_mat_freqs_asm_BE = np.zeros((nfreq, ndim, ndim), dtype=complex)
calib_mat_freqs_asm_BE[:, 0:3, 0:3] = B_cal[F][index_freq[F], ...]
calib_mat_freqs_asm_BE[:, 3:5, 3:5] = E_cal[F][index_freq[F], ...]
asm_unit_calibrated = calibrate_sm(asm_unit['asm_idl'], calib_mat_freqs_asm_BE, echo=True)
calibrated sm shape: (16, 88, 5, 5)
i1 = index_from_date(3168, asm_unit['freq'], epsi=None)
i2 = index_from_date(3168, asm_default['freq'], epsi=None)
i1, asm_unit['freq'][i1], i2, asm_default['freq'][i1]
(16, 3168.0, 16, 3168.0)
np.array(asm_unit['freq'])
array([1632., 1728., 1824., 1920., 2016., 2112., 2208., 2304., 2400.,
2496., 2592., 2688., 2784., 2880., 2976., 3072., 3168., 3264.,
3360., 3456., 3552., 3648., 3744., 3840., 3936., 4032., 4128.,
4224., 4320., 4416., 4512., 4608., 4704., 4800., 4896., 4992.,
5088., 5184., 5280., 5376., 5472., 5568., 5664., 5760., 5856.,
5952., 6048., 6144., 6240., 6336., 6432., 6528., 6624., 6720.,
6816., 6912., 7008., 7104., 7200., 7296., 7392., 7488., 7584.,
7680., 7776., 7872., 7968., 8064., 8160., 8256., 8352., 8448.,
8544., 8640., 8736., 8832., 8928., 9024., 9120., 9216., 9312.,
9408., 9504., 9600., 9696., 9792., 9888., 9984.])
for i in range(5):
for j in range(5):
if i <= j:
print("i=%d, j=%d :"%(i, j))
print(np.abs(asm_unit['asm_idl'][:, i1, i, j]))
print(np.angle(asm_unit['asm_idl'][:, i1, i, j], deg=True))
i=0, j=0 :
[1.7031250e+00 2.0179524e+07 2.0183552e+07 2.0189182e+07 2.0192142e+07
2.0189408e+07 2.0183824e+07 2.0179234e+07 2.0180634e+07 2.0185390e+07
2.0190542e+07 2.0191966e+07 2.0188116e+07 2.0181438e+07 2.0179392e+07
2.0181008e+07]
[0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.]
i=0, j=1 :
[5.9302425e-01 2.0704750e+07 2.0709092e+07 2.0715752e+07 2.0718586e+07
2.0716286e+07 2.0709424e+07 2.0704552e+07 2.0706194e+07 2.0711516e+07
2.0717270e+07 2.0718516e+07 2.0714076e+07 2.0707248e+07 2.0704624e+07
2.0707044e+07]
[-10.37236 0.36110428 0.36207673 0.3620051 0.3639064
0.36444402 0.36382058 0.36269024 0.36065224 0.363144
0.36456618 0.36435312 0.3650324 0.36174795 0.3609276
0.36296123]
i=0, j=2 :
[5.7709765e-01 1.9743272e+07 1.9766876e+07 1.9783888e+07 1.9760480e+07
1.9781506e+07 1.9769554e+07 1.9752470e+07 1.9757908e+07 1.9791438e+07
1.9777100e+07 1.9760246e+07 1.9815364e+07 1.9765490e+07 1.9803520e+07
1.9770006e+07]
[ -33.403233 31.458881 89.91196 148.39124 -153.46336 -95.138
-36.7243 21.807304 80.48189 139.30005 -163.25691 -104.54899
-46.868984 10.970469 68.92569 126.64537 ]
i=0, j=3 :
[9.9144906e-01 1.9214894e+07 1.9222062e+07 1.9247918e+07 1.9239050e+07
1.9237340e+07 1.9251444e+07 1.9250224e+07 1.9246500e+07 1.9278356e+07
1.9276824e+07 1.9277470e+07 1.9285338e+07 1.9269530e+07 1.9294014e+07
1.9245966e+07]
[ 12.283956 38.128067 113.342514 -171.49811 -96.42323
-21.454256 53.442364 127.913284 -157.26578 -82.64235
-8.607178 65.26556 139.28949 -146.63483 -72.948364
1.0895917]
i=0, j=4 :
[8.0332106e-01 1.9280230e+07 1.9287762e+07 1.9314208e+07 1.9304714e+07
1.9302952e+07 1.9317626e+07 1.9315810e+07 1.9311984e+07 1.9344248e+07
1.9342822e+07 1.9342960e+07 1.9351404e+07 1.9335162e+07 1.9360106e+07
1.9311410e+07]
[ -7.2632456 38.076633 113.29046 -171.54974 -96.47696
-21.507225 53.390293 127.86016 -157.31903 -82.69299
-8.659789 65.213036 139.23843 -146.687 -72.99957
1.0375913]
i=1, j=1 :
[8.7760419e-01 2.1243642e+07 2.1248314e+07 2.1256056e+07 2.1258758e+07
2.1256920e+07 2.1248712e+07 2.1243544e+07 2.1245446e+07 2.1251360e+07
2.1257730e+07 2.1258804e+07 2.1253746e+07 2.1246756e+07 2.1243530e+07
2.1246796e+07]
[0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.]
i=1, j=2 :
[2.2729467e-01 2.0257190e+07 2.0281540e+07 2.0299876e+07 2.0275634e+07
2.0297802e+07 2.0284292e+07 2.0266642e+07 2.0272500e+07 2.0307372e+07
2.0293146e+07 2.0275522e+07 2.0331568e+07 2.0280500e+07 2.0318970e+07
2.0285440e+07]
[ -18.01977 31.097504 89.54979 148.02888 -153.82732 -95.50284
-37.08816 21.444685 80.12145 138.93723 -163.62178 -104.91341
-47.23393 10.608889 68.5645 126.282646]
i=1, j=3 :
[4.8516539e-01 1.9715064e+07 1.9722542e+07 1.9749920e+07 1.9740598e+07
1.9739456e+07 1.9752678e+07 1.9751312e+07 1.9747778e+07 1.9780926e+07
1.9779846e+07 1.9780144e+07 1.9787720e+07 1.9771636e+07 1.9796206e+07
1.9747772e+07]
[ 17.170563 37.76664 112.980316 -171.86058 -96.78721
-21.819237 53.078495 127.5507 -157.62617 -83.00507
-8.972123 64.901115 138.92456 -146.99635 -73.30964
0.72694427]
i=1, j=4 :
[3.8024402e-01 1.9782106e+07 1.9789954e+07 1.9817940e+07 1.9807976e+07
1.9806784e+07 1.9820582e+07 1.9818608e+07 1.9814966e+07 1.9848540e+07
1.9847566e+07 1.9847348e+07 1.9855504e+07 1.9838972e+07 1.9864014e+07
1.9814920e+07]
[ 0.7848246 37.71519 112.92827 -171.91219 -96.84093
-21.872204 53.026417 127.49756 -157.67941 -83.05571
-9.024735 64.84861 138.87352 -147.04854 -73.36085
0.67494255]
i=2, j=2 :
[7.5260419e-01 2.1127778e+07 2.1138100e+07 2.1132304e+07 2.1126594e+07
2.1136940e+07 2.1135598e+07 2.1126382e+07 2.1135346e+07 2.1136678e+07
2.1126442e+07 2.1132216e+07 2.1138836e+07 2.1128024e+07 2.1130430e+07
2.1138548e+07]
[0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.]
i=2, j=3 :
[4.4750011e-01 2.1090016e+07 2.1095094e+07 2.1088670e+07 2.1091348e+07
2.1093724e+07 2.1094980e+07 2.1097010e+07 2.1100818e+07 2.1095634e+07
2.1097174e+07 2.1104684e+07 2.1091988e+07 2.1100152e+07 2.1093414e+07
2.1096636e+07]
[ 45.707317 6.668669 23.432804 40.114117 57.0336 73.681725
90.16245 106.10529 122.252106 138.05653 154.63943 169.82352
-173.84578 -157.61221 -141.869 -125.55751 ]
i=2, j=4 :
[3.5563803e-01 2.1161804e+07 2.1167162e+07 2.1161270e+07 2.1163420e+07
2.1165624e+07 2.1167504e+07 2.1168880e+07 2.1172702e+07 2.1167720e+07
2.1169364e+07 2.1176380e+07 2.1164114e+07 2.1172056e+07 2.1165446e+07
2.1168458e+07]
[ 29.38061 6.6171823 23.380793 40.06244 56.98017
73.62894 90.11077 106.052345 122.19905 138.00569
154.5871 169.77104 -173.89696 -157.66475 -141.92021
-125.60975 ]
i=3, j=3 :
[1.7682291e+00 2.1197050e+07 2.1206518e+07 2.1196406e+07 2.1199898e+07
2.1206410e+07 2.1194606e+07 2.1203302e+07 2.1204554e+07 2.1194908e+07
2.1205346e+07 2.1202694e+07 2.1195884e+07 2.1206712e+07 2.1199554e+07
2.1197814e+07]
[0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.]
i=3, j=4 :
[5.9676337e-01 2.1269218e+07 2.1278954e+07 2.1269370e+07 2.1272366e+07
2.1278680e+07 2.1267482e+07 2.1275536e+07 2.1276818e+07 2.1267322e+07
2.1277902e+07 2.1274724e+07 2.1268330e+07 2.1278986e+07 2.1271902e+07
2.1270000e+07]
[-22.312138 -0.05151056 -0.05198732 -0.05172055 -0.0533449
-0.05274652 -0.0515978 -0.0528972 -0.05299045 -0.05089807
-0.05226056 -0.05242247 -0.05121182 -0.05263494 -0.05121881
-0.0522968 ]
i=4, j=4 :
[1.2630209e+00 2.1341638e+07 2.1351638e+07 2.1342586e+07 2.1345084e+07
2.1351202e+07 2.1340602e+07 2.1348012e+07 2.1349320e+07 2.1339984e+07
2.1350706e+07 2.1346996e+07 2.1341030e+07 2.1351510e+07 2.1344498e+07
2.1342438e+07]
[0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.]
for i in range(5):
for j in range(5):
if i <= j:
print("i=%d, j=%d :"%(i, j))
print(np.abs(asm_default['asm_idl'][:, i1, i, j]))
print(np.angle(asm_default['asm_idl'][:, i1, i, j], deg=True))
i=0, j=0 :
[6.0300684e-01 2.7208230e+06 2.9792838e+06 3.5156168e+06 3.8526488e+06
3.6913698e+06 3.1736012e+06 2.7585242e+06 2.8355200e+06 3.3397178e+06
3.7951772e+06 3.7894180e+06 3.3292568e+06 2.8351000e+06 2.7570028e+06
3.1619358e+06]
[0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.]
i=0, j=1 :
[1.8670937e-01 2.2111082e+06 5.6775295e+06 8.3674305e+06 8.9717640e+06
6.9334340e+06 3.1790042e+06 8.4584512e+05 4.4857480e+06 7.6477750e+06
9.0760120e+06 7.7786495e+06 4.2989790e+06 2.9236531e+05 3.4745610e+06
6.7900705e+06]
[-170.60301 88.499405 53.263477 24.17052 -1.219734 -25.631691
-53.935574 106.99503 64.81831 33.07858 6.535297 -18.087109
-45.3228 -104.8714 74.95839 42.289257]
i=0, j=2 :
[3.1509805e-01 9.5905920e+06 8.5894420e+06 6.8007180e+06 6.0088140e+06
7.8600685e+06 9.5320770e+06 9.7526180e+06 9.0234610e+06 7.4164825e+06
5.9302815e+06 7.1595760e+06 9.2105010e+06 9.8196800e+06 9.3117390e+06
8.0239040e+06]
[-178.20557 162.95738 113.55418 57.429207 -17.206724
-82.57788 -131.75385 -179.18323 131.34766 77.3231
7.6362095 -65.45099 -118.76032 -165.7516 145.67683
94.9096 ]
i=0, j=3 :
[1.3396013e-02 2.7960084e+05 2.9221272e+05 3.1655553e+05 3.3110522e+05
3.2409381e+05 3.0111512e+05 2.8144016e+05 2.8513441e+05 3.0876819e+05
3.2873969e+05 3.2850875e+05 3.0823666e+05 2.8517119e+05 2.8124384e+05
3.0057116e+05]
[ 174.0997 -167.85533 -154.94177 -137.14563 -115.35638
-92.478676 -71.74305 -56.71396 -43.329956 -26.500013
-5.5114694 16.846052 38.07099 54.785686 68.54809
84.31246 ]
i=0, j=4 :
[2.2647202e-02 4.1843472e+05 4.3731903e+05 4.7374250e+05 4.9552050e+05
4.8503328e+05 4.5063659e+05 4.2118903e+05 4.2672609e+05 4.6209225e+05
4.9198703e+05 4.9165325e+05 4.6129547e+05 4.2678400e+05 4.2089947e+05
4.4982069e+05]
[-167.13737 -167.82002 -154.90683 -137.11162 -115.32064 -92.44356
-71.70733 -56.677975 -43.295086 -26.465265 -5.477606 16.880732
38.10577 54.821476 68.58252 84.34646 ]
i=1, j=1 :
[4.0986109e-01 2.2805622e+06 1.1280543e+07 2.0310022e+07 2.1243712e+07
1.3380788e+07 3.6223288e+06 7.6218031e+05 7.5938600e+06 1.7943674e+07
2.2068472e+07 1.6338912e+07 5.9693985e+06 5.1593788e+05 4.8815460e+06
1.5009232e+07]
[0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.]
i=1, j=2 :
[3.2522708e-01 8.1725965e+06 1.6907216e+07 1.6898072e+07 1.4678340e+07
1.5142503e+07 9.7276050e+06 3.4785425e+06 1.4767572e+07 1.7665464e+07
1.4932098e+07 1.5193444e+07 1.2119222e+07 1.7067485e+06 1.2173139e+07
1.7820406e+07]
[-19.53963 67.82068 57.58756 31.882153 -15.048239 -54.51625
-72.54823 57.491364 62.971928 42.35793 1.2157724 -45.123344
-69.45898 -28.914188 66.20603 50.468258 ]
i=1, j=3 :
[9.1188028e-03 2.2462816e+05 5.4837281e+05 7.4219925e+05 7.6014062e+05
5.9868000e+05 2.9452800e+05 8.5648734e+04 4.4365341e+05 6.9620819e+05
7.7505231e+05 6.6401050e+05 3.8989469e+05 2.3840080e+04 3.4881959e+05
6.3531138e+05]
[ -30.605593 105.20829 152.11131 -161.35605 -114.38498
-67.40921 -19.29296 -158.60484 -107.502914 -59.531944
-12.229584 34.50467 82.34084 145.71976 -5.3792906
42.19379 ]
i=1, j=4 :
[1.54731460e-02 3.36163750e+05 8.20683562e+05 1.11074125e+06
1.13759825e+06 8.95974375e+05 4.40776469e+05 1.28182883e+05
6.63963625e+05 1.04192075e+06 1.15993088e+06 9.93773062e+05
5.83501250e+05 3.56806406e+04 5.22032500e+05 9.50777688e+05]
[ -12.373068 105.24408 152.14607 -161.32217 -114.34943 -67.37406
-19.256996 -158.56961 -107.46829 -59.497307 -12.195631 34.539394
82.3759 145.75511 -5.344852 42.227795]
i=2, j=2 :
[1.0608796e+00 3.5867556e+07 2.6728834e+07 1.4838472e+07 1.0867687e+07
1.8261544e+07 3.0496666e+07 3.6622932e+07 3.0836150e+07 1.8305264e+07
1.0816322e+07 1.5110531e+07 2.7264250e+07 3.6082300e+07 3.3593476e+07
2.2186218e+07]
[0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.]
i=2, j=3 :
[2.2601822e-02 9.7180388e+05 8.2772512e+05 5.9291775e+05 4.9288581e+05
6.6936306e+05 8.8665069e+05 9.7951650e+05 8.9127150e+05 6.6810500e+05
4.9056019e+05 5.9929994e+05 8.3459231e+05 9.7144519e+05 9.3343594e+05
7.4622356e+05]
[ -2.7997139 29.121058 90.723335 164.22049 -98.02572
-8.865092 60.804245 122.66501 -175.26044 -105.004234
-13.6365385 83.25901 157.7302 -139.07834 -77.48632
-11.620055 ]
i=2, j=4 :
[3.8782604e-02 1.4543400e+06 1.2387609e+06 8.8733612e+05 7.3763038e+05
1.0017591e+06 1.3269279e+06 1.4659005e+06 1.3338645e+06 9.9986475e+05
7.3416362e+05 8.9692769e+05 1.2490254e+06 1.4538541e+06 1.3969476e+06
1.1167634e+06]
[ 10.9002285 29.156284 90.75823 164.25421 -97.99041
-8.829959 60.8402 122.70137 -175.22556 -104.96979
-13.602385 83.29366 157.76508 -139.04272 -77.45175
-11.585998 ]
i=3, j=3 :
[1.8012272e-03 2.8830734e+04 2.8843389e+04 2.8830459e+04 2.8833572e+04
2.8843350e+04 2.8829178e+04 2.8838201e+04 2.8842193e+04 2.8828080e+04
2.8839258e+04 2.8840398e+04 2.8828271e+04 2.8840025e+04 2.8840398e+04
2.8828742e+04]
[0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.]
i=3, j=4 :
[1.7697049e-03 4.3146496e+04 4.3166348e+04 4.3146352e+04 4.3151387e+04
4.3166441e+04 4.3144562e+04 4.3157746e+04 4.3164566e+04 4.3143184e+04
4.3160426e+04 4.3163105e+04 4.3143266e+04 4.3161637e+04 4.3161496e+04
4.3143723e+04]
[11.466292 0.03533728 0.03499222 0.03410931 0.03582298 0.03506905
0.03563844 0.03589203 0.03491805 0.03482252 0.03380174 0.03464787
0.03468239 0.03582527 0.03441754 0.03397954]
i=4, j=4 :
[3.7607653e-03 6.4570676e+04 6.4601777e+04 6.4570859e+04 6.4578957e+04
6.4602145e+04 6.4568387e+04 6.4587652e+04 6.4599105e+04 6.4566719e+04
6.4593281e+04 6.4598750e+04 6.4566547e+04 6.4595207e+04 6.4593941e+04
6.4566844e+04]
[0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.]
for i in range(5):
for j in range(5):
if i <= j:
print("i=%d, j=%d :"%(i, j))
print(np.abs(asm_unit_calibrated[:, i1, i, j]))
print(np.angle(asm_unit_calibrated[:, i1, i, j], deg=True))
i=0, j=0 : [5.31237704e-01 2.72449918e+06 3.01322350e+06 3.56766692e+06 3.85853227e+06 3.61564156e+06 3.06294612e+06 2.72860867e+06 2.93595543e+06 3.48621795e+06 3.84603185e+06 3.68660161e+06 3.15833709e+06 2.75249759e+06 2.85288917e+06 3.36454568e+06] [-2.99353252e-15 -4.89639098e-15 0.00000000e+00 1.49567922e-14 -3.45732841e-15 1.10687520e-14 0.00000000e+00 0.00000000e+00 4.54373833e-15 3.82655743e-15 -6.93713092e-15 1.08556996e-14 -1.26714276e-14 4.84658489e-15 1.40281088e-14 -3.96493747e-15] i=0, j=1 : [1.95869509e-01 2.31573407e+06 5.89368537e+06 8.54252927e+06 8.83414230e+06 6.35085433e+06 2.33270179e+06 1.69362409e+06 5.34702056e+06 8.25289678e+06 8.99373476e+06 6.89944019e+06 3.06224246e+06 9.85872114e+05 4.66114192e+06 7.75520439e+06] [-1.56101804e+02 8.66691862e+01 5.10006261e+01 2.14344496e+01 -4.24182456e+00 -3.00296783e+01 -6.08668689e+01 9.37335197e+01 5.62942312e+01 2.57051407e+01 -3.12770817e-02 -2.56755356e+01 -5.48834927e+01 1.03664625e+02 6.31282968e+01 3.18138043e+01] i=0, j=2 : [3.10301511e-01 9.52927684e+06 8.46618057e+06 6.60171233e+06 6.10376425e+06 8.21090871e+06 9.67404597e+06 9.63150226e+06 8.69883760e+06 6.89995530e+06 5.95225446e+06 7.83379378e+06 9.56005592e+06 9.72244701e+06 8.97004890e+06 7.32395433e+06] [ 176.11568126 161.47907274 110.05771103 50.70287435 -26.81862206 -91.83793185 -141.42085129 169.75731992 118.56637234 61.12150428 -13.46870278 -82.94914042 -133.09564857 178.97531658 128.81153395 74.72278889] i=0, j=3 : [1.12032263e-02 2.79853965e+05 2.93914182e+05 3.19110429e+05 3.31670810e+05 3.21227369e+05 2.96275848e+05 2.80137801e+05 2.90298834e+05 3.15608299e+05 3.31263386e+05 3.24541044e+05 3.00674823e+05 2.81352521e+05 2.86203283e+05 3.10166557e+05] [-176.00336982 -168.17985209 -155.62639294 -138.62565866 -117.85282955 -97.16313958 -79.90086235 -67.85161393 -56.39541889 -41.04401808 -21.00787382 -1.53187812 16.31623687 29.40027026 40.15874017 55.02137462] i=0, j=4 : [1.92435339e-02 4.18831600e+05 4.39870535e+05 4.77570719e+05 4.96376725e+05 4.80749690e+05 4.43398439e+05 4.19255618e+05 4.34462600e+05 4.72337501e+05 4.95765347e+05 4.85712783e+05 4.49987644e+05 4.21073098e+05 4.28331009e+05 4.64197371e+05] [-161.23210232 -168.14556436 -155.59175245 -138.59123693 -117.81723908 -97.12799582 -79.86645638 -67.81637499 -56.36011129 -41.01015865 -20.97300702 -1.49695992 16.35030969 29.4352867 40.19285846 55.05614089] i=1, j=1 : [4.33183076e-01 2.51546363e+06 1.20137807e+07 2.08572731e+07 2.06075537e+07 1.15550459e+07 2.25278644e+06 1.59168473e+06 1.02403894e+07 1.99460736e+07 2.14070568e+07 1.33130956e+07 3.40949061e+06 8.82405971e+05 8.10469585e+06 1.83090447e+07] [ 0.00000000e+00 -7.95492313e-15 -2.22081850e-15 0.00000000e+00 0.00000000e+00 -4.61796977e-15 -8.88247526e-15 8.38119066e-15 7.81623399e-15 0.00000000e+00 -2.49267582e-15 0.00000000e+00 1.17380114e-14 3.77949993e-15 3.29197134e-15 2.91445314e-15] i=1, j=2 : [3.71629242e-01 8.52686215e+06 1.71454070e+07 1.65595586e+07 1.46871975e+07 1.47876930e+07 7.53497589e+06 6.40055980e+06 1.64024838e+07 1.70560856e+07 1.46631744e+07 1.50835481e+07 9.44473790e+06 3.96018703e+06 1.51514213e+07 1.75811193e+07] [-24.48846382 67.59608952 56.2866683 28.01214199 -21.23018327 -58.90388978 -73.09789779 66.44640095 59.16986935 33.91244351 -12.56003082 -54.53169046 -72.74660652 60.10449607 62.29996516 41.05345096] i=1, j=3 : [1.14882735e-02 2.35178611e+05 5.66164975e+05 7.53451259e+05 7.48872844e+05 5.54908053e+05 2.19828222e+05 1.72422618e+05 5.21076096e+05 7.37019080e+05 7.64356528e+05 5.98358024e+05 2.85011789e+05 1.00406085e+05 4.60986954e+05 7.04872922e+05] [ -27.72491025 106.60756093 153.62806354 -160.1365426 -113.87585172 -67.71650072 -20.76150783 -159.55915909 -112.37837322 -66.79324991 -21.20024833 23.67749195 69.85698273 -70.56359623 -22.51530318 23.21606235] i=1, j=4 : [1.92240978e-02 3.51970216e+05 8.47319496e+05 1.12759120e+06 1.12075987e+06 8.30475514e+05 3.28986631e+05 2.58049686e+05 7.79847012e+05 1.10301789e+06 1.14392710e+06 8.95510835e+05 4.26546215e+05 1.50265759e+05 6.89909023e+05 1.05491849e+06] [ -13.3667651 106.64160967 153.6627435 -160.10214424 -113.84015279 -67.6813282 -20.72694029 -159.52380114 -112.34306869 -66.7594557 -21.16531025 23.71244336 69.89071847 -70.52906266 -22.48093069 23.25073203] i=2, j=2 : [1.12164946e+00 3.56617108e+07 2.58591592e+07 1.39332395e+07 1.12827341e+07 2.03506847e+07 3.25847560e+07 3.63011367e+07 2.79141245e+07 1.53996208e+07 1.08132133e+07 1.83551668e+07 3.08150211e+07 3.65977656e+07 3.02885997e+07 1.77899040e+07] [ 0.00000000e+00 0.00000000e+00 0.00000000e+00 -3.82975207e-15 -4.72942572e-15 0.00000000e+00 -3.27520346e-15 2.93989984e-15 0.00000000e+00 -6.93015156e-15 0.00000000e+00 0.00000000e+00 -3.46330140e-15 0.00000000e+00 0.00000000e+00 0.00000000e+00] i=2, j=3 : [2.27247779e-02 9.65308771e+05 8.11441535e+05 5.71751455e+05 5.02057850e+05 7.10916425e+05 9.20931225e+05 9.75861778e+05 8.47036462e+05 6.07683508e+05 4.90557937e+05 6.71821772e+05 8.94567513e+05 9.80717755e+05 8.86873050e+05 6.59150320e+05] [ 1.76475242 30.2534173 93.4735016 169.52342259 -90.68405831 -4.3062831 62.18541925 122.43071613 -175.66964162 -103.30322698 -7.50944855 82.38484256 150.14568253 -149.39847615 -89.20986037 -20.8348469 ] i=2, j=4 : [3.96059307e-02 1.44469091e+06 1.21439814e+06 8.55664750e+05 7.51380381e+05 1.06395880e+06 1.37824256e+06 1.46047893e+06 1.26768121e+06 9.09455910e+05 7.34162776e+05 1.00545802e+06 1.33879782e+06 1.46774725e+06 1.32728378e+06 9.86492748e+05] [ 14.53354572 30.28769 93.50814399 169.55776199 -90.64822613 -4.27098596 62.22003914 122.46606858 -175.63416809 -103.26953925 -7.47431638 82.41981706 150.17973358 -149.36366568 -89.17583692 -20.80016892] i=3, j=3 : [1.70624083e-03 2.88308569e+04 2.88443661e+04 2.88321376e+04 2.88355122e+04 2.88437771e+04 2.88294574e+04 2.88394677e+04 2.88412347e+04 2.88286225e+04 2.88431071e+04 2.88380951e+04 2.88300356e+04 2.88441932e+04 2.88347205e+04 2.88319377e+04] [ 0.00000000e+00 -3.61489137e-15 -3.61319834e-15 0.00000000e+00 0.00000000e+00 -3.61327212e-15 0.00000000e+00 0.00000000e+00 0.00000000e+00 3.61517155e-15 0.00000000e+00 0.00000000e+00 3.61499436e-15 0.00000000e+00 0.00000000e+00 0.00000000e+00] i=3, j=4 : [1.62797119e-03 4.31484486e+04 4.31683559e+04 4.31492944e+04 4.31550268e+04 4.31676888e+04 4.31454037e+04 4.31612869e+04 4.31639023e+04 4.31447699e+04 4.31663053e+04 4.31595062e+04 4.31468382e+04 4.31683129e+04 4.31540062e+04 4.31500457e+04] [10.83642442 0.03430144 0.03461879 0.03444085 0.03552279 0.03512444 0.03435915 0.0352248 0.03528689 0.03389344 0.03480068 0.03490878 0.03410236 0.03505014 0.03410707 0.034825 ] i=4, j=4 : [3.56111792e-03 6.45762567e+04 6.46055771e+04 6.45759122e+04 6.45855179e+04 6.46049111e+04 6.45702599e+04 6.45953877e+04 6.45992527e+04 6.45702431e+04 6.46022579e+04 6.45931290e+04 6.45732766e+04 6.46058391e+04 6.45842312e+04 6.45786154e+04] [-3.48880380e-15 -3.22782463e-15 3.22635972e-15 0.00000000e+00 -3.22736177e-15 0.00000000e+00 -6.45624880e-15 3.22686865e-15 3.22667559e-15 -3.22812524e-15 -3.22652549e-15 0.00000000e+00 -3.22797359e-15 0.00000000e+00 0.00000000e+00 -3.22770673e-15]
i = 0
j = 1
print(np.abs(asm_unit_calibrated[:, i1, i, j] / asm_default['asm_idl'][:, i2, i, j]))
print(np.angle(asm_unit_calibrated[:, i1, i, j] / asm_default['asm_idl'][:, i2, i, j], deg=True))
[1.04906095 1.04731825 1.0380722 1.02092622 0.98466059 0.91597529 0.73378378 2.00228623 1.19200201 1.07912388 0.99093469 0.88697144 0.71231856 3.3720556 1.34150529 1.14213903] [ 14.50119637 -1.83021932 -2.26285199 -2.73607096 -3.02209065 -4.39798843 -6.93129566 -13.26150887 -8.52407923 -7.37343643 -6.56657411 -7.58842702 -9.5606926 -151.4639766 -11.83009609 -10.47545176]
k0 = 1
for i in range(5):
for j in range(5):
if i <= j:
print("i=%d, j=%d => %8.5f %8.5f %5.3f %5.3f"%(i, j,
np.min(np.abs(asm_unit_calibrated[k0:, i1, i, j] / asm_default['asm_idl'][k0:, i2, i, j])),
np.max(np.abs(asm_unit_calibrated[k0:, i1, i, j] / asm_default['asm_idl'][k0:, i2, i, j])),
np.min(abs(np.angle(asm_unit_calibrated[k0:, i1, i, j] / asm_default['asm_idl'][k0:, i2, i, j], deg=True))),
np.max(abs(np.angle(asm_unit_calibrated[k0:, i1, i, j] / asm_default['asm_idl'][k0:, i2, i, j], deg=True))))
)
i=0, j=0 => 0.94866 1.06408 0.000 0.000 i=0, j=1 => 0.71232 3.37206 1.830 151.464 i=0, j=2 => 0.91277 1.09417 1.478 21.105 i=0, j=3 => 0.97547 1.03192 0.325 29.291 i=0, j=4 => 0.97549 1.03196 0.326 29.290 i=1, j=1 => 0.57116 2.08833 0.000 0.000 i=1, j=2 => 0.77460 2.32031 0.225 89.019 i=1, j=3 => 0.73100 4.21165 0.307 143.717 i=1, j=4 => 0.73101 4.21141 0.307 143.716 i=2, j=2 => 0.80184 1.21473 0.000 0.000 i=2, j=3 => 0.88331 1.12101 0.234 11.724 i=2, j=4 => 0.88335 1.12100 0.235 11.724 i=3, j=3 => 0.99980 1.00014 0.000 0.000 i=3, j=4 => 0.99983 1.00015 0.000 0.001 i=4, j=4 => 0.99985 1.00018 0.000 0.000
case = 'F1-912Hz-2Vpp_on_E1E2B1B2B3-UnitKCOEFFMatrices'
dir_Packets = dir_DATA + '/' + case + '/decom'
files = glob.glob(dir_Packets + '/*.af1')
files.sort()
print(len(files))
for file in files:
print(os.path.basename(file))
asm_unit = get_lfr_decom_l1_asm(file)
1 2022_03_10_14_58_28_packet_record_NORMAL.af1
case = 'F1-912Hz-2Vpp_on_E1E2B1B2B3-DefaultKCOEFFMatrices'
dir_Packets = dir_DATA + '/' + case + '/decom'
files = glob.glob(dir_Packets + '/*.af1')
files.sort()
print(len(files))
for file in files:
print(os.path.basename(file))
asm_default = get_lfr_decom_l1_asm(file)
1 2022_03_10_14_54_31_packet_record_NORMAL.af1
print('asm mode :', asm_unit['mode'], '@F%d'%(LFR_Fs2F(asm_unit['fe'])))
print('asm time :', asm_unit['relative_time'])
print('asm shape:', asm_unit['asm_idl'].shape)
nspec, nfreq, ndim, ndim = asm_unit['asm_idl'].shape
asm mode : NORMAL @F1 asm time : [ 0. 3.99992371 7.99984741 11.99977112 15.99969482 19.99961853 23.99954224 27.99946594 31.99938965 35.99931335 39.99923706 43.99916077 47.99908447 51.99900818 55.99893188 59.99885559] asm shape: (16, 104, 5, 5)
print('asm mode :', asm_default['mode'], '@F%d'%(LFR_Fs2F(asm_default['fe'])))
print('asm time :', asm_default['relative_time'])
print('asm shape:', asm_default['asm_idl'].shape)
asm mode : NORMAL @F1 asm time : [ 0. 3.99992371 7.99984741 11.99977112 15.99969482 19.99961853 23.99954224 27.99946594 31.99938965 35.99931335 39.9992218 43.99914551 47.99906921 51.99899292 55.99890137 59.99882507 63.99874878] asm shape: (17, 104, 5, 5)
F = 1
calib_mat_freqs_asm_BE = np.zeros((nfreq, ndim, ndim), dtype=complex)
calib_mat_freqs_asm_BE[:, 0:3, 0:3] = B_cal[F][index_freq[F], ...]
calib_mat_freqs_asm_BE[:, 3:5, 3:5] = E_cal[F][index_freq[F], ...]
asm_unit_calibrated = calibrate_sm(asm_unit['asm_idl'], calib_mat_freqs_asm_BE, echo=True)
calibrated sm shape: (16, 104, 5, 5)
i1 = index_from_date(912, asm_unit['freq'], epsi=None)
i2 = index_from_date(912, asm_default['freq'], epsi=None)
i1, asm_unit['freq'][i1], i2, asm_default['freq'][i1]
(51, 912.0, 51, 912.0)
np.array(asm_unit['freq'])
array([ 96., 112., 128., 144., 160., 176., 192., 208., 224.,
240., 256., 272., 288., 304., 320., 336., 352., 368.,
384., 400., 416., 432., 448., 464., 480., 496., 512.,
528., 544., 560., 576., 592., 608., 624., 640., 656.,
672., 688., 704., 720., 736., 752., 768., 784., 800.,
816., 832., 848., 864., 880., 896., 912., 928., 944.,
960., 976., 992., 1008., 1024., 1040., 1056., 1072., 1088.,
1104., 1120., 1136., 1152., 1168., 1184., 1200., 1216., 1232.,
1248., 1264., 1280., 1296., 1312., 1328., 1344., 1360., 1376.,
1392., 1408., 1424., 1440., 1456., 1472., 1488., 1504., 1520.,
1536., 1552., 1568., 1584., 1600., 1616., 1632., 1648., 1664.,
1680., 1696., 1712., 1728., 1744.])
for i in range(5):
for j in range(5):
if i <= j:
print("i=%d, j=%d :"%(i, j))
print(np.abs(asm_unit['asm_idl'][:, i1, i, j]))
print(np.angle(asm_unit['asm_idl'][:, i1, i, j], deg=True))
i=0, j=0 : [3.2812500e-01 2.0768036e+07 2.0764050e+07 2.0764100e+07 2.0762492e+07 2.0760184e+07 2.0761106e+07 2.0761324e+07 2.0762452e+07 2.0762996e+07 2.0764460e+07 2.0764748e+07 2.0767588e+07 2.0767156e+07 2.0769392e+07 2.0771470e+07] [0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.] i=0, j=1 : [0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.] [0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.] i=0, j=2 : [0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.] [0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.] i=0, j=3 : [0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.] [0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.] i=0, j=4 : [0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.] [0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.] i=1, j=1 : [7.8125000e-01 2.0622732e+07 2.0621126e+07 2.0619828e+07 2.0619380e+07 2.0618884e+07 2.0615356e+07 2.0619282e+07 2.0618502e+07 2.0620082e+07 2.0620030e+07 2.0621776e+07 2.0621248e+07 2.0624208e+07 2.0625696e+07 2.0627756e+07] [0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.] i=1, j=2 : [0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.] [0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.] i=1, j=3 : [0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.] [0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.] i=1, j=4 : [0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.] [0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.] i=2, j=2 : [6.0937500e-01 2.0760014e+07 2.0757918e+07 2.0755192e+07 2.0755950e+07 2.0757776e+07 2.0761366e+07 2.0764350e+07 2.0764860e+07 2.0764542e+07 2.0764238e+07 2.0762556e+07 2.0757160e+07 2.0756748e+07 2.0755936e+07 2.0756840e+07] [0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.] i=2, j=3 : [0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.] [0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.] i=2, j=4 : [0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.] [0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.] i=3, j=3 : [8.7500000e-01 2.0512638e+07 2.0508508e+07 2.0507970e+07 2.0508914e+07 2.0510494e+07 2.0516198e+07 2.0518776e+07 2.0519032e+07 2.0518322e+07 2.0516156e+07 2.0511712e+07 2.0510064e+07 2.0509760e+07 2.0510830e+07 2.0515150e+07] [0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.] i=3, j=4 : [0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.] [0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.] i=4, j=4 : [8.9062500e-01 2.0643796e+07 2.0639256e+07 2.0637216e+07 2.0640184e+07 2.0642284e+07 2.0646068e+07 2.0648228e+07 2.0647988e+07 2.0646124e+07 2.0645436e+07 2.0641624e+07 2.0638752e+07 2.0640480e+07 2.0640284e+07 2.0644842e+07] [0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.]
for i in range(5):
for j in range(5):
if i <= j:
print("i=%d, j=%d :"%(i, j))
print(np.abs(asm_default['asm_idl'][:, i1, i, j]))
print(np.angle(asm_default['asm_idl'][:, i1, i, j], deg=True))
i=0, j=0 : [4.5163870e-01 1.4301202e+07 1.4300044e+07 1.4298617e+07 1.4300213e+07 1.4299640e+07 1.4300546e+07 1.4301526e+07 1.4301576e+07 1.4302438e+07 1.4302972e+07 1.4303068e+07 1.4303181e+07 1.4304150e+07 1.4304558e+07 1.4305916e+07 1.4306120e+07] [0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.] i=0, j=1 : [1.0192572e-02 1.3343892e+06 1.3349994e+06 1.3343700e+06 1.3338928e+06 1.3331784e+06 1.3325029e+06 1.3324950e+06 1.3320578e+06 1.3331691e+06 1.3348896e+06 1.3352529e+06 1.3365322e+06 1.3366189e+06 1.3360264e+06 1.3357490e+06 1.3358272e+06] [ 62.304493 168.16464 168.17143 168.16777 168.1605 168.15315 168.14497 168.1439 168.13995 168.15013 168.16815 168.17259 168.18562 168.18593 168.17908 168.1754 168.17474 ] i=0, j=2 : [8.2495503e-02 1.1817081e+06 1.1809496e+06 1.1824476e+06 1.1816476e+06 1.1809582e+06 1.1806554e+06 1.1806751e+06 1.1813706e+06 1.1809331e+06 1.1819279e+06 1.1828646e+06 1.1829014e+06 1.1833712e+06 1.1829486e+06 1.1835934e+06 1.1821974e+06] [-1.404121 -0.91929543 -0.9188002 -0.91981083 -0.9193667 -0.9190967 -0.9189866 -0.9189835 -0.9194488 -0.9190009 -0.919282 -0.919765 -0.919579 -0.9198071 -0.91965127 -0.9200157 -0.91918147] i=0, j=3 : [0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.] [0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.] i=0, j=4 : [0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.] [0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.] i=1, j=1 : [3.8983741e-01 1.3172003e+07 1.3171310e+07 1.3168860e+07 1.3171150e+07 1.3171230e+07 1.3172486e+07 1.3173452e+07 1.3173055e+07 1.3174049e+07 1.3173516e+07 1.3172836e+07 1.3172690e+07 1.3173288e+07 1.3174124e+07 1.3175042e+07 1.3176292e+07] [0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.] i=1, j=2 : [2.4879636e-02 3.6960219e+05 3.7000619e+05 3.6968022e+05 3.6932884e+05 3.6890431e+05 3.6845906e+05 3.6842275e+05 3.6816094e+05 3.6880866e+05 3.6985419e+05 3.7008259e+05 3.7086375e+05 3.7088959e+05 3.7050928e+05 3.7030234e+05 3.7032959e+05] [161.41367 33.17491 33.126328 33.15659 33.203312 33.251648 33.306526 33.31372 33.343025 33.271545 33.151104 33.123554 33.034733 33.03393 33.079266 33.10601 33.10596 ] i=1, j=3 : [0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.] [0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.] i=1, j=4 : [0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.] [0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.] i=2, j=2 : [4.8520991e-01 1.4419251e+07 1.4417640e+07 1.4417179e+07 1.4418298e+07 1.4417400e+07 1.4418154e+07 1.4419114e+07 1.4419602e+07 1.4420093e+07 1.4421049e+07 1.4421659e+07 1.4421686e+07 1.4422890e+07 1.4423087e+07 1.4424798e+07 1.4424179e+07] [0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.] i=2, j=3 : [0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.] [0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.] i=2, j=4 : [0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.] [0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.] i=3, j=3 : [1.0303510e-03 2.6169168e+04 2.6167914e+04 2.6168580e+04 2.6169234e+04 2.6174865e+04 2.6177635e+04 2.6180299e+04 2.6179734e+04 2.6174576e+04 2.6170139e+04 2.6167104e+04 2.6169406e+04 2.6170021e+04 2.6174547e+04 2.6179324e+04 2.6182246e+04] [0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.] i=3, j=4 : [5.1517552e-04 1.3084584e+04 1.3083957e+04 1.3084290e+04 1.3084617e+04 1.3087433e+04 1.3088817e+04 1.3090149e+04 1.3089867e+04 1.3087288e+04 1.3085069e+04 1.3083552e+04 1.3084703e+04 1.3085011e+04 1.3087273e+04 1.3089662e+04 1.3091123e+04] [ 0.0000000e+00 0.0000000e+00 -5.3455636e-07 0.0000000e+00 5.3452942e-07 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 -5.3457296e-07 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 -5.3426379e-07] i=4, j=4 : [1.4464543e-03 3.2542377e+04 3.2540777e+04 3.2540324e+04 3.2543318e+04 3.2549043e+04 3.2554293e+04 3.2555883e+04 3.2556607e+04 3.2551553e+04 3.2545928e+04 3.2541031e+04 3.2542082e+04 3.2548197e+04 3.2551998e+04 3.2556898e+04 3.2559539e+04] [0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.]
for i in range(5):
for j in range(5):
if i <= j:
print("i=%d, j=%d :"%(i, j))
print(np.abs(asm_unit_calibrated[:, i1, i, j]))
print(np.angle(asm_unit_calibrated[:, i1, i, j], deg=True))
i=0, j=0 : [3.98423199e-01 1.43050180e+07 1.43032357e+07 1.43024348e+07 1.43020391e+07 1.43016351e+07 1.43014936e+07 1.43031274e+07 1.43032782e+07 1.43037949e+07 1.43040932e+07 1.43043638e+07 1.43040209e+07 1.43046593e+07 1.43054881e+07 1.43067223e+07] [ 3.99142757e-15 0.00000000e+00 3.73068402e-15 0.00000000e+00 0.00000000e+00 3.73110156e-15 0.00000000e+00 0.00000000e+00 3.73067295e-15 0.00000000e+00 -3.73046038e-15 0.00000000e+00 3.73047924e-15 0.00000000e+00 0.00000000e+00 0.00000000e+00] i=0, j=1 : [4.19949671e-02 1.33472483e+06 1.33406368e+06 1.33466936e+06 1.33405267e+06 1.33300665e+06 1.33246736e+06 1.33188915e+06 1.33209022e+06 1.33231715e+06 1.33279258e+06 1.33324678e+06 1.33521869e+06 1.33519810e+06 1.33600569e+06 1.33639580e+06] [ 10.31078018 168.16339335 168.15877819 168.16571616 168.15980456 168.1494467 168.14234131 168.1355806 168.13701315 168.13928183 168.14366504 168.14877632 168.16941223 168.16961063 168.17726737 168.18023228] i=0, j=2 : [1.46595686e-01 1.18086841e+06 1.18124532e+06 1.18112487e+06 1.18113865e+06 1.18111473e+06 1.17920783e+06 1.18019996e+06 1.17965727e+06 1.18017787e+06 1.17994810e+06 1.18075183e+06 1.18078277e+06 1.18198261e+06 1.18224484e+06 1.18251584e+06] [-1.55730378 -0.91862346 -0.91901226 -0.91886956 -0.91899127 -0.91916098 -0.91813094 -0.9187551 -0.9183977 -0.91865006 -0.91842732 -0.9188174 -0.91852593 -0.9192134 -0.91920841 -0.91926295] i=0, j=3 : [0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.] [0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.] i=0, j=4 : [0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.] [0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.] i=1, j=1 : [2.94494719e-01 1.31763983e+07 1.31744486e+07 1.31736290e+07 1.31733342e+07 1.31731375e+07 1.31745405e+07 1.31755273e+07 1.31760535e+07 1.31761344e+07 1.31765211e+07 1.31760998e+07 1.31753761e+07 1.31751083e+07 1.31755924e+07 1.31765510e+07] [ 0.00000000e+00 0.00000000e+00 -8.10065824e-15 4.05058113e-15 4.05067176e-15 0.00000000e+00 0.00000000e+00 4.04999752e-15 0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00 -4.05012632e-15 0.00000000e+00 4.04968286e-15] i=1, j=2 : [4.78268833e-02 3.69677277e+05 3.69333369e+05 3.69728704e+05 3.69363962e+05 3.68736696e+05 3.68393062e+05 3.67997096e+05 3.68109881e+05 3.68237043e+05 3.68515858e+05 3.68792675e+05 3.70010923e+05 3.69989459e+05 3.70460250e+05 3.70662353e+05] [170.53240081 33.17968581 33.21300538 33.16528947 33.20594341 33.27704469 33.3197317 33.36902082 33.3573595 33.34330071 33.31233562 33.27974601 33.138465 33.1407975 33.08905286 33.06933565] i=1, j=3 : [0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.] [0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.] i=1, j=4 : [0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.] [0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.] i=2, j=2 : [4.73837804e-01 1.44224277e+07 1.44209749e+07 1.44200813e+07 1.44197558e+07 1.44194348e+07 1.44182335e+07 1.44204373e+07 1.44202518e+07 1.44210360e+07 1.44211531e+07 1.44218454e+07 1.44213734e+07 1.44226899e+07 1.44235795e+07 1.44248999e+07] [ 0.00000000e+00 0.00000000e+00 0.00000000e+00 3.70045437e-15 0.00000000e+00 3.70062029e-15 0.00000000e+00 0.00000000e+00 0.00000000e+00 -3.70020939e-15 3.70017935e-15 0.00000000e+00 0.00000000e+00 -3.69978509e-15 0.00000000e+00 0.00000000e+00] i=2, j=3 : [0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.] [0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.] i=2, j=4 : [0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.] [0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.] i=3, j=3 : [1.12942321e-03 2.61788995e+04 2.61731422e+04 2.61705552e+04 2.61743190e+04 2.61769821e+04 2.61817807e+04 2.61845198e+04 2.61842155e+04 2.61818517e+04 2.61809792e+04 2.61761451e+04 2.61725031e+04 2.61746944e+04 2.61744458e+04 2.61802260e+04] [ 0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00 -9.95069108e-16 0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00 9.95431066e-16 0.00000000e+00 -9.95220744e-16] i=3, j=4 : [5.64711606e-04 1.30894498e+04 1.30865711e+04 1.30852776e+04 1.30871595e+04 1.30884910e+04 1.30908903e+04 1.30922599e+04 1.30921077e+04 1.30909258e+04 1.30904896e+04 1.30880726e+04 1.30862515e+04 1.30873472e+04 1.30872229e+04 1.30901130e+04] [ 0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00 -9.95069108e-16 0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00 9.95431066e-16 0.00000000e+00 -9.95220744e-16] i=4, j=4 : [1.39196457e-03 3.25572997e+04 3.25506231e+04 3.25492941e+04 3.25514321e+04 3.25541015e+04 3.25625345e+04 3.25664886e+04 3.25667371e+04 3.25652458e+04 3.25622809e+04 3.25554369e+04 3.25524365e+04 3.25525988e+04 3.25538935e+04 3.25608169e+04] [-5.57846531e-16 -8.00284550e-16 0.00000000e+00 0.00000000e+00 -8.00428806e-16 0.00000000e+00 0.00000000e+00 0.00000000e+00 -8.00052638e-16 0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00]
i = 0
j = 1
print(np.abs(asm_unit_calibrated[:, i1, i, j] / asm_default['asm_idl'][:-1, i2, i, j]))
print(np.angle(asm_unit_calibrated[:, i1, i, j] / asm_default['asm_idl'][:-1, i2, i, j], deg=True))
[4.12015402 1.00025145 0.99929914 1.00022438 1.00011987 0.99987116 0.99997337 0.99954534 1.00002436 0.99936096 0.99842905 0.99849756 0.99901718 0.99893701 0.99998456 1.00048421] [-5.19937144e+01 -1.24250979e-03 -1.26485355e-02 -2.04690602e-03 -6.91030572e-04 -3.70177752e-03 -2.62680739e-03 -8.32804484e-03 -2.94747828e-03 -1.08569439e-02 -2.44897086e-02 -2.38113952e-02 -1.62150233e-02 -1.63172986e-02 -1.81132975e-03 4.83185565e-03]
k0 = 1
k1= -1
for i in range(5):
for j in range(5):
if i <= j:
print("i=%d, j=%d => %8.5f %8.5f %5.3f %5.3f"%(i, j,
np.min(np.abs(asm_unit_calibrated[k0:, i1, i, j] / asm_default['asm_idl'][k0:k1, i2, i, j])),
np.max(np.abs(asm_unit_calibrated[k0:, i1, i, j] / asm_default['asm_idl'][k0:k1, i2, i, j])),
np.min(abs(np.angle(asm_unit_calibrated[k0:, i1, i, j] / asm_default['asm_idl'][k0:k1, i2, i, j], deg=True))),
np.max(abs(np.angle(asm_unit_calibrated[k0:, i1, i, j] / asm_default['asm_idl'][k0:k1, i2, i, j], deg=True))))
)
i=0, j=0 => 1.00004 1.00027 0.000 0.000 i=0, j=1 => 0.99843 1.00048 0.001 0.024 i=0, j=2 => 0.99821 1.00025 0.000 0.001 i=0, j=3 => nan nan nan nan i=0, j=4 => nan nan nan nan i=1, j=1 => 1.00011 1.00036 0.000 0.000 i=1, j=2 => 0.99638 1.00097 0.003 0.161 i=1, j=3 => nan nan nan nan i=1, j=4 => nan nan nan nan i=2, j=2 => 0.99998 1.00023 0.000 0.000 i=2, j=3 => nan nan nan nan i=2, j=4 => nan nan nan nan i=3, j=3 => 1.00000 1.00041 0.000 0.000 i=3, j=4 => 1.00000 1.00041 0.000 0.000 i=4, j=4 => 1.00006 1.00050 0.000 0.000
/tmp/ipykernel_21089/3168270230.py:7: RuntimeWarning: invalid value encountered in true_divide np.min(np.abs(asm_unit_calibrated[k0:, i1, i, j] / asm_default['asm_idl'][k0:k1, i2, i, j])), /tmp/ipykernel_21089/3168270230.py:8: RuntimeWarning: invalid value encountered in true_divide np.max(np.abs(asm_unit_calibrated[k0:, i1, i, j] / asm_default['asm_idl'][k0:k1, i2, i, j])), /tmp/ipykernel_21089/3168270230.py:9: RuntimeWarning: invalid value encountered in true_divide np.min(abs(np.angle(asm_unit_calibrated[k0:, i1, i, j] / asm_default['asm_idl'][k0:k1, i2, i, j], deg=True))), /tmp/ipykernel_21089/3168270230.py:10: RuntimeWarning: invalid value encountered in true_divide np.max(abs(np.angle(asm_unit_calibrated[k0:, i1, i, j] / asm_default['asm_idl'][k0:k1, i2, i, j], deg=True))))
case = 'F2-64Hz-2Vpp_on_E1E2B1B2B3-UnitKCOEFFMatrices'
dir_Packets = dir_DATA + '/' + case + '/decom'
files = glob.glob(dir_Packets + '/*.af2')
files.sort()
print(len(files))
for file in files:
print(os.path.basename(file))
asm_unit = get_lfr_decom_l1_asm(file)
1 2022_03_10_15_05_10_packet_record_NORMAL.af2
case = 'F2-64Hz-2Vpp_on_E1E2B1B2B3-DefaultKCOEFFMatrices'
dir_Packets = dir_DATA + '/' + case + '/decom'
files = glob.glob(dir_Packets + '/*.af2')
files.sort()
print(len(files))
for file in files:
print(os.path.basename(file))
asm_default = get_lfr_decom_l1_asm(file)
1 2022_03_10_15_01_34_packet_record_NORMAL.af2
print('asm mode :', asm_unit['mode'], '@F%d'%(LFR_Fs2F(asm_unit['fe'])))
print('asm time :', asm_unit['relative_time'])
print('asm shape:', asm_unit['asm_idl'].shape)
nspec, nfreq, ndim, ndim = asm_unit['asm_idl'].shape
asm mode : NORMAL @F2 asm time : [ 0. 3.99992371 7.99984741 11.99978638 15.99971008 19.99963379 23.9995575 27.9994812 31.99940491 35.99932861 39.99925232 43.99919128 47.99911499 51.9990387 55.9989624 59.99888611] asm shape: (16, 96, 5, 5)
print('asm mode :', asm_default['mode'], '@F%d'%(LFR_Fs2F(asm_default['fe'])))
print('asm time :', asm_default['relative_time'])
print('asm shape:', asm_default['asm_idl'].shape)
asm mode : NORMAL @F2 asm time : [ 0. 3.99992371 7.99984741 11.99977112 15.99969482 19.99961853 23.99954224 27.99946594 31.99938965 35.99931335 39.99923706 43.99916077 47.99908447 51.99900818 55.99893188 59.99887085 63.99879456 67.99871826 71.99864197 75.99858093] asm shape: (20, 96, 5, 5)
F = 2
calib_mat_freqs_asm_BE = np.zeros((nfreq, ndim, ndim), dtype=complex)
calib_mat_freqs_asm_BE[:, 0:3, 0:3] = B_cal[F][index_freq[F], ...]
calib_mat_freqs_asm_BE[:, 3:5, 3:5] = E_cal[F][index_freq[F], ...]
asm_unit_calibrated = calibrate_sm(asm_unit['asm_idl'], calib_mat_freqs_asm_BE, echo=True)
calibrated sm shape: (16, 96, 5, 5)
i1 = index_from_date(64, asm_unit['freq'], epsi=None)
i2 = index_from_date(64, asm_default['freq'], epsi=None)
i1, asm_unit['freq'][i1], i2, asm_default['freq'][i1]
(57, 64.0, 57, 64.0)
np.array(asm_unit['freq'])
array([ 7., 8., 9., 10., 11., 12., 13., 14., 15., 16., 17.,
18., 19., 20., 21., 22., 23., 24., 25., 26., 27., 28.,
29., 30., 31., 32., 33., 34., 35., 36., 37., 38., 39.,
40., 41., 42., 43., 44., 45., 46., 47., 48., 49., 50.,
51., 52., 53., 54., 55., 56., 57., 58., 59., 60., 61.,
62., 63., 64., 65., 66., 67., 68., 69., 70., 71., 72.,
73., 74., 75., 76., 77., 78., 79., 80., 81., 82., 83.,
84., 85., 86., 87., 88., 89., 90., 91., 92., 93., 94.,
95., 96., 97., 98., 99., 100., 101., 102.])
for i in range(5):
for j in range(5):
if i <= j:
print("i=%d, j=%d :"%(i, j))
print(np.abs(asm_unit['asm_idl'][:, i1, i, j]))
print(np.angle(asm_unit['asm_idl'][:, i1, i, j], deg=True))
i=0, j=0 : [ 0. 19569292. 19567126. 19568844. 19569488. 19569428. 19568360. 19566660. 19571348. 19569760. 19568978. 19569054. 19566608. 19569090. 19567220. 19568286.] [0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.] i=0, j=1 : [ 0. 19449690. 19447414. 19449352. 19449048. 19450874. 19448468. 19448326. 19449620. 19448744. 19448992. 19449074. 19446880. 19449558. 19447298. 19450496.] [ 0. 0.00977211 -0.00378954 0.00543076 -0.00094565 0.00067971 0.0048035 0.00336218 -0.00071805 0.00154738 -0.00336428 0.00324421 0.00269584 0.00208715 0.00567587 0.00223507] i=0, j=2 : [ 0. 19529524. 19528498. 19529714. 19530704. 19529444. 19528430. 19527610. 19530370. 19529268. 19528852. 19528724. 19527758. 19528250. 19526892. 19528632.] [ 0. 0.6479112 1.8407507 3.0333145 4.217862 5.4025974 6.578345 7.752088 8.931496 10.121181 11.30142 12.494687 13.6476145 14.811148 15.981337 17.156614 ] i=0, j=3 : [ 0. 19398270. 19397292. 19397690. 19397544. 19399638. 19397284. 19397256. 19398876. 19398064. 19397000. 19398094. 19395184. 19395760. 19397288. 19396122.] [ 0. 0.7747546 2.3025057 3.857883 5.40566 6.929376 8.437868 9.946199 11.45015 12.957338 14.466389 15.977189 17.471016 18.984085 20.497 22.025557 ] i=0, j=4 : [ 0. 19456838. 19456580. 19456174. 19457690. 19456736. 19456248. 19454642. 19457754. 19455152. 19456510. 19455592. 19456500. 19456742. 19453608. 19454778.] [ 0. 0.7748144 2.3043892 3.8560944 5.401253 6.9308887 8.436542 9.9428215 11.447791 12.958006 14.470352 15.981724 17.471436 18.99007 20.494686 22.022432 ] i=1, j=1 : [ 0. 19330820. 19328432. 19330588. 19329348. 19333036. 19329308. 19330708. 19328646. 19328476. 19329744. 19329828. 19327884. 19330760. 19328110. 19333412.] [0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.] i=1, j=2 : [ 0. 19410172. 19409026. 19410462. 19410506. 19411136. 19408774. 19409502. 19408894. 19408496. 19409122. 19409004. 19408268. 19408966. 19407224. 19411080.] [ 0. 0.63811386 1.844556 3.0278811 4.2188106 5.401916 6.5735593 7.748743 8.9322195 10.119621 11.304823 12.491437 13.644886 14.809051 15.975668 17.154404 ] i=1, j=3 : [ 0. 19279722. 19278622. 19279242. 19278166. 19282116. 19278430. 19279936. 19278220. 19278100. 19278080. 19279180. 19276504. 19277284. 19278416. 19279366.] [ 0. 0.764951 2.3063164 3.8524513 5.406609 6.928696 8.433087 9.942861 11.450873 12.955778 14.469806 15.973936 17.468279 18.981985 20.491335 22.023352 ] i=1, j=4 : [ 0. 19337928. 19337546. 19337370. 19337944. 19338868. 19337032. 19336976. 19336730. 19334838. 19337222. 19336326. 19337444. 19337890. 19334394. 19337670.] [ 0. 0.7650108 2.3082004 3.8506613 5.402202 6.9302077 8.43176 9.939482 11.448515 12.956447 14.473768 15.978472 17.4687 18.987974 20.489017 22.02023 ] i=2, j=2 : [ 0. 19490526. 19490592. 19491328. 19492652. 19490196. 19489202. 19489274. 19490104. 19489526. 19489468. 19489122. 19489602. 19488122. 19487298. 19489700.] [0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.] i=2, j=3 : [ 0. 19359708. 19359838. 19359764. 19359946. 19360836. 19358492. 19359354. 19359062. 19358750. 19358064. 19358936. 19357474. 19356090. 19358138. 19357634.] [0. 0.12686585 0.46176744 0.8245694 1.1878387 1.5267832 1.8595581 2.1941137 2.5186713 2.836186 3.164993 3.482489 3.8233678 4.172932 4.515667 4.868936 ] i=2, j=4 : [ 0. 19418160. 19419026. 19418132. 19419974. 19417822. 19417346. 19416634. 19417822. 19415718. 19417448. 19416318. 19418670. 19416942. 19414346. 19416186.] [0. 0.12692213 0.4636373 0.82278866 1.1834147 1.5282868 1.8582442 2.1907296 2.516293 2.8368587 3.1689372 3.4870286 3.823814 4.1788926 4.51333 4.865806 ] i=3, j=3 : [ 0. 19229808. 19230018. 19229154. 19228204. 19232384. 19228706. 19230344. 19228948. 19228894. 19227596. 19229668. 19226296. 19225006. 19229886. 19226518.] [0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.] i=3, j=4 : [ 0. 19287868. 19288812. 19287128. 19287820. 19288992. 19287168. 19287246. 19287310. 19285478. 19286578. 19286664. 19287076. 19285448. 19285722. 19284676.] [ 0.0000000e+00 5.5698009e-05 1.8654208e-03 -1.7779486e-03 -4.4276314e-03 1.5007883e-03 -1.3108074e-03 -3.3858055e-03 -2.3832012e-03 6.7365810e-04 3.9392267e-03 4.5415279e-03 4.5302912e-04 5.9530078e-03 -2.3440330e-03 -3.1322318e-03] i=4, j=4 : [ 0. 19346102. 19347788. 19345272. 19347620. 19345768. 19345804. 19344316. 19345850. 19342228. 19345740. 19343832. 19348044. 19346076. 19341720. 19343008.] [0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.]
for i in range(5):
for j in range(5):
if i <= j:
print("i=%d, j=%d :"%(i, j))
print(np.abs(asm_default['asm_idl'][:, i1, i, j]))
print(np.angle(asm_default['asm_idl'][:, i1, i, j], deg=True))
i=0, j=0 : [0.0000000e+00 1.7585013e+08 1.7596922e+08 1.7603837e+08 1.7620131e+08 1.7636691e+08 1.7653709e+08 1.7669237e+08 1.7689288e+08 1.7709459e+08 1.7725837e+08 1.7747774e+08 1.7772405e+08 1.7794077e+08 1.7820194e+08 1.7839707e+08 1.7862851e+08 1.7894547e+08 1.7922597e+08 1.7952790e+08] [0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.] i=0, j=1 : [0.0000000e+00 5.8858452e+07 5.9817912e+07 6.1123452e+07 6.2773008e+07 6.4768452e+07 6.7000608e+07 6.9469832e+07 7.2207248e+07 7.5161976e+07 7.8253392e+07 8.1547536e+07 8.4988120e+07 8.8516456e+07 9.2144120e+07 9.5799384e+07 9.9605528e+07 1.0358768e+08 1.0758340e+08 1.1157708e+08] [ 0. 6.918542 10.802507 14.537177 18.041357 21.3602 24.364206 27.105206 29.608511 31.883043 33.95451 35.778942 37.431404 38.881058 40.159786 41.27864 42.27243 43.18143 43.95068 44.625343] i=0, j=2 : [0.0000000e+00 4.8256467e+08 4.8258906e+08 4.8249498e+08 4.8249168e+08 4.8241754e+08 4.8230582e+08 4.8212022e+08 4.8198870e+08 4.8178778e+08 4.8148906e+08 4.8121005e+08 4.8092928e+08 4.8056803e+08 4.8021382e+08 4.7972090e+08 4.7925184e+08 4.7879040e+08 4.7829165e+08 4.7773594e+08] [ 0. 176.71169 175.75185 174.79433 173.83789 172.86798 171.91197 170.96484 170.02899 169.06853 168.11949 167.17078 166.21928 165.28162 164.34167 163.40955 162.47229 161.51752 160.57056 159.63202] i=0, j=3 : [ 0. 2147293.8 2148111.5 2148609.5 2149546.5 2150427.8 2151623. 2152312. 2153704.8 2154977.2 2155919.8 2157238.5 2158639.8 2160056.5 2161501. 2162954.2 2164165.2 2166180.5 2167880.5 2169797.8] [ 0. 89.807396 90.03364 90.27439 90.51939 90.730484 90.96594 91.215836 91.45412 91.6791 91.88999 92.10243 92.31521 92.52901 92.74457 92.991295 93.25113 93.48178 93.71765 93.957405] i=0, j=4 : [ 0. 3214468.8 3215421.2 3216235.2 3217623.5 3219096.8 3220782.2 3222060.5 3223879.5 3225820.8 3227237.5 3229264. 3231350. 3233360.2 3235862.5 3237735.8 3239711.5 3242727. 3245289.8 3247904.8] [ 0. 89.806755 90.034096 90.275665 90.523125 90.73613 90.96801 91.21566 91.4592 91.677956 91.89419 92.103035 92.31417 92.53261 92.7457 92.99004 93.24879 93.48449 93.716774 93.955376] i=1, j=1 : [ 0. 19707562. 20341828. 21230518. 22370932. 23792976. 25436260. 27320524. 29482004. 31907602. 34553256. 37477276. 40648756. 44039668. 47652708. 51451104. 55548920. 59972012. 64585972. 69352240.] [0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.] i=1, j=2 : [0.0000000e+00 1.6150456e+08 1.6403376e+08 1.6751626e+08 1.7187754e+08 1.7714837e+08 1.8303549e+08 1.8954250e+08 1.9673547e+08 2.0446768e+08 2.1255066e+08 2.2109659e+08 2.2997211e+08 2.3904912e+08 2.4829909e+08 2.5760290e+08 2.6722899e+08 2.7715411e+08 2.8709584e+08 2.9690768e+08] [ 0. 169.7925 164.94824 160.25552 155.7947 151.5057 147.54538 143.85735 140.41794 137.1828 134.16234 131.38913 128.78526 126.39781 124.1793 122.12849 120.19719 118.3336 116.617424 115.00435 ] i=1, j=3 : [ 0. 718731.7 730229.8 746048. 765805.1 789728.7 816613.2 846232.44 879148.06 914619. 951772.6 991215.5 1032274.25 1074521.2 1117667. 1161512.8 1206768.8 1253958.1 1301309.4 1348536.2 ] [ 0. 82.88917 79.23158 75.737816 72.47873 69.37098 66.6025 64.111565 61.846436 59.796757 57.936226 56.32432 54.88469 53.64855 52.58558 51.71355 50.97953 50.3012 49.767647 49.332775] i=1, j=4 : [ 0. 1075933.1 1093050.6 1116752. 1146322.6 1182189.6 1222393.9 1266829.1 1315996.2 1369107.1 1424726.5 1483794.4 1545250.8 1608437.9 1673196.8 1738674. 1806509.8 1877149.1 1948044.5 2018581.8] [ 0. 82.88856 79.23204 75.73898 72.482414 69.37665 66.60456 64.11134 61.851498 59.795624 57.940414 56.324898 54.88363 53.652214 52.58666 51.712273 50.977215 50.3039 49.76676 49.330753] i=2, j=2 : [0.0000000e+00 1.3242714e+09 1.3235127e+09 1.3224759e+09 1.3212342e+09 1.3195877e+09 1.3177047e+09 1.3155343e+09 1.3133267e+09 1.3107379e+09 1.3079022e+09 1.3047726e+09 1.3014440e+09 1.2979082e+09 1.2940946e+09 1.2900264e+09 1.2858387e+09 1.2810889e+09 1.2764206e+09 1.2713152e+09] [0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.] i=2, j=3 : [ 0. 5892531. 5891089. 5888985. 5886064.5 5882047. 5878278.5 5872734.5 5868272.5 5862611.5 5856120.5 5849073.5 5841345. 5833684. 5824723. 5816284.5 5806322. 5795846.5 5785287. 5773948.5] [ 0. -86.90433 -85.71828 -84.51998 -83.318565 -82.13755 -80.94606 -79.74911 -78.57494 -77.389465 -76.229546 -75.068436 -73.90419 -72.75269 -71.597206 -70.41837 -69.221306 -68.03587 -66.85298 -65.67474 ] i=2, j=4 : [ 0. 8821037. 8818135. 8815172. 8810762. 8805169. 8799246. 8791621. 8784214. 8775841. 8766139. 8755735. 8744133. 8732360. 8719870. 8706423. 8691943. 8676262. 8660502. 8642846.] [ 0. -86.90497 -85.71781 -84.518684 -83.31481 -82.131905 -80.944 -79.749275 -78.569855 -77.39062 -76.22535 -75.06783 -73.90522 -72.749084 -71.59606 -70.41962 -69.223625 -68.03315 -66.853874 -65.67679 ] i=3, j=3 : [6.7785347e-04 2.6220486e+04 2.6222699e+04 2.6224576e+04 2.6223156e+04 2.6220016e+04 2.6223881e+04 2.6217625e+04 2.6221801e+04 2.6222900e+04 2.6221582e+04 2.6221240e+04 2.6218930e+04 2.6221355e+04 2.6217961e+04 2.6224525e+04 2.6219891e+04 2.6222203e+04 2.6222264e+04 2.6224504e+04] [0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.] i=3, j=4 : [3.3892679e-04 3.9251715e+04 3.9251703e+04 3.9255332e+04 3.9253047e+04 3.9250227e+04 3.9254742e+04 3.9248391e+04 3.9251402e+04 3.9253484e+04 3.9251586e+04 3.9251715e+04 3.9248105e+04 3.9250398e+04 3.9249445e+04 3.9255613e+04 3.9250633e+04 3.9254090e+04 3.9254402e+04 3.9254656e+04] [ 0. -0.00063425 0.0004625 0.00126345 0.00373699 0.00564609 0.00207326 -0.00018061 0.005082 -0.00113831 0.00419478 0.00060307 -0.00103658 0.00360574 0.00112859 -0.00125553 -0.00234021 0.00271245 -0.0008904 -0.00202932] i=4, j=4 : [1.6946340e-04 5.8759289e+04 5.8754277e+04 5.8760953e+04 5.8757301e+04 5.8755891e+04 5.8760742e+04 5.8755750e+04 5.8755402e+04 5.8759184e+04 5.8756453e+04 5.8757602e+04 5.8751969e+04 5.8753402e+04 5.8758148e+04 5.8761898e+04 5.8757391e+04 5.8762555e+04 5.8763352e+04 5.8759102e+04] [0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.]
for i in range(5):
for j in range(5):
if i <= j:
print("i=%d, j=%d :"%(i, j))
print(np.abs(asm_unit_calibrated[:, i1, i, j]))
print(np.angle(asm_unit_calibrated[:, i1, i, j], deg=True))
i=0, j=0 : [0.00000000e+00 1.75901690e+08 1.76022456e+08 1.76139086e+08 1.76302506e+08 1.76381645e+08 1.76540207e+08 1.76670734e+08 1.76908628e+08 1.77064671e+08 1.77238120e+08 1.77419126e+08 1.77619424e+08 1.77809584e+08 1.78019617e+08 1.78239926e+08] [ 0.00000000e+00 -3.03356112e-16 1.51573993e-16 3.02947257e-16 -3.02666445e-16 1.51265323e-16 -1.20903569e-15 -1.51017805e-15 1.20651782e-15 -1.50681817e-15 -2.55908406e-15 -1.95495013e-15 1.50211198e-15 -9.00303319e-16 5.09569963e-15 4.04158332e-15] i=0, j=1 : [ 0. 58888689.38286138 59788124.16481612 61002554.19880284 62532651.80364019 64313711.22706834 66367184.13273928 68654706.89681692 71165944.32808751 73893728.94954431 76774324.13150269 79838689.99646686 82947994.46373634 86174850.96068308 89545001.45750031 93033680.24013746] [ 0. 6.93898824 10.61524846 14.17092094 17.50388883 20.64680526 23.549324 26.2141078 28.6843912 30.94784275 32.97727862 34.84531395 36.45564185 37.93435258 39.25875177 40.4356241 ] i=0, j=2 : [0.00000000e+00 4.82693528e+08 4.82703845e+08 4.82706688e+08 4.82709006e+08 4.82566962e+08 4.82451936e+08 4.82287794e+08 4.82226055e+08 4.81994445e+08 4.81758441e+08 4.81481522e+08 4.81187864e+08 4.80868898e+08 4.80490647e+08 4.80144594e+08] [ 0. 176.72301254 175.79133842 174.89402586 173.98037323 173.07724789 172.18670063 171.28809773 170.38442238 169.48013438 168.57227436 167.67229796 166.79296437 165.90470302 165.01889541 164.11626438] i=0, j=3 : [ 0. 2147902.0057427 2148732.03434523 2149303.93888553 2150431.20025093 2150811.24695486 2151779.93302813 2152492.18036445 2154026.63809148 2154774.76567744 2156025.63467687 2157019.16570584 2158471.84090442 2159516.7679075 2160548.4551151 2161957.04057941] [ 0. 89.8205619 90.07140954 90.36836637 90.65226486 90.92781313 91.19350956 91.45533234 91.71069191 91.96516423 92.22478924 92.48258858 92.7549085 93.04525634 93.32022216 93.60435427] i=0, j=4 : [ 0. 3215386.94269125 3216548.6783564 3217495.5539551 3218999.89972711 3219904.92402757 3221147.88052506 3222388.3214954 3224521.72255308 3225840.08830453 3227444.47863004 3229153.99477712 3230905.63976018 3232506.963528 3234567.58343702 3236415.79081548] [ 0. 89.82052395 90.0701643 90.36955135 90.65521239 90.92681194 91.19438345 91.45758603 91.71227646 91.96471631 92.22216304 92.47956654 92.75461128 93.0412893 93.32177841 93.60644017] i=1, j=1 : [ 0. 19722001.97452998 20314453.8065901 21133979.15636528 22186472.86953313 23457361.72481325 24955891.20866026 26685913.8261834 28634696.02509005 30844544.29868944 33263272.78071916 35934105.7623807 38742851.53059809 41770740.63968249 45048407.13247438 48566177.0178309 ] [ 0.00000000e+00 -2.02923819e-15 6.56685794e-16 -7.57465300e-15 -1.20255376e-15 1.19427104e-14 -2.67275833e-15 1.29973268e-14 -6.98813768e-15 4.32498308e-16 -8.02098658e-16 1.78195678e-14 -2.30699149e-14 7.66481782e-15 -2.45788424e-14 -2.47213033e-14] i=1, j=2 : [0.00000000e+00 1.61583376e+08 1.63943452e+08 1.67163771e+08 1.71199557e+08 1.75945637e+08 1.81358189e+08 1.87407376e+08 1.93977273e+08 2.01138321e+08 2.08674338e+08 2.16657945e+08 2.24705309e+08 2.33043429e+08 2.41682042e+08 2.50607763e+08] [ 0. 169.78347813 165.17500744 160.7217592 156.47486246 152.4286187 148.63530977 145.0717375 141.69778673 138.52988072 135.59246532 132.82465223 130.33509876 127.96800678 125.75773073 123.67818471] i=1, j=3 : [ 0. 719094.48169072 729862.15640691 744389.26055765 762751.48858677 784261.40975605 808935.12760713 836476.30476164 866523.19804955 899250.95992186 933937.89813298 970669.99854002 1008011.27608186 1046611.32607033 1086778.13102053 1128457.14480295] [ 0. 82.88202842 79.45662455 76.19807094 73.14911023 70.28173628 67.6449689 65.24194955 63.02711727 61.0180997 59.24832078 57.63817622 56.30019646 55.11173114 54.06233025 53.16957189] i=1, j=4 : [ 0. 1076476.98352462 1092566.93511222 1114346.9490173 1141769.48040806 1174090.0443137 1210950.52257079 1252246.07891935 1297161.83009213 1346238.58986716 1398050.4388984 1453136.66355259 1508841.32719872 1566636.43331837 1627019.86862508 1689281.92539916] [ 0. 82.88199751 79.45538449 76.19924963 73.1520933 70.28074107 67.64580304 65.24419458 63.02871526 61.01766697 59.24573384 57.63515089 56.29988317 55.10780043 54.06390009 53.17162345] i=2, j=2 : [0.00000000e+00 1.32459045e+09 1.32373660e+09 1.32287691e+09 1.32166273e+09 1.32029196e+09 1.31847821e+09 1.31660869e+09 1.31449986e+09 1.31208246e+09 1.30951298e+09 1.30667198e+09 1.30360752e+09 1.30048914e+09 1.29691094e+09 1.29344379e+09] [ 0.00000000e+00 3.22278348e-16 0.00000000e+00 0.00000000e+00 0.00000000e+00 -3.23327595e-16 -3.23772376e-16 -3.24232117e-16 0.00000000e+00 0.00000000e+00 3.25988998e-16 -3.26697772e-16 0.00000000e+00 6.56501942e-16 -3.29156621e-16 0.00000000e+00] i=2, j=3 : [ 0. 5894045.94053436 5892393.86758366 5890099.00083486 5887756.58135258 5884421.74594072 5880386.96782844 5875987.99150448 5871516.04096591 5865568.25039971 5860350.41184423 5853702.95272265 5847469.1567777 5840172.66296005 5831476.10251911 5823866.03807478] [ 0. -86.90252364 -85.7199475 -84.52571068 -83.3281678 -82.14948017 -80.99326332 -79.832804 -78.67378484 -77.5150698 -76.34750984 -75.18977231 -74.03816243 -72.85951728 -71.69873185 -70.51196815] i=2, j=4 : [ 0. 8823327.58632295 8820634.64713217 8817443.04237343 8813434.97321468 8809364.03956107 8802758.86069069 8796648.67972678 8789507.34498179 8781140.94378513 8772602.17013844 8763254.61847974 8752776.14853051 8741955.42115925 8730331.21550231 8718239.06470973] [ 0. -86.90256268 -85.72119888 -84.52452072 -83.32522973 -82.15048533 -80.99238175 -79.83055283 -78.6722098 -77.51551635 -76.35014578 -75.19279335 -74.03844808 -72.86349968 -71.69718585 -70.50988521] i=3, j=3 : [ 0. 26227.64957744 26229.93529976 26226.52434047 26229.70754095 26227.19677124 26227.2455767 26225.2282844 26227.30793923 26222.39756779 26227.15881144 26224.57212212 26230.2823608 26227.61432905 26221.70886853 26223.45501939] [ 0.00000000e+00 0.00000000e+00 0.00000000e+00 -1.24182982e-16 0.00000000e+00 0.00000000e+00 0.00000000e+00 1.24189119e-16 0.00000000e+00 1.24202525e-16 0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00] i=3, j=4 : [ 0. 39262.52611367 39264.94877164 39260.96028024 39263.49007347 39263.82353142 39261.37512282 39260.47224768 39261.59882501 39256.65996986 39260.53190954 39259.35516925 39262.76878579 39259.22777205 39256.64642615 39256.10144331] [ 0.00000000e+00 -3.70946746e-05 -1.24234764e-03 1.18410924e-03 2.94870694e-03 -9.99544801e-04 8.72986934e-04 2.25497958e-03 1.58719496e-03 -4.48665299e-04 -2.62346864e-03 -3.02469675e-03 -3.01701489e-04 -3.96451728e-03 1.56117770e-03 2.08605125e-03] i=4, j=4 : [ 0. 58775.6026848 58777.73860148 58773.43151974 58773.8775474 58780.50559986 58773.05870026 58774.88077073 58773.59488246 58769.81042584 58770.73230672 58773.01125051 58770.4259169 58765.80295538 58771.31390511 58765.76635129] [ 0.00000000e+00 0.00000000e+00 1.10820459e-16 1.10828581e-16 0.00000000e+00 0.00000000e+00 0.00000000e+00 -1.10825848e-16 -1.10828272e-16 0.00000000e+00 -1.10833671e-16 1.10829373e-16 0.00000000e+00 0.00000000e+00 0.00000000e+00 1.10843037e-16]
i = 0
j = 3
k1= -4
print(np.abs(asm_unit_calibrated[:, i1, i, j] / asm_default['asm_idl'][:k1, i2, i, j]))
print(np.angle(asm_unit_calibrated[:, i1, i, j] / asm_default['asm_idl'][:k1, i2, i, j], deg=True))
[ nan 1.00028332 1.00028882 1.00032326 1.00041154 1.00017833 1.00007293 1.00008368 1.00014943 0.99990601 1.00004907 0.99989836 0.99992218 0.99975008 0.99955936 0.99953892] [ nan 0.01317233 0.03777593 0.09397279 0.13287955 0.19732646 0.22756568 0.23949089 0.25657409 0.28606967 0.33479749 0.38015874 0.43970305 0.51625516 0.57565067 0.61305688]
/tmp/ipykernel_21089/450759236.py:4: RuntimeWarning: invalid value encountered in true_divide print(np.abs(asm_unit_calibrated[:, i1, i, j] / asm_default['asm_idl'][:k1, i2, i, j])) /tmp/ipykernel_21089/450759236.py:5: RuntimeWarning: invalid value encountered in true_divide print(np.angle(asm_unit_calibrated[:, i1, i, j] / asm_default['asm_idl'][:k1, i2, i, j], deg=True))
k0 = 1
k1= -4
for i in range(5):
for j in range(5):
if i <= j:
print("i=%d, j=%d => %8.5f %8.5f %5.3f %5.3f"%(i, j,
np.min(np.abs(asm_unit_calibrated[k0:, i1, i, j] / asm_default['asm_idl'][k0:k1, i2, i, j])),
np.max(np.abs(asm_unit_calibrated[k0:, i1, i, j] / asm_default['asm_idl'][k0:k1, i2, i, j])),
np.min(abs(np.angle(asm_unit_calibrated[k0:, i1, i, j] / asm_default['asm_idl'][k0:k1, i2, i, j], deg=True))),
np.max(abs(np.angle(asm_unit_calibrated[k0:, i1, i, j] / asm_default['asm_idl'][k0:k1, i2, i, j], deg=True))))
)
i=0, j=0 => 0.99898 1.00057 0.000 0.000 i=0, j=1 => 0.97113 1.00051 0.020 0.977 i=0, j=2 => 1.00024 1.00088 0.011 0.707 i=0, j=3 => 0.99954 1.00041 0.013 0.613 i=0, j=4 => 0.99959 1.00043 0.014 0.616 i=1, j=1 => 0.94393 1.00073 0.000 0.000 i=1, j=2 => 0.97285 1.00049 0.009 1.578 i=1, j=3 => 0.97154 1.00050 0.007 1.477 i=1, j=4 => 0.97159 1.00051 0.007 1.477 i=2, j=2 => 1.00017 1.00265 0.000 0.000 i=2, j=3 => 1.00019 1.00130 0.002 0.134 i=2, j=4 => 1.00026 1.00136 0.002 0.133 i=3, j=3 => 0.99996 1.00043 0.000 0.000 i=3, j=4 => 1.00001 1.00037 0.000 0.008 i=4, j=4 => 1.00007 1.00042 0.000 0.000