| Christophe Favergeon | 0c9cd47 | 2021-07-20 10:51:46 +0200 | [diff] [blame] | 1 | import cmsisdsp as dsp |
| 2 | import numpy as np |
| Christophe Favergeon | e7de243 | 2021-09-07 13:28:35 +0200 | [diff] [blame] | 3 | import cmsisdsp.fixedpoint as f |
| Christophe Favergeon | 0c9cd47 | 2021-07-20 10:51:46 +0200 | [diff] [blame] | 4 | |
| 5 | # Test vlog q31 and q15 |
| 6 | x = np.array([0.9,0.5,2**-16]) |
| 7 | |
| 8 | r=dsp.arm_vlog_q15(f.toQ15(x)) |
| Christophe Favergeon | 0c9cd47 | 2021-07-20 10:51:46 +0200 | [diff] [blame] | 9 | print(f.Q15toF32(r)*16.0) |
| 10 | |
| Christophe Favergeon | 19dc952 | 2021-09-22 15:37:22 +0200 | [diff] [blame^] | 11 | r=dsp.arm_vlog_q31(f.toQ31(x)) |
| 12 | print(f.Q31toF32(r)*32.0) |
| 13 | |
| Christophe Favergeon | 0c9cd47 | 2021-07-20 10:51:46 +0200 | [diff] [blame] | 14 | print(np.log(x)) |
| 15 | |
| 16 | print("") |
| 17 | # Test sin_cos |
| 18 | t=20 |
| 19 | |
| 20 | sinRef=np.sin(t * np.pi / 180) |
| 21 | cosRef=np.cos(t * np.pi / 180) |
| 22 | print(sinRef) |
| 23 | print(cosRef) |
| 24 | |
| 25 | s,c=dsp.arm_sin_cos_f32(t) |
| 26 | print(s) |
| 27 | print(c) |
| 28 | |
| 29 | s,c=dsp.arm_sin_cos_q31(f.toQ31(t/180.0)) |
| 30 | print(f.Q31toF32(s)) |
| 31 | print(f.Q31toF32(c)) |
| 32 | |
| 33 | print("") |
| 34 | # Test sqrt |
| 35 | a=0.6 |
| 36 | print(np.sqrt(a)) |
| 37 | |
| 38 | err,r=dsp.arm_sqrt_f32(a) |
| 39 | print(err,r) |
| 40 | |
| 41 | err,r=dsp.arm_sqrt_q31(f.toQ31(a)) |
| 42 | print(err,f.Q31toF32(r)) |
| 43 | |
| 44 | err,r=dsp.arm_sqrt_q15(f.toQ15(a)) |
| 45 | print(err,f.Q15toF32(r)) |
| 46 | |
| 47 | err,r=dsp.arm_sqrt_f32(-a) |
| 48 | print(err,r) |
| 49 | |
| 50 | err,r=dsp.arm_sqrt_q31(f.toQ31(-a)) |
| 51 | print(err,f.Q31toF32(r)) |
| 52 | |
| 53 | err,r=dsp.arm_sqrt_q15(f.toQ15(-a)) |
| 54 | print(err,f.Q15toF32(r)) |
| 55 | |