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
