茨城の山奥に行ったりとか，書類書いてたりとかして1週間も経ってしまいましたが書きます．
"CilkPlus, TBB, OpenMP 〜 並列プログラミングを支える技術" という大層なタイトルで発表してましたが，つまるところこの辺のライブラリとかコンパイラ拡張の紹介です．
発表スライド: Boost study.tokyo11
CilkPlus, TBB については以下のようなページがありますので，こちらを見ると良いと思います．
後は CG 法自体についての説明ですが... まあ興味ある方は調べてみてください．簡単に言うと線形代数で連立1次方程式を解くための解法です．
このスライドでは，ベクトル演算中心で並列化が容易だったので，それぞれのライブラリやコンパイラ拡張だけで差が出たら良いなあと思っていただけですので．
今回解いた問題で用いたのは Givens 行列という行列なんだそうです．
未知ベクトルを全て 1 として既知ベクトルを求めた後，未知ベクトルを全て 1 以外で初期化し，CG 法で解けば，未知ベクトルは 1 に近似された値が入る，となります．
ああ，後，行列ベクトル積のそれぞれのサンプルコードを載せて「簡単に書けますね」しか言わなかったのですが，もちろんご指摘がありましたように「簡単に書ける = パフォーマンスが良い」という風にはもちろん思ってないです．
実際にはパフォーマンスを調べてみないといけない訳ですが，一般的なアプリケーションにおいては簡単に書けて，かつスケーリングに強いというの方が良いと思っています．(もちろん，パフォーマンスが逐次計算時よりも良くなっているのが前提)
コンテナのパフォーマンスに関する発表では，会場に居ない人の代わりというわけでもなく色々指摘しましたけど，ちゃんと思ったことは言った方が良いですよ．
勉強会で発表してる人って，「僕何か間違ってないかな...」とか思いながら発表してる部分もきっとあると思います．(間違いに気づいてない人もいますけど) そういうのは直接言ってあげないと気がつかないと思います．僕もそんな感じなので．
懇親会も楽しかったです．ラズベリーパイとバイクと某氏の妹の話しか覚えてないですが．
#12 大阪は行きたいのですが，まあ色々予定が辛いので難しいですね...
闇鍋は夏が終わるまで無理です．暑いし．
]]>コード書くの疲れたので休憩も兼ねて書きます．
知っている人も多いと思いますが，CilkPlus は並列処理を容易にするため，Intel が自社のコンパイラに用意している C/C++ 用の言語拡張です．
元々は MIT で作られていた C 言語用の拡張だったとかなんだとか．
まあそれはともかく，CilkPlus はキーワードが3つしかないとてもシンプルな設計です．
の3つだけです．
後は reducer などが用意されています．
他には，ベクトル演算を簡略化する Array Notation，関数をベクタライズする Elemental Functions，ループのベクタライズを明示的に許可する #pragma simd などがあります．(これらも CilkPlus の一部だったのか...)
これらの定義は，全て <cilk/cilk.h> にあります *1．
CilkPlus のことについては以下のサイトが良いと思います．
現在のところ，Intel コンパイラなら使用可能ですが，それ以外のコンパイラでは GCC は 4.8 で CilkPlus が使えるブランチがあり，4.9 では development branch があります．また LLVM/Clang でも実装されており，github に置かれています．詳しい話は以下．
cilk_for は並列 for 文です．OpenMP だと pragma omp parallel for に該当します．
#include <cilk/cilk.h> template<class T> void add(std::vector<T> & out, std::vector<T> const& x, std::vector<T> const& y) { const std::size_t n = x.size(); cilk_for(std::size_t i = 0 ; i < n ; ++i) { out[i] = x[i] + y[i]; } }
cilk_spawn は関数を非同期実行させる命令で，cilk_sync は cilk_spawn で実行した全ての処理が完了するのを待機します．
#include <cilk/cilk.h> int fib(int n) { if (n < 2) { return n; } int x = cilk_spawn fib(n-1); int y = fib(n-2); cilk_sync; return x + y; }
reducer は並行プログラミングでは非常に困る，複数のスレッド間で同期しなければならない計算を助けるための機能です．
例えば総和計算の関数は，全スレッドが参照する総和結果を返す変数をロックしなければいけませんが，reducer はそのようなロックをうまく隠してくれます．
#include <cilk/cilk.h> template<class T> T sum(std::vector<T> const& x) { cilk::reducer_opadd<T> result; // 自動でゼロ初期化されているらしい． cilk_for(std::size_t i = 0 ; i < x.size() ; ++i) { result += x[i]; } return result.get_value(); }
OpenMP と比べると，#ifdef とかない分書きやすいと思います．
CilkPlus はオプションの指定が必要ないため，使える環境ならどこでも使えます．
後，OpenMP ではリダクションの指定が色々制限が強いのですが，CilkPlus なら reducer を自前で実装することで，ユーザ定義型でも使用できるようです.
Tutorial Cilk Plus Reducers | CilkPlus
そういえば OpenMP と CilkPlus の比較記事をあまり見たことがない．
どこかでやっているような気はするんですが...
*1：cilk_for, cilk_spawn, cilk_sync は実際にはマクロで，実際のキーワードは _Cilk_for, _Cilk_spawn, _Cilk_sync なのですが，cilk_for, cilk_spawn, cilk_sync を使うのが推奨されています
無視する．
タイトルが学籍番号と名前で，目的は何で，何があったのか不明，「出席管理サーバに登録してください．」だけで「よろしくお願いいたします」と言ったしめもなし．*1
おまけにメールフッターはついてない．
友達気分でメールし過ぎじゃないかと思いましたね．
僕はお前の友達でもなんでも無いんだが．
目上の人間にメール送るって意識が足りなさすぎですね．
というか失礼過ぎ．
今回はそもそも問い合わせ先が変わってるのでそっちに送ってねで対応しましたが，そうじゃなかったら完全に無視してましたね．
*1：本当は学生が自分で登録する
最近 TeX に載せるグラフを書くために Numbers とか使って EPS に変換したりとかしてたのですが, よくある2段組み用に使うと汚くてかなり見づらい.
今日とうとう諦めて PLplot を使って描画プログラムを書くことにしました.
PLplot は pgplot みたいなグラフ生成用のミドルウェアです.
バインディングは C, C++, Python, Fortran など.
EPS で出力すると結構奇麗な表示になって大変良いと思います.
ただ C++ の example が大変残念だったり, 凡例表示用関数が引数30超えてるとか*1まあ色々不満はあるのですが, ベタで書くよりは楽だろうと思います.
僕はいちいちデータファイルを変更するためにコンパイルし直したくないのでコマンドライン引数に渡したかったのですが, PLplot にオプション解析を乗っ取られているのでどうすれば良いのかと思ってたのですが, なんとなく調べても出てこなかったので一応メモ.
plparseopts 関数が引数を解析しますが, サンプルだと PL_PARSE_FULL だけが渡されているのが多いと思います.
ここに | で PL_PARSE_SKIP を指定すれば PLplot が解析できない引数をスキップできます. スキップしないとエラー出力で落ちます.
plparseopts 関数を通った後, argc, argv の内容は変更され, PLplot が認識した引数は全て削除されます. なので PLplot のオプションだけを渡せば, argc は 1 になります.
グラフに使用するデータファイルをコマンドライン引数に指定させようとすると, こういう書き方になります.
int main(int argc, const char** argv) { plparseopts(&argc, argv, PL_PARSE_FULL | PL_PARSE_SKIP); if(argc != 2) { std::fprintf(stderr, "[Usage] %s <PLplot options...> data-file\n", argv[0]); return -1; } const std::string dataFile = argv[1]; ... }
PLplot はちゃんと設定すれば Numbers で描画したグラフと遜色ないぐらい奇麗なグラフを EPS で出力できるので, TeX 使う人には便利だと思います.
Ruby のバインドがあれば良いなあとは思いました.
]]>これをやるきっかけは, 大学の1年生が主に使っているシステム (がマウントしているファイルサーバだったかな) が今年度でサポート完全停止でいつ壊れても修理できないそうで, 今学内全体で使うためのクラウドシステムの準備中らしく, その期間安全に使えるシステムがないとかなんとかで, 出席管理が欲しいとかそういう理由だったと思います.
FreeBSD は別にいいんですが, Moodle は説明します.
"Moodleはインターネット上で授業用のWebページを作るためのソフトです。 教育学でいう社会的構築主義の考え方に基づいて作られており，日々改良が行われています。"
とのことで, 簡単に言うと e ラーニングシステムみたいなものだと思います.
先生にその名前を聞くまで知らなかったのですが, セットアップして使えるようになるまで忘れてました. 僕も微積分の講義で使ってました. というかデザインを見て思い出したので, 「あああれは Moodle だったのか」という感じです.
Moodle の設定は上のページのように日本語で公式のドキュメントもあるので (若干ふるい気がする) そこまで難しくないと思います.
必要なのは SQL サーバの設定と, PHP の設定, Apache の設定なので.
FreeBSD はサーバをたてるかって時に「Ubuntu はサーバとして信用ならない」という先生が居たので FreeBSD にしました. あと Ports に Moodle も入ってたからです.
インストール自体は公式のドキュメント Moodleをインストールする - MoodleDocs を見れば大体できます.
問題は, Ports でインストールするとディレクトリが /usr/local/www/moodle, /usr/local/www/moodledata になると思います.
Apache のドキュメントルートは /usr/local/www/apache22/data 辺りになるので, 修正したいのですが, ドキュメントルート自体はそのままで, /moodle とウェブでアクセスするときは /usr/local/www/moodle を使用するようにしたい.
それ自体は簡単で Apache の Alias モジュールを使えば解決.
mod_alias - Apache HTTP Server
なのでかなり簡単です.
問題は, 僕にウェブの知識がほとんどないので詰んでる.
既に設定は全て終わってるので, 後は実際に使う先生方に Moodle の管理者権限を渡せば完了.
]]>