为了和后续的HSP加速运算进行对比，先测试一下U3C5的FPU进行简单的正弦函数FFT计算的状况。
2 ~9 n1 J1 V1 L3 y- [还是先从Cubemx这个工具开始。有两种方式使用DSP library。一种是下载库，然后移植到项目里，一种是在Cubemx里配置DSP library。如图1所示：
$ Z2 ^- ]: T( X
图1
  E( b3 D- v! }7 W
! q/ O1 S: J4 n, i& w左边栏中X-CUBE-ALGOBUILD，选择，并勾选DSP library。生成项目。
9 a9 V+ B# ^9 r5 j/ A1 V, m然后使用VSCODE打开项目，编写FFT实现函数。主要代码如下所示：
3 X. f6 _6 |) t1 B9 X" o

1.   arm_rfft_fast_init_f32(&fft_handler, FFT_SIZE);
   
2.   GenerateSineWave();
   
3.   uint32_t start = DWT->CYCCNT;
   
4.   arm_rfft_fast_f32(&fft_handler, input, fft_output, 0);
   ' I: e( ?8 y* T) R
5.   arm_cmplx_mag_f32(fft_output, magnitude, FFT_SIZE / 2);
   7 g: ^, K( R+ }+ H- S
6.   uint32_t end = DWT->CYCCNT;
   
7.   volatile uint32_t cycles = end - start;
   
8.   printf("FFT computation took %lu cycles\r\n", cycles);

复制代码

在gcc-arm-none-eabi.cmake 中有


9 m$ k2 n# _" N6 `/ ^set(TARGET_FLAGS "-mcpu=cortex-m33 -mfpu=fpv5-sp-d16 -mfloat-abi=hard ")
: I0 i( n+ _9 Y

, n2 A% u0 J. [3 S表示编译器已经按硬件FPU方式生成浮点指令。
8 S& r* x9 U" ^- K7 ^$ [0 |编译，Debug,观察左边Variables 中时间变量circles 的时间数值，包括观察FFT后峰值是否在bin 10的位置。测试结果如图2所示：

图2
从测试结果看和预期相同。
3 f' _4 D) \7 D
; R: H3 ~6 y& Z) F4 V9 N