Visual Studio 2008 Part 22

1 名前：デフォルトの名無しさん mailto:sage [2014/10/13(月) 15:13:20.18 ID:yFRqmPNp.net] ■公式 ttp://www.microsoft.com/japan/msdn/vstudio/default.aspx ■前スレ Visual Studio 2008 Part 21 peace.2ch.net/test/read.cgi/tech/1290969016/ ■リンク Visual Studio 2008に搭載された17の新機能 ttp://www.atmarkit.co.jp/fdotnet/special/visualstudio2008_01/visualstudio2008_01_01.html 5000個のバグと戦った、MSが「Visual Studio 2008」RTM出荷 ttp://www.atmarkit.co.jp/news/200711/20/vs.html ■関連 Visual Studio 2013 part4 peace.2ch.net/test/read.cgi/tech/1404333757/ Visual Studio 2012 Part8 peace.2ch.net/test/read.cgi/tech/1392639689/ Visual Studio 2010 Part21 peace.2ch.net/test/read.cgi/tech/1412136476/ Visual Studio 2005 Part 27 peace.2ch.net/test/read.cgi/tech/1291513609/ その他テンプレ>>2-5 196 名前：デフォルトの名無しさん mailto:sage [2018/09/15(土) 20:42:48.06 ID:UR1d6CKz.net] オプション等(コマンドライン）：（全てデフォのままのはずだが一応） C/C++： 　/GL /D "WIN32" /D "NDEBUG" /D "_UNICODE" /D "UNICODE" /FD /EHa /MD /Yu"stdafx.h" /Fp"Release\test_floating_error4.pch" /Fo"Release\\" /Fd"Release\vc90.pdb" /W3 /nologo /c /Zi /clr /TP /errorReport:prompt /FU "c:\Windows\Microsoft.NET\Framework\v2.0.50727\System.dll" /FU "c:\Windows\Microsoft.NET\Framework\v2.0.50727\System.Data.dll" /FU "c:\Windows\Microsoft.NET\Framework\v2.0.50727\System.XML.dll" リンカ： /OUT:"MYPATH\test_floating_error4\Release\test_floating_error4.exe" /INCREMENTAL:NO /NOLOGO /MANIFEST /MANIFESTFILE:"Release\test_floating_error4.exe.intermediate.manifest" /MANIFESTUAC:"level='asInvoker' uiAccess='false'" /DEBUG /PDB:"MYPATH\test_floating_error4\Release\test_floating_error4.pdb" /LTCG /DYNAMICBASE /FIXED:No /NXCOMPAT /MACHINE:X86 /ERRORREPORT:PROMPT パスは$MYPATHと書き換えた。また、投稿の為に適宜改行を入れた。 このプロジェクト名は見れば分かるとおりtest_floating_error4。 197 名前：デフォルトの名無しさん mailto:sage [2018/09/15(土) 20:43:03.84 ID:UR1d6CKz.net] 備考2： なおこの方法で見える calc_norm_and_regulate 関数の『逆アセンブル』結果は 俺の環境での物とコールアドレス以外は一致していることを確認している。 一応、diff結果は以下。 8c8 < 0000000c cmp dword ptr ds:[00752E14h],0 --- > 0000000c cmp dword ptr ds:[007D2E14h],0 10c10 < 00000015 call 676F58B9 --- > 00000015 call 683C58B9 44c44 < 0000006a call FF0455E8 --- > 0000006a call FFD955E8 198 名前：デフォルトの名無しさん mailto:sage [2018/09/16(日) 02:47:43.28 ID:wIV2HUNW.net] >>191 C++/CLR では、.Net を使っているから、起動方法が違うだけでも、 fpu control word の値や、使うCPU命令がx87 FPUなのか、SSE なのかが違ってくる可能性がある� 199 名前：ｩもしれない。 fpu control word は、main()関数に入る前の start up codeの中で 初期化される。 [] [ここ壊れてます] 200 名前：デフォルトの名無しさん mailto:sage [2018/09/16(日) 03:19:36.22 ID:wIV2HUNW.net] >>192 // Release build // 0.000000, 0x1ff68ddfb62221dd from IDE // 0.000000, 0x1ff68ddfb62221de from command prompt それにしても、随分小さな値だね。ちなみに、浮動小数点表示 の場合の有効数字の桁数を上げたら、どのようになる？ 1.xxxe-yy 表示にして。 201 名前：デフォルトの名無しさん mailto:sage [2018/09/16(日) 03:40:33.66 ID:wIV2HUNW.net] >>192 試しに、ソースの冒頭に #include <stdio.h> を追加してから、 Console::Write(String::Format("{0:F6}, 0x{1:x16}\r\n",norm, *(__int64*)&norm)); の部分を、 printf( "%30.30e, 0x%016X\n", norm, *(__int64*)&norm) ); としてみるとどうなる？ 202 名前：デフォルトの名無しさん mailto:sage [2018/09/16(日) 03:42:01.82 ID:wIV2HUNW.net] >>197 誤： printf( "%30.30e, 0x%016X\n", norm, *(__int64*)&norm) ); 正： printf( "%30.30e, 0x%016X\n", norm, *(__int64*)&norm ); 203 名前：デフォルトの名無しさん [2018/09/16(日) 07:27:22.50 ID:SOVIz+sV.net] > 0x1ff68ddfb62221dd（Debug） > 0x1ff68ddfb62221de（Release） VS 2010 VC++ Express でも再現した 204 名前：デフォルトの名無しさん [2018/09/16(日) 07:38:10.70 ID:SOVIz+sV.net] ↓このﾙｰﾌﾟを抜けたあと、すでにReleaseﾋﾞﾙﾄﾞとDebugﾋﾞﾙﾄﾞでは 　normの値に差異が発生してることが確認できた for (int i=0;i<num;i++) norm += (double)r[i] * (double)r[i]; ↓この下に（ﾙｰﾌﾟ内に）fprintf文を入れるだけで 　ReleaseﾋﾞﾙﾄﾞとDebugﾋﾞﾙﾄﾞが同じ実行結果になることが確認できた norm += (double)r[i] * (double)r[i]; とりあえずまずこれだけは分かったから 低学歴知恵遅れが書いたｳﾝｺｺｰﾄﾞの問題箇所を限定する 205 名前：デフォルトの名無しさん [2018/09/16(日) 07:51:03.41 ID:SOVIz+sV.net] �@-1 ﾃﾞﾌｫﾙﾄ設定（Release） 【ｺｰﾄﾞ】 #include "stdafx.h" #include <stdio.h> #include <stdint.h> #include <math.h> int main(array<System::String ^> ^args) { 　　__int64 inputs_hex[16] = { 　　　　0x0000000000000000, 0x0000000000000000, 0x0000000000000000, 0x1fedb1530240aa54, 　　　　0x0000000000000000, 0x0000000000000000, 0x0000000000000000, 0x1ff0af0d95025bc3, 　　　　0x1fc9353df6af376b, 0x0000000000000000, 0x0000000000000000, 0x0000000000000000, 　　　　0x0000000000000000, 0x0000000000000000, 0x0000000000000000, 0x0000000000000000 　　}; 　　double* r = (double*)inputs_hex; 　　double norm = 0; 　　for (int i = 0; i < 16; i++) { 　　　　norm += (double)r[i] * (double)r[i]; //　　　　fprintf(stdout, "[1]0x%016llX:%.19lg\n", *(uint64_t*)&norm, norm); 　　} 　　fprintf(stdout, "0x%016llX:%.19lg\n", *(uint64_t*)&norm, norm); 　　return 0; } 206 名前：デフォルトの名無しさん [2018/09/16(日) 07:51:56.62 ID:SOVIz+sV.net] �@-1 ﾃﾞﾌｫﾙﾄ設定（Release） ↓このｺｰﾄﾞの逆ｱｾﾝﾌﾞﾙｺｰﾄﾞ https://ideone.com/Ryyxkl 【実行結果】 0x0007F2C44DFFF8F2:1.1053482540585106e-308 207 名前：デフォルトの名無しさん [2018/09/16(日) 07:53:20.79 ID:SOVIz+sV.net] �@-2 ﾃﾞﾌｫﾙﾄ設定（Release） #include "stdafx.h" #include <stdio.h> #include <stdint.h> #include <math.h> int main(array<System::String ^> ^args) { 　　__int64 inputs_hex[16] = { 　　　　0x0000000000000000, 0x0000000000000000, 0x0000000000000000, 0x1fedb1530240aa54, 　　　　0x0000000000000000, 0x0000000000000000, 0x0000000000000000, 0x1ff0af0d95025bc3, 　　　　0x1fc9353df6af376b, 0x0000000000000000, 0x0000000000000000, 0x0000000000000000, 　　　　0x0000000000000000, 0x0000000000000000, 0x0000000000000000, 0x0000000000000000 　　}; 　　double* r = (double*)inputs_hex; 　　double norm = 0; 　　for (int i = 0; i < 16; i++) { 　　　　norm += (double)r[i] * (double)r[i]; 　　　　fprintf(stdout, "[1]0x%016llX:%.19lg\n", *(uint64_t*)&norm, norm); 　　} 　　fprintf(stdout, "0x%016llX:%.19lg\n", *(uint64_t*)&norm, norm); 　　return 0; } 208 名前：デフォルトの名無しさん [2018/09/16(日) 08:02:12.99 ID:SOVIz+sV.net] �@-2 ﾃﾞﾌｫﾙﾄ設定(Release) ↓このｺｰﾄﾞの逆ｱｾﾝﾌﾞﾙｺｰﾄﾞ https://ideone.com/4j12ib 【実行結果】 ↓実行結果を書き込めないからこっちに書き込んどいた https://ideone.com/gOqtki 0x0007F2C44DFFF8F1:1.1053482540585101e-308 209 名前：デフォルトの名無しさん [2018/09/16(日) 08:05:50.92 ID:SOVIz+sV.net] �A-1 最適化無効 (/Od)（Release） ※ ｺｰﾄﾞは�@-1（>>201）と同じ ↓このｺｰﾄﾞの逆ｱｾﾝﾌﾞﾙｺｰﾄﾞ https://ideone.com/cONCAx 【実行結果】 0x0007F2C44DFFF8F1:1.1053482540585101e-308 210 名前：デフォルトの名無しさん [2018/09/16(日) 08:09:59.78 ID:SOVIz+sV.net] �A-2 最適化無効 (/Od)（Release） ※ ｺｰﾄﾞは�A-2（>>203）と同じ ↓このｺｰﾄﾞの逆ｱｾﾝﾌﾞﾙｺｰﾄﾞ https://ideone.com/kuRt3w 【実行結果】 ※ �@-2（>>204）と同じ 0x0007F2C44DFFF8F1:1.1053482540585101e-308 211 名前：デフォルトの名無しさん [2018/09/16(日) 08:20:27.67 ID:SOVIz+sV.net] �@-1、�@-2の逆ｱｾﾝﾌﾞﾙの出力結果を比較すると原型をとどめてないぐらいｸﾞﾁｮｸﾞﾁｮに違う（最適化のせいと考えられる） �A-1、�A-2の逆ｱｾﾝﾌﾞﾙの出力結果を比較すると差異はほとんどない（�Aは両方ともまったく最適化されてないから当然） �@-1と�A-1の逆ｱｾﾝﾌﾞﾙの出力結果を比較すると原型をとどめてないぐらいｸﾞﾁｮｸﾞﾁｮに違う（�@-1のｺｰﾄﾞ(>>201)ははげしくｳﾝｺ最適化されてると考えられる） �@-2と�A-2の逆ｱｾﾝﾌﾞﾙの出力結果を比較すると差異はほとんどない（�@-2のｺｰﾄﾞ(>>203)はあまり最適化されてないと考えられる） はっきりいって、これ以上見る気もしないしﾃｷﾄｰだが ｳﾝｺみたいな最適化で演算の順序が入れ替わったせいで、誤差が発生しているものと考えられる 212 名前：デフォルトの名無しさん mailto:sage [2018/09/16(日) 08:28:37.25 ID:zL1WUjLu.net] >>199以降、 すまん、入れ替わりになるかもしれんが後で確認する。 まず>>195その他について回答する。 >>195 SSEは /arch:SSE または /arch:SSE2 でないと出ないことになっており、勿論設定はしていない。 また、逆アセンブル結果では x87 命令のみであるのも確認している。 ただ今回の問題は、本当にReleaseビルドのバイナリを逆アセンブルしているか怪しい事だが。 >>196 小さい値なのは偶々だ。 辿って行ってそれが1回目にヒットする入力データだっただけのこと。 >>197 .NETの書式指定はググり難いが以下。 https://docs.microsoft.com/en-us/dotnet/standard/base-types/standard-numeric-format-strings 概ねprintfと同じで、自動的にやってくれるのが増えている。 30桁欲しければ以下。 Console::Write(String::Format("{0:E6}, {0:E30}\r\n",norm)); // 出力は 1.051355E-154, 1.051355436595307800000000000000E-154 なおdoubleは16桁な。(15.9=53*log(10)2) ただ当たり前だが、書式を変えたところで計算結果は変わらないし、 精度の問題には関係ない。 213 名前：デフォルトの名無しさん mailto:sage [2018/09/16(日) 08:29:03.14 ID:zL1WUjLu.net] >>195 FPU control registor については何故か安定した結果を得られていない。 インラインアセンブラは以下の通りで、 #pragma unmanaged inline void fpu_getcw(unsigned short* cw) { __asm{ fnstcw [cw]; } } #pragma managed これを norm = calc_norm_and_regulate( ... ) の直前/直後に配置して読み出し、 同様にコンソール出力すると、以下となる。 また、IDEで起動した場合は、「レジスタ」で見れる。 なお定義は以下の通り。 [9:8]に対し、 0x00 : 単精度(24bit) 0x01 : reserved 0x10 : 倍精度(53bit) 0x11 : 拡張倍精度(64bit) [11:10]に対し 0x00 : 最近値 0x01 : 切り捨て 0x10 : 切り上げ 0x11 : 0方向への切り捨て 214 名前：デフォルトの名無しさん mailto:sage [2018/09/16(日) 08:29:57.53 ID:zL1WUjLu.net] >>195 直後のみに配置： 0x027F (倍精度） ＝　Debug(IDE起動）のIDE内表示、Release(IDE起動）のIDE内表示、 0x03a5 (拡張倍精度)　＝　Debug(IDE起動）、Release(IDE起動）、 0x3fdc (拡張倍精度） ＝　Debug（コマンドプロンプト）、 0xf280, 0xf290, 0xf160, 0xf010等、不安定　＝　Release（コマンドプロンプト） 直前のみに配置： 直後のみと同じ結果。（つまり『何故か』安定している） Release（コマンドプロンプト）は不安定なのも同じ。 直前と直後に配置： 直前側は当然不安定になる。 直後側は「直後のみ」の結果と同じ。(Release（コマンドプロンプト）は不安定なのも同じ) 雰囲気からすると、IDE内表示は当てにならず、 命令自体は rdtsc と同じで非同期に実行されている雰囲気だが、 rdtsc命令の注意書きにある「シリアル化命令ではない」という但し書きが無く、状況は不明。 正直、正しく読み出せているか怪しい。（あてにならない） これらから推測すると、暫定的には以下。 拡張倍精度　＝　Debug(IDE起動）、Release(IDE起動）、Debug（コマンドプロンプト）、 不明　＝　Release（コマンドプロンプト） 以上が>>195その他に対する回答。 これから>>199その他について確認する。 215 名前：デフォルトの名無しさん mailto:sage [2018/09/16(日) 12:52:37.56 ID:haV9TZ8e.net] >>210 興味深い結果だ。 216 名前：デフォルトの名無しさん mailto:sage [2018/09/16(日) 13:21:26.16 ID:haV9TZ8e.net] >>210 >命令自体は rdtsc と同じで非同期に実行されている雰囲気だが、 >rdtsc命令の注意書きにある「シリアル化命令ではない」という但し書きが無く、状況は不明。 >正直、正しく読み出せているか怪しい。（あてにならない） インラインアセンブラを使わずに、 _control87(), _controlfp() : Get and set the floating-point control word. unsigned int _control87( unsigned int new, unsigned int mask ); unsigned int _controlfp( unsigned int new, unsigned int mask ); を使ってみたらどうなる？ 217 名前：デフォルトの名無しさん mailto:sage [2018/09/16(日) 13:31:15.50 ID:Q5j4SiHR.net] win32コンソールなら結果が同じ。 もう理由は分かったのに何が問題なんだ？こんなの何の影響もないだろう。 218 名前：デフォルトの名無しさん mailto:sage [2018/09/16(日) 13:33:14.92 ID:zL1WUjLu.net] >>199 再現実験ありがとう。 しかし色々問題がある。 1. 俺は起動方法による違いについてフォーカスしているが、 　君はRelease/Debugの違いにフォーカスしている。 2. VC++2008では再現しない。(VC++2010では再現する） 3. ソース改変しすぎ。それでは意味がない。 4. >>207の結論は間違い。 まず問題なのはソースの改変だ。 ループ回数を16回と決め打ちしたことで 8*2 に展開されている。 その結果、元のソース（俺が遭遇した状況）では発生しえないことが発生している。 これでは意味がない。 そして、君の結論は間違いだ。 × > ｳﾝｺみたいな最適化で演算の順序が入れ替わったせいで、誤差が発生しているものと考えられる 逆アセンブルを追えば分かるが、演算順序は入れ替わっていない。 原因は、Debugでは fld/fmul/fadd/fstp と毎回64bitに整形されるのに対し、 Releaseでは (fld/fmul/fadd)*8 + fstp と整形が8回に1回と減り、 8回は80bit(拡張倍精度）で演算されるからだ。 (こうなったのは君が16回ループ決め打ちコードに改変したから) ただしIDE上の fpu control registor の値は 相変わらず0x027F(倍精度）となっており、 IDEのこの表示が当てにならない事は分かる。 なおVC++2008では再現しなかった。 俺の環境では、16回決め打ちコードでも 8*2 に展開されず、Debugと同じコードだからだ。 勿論結果も同じだった。 219 名前：デフォルトの名無しさん mailto:sage [2018/09/16(日) 13:33:48.22 ID:zL1WUjLu.net] >>199 問題は、俺の環境で俺が提供したコード>>192だと、 同様に展開されないにも関わらず、『起動方法によって』結果が異なってしまう点だ。 俺の環境でのRelease/Debugの逆アセンブル結果のdiffは以下。 17c17 < 0000000c cmp dword ptr ds:[001C2E14h],0 --- > 0000000c cmp dword ptr ds:[00702E14h],0 19c19 < 00000015 call 68302BA9 --- > 00000015 call 683A5AB1 93c93 < 0000015a call FF6C3098 --- > 0000015a call FFCA57E8 98c98 < 0000016f push 0B5311Ch --- > 0000016f push 0D03188h 104,105c104,105 < 00000183 push 4F9D68h < 00000188 call FF6C30A4 --- > 00000183 push 2B71C0h > 00000188 call FFCA57F4 アドレスの変更だけであり、君の結果 「ループ回数を決め打ちしたことによりアンローリングされ、一部の演算がx87精度で計算される」には該当しない。 そして、この状況でも結果が異なってしまうことが問題なのだ。 君は君が勝手に新しく作り込んだ問題に対し、間違った結論でお茶を濁したにすぎない。 君が知っているFPU関連のことはこちらも知っている上で、質問している。 220 名前：デフォルトの名無しさん mailto:sage [2018/09/16(日) 13:53:17.67 ID:haV9TZ8e.net] >>214 なるほど、全く別の2つの理由で、精度が変わっている可能性（というより多分、確実）があると。 それは以下の2つ： 1. Debug版とRelease版では、最適化の結果、x87 FPU命令の使われ方が変わる。 　 x87では、メモリに書き戻さずに st(0)〜st(7)レジスタに入っている途中では、 　拡張倍精度の80BITで計算されるが、書き戻すとdoubleの場合でも64BITに丸め 　られる。なるべくメモリに書き戻さずにレジスタった方が高速なので、Release版 　では、80BITで計算される「期間」が多くなる。そのため、Debug版とRelese版では 　結果が僅かに違ってくる。 2. fpu control word が違っていて、st(0)〜st(7)に入っていても、計算が 　80BITか、64BIT、32BITのどれで計算されるか異なったり、丸め方が四捨五入 　か、正負二種類の方向の切り捨てなどが変わっている。 そして、IDEから起動した場合と、コマンドラインから起動した場合で結果が違う のは、「2.」の理由によるものではないかと。そして、実際に、インラインアセンブラ で fpu control word を読み取ってみると、不安定な値が読み出されたと。 221 名前：デフォルトの名無しさん [2018/09/16(日) 13:58:15.86 ID:haV9TZ8e.net] まず、 __asm{ 　fnstcw [cw]; } ではなく、_control87() を使ってみて欲しい。 インラインアセンブラは、独立した *.asm で書くより危険な場合が あるかも知れないので。特にC関数の引数、今の場合は、「cw」を インライン・アセンブラで用いた場合、正しいコードが出ているかどうか は注意が必要。 222 名前：デフォルトの名無しさん mailto:sage [2018/09/16(日) 14:05:28.41 ID:haV9TZ8e.net] >>209 よく見ると、それは、かなり複雑な事情が絡みそうなコード。 以下のようにした方が安心。なお、「cw」という短すぎ� 223 名前：驤�数名 も長年のプログラミング経験からすると、インラインアセンブラでは 怖い。また、 TTTT reg,引数名 や TTTT 引数名 は大丈夫でも、 TTTT reg,[引数名] や TTTT [引数名] は1命令では不可能な事をコンパイラに指示している事になるので ちょっと怖い。間接の間接、つまり、[[ebp+8]]みたいな事を要求 しているが、そんなオペランドが使えるアセンブリ命令はx86/x64 では存在しないので。 #pragma unmanaged inline void fpu_getcw(unsigned short *pCW) { 　__asm{ 　　mov edx,pCW 　　fnstcw [edx]; 　} } #pragma managed [] [ここ壊れてます] 224 名前：デフォルトの名無しさん mailto:sage [2018/09/16(日) 14:17:40.34 ID:haV9TZ8e.net] あ、後、インライン・アセンブラで実験する場合は、関数名の inline は 「取った」方がいい。つまり、以下の方が安心： #pragma unmanaged void fpu_getcw(unsigned short *pCW) { 　__asm{ 　　mov edx,pCW 　　fnstcw [edx]; 　} } #pragma managed 225 名前：デフォルトの名無しさん mailto:sage [2018/09/16(日) 15:23:18.78 ID:h8nMbN0G.net] また基本に戻るが、>>193で/MDになってるので /MTや/MTdでも発生するかしてみた方がいい 226 名前：デフォルトの名無しさん mailto:sage [2018/09/16(日) 15:37:02.44 ID:Q5j4SiHR.net] .netはx87コンテキストをすべて保持しませんって分かったんだからもう十分。 win32かx64にすれば解決。そもそも問題になる仕様バグじゃない。 227 名前：デフォルトの名無しさん mailto:sage [2018/09/16(日) 15:40:56.17 ID:haV9TZ8e.net] >>220 ホントだ。/MDだと、Runtime Library として DLL のものを使ってしまう。 これは、今回の現象に非常に重要な影響を与えているかも知れない。 228 名前：デフォルトの名無しさん mailto:sage [2018/09/16(日) 15:42:36.82 ID:haV9TZ8e.net] >>221 いや、まだまだ興味深い。 ここで終わらせずにもっと徹底して追及すべきだ。 229 名前：デフォルトの名無しさん mailto:sage [2018/09/16(日) 15:48:45.63 ID:zL1WUjLu.net] >>212 それはどうやらclrでは使えないらしい。 > These functions are ignored when you use /clr (Common Language Runtime Compilation) or /clr:pure to compile > because the common language runtime (CLR) only supports the default floating-point precision. > https://msdn.microsoft.com/en-us/library/e9b52ceh.aspx とはいえ無理矢理やってみた。警告は出るがコンパイルは通る。 結果は、どこに置いても、Debug/Releaseでも、常に 0x9001f が読み出される。 ただし、これは上記の仕様からして、当てにならない。 230 名前：デフォルトの名無しさん mailto:sage [2018/09/16(日) 15:49:44.99 ID:zL1WUjLu.net] >>219 218のコードで試してみた結果、209で言った不安定さはなくなり、 全てにおいて 0x027f が安定して読み出せるようになった。 ただしその過程で気づいたが、 IDEから起動した場合はReleaseビルドであっても、「未初期化のスタック値」も0x00が読み出せるようだ。 どうやらこれが原因の可能性が出てきた。（はっきり言って俺のバグだが） コードは以下の通りだが、 unsigned short fpu_cw, fpu_cw_after; // fpu_getcw(&fpu_cw); double norm = calc_norm_and_regulate(count, inputs, false); fpu_getcw(&fpu_cw_after); Console::Write(String::Format("{0:D}, 0x{0:x4}\r\n",fpu_cw)); Console::Write(String::Format("{0:D}, 0x{0:x4}\r\n",fpu_cw_after)); 読み出しと書き出し(Console::Write)を両方ともコメントアウトするのが面倒なので、 色々試す際、読み出しだけコメントアウトし、不定を表示させて脳内で省略していたのだが、 IDEから起動した場合はReleaseビルドであっても必ず0x0000が表示される事に気づいた。 上記『初期化していない』 fpu_cw を Releaseビルドをコマンドプロンプトから実行： 不定 ReleaseビルドをIDEから実行： 常に0x0000 となる。 実行前にあらかじめスタック領域を0fillでもしているのか？ まあこれに当たっているのなら確実に俺のバグだし、これなら辻褄は合ってしまうのだが。 231 名前：デフォルトの名無しさん mailto:sage [2018/09/16(日) 15:51:01.50 ID:zL1WUjLu.net] >>219 なお、逆アセンブルでコードバイトを表示させて確かめることは出来る。 正しいコードは出ている。(ただし不安定） inline void fpu_getcw(unsigned short* cw) { 00DA1540 55 push ebp 00DA1541 8B EC mov ebp,esp __asm{ fnstcw [cw]; 00DA1543 D9 7D 08 fnstcw word ptr [cw] } } 00DA1546 5D pop ebp 00DA1547 C3 ret fnstcwは D9 /7 で 7D なら [EBP+disp8] となり、 7D 08 は [EBP+08] となる。 つまりスタックポインタ+8の領域に書き戻せ、となる。 [ebp+0]は元のebpが入っているから、(pushしているので） [ebp+4]にcallの戻り値アドレス [ebp+8]にcw（第一引数）が入っていることになる。 これは正しいコードだ。 しかし再度試したが、確かに不安定だ。何故かは分からん。 inline取ってみても不安定のまま。 > そんなオペランドが使えるアセンブリ命令はx86/x64 > では存在しないので。 正直、/7の意味が分からないのだが、説明は > /digit − 0 から7 までの数字で、命令のModR/M バイトがr/m（レジスタまたはメモリ）オペランドだけを使用することを示す。 > reg フィールドには、命令のオペコードを拡張する数字が入っている。(Intelのマニュアルより） となっているのだが、これはどういう意味だ？ ModR/Mバイトが全部使えるとすると [ebp+disp8]出来ることになる。そしてそのコードは出ている。 ただし、動作は怪しいのも事実。 ModR/Mの一部しか使えない、ということか？ 232 名前：デフォルトの名無しさん mailto:sage [2018/09/16(日) 15:51:25.16 ID:zL1WUjLu.net] >>219 218のコードだと、 00381002 EC in al,dx __asm{ mov edx,pCW 00381003 8B 55 08 mov edx,dword ptr [pCW] fnstcw [edx]; 00381006 D9 3A fnstcw word ptr [edx] } } 00381008 5D pop ebp 00381009 C3 ret D9 3A ならまんま fnstcw [edx] だ。 理由は分からんがこちらだと安定しているので、結果としてはこのやり方が正しい。 233 名前：デフォルトの名無しさん mailto:sage [2018/09/16(日) 16:02:47.64 ID:haV9TZ8e.net] >>226 をを。やはり、ある意味ではVCが間違ったアセンブリコードを出していたよ。 それだと、 fnstcw [EBP+08] という意味になってしまって、 fnstcw pCW の意味になっている。つまり： pCW = control_word; あなたが、やりたいのは、 *pCW = control_word; だったのだから、アセンブリ・コードが間違ってる。 あなたが指示したのは、 fnstcw [pCW] だった。実際に生成されたコードは、 fnstcw pCW だった。 VC のインラインアセンブラは、エラーも出さずに間違ったコードを 出すことが証明された。 これと、精度が不安定な問題とは全く別ではあるけれど。 234 名前：デフォルトの名無しさん mailto:sage [2018/09/16(日) 16:06:05.39 ID:zL1WUjLu.net] すまん、間違いの修正 >>225 × > どうやらこれが原因の可能性が出てきた。（はっきり言って俺のバグだが） × > まあこれに当たっているのなら確実に俺のバグだし、これなら辻褄は合ってしまうのだが。 今回は俺はあくまで俺の本番コードのデバッグを念頭に置いていて、この発言だった。 ただし>>192の再現コードで『不定スタック領域』を掴んでいるわけもなく、 一応IDE起動とコマンドプロンプト起動での挙動の違いを再現出来ているわけだから、 これだけが問題ではないのも事実だ。 俺にとっては一つ新しい知見として、 ・IDEから起動した場合、スタックが初期化されるっぽい ということが分かった。とはいえOSは0fillしてから各プロセスにメモリを与えるので、実際は、 ・コマンドプロンプト起動ならmain前に設定した続きでそのまま実行、 ・IDE起動ならmain前に色々やって0fillして実行、 　或いはmain前に色々やることが多く、スタックが進み、(例えばデバッガをアタッチする為) 　結果的にOSが初期化済みの領域から始動 となって違いが発生するというところか。 235 名前：デフォルトの名無しさん mailto:sage [2018/09/16(日) 16:11:52.18 ID:haV9TZ8e.net] >>226 >正直、/7の意味が分からないのだが、 ModRM とは、 mod reg r/m 76 543 210 のようなオペランドを指しているのだけど、/7 は、regの部分を2進数の111、 10進数の「7」にするという意味。このタイプのマシン語は、 mod ttt r/m とも書かれる。tttの部分は、命令の主幹部分(ニモニック部分)によって変わる。 普通は、レジスタ番号を入れるところに、命令の種類を表す3BITの値を入れる 仕様になっている。 あなたがインラインアセンブラでVCに出させたかったコードは、意味的には、 fnstcw [[EBP+08]] なのだが、[ ] を二重にしたようなそんなx86/x64命令は存在しないので VC がエラーも出さずに勝手に一重の fnstcw [EBP+08] にしてしまった、という事。本当は、 mov edx,[EBP+08] fnstcw [edx] 236 名前： というコードにしなくてはならなかったのに、VCがある意味では間違った。 これが、単独の *.asm ではなく、VC の asm {・・・}　が危険な理由。 VC の asm は特に危険。 [] [ここ壊れてます] 237 名前：デフォルトの名無しさん mailto:sage [2018/09/16(日) 16:17:02.07 ID:haV9TZ8e.net] >>229 >ただし>>192の再現コードで『不定スタック領域』を掴んでいるわけもなく、 >一応IDE起動とコマンドプロンプト起動での挙動の違いを再現出来ているわけだから、 >これだけが問題ではないのも事実だ。 そうだよ。精度が変わるのはあなたの間違いではない。スタック領域が0クリア されようがれまいが、あなたのコード自体には特に不安定さはない。 非初期化領域を参照しているコードは見当たらないし。 238 名前：デフォルトの名無しさん mailto:sage [2018/09/16(日) 16:20:22.27 ID:haV9TZ8e.net] 逆アセンブラ結果を見てないで言うけど、もし、sqrt() が call文で関数呼び出し されているんだったら、そこで精度の違いが出てるかもしれない。 239 名前：デフォルトの名無しさん mailto:sage [2018/09/16(日) 16:23:20.75 ID:zL1WUjLu.net] >>228 なるほど、了解した。 つまり、>>210は全面的に間違いで、正しくは、 ・fpu control register は 0x027F で、IDEからも正しく読めている だな。 俺がやるべきだったのは fnstcw [[cw]] なのだと思うが、これはSyntaxErrorだ。 そして、こんな命令はないから、 []内に変数を書かず、レジスタ名にしろ、ということだったのだな。 全くもって了解だ。 VCの問題ではなくて、 俺が fnstcw [cw] と書いたのが間違いで、それをそのままコードにされてしまっただけだな。 正しく書けばSyntaxErrorだったのだし。 なお fnstcw [*cw] もSyntaxErrorだ。手動で一旦レジスタに移さないと駄目だな。 全くもって>>219のコードが正しい。 240 名前：デフォルトの名無しさん mailto:sage [2018/09/16(日) 16:35:08.59 ID:LrdaMWHl.net] >>233 >俺がやるべきだったのは fnstcw [[cw]] なのだと思うが、これはSyntaxErrorだ。 ちょっと違う。あなたはやるべきことをちゃんと正しく、 fnstcw [cw] と書いた。しかし、cw=[ebp+8]なので、これは、 fnstcw [[ebp+8]] という「意味」になる。でも、x86/x64のマシン語にはこんな[ ]を二重にした オペランドは存在しないので、VCが無断で勝手に[ ]を一重にして、 fnstcw [ebp+8] に改変してしまった。 **(ebp+8) = control_word; としなくてはならないのに、VCが勝手に、 *(ebp+8) = control_word; としたということ。 241 名前：デフォルトの名無しさん mailto:sage [2018/09/16(日) 16:36:04.13 ID:zL1WUjLu.net] >>230-231 了解だ。ありがとう。 >>232 その部分の逆アセンブラは以下の通り。 普通にcallされている。(行数オーバーなので切るが) ただし、 > そこで精度の違いが出てるかもしれない との繋がりがよくからない。 sqrt()でcallされると、スタックが改変される。おそらくデータ依存か？ なら未初期化のスタックを掴みに行っているコードが有ればバグる。 ただし今回の『再現コード』はこの限りではない。 （俺の本番コードはさておき） 242 名前：デフォルトの名無しさん mailto:sage [2018/09/16(日) 16:37:19.30 ID:zL1WUjLu.net] >>232 逆アセンブラ 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] 243 名前：デフォルトの名無しさん mailto:sage [2018/09/16(日) 16:37:34.94 ID:zL1WUjLu.net] >>232 逆アセンブラ（続き） 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] 244 名前：デフォルトの名無しさん mailto:sage [2018/09/16(日) 16:40:36.14 ID:LrdaMWHl.net] >>235 よく見ると、最小(?)の実験コードでは sqrt() が使われていなかった。 スマン。 245 名前：デフォルトの名無しさん mailto:sage [2018/09/16(日) 16:42:54.99 ID:zL1WUjLu.net] >>234 ああ、なるほど、了解。 246 名前：デフォルトの名無しさん mailto:sage [2018/09/16(日) 16:49:54.65 ID:zL1WUjLu.net] >>238 いや、俺が提供した>>192のソースなら使われてるぞ。 >>201のソースでは使われてないが。 ただまあ、彼(200)がsqrtを落としたのも分からなくはない。 誤差が生じる＝通常は桁落ちだから、この場合は当然積和部分が怪しい。 あらかじめ彼はそうなると分かっていてそれを落とし、予定調和的な結論にたどり着いてしまった。 それが彼の間違いだった、ということ。 俺は出来るだけ元のソースのままで追跡しようとしている。 元のソースの該当ケースと離れてしまっては意味がないから。 そして元ソースではsqrtを使っている。 247 名前：デフォルトの名無しさん mailto:sage [2018/09/16(日) 16:53:58.12 ID:/oSJzlqn.net] たぶん2008の最適化ミスだと思う。 static double norm = 0;// ←"static"を追加する にするとか、最適化オプションをいじると Release/コマンドプロンプトからの起動でも 0x1ff68ddfb62221ddになる 248 名前：デフォルトの名無しさん mailto:sage [2018/09/16(日) 16:54:10.71 ID:LrdaMWHl.net] >>238 ああ。また訂正。 sqrt()が使われていないのは、>>201, >>202, >>203, >>204 の場合で、 それは、ループ内にfprintf()を入れた場合と入れない場合とで、 x87 fpuレジスタのst(0)〜st(7)を使う「期間」が変わるために 80BITから 64BITへの書き戻し丸めの問題のために精度が変わっているだけだった。 一方、あなたが指摘した >>192 では、ちゃんと sqrt() 関数が使われていて、 それだと、IDEからの起動とコマンド・プロンプトからの起動とで、精度が変 わってくると。そして、その場合の逆アセンブル結果は >>236 のように sqrt() 関数がその場で x87 fpu の fsqrt 命令を使わずに、call 文によって 実際に本当のサブ・ルーチンを呼び出していると。 これはとても興味深い。そのサブ・ルーチンの中が、時と場合によって 精度が変わってくるような書き方をされている可能性が見えてきた。 249 名前：デフォルトの名無しさん mailto:sage [2018/09/16(日) 16:56:42.67 ID:LrdaMWHl.net] >>240 >いや、俺が提供した>>192のソースなら使われてるぞ。 > >>201のソースでは使われてないが。 了解。 問題を切り分けるため、sqrt() を使わなかった場合の Release版での、 IDE起動とコマンドｒライン起動の精度の違いを実験してみて欲しい。 250 名前：デフォルトの名無しさん mailto:sage [2018/09/16(日) 17:02:35.56 ID:LrdaMWHl.net] ちょっとしばらく、ここを離れる。 251 名前：デフォルトの名無しさん mailto:sage [2018/09/16(日) 17:22:07.21 ID:zL1WUjLu.net] >>241 現象確認した。こちらでも再現した。 逆アセンブルは、以下。(肝心のループ部分は次レス内） 正直、fld/fmul/fadd/fstpのループ部分は変わらず、 normのアドレスが [ebp-10h]（つまりローカル）から ds:[00A4AD40h]（つまりグローバル）に変わっただけであり、 これで結果が変わるのはかなり奇妙な気もするが、何か見落としがあるのかも。 >>241逆アセンブル(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 252 名前：デフォルトの名無しさん mailto:sage [2018/09/16(日) 17:22:29.32 ID:zL1WUjLu.net] >>241逆アセンブル（続き）(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] 253 名前：デフォルトの名無しさん mailto:sage [2018/09/16(日) 17:22:46.06 ID:zL1WUjLu.net] >>241逆アセンブル（続き）(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 254 名前：デフォルトの名無しさん mailto:sage [2018/09/16(日) 17:35:35.70 ID:zL1WUjLu.net] >>243 まだ異なった出力が得られた。 この意味では200がsqrtを外した判断は正しかった。 (彼はそこからさらにループ回数を固定してしまったのが間違いだった) 191ソースを以下に変更した。(sqrtをコメントアウト） ついでに Console::Write(String::Format("{0:E6}, {0:E30}\r\n",norm)); の出力も付けておく。 ソース： template<typename T> static double calc_norm_and_regulate(int num, T* r, bool regulate){ // <float> for debug. double norm = 0; for (int i=0;i<num;i++) norm += (double)r[i] * (double)r[i]; // norm = sqrt(norm); if (regulate) for (int i=0;i<num;i++) r[i] = (T)(r[i]/norm); return norm; } 結果：（Releaseビルド/コマンドプロンプトからの起動） 0.000000, 0x0007f2c44dfff8f2 1.105348E-308, 1.105348254058510600000000000000E-308 結果：（Releaseビルド/IDEからの起動、Debugビルドは起動方法によらずこちら) 0.000000, 0x0007f2c44dfff8f1 1.105348E-308, 1.105348254058510100000000000000E-308 >>244 了解。いずれにしても助かってる。 こちらも後30分くらいでちょっと離れる予定。 255 名前：デフォルトの名無しさん [2018/09/16(日) 18:30:34.19 ID:HF0YmRsW.net] >>213 ほんそれ 256 名前：デフォルトの名無しさん mailto:sage [2018/09/16(日) 20:54:27.39 ID:zL1WUjLu.net] >>241 さて再見したが、やはりstaticだけで直る理由は分からない。 なお、最適化ミスの場合は、逆アセンブラを読めば分かる。 今のところそれではない。 一応、>>192ソースのtemplate部の逆アセンブルを上げておく。（ただし重複するので頭のみ） 頭はこれ。続きが>>236,237。 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] 257 名前：デフォルトの名無しさん mailto:sage [2018/09/16(日) 21:25:23.90 ID:zL1WUjLu.net] >>220 >>222 /MTと/clrは同時に指定出来ないらしい。(error D8016) /MTdも同じく無理。 もう一つ /MDd ってのがあるから試してみた。 /MDdの結果： Releaseビルドでコマンドプロンプト起動の時のみ ****de、 ReleaseビルドでIDEからの起動だと ***dd。（Debugビルドは起動方法を問わずこっち） (/MDと全く挙動は同じ) これで有効な指摘については全て回答してるかな？ 見落としが有れば指摘よろしく。 （規制に引っかかったので遅くなってすまん） 今のところ、可能性があるのは以下か？ ・Releaseビルドをコマンドプロンプトから起動したときのみなぜか精度が高い 　(>>201から結果的に検出された。今のところ精度が高いときと同じ挙動をしている為） ・ReleaseビルドもIDEから起動すれば結果的にスタックが0初期化されている状態になっており、 　俺の本番プログラムに関してはここに当たるバグがある？(>>229) 　(ただしこれは>>192には該当しない） 258 名前：240 mailto:sage [2018/09/16(日) 21:43:04.47 ID:/oSJzlqn.net] .netの場合、デバッガ配下では(デバッグのため)違うコードを実行しているような気がする。 デバッガの逆アセンブル表示とかasm出力はあまり当てにならないような気もする。 ループ部分だけど、レジスタのみで処理するか、メモリを使用するかで精度が変わるのかも。 そもそも、どっちが正しいのかよくわからんけど... ループ部分の関数を#pragma unmanagedすると結果が変わるでそれが正しいのかも。 259 名前：240 mailto:sage [2018/09/16(日) 21:43:27.78 ID:/oSJzlqn.net] 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 260 名前：デフォルトの名無しさん mailto:sage [2018/09/16(日) 22:27:34.58 ID:zL1WUjLu.net] >>252 とりあえず落ち着け。一つずつ行こう。 > ループ部分の関数を#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 ec 261 名前：x,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) [] [ここ壊れてます] 262 名前：デフォルトの名無しさん mailto:sage [2018/09/16(日) 22:33:32.07 ID:zL1WUjLu.net] >>253 そちらの逆アセンブルは以下の違いが出てるだろ。 static版： fld/fmul/fadd/fstp 非static版： fld/fmul/faddp (fstpが無い） この非static版の場合、拡張倍精度(80bit)で演算されるから精度が高いことになり、 static版との演算結果に違いが出るのも仕様通りなんだよ。(これは>>201と同じ間違い) 一応、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でどう書くのかは知らん） だから>>253の場合の誤差なら、仕様通りなんだよ。(片方が倍精度、もう片方は拡張倍精度） ただし、>>192は逆アセンブル(>>236)を見る限りそれに該当しないし、(両方とも倍精度） 今回の俺の上記逆アセンブル（>>254、中身は君の指摘通りunmanagedにしただけ）も該当しない。(両方とも倍精度） そして253は何故か直ってしまった。 263 名前：デフォルトの名無しさん mailto:sage [2018/09/16(日) 22:34:16.44 ID:zL1WUjLu.net] >>253 > .netの場合、デバッガ配下では(デバッグのため)違うコードを実行しているような気がする。 > デバッガの逆アセンブル表示とかasm出力はあまり当てにならないような気もする。 これは俺も相当疑っているのだが、今のところ尻尾を掴めない。 ILspyだっけ？外部の逆アセンブルツール使えばチェック出来るのかな？ いずれにしても、>>252の指摘 ・unmanagedにすれば直る のも事実だし、逆アセンブルを見る限り、これを説明出来る理由もないのも事実。 264 名前：240 mailto:sage [2018/09/16(日) 23:13:48.62 ID:/oSJzlqn.net] > だから>>253の場合の誤差なら、仕様通りなんだよ。(片方が倍精度、もう片方は拡張倍精度） そうなの？　これが仕様通りならstatic版での違いは仕様通りということになる。 252はRelease版をコンソールで実行したときの逆アセンブル結果。 よって、Release版をコンソールで実行したときのみ(たまたま)レジスタ(80ビット)での演算になるので、 計算結果が変わるのはやむを得ないという結論になるのだが... ちなみに、235はDebugモードでコンパイルし、デバッガ配下の逆アセンブル結果でしょ。 265 名前：デフォルトの名無しさん mailto:sage [2018/09/16(日) 23:24:30.83 ID:zL1WUjLu.net] >>257 > 252はRelease版をコンソールで実行したときの逆アセンブル結果。 それはどうやって得たの？俺はそれが出来ないから困ってる。 > ちなみに、235はDebugモードでコンパイルし、デバッガ配下の逆アセンブル結果でしょ。 235は、IDE上でReleaseモードでF5で起動し、ブレークポイントを当てて止めて逆アセンブルした結果。 俺が貼ってる逆アセンブル結果は全てこの方法で、IDEで表示されているもの。 だからIDEの表示がおかしかったら話が全部おかしくなる。 君がIDEから独立して逆アセンブル出来ているのなら、その方法を知りたい。 こちらでも試す。 なおILSpy、グダグダ言わずに試してみたが、 当たり前だがmanaged code だとILが出る(x86ではない）ので、 俺って根本的に間違ってたかも？ 今までx86のアセンブラで議論してたけど、これって .NET アプリには同梱されていないというオチ？ (まあその場合は君がやっている外部逆アセンブルが単純には出来ないはずなのだが） 266 名前：240 mailto:sage [2018/09/16(日) 23:35:36.48 ID:/oSJzlqn.net] >> 252はRelease版をコンソールで実行したときの逆アセンブル結果。 >それはどうやって得たの？俺はそれが出来ないから困ってる。 calc_norm_and_regulateの次の行に System::Diagnostics::Debugger::Launch(); を入れてコンソールから実行すると just in time デバッグできるので、デバッガを選んだ後、 Visual Studioの 呼び出し履歴から calc_norm_and_regulate を探して移動する >> ちなみに、235はDebugモードでコンパイルし、デバッガ配下の逆アセンブル結果でしょ。 >236は、IDE上でReleaseモードでF5で起動し、ブレークポイントを当てて止めて逆アセンブルした結果。 あれ？ だとすると最適化していないのでは？ こちらの結果と違うのだが。 267 名前：デフォルトの名無しさん [2018/09/16(日) 23:49:32.68 ID:SOVIz+sV.net] とりあえず、>>199に書いてあるとおり、元のｺｰﾄﾞで再現はしてるからな その再現したｺｰﾄﾞと逆ｱｾﾝﾌﾞﾙした結果はあげとく �B-1 ﾃﾞﾌｫﾙﾄ設定(Release)　【ｺｰﾄﾞ】 #include "stdafx.h" #include <stdio.h> #include <stdint.h> #include <math.h> using namespace System; template<typename T> static double calc_norm_and_regulate(int num, T* r, bool regulate){ // <float> for debug. 　　　double norm = 0; 　　　for (int i = 0; i < num; i++) { 　　　　　　norm += (double)r[i] * (double)r[i]; //　　　　　　fprintf(stdout, "[0]0x%016llX:%.19lg\n", *(uint64_t*)&norm, norm); 　　　} 　　　fprintf(stdout, "[1]0x%016llX:%.19lg\n", *(uint64_t*)&norm, norm); 　　　norm = sqrt(norm); 　　　if (regulate) 　　　　　　for (int i=0;i<num;i++) 　　　　　　　　　r[i] = (T)(r[i]/norm); 　　　return norm; } 268 名前：デフォルトの名無しさん [2018/09/16(日) 23:50:40.20 ID:SOVIz+sV.net] �B-1 ﾃﾞﾌｫﾙﾄ設定(Release)　【ｺｰﾄﾞ】（その2） int main(array<System::String ^> ^args) { 　　　int count = 16; 　　　__int64 inputs_hex[16] = { 　　　　　　0x0000000000000000, 0x0000000000000000, 0x0000000000000000, 0x1fedb1530240aa54, 　　　　　　0x0000000000000000, 0x0000000000000000, 0x0000000000000000, 0x1ff0af0d95025bc3, 　　　　　　0x1fc9353df6af376b, 0x0000000000000000, 0x0000000000000000, 0x0000000000000000, 　　　　　　0x0000000000000000, 0x0000000000000000, 0x0000000000000000, 0x0000000000000000 　　　}; 　　　double* inputs = (double*)inputs_hex; 　　　double norm = calc_norm_and_regulate(count, inputs, false); 　　　fprintf(stdout, "[2]0x%016llX:%.19lg\n", *(uint64_t*)&norm, norm); 　　　// Console::Write(String::Format("{0:F6}, 0x{1:x16}\r\n",norm, *(__int64*)&norm)); 　　　// Release build 　　　// 0.000000, 0x1ff68ddfb62221dd from IDE 　　　// 0.000000, 0x1ff68ddfb62221de from command prompt 　　　return 0; } 269 名前：デフォルトの名無しさん [2018/09/16(日) 23:51:14.21 ID:SOVIz+sV.net] �B-1 ﾃﾞﾌｫﾙﾄ設定（Release）　【実行結果】 ↓このｺｰﾄﾞの逆ｱｾﾝﾌﾞﾙｺｰﾄﾞ https://ideone.com/Gf4qUQ [1]0x0007F2C44DFFF8F2:1.1053482540585106e-308 [2]0x1FF68DDFB62221DE:1.051355436595308e-154 270 名前：デフォルトの名無しさん [2018/09/16(日) 23:54:13.85 ID:SOVIz+sV.net] で、>>200に書いてあるとおり↓Debugﾋﾞﾙﾄﾞと同じ結果が再現された > ↓この下に（ﾙｰﾌﾟ内に）fprintf文を入れるだけで > 　ReleaseﾋﾞﾙﾄﾞとDebugﾋﾞﾙﾄﾞが同じ実行結果になることが確認できた > norm += (double)r[i] * (double)r[i]; �B-2 ﾃﾞﾌｫﾙﾄ設定(Release)　【ｺｰﾄﾞ】（その1） #include "stdafx.h" #include <stdio.h> #include <stdint.h> #include <math.h> using namespace System; template<typename T> static double calc_norm_and_regulate(int num, T* r, bool regulate){ // <float> for debug. 　　　double norm = 0; 　　　for (int i = 0; i < num; i++) { 　　　　　　norm += (double)r[i] * (double)r[i]; 　　　　　　fprintf(stdout, "[0]0x%016llX:%.19lg\n", *(uint64_t*)&norm, norm); 　　　} 　　　fprintf(stdout, "[1]0x%016llX:%.19lg\n", *(uint64_t*)&norm, norm); 　　　norm = sqrt(norm); 　　　if (regulate) 　　　　　　for (int i=0;i<num;i++) 　　　　　　　　　r[i] = (T)(r[i]/norm); 　　　return norm; } 271 名前：デフォルトの名無しさん [2018/09/16(日) 23:54:58.60 ID:SOVIz+sV.net] �B-2 ﾃﾞﾌｫﾙﾄ設定(Release)　【ｺｰﾄﾞ】（その2） int main(array<System::String ^> ^args) { 　　　int count = 16; 　　　__int64 inputs_hex[16] = { 　　　　　　0x0000000000000000, 0x0000000000000000, 0x0000000000000000, 0x1fedb1530240aa54, 　　　　　　0x0000000000000000, 0x0000000000000000, 0x0000000000000000, 0x1ff0af0d95025bc3, 　　　　　　0x1fc9353df6af376b, 0x0000000000000000, 0x0000000000000000, 0x0000000000000000, 　　　　　　0x0000000000000000, 0x0000000000000000, 0x0000000000000000, 0x0000000000000000 　　　}; 　　　double* inputs = (double*)inputs_hex; 　　　double norm = calc_norm_and_regulate(count, inputs, false); 　　　fprintf(stdout, "[2]0x%016llX:%.19lg\n", *(uint64_t*)&norm, norm); 　　　// Console::Write(String::Format("{0:F6}, 0x{1:x16}\r\n",norm, *(__int64*)&norm)); 　　　// Release build 　　　// 0.000000, 0x1ff68ddfb62221dd from IDE 　　　// 0.000000, 0x1ff68ddfb62221de from command prompt 　　　return 0; } 272 名前：デフォルトの名無しさん [2018/09/16(日) 23:56:17.88 ID:SOVIz+sV.net] �B-2 ﾃﾞﾌｫﾙﾄ設定（Release）　【実行結果】 ↓このｺｰﾄﾞの逆ｱｾﾝﾌﾞﾙｺｰﾄﾞ https://ideone.com/E3Nxt8 ↓実行結果を書き込めないからこっちに書き込んどいた https://ideone.com/1cky6N [1]0x0007F2C44DFFF8F1:1.1053482540585101e-308 [2]0x1FF68DDFB62221DD:1.0513554365953078e-154 以上だ 273 名前：デフォルトの名無しさん mailto:sage [2018/09/16(日) 23:58:42.66 ID:zL1WUjLu.net] >>259 おお、そのやり方は知らなかった。大変助かった。ありがとう。 で、結果だが、>>253とは微妙に違うが、確かに拡張倍精度で計算されている。 逆アセンブル結果は、以下。 0000000e D9 EE fldz for (int i=0;i<num;i++) norm += (double)r[i] * (double)r[i]; 00000010 33 C9 xor ecx,ecx 00000012 EB 01 jmp 00000015 00000014 41 inc ecx 00000015 3B CE cmp ecx,esi 00000017 7D 0B jge 00000024 00000019 DD 04 CF fld qword ptr [edi+ecx*8] 0000001c D9 C0 fld st(0) 0000001e DE C9 fmulp st(1),st 00000020 DE C1 faddp st(1),st 00000022 EB F0 jmp 00000014 norm = sqrt(norm); 00000024 83 EC 08 sub esp,8 00000027 DD 1C 24 fstp qword ptr [esp] 0000002a E8 49 7C F2 FF call FFF27C78 とにかく、Releaseビルドをコンソールから起動した場合は拡張倍精度になってるのは分かった。 なら、ReleaseビルドをIDEから起動した場合は何を起動してるんだこれは？ Debugビルドとも微妙にアドレス等が違うんだが。 とはいえ、これは「そもそも色々間違っている」可能性が出てきたので、もう一度全体を見直す。 明日（だけで済むとも思えないが）確認し、整理してまた投稿する。 とにかくありがとう。これはだいぶインパクトがある。（はず) 274 名前：デフォルトの名無しさん mailto:sage [2018/09/17(月) ] [ここ壊れてます] 275 名前： 00:49:10.27 ID:gBxGzBbQ.net mailto: >>266 ideとconsoleで微妙に違うバージョンのvcruntimeがロードされてるとか？ process monitorで何が実行されて何がロードされてるか調べるとか？ あとはwindbgから起動したら、ideとconsoleのどちらの結果と一致するのか気になる。 [] [ここ壊れてます] 276 名前：デフォルトの名無しさん mailto:sage [2018/09/17(月) 01:06:37.11 ID:+dwRu2dr.net] >>262 だからそれは>>201と同じなんだよ。 その逆アセンブルでいうと、以下部分がメモリに出力されず、拡張倍精度で動作してるだろ。 00000281 fld qword ptr [ebp+FFFFFF14h] 00000287 fmul st,st(0) 00000289 fadd qword ptr [ebp+FFFFFF70h] 0000028f fld qword ptr [ebp+FFFFFF1Ch] 00000295 fmul st,st(0) 00000297 faddp st(1),st 00000299 fld qword ptr [ebp+FFFFFF24h] 0000029f fmul st,st(0) 000002a1 faddp st(1),st 000002a3 fld qword ptr [ebp+FFFFFF2Ch] 000002a9 fmul st,st(0) 000002ab faddp st(1),st 000002ad fld qword ptr [ebp+FFFFFF34h] 000002b3 fmul st,st(0) 000002b5 faddp st(1),st 000002b7 fld qword ptr [ebp+FFFFFF3Ch] 000002bd fmul st,st(0) 000002bf faddp st(1),st 000002c1 fld qword ptr [ebp+FFFFFF44h] 000002c7 fmul st,st(0) 000002c9 faddp st(1),st 000002cb fld qword ptr [ebp+FFFFFF4Ch] 000002d1 fmul st,st(0) 000002d3 faddp st(1),st 000002d5 fstp qword ptr [ebp+FFFFFF70h] 277 名前：デフォルトの名無しさん mailto:sage [2018/09/17(月) 01:07:41.74 ID:+dwRu2dr.net] >>262（続き) これは少なくとも「ループ回数が8の倍数である」事がコンパイラに見えないと出来ない最適化だ。 そうでなければ、例えばループ回数が6回や14回の時に、 最初の１回だけ 0299 に飛び込んで始める（頭の2回をスキップする）コードが必要になるが、 それは出てないだろ。 (そもそもこのアンローリングがx86的に意味があるのかも疑問だが） 一般的に、可変回数ループを展開すると、必ず上記の端切れ処理（キリが良くないときの処理）が必要になる。 だから「可変」だと確定しているのなら普通は展開しない。 つまり、一般的には、別関数でループ回数が引数で与えられてたら、その最適化はかからない。 今回ヒットするデータが偶々16回ループだっただけで、 俺の本番コードは可変で用いている為、 このようなアンローリングはされてない。（と思ってる。ただし265の通りもう一度確認必要） こちらではデフォでそこまで最適化がかからないのでそちらの状況はよく分からないが、 その最適化がかかってるのなら、俺の本番コードとは違う状況で動いていることになる。(はず） だから俺はその最適化がかからない範囲で議論してる。 それが俺が最初から別関数にしてた理由。（というか元のコードが別関数だからだが） てゆうかマジでそれデフォのままか？ もしかして俺の環境がおかしくて、全く最適化されてないコードが出てる？ 278 名前：247 mailto:sage [2018/09/17(月) 01:12:11.55 ID:yaPtorLJ.net] 今、戻った。 席を離れて思ったが、多分、C++/CLI の場合、IDEでコンパイル後、 exeファイルになっても、unmanaged コード部分以外は、本質的には .NETの共通中間言語(CIL、MSIL)になるだけで、x86のマシン語には なってないと思う。一方、昔ながらの本当のC++では、exeの中身は、 本質的には、x86のマシン語が入っていた。 C++/CLIの場合、作成されたexeは、起動時に、中間言語がx86の マシン語に直されてから実行される。その際、最適化の関係で、 どう直されるかが必ずしも一定していない。だから、IDEからの起動と、 コマンドラインからの起動で異なったx86命令に直されてしまう のかも知れない。 出来たx86コードが、良く最適化された場合、x87 fpu命令を長きに 渡ってメモリに書き戻さずに、st(0)〜st(7)のレジスタ上で演算し 続けられることになる。その場合、80BITの精度で計算される 期間が長くなる。 一方、x87 fpu命令に関する最適化が少し悪いか、悪くなくても命令の 使い方が変わってしまうと、80BIT(拡張double)から64BIT(double)へ 直されるタイミングや回数が変わる。 この両者では結果の精度が変わってくる。前者の方が後者より僅か に精度が高くなる。 279 名前：243=269 mailto:sage [2018/09/17(月) 01:14:28.51 ID:yaPtorLJ.net] >>270 名前欄は、「247」ではなく「243」の間違いだった。 280 名前：デフォルトの名無しさん mailto:sage [2018/09/17(月) 01:24:06.04 ID:yaPtorLJ.net] >>258 >なおILSpy、グダグダ言わずに試してみたが、 >当たり前だがmanaged code だとILが出る(x86ではない）ので、 >俺って根本的に間違ってたかも？ >今までx86のアセンブラで議論してたけど、これって .NET アプリには同梱されていないというオチ？ やはり、managed code部分は、x86命令では無く、ILにコンパイルされていて、 普通のC++とは違ってたんだ。 281 名前：デフォルトの名無しさん mailto:sage [2018/09/17(月) 01:30:31.16 ID:yaPtorLJ.net] >calc_norm_and_regulateをunmanaged関数にすると、違いはなくなる。 >(Releaseビルドの｀をコマンドプロンプトで起動した際にも、****ddの結果となる） やはり。 unmanaged関数の場合は、CIL(MSIL)ではなく、exeの段階で既に x86マシン語に直されたものが格納されるんだろう。だとすると、起動方法に 関係なく、少なくともその部分に関しては、x87 fpu命令の使われ方が 全く同じになる。callしたsqrt()関数の中は除いて。 282 名前：デフォルトの名無しさん [2018/09/17(月) 01:37:20.55 ID:dj7qSZnZ.net] >>260-261のｺｰﾄﾞも>>263-264 も同じ条件でﾃﾞﾌｫﾙﾄでｺﾝﾊﾟｲﾙしてる ｺﾚは間違いない 一行ｺﾒﾝﾄはずしてｺﾝﾊﾟｲﾙしなおすだけだからな で、>>262、>>265みたいな結果になる >>265の実行結果はDebugﾋﾞﾙﾄﾞとまったく同じになる そのまんま 283 名前：デフォルトの名無しさん mailto:sage [2018/09/17(月) 02:03:16.26 ID:yaPtorLJ.net] はっきり書いてある。managed code は、起動時にJITコンパイルされる。 だから、どんなマシン語に置き換わるかが、コンパイルしただけでは まだ完全には決定されてない。 https://en.wikipedia.org/wiki/Managed_code Drawbacks include slower startup speed (the managed code must be JIT compiled by the VM) and generally increased use of system resources on any machine that is executing the code. managed code は、VMによって、JIT コンパイルされないといけないので、 起動速度が遅くなり、かつ、一般的に、システム・リソースの使用が増える。 284 名前：デフォルトの名無しさん [2018/09/17(月) 02:49:02.18 ID:dj7qSZnZ.net] 普通のｺﾝｿｰﾙｱﾌﾟﾘｹｰｼｮﾝ(CLRじゃないほう)でも 同じｺｰﾄﾞで普通にまったく同じように再現するﾜｹだが ﾎﾝﾄなｷﾐﾗはなにを頭悪いことばっかりいってんの？ 実行結果だけははっといてやるが 285 名前：デフォルトの名無しさん [2018/09/17(月) 02:50:12.89 ID:dj7qSZnZ.net] �C-1 ﾃﾞﾌｫﾙﾄ設定(Release)　【ｺｰﾄﾞ】（その1） #include "stdafx.h" #include <stdio.h> #include <stdint.h> #include <math.h> template<typename T> static double calc_norm_and_regulate(int num, T* r, bool regulate){ // <float> for debug. 　　　double norm = 0; 　　　for (int i = 0; i < num; i++) { 　　　　　　norm += (double)r[i] * (double)r[i]; //　　　　　　fprintf(stdout, "[0]0x%016llX:%.19lg\n", *(uint64_t*)&norm, norm); 　　　} 　　　fprintf(stdout, "[1]0x%016llX:%.19lg\n", *(uint64_t*)&norm, norm); 　　　norm = sqrt(norm); 　　　if (regulate) 　　　　　　for (int i=0;i<num;i++) 　　　　　　　　　r[i] = (T)(r[i]/norm); 　　　return norm; } 286 名前：デフォルトの名無しさん [2018/09/17(月) 02:50:37.91 ID:dj7qSZnZ.net] �C-1 ﾃﾞﾌｫﾙﾄ設定(Release)　【ｺｰﾄﾞ】（その2） int _tmain(int argc, _TCHAR* argv[]) { 　　　int count = 16; 　　　__int64 inputs_hex[16] = { 　　　　　　0x0000000000000000, 0x0000000000000000, 0x0000000000000000, 0x1fedb1530240aa54, 　　　　　　0x0000000000000000, 0x0000000000000000, 0x0000000000000000, 0x1ff0af0d95025bc3, 　　　　　　0x1fc9353df6af376b, 0x0000000000000000, 0x0000000000000000, 0x0000000000000000, 　　　　　　0x0000000000000000, 0x0000000000000000, 0x0000000000000000, 0x0000000000000000 　　　}; 　　　double* inputs = (double*)inputs_hex; 　　　double norm = calc_norm_and_regulate(count, inputs, false); 　　　fprintf(stdout, "[2]0x%016llX:%.19lg\n", *(uint64_t*)&norm, norm); 　　　return 0; } 287 名前：デフォルトの名無しさん [2018/09/17(月) 02:51:12.35 ID:dj7qSZnZ.net] �C-1 ﾃﾞﾌｫﾙﾄ設定（Release）　【実行結果】 ↓このｺｰﾄﾞの逆ｱｾﾝﾌﾞﾙｺｰﾄﾞ https://ideone.com/Dqqn6J [1]0x0007F2C44DFFF8F2:1.1053482540585106e-308 [2]0x1FF68DDFB62221DE:1.051355436595308e-154 288 名前：デフォルトの名無しさん [2018/09/17(月) 02:51:53.99 ID:dj7qSZnZ.net] �C-2 ﾃﾞﾌｫﾙﾄ設定(Release)　【ｺｰﾄﾞ】（その1） #include "stdafx.h" #include <stdio.h> #include <stdint.h> #include <math.h> template<typename T> static double calc_norm_and_regulate(int num, T* r, bool regulate){ // <float> for debug. 　　　double norm = 0; 　　　for (int i = 0; i < num; i++) { 　　　　　　norm += (double)r[i] * (double)r[i]; //　　　　　　fprintf(stdout, "[0]0x%016llX:%.19lg\n", *(uint64_t*)&norm, norm); 　　　} 　　　fprintf(stdout, "[1]0x%016llX:%.19lg\n", *(uint64_t*)&norm, norm); 　　　norm = sqrt(norm); 　　　if (regulate) 　　　　　　for (int i=0;i<num;i++) 　　　　　　　　　r[i] = (T)(r[i]/norm); 　　　return norm; } 289 名前：デフォルトの名無しさん [2018/09/17(月) 02:52:10.55 ID:dj7qSZnZ.net] �C-2 ﾃﾞﾌｫﾙﾄ設定(Release)　【ｺｰﾄﾞ】（その2） int _tmain(int argc, _TCHAR* argv[]) { 　　　int count = 16; 　　　__int64 inputs_hex[16] = { 　　　　　　0x0000000000000000, 0x0000000000000000, 0x0000000000000000, 0x1fedb1530240aa54, 　　　　　　0x0000000000000000, 0x0000000000000000, 0x0000000000000000, 0x1ff0af0d95025bc3, 　　　　　　0x1fc9353df6af376b, 0x0000000000000000, 0x0000000000000000, 0x0000000000000000, 　　　　　　0x0000000000000000, 0x0000000000000000, 0x0000000000000000, 0x0000000000000000 　　　}; 　　　double* inputs = (double*)inputs_hex; 　　　double norm = calc_norm_and_regulate(count, inputs, false); 　　　fprintf(stdout, "[2]0x%016llX:%.19lg\n", *(uint64_t*)&norm, norm); 　　　return 0; } 290 名前：デフォルトの名無しさん [2018/09/17(月) 02:52:25.20 ID:dj7qSZnZ.net] �C-2 ﾃﾞﾌｫﾙﾄ設定（Release）　【実行結果】 ↓このｺｰﾄﾞの逆ｱｾﾝﾌﾞﾙｺｰﾄﾞ https://ideone.com/5OaUe6 ↓実行結果を書き込めないからこっちに書き込んどいた https://ideone.com/FFW0P0 [1]0x0007F2C44DFFF8F1:1.1053482540585101e-308 [2]0x1FF68DDFB62221DD:1.0513554365953078e-154 291 名前：デフォルトの名無しさん [2018/09/17(月) 02:53:45.62 ID:dj7qSZnZ.net] �C-2 ﾃﾞﾌｫﾙﾄ設定(Release)　【ｺｰﾄﾞ】（その1） #include "stdafx.h" #include <stdio.h> #include <stdint.h> #include <math.h> template<typename T> static double calc_norm_and_regulate(int num, T* r, bool regulate){ // <float> for debug. 　　　double norm = 0; 　　　for (int i = 0; i < num; i++) { 　　　　　　norm += (double)r[i] * (double)r[i]; 　　　　　　fprintf(stdout, "[0]0x%016llX:%.19lg\n", *(uint64_t*)&norm, norm); 　　　} 　　　fprintf(stdout, "[1]0x%016llX:%.19lg\n", *(uint64_t*)&norm, norm); 　　　norm = sqrt(norm); 　　　if (regulate) 　　　　　　for (int i=0;i<num;i++) 　　　　　　　　　r[i] = (T)(r[i]/norm); 　　　return norm; } 292 名前：デフォルトの名無しさん [2018/09/17(月) 02:54:34.82 ID:dj7qSZnZ.net] （正）>>283 （誤）>>280 293 名前：デフォルトの名無しさん [2018/09/17(月) 02:57:46.48 ID:dj7qSZnZ.net] >>277-278のｺｰﾄﾞも>>283,281のｺｰﾄﾞも も同じ条件でﾃﾞﾌｫﾙﾄでｺﾝﾊﾟｲﾙしてる 一行ｺﾒﾝﾄはずしてｺﾝﾊﾟｲﾙしなおすだけだからな で、>>279、>>282みたいな結果になる >>282の実行結果はDebugﾋﾞﾙﾄﾞとまったく同じになる そのまんま CLRのとき（>>274）の動作とまったく同じ ｷﾐﾗはずーっとなにやってるわけ？ 294 名前：デフォルトの名無しさん mailto:sage [2018/09/17(月) 09:44:34.76 ID:yu1Dprt2.net] >>285 同一のreleaseをコンソールで実行するかデバッガで実行するかで結果が異なるのはなぜだろう。 という話をしていたのであって、 debug/releaseで別の結果になることを問題にしているのではないです。 295 名前：デフォルトの名無しさん mailto:sage [2018/09/17(月) 10:25:22.39 ID:+dwRu2dr.net] >>192がコンソール起動とIDE起動で挙動が異なる理由は分かりました。 ありがとう。 結論はつまり以下だ。 > JIT の最適化とデバッグ（抜粋） > マネージ アプリケーションをデバッグするとき、Visual Studio では、既定で、 > ジャスト イン タイム (JIT: Just-In-Time) コードの最適化が省略されています。 > 最適化されたコードをデバッグするのは困難であるため、 > 最適化されたコードで発生するバグが、非最適化バージョンでは再現しないときにのみお勧めします。 > JIT 最適化は、Visual Studio の [モジュールの読み込み中に JIT 最適化を抑制する] オプションで制御されます。 > 実行中のプロセスにアタッチする場合、既に読み込まれ、JIT でコンパイルされ、 > 最適化されているコードが含まれることがあります。 > このようなコードの場合、[モジュールの読み込み中に JIT 最適化を抑制する] オプションの影響はありません。 > https://msdn.microsoft.com/ja-jp/library/ms241594.aspx 確かにこのオプションで直った。 296 名前：デフォルトの名無しさん mailto:sage [2018/09/17(月) 10:25:48.51 ID:+dwRu2dr.net] その他諸々、話を整理すると、以下となる。（ソースは>>192参照） 1. managedコードではMSILが出力され、x86コードは含まれていない。 2. 起動時、MSILはJITされ、x86コードに落とされる。 3. このため、mainの1行目でブレークポイントで止め、calc_norm_and_regulateの逆アセンブルを見ようとしても、 　IDE上で「逆アセンブルを表示できません。式がまだネイティブ マシン コードに翻訳されていません。」と出る。 　これはmainの1行目に System::Diagnostics::Debugger::Launch(); を入れたときも同様。 4. そしてこのJITに関して、上記IDE中の 『[モジュールの読み込み中に JIT 最適化を抑制する] オプション』 が効いてくる。 　規定ではオフ、つまり、ReleaseビルドでもIDE起動ならJIT最適化は抑制される。 　これがfld/fmul/fadd/fstpのループコードになる理由。 　これをオンにすれば、確かにReleaseビルドIDE起動でも、 　fld/fmul/faddのループコードとなり、コマンドプロンプト起動と同じ結果になることは確認した。 5. 上記では表現が微妙だが、JIT最適化をするかどうかは読み込まれるときに決まるらしい。 　したがって、Releaseビルドを起動後にアタッチした場合は通常通り最適化され、 　IDEからReleaseビルドを起動した場合は『既定では』最適化が抑制されてしまう。 　これがIDE起動とコマンドプロンプト起動で挙動が異なった原因。 　上記、『[モジュールの読み込み中に JIT 最適化を抑制する]』のチェックを外せば、直った。 6. おそらくこのオプションはソリューション毎ではなく、IDEのインストール毎なんだと思う。 　（ソリューション毎のオプションはプロジェクトのプロパティにあり、場所が違う) 　だからその人の環境によっては最初からオフにしている人がいたかも？ 　これが再現実験をしてくれた人たちと微妙に結果が異なったりした原因か？ これで>>192についての疑問は解消した。（はず） 俺の本番コードについては再度確認し、また報告する。 297 名前：デフォルトの名無しさん [2018/09/17(月) 10:28:21.54 ID:dj7qSZnZ.net] 同じﾘﾘｰｽﾋﾞﾙﾄﾞで 結果がかわってんのになにいってんの？ fprintf入れるだけで ﾘﾘｰｽﾋﾞﾙﾄﾞで結果が変わる 298 名前：デフォルトの名無しさん [2018/09/17(月) 10:31:16.66 ID:dj7qSZnZ.net] やはり低学歴知恵遅れは 結果が意味するところが分かってないわ。。。 299 名前：デフォルトの名無しさん [2018/09/17(月) 10:32:49.39 ID:dj7qSZnZ.net] ちなみにｵﾚがあげた結果は すべてﾘﾘｰｽﾋﾞﾙﾄﾞの結果だからな ﾃﾞﾊﾞｯｸﾞﾋﾞﾙﾄﾞの結果なんかあげても 意味ないからな 300 名前：デフォルトの名無しさん [2018/09/17(月) 10:38:00.02 ID:dj7qSZnZ.net] CLRのｹｰｽもCLRでない普通のexeのｹｰｽでも 結果はまったく同じだからな しかもすべてﾘﾘｰｽﾋﾞﾙﾄﾞで おきてる誤差までぴったり一致してる 301 名前：デフォルトの名無しさん mailto:sage [2018/09/17(月) 10:52:24.35 ID:+dwRu2dr.net] >>289 お前は相変わらず理解してないな。 80bit（拡張倍精度）と64bit(倍精度)の演算で桁落ちが異なり、結果が異なるのは当然なんだよ。 問題は同じバイナリの癖に何故起動方法によって異なるのか？だったんだ。 理由はMSILだからだ。 MSILはCLR上でJITされ、x86コードに落とされる。 このときにJIT最適化がかかれば、拡張倍精度（を保ったまま）のコードになるし、 最適化がかからず毎回メモリに書き戻していれば、倍精度のコードになる。 .NETにおける同一バイナリってのは、 302 名前：同一MSILという意味であって、同一x86機械語という意味ではない。 だから、確かに同一バイナリを掴んでいたが、起動方法によって結果が異なっていたんだよ。 (VSがデバッグ用に意図的にそういう仕様にしていただけ。俺はそれを知らなかった） 君のコードについては、 いちいちfprintする場合はdouble(倍精度）が毎回必要になるから、 コンパイラはその部分での拡張倍精度でのループを断念し、 結果的にそのループが倍精度で回っているだけのこと。 もし仮にCが拡張倍精度型 ExDouble を持っていたとして、printfもそれに対応していれば、 もしかするとその毎回printfするコードでも拡張倍精度で回っていたかもしれん。 勿論手動でそういうコードのすることも可能だ。 そこはコンパイラがどう判断したのかでしかなく、あまり詰めても意味がない。 [] [ここ壊れてます] 303 名前：デフォルトの名無しさん mailto:sage [2018/09/17(月) 10:54:19.52 ID:+dwRu2dr.net] × fprint ○ printf まあ、分かると思いますが 304 名前：デフォルトの名無しさん [2018/09/17(月) 10:57:31.99 ID:dj7qSZnZ.net] で、最適化されてるかされてないかすら いまのいままで気付くことすらできない そして気付く方法すらわからなかったわけか うあわ 頭わるう。。。 ｷﾐ、ﾌﾟﾛｸﾞﾗﾑくむの向いてないわ 305 名前：デフォルトの名無しさん [2018/09/17(月) 10:58:47.63 ID:dj7qSZnZ.net] 問題の切り分けができない人は ﾌﾟﾛｸﾞﾗﾑはくめない ｺﾚは定説だからな 306 名前：デフォルトの名無しさん mailto:sage [2018/09/17(月) 11:09:05.96 ID:ivGPGa/P.net] >>287 なるほど面白いね。 レスが膨大過ぎて最初の方しか読んでなかったけど、ネイティブコードの話と思い込んでたら .NETの話だったのかw 307 名前：デフォルトの名無しさん mailto:sage [2018/09/17(月) 12:34:14.16 ID:yaPtorLJ.net] 個人的には、C++やx87 FPUは割と知識があったけど、.NET関連は余り追いかけてなかったので、気づくのが遅れた。 managedコード、unmanagedコードについて、今回初めて調べてみたくらいだし。 308 名前：デフォルトの名無しさん mailto:sage [2018/09/17(月) 12:53:34.77 ID:yaPtorLJ.net] >>293 「ID:dj7qSZnZ」の人は理解していないね。 彼は人の事馬鹿にしてるけど・・・。 309 名前：デフォルトの名無しさん mailto:sage [2018/09/17(月) 13:39:23.91 ID:F2vzl5VC.net] 最初に指摘されたことだろうに。 アセンブラレベルで精度や効率に介入したきゃ.netなんて使うな、なんて分かりきったこと。 310 名前：デフォルトの名無しさん mailto:sage [2018/09/17(月) 14:58:46.52 ID:zCVYDMXL.net] エディタ使ってるとたまに Intellisense機能が効かないときがあるんだが あれなんなの？ 中間ファイルとか消せば直るの？ 311 名前：デフォルトの名無しさん mailto:sage [2018/09/17(月) 17:03:46.83 ID:+dwRu2dr.net] >>301 はい。よく壊れます。 > [C++] There is an issue with the .ncb file > Close the solution. > Delete the . ncb file. > Reopen the solution. > Reopening the solution creates a new . ncb file. > https://msdn.microsoft.com/ja-jp/library/ks1ka3t6(v=vs.110) 312 名前：デフォルトの名無しさん mailto:sage [2018/09/17(月) 17:17:56.66 ID:zCVYDMXL.net] >>302 なるほど .ncbファイルね 今度消してやってみるわ 313 名前：デフォルトの名無しさん mailto:sage [2018/09/17(月) 18:30:58.83 ID:+dwRu2dr.net] さて俺の本番コード、以下のようだ。 疑問は解消した。協力してくれた皆様ありがとう。 ◎：拡張倍精度、○：倍精度、として、(ソースは>>192参照) ・Releaseビルドをコマンドプロンプトから起動→◎積和、◎平方根 ・Debugビルドをコマンドプロンプトから起動→◎積和、○平方根 ・IDEから起動→○積和、○平方根 これで3種類出来上がってた。 (なお、>>167内バイナリをアタッチした際の「AまたはC」は、「AまたはB」の間違い) そしてIDE上で『[モジュールの読み込み中に JIT 最適化を抑制する]』を変更すると、 確かにRelease/Debugの2種類に絞れる。 Debugだからといって、全く最適化がかからないわけでもないようだ。 （1行内なら最適化がかかる？） 参考に、Releaseビルドの該当部分の逆アセンブルは以下。 積和が拡張倍精度で行われ、そのまま fsqrt で平方根が取られる。 （関数ごとインライン化されているのでアドレスが中途半端だが） double retval = calc_norm_and_regulate(count, vec, false); 0000003e fldz 00000040 xor edx,edx 00000042 test esi,esi 00000044 jle 00000056 00000046 lea eax,[esp+28h] 0000004a fld qword ptr [eax+edx*8] 0000004d fmul st(0),st 0000004f faddp st(1),st 00000051 inc edx 00000052 cmp edx,esi 00000054 jl 00000046 00000056 fsqrt 00000058 fstp qword ptr [esp+10h] 314 名前：デフォルトの名無しさん [2018/09/19(水) 17:24:44.24 ID:GzaYiinu.net] VC++2008だけど、突然一部のファイルだけブレークポイントが入らなくなった・・・ .ncb消したり、そのプロジェクトだけリビルドしたけど直らない 全リビルドすればおｋ？ 315 名前：デフォルトの名無しさん mailto:sage [2018/09/19(水) 18:59:08.26 ID:KTUM/Tgt.net] >>305 根本的な解決策とは違うが 「ブレークポイントは現在の設定ではヒットしません。ソースコードが元のバージョンと異なります。」 なら [オプション] の[デバッグ]から[元のバージョンと完全に一致するソース ファイルを必要とする] をオフで 316 名前：デフォルトの名無しさん mailto:sage [2018/09/20(木) 12:44:39.12 ID:rLmEYUhk.net] >>306 ありがとうございます。 結局、全リビルドでも直らなかったので、.sln と .vcproj 以外全部消して Windowsも再起動して完全にまっさらな状態にして全ビルドしたら直りましたが、 そのうちまたブレークポイントが入らなくなったので、>>306の方法で 回避しました。 317 名前：デフォルトの名無しさん [2019/06/19(水) 14:12:16.93 ID:Yoy0IPRe.net] お、おう 318 名前：デフォルトの名無しさん mailto:sage [2019/10/04(金) 14:14:25.88 ID:gozwERy2.net] C++　CLRのWindowsフォームアプリ作ってると、 イベントとか記入するときインデントがおかしい（スペースが1つ付く）んですが、 これを直す方法ってありますか？ 319 名前：デフォルトの名無しさん mailto:sage [2019/10/06(日) 21:09:30.88 ID:xhkeezXX.net] テンプレート特殊化ってVisualC++2008 SP1 ではできないですか？ 320 名前：デフォルトの名無しさん mailto:sage [2019/10/06(日) 21:48:55.66 ID:xhkeezXX.net] >310 https://programming-place.net/ppp/contents/cpp/language/023.html こちらでできました、スレ汚しすみません 321 名前：デフォルトの名無しさん [2021/04/20(火) 20:00:45.05 ID:rVJPu0WM.net] わからないので教えてください。 シリアル通信をするプログラムを作りたく、 VS2008にて、「SerialPort」を使え、というのが サイトにあったのですが、 ツールボックスにSerialPortがありません。 追加する方法を教えてください。 使用環境：VS2008 C++ MFCアプリケーション 322 名前：デフォルトの名無しさん mailto:sage [2021/04/20(火) 20:36:15.98 ID:mPkgRyxP.net] ねぇねぇ今日は何処まで逝こうかな〜 COMはネイティブなのかね 323 名前：デフォルトの名無しさん mailto:sage [2021/04/23(金) 21:39:54.21 ID:OsunfFOw.net] >>312 SerialPortは使ったこと無いから知らんが、 その程度のことを自力で解決出来ない奴が今更VS2008でしかもMFCとか無理だ。 サイトを参考にするのもいいが、日付は必ず見るようにしろ。 .NET serial port で検索すると以下が当たるし、今なら普通にこれだと思うが。 https://docs.microsoft.com/ja-jp/dotnet/api/system.io.ports.serialport?view=dotnet-plat-ext-5.0 どうせ見えた展開だし、面倒だから先に言っておくが、 お前のやってることは初心者のあるある ・そもそもやり方を間違っているのに、それ以外の解決方法を認めない ・正しいやり方を教えても、それは聞いてないとして受け付けない だ。だからどうせこれから文句も言うのだろうが、重ねて言うが、 今更その環境(VS2008+MFC)はあり得ないくらいの異常さだ。 お前が初心者で何も知らないだけなら、常に最新の環境を使うように心がけろ。それだけで無駄なことに嵌りにくくなる。 その環境に拘る何らかの理由があり、それを強いてきた上司等が居るのなら、まずそいつに聞け。 聞く人が居ないのなら、今のお前にその環境でやりきる能力はないから諦めろ。 324 名前：デフォルトの名無しさん mailto:sage [2022/02/16(水) 02:27:10.00 ID:ky7QDxn6.net] 記念カキコ 325 名前：過去ログ ★ [[過去ログ]] ■ このスレッドは過去ログ倉庫に格納されています

---

---

[ 新着レスの取得/表示 (agate) ] / [ 携帯版 ]

read.cgi ver5.27 [feat.BBS2 +1.6] / e.0.2 (02/09/03) / eucaly.net products.
担当:undef