POPPUR爱换

标题: INTEL：新处理器快过CELL 5000倍，下代游戏主机无须手柄 [打印本页]
作者: CC9K    时间: 2007-12-17 10:55
标题: INTEL：新处理器快过CELL 5000倍，下代游戏主机无须手柄
著名芯片生产商Intel公司日前向索尼、微软和任天堂三家家用机生产商，就下一代芯片将如何进一步拓展家用机性能的问题进行了说明。Intel相信，新一代家用机霸主Wii的继任者将不再使用传统的游戏手柄，而是完全依靠身体动作感应进行游戏操作。    “我们相信在不久的将来，任天堂Wii主机的新一代后继者将不再使用手柄进行操作。”Intel公司首席技术总监Justin Rattner对媒体表示：“你只需要将房间中的摄像头设置好，然后挥挥手就可以模拟网球等游戏操作。”

    根据美国商业周刊的报道，Intel公司对于家用机芯片的业务一直十分看好。虽然在与IBM的次世代主机芯片竞赛中落败，但Intel显然没有放弃下一代家用机芯片业务中潜藏的巨大价值并已经开始着手进行准备。

    尽管IBM公司为索尼PS3主机提供的CELL处理器，运算能力可以达到每秒2亿次以上，但Rattner表示，Intel的下一代微型处理器将可以为家用游戏主机提供每秒1万亿次的运算能力，性能将远远超过IBM的CELL处理器。
作者: jiuzhege    时间: 2007-12-17 11:19
你只需要将房间中的摄像头设置好，然后挥挥手就可以模拟网球等游戏操作。

eyetoy 么。。。


主机怎样 不是单纯的技术领先就行了
作者: librae    时间: 2007-12-17 11:32
还没出来就开始吹。。
作者: 紫色    时间: 2007-12-17 11:37
intel必须吹。你见过街上卖东西的小贩不吆喝的么？

你必须告诉别人，你必须吸引别人!

哪怕忽悠大众、吸引别人来骂你都行啊

最怕的就是自己成了粪坑里的石头-又臭又硬 ，没人理！
作者: EviL    时间: 2007-12-17 11:50
干脆说以后的游戏机就是一跟线头，直接插进脑子里
作者: HerculesVR    时间: 2007-12-17 12:05
支持任魔出品Hgame···:loveliness:
作者: itany    时间: 2007-12-17 12:14
20G flops变成2亿次了？作者真是挨驴踢了……
Intel只是鼓吹一下Larrabee的浮点能力而以，1T flops

没准下一代家用机会采用Larrabee+Gesher的异构多核心呢
作者: boris_lee    时间: 2007-12-17 12:25
国内的翻译，呵呵
作者: bwdinlife    时间: 2007-12-17 13:54
原帖由 紫色 于 2007-12-17 11:37 发表
intel必须吹。你见过街上卖东西的小贩不吆喝的么？

你必须告诉别人，你必须吸引别人!

哪怕忽悠大众、吸引别人来骂你都行啊

最怕的就是自己成了粪坑里的石头-又臭又硬 ，没人理！



AMD必须吹。你见过街上卖东西的小贩不吆喝的么？

你必须告诉别人，你必须吸引别人!

哪怕忽悠大众、吸引别人来骂你都行啊

最怕的就是自己成了粪坑里的石头-又臭又硬 ，没人理！
作者: airpotato    时间: 2007-12-17 17:05
原帖由 紫色 于 2007-12-17 11:37 发表
intel必须吹。你见过街上卖东西的小贩不吆喝的么？

你必须告诉别人，你必须吸引别人!

哪怕忽悠大众、吸引别人来骂你都行啊

最怕的就是自己成了粪坑里的石头-又臭又硬 ，没人理！

i是小贩，a是什么？乞 丐？:w00t):
作者: Asuka    时间: 2007-12-17 17:14
原帖由 紫色 于 2007-12-17 11:37 发表
intel必须吹。你见过街上卖东西的小贩不吆喝的么？

你必须告诉别人，你必须吸引别人!

哪怕忽悠大众、吸引别人来骂你都行啊

最怕的就是自己成了粪坑里的石头-又臭又硬 ，没人理！


看AFAN现在一个个粪世嫉俗的，真是可怜啊 :w00t):
作者: the_god_of_pig    时间: 2007-12-17 18:06
原帖由 紫色 于 2007-12-17 11:37 发表
intel必须吹。你见过街上卖东西的小贩不吆喝的么？

你必须告诉别人，你必须吸引别人!

哪怕忽悠大众、吸引别人来骂你都行啊

最怕的就是自己成了粪坑里的石头-又臭又硬 ，没人理！




哪怕忽悠大众、吸引别人来骂你都行啊

最怕的就是自己成了粪坑里的石头-又臭又硬 ，没人理！


您这话说得太对了:p
作者: 北欧候鸟    时间: 2007-12-17 19:25
不知道啥时候才能真的出来
作者: bwdinlife    时间: 2007-12-17 20:17
原帖由 EviL 于 2007-12-17 11:50 发表
干脆说以后的游戏机就是一跟线头，直接插进脑子里


