melancholic afternoon http://homepage1.nifty.com/herumi/diary/latest.html With Wnz ja Sun, 08 Nov 2009 22:24:43 GMT tekito ver. 2.0 11/8 float xに関するexp(x)の実... http://homepage1.nifty.com/herumi/diary/0911.html#8 Sun, 8 Nov 2009 14:59:59 GMT
  • e-x = 1 / exだからx >= 0とする.
  • ex = Σxi / i!を使うのだけど収束に時間がかかるのでx = n + d(n∈Z, d∈[0, 1))と分解しよう. ex = en ed. n～128程度でfloatの最大値になるからそちらはテーブル引きすればいい.
  • その場合11!=39916800～4.0e7だから1e-7の精度を出すにはΣはi = 11ぐらいまでとる必要がある. これは大きすぎる. d >= 0.5の場合はd = 0.5 + d1, d1∈[0, 0.5)と更に分解してみよう. ex = en e0.5 ed1
  • この場合0.58/8!～1e-6ぐらい. うーん, まだ大きい.d1 = 0.25 + d2, d2∈[0, 0.25)ともう一つ分けるか. それでも0.256/6!～3e-7. いまいち. 別の方法を考える.
  • d1, d2は要は2進分数なんだからそちらもテーブル引きすればよさそう.
  • x = n + (2進数8桁のd1) + (次の2進数8桁のd2) + ...として全部テーブル引きでいいか. なら級数展開は不要だなあ.
]]> 11/7 昔作ってたlogの計算を... http://homepage1.nifty.com/herumi/diary/0911.html#7 Sat, 7 Nov 2009 14:59:59 GMT logの計算をintrinsicでやってみた(log.cpp). とりあえずasmを焼き直しただけなので改良の余地は多々あるだろうがCore2Duo 1.8GHzでベンチ.
32bit VC9 32bit VC9 /fp:fast 64bit VC 32bit gcc 4.1.2 64bit gcc 4.1.2
log 0.687 0.375 0.188 0.48 0.330
my_log 0.156 0.125 0.130 0.110 0.100
計測がいいかげんなので有効数字2桁なさそうやけど32bitで概ね3～4倍. 64bitのVCのlogはがんばってるなあ. 昔はdivpsよりもrcpps+ニュートン補完が速かったんやけど今は違うのか?(要検証) まあ, 昔の最適化が現在でもそれなりに意味はあると分かったので次はexpを考えたい. ]]> 11/6 from_time_tの問題は1.41.0で... http://homepage1.nifty.com/herumi/diary/0911.html#6 Fri, 6 Nov 2009 14:59:59 GMT 対応するっぽい. と思ったらまだ駄目だった. ]]> 11/5 ・第二回カーネル/VM探... http://homepage1.nifty.com/herumi/diary/0911.html#5 Thu, 5 Nov 2009 14:59:59 GMT 第二回カーネル/VM探検隊
ほどほどに面白かった. ただ16時～20時半の長丁場は辛かった. 途中でおなか空きすぎて頭が半分回ってなかった. MINIX SMPのところで66h, 67hをnopに置き換えてる理由がわからないと言ってはったので32bitモードで16bitコード生成したいからとちゃうって適当に答えてしまったが, ソースのコメントに書いてあった.
・副都心線初めて乗った. 渋谷で上に出ようとして階段使ったら100段あった. なんて深いところにあるんや. ]]> 11/4 ・CANS2009で発表するF212... http://homepage1.nifty.com/herumi/diary/0911.html#4 Wed, 4 Nov 2009 14:59:59 GMT 21223, F3509上のTateペアリングのソースコードが公開された. 私は途中参加だったのでC++ならもちっと楽できるんだけどなあと思いつつ郷に従う. なにはともあれ, こういう最速狙うコードで公開されてるのは殆ど無いのでめでたい.
・下記のposix_time::from_time_tの件, 2007年に既に指摘されてるのに直ってないということは直す気ないんやろか. trackが見えんでよう分からん. こっちの人もfrom_time_tについて何か言ってるなあ. これも2年前だ. ]]> 11/2 boost::posix_time::from_time_t... http://homepage1.nifty.com/herumi/diary/0911.html#2 Mon, 2 Nov 2009 14:59:59 GMT #include <stdio.h> #include <time.h> #include <iostream> #include <boost/date_time/posix_time/posix_time.hpp> #include <boost/date_time/gregorian/gregorian.hpp> void Put(time_t t) { struct tm tm; char buf[64]; #ifdef _WIN32 _gmtime64_s(&tm, &t); asctime_s(buf, &tm); #else gmtime_r(&t, &tm); asctime_r(&tm, buf); #endif printf("time=%s", buf); } int main() { time_t t; boost::posix_time::ptime bt; time(&t); printf("sizeof=%d\n", sizeof(t)); Put(t); bt = boost::posix_time::from_time_t(t); std::cout << "boost:" << bt << std::endl; std::cout << "boost:" << (bt + boost::gregorian::years(29)) << std::endl; t += 84600 * 365 * 29; Put(t); bt = boost::posix_time::from_time_t(t); std::cout << "boost:" << bt << std::endl; } sizeof=8 time=Mon Nov 02 13:39:33 2009 boost:2009-Nov-02 13:39:33 boost:2038-Nov-02 13:39:33 time=Sat Mar 20 01:09:33 2038 boost:1902-Feb-11 18:41:17 ; おかしい sizeof(time_t) = 4な環境で駄目なのはしかたないけど, そうじゃない環境では対応しておいて欲しいなあ. 言ってみるか. 
  ptime from_time_t(std::time_t t)
  {
    ptime start(gregorian::date(1970,1,1));
    return start + seconds(static_cast<long>(t));
  }
longにキャストしてるからか? でも64-bit Linuxでも駄目なんだよなあ. ]]>