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

图1

% {; I/ o* ~: M6 y; U左边栏中X-CUBE-ALGOBUILD，选择，并勾选DSP library。生成项目。
然后使用VSCODE打开项目，编写FFT实现函数。主要代码如下所示：
% P+ a  h+ V2 R- R9 A& h% v8 w

1.   arm_rfft_fast_init_f32(&fft_handler, FFT_SIZE);
   
2.   GenerateSineWave();
   
3.   uint32_t start = DWT->CYCCNT;
   ! Y2 H2 J. I3 @. g$ ~. |/ U* L. ^; |- r
4.   arm_rfft_fast_f32(&fft_handler, input, fft_output, 0);
   6 B6 u  k: h7 R* `. x
5.   arm_cmplx_mag_f32(fft_output, magnitude, FFT_SIZE / 2);
   7 t( l2 {$ ?( E! _0 P9 P
6.   uint32_t end = DWT->CYCCNT;
   
7.   volatile uint32_t cycles = end - start;
   5 _. x+ N; K9 H) H+ ~# z- j1 \  n( }" ~
8.   printf("FFT computation took %lu cycles\r\n", cycles);

复制代码

$ A0 ?( D) D5 a在gcc-arm-none-eabi.cmake 中有
2 b; g+ U6 S( s. e( A
" T+ H, k7 [/ F; _
+ W# Q% b# o) W/ d4 L2 xset(TARGET_FLAGS "-mcpu=cortex-m33 -mfpu=fpv5-sp-d16 -mfloat-abi=hard ")
0 y( C8 A# c0 y0 g; ^

表示编译器已经按硬件FPU方式生成浮点指令。
编译，Debug,观察左边Variables 中时间变量circles 的时间数值，包括观察FFT后峰值是否在bin 10的位置。测试结果如图2所示：

图2
, m9 S6 h' R, Z) P. @, e$ C从测试结果看和预期相同。
& L& o8 x; Q0 ?, D! S