貌似ps2已经说过了
作者: bessel    时间: 2007-12-17 20:29
大哥，您从哪里转来的这篇报道？
原帖由 CC9K 于 2007-12-17 10:55 发表
著名芯片生产商Intel公司日前向索尼、微软和任天堂三家家用机生产商，就下一代芯片将如何进一步拓展家用机性能的问题进行了说明。Intel相信，新一代家用机霸主Wii的继任者将不再使用传统的游戏手柄，而是完全依靠身体 ...
作者: IAMPG    时间: 2007-12-17 20:38
完了，几大护法又把火烧到AMD身上去了,:(
作者: yzsxlijunjun    时间: 2007-12-17 20:42
S饭不少么。
i是该吹了。
当年s吹ee浮点是i的不知道几倍，结果s的是32b精度，i的是80b精度。
作者: OCbaichi    时间: 2007-12-17 20:43
“你只需要将房间中的摄像头设置好，然后挥挥手就可以模拟网球等游戏操作。”
不知道有没有18x游戏?
作者: 路路步兵    时间: 2007-12-17 20:47
提示: 作者被禁止或删除 内容自动屏蔽
作者: D65    时间: 2007-12-17 21:02
AMD吹牛就算了，还编些假事件，假数据，泡制假新闻，弄出一付Intel也爱吹牛的模样。俗称抹黑战术。我黑了，也要搞黑对手。

和民进党非常像，实在是又臭又硬。
作者: naze    时间: 2007-12-17 22:33
原帖由 D65 于 2007-12-17 21:02 发表
AMD吹牛就算了，还编些假事件，假数据，泡制假新闻，弄出一付Intel也爱吹牛的模样。俗称抹黑战术。我黑了，也要搞黑对手。

和民进党非常像，实在是又臭又硬。


同求新闻来源链接...
不会是某某小编梦中呓语得到的吧
作者: Jason21    时间: 2007-12-18 08:03
原帖由 紫色 于 2007-12-17 11:37 发表
intel必须吹。你见过街上卖东西的小贩不吆喝的么？

你必须告诉别人，你必须吸引别人!

哪怕忽悠大众、吸引别人来骂你都行啊

最怕的就是自己成了粪坑里的石头-又臭又硬 ，没人理！

越看越像AMD:whistling:
作者: jiuzhege    时间: 2007-12-18 09:11
原帖由 bwdinlife 于 2007-12-17 20:17 发表


貌似ps2已经说过了


ps2 说脑后插管？  搞笑
作者: 紫色    时间: 2007-12-18 09:15
有些东西是必然的，小贩总要叫卖、鸡总是要给旅客打电话。
商人必须叫卖自己的产品，必须产品有了点眉目就放出风声去吸引人，因为让公众作出反应、去猜测、去研究从而做好购买的准备都需要时间，这样产品一出来就可以顺利地被大家接受、购买。
所以有的时候就会发生跳票事件：nb吹过头了，公众的胃口是吊起来了，自己却又没有能耐实现~~~所以对产品上市前的风声应半信半疑，不要激动，只有见到了实物、见到了独立可信网站的实测才认可。
作者: 紫色    时间: 2007-12-18 09:22
继续。
还有一点值得注意，厂商总是要吹嘘自家的产品的，因为没人嫌钱多。结论就是：自家的产品即使不好也必须吹、落后了也要吹。ia都是如此，谴责他们没用。
重要的是自己要有独立思考、独立判断的能力，不要盲目相信别人，这一点对一个人来说实在是太重要了。
作者: 紫色    时间: 2007-12-18 09:29
说教继续。
什么是事实？事实是不以你我意志为转移的“独立“的东西吗？错了，事实是人的意志的一种状态。我在看叔本华“作为意志和表象的世界”这本书之前就悟出了这一点。:a)
当我判断哪种U好的时候，那就是事实（对我）:p。
作者: itany    时间: 2007-12-18 10:47
原帖由 紫色 于 2007-12-18 09:29 发表
说教继续。
什么是事实？事实是不以你我意志为转移的“独立“的东西吗？错了，事实是人的意志的一种状态。我在看叔本华“作为意志和表象的世界”这本书之前就悟出了这一点。:a)
当我判断哪种U好的时候，那就是事实 ...


你可以认为别人都是纱布，就像别人都认为你是纱布一样
你可以认为疯人院外边的人都是疯子，就像你被关进疯人院一样
你可以认为AMD很好，就像认为那一陀很香一样……
作者: shike_cuke    时间: 2007-12-18 10:49
等实际产品出来以后才知道分晓啊!
作者: 紫色    时间: 2007-12-19 10:09
嗬嗬，无所谓。
ibm为cell体系优化了其编程语言，在c,fortran中引入vector数据类型很大程度上就是为了提供针对对cell spu的vector intrinsics，也算是踏踏实实干点实事吧。
即使intel实现诺言，其新处理器有强大的峰值性能，仅提供编程工具也至少要向后推迟1-2年，除非他现在就动手。
ifan们欢呼是不是早了点呢？
http://publib.boulder.ibm.com/in ... /vectordatatype.htm
作者: potomac    时间: 2007-12-19 11:36
提示: 作者被禁止或删除 内容自动屏蔽
作者: itany    时间: 2007-12-19 11:38
原帖由 紫色 于 2007-12-19 10:09 发表
嗬嗬，无所谓。
ibm为cell体系优化了其编程语言，在c,fortran中引入vector数据类型很大程度上就是为了提供针对对cell spu的vector intrinsics，也算是踏踏实实干点实事吧。
即使intel实现诺言，其新处理器有强大的 ...


