CでRの拡張したら速すぎて(40〜50倍)吹いたwww - yasuhisa's blog

昨日Gibbs Sampler Algorithmをやってみたわけだが、Rの中でfor文を書いていて必要となるサンプル数が多くなると非常につらくなってくることは目に見えている。しかも、MCMCでは初期値依存となる期間のサンプルを捨てないといけない。そういうわけでじゃんじゃんサンプルを作っても大丈夫なような速度が必要。 Rで速度を上げようと思ったらapplyファミリーを使うとかベクトル単位での処理をするetcが常套手段*1。が、今回は本質的にfor文が必要なケースである。 で、困るわけだがRにはC、C++、fortranを使って拡張する機能がある。詳しくはこの辺に載っている。そういうわけでCのポインタもアドレスも理解していないid:syou6162がRが好きすぎたためにCを書いてみたという感じの内容。 #include <R.h> #include <Rinternals.h> SEXP r

[syou6162]

これから拡張書きまくってRの中をもうちょい知ってみる予定

[MLAB]

Rの拡張

[ichan]

「はじめてのC♪　RとC♪ I will give you all my love 速度が出ちゃう Rのくせに♪ be in love with you」

[soulhack]

R.hとかあるのか，知らなかった～。

[kkobayashi]

Rの拡張。Rレベルでのチューニングの余地もあるかもしれない。cf. http://www.okada.jp.org/RWiki/?R%A5%B3%A1%BC%A5%C9%BA%C7%C5%AC%B2%BD%A4%CE%A5%B3%A5%C4%A4%C8%BC%C2%CE%E3%BD%B8

[Itisango]

R便利そう