def dpss0Fast(N, W): import dpss sines = numpy.arange(N) Ntapers = 0 #with this code all tapers can be generated very fast v, sig, totit, ifault = dpss.dpss(N, Ntapers, W, sines, sines.copy(), sines.copy(), sines.copy()) return v
def dpss0Fast(N, W): import dpss sines = numpy.arange(N) Ntapers = 0 # with this code all tapers can be generated very fast v, sig, totit, ifault = dpss.dpss(N, Ntapers, W, sines, sines.copy(), sines.copy(), sines.copy()) return v
def dpssFast(N,W,K): import dpss sines = numpy.zeros(N) taperOrder = K v,sig,totit,ifault = dpss.dpss(N,taperOrder,W,sines,sines.copy(),sines.copy(),sines.copy()) sig += 1. return v, sig
def dpssFast(N, W, K): import dpss sines = numpy.zeros(N) taperOrder = K v, sig, totit, ifault = dpss.dpss(N, taperOrder, W, sines, sines.copy(), sines.copy(), sines.copy()) sig += 1. return v, sig
VERSION=1 SUBSAMPLE=1 NWORKERS=cpu_count() FS=int(FS/SUBSAMPLE); NFFT=int(nextpow2(NFFT*FS)) # convert to ticks NWINDOWS_PER_WORKER=int(12*256*1000/NFFT) # NFFT/2 ticks FIRST_MT=float('nan') LAST_MT=float('nan') FRACTION_MT=float('nan') tapers,eig = dpss(NFFT, NW, K) tapers = np.array(tapers, dtype=np.float32) #tapers = tapers * np.sqrt(FS) f=np.array(range(0,NFFT//2+1))*FS/NFFT df=f[1]-f[0]; DIROUT=os.path.dirname(FILEIN); FILEINs=sorted(glob.glob(FILEIN+'.ch*')); FILE_TYPE=1 if (len(FILEINs)==0): FILEINs=sorted(glob.glob(FILEIN+'*.wav')); FILE_TYPE=2 if (len(FILEINs)==0): print(["can't find any .wav or .ch files with basename '"+FILEIN]); exit