前几天还讨论过Larrabee的体系结构和CT……
忘了？
作者: Prescott    时间: 2007-12-19 11:45
原帖由 紫色 于 2007-12-19 10:09 发表
嗬嗬，无所谓。
ibm为cell体系优化了其编程语言，在c,fortran中引入vector数据类型很大程度上就是为了提供针对对cell spu的vector intrinsics，也算是踏踏实实干点实事吧。
即使intel实现诺言，其新处理器有强大的 ...

Cell？你不是很鄙视SIMD吗？Cell还能活多久都是个问题，那个双精度加强版Cell去哪里了？

至于Intel编译器和MS C++编译器的SIMD intrinsic，你总不能选择性请注意文明用词吧。Intel的我就懒得说了，直接帖MSVC的吧。
在Intel编译器和MS编译器中都已经用了几十年的东西，你也好意思替IBM拿出来吹？忘了说了，Intel的intrinsic和MSVC的完全兼容，不过Intel的出来的更早些，明白怎么回事了吗？
不懂就虚心学，别什么都不懂就知道整天到处攻击。

http://msdn2.microsoft.com/en-us/library/y0dh78ez.aspx

Streaming SIMD Extensions 2 Instructions

This section describes the C/C++ language-level features supporting the Streaming SIMD Extensions 2 (SSE2) instructions:  
Other topics discussed in this section include:  
Theemmintrin.h header file contains the declarations for the SSE2instructions intrinsics. The file dvec.h contains operator overloadsfor some of the SSE2 instructions intrinsics, which are available foruse in C++ programs.
SSE2 intrinsics use the __m128, __m128i, and __m128d data types, which are not supported on Itanium Processor Family (IPF) processors.  Any SSE2 intrinsics that use the __m64 data type are not supported on x64 processors.

Floating-Point Intrinsics Using Streaming SIMD Extensions 2 Instructions

The following topics list the floating-point and integer intrinsics brokeninto groups by the nature of the operation. Each intrinsic entry has aninformal pseudo-code, and it is followed with a correspondinginstruction name in uppercase letters; for example, ADDSD is the name of the first instruction listed in this section. The variable ris generally used for the intrinsic's return value. A number appendedto a variable name indicates the element of a packed object. Forexample, r0 is the lowest double of r. Some intrinsics are compositesbecause they require more than one instruction to implement them. Formore details, refer to the Streaming SIMD Extensions 2 (SSE2)instructions external architecture specification (EAS). You should befamiliar with the hardware features provided by the SSE2 instructionswhen writing programs with the intrinsics. The following are threeimportant issues to keep in mind:
This section contains the following topics:  

后面的自己看吧

[ 本帖最后由 Prescott 于 2007-12-19 11:50 编辑 ]
作者: 紫色    时间: 2007-12-19 12:13
谁说我鄙视simd了，任意翻我的帖子你找不到我犯了这条罪。我鄙视的是那些盲目YY理论性能的，cell和sse都是如此，我很希望能体验到simd的好处。intel编译器带simd intrinsic我知道，一个例子是intel确实用sse2实现了sin函数（libimf.so库里的sin_sse2.o），不过请P大研究研究它的速度，好像慢的很，比fsin还慢。在我看来真不如象sqrt那样，设置sinpd,sinps,sinsd,sinss之类的sse指令对我更实惠些。
老话题不要再讨论了。跑题了已经，我上面的帖子说的是“万亿处理器”编程工具落后于cell，落后了几年那是事实。

[ 本帖最后由 紫色 于 2007-12-19 12:53 编辑 ]
作者: Prescott    时间: 2007-12-19 13:46
原帖由 紫色 于 2007-12-19 12:13 发表
谁说我鄙视simd了，任意翻我的帖子你找不到我犯了这条罪。我鄙视的是那些盲目YY理论性能的，cell和sse都是如此，有了编程工具让我体验到好处，我为啥鄙视它。intel编译器带simd intrinsic我知道，我琢磨了libimf.so库 ...

谈万亿次？Cell只是个骗Sony的玩具而已。

难道IBM的intrinsic就是为万亿次作准备，Intel/MS的intrinsic就什么都不是？你也太有色眼睛了吧。难道就因为Intel/MS的intrinsic出来的太早，那个时候还没有万亿次的概念？

sin sse2实现比fsin慢？您用的是AMD的处理器吧。:whistling:
作者: 紫色    时间: 2007-12-19 15:15
也许吧。intel的simd内置函数出来的时候有谁联想到了万亿次？
用sse实现sin就涉及到算法，还要考虑到那些sse指令的延迟/吞吐，弄得不好就会慢了点，也不算奇怪。除了少数几个函数以外，大多数操作都是sse会比x87快点。

