| Christophe Favergeon | 936a5ea | 2022-05-25 10:39:00 +0200 | [diff] [blame] | 1 | import cmsisdsp as dsp |
| 2 | import numpy as np |
| 3 | |
| 4 | recipQ15=np.array([0x7F03, 0x7D13, 0x7B31, 0x795E, 0x7798, 0x75E0, |
| 5 | 0x7434, 0x7294, 0x70FF, 0x6F76, 0x6DF6, 0x6C82, |
| 6 | 0x6B16, 0x69B5, 0x685C, 0x670C, 0x65C4, 0x6484, |
| 7 | 0x634C, 0x621C, 0x60F3, 0x5FD0, 0x5EB5, 0x5DA0, |
| 8 | 0x5C91, 0x5B88, 0x5A85, 0x5988, 0x5890, 0x579E, |
| 9 | 0x56B0, 0x55C8, 0x54E4, 0x5405, 0x532B, 0x5255, |
| 10 | 0x5183, 0x50B6, 0x4FEC, 0x4F26, 0x4E64, 0x4DA6, |
| 11 | 0x4CEC, 0x4C34, 0x4B81, 0x4AD0, 0x4A23, 0x4978, |
| 12 | 0x48D1, 0x482D, 0x478C, 0x46ED, 0x4651, 0x45B8, |
| 13 | 0x4521, 0x448D, 0x43FC, 0x436C, 0x42DF, 0x4255, |
| 14 | 0x41CC, 0x4146, 0x40C2, 0x4040]) |
| 15 | |
| 16 | s,v=dsp.arm_recip_q15(int(0x2000),recipQ15) |
| 17 | print("1 / 0.25") |
| 18 | print(s) |
| 19 | print("%04X -> %f" % (v,((v<<s)/(1<<15)))) |
| 20 | print("----\n") |
| 21 | |
| 22 | s,v=dsp.arm_recip_q15(1,recipQ15) |
| 23 | print("1 / (1 << 15)") |
| 24 | print(s) |
| 25 | print("%04X -> %f" % (v,((v<<s)/(1<<15)))) |
| 26 | print(1<<15) |
| 27 | print("----\n") |
| 28 | |