CMSIS-DSP/PythonWrapper/testdsp3.py

import cmsisdsp as dsp
import numpy as np
import cmsisdsp.fixedpoint as f

# Test vlog q31 and q15
x = np.array([0.9,0.5,2**-16])

r=dsp.arm_vlog_q15(f.toQ15(x))

print(f.Q15toF32(r)*16.0)

print(np.log(x))

print("")
# Test sin_cos 
t=20

sinRef=np.sin(t * np.pi / 180)
cosRef=np.cos(t * np.pi / 180)
print(sinRef)
print(cosRef)

s,c=dsp.arm_sin_cos_f32(t)
print(s)
print(c)

s,c=dsp.arm_sin_cos_q31(f.toQ31(t/180.0))
print(f.Q31toF32(s))
print(f.Q31toF32(c))

print("")
# Test sqrt
a=0.6
print(np.sqrt(a))

err,r=dsp.arm_sqrt_f32(a)
print(err,r)

err,r=dsp.arm_sqrt_q31(f.toQ31(a))
print(err,f.Q31toF32(r))

err,r=dsp.arm_sqrt_q15(f.toQ15(a))
print(err,f.Q15toF32(r))

err,r=dsp.arm_sqrt_f32(-a)
print(err,r)

err,r=dsp.arm_sqrt_q31(f.toQ31(-a))
print(err,f.Q31toF32(r))

err,r=dsp.arm_sqrt_q15(f.toQ15(-a))
print(err,f.Q15toF32(r))