sin_sse.o我反汇编了给各位看看，P大认为很简单嘛？：
sin_sse2.o:     file format elf32-i386

Disassembly of section .text:

00000000 <__libm_sse2_sin>:
   0:   55                      push   %ebp
   1:   8b ec                   mov    %esp,%ebp
   3:   83 ec 68                sub    $0x68,%esp
   6:   66 0f c5 c0 03          pextrw $0x3,%xmm0,%eax
   b:   25 ff 7f 00 00          and    $0x7fff,%eax
  10:   2d 30 30 00 00          sub    $0x3030,%eax
  15:   3d c5 10 00 00          cmp    $0x10c5,%eax
  1a:   0f 87 37 01 00 00       ja     157 <__libm_sse2_sin+0x157>
  20:   f2 0f 10 0d 70 08 00    movsd  0x870,%xmm1
  27:   00
  28:   f2 0f 59 c8             mulsd  %xmm0,%xmm1
  2c:   f2 0f 10 15 80 08 00    movsd  0x880,%xmm2
  33:   00
  34:   f2 0f 2d d1             cvtsd2si %xmm1,%edx
  38:   f2 0f 58 ca             addsd  %xmm2,%xmm1
  3c:   f2 0f 10 1d 50 08 00    movsd  0x850,%xmm3
  43:   00
  44:   f2 0f 5c ca             subsd  %xmm2,%xmm1
  48:   66 0f 28 15 40 08 00    movapd 0x840,%xmm2
  4f:   00
  50:   f2 0f 59 d9             mulsd  %xmm1,%xmm3
  54:   66 0f 14 c9             unpcklpd %xmm1,%xmm1
  58:   81 c2 00 76 1c 00       add    $0x1c7600,%edx
  5e:   66 0f 28 e0             movapd %xmm0,%xmm4
  62:   83 e2 3f                and    $0x3f,%edx
  65:   66 0f 28 2d 30 08 00    movapd 0x830,%xmm5
  6c:   00
  6d:   8d 05 00 00 00 00       lea    0x0,%eax
  73:   c1 e2 05                shl    $0x5,%edx
  76:   03 c2                   add    %edx,%eax
  78:   66 0f 59 d1             mulpd  %xmm1,%xmm2
  7c:   f2 0f 5c c3             subsd  %xmm3,%xmm0
  80:   f2 0f 59 0d 60 08 00    mulsd  0x860,%xmm1
  87:   00
  88:   f2 0f 5c e3             subsd  %xmm3,%xmm4
  8c:   f2 0f 10 78 08          movsd  0x8(%eax),%xmm7
  91:   66 0f 14 c0             unpcklpd %xmm0,%xmm0
  95:   66 0f 28 dc             movapd %xmm4,%xmm3
  99:   f2 0f 5c e2             subsd  %xmm2,%xmm4
  9d:   66 0f 59 e8             mulpd  %xmm0,%xmm5
  a1:   66 0f 5c c2             subpd  %xmm2,%xmm0
  a5:   66 0f 28 35 10 08 00    movapd 0x810,%xmm6
  ac:   00
  ad:   f2 0f 59 fc             mulsd  %xmm4,%xmm7
  b1:   f2 0f 5c dc             subsd  %xmm4,%xmm3
  b5:   66 0f 59 e8             mulpd  %xmm0,%xmm5
  b9:   66 0f 59 c0             mulpd  %xmm0,%xmm0
  bd:   f2 0f 5c da             subsd  %xmm2,%xmm3
  c1:   66 0f 28 10             movapd (%eax),%xmm2
  c5:   f2 0f 5c cb             subsd  %xmm3,%xmm1
  c9:   f2 0f 10 58 18          movsd  0x18(%eax),%xmm3
  ce:   f2 0f 58 d3             addsd  %xmm3,%xmm2
  d2:   f2 0f 5c fa             subsd  %xmm2,%xmm7
  d6:   f2 0f 59 d4             mulsd  %xmm4,%xmm2
  da:   66 0f 59 f0             mulpd  %xmm0,%xmm6
  de:   f2 0f 59 dc             mulsd  %xmm4,%xmm3
  e2:   66 0f 59 d0             mulpd  %xmm0,%xmm2
  e6:   66 0f 59 c0             mulpd  %xmm0,%xmm0
  ea:   66 0f 58 2d 20 08 00    addpd  0x820,%xmm5
  f1:   00
  f2:   f2 0f 59 20             mulsd  (%eax),%xmm4
  f6:   66 0f 58 35 00 08 00    addpd  0x800,%xmm6
  fd:   00
  fe:   66 0f 59 e8             mulpd  %xmm0,%xmm5
