In [57]:

    

from __future__ import division, print_function

# Third-party
import astropy.units as u
import matplotlib.pyplot as plt
import numpy as np
%matplotlib inline

# Custom
import gary.dynamics as gd
import gary.integrate as gi
import gary.io as io
import gary.potential as gp
from gary.units import galactic


    

In [74]:

    

# w0 = np.load("/Users/adrian/projects/morphology/output/freqmap/triaxial-NFW/E-0.135_tube_grid_xz/w0.npy")
# lyap = np.memmap("/Users/adrian/projects/morphology/output/freqmap/triaxial-NFW/E-0.135_tube_grid_xz/alllyap.dat", 
#                  shape=(len(w0),2), mode='r', dtype=float)

w0 = np.load("/Users/adrian/projects/morphology/output/freqmap/triaxial-NFW-yz/E-0.135_tube_grid_xz/w0.npy")
lyap = np.memmap("/Users/adrian/projects/morphology/output/freqmap/triaxial-NFW-yz/E-0.135_tube_grid_xz/alllyap.dat", 
                 shape=(len(w0),2), mode='r', dtype=float)


    

In [75]:

    

donezo = lyap[:,1] == 1.


    

In [76]:

    

lyap[donezo,0].max()


    

    
Out[76]:





memmap(0.00027137784360404356)

In [79]:

    

plt.figure(figsize=(8,8))
plt.scatter(w0[donezo,0], w0[donezo,2], c=np.log10(lyap[donezo,0]), marker='s', s=5, cmap='Greys_r')
plt.scatter(w0[~donezo,0], w0[~donezo,2], c='k', marker='s', s=5)
plt.xlim(-1,50)
plt.ylim(-1,50)


    

    
Out[79]:





(-1, 50)






    

In [65]:

    




    

In [ ]: