三浦です

私のメールでプログラムまで作ってもらってありがとうございます。
ロールじゃなくてローテートですね。間違えて覚えていました。

ローテート命令はインラインアセンブラで直接書かないと出てこないかな
と思っていたのですが、出てくるんですね。びっくりしました。
それでも、アセンブルリストを見ると何行か削れそうな部分が出てきて
手を入れたくなります。Rubyではアセンブリ言語はご法度なんでしょ
うが。

昔、Green Hills Cで
 (x > 0) ? x : -x
みたいなコードが、
  (x >> 31) ^ x + ((unsigned)x >> 31)
みたいにコンパイルされてびっくりしたことがあります。


----- Original Message ----- 
From: "SASADA Koichi" <ko1 atdot.net>
To: <yarv-dev quickml.atdot.net>
Sent: Wednesday, November 21, 2007 7:44 PM
Subject: [yarv-dev:1145] Re: 考えている案


> 　ささだです．
>
> 　なるほど，rotate を使うと，うまいこと削りたいところが頭出しされる
> のですね．すごい．頭が悪いので，こういう bit 演算が苦手なんです．
>
> 　というわけで，ちょっと作ってみました．
>
> typedef unsigned long long VALUE;
>
> VALUE
> d2v(double d)
> {
>    union {
> double d;
> VALUE v;
>    } u;
>    u.d = d;
>    return u.v;
> }
>
> double
> v2d(VALUE v)
> {
>    union {
> double d;
> VALUE v;
>    } u;
>    u.v = v;
>    return u.d;
> }
>
> #define ROTL(v, n) (((v) << (n)) | ((v) >> ((sizeof(v) * 8) - n)))
> #define ROTR(v, n) (((v) >> (n)) | ((v) << ((sizeof(v) * 8) - n)))
>
> VALUE
> DOUBLE2VALUE(double d)
> {
>    VALUE v1 = d2v(d);
>    VALUE v2 = ROTL(v1, 3);
>
>    if ((v2 & 0x01) ^ ((v2 >> 1) & 0x01)) {
> return v2 | 0x01;
>    }
>    else {
> /* from heap */
> return -1;
>    }
> }
>
> double
> VALUE2DOUBLE(VALUE v)
> {
>    VALUE v1 = v ^ ((v >> 1) & 0x01);
>    VALUE v2 = ROTR(v1, 3);
>    return v2d(v2);
> }
>
> 　x64 用 gcc -O3 の結果です．
>
> DOUBLE2VALUE:
> .LFB16:
> movd %xmm0, %rcx
> movq $-1, %rax
> movq %xmm0, -8(%rsp)  # これ無駄では...
> rolq $3, %rcx
> movq %rcx, %rdx
> movq %rcx, %rsi
> orq $1, %rcx
> shrq %rdx
> andl $1, %esi
> andl $1, %edx
> cmpq %rdx, %rsi
> cmovne %rcx, %rax
> ret
>
> VALUE2DOUBLE:
> .LFB17:
> movq %rdi, %rax
> shrq %rax
> andl $1, %eax
> xorq %rax, %rdi
> rorq $3, %rdi
> movq %rdi, -8(%rsp)
> movlpd -8(%rsp), %xmm1
> movsd %xmm1, %xmm0
> ret
>
>
> 感動したところ：
> ・jmp 命令が無い！（まぁ，コード小さいから）
> ・rol, ror 命令使ってるよ！
> ・短い！
>
> 　いやぁ，凄い．実は，rotate 命令を使うかどうか，昨日ちょっと調べて
> みると，gcc 2.95 でも -O0 で rotate 命令を使うようなコードを出してた
> んですよね．凄い．
>
> 　実際に速度比較やってみないとなんとも言えませんが，メモリアクセスが
> なくなる分，速くなったりしないかなー．
>



--
ML: yarv-dev quickml.atdot.net
Info: http://www.atdot.net/~ko1/quickml