102:   66 0f 28 c3             movapd %xmm3,%xmm0
106:   f2 0f 58 58 08          addsd  0x8(%eax),%xmm3
10b:   66 0f 59 cf             mulpd  %xmm7,%xmm1
10f:   66 0f 28 fc             movapd %xmm4,%xmm7
113:   f2 0f 58 e3             addsd  %xmm3,%xmm4
117:   66 0f 58 f5             addpd  %xmm5,%xmm6
11b:   f2 0f 10 68 08          movsd  0x8(%eax),%xmm5
120:   f2 0f 5c eb             subsd  %xmm3,%xmm5
124:   f2 0f 5c dc             subsd  %xmm4,%xmm3
128:   f2 0f 58 48 10          addsd  0x10(%eax),%xmm1
12d:   66 0f 59 f2             mulpd  %xmm2,%xmm6
131:   f2 0f 58 e8             addsd  %xmm0,%xmm5
135:   f2 0f 58 df             addsd  %xmm7,%xmm3
139:   f2 0f 58 cd             addsd  %xmm5,%xmm1
13d:   f2 0f 58 cb             addsd  %xmm3,%xmm1
141:   f2 0f 58 ce             addsd  %xmm6,%xmm1
145:   66 0f 15 f6             unpckhpd %xmm6,%xmm6
149:   f2 0f 58 ce             addsd  %xmm6,%xmm1
14d:   f2 0f 58 e1             addsd  %xmm1,%xmm4
151:   66 0f 28 c4             movapd %xmm4,%xmm0
155:   eb 72                   jmp    1c9 <__libm_sse2_sin+0x1c9>
157:   7f 2e                   jg     187 <__libm_sse2_sin+0x187>
159:   c1 e8 04                shr    $0x4,%eax
15c:   3d fd fc ff 0f          cmp    $0xffffcfd,%eax
161:   75 0a                   jne    16d <__libm_sse2_sin+0x16d>
163:   f2 0f 59 05 b0 08 00    mulsd  0x8b0,%xmm0
16a:   00
16b:   eb 5c                   jmp    1c9 <__libm_sse2_sin+0x1c9>
16d:   f2 0f 10 1d 90 08 00    movsd  0x890,%xmm3
174:   00
175:   f2 0f 59 d8             mulsd  %xmm0,%xmm3
179:   f2 0f 5c d8             subsd  %xmm0,%xmm3
17d:   f2 0f 59 1d a0 08 00    mulsd  0x8a0,%xmm3
184:   00
185:   eb 42                   jmp    1c9 <__libm_sse2_sin+0x1c9>
187:   66 0f c5 c0 03          pextrw $0x3,%xmm0,%eax
18c:   25 f0 7f 00 00          and    $0x7ff0,%eax
191:   3d f0 7f 00 00          cmp    $0x7ff0,%eax
196:   74 29                   je     1c1 <__libm_sse2_sin+0x1c1>
198:   83 ec 20                sub    $0x20,%esp
19b:   f2 0f 11 04 24          movsd  %xmm0,(%esp)
1a0:   8d 44 24 20             lea    0x20(%esp),%eax
1a4:   89 44 24 08             mov    %eax,0x8(%esp)
1a8:   b8 02 00 00 00          mov    $0x2,%eax
1ad:   89 44 24 0c             mov    %eax,0xc(%esp)
1b1:   e8 fc ff ff ff          call   1b2 <__libm_sse2_sin+0x1b2>
1b6:   83 c4 20                add    $0x20,%esp
1b9:   f3 0f 7e 44 24 08       movq   0x8(%esp),%xmm0
1bf:   eb 08                   jmp    1c9 <__libm_sse2_sin+0x1c9>
1c1:   f2 0f 59 05 c0 08 00    mulsd  0x8c0,%xmm0
1c8:   00
1c9:   8b e5                   mov    %ebp,%esp
1cb:   5d                      pop    %ebp
1cc:   c3                      ret   
1cd:   90                      nop   
1ce:   90                      nop   
1cf:   90                      nop
作者: Prescott    时间: 2007-12-19 15:21
哈哈，又不要你写，水平不够就老老实实用Intel给的函数啊。

你觉得这么一大堆指令会比三条fld/fsin/fstp慢？？哈哈哈。不是指令多就慢的。小朋友。
作者: sunboy77_2000    时间: 2007-12-19 15:23
英特尔好像是说的希望下代游戏机不需要手柄。
作者: 紫色    时间: 2007-12-19 15:28
P大屡屡说gcc弱智，我看了gcc的论坛，关于simd intrinsic讨论的太多了。我认为真要解决它，他们有那个水平 ，gcc没有实现更多的是考虑ABI接口兼容性。目前gcc确实是没有simd内置函数实现（对32位），加-mno-80387编译libm就肯定失败。我很期待能早点实现。

