Visual Studio 2008 Part 22

[過去ﾛｸﾞ] Visual Studio 2008 Part 22 (314ﾚｽ)
上下前次 1-新
通常表示 512ﾊﾞｲﾄ分割ﾚｽ栞

このｽﾚｯﾄﾞは過去ﾛｸﾞ倉庫に格納されています｡
次ｽﾚ検索歴削→次ｽﾚ栞削→次ｽﾚ過去ﾛｸﾞﾒﾆｭｰ

235: デフォルトの名無しさん [sage] 2018/09/16(日) 16:37:19.30 ID:zL1WUjLu

>>231
逆アセンブラ

for (int i=0;i<num;i++) norm += (double)r[i] * (double)r[i];
00000033 33 D2            xor         edx,edx
00000035 89 55 E8         mov         dword ptr [ebp-18h],edx
00000038 90               nop
00000039 EB 03            jmp         0000003E
0000003b FF 45 E8         inc         dword ptr [ebp-18h]
0000003e 8B 45 E8         mov         eax,dword ptr [ebp-18h]
00000041 3B 45 FC         cmp         eax,dword ptr [ebp-4]
00000044 7D 1B            jge         00000061
00000046 8B 45 F8         mov         eax,dword ptr [ebp-8]
00000049 8B 55 E8         mov         edx,dword ptr [ebp-18h]
0000004c DD 04 D0         fld         qword ptr [eax+edx*8]
0000004f 8B 45 F8         mov         eax,dword ptr [ebp-8]
00000052 8B 55 E8         mov         edx,dword ptr [ebp-18h]
00000055 DC 0C D0         fmul        qword ptr [eax+edx*8]
00000058 DC 45 F0         fadd        qword ptr [ebp-10h]
0000005b DD 5D F0         fstp        qword ptr [ebp-10h]
0000005e 90               nop
0000005f EB DA            jmp         0000003B
norm = sqrt(norm);
00000061 DD 45 F0         fld         qword ptr [ebp-10h]
00000064 83 EC 08         sub         esp,8
00000067 DD 1C 24         fstp        qword ptr [esp]
0000006a E8 0D 50 7B FF   call        FF7B507C
0000006f DD 5D D8         fstp        qword ptr [ebp-28h]
00000072 DD 45 D8         fld         qword ptr [ebp-28h]
00000075 DD 5D F0         fstp        qword ptr [ebp-10h]

http://mevius.5ch.net/test/read.cgi/tech/1413180800/235

236: デフォルトの名無しさん [sage] 2018/09/16(日) 16:37:34.94 ID:zL1WUjLu

>>231
逆アセンブラ（続き）

if (regulate) for (int i=0;i<num;i++) r[i] = (T)(r[i]/norm);
00000078 0F B6 45 08      movzx       eax,byte ptr [ebp+8]
0000007c 85 C0            test        eax,eax
0000007e 74 25            je          000000A5
00000080 33 D2            xor         edx,edx
00000082 89 55 EC         mov         dword ptr [ebp-14h],edx
00000085 90               nop
00000086 EB 03            jmp         0000008B
00000088 FF 45 EC         inc         dword ptr [ebp-14h]
0000008b 8B 45 EC         mov         eax,dword ptr [ebp-14h]
0000008e 3B 45 FC         cmp         eax,dword ptr [ebp-4]
00000091 7D 12            jge         000000A5
00000093 8B 45 F8         mov         eax,dword ptr [ebp-8]
00000096 8B 55 EC         mov         edx,dword ptr [ebp-14h]
00000099 DD 45 F0         fld         qword ptr [ebp-10h]
0000009c DC 3C D0         fdivr       qword ptr [eax+edx*8]
0000009f DD 1C D0         fstp        qword ptr [eax+edx*8]
000000a2 90               nop
000000a3 EB E3            jmp         00000088
return norm;
000000a5 DD 45 F0         fld         qword ptr [ebp-10h]
000000a8 DD 5D E0         fstp        qword ptr [ebp-20h]

http://mevius.5ch.net/test/read.cgi/tech/1413180800/236

244: デフォルトの名無しさん [sage] 2018/09/16(日) 17:22:07.21 ID:zL1WUjLu

