「.NETのfloat.Epsilon定数は計算機イプシロンではない」の版間の差分
imported>Administrator |
imported>Administrator 編集の要約なし |
||
| 13行目: | 13行目: | ||
== 解決 == | == 解決 == | ||
「奥村晴彦『C言語による最新アルゴリズム事典』技術評論社，1991年，ISBN4-87408-414-1，2400円」という書籍に[[計算機イプシロン]]を動的に求める[[アルゴリズム]]が詳細に解説されており、また[[ソースコード]]は公式サイトにおいても配布されている。この本は非常に面白いので迷わず買え。個人的にはこの手の技術書を買う際に目次と技術評論社という名前だけで信頼して買うようになった一冊である。 | 「奥村晴彦『C言語による最新アルゴリズム事典』技術評論社，1991年，ISBN4-87408-414-1，2400円」という書籍に[[計算機イプシロン]]を動的に求める[[アルゴリズム]]が詳細に解説されており、また[[ソースコード]]は公式サイトにおいても配布されている。この本は非常に面白いので迷わず買え。個人的にはこの手の技術書を買う際に目次と技術評論社という名前だけで信頼して買うようになった一冊である。 | ||
https://oku.edu.mie-u.ac.jp/~okumura/algo/ | https://oku.edu.mie-u.ac.jp/~okumura/algo/ | ||
速攻で[[C言語]]からC#に[[移植]]した（ほぼ[[コピペ]]しただけ）。<source lang="csharp"> | |||
<source lang="csharp"> | |||
using System; | using System; | ||
| 48行目: | 48行目: | ||
</source> | </source> | ||
計算機イプシロンの定番の[[アルゴリズム]] | 計算機イプシロンの定番の[[アルゴリズム]]といえば、変数をひたすら2で割り続け、ゼロになった瞬間のひとつ前を捉えるというもの。「1以上」という条件はコンピューターの指数と仮数の組み合わせで表現される浮動小数点数では、1未満だと指数部だけを使ってどこまでも小さくなってしまうためであり、なんとしても仮数部を使おうということらしい。<source> | ||
<source> | |||
static double MachineEpsilon() | static double MachineEpsilon() | ||
{ | { | ||