[ 本帖最后由 紫色 于 2007-12-24 20:55 编辑 ]
作者: Prescott    时间: 2007-12-19 15:29
  1. #include <stdio.h>
  2. #include <math.h>
  3. #include <sys/time.h>


  6. int main()
  7. {
  8.     struct timeval tv1, tv2;
  9.     int timediff;

  11.     gettimeofday(&tv1, NULL);
  12.     double x=0.654625453;
  13.     double y=0.754656857;
  14.     double z1, z2;
  15.     int i;
  16.     for(i=1;i<20000000;i++){
  17.         z1 = sin(x);
  18.         z2 = sin(y);
  19.         x += z1 + 1;
  20.         y += z2 + 1;
  21.     }
  22.     gettimeofday(&tv2, NULL);
  23.     timediff = (tv2.tv_sec - tv1.tv_sec) * 1000000 + (tv2.tv_usec - tv1.tv_usec);
  24.     printf("sse2: x=%e, y=%e, time = %dn", x, y, timediff);

  26.     gettimeofday(&tv1, NULL);
  27.     x=0.654625453;
  28.     y=0.754656857;
  29.     asm("finit");
  30.     for(i=1;i<20000000;i++){
  31.         asm("fsin" : "=t" (z1) : "0" (x));
  32.         asm("fsin" : "=t" (z2) : "0" (y));
  33.         x += z1 + 1;
  34.         y += z2 + 1;
  35.     }
  36.     gettimeofday(&tv2, NULL);
  37.     timediff = (tv2.tv_sec - tv1.tv_sec) * 1000000 + (tv2.tv_usec - tv1.tv_usec);
  38.     printf("fsin: x=%e, y=%e, time = %dn", x, y, timediff);

  40.     return 0;
  41. }

  43. [xxxx@xxxx-desk test]$ icc -g -fast -static sin_test1.cpp -lstdc++ -lm -o sin.icc
  44. ipo: remark #11001: performing single-file optimizations
  45. ipo: remark #11006: generating assembly file /tmp/iccPOwkIPas_.s
  46. [xxxx@xxxx-desk test]$ ./sin.icc
  47. sse2: x=4.712389e+00, y=4.712389e+00, time = 758399
  48. fsin: x=4.712389e+00, y=4.712389e+00, time = 1254399
复制代码
:whistling: :whistling:

[ 本帖最后由 Prescott 于 2007-12-19 15:34 编辑 ]
作者: 紫色    时间: 2007-12-19 15:50
能不能跑一下简单点的。
   double x=0.0,tmp
    for(i=1;i<20000000;i++){
        tmp = sin(i);
        x=x+tmp;
     }

我就是跑这个的时候发现fsin快些。
cpu : p4 2.4B (478)，支持sse ,sse2
作者: Prescott    时间: 2007-12-19 16:03
  1. #include <stdio.h>
  2. #include <math.h>
  3. #include <sys/time.h>


  6. int main()
  7. {
  8.     struct timeval tv1, tv2;
  9.     int timediff;
  10.     int i;
  11.     double x = 0.0,tmp;
  12.     gettimeofday(&tv1, NULL);
  13.     for(i=1;i<20000000;i++){
  14.         tmp = sin(i);
  15.         x += tmp;
  16.     }
  17.     gettimeofday(&tv2, NULL);
  18.     timediff = (tv2.tv_sec - tv1.tv_sec) * 1000000 + (tv2.tv_usec - tv1.tv_usec);
  19.     printf("sse2: x=%e, time = %dn", x, timediff);

  21.     x=0.0;
  22.     gettimeofday(&tv1, NULL);
  23.     asm("finit");
  24.     for(i=1;i<20000000;i++){
  25.         tmp = (double)i;
  26.         asm("fsin" : "=t" (tmp) : "0" (tmp));
  27.         x += tmp;
  28.     }
  29.     gettimeofday(&tv2, NULL);
  30.     timediff = (tv2.tv_sec - tv1.tv_sec) * 1000000 + (tv2.tv_usec - tv1.tv_usec);
  31.     printf("fsin: x=%e, time = %dn", x, timediff);

  33.     return 0;
  34. }

  36. [xxxxxx@xxxxxx-desk test]$ icc -O3 -xT -static sin_test2.cpp -o sin2.icc
  37. sin_test2.cpp(13): (col. 5) remark: LOOP WAS VECTORIZED.

  39. [xxxxxx@xxxxx-desk test]$ ./sin2.icc
  40. sse2: x=7.052914e-01, time = 387135
  41. fsin: x=7.052914e-01, time = 786659

  43. [xxx@xx-desk ~]$ cat /proc/cpuinfo
  44. ............
  45. model name      : Intel(R) Core(TM)2 CPU          6700  @ 2.66GHz
  46. ............
复制代码
:whistling:

[ 本帖最后由 Prescott 于 2007-12-19 16:06 编辑 ]
作者: 紫色    时间: 2007-12-19 16:27
我也跑了几遍。第一个程序确实是sse快。而我要求改的那个就反过来。
结果是这样（gfortran是必然要调用fsin的，所以使用x87或者sse运行时间几乎完全相同，已经验证）：
第一个：
编译选项                         运行时间
ifort -O3 1.f90                 1.836s
ifort -O3 -xN 1.f90            1.468s
gfortran 1.f90 -O3 ...         4.240s

第二个：
ifort -O3 2.f90                 2.436s
ifort -O3 -xN 2.f90            6.280s   :funk:  ????
gfortran 2.f90 -O3 ...         1.872s

第二个怎么解释，矢量sse那么慢
作者: 慕尼黑    时间: 2007-12-19 16:36
:funk: :funk: :funk::wacko: :wacko: :wacko:
作者: Prescott    时间: 2007-12-19 16:49
原帖由 紫色 于 2007-12-19 16:27 发表
我也跑了几遍。第一个程序确实是sse快。而我要求改的那个就反过来。
结果是这样（gfortran是必然要调用fsin的，所以使用x87或者sse运行时间几乎完全相同，已经验证）：
第一个：
编译选项                        ...


