def VU_vfrdB(a, rng):
N_len = int(a.length)
ca = a.otherEmpty.fill(0.0)
fftType = {'cvview_d': 'ccfftip_d', 'cvview_f': 'ccfftip_f'}[ca.type]
fft = pv.FFT(fftType, N_len, 1, pv.VSIP_FFT_FWD, 0, 0)
ra = ca.realview
ia = ca.imagview
ta = a.cloneview
s = ta.stride
pv.copy(a, ra)
fft.dft(ca)
pv.cmagsq(ca, ra)
maxv = ra.maxval
minv = maxv * rng
pv.clip(ra, minv, maxv, minv, maxv, ra)
Nlen = int(N_len / 2)
if N_len % 2:
ta.putlength(Nlen + 1)
ra.putlength(Nlen + 1)
ta.putoffset(Nlen * s)
pv.copy(ra, ta)
ra.putlength(Nlen)
ta.putlength(Nlen)
ta.putoffset(a.offset)
ra.putoffset(Nlen + 1)
pv.copy(ra, ta)
else:
pv.copy(ra, ta)
ta.putlength(Nlen)
a.putlength(Nlen)
ta.putoffset(Nlen * s)
pv.swap(ta, a)
a.putlength(N_len)
a.log10
a *= 10.0
def VU_vfrdB(a,rng):
N_len=int(a.length)
ca=a.otherEmpty.fill(0.0)
fftType={'cvview_d':'ccfftip_d','cvview_f':'ccfftip_f'}[ca.type]
fft = pv.FFT(fftType,N_len,1,pv.VSIP_FFT_FWD,0,0)
ra = ca.realview
ia = ca.imagview
ta = a.cloneview
s = ta.stride
pv.copy(a,ra)
fft.dft(ca)
pv.cmagsq(ca,ra)
maxv = ra.maxval
minv = maxv * rng
pv.clip(ra,minv,maxv,minv,maxv,ra)
Nlen = int(N_len/2)
if N_len%2:
ta.putlength(Nlen+1)
ra.putlength(Nlen+1)
ta.putoffset(Nlen * s)
pv.copy(ra,ta)
ra.putlength(Nlen)
ta.putlength(Nlen)
ta.putoffset(a.offset)
ra.putoffset(Nlen+1)
pv.copy(ra,ta)
else :
pv.copy(ra,ta)
ta.putlength(Nlen)
a.putlength(Nlen)
ta.putoffset(Nlen * s)
pv.swap(ta,a)
a.putlength(N_len)
a.log10
a *= 10.0
import pyJvsip as pjv
from math import pi as pi
n = 1000
x = pjv.create('vview_d', n).ramp(0, 2 * pi / float(n))
c = pjv.cos(x, x.empty)
cClip = pjv.clip(c, -.9, .9, -.8, .8, c.empty)
import pyJvsip as pjv
from math import pi as pi
n=1000
x=pjv.create('vview_d',n).ramp(0,2*pi/float(n))
c=pjv.cos(x,x.empty)
cClip=pjv.clip(c,-.9,.9,-.8,.8,c.empty)
