为了和后续的HSP加速运算进行对比，先测试一下U3C5的FPU进行简单的正弦函数FFT计算的状况。
还是先从Cubemx这个工具开始。有两种方式使用DSP library。一种是下载库，然后移植到项目里，一种是在Cubemx里配置DSP library。如图1所示：

. [0 Z$ M9 o5 ~3 J图1

左边栏中X-CUBE-ALGOBUILD，选择，并勾选DSP library。生成项目。
3 \& n% b) f% H然后使用VSCODE打开项目，编写FFT实现函数。主要代码如下所示：
: B4 [4 e/ Y" N7 u* R

1.   arm_rfft_fast_init_f32(&fft_handler, FFT_SIZE);
   9 e( N8 c9 X+ V
2.   GenerateSineWave();
   
3.   uint32_t start = DWT->CYCCNT;
   $ n  }' N' a+ S3 I! Y  M* {( {
4.   arm_rfft_fast_f32(&fft_handler, input, fft_output, 0);
   # G, r" Z/ m2 M6 c) X3 Q
5.   arm_cmplx_mag_f32(fft_output, magnitude, FFT_SIZE / 2);
   
6.   uint32_t end = DWT->CYCCNT;
   
7.   volatile uint32_t cycles = end - start;
   - k# G$ H0 \' i/ x
8.   printf("FFT computation took %lu cycles\r\n", cycles);

复制代码

在gcc-arm-none-eabi.cmake 中有
6 u# p6 ]% S  a; R' H& w
, i1 h5 c. L% T# {$ }
set(TARGET_FLAGS "-mcpu=cortex-m33 -mfpu=fpv5-sp-d16 -mfloat-abi=hard ")
& h# u) D% V+ S

7 v7 n$ L) j: g表示编译器已经按硬件FPU方式生成浮点指令。
编译，Debug,观察左边Variables 中时间变量circles 的时间数值，包括观察FFT后峰值是否在bin 10的位置。测试结果如图2所示：
6 N+ T: j. n5 S8 v, \5 N
图2
从测试结果看和预期相同。

2 @4 `; Y7 P$ I3 {& r; Y