DONLP2 の下調べ

今日は、ほとんどの時間を DONLP2 の下調べに利用。ユーザマニュアルがあるにはあるけど、これが読みにくく苦戦。まだ、ソースの中身はチェックできていない。

あと、他にどんなソフトがあるのかも下調べ。有力そうなものは Ipopt というもので、COIN-OR の一部にもなっている。(Ipopt へのリンク)

ということで、明日は Ipopt をざっとみるのと、DONLP2 の詳細チェックになりそうだ。

今日の作業内容: DONLP2 チェック 5h
今日のBGM: FF12 OST [1,2,4]
今日のランチ: らく うな重
明日の予測作業時間: 3h

計算サーバチェック

今日は、計算サーバとしてどんなスペックがあったらいいか、を検討してみた。
すでに、大きなＰＣクラスタなどは他のところにもあって、似たようなのを目指すのは効率がよくないので、ここは思い切って、いろんな実験ができるようなものを探してみることにする。

たとえば、グラフィックスボードで BLAS の計算ができる、という CUDA は、名前は聞いたことのあるものの、実際に使ったことがないので、CUDA の BLAS を SDPA に組み込めたらどうなるか、などは確認してみたいところである。

今きがついたが、今日はうっかり DONLP2 を調べる時間を取ってなかった。これは、明日に回そう。


今日の作業内容: 計算サーバチェック 3h + 書類処理 2h
今日のBGM: Madlax OST [1-2]
今日のランチ: 味庵 かにとレタスのチャーハン
明日の予測作業時間: 3h

Python と convex な日

今日は、まず Mangasarian の本に目を通してみた。ひとつひとつの証明がなかなかに面白いのだが、全部読んでいるときりがないので、定理で言っていることだけを確認。それにあわせて、そもそも何でこの本を読んでいたかを確認。

もともとの動機としては、ある集合 F \in \Real^n に対して
$ W \in conv\{vv^T : v \in F\}$
が簡単に判断できるか？ということであった（のを思い出した）。もちろん、v \in F を判断するのと同じ程度か n 倍ぐらいの計算時間で判断できるのか？ということ。

これが案外難しいらしく、簡単には見つからない。Rockefeler の convex analysis, Hiriart-Urruty & Lemarchal の Fundamentals of Convex Analysis もチェックしてみたが、あまり重要なことは書いていない。ただ、Caratheodory の定理は、使い道があるかもしれない。

あと、今日は Python について、少し調べてみた。
CVXOPT と CHOMPACK というソフトがあるようで、特に CHOMPACK は sparse Cholesky factorization ができるということが書いてあるので、要チェック。ただし、両方とももちろん Python な感じである。そうそう、CVXOPT は Ubuntu では標準パッケージがあった。SDPA にも Ubuntu のパッケージがあったらいいかもしれない。

そういえば、今日は DONLP2 が進まなかったので、明日もうちょっと確認してみようと思う。

今日の作業時間: Mangasarian チェック 3h + Python チェック 1h + その他 1h
今日のBGM: Macross F OST [1-2], Chrono Trigger OST [1-3]
今日のランチ: いろは 豚アスパラチーズ焼き
明日の作業時間予測: 6h

TSUBAME見学と DONLP2

今日は、午前中に TSUBAME の見学をさせてもらうことができた。ラッキーである。

普段はネットワーク越しにログインしているだけなので、実際にどんな風になっているのかはコマンドラインからだけではよく解からなかった。見学をしてみて、各ノードに貼ってある各ノードの番号を示す小さなシールを見ると、自分がログインしているノードに似たような番号が貼ってあって、「おぉ、やっぱり本物だ」という感慨があった。
ファンが回っているので、騒音が大きく長時間滞在できるわけではないが、こうやって実際のノードを見てみるというのも大事なことだと改めて思う。

そのあと、いろいろと話を聞いたけど、
DONLP2というソフトは勉強になるようなので、明日あたりもう少し詳しく調べてみたい。
DONLP2へのリンク
ただ、DONLP2 は AMPL をインターフェースの一種として使っているので、AMPLについても勉強が必要だ。
AMPLへのリンク

