In [11]:

    

%load_ext autoreload
%autoreload 2
%matplotlib inline

import sys
sys.path.append("..")
import splitwavepy as sw

import numpy as np
import matplotlib.pyplot as plt


    

    




The autoreload extension is already loaded. To reload it, use:
  %reload_ext autoreload

In [36]:

    

a = sw.Trio(split=(30,1.2),ray=(30,45),delta=0.05)
a.plot()
print(a.rayvecs)
print(a.cmpvecs)


    

    













    




[[-0.61237244  0.5         0.61237244]
 [-0.35355339 -0.8660254   0.35355339]
 [ 0.70710678  0.          0.70710678]]
[[ 1.  0.  0.]
 [ 0.  1.  0.]
 [ 0.  0.  1.]]

In [45]:

    

a.rotate2ray()
a.set_window(15,35,tukey=0.1)
a.plot()
print(a.rayvecs)
print(a.cmpvecs)


    

    













    




[[  1.00000000e+00  -3.08148791e-33   1.11022302e-16]
 [ -3.08148791e-33   1.00000000e+00  -1.11022302e-16]
 [  1.11022302e-16  -1.11022302e-16   1.00000000e+00]]
[[ -6.12372436e-01  -3.53553391e-01   7.07106781e-01]
 [  5.00000000e-01  -8.66025404e-01  -3.92523115e-17]
 [  6.12372436e-01   3.53553391e-01   7.07106781e-01]]

In [46]:

    

b = sw.Eig3dM(a)


    

In [47]:

    

b.plot()


    

In [61]:

    

plt.contourf(b.lags,b.degs,b.lam1/b.lam2)


    

    
Out[61]:





<matplotlib.contour.QuadContourSet at 0x108cf28d0>






    

In [52]:

    

np.std(b.lam3)


    

    
Out[52]:





2.3347220871022585e-09

In [ ]: