i using code spectral estimation of geophysical time series.
import matplotlib.pyplot plt plt.style.use("ggplot") import numpy np mtspec import mtspec mtspec.util import _load_mtdata guppy import hpy h = hpy() print h.heap() data = np.loadtxt('500.dat') spec,freq,jackknife,f_statistics,degrees_of_f = mtspec(data=data, delta= 40.0, time_bandwidth=4 ,number_of_tapers=7, nfft= 512, statistics=true) i got this:
milenko@milenko-hp-compaq-6830s:~$ python mt4.py partition of set of 214843 objects. total size = 35303944 bytes. index count % size % cumulative % kind (class / dict of class) 0 15461 7 8538896 24 8538896 24 unicode 1 77241 36 7769768 22 16308664 46 str 2 43700 20 3870600 11 20179264 57 tuple 3 582 0 1443984 4 21623248 61 dict of pyqt4.qtcore.pyqtwrappertype 4 352 0 1395712 4 23018960 65 dict of module 5 10118 5 1295104 4 24314064 69 types.codetype 6 858 0 1265136 4 25579200 72 dict (no owner) 7 9893 5 1187160 3 26766360 76 function 8 1123 1 1112392 3 27878752 79 dict of type 9 1123 1 1012552 3 28891304 82 type <900 more rows. type e.g. '_.more' view.> if take time series 1000points(500 in previous case),i got same:
milenko@milenko-hp-compaq-6830s:~$ python mt5.py partition of set of 214843 objects. total size = 35303944 bytes. index count % size % cumulative % kind (class / dict of class) 0 15461 7 8538896 24 8538896 24 unicode 1 77241 36 7769768 22 16308664 46 str 2 43700 20 3870600 11 20179264 57 tuple 3 582 0 1443984 4 21623248 61 dict of pyqt4.qtcore.pyqtwrappertype 4 352 0 1395712 4 23018960 65 dict of module 5 10118 5 1295104 4 24314064 69 types.codetype 6 858 0 1265136 4 25579200 72 dict (no owner) 7 9893 5 1187160 3 26766360 76 function 8 1123 1 1112392 3 27878752 79 dict of type 9 1123 1 1012552 3 28891304 82 type <900 more rows. type e.g. '_.more' view.> so means doubling time series nothing(that can see) changes.how interpret this?is guppy right choice me?
Comments
Post a Comment