float vs. double
å”çªã§ã™ãŒ, 下記ã®ã‚³ãƒ¼ãƒ‰ã‚’実行ã—ãŸã¨ã速ã„ã®ã¯ã©ã¡ã‚‰ã§ã—ょã†. ã¾ãŸæ¼”ç®—çµæžœã¯ã©ã†ãªã‚‹ã¨äºˆæƒ³ã•ã‚Œã¾ã™ã‹. double calcd(int n) { double ret = 0; for (int i = 1; i < n; i++) { ret += 1 / double(i); } return ret; } float calcf(int n) { float ret = 0; for (int i = 1; i < n; i++) { ret += 1 / float(i); } return ret; } int main() { printf("calcf=%.15f\n", calcf(100000000)); printf("calcd=%.15f\n", calcd(100000000)); } floatåž‹ã¯é…ã„? ã„ãˆ
後藤弘茂ã®Weekly海外ニュース - GeForce GTX 280ã®å€ç²¾åº¦æµ®å‹•å°æ•°ç‚¹æ¼”ç®—
â—å€ç²¾åº¦æµ®å‹•å°æ•°ç‚¹æ¼”ç®—ã®4ã¤ã®ãƒã‚¤ãƒ³ãƒˆ NVIDIAã¯ã€CUDAã¨G80アーã‚テクãƒãƒ£ã«ã‚ˆã£ã¦ã€ãƒã‚¤ãƒ‘フォーマンスコンピューティング(HPC)ã§ã‚る程度ã®æˆåŠŸã®è¶³ãŒã‹ã‚Šã‚’掴んã 。ãã†ã—ãŸNVIDIAã«ã¨ã£ã¦ã€å€ç²¾åº¦æµ®å‹•å°æ•°ç‚¹æ¼”ç®—ã®ã‚µãƒãƒ¼ãƒˆã¯ã€æ¬ ã‹ã›ãªã„è¦ç´ ã 。HPCã®ã‚¢ãƒ—リケーションã§ã¯ã€å€ç²¾åº¦ãŒå¿…è¦ã¨ãªã‚‹å±€é¢ãŒã‚ã‚‹ã‹ã‚‰ã 。GPUコンピューティングå‘ã‘ã®Tesla製å“ã§ã¯ã€å€ç²¾åº¦æ¼”ç®—ã¯ã‚«ã‚®ã¨ãªã‚‹ã¨è¨€ã£ã¦ã‚‚ã„ã„。「ã“ã‚Œã¾ã§ã‚‚ã€å€ç²¾åº¦ã®å£ã«ã‚ˆã£ã¦ç§»æ¤ã§ããªã‹ã£ãŸã‚¢ãƒ—リケーションãŒã‹ãªã‚Šã‚ã£ãŸã€ã¨NVIDIAã¯å¿…è¦æ€§ã‚’強調ã™ã‚‹ã€‚ å•é¡Œã¯ã€ç¾çŠ¶ã®ãƒªã‚¢ãƒ«ã‚¿ã‚¤ãƒ グラフィックスã§ã¯ã€å˜ç²¾åº¦(32-bit)ã¾ã§ã®æµ®å‹•å°æ•°ç‚¹æ¼”ç®—ã—ã‹å¿…è¦ã¨ã•ã‚Œãªã„ã“ã¨ã€‚ãã®ãŸã‚ã€ã“ã‚Œã¾ã§ã®GPUã¯å˜ç²¾åº¦æ¼”算ユニットã—ã‹å®Ÿè£…ã—ã¦ãŠã‚‰ãšã€GPUã§ã®å€ç²¾åº¦æ¼”ç®—ã®ã‚µãƒãƒ¼ãƒˆã«ã¯ã€ã„ãã¤ã‹ã®ãƒã‚¤ãƒ³ãƒˆã¨ç–‘å•ç‚¹ãŒã‚る。 (1)IEEE 75
Gauche:拡張浮動å°æ•°ç‚¹æ¼”ç®—ã®è¬Ž
SCMã®è¬Ž SCMã®æµ®å‹•å°æ•°ç‚¹æ›¸ã出ã—アルゴリズムã¯ã€Steele&White ("How to print floating-point numbers accurately, SIGPLAN '90, pp.112-123) ã«ç¤ºã•ã‚Œã¦ã„ã‚‹FP^3アルゴリズムã¨åŸºæœ¬çš„ã«åŒã˜ã§ã‚る。 ã™ãªã‚ã¡ã€æ›¸ã出ã—ãŸã„数値v (v ＞ 0) ã«å¯¾ã—ã¦ã€ ( m - ε) × 10^e ≦ v ≦ ( m + ε) × 10^e 1 ≦ m ＜ 10 ε ã¯ä¸Šã®å¼ã®æˆç«‹ã•ã›ã‚‹æœ€ã‚‚å°ã•ã„æ£ã®æµ®å‹•å°æ•°ç‚¹æ•° ã¨ãªã‚‹m, ε, e を求ã‚〠1. m ã®æ•´æ•°éƒ¨ã‚’書ã出ã™ã€‚ 2. ε ≦ mã®å°æ•°éƒ¨ ≦ 1.0 - ε ãªã‚‰ m ↠(mã®å°æ•°éƒ¨)×10.0 ε ↠ε × 10.0 ã¨ã—ã¦ã€1.ã‹ã‚‰ç¹°ã‚Šè¿”ã™ 3. ãã†ã§ãªã‘ã‚Œã°ã€ mã®å°æ•°éƒ¨ ＜ 0.5 ãªã‚‰ ãã®ã¾ã¾çµ‚了 mã®å°æ•°éƒ¨ ＞ 0.5 ãªã‚‰æœ€å¾Œã®
浮動å°æ•°ç‚¹æ¼”ç®—ã§ã¯ã¾ã£ãŸè©± - bkブãƒã‚°
浮動å°æ•°ç‚¹æ¼”ç®—ã§ã¯ã¾ã£ãŸè©± 浮動å°æ•°ç‚¹æ¼”ç®—ã®ã‚ã‚ŠãŒã¡ãªå•é¡Œã§ã¯ã¾ã‚Šã¾ã—ãŸã€‚ ã„ã‚ã„ã‚調ã¹ã¦ã„ã‚‹ã†ã¡ã« x86 特有ã®ã¡ã‚‡ã£ã¨ãŠã‚‚ã—ã‚ã„ ç¾è±¡ã«éé‡ã—ãŸã®ã§ç´¹ä»‹ã—ãŸã„ã¨æ€ã„ã¾ã™ã€‚ パーセンテージã®è¨ˆç®— ç°¡å˜ãª C ã®ãƒ—ãƒã‚°ãƒ©ãƒ ã§ãƒ‘ーセンテージを計算ã—よã†ã¨æ€ã„〠次ã®ã‚ˆã†ãªã‚³ãƒ¼ãƒ‰ã‚’書ãã¾ã—ãŸã€‚ int x, y; ... int a = (double)x / y * 100; int a = x * 100 / y ã¨ã—ãªã„ã®ã¯ã€ x ãŒå¤§ãã„å ´åˆã« x * 100 ㌠オーãƒãƒ¼ãƒ•ãƒãƒ¼ã‚’èµ·ã“ã™ (INT_MAX を越ãˆã‚‹) ãŸã‚ã§ã™ã€‚ ã“ã®ã‚³ãƒ¼ãƒ‰ã¯ä¸€è¦‹ã€æœŸå¾…通りã«å‹•ã„ã¦ã„ãŸã®ã§ã™ãŒã€ ã—ã°ã‚‰ã使ã£ã¦ã„ã‚‹ã¨ã€æ‰‹å…ƒã®ç’°å¢ƒã§ã¯ x = 53, y = 100 ã®ã¨ãã« a 㯠53 ã§ã¯ãªã 52 ã«ãªã‚‹ã“ã¨ã«æ°—ã¥ãã¾ã—ãŸã€‚ ã“ã‚Œã¯æ¬¡ã®ç†ç”±ã«ã‚ˆã‚Šã¾ã™ã€‚ å¼ã®æœ€åˆã® (double)53 / 10
1
リリースã€éšœå®³æƒ…å ±ãªã©ã®ã‚µãƒ¼ãƒ“スã®ãŠçŸ¥ã‚‰ã›
最新ã®äººæ°—エントリーã®é…ä¿¡
j次ã®ãƒ–ックマーク
kå‰ã®ãƒ–ックマーク
lã‚ã¨ã§èªã‚€
eコメント一覧を開ã
oページを開ã
Manpage of FENV
インテル(R) アーã‚テクãƒãƒ£ (IA) 浮動å°æ•°ç‚¹ãƒ¦ãƒ‹ãƒƒãƒˆ (FPU)ã€ã‚¹ãƒˆãƒªãƒ¼ãƒŸãƒ³ã‚° SIMD 拡張命令 (SSE)ã€ã‚¹ãƒˆãƒªãƒ¼ãƒŸãƒ³ã‚° SIMD 拡張命令2 (SSE2) を使用ã—ãŸæµ®å‹•å°æ•°ç‚¹ç®—術演算
{{#tags}}- {{label}}
{{/tags}}