sorry, 不会fortran，帖完整程序吧
作者: 紫色    时间: 2007-12-19 17:39
我的c是菜鸟水平，平时都用fortran。两个程序其实就是你的两个c程序，你一看就懂。
当然，去掉了内联汇编，因为fortran不支持那玩意。
看运行时间用time命令。

1.f90
--------------------
program test1
implicit none
integer,parameter :: N=10000000
real(kind=8) :: x=0.6333333333_8,y=0.7222222222_8
real(kind=8) :: z1,z2
integer :: i
do i=1,N
z1 = sin(x)
z2 = sin(y)
x = z1 + 1.0_8
y = z2 + 1.0_8
enddo
print *, x
end program test2


2.f90
--------------------------
program test2
implicit none
integer,parameter :: N=10000000
real(kind=8) :: x=0.0_8,tmp=0.0_8
integer :: i
do i=1,N
tmp=sin(dble(i))
x=x+tmp
enddo
print *, x
end program test2

已经很有收获了。第二种情况可能是dble花时间（程序中的dble是不可省的）, 或者数据在通用寄存器与sse寄存器之间频繁交换导致延迟。我会逐一测试其他常用函数。

[ 本帖最后由 紫色 于 2007-12-19 17:47 编辑 ]
作者: Prescott    时间: 2007-12-19 17:47
  1. [xxx@xxx-desk test]$ gfortran -O3 2.f90 -o sin.gfort
  2. [xxx@xxx-desk test]$ time ./sin.gfort
  3.    1.95589140854128

  5. real    0m0.654s
  6. user    0m0.653s
  7. sys     0m0.001s
  8. [xxx@xxx-desk test]$ ifort -fast 2.f90 -o sin.ifort
  9. ipo: remark #11001: performing single-file optimizations
  10. ipo: remark #11005: generating object file /tmp/ipo_ifortVMo2Ue.o
  11. 2.f90(6): (col. 1) remark: LOOP WAS VECTORIZED.
  12. [xxx@xxx-desk test]$ time ./sin.ifort
  13.    1.95589140854105

  15. real    0m0.198s
  16. user    0m0.197s
  17. sys     0m0.001s

  19. [xxx@xxx-desk test]$ ifort --version
  20. ifort (IFORT) 10.1 20070913
  21. Copyright (C) 1985-2007 Intel Corporation.  All rights reserved.
复制代码
在我这里，E6700/EL5/ifort ver10.1.008下，ifort比gfortran快3倍多。

[ 本帖最后由 Prescott 于 2007-12-19 17:49 编辑 ]
作者: 紫色    时间: 2007-12-19 17:52
看见了。

只有core2才能跑-fast吧，
能不能分别-O3 和 -O3 -xN看看结果。

我的ifort版本更新些呢
Intel(R) Fortran Compiler for applications running on IA-32, Version 10.1    Build 20071116 Package ID: l_fc_p_10.1.011

[ 本帖最后由 紫色 于 2007-12-19 17:53 编辑 ]
作者: Prescott    时间: 2007-12-19 17:55
原帖由 紫色 于 2007-12-19 17:52 发表
看见了。

只有core2才能跑-fast吧，
能不能分别-O3 和 -O3 -xN看看结果。

我的ifort版本更新些呢
Intel(R) Fortran Compiler for applications running on IA-32, Version 10.1    Build 20071116 Package I ...
  1. [xxx@xxx-desk test]$ ifort -O3 -xW 2.f90 -o sin2.ifort
  2. 2.f90(6): (col. 1) remark: LOOP WAS VECTORIZED.
  3. [xxx@xxx-desk test]$ time ./sin2.ifort
  4.    1.95589140854105

  6. real    0m0.195s
  7. user    0m0.193s
  8. sys     0m0.002s

  10. [xxx@xxx-desk test]$ ifort -O3 2.f90 -o sin2.ifort
  11. 2.f90(6): (col. 1) remark: LOOP WAS VECTORIZED.
  12. [xxx@xxx-desk test]$ time ./sin2.ifort
  13.    1.95589140854105

  15. real    0m0.194s
  16. user    0m0.194s
  17. sys     0m0.000s
复制代码
-xW等同于-xN，新版的ifort取消了-xN。
最后的那一点点性能差别，应该是speedstep造成的，我的处理器自动降频到1.6GHz

哦，你是32位的 :lol:

