为什么预计计算速度的单位都是GFlops/sec ?

lch=gta · 发表于 2011-8-2 03:39:16

除了数学项目之外，基本也就没有用整数的了吧……

引用一段Ap的源代码

CPUComputation::CPUComputation ()
14 {
15 m_bInterleaved = false;
16 if ( m_bInterleaved )
17 {
18 m_pfInput =
19 (float *)fftwf_malloc ( 4 * sciFFTLength * sizeof ( float ) );
20 m_pfOutput =
21 (float *)fftwf_malloc ( 4 * sciFFTLength * sizeof ( float ) );
22 m_pfInputInv =
23 (float *)fftwf_malloc ( 4 * sciFFTLength * sizeof ( float ) );
24 m_pfOutputInv =
25 (float *)fftwf_malloc ( 4 * sciFFTLength * sizeof ( float ) );
26
27 m_fpForwardPlan = fftwf_plan_many_dft ( 1, &sciFFTLength, 2,
28 (fftwf_complex *)m_pfInput,
29 NULL, 2, 1,
30 (fftwf_complex *)m_pfOutput,
31 NULL, 2, 1,
32 FFTW_FORWARD, FFTW_MEASURE );
33 m_fpInversePlan = fftwf_plan_many_dft ( 1, &sciFFTLength, 2,
34 (fftwf_complex *)m_pfInputInv,
35 NULL, 2, 1,
36 (fftwf_complex *)m_pfOutputInv,
37 NULL, 2, 1,
38 FFTW_BACKWARD, FFTW_MEASURE );
39 }
40 else
41 {
42 m_pfInput =
43 (float *)fftwf_malloc ( 2 * sciFFTLength * sizeof ( float ) );
44 m_pfOutput =
45 (float *)fftwf_malloc ( 2 * sciFFTLength * sizeof ( float ) );
46 m_pfInputInv =
47 (float *)fftwf_malloc ( 2 * sciFFTLength * sizeof ( float ) );
48 m_pfOutputInv =
49 (float *)fftwf_malloc ( 2 * sciFFTLength * sizeof ( float ) );
50
51 m_fpForwardPlan = fftwf_plan_dft_1d ( sciFFTLength,
52 (fftwf_complex *)m_pfInput,
53 (fftwf_complex *)m_pfOutput,
54 FFTW_FORWARD, FFTW_MEASURE );
55 m_fpInversePlan = fftwf_plan_dft_1d ( sciFFTLength,
56 (fftwf_complex *)m_pfInputInv,
57 (fftwf_complex *)m_pfOutputInv,
58 FFTW_BACKWARD, FFTW_MEASURE );
59 }
60 }
61
62 CPUComputation::~CPUComputation ()
63 {
64 fftwf_destroy_plan ( m_fpForwardPlan );
65 fftwf_destroy_plan ( m_fpInversePlan );
66
67 fftwf_free ( m_pfInput );
68 fftwf_free ( m_pfOutput );
69 fftwf_free ( m_pfInputInv );
70 fftwf_free ( m_pfOutputInv );
71 }
72
73 float CPUComputation::Execute ( float *pfInput, float *pfOutput )
74 {
75
76 int i;
77
78 if ( m_bInterleaved )
79 {
80 for ( i = 0; i < sciFFTQuantity; i++ )
81 {
82 memcpy ( m_pfInput, pfInput + 4 * i * sciFFTLength,
83 4 * sciFFTLength * sizeof ( float ) );
84
85 fftwf_execute ( m_fpForwardPlan );
86
87 memcpy ( m_pfInputInv, m_pfOutput,
88 4 * sciFFTLength * sizeof ( float ) );
89
90 fftwf_execute ( m_fpInversePlan );
91
92 memcpy ( pfOutput + 4 * i * sciFFTLength, m_pfOutputInv,
93 4 * sciFFTLength * sizeof ( float ) );
94 }
95 }
96 else
97 {
98 for ( i = 0; i < 2 * sciFFTQuantity; i++ )
99 {
100 memcpy ( m_pfInput, pfInput + 2 * i * sciFFTLength,
101 2 * sciFFTLength * sizeof ( float ) );
102
103 fftwf_execute ( m_fpForwardPlan );
104 /*
105 memcpy ( m_pfInputInv, m_pfOutput,
106 2 * sciFFTLength * sizeof ( float ) );
107
108 fftwf_execute ( m_fpInversePlan );
109
110 memcpy ( pfOutput + 2 * i * sciFFTLength, m_pfOutputInv,
111 2 * sciFFTLength * sizeof ( float ) );
112 */
113
114 memcpy ( pfOutput + 2 * i * sciFFTLength, m_pfOutput,
115 2 * sciFFTLength * sizeof ( float ) );
116
117 }
118 }

复制代码

除了for循环之外就全是float了……

nekoko · 发表于 2011-8-2 04:53:26

那就麻烦您给解释一下这段东西的意义~~

nekoko · 发表于 2011-8-2 04:58:14

本帖最后由 nekoko 于 2011-8-2 04:59 编辑

回复 15# lch=gta

我不举别的例子，SmithWaterman就是一个纯整数的算法……

争论这个没有意义，我只是认为您所谓的“科学计算基本上都是浮点”的概念有偏颇之嫌

而且实际上根据AMD的调查，在科学计算服务器上90%左右的操作都是整数操作，这也是他们做出Bulldozer架构加强整数线程的原因之一

lch=gta · 发表于 2011-8-2 10:34:20

那就麻烦您给解释一下这段东西的意义~~
nekoko 发表于 2011-8-2 04:53

这段主要是对内存的存取，仅仅是代码中的一段，没有上下文具体他要干嘛是很难理解的，但是可以看出它操作的变量都是浮点

可能也是我说的不严谨

但是对于科学计算，肯定还是浮点要多一些，当然了比较算法，排序算法这些还是要涉及到整数

但是做数值运算的话，像seti einstein这种，都要用浮点的了
不然为什么大名鼎鼎的ScienceMark都是在跑浮点

90%整数操作这个悬了，推土机加强整数不能说明问题，那现在世界上越来越多的超级计算机GPU加速，是否可以说明90%都是浮点操作？

lch=gta · 发表于 2011-8-2 10:38:37

而且BOINC官方给出的help也说了，credit就是根据浮点操作给的
如果绝大多数项目都是整数，那为什么不用整数操作呢？

jpegreader · 发表于 2011-8-2 19:17:03

讨论科学计算是浮点多还是整数多是很重要的哦！

		自动登录	找回密码
密码			新注册用户

为什么预计计算速度的单位都是GFlops/sec ?

浏览过的版块