メガネの歌のうまい学生

昨日、またもやGotoBLASについてコメントをいただいたので、「メガネの歌のうまい学生」に連絡を取ってみようと思ったところ、予想外にもメールアドレスをまだ入手していなかったことに気づいた。
すでに教えてもらってあると勘違いしていた。うっかり。

ということで、他の人にメールアドレスを教えてもらうか、Blog 経由で連絡にするか、それとも mixi 経由にするか、であるが、どれが一番面白いだろうか、を現在検討中。

ところで、投稿用の文章は今日も１ページ程度進んだ。そろそろ数値実験結果が入ってこないと文章が組み立てづらくなってきたので、数値実験もあわせて進めていきたい。

今日の作業内容:投稿用準備 4h
今日のBGM: SKY OST [1-2], DewPrism OST [1-2]
今日のランチ: シッダルータ ほうれんそうとチキンのカレー
明日の予測作業時間: 5h

投稿用、まだまだ続く

今日はさらに１ページ追加して、図などを追加。あと、参考文献も追加したけど、こちらはフォーマットが違うので、明日調整する必要あり。

数値実験も準備段階は終わりつつあって、どれくらいのサイズの問題を解けばいいのか解かってきたので、それを使って実験をスタートしようと思う。

今日の作業内容: 論文作成 4h
今日のBGM: MADLAX OST [1-2], FF6 OST [1,2]
今日のランチ: らく 焼魚定食
明日の予測作業時間: 4h

投稿用の続き

今日は投稿用の文章の続き。１ページ分追加。ただ、今のペースで書いていくと、６ページ持たないかもしれない。自分のやったところをバランスよく増やして、6ページになるようにペース配分をしたい。明日は、できればもう１ページ追加しておきたい。

これと並行して数値実験。制御の問題で、特定のサイズで乱数で発生させている問題が pdINF でうまく解けていなかったが、初期点のパラメータを変更したところ、それなりに解けるようになった。どうやら中心パスから、かなり遠いところでスタートしていたようだ。

あと、昨日いただいたコメントに基づいて GotoBLAS2 の起動にかかる時間をチェック。Makefile.rule で、一度
NO_WARMUP = 1
を設定したが、これでは変わらず。make clean をしたかどうか覚えていないので、いっそのこと tar+gz を展開するところからやり直し。これで起動にかかる時間が正常になった。
ところで、このWARMUP って何だろう？今度時間があるときに何をしていて、それでどうして遅くなっていたのかを調べることにしよう。


今日の作業内容: 投稿用続き 4h + 論文対応 1h
今日のBGM: FF10 OST[1,2,4]
今日のランチ: しばた さけいくら丼
明日の予測作業時間: 4h

投稿用の文章の作成

今日は、１月〆切の投稿用の文章の作成。6ページ or 8ページのあまり長くない文章ではあるけど、正月をはさむことを考えるとあまり時間に余裕がないので、とりあえず書きはじめ。ただ、今日は英語の細かいところを気にせずにラフに書いて１ページ分しか進まなかった。思っていたよりも時間がかかりそうなので、これを早めに進めることにしよう。

これにあわせて行っている数値実験だけど、SDP の問題生成ルーチンが乱数で生成している関係で、変なサイズだと pdINF で終わってしまうことがある。まずは、別の乱数で生成しなおしている。

あと、別に投稿していた論文の refree からのコメントへの対応策を考え中。こっちは、調べれば調べるほど深みにはまりそうな展開で、さらっと流すほうがいいのかもしれない。

ところで、今日、このブログに初めてコメントがついた。どうやら、このブログにもコメントをつける機能があったらしい。初めて知った。GotoBLAS のコンパイルに関する内容で、かなり細かいところを教えてもらったので、明日以降文章を書くのと並行してコンパイルしてみようと思う。

今日の作業内容: 文章作成 3h + 論文対応 1h
今日のBGM: MADLAX OST [1-2]
今日のランチ: 信華園 鶏肉両面焼きそば
明日の予測作業時間: 4h

SDPA-Python Interface を実装

昨日からの続きの SDPA-Python Interface だが、今日のあまりある時間をフルに投入して、ファイル書き込みと sedumiwrap を実装した。これにより、例えば control7.dat-s などは、

import sdpap;
At,b,c,K,blockStruct = sdpap.fromsdpa('control7.dat-s');
x,y,info = sdpap.sedumiwrap(At,b,c,K);

のように sedumi っぽく書くことができる。もちろん、中身で動くのは SDPA のアルゴリズムそのもので、インターフェースだけが変わっていることになる。(SDPA の callable library にリンクしている。)
とりあえず自分のところではコンパイルできるようになったので、ほかのところ(Macなども)でもできるのか、確認してみたいところだ。

あと、SIAM Review の最新号が到着。 Lassere の論文が、この前投稿した論文に似たようなタイトルだったので、ちょっと焦ったが、内容はかなりちがかったので、ほっと一安心。

今日は、あとは SDPA7 の論文の変更すべき点の確認と、他の論文のフォーマットの確認。
明日は午前中に論文構成を練ることにして、午後に SDPA7 の論文をチェックしよう。

今日の作業内容: Python 実装 5h + 論文確認 1h + フォーマット確認 1h
今日のBGM: FF12 OST [1-4]
今日のランチ: いろは チキンステーキ
明日の予測作業時間: 6h

Python で SDP を標準形に

今朝はうっかり２度寝をしてしまい、１時間ぐらいのロス。

今日は Python の続きで、SDPA sparse format のファイル出力に先立って、標準形に変更する (free 変数をふたつの非負変数にする) ところを書いた。ある程度 Matlab に近い形で書けるので、かなり便利だ。

あと、Python で型を調べるときには、type などもあるが、どうやら isinstance を使うのがいいようだ。type だと自分で作った class の instance を、class との名前で判断できないときがある。

たとえば、
def myclass(object):
　　クラスの中身
であるときに、
>>> myclassI = myclass();
とすると、
type(myclassI) is myclass
は False になるが
isinstance(myclassI,myclass)
は True になる。

明日は、ファイル出力関数を C/API で書く。
また、SDPA の論文についても修正点を整理しないと。

今日の作業内容: Python 3h + その他 1h
今日の BGM: Chrono Trigger OST [1-3]
今日のランチ: シッダルータ　チキンとひよこ豆のカレー
明日の予測作業時間: 6h

Python で fromsdpa

今日は、SDPA フォーマットからファイルを読み込む fromsdpa() を作成。ざっと 400 行ぐらいのコードでできた。アルゴリズムは mex と同じで済むのが有り難い。

ただ、まだ Python になれていないので、変なバグに遭遇すると時間がかかる。今日も、help(sdpa) としたら、「配列の要素が越えてるぞ！」という意味不明のバグに遭遇し、30分消費。結局、pythonを起動しなおしたら、バグが取れた。いまひとつ謎。

あとは、SDPARA の数値実験開始。思っていたよりも数値精度がでないときがあるので、明日結果を見て判断しよう。制御の問題を生成するプログラムも作ったし、これの結果も確認していきたい。

SDPA の論文は、まだまだ英語がこなれていないので、英語を頑張ろうと思う。英語は、単語の微妙なニュアンスが難しいけど、そこが頑張りどころ。

今日の作業内容: python 5h + SDPARA 数値実験 1h + SDPA 論文 1h
今日のBGM: EVA OST [1-3]
今日のランチ: 味庵 回鍋肉
明日の予測作業時間: 3h