>>240
現象確認した。こちらでも再現した。
逆アセンブルは、以下。(肝心のループ部分は次レス内）

正直、fld/fmul/fadd/fstpのループ部分は変わらず、
normのアドレスが [ebp-10h]（つまりローカル）から
ds:[00A4AD40h]（つまりグローバル）に変わっただけであり、
これで結果が変わるのはかなり奇妙な気もするが、何か見落としがあるのかも。

>>240逆アセンブル(static付加版）
template<typename T> static double calc_norm_and_regulate(int num, T* r, bool regulate){ // <float> for debug.
static double norm = 0;
for (int i=0;i<num;i++) norm += (double)r[i] * (double)r[i];
00000000 55               push        ebp
00000001 8B EC            mov         ebp,esp
00000003 83 EC 20         sub         esp,20h
00000006 89 4D FC         mov         dword ptr [ebp-4],ecx
00000009 89 55 F8         mov         dword ptr [ebp-8],edx
0000000c 83 3D 14 2E 38 00 00 cmp         dword ptr ds:[00382E14h],0
00000013 74 05            je          0000001A
00000015 E8 FF 52 30 68   call        68305319
0000001a 33 D2            xor         edx,edx
0000001c 89 55 F0         mov         dword ptr [ebp-10h],edx
0000001f 33 D2            xor         edx,edx
00000021 89 55 F4         mov         dword ptr [ebp-0Ch],edx
00000024 D9 EE            fldz
00000026 DD 5D E8         fstp        qword ptr [ebp-18h]
00000029 33 D2            xor         edx,edx
0000002b 89 55 F0         mov         dword ptr [ebp-10h],edx
0000002e 90               nop
0000002f EB 03            jmp         00000034

http://mevius.5ch.net/test/read.cgi/tech/1413180800/244

245: デフォルトの名無しさん [sage] 2018/09/16(日) 17:22:29.32 ID:zL1WUjLu

>>240逆アセンブル（続き）(static付加版）

00000031 FF 45 F0         inc         dword ptr [ebp-10h]
00000034 8B 45 F0         mov         eax,dword ptr [ebp-10h]
00000037 3B 45 FC         cmp         eax,dword ptr [ebp-4]
0000003a 7D 21            jge         0000005D
0000003c 8B 45 F8         mov         eax,dword ptr [ebp-8]
0000003f 8B 55 F0         mov         edx,dword ptr [ebp-10h]
00000042 DD 04 D0         fld         qword ptr [eax+edx*8]
00000045 8B 45 F8         mov         eax,dword ptr [ebp-8]
00000048 8B 55 F0         mov         edx,dword ptr [ebp-10h]
0000004b DC 0C D0         fmul        qword ptr [eax+edx*8]
0000004e DC 05 40 AD A4 00 fadd        qword ptr ds:[00A4AD40h]
00000054 DD 1D 40 AD A4 00 fstp        qword ptr ds:[00A4AD40h]
0000005a 90               nop
0000005b EB D4            jmp         00000031
norm = sqrt(norm);
0000005d DD 05 40 AD A4 00 fld         qword ptr ds:[00A4AD40h]
00000063 83 EC 08         sub         esp,8
00000066 DD 1C 24         fstp        qword ptr [esp]
00000069 E8 0E 50 88 FF   call        FF88507C
0000006e DD 5D E0         fstp        qword ptr [ebp-20h]
00000071 DD 45 E0         fld         qword ptr [ebp-20h]
00000074 DD 1D 40 AD A4 00 fstp        qword ptr ds:[00A4AD40h]

http://mevius.5ch.net/test/read.cgi/tech/1413180800/245

246: デフォルトの名無しさん [sage] 2018/09/16(日) 17:22:46.06 ID:zL1WUjLu

>>240逆アセンブル（続き）(static付加版）
if (regulate) for (int i=0;i<num;i++) r[i] = (T)(r[i]/norm);
0000007a 0F B6 45 08      movzx       eax,byte ptr [ebp+8]
0000007e 85 C0            test        eax,eax
00000080 74 28            je          000000AA
00000082 33 D2            xor         edx,edx
00000084 89 55 F4         mov         dword ptr [ebp-0Ch],edx
00000087 90               nop
00000088 EB 03            jmp         0000008D
0000008a FF 45 F4         inc         dword ptr [ebp-0Ch]
0000008d 8B 45 F4         mov         eax,dword ptr [ebp-0Ch]
00000090 3B 45 FC         cmp         eax,dword ptr [ebp-4]
00000093 7D 15            jge         000000AA
00000095 8B 45 F8         mov         eax,dword ptr [ebp-8]
00000098 8B 55 F4         mov         edx,dword ptr [ebp-0Ch]
0000009b DD 05 40 AD A4 00 fld         qword ptr ds:[00A4AD40h]
000000a1 DC 3C D0         fdivr       qword ptr [eax+edx*8]
000000a4 DD 1C D0         fstp        qword ptr [eax+edx*8]
000000a7 90               nop
000000a8 EB E0            jmp         0000008A
return norm;
000000aa DD 05 40 AD A4 00 fld         qword ptr ds:[00A4AD40h]
000000b0 DD 5D E8         fstp        qword ptr [ebp-18h]
}
000000b3 DD 45 E8         fld         qword ptr [ebp-18h]
000000b6 8B E5            mov         esp,ebp
000000b8 5D               pop         ebp
000000b9 C2 04 00         ret         4

http://mevius.5ch.net/test/read.cgi/tech/1413180800/246

249: デフォルトの名無しさん [sage] 2018/09/16(日) 20:54:27.39 ID:zL1WUjLu

>>240
さて再見したが、やはりstaticだけで直る理由は分からない。
なお、最適化ミスの場合は、逆アセンブラを読めば分かる。
今のところそれではない。

一応、>>191ソースのtemplate部の逆アセンブルを上げておく。（ただし重複するので頭のみ）
頭はこれ。続きが>>235,236。

template<typename T> static double calc_norm_and_regulate(int num, T* r, bool regulate){ // <float> for debug.
double norm = 0;
00000000 55               push        ebp
00000001 8B EC            mov         ebp,esp
00000003 83 EC 28         sub         esp,28h
00000006 89 4D FC         mov         dword ptr [ebp-4],ecx
00000009 89 55 F8         mov         dword ptr [ebp-8],edx
0000000c 83 3D 14 2E 76 00 00 cmp         dword ptr ds:[00762E14h],0
00000013 74 05            je          0000001A
00000015 E8 FF 52 1B 68   call        681B5319
0000001a 33 D2            xor         edx,edx
0000001c 89 55 E8         mov         dword ptr [ebp-18h],edx
0000001f 33 D2            xor         edx,edx
00000021 89 55 EC         mov         dword ptr [ebp-14h],edx
00000024 D9 EE            fldz
00000026 DD 5D F0         fstp        qword ptr [ebp-10h]
00000029 D9 EE            fldz
0000002b DD 5D E0         fstp        qword ptr [ebp-20h]
0000002e D9 EE            fldz
00000030 DD 5D F0         fstp        qword ptr [ebp-10h]

http://mevius.5ch.net/test/read.cgi/tech/1413180800/249

252: 240 [sage] 2018/09/16(日) 21:43:27.78 ID:/oSJzlqn

static版
0000000e 33 C0            xor         eax,eax
00000010 85 F6            test        esi,esi
00000012 7E 16            jle         0000002A
00000014 DD 04 C7         fld         qword ptr [edi+eax*8]
00000017 DC C8            fmul        st(0),st
00000019 DC 05 00 30 CC 00 fadd        qword ptr ds:[00CC3000h]
0000001f DD 1D 00 30 CC 00 fstp        qword ptr ds:[00CC3000h]
00000025 40               inc         eax
00000026 3B C6            cmp         eax,esi
00000028 7C EA            jl          00000014
0000002a DD 05 00 30 CC 00 fld         qword ptr ds:[00CC3000h]

非static版
0000000e D9 EE            fldz
00000010 33 C0            xor         eax,eax
00000012 85 F6            test        esi,esi
00000014 7E 0C            jle         00000022
00000016 DD 04 C7         fld         qword ptr [edi+eax*8]
00000019 DC C8            fmul        st(0),st
0000001b DE C1            faddp       st(1),st
0000001d 40               inc         eax
0000001e 3B C6            cmp         eax,esi
00000020 7C F4            jl          00000016

http://mevius.5ch.net/test/read.cgi/tech/1413180800/252

253: デフォルトの名無しさん [sage] 2018/09/16(日) 22:27:34.58 ID:zL1WUjLu

>>251
とりあえず落ち着け。一つずつ行こう。

> ループ部分の関数を#pragma unmanagedすると結果が変わるでそれが正しいのかも。
こちらでも確認した。
calc_norm_and_regulateをunmanaged関数にすると、違いはなくなる。
(Releaseビルドの｀をコマンドプロンプトで起動した際にも、****ddの結果となる）

ただしこちらの逆アセンブル結果は以下だ。(fld/fmul/fadd/fstpであることに注意）
for (int i=0;i<num;i++) norm += (double)r[i] * (double)r[i];
0007272C C7 45 F4 00 00 00 00 mov         dword ptr [i],0
00072733 EB 09            jmp         `anonymous namespace'::calc_norm_and_regulate<double>+1Eh (7273Eh)
00072735 8B 45 F4         mov         eax,dword ptr [i]
00072738 83 C0 01         add         eax,1
0007273B 89 45 F4         mov         dword ptr [i],eax
0007273E 8B 4D F4         mov         ecx,dword ptr [i]
00072741 3B 4D 08         cmp         ecx,dword ptr [num]
00072744 7D 1A            jge         `anonymous namespace'::calc_norm_and_regulate<double>+40h (72760h)
00072746 8B 55 F4         mov         edx,dword ptr [i]
00072749 8B 45 0C         mov         eax,dword ptr [r]
0007274C 8B 4D F4         mov         ecx,dword ptr [i]
0007274F 8B 75 0C         mov         esi,dword ptr [r]
00072752 DD 04 D0         fld         qword ptr [eax+edx*8]
00072755 DC 0C CE         fmul        qword ptr [esi+ecx*8]
00072758 DC 45 F8         fadd        qword ptr [norm]
0007275B DD 5D F8         fstp        qword ptr [norm]
0007275E EB D5            jmp         `anonymous namespace'::calc_norm_and_regulate<double>+15h (72735h)

http://mevius.5ch.net/test/read.cgi/tech/1413180800/253

254: デフォルトの名無しさん [sage] 2018/09/16(日) 22:33:32.07 ID:zL1WUjLu

>>252
そちらの逆アセンブルは以下の違いが出てるだろ。
static版： fld/fmul/fadd/fstp
非static版： fld/fmul/faddp (fstpが無い）
この非static版の場合、拡張倍精度(80bit)で演算されるから精度が高いことになり、
static版との演算結果に違いが出るのも仕様通りなんだよ。(これは>>200と同じ間違い)

一応、fstpにも80bit版はあって、Intelのマニュアルによると以下。
> オペコード命令説明
> D9 /2 FST m32fp ST(0) をm32fp にコピーする。
> DD /2 FST m64fp ST(0) をm64fp にコピーする。
> DD D0+i FST ST(i) ST(0) をST(i) にコピーする。
> D9 /3 FSTP m32fp ST(0) をm32fp にコピーし、レジスタスタックをポップする。
> DD /3 FSTP m64fp ST(0) をm64fp にコピーし、レジスタスタックをポップする。
> DB /7 FSTP m80fp ST(0) をm80fp にコピーし、レジスタスタックをポップする。
> DD D8+i FSTP ST(i) ST(0) をST(i) にコピーし、レジスタスタックをポップする。
つまり君のstatic版
> 0000001f DD 1D 00 30 CC 00 fstp qword ptr ds:[00CC3000h]
では FSTP /3 m64fp [disp32] であり、そこで64bit（倍精度）に丸められてる。
だからレジスタ(80bit=拡張倍精度）で演算される非static版と結果が異なる。
static版のsftpが DB /7 m80fp なら誤差は出ないはずなんだよ。（Cでどう書くのかは知らん）

だから>>252の場合の誤差なら、仕様通りなんだよ。(片方が倍精度、もう片方は拡張倍精度）
ただし、>>191は逆アセンブル(>>235)を見る限りそれに該当しないし、(両方とも倍精度）
今回の俺の上記逆アセンブル（>>253、中身は君の指摘通りunmanagedにしただけ）も該当しない。(両方とも倍精度）
そして253は何故か直ってしまった。

http://mevius.5ch.net/test/read.cgi/tech/1413180800/254

上下前次 1-新書関写板覧索設栞歴

あと 59 ﾚｽあります
ｽﾚ情報赤ﾚｽ抽出画像ﾚｽ抽出歴の未読ｽﾚ AAｻﾑﾈｲﾙ

ぬこの手ぬこTOP 0.016s