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

, M( {7 \" m4 ^0 L0 W左边栏中X-CUBE-ALGOBUILD，选择，并勾选DSP library。生成项目。
然后使用VSCODE打开项目，编写FFT实现函数。主要代码如下所示：
5 _, C& m1 g8 t; p7 K

1.   arm_rfft_fast_init_f32(&fft_handler, FFT_SIZE);
   + R! x) ^2 C/ U1 i6 C7 }, F
2.   GenerateSineWave();
   $ b- {0 m$ m1 X/ i7 v+ F4 ~' z( F. L7 N, F
3.   uint32_t start = DWT->CYCCNT;
   
4.   arm_rfft_fast_f32(&fft_handler, input, fft_output, 0);
   9 x6 t) }" s, M) e: p
5.   arm_cmplx_mag_f32(fft_output, magnitude, FFT_SIZE / 2);
   9 M3 b' [/ g$ K, J
6.   uint32_t end = DWT->CYCCNT;
   
7.   volatile uint32_t cycles = end - start;
   7 A: |' {1 ]2 e9 E. H: _% p
8.   printf("FFT computation took %lu cycles\r\n", cycles);

复制代码

* @) f3 [" X2 l1 M2 o2 u. r在gcc-arm-none-eabi.cmake 中有
% A* z8 C4 [+ ~3 ~
  o8 ~/ u/ e& j/ `# Z8 [- a2 [
set(TARGET_FLAGS "-mcpu=cortex-m33 -mfpu=fpv5-sp-d16 -mfloat-abi=hard ")


; W! s3 Q" ~* o  b0 K' t% }表示编译器已经按硬件FPU方式生成浮点指令。
  n$ e3 K$ b, P" k/ G/ s6 i( x编译，Debug,观察左边Variables 中时间变量circles 的时间数值，包括观察FFT后峰值是否在bin 10的位置。测试结果如图2所示：

图2
从测试结果看和预期相同。