ほかに今日得た知識では、

1. QAP には、positive definite と negative definite に分けるという今までと違った定式化があって、これはひょっとしたら調べる価値があるかもしれない
2. kissing number は SDP 緩和で近似解が求まる（ただ、詳しい定式化は確認しなおしたほうがいい）

などなど、今日はとても勉強になった。

今日の作業時間: DONLP2 チェック 2h
今日のBGM: なし
今日のランチ: シッダルータ シーフードカレー
明日の予測作業時間: 6h

PDSDPのコンパイル

今日は、PDSDP のコンパイルをしてみた。久しぶりのコンパイルだったので、なかなかてこずったけど、ソースを微修正してコンパイルを成功した。やはり、メンテナンスされていないので、細かい綻びがある。で、ちょっと解いてみると、SDPARA, PCSDP に比較するとかなり遅い。dual interior-point method だと primal-dual interior-point method のようにスピードが出にくいようだ。

あと、連休の間に他の研究内容に関係ありそうな論文をチェックした。
On the Minimum Volume Covering Ellipsoid of Ellipsoids, SIAM Journal on Optimization, 17(3) , 621-641, 2006

やっている内容は面白いんだけど、証明などに一部きちんとできていないところがある。
とりあえず、明日はMangasarian の本を読むのを進めたほうがよさそうだ。

SDPA の論文は、一通り校正が終わったので、これで他の人の校正をお願いしたところ。ついに終わりが見えてきた（という気が、個人的にはする）。

今日の作業内容: PDSDP コンパイル 2h + SDPA 論文校正 2h + 計算機不具合解消 1h
今日のBGM: FF7 Reunion, ARIA the Animation, ARIA the Natural, ARIA the Origination
今日のランチ: さか本 ねぎとろとろろ丼
明日の作業時間予測: 3h

ディスカッション

今日は、英語でのディスカッションをしてみた。研究内容になると英単語のボキャブラリが狭くなるので、かなり楽にはなるけど、英語についてはまだまだ改善の余地が大きい。

教えてもらったことの要点をまとめておくと、

1. CPU の cache size は性能に大きく影響するので、計算時間予測に使ったほうがいい。(SDPT3はひょっとしたら cache size を考慮しているのかも）
2. 似たような研究をしている人に、数値結果を送れば興味を持ってくれるので、連絡してみるのがいい。
3. どういった問題が並列計算に向いているか気がついたら、改めて教えてくれる、とのこと。

英語は大変だったけど、とってもためになった。

今日の作業内容: ディスカッション準備 2h + ディスカッション 1.5h + その他雑多 1h
今日のBGM: Madlax OST [1]
今日のランチ: たちばな シマアジ丼
明日の予測作業時間: 3h

お迎え

今日は、朝一に品川までお迎えに行ってきた。そのあと、一緒に職場まで移動して勤務開始。
ただ、そのおかげで雑多なことが思っていた以上に発生したので、今日は論文探索しか進んでいない。

今回読んでいた
Support Vector Data Description, Machine Learning, 54, 45-66, 2004
に元になっている論文があるというので、この作者の論文をチェックする。
Support Vector Data Description, Proceedings of the Fifth Annual Confenrence of the ASCI, 398-405, 1999
これにも元になっている論文があるというので、さらに孫引。
Data domain description using support vectors, the special issue of Pattern Recognition Letters, proceedings of Pattern Recognition in Practice VI, 1999

ただ、これらの内容はかなり似ていて、どうやら下２本の会議の内容を論文として投稿したのが1本目のようだ。
ということで、1本めの論文で、もう一個重要だと書かれていた論文を明日チェックすることにする。

ちなみに、作業として残っているのは、
1. 論文探索 2. 論文校正 3. Mangasarian の本をチェック
といったところ。

今日の作業内容: SVDD チェック 2h + 雑多な仕事 2h
今日のBGM: 電脳コイル OST [1-2]
今日のランチ: らく うな重
明日の作業時間予測: 4h