[ 本帖最后由 Prescott 于 2007-12-19 18:06 编辑 ]
作者: 紫色    时间: 2007-12-19 18:01
满点高兴。把gfortran版本告诉我。按我的选项跑跑看，时间就不是3倍了。
作者: Prescott    时间: 2007-12-19 18:04
原帖由 紫色 于 2007-12-19 18:01 发表
满点高兴。把gfortran版本告诉我。按我的选项跑跑看，时间就不是3倍了。
  1. [xxx@xxx-desk test]$ ifort --version
  2. ifort (IFORT) 10.1 20070913
  3. Copyright (C) 1985-2007 Intel Corporation.  All rights reserved.

  5. [xxx@xxx-desk test]$ gfortran --version
  6. GNU Fortran 95 (GCC) 4.1.1 20070105 (Red Hat 4.1.1-52)
  7. Copyright (C) 2006 Free Software Foundation, Inc.

  9. GNU Fortran comes with NO WARRANTY, to the extent permitted by law.
  10. You may redistribute copies of GNU Fortran
  11. under the terms of the GNU General Public License.
  12. For more information about these matters, see the file named COPYING

  14. [xxx@xxx-desk test]$ uname -a
  15. Linux xxx-desk 2.6.18-8.el5 #1 SMP Fri Jan 26 14:15:14 EST 2007 x86_64 x86_64 x86_64 GNU/Linux
复制代码
作者: 紫色    时间: 2007-12-19 18:10
gfortran 2.f90 -Wall -g -fverbose-asm -fbounds-check -O3 -ftree-vectorize -funroll-all-loops -mfpmath=sse -march=pentium4 -mfpmath=sse -ffixed-line-length-none -ffree-line-length-none -fimplicit-none -frange-check -Wimplicit-interface -Wunused-parameter -Wconversion -Wunderflow -ftree-vectorizer-verbose=5 -ffpe-trap=invalid,zero,overflow,underflow

这就是一行，不要断行。请看看时间吧。
作者: 紫色    时间: 2007-12-19 18:16
好像也差不了多少。
作者: Prescott    时间: 2007-12-19 18:22
原帖由 紫色 于 2007-12-19 18:10 发表
gfortran 2.f90 -Wall -g -fverbose-asm -fbounds-check -O3 -ftree-vectorize -funroll-all-loops -mfpmath=sse -march=pentium4 -mfpmath=sse -ffixed-line-length-none -ffree-line-length-none -fimplicit-none  ...
  1. [xxx@xxx-desk test]$ gfortran 2.f90 -Wall -g -fverbose-asm -fbounds-check -O3 -ftree-vectorize -funroll-all-loops -mfpmath=sse -march=pentium4 -mfpmath=sse -ffixed-line-length-none -ffree-line-length-none -fimplicit-none -frange-check -Wimplicit-interface -Wunused-parameter -Wconversion -Wunderflow -ftree-vectorizer-verbose=5 -ffpe-trap=invalid,zero,overflow,underflow
  2. 2.f90:0: error: CPU you selected does not support x86-64 instruction set
  3. 2.f90:0: error: CPU you selected does not support x86-64 instruction set

  5. 只好加个-m32了

  7. [xxx@xxx-desk test]$ gfortran 2.f90 -m32 -Wall -g -fverbose-asm -fbounds-check -O3 -ftree-vectorize -funroll-all-loops -mfpmath=sse -march=pentium4 -mfpmath=sse -ffixed-line-length-none -ffree-line-length-none -fimplicit-none -frange-check -Wimplicit-interface -Wunused-parameter -Wconversion -Wunderflow -ftree-vectorizer-verbose=5 -ffpe-trap=invalid,zero,overflow,underflow

  9. 2.f90:6: note: not vectorized: unsupported use in stmt.
  10. 2.f90:6: note: vectorized 0 loops in function.
  11. [xxx@xxx-desk test]$ time ./a.out
  12.    1.95589140854113

  14. real    0m0.450s
  15. user    0m0.450s
  16. sys     0m0.001s

  18. 原来gfortran生成的64位代码还要更慢。。。。。。嘿嘿，说他弱智不算冤枉吧。

  20. [xxx@xxx-desk test]$ gfortran 2.f90 -Wall -g -fverbose-asm -fbounds-check -O3 -ftree-vectorize -funroll-all-loops -mfpmath=sse -march=nocona -mfpmath=sse -ffixed-line-length-none -ffree-line-length-none -fimplicit-none -frange-check -Wimplicit-interface -Wunused-parameter -Wconversion -Wunderflow -ftree-vectorizer-verbose=5 -ffpe-trap=invalid,zero,overflow,underflow

  22. 2.f90:6: note: not vectorized: unsupported use in stmt.
  23. 2.f90:6: note: vectorized 0 loops in function.
  24. [xxx@xxx-desk test]$ time ./a.out
  25.    1.95589140854128

  27. real    0m0.691s
  28. user    0m0.689s
  29. sys     0m0.001s
复制代码
作者: Prescott    时间: 2007-12-19 18:26
原帖由 紫色 于 2007-12-19 18:16 发表
好像也差不了多少。

也不是很多，0.45/0.19 = 2.36
确实不多。。。:sweatingbullets:
作者: 紫色    时间: 2007-12-19 18:40
2.36:1？我跑第二个程序连1:1都没有。不过今天我已经很高兴了，因为看到1.f90确实是ifort快，那也是有sin的程序。

非常感谢P大！

[ 本帖最后由 紫色 于 2007-12-24 20:58 编辑 ]



欢迎光临 POPPUR爱换 (https://we.poppur.com/) Powered by Discuz! X3.4