Open Diary - June 2008

6/30 (Mon)

[phys]ryuon
- 高速な implicit algorithm の実装：
  - 前回のテストの続き：
    - N=100 で NITSOL がこける問題、解決。
      - 問題は、またもや config の設定（else if すべき所が一ヶ所 if だけになっていた、というおバカなもの）。
    - 以上で no-HI 版のテストは終了とする。
  - HI 化：
    - ノートに戻って、定式化を確認。
      - ふむ、 right-hand side vector のみの変更で対応できる。
    - しかるべく、実装。
      - できた。 straight chain の F の決め打ちで、しかも non periodic だけ。
      - ついでなので corrector step の初期と update の統合もしておく。（詳しい理由は知らないけど、微妙に違うんだよな。この違いについては、後でちょっと調べておこう。）
      - これまでのチェックで通ることを確認ののち、古い（非統合）コードは破棄。
    - これで、一応 HI 化はできたはす。
      - 不十分な点がもう一つ、 connector から position を作るのに、 COM を手で与えている。
    - テスト：
      - 既存のテストに flag_noHI を追加。一通り走るようになった！！
      - と思っていたら、何かの加減で check_fastSI_solve_cubic() が文句を言うようになった。
        
        見ると、いろいろ下らない bug や、思い違いがあった。
        
        それらを修正後も、まだどこかで nan になっている模様……。
      - 7/3/2008: テストの続き。
    - 今月はここまで。
[monton] モンちゃんの散髪。
- をするために、バリカンを買って帰る。（あと、モンちゃんの暑さ対策として氷嚢を１つ。）
- いくつか店を回って、結局４０数ドルの一番安い light duty という奴を買った。他の店でみたものよりも、 quality はよさそうだったので。
- で、さっそく試してみた。が、これがなかなか思うように刈れない。本当は目の前、鼻の上の毛を（病院で刈られたように）切りたかったのだが（モンちゃんは、おめめが大きいので、ハサミだとちょっと怖い）、そんなことは不可能なほど、弱い。刈りやすそうな背中から試したが、きちんと虎刈りになってしまった。ちょっと値が張るが、これからずっと使うわけだし、返品して heavy duty のものを買ってこよう。
- 当のトラ・モンはまったく気にしていない様子だし、明らかに涼しくなっただろうから、それはそれでよし、と。
  - 7/13/2008: 虎刈り？
- 7/4/2008: heavy duty に交換してきた。
- 今日の写真、氷嚢で遊ぶモンの図。
  
  (2008/7/1撮影)
  背中の虎刈りはちょっと見えないな。ちなみに、この氷嚢は数日後に穴があいた……。
自己表現について。
- というか、あれだ、「見てくれ」というか、ペルソナというか、自分をどう見せたいかということについて。
  - 8/13/2002: 研究者の経営感覚。
  - 8/19/2002: 手品師のような教育。
- 思わせぶりな発言（本質的なことは何も言わないでおいて）だとか、やればできる（けどやらない）とか、やればできる（と思って、とりあえずできたことにしておいて、良いの悪いのいう）とか、そういうことはしてはいけないな、というか、してもいいけど、みっともないな、と思った。広い意味での「はったり」なのかな。まあ、多かれ少なかれ、みんなしていることだとは思うけれど。
  - 2/2/2004: 自己宣伝が過剰な日記を見る気がしなくなった。
  - 6/8/2006: ハッタリをかまさないという生き方。
  - 7/28/2008: ハッタリではなくて、「有言実行」の有言として（かな）。
  - 10/3/2008: カズがあの年でケガが少ないということは、ハッタリでできることではない。
- そういうことをしないで、では、どうすればいいのか：
  - 「思わせぶりな発言」は、まず、しない。発言するなら、隠さずに言う。隠さずに言うとボロがでそうなら言わない。隠さざるを得ない内容なら言わない。
    - 定理「囲い込まなければいけないような情報には、即物的で小さな価値しかない」 (江島健太郎)
  - 「やればできるけどやらない」は、「自分にはできないんだ」と認識する。
    - 2/10/2007: 「できるけどやらない」は負け惜しみ。
  - 良いの悪いの言いたいのならば、自分できちんとできた後にする。できないのならば、言わない。
    - 評論家嫌い。（例えば4/28/2004）
- 以上が文句。つまりは、あれだ、「良いの悪いの」言っているわけだ。
- 「じゃあお前は自分を人にどう示したいんだ」ということだな。
  - 普通の（ハッタリ半分な人間の）主張を５割引きで聞くところを、割引きなしで（せめて１割引き程度で）聞いてもらえるような、そういう人間になりたい。信用してもらえる人間になりたい、ということだろうか。
  - さて、そうするためには、具体的に、どうすればいいのだろう？（きちんと、日々、飢えることなく生きていくと同時に。）
  - 8/24/2008: 「これまでの研究業績」という自己解説。
- 3/18/2007: 対人インターフェイス -- 関連事項として。

6/29 (Sun)

[music] ふと魔が差して、とうとうアンプをポチッとしてしまう。
- 一応、うちの人の許可はとった。誕生日プレゼントだって。はやく来ないかなあ。しかし何か中学生のガキみたいだな。
- 7/2/2008: 結局、直接、店に行って買ってきた。
[music] maceo parker band のライブにでかける。
- 今年の eijf のトリを飾る、 maceo parker. 6/24/2008にオーダーし、 6/27/2008にピックアップしていたもの。
- 昨日 youtube で勉強したら、 JB のバンドに居た人で、ジャズというよりも完璧にそっち系の人だ。
- 場所は whyte 沿いなので、最初、バスで行くかと言っていたけど、時間（日曜の１０時開始）を考えると、どうも帰りの便がなさそうなので、車で行く。幸い２ブロック程度はなれた道路に停められた。会場は、 baltimore の funkbox （9/29/2003）の二階テラスなし、という感じ。 dance floor とそのまわりを囲むように小さなテーブルがちょこちょことある。
- 最初、テーブルに陣取ったが、始まってすぐに、もっと近くで見たい、聞きたい、踊りたい、と言う人の強力な要請により、椅子とテーブルに別れを告げて、前方に移動。 first set を耐える。何というか、まあシャイなんだな、あとまわりの目をやっぱり気にしてしまうんだな、体は多少揺すっていたけれど（そうしないと、逆に目立つ）、ちょっとつらいものがありました。休憩で私は水、もう一人の人はビールを飲んで、後半は私は後ろで休ませてもらうことにした。
- ジャズフェスの volunteer の人たちが「ご苦労さま」の意味で招待されていたのかな？そういう（年が我々よりもかなーり上の）人たちも結構いて、その休憩中になんか抽選会みたいなことをやっていた。そういう感じで、観客はいろいろで、人目を気にせずに自分なりにみんな楽しんでいた。（この「人目を気にせずに自分なりに楽しむ」という姿勢が、なんというか培われていないんだな、自分に（というか、おそらく日本人に）。）
- 演奏は、何というか、そうそう guitar の人はよかった。あと maceo parker は、 enternainer なんだな。一曲目、みんな何かまだ乗れていない状況で、ステージから消えたと思ったら floor の最前列をサックス片手にカニさん歩きで横切り、盛り上げていた（candy dulfer がテーブルの間を歩き回ってブロウしていた光景を思い出した）。しかし正直なところ、しばらく funk 系は聞かなくてもいいや、という気持ち。多分、個人的には david sanborn （8/20/1999）とか candy dulfer （5/25/1998）とかがギリギリの線で、今回の maceo parker はボーダーのあっち側だったみたい。
- 7/7/2008: この日手配した別のチケットが届く。

6/28 (Sat)

[monton] 最近のモンちゃん。
- 一階すべてを開放されたモンちゃんだが、実は最近、少しトイレの精度が悪くなってきている。枠から少しはみ出ている。（6/14/2008）ウンチは前からはみ気味で、まあいいかと思っていたのだが、最近はおしっこも。前足は入っているようなんだが、後ろ足が出ることがある模様で、しかも最近、その頻度が高くなってきている。
- ということでダンボールで、トイレに高さ５センチくらいの外枠を付けてみる。どうなるかな？
- 7/23/2008: おしっこ問題でモンちゃんを病院へ連れていく。

6/27 (Fri)

[music] 出掛けにチケットをピックアップしてきた。
- 6/24/2008にオンラインでオーダーしていたもの。
- 6/29/2008: ライブ当日。
[phys]ryuon
- 高速な implicit algorithm の実装：
  - テストのテストの続きの続き。
  - ３次方程式の解法。
    - 自分のコードにやっとこさ間違いを発見。正直に書くと、 x^y の計算で、結果を代入していなかった。つまり
      x = pow (x, y);
      と書くべき所を
      pow (x, y);
      としただけで、のんびりしていた。 bug なんてこういうものだよね、という教訓のために明記しておこう。
    - その後、 GSL と速度比べしたら、あれれ倍近く遅いや。
    - GSL のコードを見ると、きちんと戻り値で解の数を返しているので（これが欲しかったのだ。しかも後で見たらきちんと document にも書かれているではないか…… 手を抜かずにきちんと読みましょう）、 GSL のルーチンをそのまま使うことにする。
      - 7/26/2008: 挙動がちょっと不安定？
    - 一つ一つ bug を潰して行って、やっと６種類で妥当な結果を得る所まできた。勝率１０割！！
      - FENE-Fraenkel では s を大きくとりすぎていて、おかしい（精度が悪い）としばらく悩んだ。論文に従って s=0.01 程度にしておく。
  - 高速 scheme と普通の非線形方程式解法（NITSOL and GSL）の比較：
    - これもなかなか骨が折れた。
      - ポイントは configuration を作る際、 FENE モデルでは bead distance を Q₀ よりも小さく取らないといけないが、そこをしくじっていた。
    - 速度を比較したところ、本当に高速 scheme がダントツだな（これを期待して頑張っていたわけだけど）。
      - N=3 だと少なくとも３倍、全般的に５倍程度は出ている。
      - N=30 だと GSL が急に performance が悪くなる。それでも NITSOL の３倍弱は出ている。
      - 調子にのって N=100 も見てみると、いくつかの場合 NITSOL がこけるなぁ……。
    - 6/30/2008: noHI のテスト、終了。

6/26 (Thu)

[phys]ryuon
- 高速な implicit algorithm の実装：
  - 昨日の続き、テストのテスト。
    - 6/27/2008: さらに続き。
  - 全然テストになっていなかった（つまり、勝率５割なんて大嘘だったということ、やれやれ）。
  - 一般的な非線形解法の debug は、終わったかな。
  - 一方、高速 scheme は全然ダメであることが判明：
    - 調べてみたら、３次方程式が正しく解けていない。
    - 自分でコードを書いているのだが（6/24/2008）、そいつがうまく動いていない模様。参ったね。
    - 結局、これで一日潰す。しかもまだ動いていない……。
    - 帰宅後、 GSL を覗くとあったよ。 gsl_poly_solve_cubic(). こいつは（当然）きちんと解けているな。

6/25 (Wed)

[phys]お勉強
- NLD 問題：
  - 気分転換に、時間を無理やりここに割く。
  - 理論の論文にざっと目を通す。よく分からない。
  - 差分法のおさらい：
    - 以前一瞬やってる（9/21/2007）。その時のテスト・コードを動くようにしておく。
    - 1D は動いていて、 2D も組んではいる模様。動くのかな？
[phys]ryuon
- 高速な implicit algorithm の実装：
  - 昨日の続き。
  - 非線形解法のレベルでテスト。
    - ３勝３敗、勝率５割だ！（悪いな -- 内二つは線形だし。）
    - 6/26/2008: 勝率５割なんて大嘘だった……。
    - 6/27/2008: さらに続き。
ここ数日、二人共体調を崩している。
- 風邪で、かつ、のどが痛い。
- 前回の風邪の時に薬が切れていたので、帰宅時に買って帰る。

6/24 (Tue)

[phys]nanoプロジェクト
- dWLC の simulations:
  - CF-cylinder:
    - 先週はじめた計算、予定通り終了（ってことは、７日）。
      - もう一本は、この調子だと倍ちょっとかかりそうだな。
    - 結果は、ふむ、 d を大きくとりすぎていて面白くない（し、分布がきれいではない == 相空間を十分走り回っていない）。
    - 計算の追加を流しておく。これも一週間の予定（１００％で）。
    - 7/2/2008: 予定通り、終了。
[phys]ryuon
- 高速な implicit algorithm の実装：
  - さて、今度こそ本当に実装するぞ！
    - 作業場にて実装開始。
    - まずは３次方程式の解から……。えっと、千里の道も一歩から、と。
    - 6/26/2008: バグっている。
  - 丸一日かかって（家に持ち帰って）、何とか書き上げた。
    - とりあえず straight chain の決め打ち
      - また、である……後で一般化するのが、また大変なのは分かっているが……。
    - あと、とりあえず noHI だけ。
    - ただし force law は ILC 以外６種類、一通り組みました。
    - 対応する非線形解法アプローチも GSL multiroot 版と NITSOL 版とも組んだ。
  - 今のところ、全て未チェック。
  - 6/25/2008: チェック開始。
[music]またチケットを手配。
- 文字通り昨日の今日ではあるが、いろいろな意味で正しいお金の使い方、ということで。ここ（edmonton）ではこういう機会を逃すと何もないし。
- 6/27/2008: ピックアップしてきた。
- 6/29/2008: ライブ当日。
[phys]一般相対性理論を勉強してみよう。
- と言ってから数日が経過した。幸い（なのか残念ながらなのかは知らないが）刺激剤というよりは入眠剤になってるようで、ベッドでパラパラ眺めるにはちょうど良い。
- しかしこの内山さんの訳が、ちょっと辛いなあ。日本語がという話では当然なくて、訳語がドイツ語系なところ。ふとみると英語版は dover から出ているみたい。そっちで見るのがいいかな？
- ちょっと調べてみたら、こいつ文庫になっているみたいだな（amazon.co.jp: 上、下）。訳もそのままなのだろうか？
- ちょっと調べてみたら^２、「第五百七十夜 2002年07月01日アルバート・アインシュタイン『わが相対性理論』」@senya(lc)を発見した。
  - まずは些細なこと。引用：
    ボルツマンは「優美にすることは靴屋と仕立屋にまかしておけばいい」と言った。
    - 私の記憶にあるのは maxwell が boltzmann の書いたものを評して言った「彼の書くものは長すぎて困る」というようなもの（出典はなにだったっけなあ）。これを見て「いかん、俺は boltzmann タイプだな」と思ったので覚えている。
  - で、本当にドンときたこと。引用：
    では最後に自慢話をひとつ。
    　ぼくは自分が科学にまったく携わってこなかった者としては、いちはやくアインシュタインの相対性理論を理解できた数少ない一人だったのではないかと思っている。これは湯川秀樹博士とのちょっとした蜜月が支えになったからなのだが（このことについてはまた書くことにするが）、明けても暮れても相対性理論に集中した濃厚な日々も役立っている。
    　ずっとのちになって、当時、大阪大学理学部教授だった内山龍雄さんが岩波新書に『相対性理論入門』（1978）を書くことになったとき、内山さんから「どういうふうに書こうか悩んでいるんで、相談にのってよ」と言われ、二人であれこれそのシナリオを練ったことがあるのだが、そのとき内山さんは「君はどうしてそんなにおもしろく相対性理論のことが考えられるの？」と聞かれた。正確ではないが、こんなふうに応えたのではないかとおもう。
    　「アインシュタインがおもしろそうに考えたということばかりが気になってましたから」。
    - いや本当に恐れ入るな、松岡さんには。「いちはやく理解できた」ことにではなく、「本物」な人の懐に飛び込む、その行動力なのか、人間的な魅力なのか、そういうものについて。
    - 4/17/2008: 本格的で本気な人に沢山会ってきた人としての松岡さん。
    - 3/3/2008: 松岡さんのアンテナ（シェイクスピアのレコーダー）。
  - ついでに、千夜千冊の今日のエントリー「第千二百四十九夜　2008年6月25日三枝充悳大乗とは何か」@senya(lc) からも引用：
    かつてぼくは「炭男」（すみおとこ）と称していたことがあるんだね。「半巡通信」というたった８ページのパーソナルメディアを、最初は３００人に、最後は４０００人に無料で送っていたころのことです。７年ほど続けたかな。その「半巡通信」のなかで、「ぼくは炭男だ」と書いた。その炭男が発火すると火元になる（笑）。
    - 考えてみると、この日記が私にとっての「パーソナルメディア」なんだな、「無料で送って」いる。しかし私には、届けたいと思う明確なターゲットがそもそも「３００人」もいないけれど。
    - そういえば、自分がやりたいことをやるメディアとしての「雑誌」というのは大塚英志をみてふうんと思ったし（4/3/2003）、「ベンチャー」は近藤淳也をみてふうんと思ったし（9/5/2005）、「（私）学校」は西郷隆盛をみてふうんと思ったな（12/27/2006）。
    - 8/4/2010: [blog]とりあえずぼくはロケットになろう@kenmogi
- 7/8/2008: 「連星パルサーから一般相対性理論を実証」@/.jp.

6/23 (Mon)

[phys]ryuon
- 高速な implicit algorithm の実装：
  - この際なので、いろいろな force laws についても、まとめておく。
    - Somasi et al (2002) は、いわゆる Werner spring,
    - Hsieh et al (2003) は Cohen's approximation,
    - Hsieh et al (2006) は FENE-Fraenkel を使っている。
    - ILC 以外は、代数的に計算できる。
  - ここまでをノートにまとめて、今日は終わった……なかなか進まないな。
[music]ジョン・スコのライブ
- のチケットをピックアップしに、出勤途中にシアターに寄る。が、 box office はお昼からだった……。先週の box office のスケジュールは、朝（9:30）から開いていたのになあ、スケジュールがその日ごとに変わるとは思っていなかったぜ。
- で出直して、４時すぎに再びピックアップしに行く。 box office は工事の関係で３階に移動していた。ということで朝と合わせて、駐車場に２ドル×２を払ったことになる。で、お金がもったいないなあということで一瞬バスで行こうかと言ったりもしたが、やっぱり不精者なので車で行く。駐車料金の上限が１５ドルだと表示にあったのを見ていたのだけど、行ってみると実は６時以降は３ドルだった。バスだと２．５ドル×２人の往復だったので、車で行ってよかった。
- 本題のジョン・スコ。会場が階段状のシアターで（しかも結構こじんまりとしたもので）、後ろから２列目というのも全然問題なかった。というか、高さ的にも、方向的にも、ちょうどよい角度であった（ジョン・スコの正面、という意味）。前半、 bill stewart のドラムのセッティングが定まらず、というか、演奏してたら、どんどんバスドラが前にすべっていって、困っていた。スタッフに助けを求め、それからは大丈夫だったみたい。しっかりドンドン踏んでいた。あと bill の頭のてっぺんにピンと立った髪が気になった（多分、そういうセットなんだろう、きっと john への当てつけだ）。それから steve swallow が集中してジッと john を見つめて弾いている様は、実は家のモンちゃんが「遊んで！」と本気で顔をグッと寄せて迫ってくる所そっくりだった。彼はピックでベースを弾いているのかな？右手が上下に動いていた。普通（なのかどうかはよく知らないが）みんなはジャコ風に指で弾いてるような気がするので、なんとなく新鮮であった。まあ見た目の話。出てくる音でそんなことは（少なくとも私には）分からない。
- とかいう話は、どうでもいいことで、演奏の方。いつもの "a go go" からの green tea とか、スタンダードからは私の好きな cole porter の "everything i love" を演った (cf. 6/16/2000. まあ正直、やっぱり bill evans の演奏の方が好きではあるが)。あと ray charles の曲とか、カントリー風の曲（カポ付けて弾いていた）とかもやったけど、これらはちょっと今一だったかな（個人的には）。あとサンプラーを駆使しようとしていて、こちらがドキドキするような瞬間もあった。まるで joe zawinul を見ている気分だった。 MC でも言ってたけど、後で amazon 見たら新しい CD が出ていた。そこからの曲がその手のものだったのかな、 "Satisfaction" と "The Low Road" という曲。個人的には、ファンク系というの？そういう奴をバリバリ弾いているのが、やっぱり好だな。
- アンコールは一曲、例の "en route" から wee を短めに演って終了。何かちょっと物足りなかったかな。しかし、やっぱりライブはよい（５年ぶりだな、ふう）。
- 外に出たら、まだ明るい。小腹も減ったので、何か食べて帰ろうか、とちょっと china town をうろつく。が、明るいので誤解してしまうが、気がついたら夜の１０時近くで、開いているのは subway くらいだった。とぼとぼとシアターの駐車場に戻ると、道で知り合いにばったり出会った。さすが大都会の edmonton. ちょうど quebec への旅行に出かけるところだった模様。
- 6/16/2008: チケットを手配。
- サマリー： john scofield
- 6/24/2008: 昨日の今日ではあるが。
- 9/20/2010: Yoshio Suzuki x 3 Pianists @ kanazawa

6/22 (Sun)

今日はすごく曇り空。
- 昨日の天気はついていたと思った。
- 食料が尽きてきたので夕方、買い物に出かける。
- どうも肉体的に昨日の疲れが響いているようで、８時すぎには本当に睡魔に襲われたので、何も考えずにすぐに寝た。

6/21 (Sat)

a day trip
- 家の人が最近参加している、とある音楽関係（？）のグループでバスをチャーターして jasper に出かけるというツアー。ここに来て（ロンドン@オンタリオも含めて）観光なんてどこにも行ってないので、もちろん初 jasper である。
- 朝７時集合、出発で片道５時間のドライブ、夜１０時帰還予定という強行軍（maps.google.ca）。予定通りに予定は遅れ、結局出発場所に戻ってきたのは真っ暗になった１１時半であった（一年で一番日が長いシーズンなので９時すぎなら楽勝で明るかったのだけど）。
- 立ち寄ったのは３ヵ所:
  - maligne canyon@explorejasper(lc)
    
    軽い昼食の後、上のレストランから下の 5th bridge まで、岩の間を流れる川に沿っての trail を歩く。途中、雪（氷）が残っているところが２ヶ所あった。後半の道が急にさみしくなって、正しい道を歩いているのか不安になった。（他の人もそう思ったみたいで、後続のグループでちょっとした行き違いがあった。）３０分程度の軽いハイキングということだったが、（この時は）晴れて日差しがあったためか、予想の外疲れた。
  - jasper town
    
    （知らなかったのだけど）ちょうど今日が national aboriginal day だったようで何かイベントが行われていた。ムースの皮の何やらをやっていた。街自体は、いわゆる観光地のお土産屋が並んでいるだけ、という感じ。どこでも食べることができるネスレのアイスを食べた。
  - athabasca falls@explorejasper(lc)
    
    最終ポイント。水量があって、結構、豪快な滝だった（落差はそんなになかったけど）。もう疲れていたので、我々は下までは降りずに帰ってきた。その後、軽い食事をとって、帰路に就く。ちなみに、食事は世話人の人が準備してくれた、いわゆるピクニック・ランチ的なもの。それが、どれも、おいしかった。
- 集合解散地点の駐車場の状況が不明だったのでバスで行ったのだが、帰り着いたのが１１時半。ちょっと不安になったが、ここのバスは（特にメインの路線は）結構遅くまで走っていて助かった（何でも２時くらいまであるとのこと）。１０分程度の待ち時間で無事に来て、乗り継ぎも待ち時間なしで、無事に帰り着いた。
- 実は出発時点から（そして予定が遅れる度に）心配していたお留守番の子。家に着いて真っ暗な家の中に向かって「ただいま」と呼びかけると、パタパタと尻尾を振る音が聞こえてきて、やっと安心した（モンはほとんど吠えない子なので）。一通りチェックすると、準備しておいたご飯（ドッグフード）には一切手をつけていなかった。しかし我々が荷物を置いたりして一息ついていたら、急に食べ始めて、あっという間に空にした。モンちゃんは結構、神経の細かい子らしい。これからは、モンちゃん込みで色々と出かけることにしよう。

6/20 (Fri)

[phys]お勉強
- 昨日のセミナーに関連して、理論について、および数値解法について、いくつか文献を印刷。
[phys]ryuon
- 高速な implicit algorithm の実装：
  - つまり Hsieh et al (2006) で使われている Somasi et al (2002) および Hsieh et al (2003) のこと。
  - 振り返ると、結構時間を費やしているな（毎度のことだが）。
    - 5/7/2007: 論文の印刷。
    - 10/12/2007: 勉強中。
    - 11/16/2007: 依然、勉強中。
    - 11/29/2007: まだまだ、勉強中。
  - ここ数日 bead-rod を陽に考えたおかげか、今度こそ、頭に入ってきた。
  - 揮発性メモリーから消える前にノートにまとめておく。
[phys]一般相対性理論を勉強してみようかな。
- と書くと、すごくチープに見えるな。まあ、しかし当然、そういう青臭い話ではなくて、地味な地道なお話として。
- 一応、重力方程式とか勉強した身なのだが（学部の頃、たしか単位は取ったと思う）、当然もうすっかり蒸発してしまっている。
- その時お世話になった「場の古典論」を探しているのだが見つからない。置いてきたか？ということで、かの wolfgang pauli が若干２０歳そこそこで書き上げたという本（内山訳）を（いや、なぜかこっちは手元にあるんだな、講談社から出てたハードカバーが）しばらく眺めてみることにしよう。
- なんでまた、ということについて。直接的には Oettinger (1994) の話がさっぱり分からなかった（6/18/2008）ということがあるのだけど、多重極展開を遊ぶのにもちょうどよいかな、という想いが半分。とはいえ、原因とか動機とかはどうでもいいことで、肝心なのは、やるかやらないかにある訳で。まあ、ぼちぼちといきましょう。
- 6/24/2008: その後。

6/19 (Thu)

[phys]ryuon
- bead-rod の勉強：
  - もうちょっとお勉強、パート２（きちんと、ぶっつかる所まで行っておく）。
    - なんというか、諦めの悪い性分だなあと思いつつ。
    - 具体的に使ってる仕事を探してみる。あんまりなかったけど、一つ見つけた。ううむ、確かに、しち面倒くさいことをやっている（やる必要がある）みたい。
  - やっぱり、尻尾を巻いておこう（少なくとも、しばらくは）。
    - うむ、確かに諦めの悪い性分だ。
- bead-spring for bead-rod の話：
  - Hsieh et al (2006) の FENE-Fraenkel を実装してみよう。
    - 第７の spring type を導入（既に６種類も実装している、ということ）。
    - こいつは３つパラメータが必要なので、 struct bonds に手を入れる。
    - guile parser も update.
  - 今日はここまで。
- 10/2/2008: 再始動。
[phys]グループ・セミナー。
- 随分と久しぶりだな。（5/2/2008以来か？）
- 具体的なナノテクのプロセスに関連した話。やっていることは結局 2D の輸送方程式を解きましょう、という話。で、解いてるモデルの出所は Flory-Huggins と。
- みんながどこに立っているのか、時々、よく分からなくなる。一方では分子の詳細を見ていて、それが見えないとつまらないという世界にあり、一方では現象論的な free energy をドンと持ってきて、それを解いている。私には、このレベル（階層）の違いがとても大きく感じられる（というか、つながってないと感じるし、大きなギャップがあると思う）。それぞれはそれぞれ、きちんとした仕事をしているけど、このグループ全体で、この両者が、本質的な所で意思疎通できているのかなと、ふと思った。私自身が、その両方に不満を感じるような立ち位置にいるというのが、この違和感の原因の一つなんだけど。
- 6/20/2008: 関連する文献を取ってくる。
- 9/12/2008: 次回のグループ・セミナー
[SF][WaoN] gWaoN にちょっと手を入れてみる。
- 6/15/2008 の「鍵盤を描こう」という話。
- とりあえず（無理やり上書きで）書くところまではできた。
- スナップショット：
  
  モードは "PV-FFT" で、鍵盤の他に半音で平均された棒グラフも書いている（ちょっと見にくいね）。
- レイアウトをもっときちんと定義して、棒グラフの色を考えて、画面と周波数の関係をきちんと debug して、それから octave removal を spectrum に課したいなあと思っている。再生時に表示している current position の indicator が消えてしまうこと（処理が追いついていない）も何とかしたい。というか、描画全体に tune up する必要がある。 spectrogram を毎回計算しているのが、まずはおバカな所だ。とか考えていたら、急に息切れしてきた……。

6/18 (Wed)

車のヘッドライトが切れた。
- ここ（カナダ）では、昼間も点灯することになっていて、実際ここの車はエンジンをかけるとヘッドライトが点く。私の車はもともとアメリカから持ってきたので、ライトのスイッチを入れないと点かなかったのだが、先日のチェックで（alberta の plate を取ったとき）交換しなければいけなくて（cf. 6/4/2007）、今はカナダ仕様になっている。
- で、先日（２、３日まえ）、ふと、左のヘッドライトが切れていることに気づいた。切れているのは一つだけなので、 fuse の問題ではないみたい。昨日、マニュアルを引っ張り出してきて見てみると、電球（bulb）は自分で交換できるみたい。ということで、今日、出勤前に canadian tire に立ち寄ってくる。「９００６」という番号がマニュアルにあったので、それを頼りに探してみると、一つ９ドルの安物から２５ドルのものまでいろいろあった（３、４種類）。中くらいの１４ドルのものを買っておく。
- 帰宅後、ガレージにてゴソゴソと交換作業をする。この作業自体はマニュアル化されてルーチン化されていて、単純なものなのだが、ヘッドライトの裏の部分に手が入らず、結構てごずった。が、何とか作業終了。切れたはずの bulb を見てみるが、（ハロゲン・ランプなので）どこがダメになっているのかはっきりせず、一瞬、別の原因か？と冷や汗をかいたが、エンジンをかけてみると、きちんと点いた。めでたしめでたし。
[phys]ryuon
- bead-rod の勉強：
  - bead-spring for bead-rod に進む前に、もうちょっとだけ勉強しておこう、と。
  - 一つの key paper として Oettinger (1994) Phys.Rev.E を（これまたもう一度）おさらい。
    - 論文には結果の式は出揃っているみたいだけど、途中がすっとばしてあり、しかも何か難解。彼の本は手に入らないし、困ったなあ……。
    - 6/20/2008: そうだ、一般相対性理論を勉強してみよう！
  - とりあえず、式だけ書き出しておく。
    - こいつら、どう実装するのか、結構悩ましいな。
    - それに多分、 stiff な bead-spring がやっぱり正解みたいだな。
  - 6/19/2008: もうちょっと、パート２。

6/17 (Tue)

[phys]nanoプロジェクト
- お約束の alex と話。
  - 基本的に私の現状報告。あと論文執筆のプランなど。
  - 8/26/2008: 会議録論文の話。
  - 8/29/2008: alex が原稿を上げてきたので、目を通す。
- dWLC の simulations:
  - CF-cylinder:
    - 先週仕込んだ結果、もう一つも終了。（ってことは７日か - 6/10/2008）
    - 結果を解析。平衡に行っているかどうかが（いつもの）懸念。前みたいに、現状ではまず、下限に注目しておくことにするか。
    - 追加計算を仕込んでおく。（予定では一本は１週間、もう一本はそれ以上。）
      - 6/24/2008: 予定通り、一本終了。もう一本は、倍ちょっとかかりそう……。
[phys]ryuon
- bead-spring for bead-rod の話：
  - bead-rod 後の、世間の流れとして。
  - Hsieh et al (2006) J.Chem.Phys. に（もう一度）目を通しておく。
    - 初見（11/5/2007）ではスルーしていた部分にある、お話の背景が分かってきた。
    - ということで、またこいつを実装してみようかな。（今回は、すでにパーツは全部そろっているはずなので、簡単だろう。）
    - 6/18/2008: 進む前に、もうちょっと。
夕食。
- 知り合いと食べにでかける。前に我々が呼ばれたことがあるので、本当ならば家に招待しなければならないのだけれど、なかなか家が片付かないので（責任は我が家の構成員全員で１／３づつかな）、どこかに食べにいこうということになったもの。彼らのスケジュールが結構たて込んでいる関係で、急遽、今日になった。別の知り合いを一人さそって（車の定員一杯）、また little india 辺りに行く。
- これまた別の知り合いに聞いた、おいしいインド・レストランを目指したのだが、結局見つけられず、先日の店と同じモールにある別の店を試す。私的にはそんなに違いは感じられなかった（ナンは今回の店の方がふっくらしいて、それが authentic なのか大衆向けなのかは知らないが、おいしかった）が、かお曰く前の店の方がおいしかったそうな。少なくとも料理のバリエーションは少し負けていた。
- 今日は夕刻、強烈な雷雨が来たのだが、それも上がって、帰る頃（９時半）には晴れていた（まだ太陽は沈んでいない）。
[job] しばらく前に出していたロンドン（オンタリオ）の公募の結果が正式に届いた。
- 残念であったという結果は、既に先月、こちらから問い合わせて informal には聞いていたけれど。
- UWO での思い出は必ずしもよいものではなかったけれど、物理だったし（一応、出身は物理なので）、ナノだったし（現職がナノテクの研究所なので）、よい意味での知り合いも居たし、「書類がとどいてませんよ」という知らせもあったし、結構、期待していたのだけど。

6/16 (Mon)

[news] 今月から何か紙がよくなってるAPS News から。
- "Physicists and Copyright-How to give away your birthright for what? (W. G. Unruh)"(lc)
- ほほう、Royal Society of London journals は copyright を要求しないのか。知らなかった。
- まあ response にあるように APS はかなり著者に自由度を与えていることは知ってる（web に載っけてもいいとか preprint server に置いてもいいとか）。
- この点、 AIP の Phys. Fluids はちょっと保守的だ。 cf. 10/30/2003
- こういう普段、あからさまな場所で話にのぼらないことを時々きちんと意識することは、大事なこと都だろうという自分への remind. （でも、普通は、 copyright transfer form はホイホイとサインして送っちゃうよな。）
[monton] 今日のモンちゃん
- 今日は注射の日。
- 暴れ度合いが減って、ずいぶんと大人になったと思う。それでも診察台にのせられたときには、相変わらずチビッたけど。獣医さん曰く「喜んでいるのではないか」と。めずらしいそうだ。
- 5/12/2008: マイクロ・チップが埋め込み、および注射。
[phys]nanoプロジェクト
- dWLC の simulations:
  - CF-cylinder:
    - 先週仕込んだ結果、一つが終了。（ってことは６日か - 6/10/2008）
    - 結果を解析。落ち着いたかなと思ったら、最後の最後にまた動き出している様子……。 endless だなあ。
    - 現状では下限を見て、そういうものだ、という結論で置いておくのがいいのかな。
      - 6/17/2008: 似たように……。
[phys]ryuon
- bead-rod の実装：
  - ここ数日、モヤモヤと考えていたことを整理。
    - 幾つかの高速化のアイデアと、取り扱いの改善案（SCM での設定項目を増やすという話）。
３時のおやつ。
- もう随分と長いこと一日２食の生活をしている（つまり、基本的に、職場で食べ物は口に入れない）。
  - 時々変な顔をされることもあるが、逆にそんなに食うと太るよね。こういうスタイルは、下手するともう１０年以上になるのではないかな。
  - まあ、いずれにしても、一日一食の森さんには負けますが…… って、そういう競争ではないね。流石に、一食には出来ないな、まだ若いし。
- 出勤時に alex に「午後、部屋にケーキを食いにこい」と誘われたので、ちょっと立ち寄る……つもりが、まあ終業時間までおしゃべりをして、一日終わったよお（そうなるだろうとは思っていましたが）。
- 明日、仕事の話をしましょう、と。
- ケーキは、自家製セロリ（別の人の家の庭でとれたもの）を使った「セロリ・ケーキ」（作ったのも彼ではない）。うまかった。
最近ここもガソリンが随分と上がりましたね。
- １リットル当たり、軒並み 133 セントのあたり。いつも入れている safeway の gas station は、買い物で７セントの割引きなので今日は 1.259 ドルで入れたことになるが、見れば去年の１０月は 0.875 ドルとかだったので、まあ随分と上がったことになる。
- ここでは通勤に使っているので、だいたい月に一回の給油という感じで、４０ドルから５０ドルといったところ。
[music] ライブのチケットを手配。
- 実は、カナダに来てはじめてかな。この週末から、いわゆる edmonton international jazz festival(lc) というのがあって、そのメイン・アクトの一つとして john scofield trio が来る(lc)ので、それをとった。一人４０ドルなり。
- 何やら
  (his) music falls somewhere between post-bop, funk-edged jazz, and R & B.
  とある（し、一般論としてそうなんだ）けど、この面子（Steve Swallow - bass, Bill Stewart - drums）で、それは無いでしょうね。このトリオの CD "en route" は持っているな。あんまり愛聴盤ではないが。
- john sco はこれで３度目だ。（9/29/2003 at baltimore, 4/24/1999 at LA）どっちも MMW というかそっち系だったけど、今回は jazz な john sco が聞けるだろう。席が後ろから２列目というのが、ちょっと悲しいが……（出遅れたので、仕方ないか）。
- 6/23/2008: ライブ当日。

6/15 (Sun)

[blog][music][WaoN] "Kurt Rosenwinkel - Untitled Improvisation"@youtube.
- これは"Kurt Rosenwinkel WBGO Studio Session - On NPR Music"@wbgo/blog(lc) および "Kurt Rosenwinkel: Jazz Guitar Alchemist"@npr(lc) だ。誰か採ってくれないかなぁ
  - 本当のはじめの部分（下の絵の左１／３くらい）：
```
E-|-----------------|---------------|------------|---------------
B-|-8-7-6-5-4-3--3--|--6------------|---7--8---?-|-??--?-??---?--
G-|-9-8-7-6-5-4-----|-7-757---7-7---|-79-97----?-|-??--?-??---?--
D-|-7-6-5-4-3-2-----|---------5-4---|----------?-|-----?------?--
A-|-8-7-------------|---------------|----------?-|-----?------?--
E-|-----6-5-4-3--0--|-------5-----5-|--------6---|----8-----8----
   
```
  - えっと、まあ、このために gwaon や pv があるんだけれども、ギタリストじゃないので、やっぱりコードが難しい（見えない……）。
  - 上の spectrum の下にも鍵盤を書いた方がいいな（というか、この状態じゃ単なる絵に過ぎないな。）
    - 6/19/2008: とりあえず描いてみた。
  - gwaon で遊びたい人は、 FreeBSD なら ports になってるし（12/16/2007）、 Ubuntu な人にも package が既にあるはずだし（4/19/2007）、その他の Linux な人たちにも（よく分かってないで作った） rpm が SF の download section に置いてあるので、どーぞ。
- それから上の WBGO/BLOG をさらに読んでいて、スゴイのを見つけた。 "Giant Steps"@youtube. 作った人は暇人だなと思うが、一度こうやって作ったものは今後ずっと人々を楽しませてくれるわけだから、正しいヒマの使い方だと思う。意図的なものではなく単にパクってきた譜面のせいだけど、省略されたピアノ・ソロに強く哀愁を感じた。
- 4/20/2011: Brad Mehldau の「Don't Be Sad」の動画。
- 12/25/2011: GarageBand で Brad のイントロの可視化。
- 12/29/2011: WaoN 0.10 デモ・ビデオとして Kurt のソロギターの transcription を作ってみた
[news][music] "Jazz champion Svensson dies at 44"@bbc(lc)
- まだ４４歳。信じられない。 wikipedia にはスキューバをしていて、とあった。スウェーデン人はタフにちょっとやんちゃしていたんだろうか、と勝手な想像。
- EST は最近知った（と思って見返したら、最初に言及しているのは 2004 年か、あんまり最近じゃないか）。 CD は持ってない。
- 個人的なベストは（実際に jukebox を確認しないで記憶だけをたよりに考えると） 2003 年の jazzbaltica での nils landgren の vocal が入った（pat metheny 込みの） "Believe, Beleft, Below" だな。（調べてみたら、このライブは一部は youtube でも見れるみたい。例えば e.s.t. Believe Beleft Belowとか。）
- サマリ：訃報
wanderlust
- 英語のメールは web で書いても構わないのだが、日本語はやっぱり emacs 上で書きたいということで、未だに wanderlust を使っている。（今時のみんなはメールをどうしてるのだろう？ web か thunderbird とかつかってんのかな？）
- ここしばらく POP server へ接続が（Reading failed というメッセージが出て）出来なくなっていたのだが、 .elmo/pop をきれいにしたら、正常に戻った。

6/14 (Sat)

[phys]ryuon
- bead-rod の実装：
  - やるならやっちまえ、ということで、 HI 系に進もう。
    - まずはノートにまとめを書く。
      - 式の上では問題なくできた（当たり前だけど）。
      - 一部詰めていない点は残っているが、とりあえず置いておく。（例の gradient 項の話）
      - no-HI で曖昧にしていた点（polydispersity）も、この一般化で全てかたが付く。
    - 既存の no-HI 系とは独立に HI 系のコードを書き下ろす。
      - もちろん使えるパーツは使って、ということだが。
      - 計算手順（計算する部分の切り分け方）を少し変えたりする。
      - とりあえず書き上げた。
    - チェック：
      - これがまた結構、手間取った。まあ今回は no-HI と one-by-one で詰めていけるので、比較的楽ではあったが。
      - 今回も、結局は下らんことであった。つまり（またしても）変数（vector）の入出力における上書きの手順の間違い（本当は、頭の中で「こう組んだ」という思い込みと、実際の実装が別であった）。トホホである（なんて、誰も言わんな）。
    - 今日は、とりあえず、目処が立った、というところまで。
  - ちなみに今の HI 系の実装は、本質的に Petera, Muthukumar (1999) と同じ（はず）だが、 connector vector に対してではなく、 position vector に対して組んでいる。
    - 6/8/2008: bead-rod の古典的な論文を読む。
[monton] 最近のモンちゃん。
- 食欲も戻ってきたみたいで、だんだんいつものモンちゃんになってきた。
- トイレも完璧になって（実は時々、はみションしてたりするけど……）、とうとうリビングも（昼間は）全面開放された。
  - 6/28/2008: はみション対策。
  - 7/24/2008: リビングを一時撤退。
- が、甘えっ子街道をまっしぐらに進んでいるみたいで、広い自由空間があっても、おにいちゃんかおねえちゃんのどっちかにペタッとしてないといけないらしい。誰かさんが時々（下で）一緒に寝ているせいではないかと、私は疑っている。誰かさんは、拾い食いも巧妙に調教したり、なかなかの trainer だからなあ。
- ご参考までに、拾い食いのトレーニング方法：
  - 犬の嫌いな匂いスプレーを染み込ませた食べ物（例えば切干大根とか）を準備する。それをさりげなく床の上に落として、犬に食べさせる。
  - （予定では）犬はそれを咥えるが、嫌な匂いが付いているので吐きだし、以後、床に落ちているものは食べなくなるはずだったのだが……
  - 注意点：犬にきちんと味覚が付いてからやりましょう！「あれ、食べちゃった」と言って、さらにパワーアップしたものを準備して、同じように落として、と繰り返していると、逆に、落ちたものを食べるという行動がどんどん強化され、家のモンちゃんのような立派な scavenger になります！！台所の友にお困りのご家庭のみなさまは、ぜひお試しを。

6/13 (Fri)

[phys]ryuon
- bead-rod の実装：
  - 一般化されたコードの debug:
    - guile parser のテストコードを書く。テスト終了。
    - constraints がきちんと満たされているかを確認。こまかい bug が一通り落ち、機能するよになった。
    - stokes3 で短い chain を走らせるが、おかしい。一つ一つ押さえていって、 position の更新順序の間違い（上書き）を発見。お粗末。
      - 6/14/2008: 似たような間違い。
  - これで bead-rod 系が走るようになった。ただし no-HI 系のみ（つまり Liu 1989 ということ）。
    - 既に定式化されていたとはいえ、 6/9/2008に思い立って一週間で動くものができたので、よしとしよう（定量的な検証は、まだ、これからだけど）。
    - 6/8/2008: bead-rod の古典的な論文を読む。
    - さて、勝負はこれからだ、と。「急がば回れ」が正しかったことをきちんと証明するためにはね、奥さん。（cf. 6/5/2008）
      - ちょっと無駄話。
      - 普通の研究者は、何か（その人にとって）新しいことをはじめるとき、どうアプローチするんだろう？風見鶏のように人のことばかり気にしてフラフラしているのはみっともないが、かと言って、人は人、自分は自分とまったく無視するのもバカだし、ということ。あさっての方向に走っていないかを確認するためにも、たまには人の振りを（意識的に）気にしてみようかな、と。
        
        9/29/2003: たまには人の振り見よう企画。
      - science ってのは、建前は、みんなが誰もまだやったことのない問題に取り組んでいるわけだけど、各人の学習プロセスで、果たしてどこまで降りているのかな。その深さが、実際にその人の発言の深さにつながるということは、経験上、だいたい分かっているが（浅い人はやっぱり浅い、という意味）。
      - 理論は、まあ式を追うというプロセスで最低限、自分の足で立てるところまで来れるとおもうが、計算をやっている人は、普通、どうしているんだろうな、と。例えば私の今回の「急がば回れ」は、理論で式の導出を追うごとく、既存の手法を自分で実装することで「自分の足で立っている」と思えるところまできたが、みんながみんなそうしているようには思われないので。
      - 理論でも、計算でも、（自分で降りることなく）テーブルの上に並んだ最新の成果をつまんで（花屋がブーケを作るが如く）というようなことをやると、まあ trend に乗ったそれらしいモノが一応できるのだろうが。計算で言えば、どっかから既存のコードを取ってきて、理解せずに使う、というのかな、よく分からんが。
      - 計算物理が理論と同様の、科学としての検証可能性を保証しようと思ったら、コードの公開は避けられないと思うが（RYUON マニュアル、おくづけ）、現状は、コードは個人（or 研究グループ）に所有されていて、何か新しいことをやりたいと思ったときに、やっぱりちょっと手が出せないし、それだけでなく、下手すると第三者による検証もできない、と。
[blog] 「長瀬弘樹という作曲家のこと」@kenn(lc)
- 以下、引用：
  長瀬「僕はもう、大衆芸術家の道をひたはしってるからｗ。でも、こびるのとはちょっと違うんだよね。また考え方が変わった。
  ……中略……
  　作曲家になりたての頃は、とにかく「大衆にうける」ってことだけしか考えてなかった。逆にいえば、大衆をばかにしてたんだよね。
  ……中略……
  　大衆というものを勝手に自分で設定して、ここで落としとけばいいやってやってて、でも、それが全然決まらなかった。３年くらいやってて、大衆なんてものを勝手に作ることは間違ってるんじゃないかと思い始めた。
  ……中略……
  　なんかね、ほんとに自分の心からでてるものじゃないなと実感するようになったの。これやればうけるだろうとか、そう思った時点で自分の価値観の範疇をでてるんだよね。
  ……中略……
  　大衆というものを客体化しちゃいけないということかな。
  ……中略……
  　自分が大衆にならなきゃいけないということだね、簡単にいえば。大衆としての自分が、大衆として、「大衆的」じゃないものを作る、そう思うようになってから徐々に認められるようになった。
  ……中略……
  　自分が大衆になって、今まででは売れそうもないものを作るというニッチの理論だね、誰もが言ってることだけど。
  
  江島「新人だと、奇抜さに走ると見透かされちゃうことが多いじゃない？目立つためにやったんだろうなって。」
  長瀬「だからね、それは目立つためだけに奇をてらってて、底があさいからだよ。
  　……中略…… なんていうかな、「あれ？これ、今までやっててもおかしくなかったし本人にすごくあってるけど新しいね」というところにもっていくということかな。なんでいままでやってなかったの？っていう。」
- 大衆つながり、ということで。
  - 6/12/2008: 大衆化で失われるもの。働くこと（by 森）は大衆化なのかな？
  - 6/8/2008: 働くということ（森博嗣）。
  - 7/15/2008: 知識と知恵とか、学習の効率とか、モヤモヤと。
  - 7/28/2008: 売れるかどうかは分からない、というのは経営感覚がないということかな。
- 「古い知り合いの成功譚」を読むと、その人がブレイクする以前の、書き手の態度がどういうものだったのかが気になる（のは、下世話なのかな）。
  - 何か条件反射的に、今まで存在も知らなかった親戚が急に名乗り出てきた、みたいなステレオタイプな光景が頭に浮かんできてしまう。
  - 例えば tatsuru さんも言っていたように（5/13/2008）
    貧しいとき、病めるときには、かつての知友は知らない顔をして通り過ぎてゆき、どこの家のドアも開かない。
    というようなこと。
  - 頑張っている人間への本当のサポートというのは、まあ難しいものではありますが。
    - おだてたり、甘やかすのが正解とは思わないけれど。
    - 「忙しい」と言い訳して今の優先順位を下げることは、結局、「何もしない」を積極的に選択していることと同じだと思うし（例えば、「現在にないものは永久にない（岡本太郎）」や「サマリー：忙しいということ」）。
  - けなすより、褒めるのは難しい（という当たり前のことは、日々実感している -- 例えばサマリー：教育について）。
  - 4/4/2009: 「有名人」が「無名人」であったころのことを綴ったものにだって意外にも多くのヒントが詰まっている(松岡正剛)。

6/12 (Thu)

[FreeBSD][coding] ATLAS-3.8.2
- The patch has been committed. (cf. 6/5/2008, 6/7/2008)
- See freebsd.org, freshports.
[phys]ryuon
- bead-rod の実装：
  - 任意の constraint への一般化：
    - bead から connector への変換だけでよい。
    - 実装は比較的 straightforward。予定通り、すんなりとは動いてくれない（動かないことが予定通り、の意）。
    - いつものようにゴニョゴニョして最後には動くところまで持ってくる。
    - で、結局 tridiagonal scheme （6/10/2008）は排除（straight のみにしか使えないので）。今は libiter を使っている。
  - コード（メインとテスト）を作業場から ryuon/libstokes へ移動。
  - 次は integrator:
    - 実はここ数日、頭の中に置いてモヤモヤと考えていた。で初手は、既存のスキームに後付けで補正を噛ます、という方向でいこうという気になっている。
    - 初期値と非拘束な最終値の両方が必要なので、上の階層で一括して、とはいかなかった。今のところ、各 scheme に処理を追加しておく。
  - guile の parser:
    - 書き上げた。
    - 途中、 BeadRod の初期化方法を変更。こういうのは、実際に使ってみる局面にならないと、何が正しいのか分からないよな。
      - こういうのを人は「頭が悪い」と言うのだろうな、と。そう言われたくない「頭の悪い」人は、そう言われないようにと、どんどん何もしない人になっていくのだろうな、と。だから、やっぱり川上弘美は正しいのだ、と。
    - stokes3 にも反映させる。今日はここまで。 integrator と parser は未チェック。
  - 帰宅後、 commit しておく。
  - 試しに、 planar flow （6/3/2008）を走らせてみるが、これまた予定通り、うまく動いていない（動かないことが予定通りということです。ってクドい。）
    - 軽く見返すと、いくつか trivial な間違いを見つけた。
    - しかし、これだけですんなり動くわけがないのが、 debug の楽しいところ（半分本音 -- サマリー：デバッグ）。
[blog][千夜千冊] 「第千二百四十七夜　2008年6月12日中西悟堂かみなりさま」@senya(lc)
- 引用：
  何度かその前を通っているうちに、ふと表札を見ると「中西悟堂」とある。引っ越して５年もたってからのことだ。それまではまったく気がつかなかった。ああ、あの鳥センセーだと懐かしく、よほど覗いてみるか、「こんにちは」とでも言ってみたかったのが、その勇気がないままになった。それからしばらくしての昭和６０年（１９８５）、８９歳で悟堂センセーは逝ってしまった。ぼくが表札を見たころは７７歳か７８歳くらいのときだったろう。
  松岡さんもこういうことをしていたのだな、と。
- 「日本野鳥の会」というみんなが表面的にだけ知っている組織の裏側に、こういう人物が存在していたのというのが、素直なおどろき。メディアって何でも表層化、矮小化、軽薄化させてしまうものなのかな。（まあ、それが大衆化なんだろうな。で森さんは「大衆化」を成した、と。（んな皮肉を言わなくてもいいのだけど、未だに何かうまく噛み砕けていないのだな……））そのプロセスで失われてしまう、しかし大事なものを感じさせる何かは、地味に（地道に）拾っていくしかないんだろうな。（4/15/2004）
  - 付記（6/13/2008）：全然関係のない話（でもないかな）。「立つ鳥後を濁さず」という言葉は、あれは下らないよな。置いていかれた人たちの側のエゴだけが聞こえてくる。「立つ鳥後を省みず」と言った松本眞さん（cf.6/1/2008）の観察の方が、多分、より真実に近いだろう。この「真実」も下らないと思うけど。こちらからは逆に、軽薄に、要領よく、都合よく、渡り歩く人間の側のエゴを感じる。
    - 6/8/2006: いわゆる「研究日記」について感じていること。
  - そうそう、松本さんのことわざ(lc)があれから２つ増えていたことを発見した。
  - 6/13/2008: 大衆芸術家の道、長瀬弘樹という作曲家のこと（kenn）。
  - 7/3/2008: わかりやすいネタに注目するな。
- ところで、「ロハス」ってなに？ wikipedia 様曰く、「LOHAS」。日本語の変な略語じゃなくて、 Lifestyles Of Health And Sustainability という英語の略語なのか。
  - なんて言っていたら、家の人に「今頃……」とか言われてしまった。まあ流行に乗り遅れるまいとアクセクする人間ではないのです。
- 以下、備忘録として。
- 「第千二百四十六夜　2008年6月06日マックス・エルンスト百頭女」(lc)
- 「第千二百四十五夜　2008年6月03日内藤湖南日本文化史研究」(lc)
[music] sufjan stevens の holland という曲。
- brad mehldau の最近のレパートリーで、私も結構気に入っている曲。てっきりオランダのことだと思ってたが、 michigan の街の holland なのだと michigan でのコンサートで彼が言っていた。
- で、 sufjan stevens って何者だろうと思ったら、いわゆるインディ系みたいで、家にいる詳しい人（特に下手な男性ボーカル系に詳しい cf. 4/7/2008）にきいてみたら、まさに彼女の favorites の一人らしい（で、以前、彼女が私に何曲か聞かしてくれたときに、何か悪く言ったらしい……）。
- 何でも彼は、全米５０州のそれぞれにアルバムを出す予定なんだそうだ。曲じゃなくてアルバム。まさに life work だな（年に一枚でも５０年！）。彼の出身が michigan で、まずは michigan は出てて、その中に holland もある、と。まあ、ツアーをやる musician なら、その土地々々の空気を感じることもできるだろう。
- youtube さまにお伺いを立てたら、いくつか引っかかって、オリジナルと思われる holland も聞いてみた（イメージ・クリップ付き）。 brad mehldau の演奏はオリジナルと同じ空気だった。
- そういえば、 beatles の有名な "black bird" という曲があるが、あれが私の頭に入ってきたのは、実は amazon の視聴で聞いた brad mehldau の "Volume 1" （これ、まだ持っていないんだな）での演奏（の、テーマの初めの半分くらい）だった。あれもアレンジ（というか、雰囲気）はオリジナルと同じだよね。そういう認識をもつ前に、実は既に、耳にタコができるほど jaco の "word of mouth" は聞いていたのだけど、あっちはずっと別物という認識だ。
- 同じ beatles の曲でも largo の "dear prudence" は、ちょっと雰囲気が違う。けど、こっちは、個人的には両方ありだ。つまり素敵な２つの曲という認識。
- 7/2/2008: sufjan stevens のCDを一枚、図書館で借りてきた。
- 7/15/2008: 今度は無事に "michigan" を借りてきた。

6/11 (Wed)

[news] 「水野晴郎、戸籍上の名前を改名していた！恐るべし『シベ超』への愛！今日都内で会見」@yahoo(lc)
- 私にとっての水野晴郎はゴールデン・タイムの映画（曜日は覚えてない）の解説者だ（淀川さんと荻さん、あと高島忠夫の時代）。うわさの「シベリア超特急」は見たことないので知らない。（日本を離れた後だったからだろうか？）
[phys]ryuon
- bead-rod の実装：
  - これまでのところをノートにまとめておく。
  - straight chain の決め打ちだったものを一般化する必要がある。
    - ノートのまとめで、基本的な構造はクリアーになってきた。原理的な困難はなさそう、ということ。
    - 文献を調べてみると、しかし closed link な時に効率が悪くなるとかで、いろいろと亜種が開発されているみたい（全部 MD の話）。
  - 6/24/2008: 「高速な implicit algorithm の実装」でも、最初は straight chain 決め打ち。

6/10 (Tue)

[phys]nanoプロジェクト
- dWLC の simulations:
  - CF-cylinder:
    - 延長してた計算が終了。解析してみると、さらに一段、緩和している模様。つまり時間が全然足りていなかったのか。
    - 別の計算を追加（d を振る、 N を振る）。
    - 6/16/2008: N=157, 50000 steps が１００％で回して６日。
    - 6/17/2008: N=200, 50000 steps が１００％で回して７日。
[phys]ryuon
- bead-rod の実装：
  - tri-diagonal の解法
    - wikipedia にあった Thomas algorithm というのを実装しておく。
    - ちなみに numerical recipes には tridag というルーチンが、 GSL には tridiag というルーチンがある。上の Thomas algorithm との関係は調べてない。
    - 6/12/2008: 一般化により、結局、不採用。
  - Liu (1989) を実装しようと結局丸一日費やすが終了せず、家に持ち帰る。
    - うまく動かず、ああでもない、こおでもない、とテストコードを書きまくる。最終的に動く所まできたが、オチはメモリ確保でのサイズ間違い……。まあこういうものだ。
    - これで Lagrange multiplier を求めるところまで完成。非線形方程式の解法には（今のところ） Liu (1989) のものと NITSOL。
  - さて本題の stokes3 にどう組み込むか、それが悩みどころだ。

6/9 (Mon)

[phys]ryuon
- bead-rod の実装：
  - 実際に手を（コードを）動かして見ないと理解できない（頭な）ので、まず実装してみよう。
  - 6/13/2008: no-HI は出来た。
家に着いたら、隣の隣の家が重機によって壊されていた。
- 先日、家の前庭を測量しているあやしいおっさんたちが居て、何してるのか聞いたら、彼らは市から派遣された人たちで、その建て替えためにやってるのだ、と（スペースが充分にあるか、とかいう類のことなのか、と勝手に解釈した）。
- こっちの家全般にそうなのだが、いわゆるツーバイフォーってやつなのかな、規格の角材を組み合わせだけの構造で（私は専門家ではないので知らないが、家の（元）専門家の人が「壁構造」とかって言っていたような気がする）、まあ紙でつくった模型みたいに、あっけなく解体された。

6/8 (Sun)

バースデーを祝って、外食に出かける。
- いろいろと悩んだ結果、インド・バフェに行く（悩んだのはかおだけど）。店に着いたのが１時半、そのとき広い店に客は２組。彼らもわたしたちが食べ始めたころには店を出て行ったので、ほとんど貸切状態。だからといって、バフェの料理がなくなってたりするようなことはなく、また、それぞれ、それなりにおいしかった。貧乏性なので、例によって例のごとく、食べ過ぎーー。まったく、学習しない人間×２である。
- 帰りにペットショップによって、モンちゃん用のトリートを買って帰ってあげる。
- 6/17/2008: インド・バフェ、パート２。
[phys]ryuon
- bead-rod の勉強（先日のつづき）。ちょっと腰を据えて、大事そうな（古典的な）論文を読む。
  - Liu (1989) J.Chem.Phys.
  - Petera, Muthukumar (1999) J.Chem.Phys.
- 6/13/2008: no-HI 系、実装完了。
- 6/14/2008: HI 系、実装開始。
携帯電話の料金を top-up しておく。
- つまりはあれから既に３ヵ月たったということだ。光陰、矢のごとし。
- 今回はだいだい２０ドルくらい残っていた計算（ということは、この３ヵ月は支払いの２５ドル分は使ったということか）。ちなみに、 expire 前に補充すれば balance は keep されることは確認した。（以前 expire させたのは、もったいなかったな、ということ。）
- 9/8/2008: おバカなことに、また expire させてしまった……。
tax return のお金が入っていた。
- 4/24/2008に仕上げた件、そのときに自分で計算（して請求）した額よりも１割弱も多いぞ、といううれしい驚き。きっと CRA の人たちがセーミツに計算してくれた結果だろう。２年前にはそのセーミツな計算によって払いが足りんと言われたが、よかったよかった。
[blog] 「働いても働いても 2008年06月05日（木曜日）」@mori(lc)
- 国語の時間。今日は自分の頭が自覚以上に悪いことを認めて、一読して全然分からなかった森さんの書いたことを、自分に分かるように噛み砕いてみる。
  - まず一番のポイントは、「働く」ことと「好きなことをする」ことは別なことだという点。なので「私は一所懸命、自分の信じることを全力でやっているのですが、ぜんぜん豊かになりません」と言い換えれば、あのように（森さんに）文句は言われない。（そのかわりに、まるでサハラ砂漠の真ん中から明後日の方向に必死に走っていく人を見るように、「大変だね」とか「それはかわいそうに」と（冷笑つきで）言われるか、もしかしたら憐れんでもらえるかもしれない。）
    - 6/13/2008: あさっての方向に走っていないか確認。
  - 「仕事」とは、そもそも需要があって成立するわけで、他者との関係以外には存在しない。だから自分の意思とか嗜好とか趣味とは関係ないし、自分を客観的に見ないとコントロールできるものではない。もし客観視しない人間がうまくいっているケースがあったとしたら、それは偶然であり、宝くじに当たったようなものだ。（まねをして出来るものではないし、その意味では天才的とも言えるかな。）
  - 自分を客観視できない人間は、問題の原因を自分と他人の関係（の不一致）に求めることができず、ある人は「自分の頑張りが足りない」と思い、より力を込めて明後日の方向に走りつづけ、ある人は原因を単に他者に押し付けることになりがち。
    - 付記：このへんはミヤダイだな。 cf. 6/14/2003
  - ある人は、仕事とやりたいことが一致していたり、似ている場合があるかもしれない。しかし、それが必ずしもハッピーな状況であるとは限らないし、むしろ逆に不幸な状況も十分にありうる。
  - 別に自己満足が悪いと（森さんは）言っているわけではない。彼は（たしか）書くことが好きなわけではないと言っていたし、毎日々々更新されている日記だって「仕事」として（お金をもらって）書いていると言っている。
- なるほど。つまり有り体にいうと、「お金が欲しければ働けばいいじゃん」ということか。
  - 別にパートだろうがバイトだろうが派遣社員だろうがポスドクだろうが、仕事は探せばあるものだ、と。
- 過去の森語録：
  - 8/27/2007: 仕事とはすべて「労働」である
  - 4/19/2008: 憧れで仕事をする人間に欠けているもの
  - サマリー：職業について
- ちなみに私は「働いても働いてもお金が貯まりません」というような愚痴は言ってないと思う。
  - まあ、「頑張ってはいるつもりなのですが、なかなか落ち着けません」とは言ってるかな。「もう引越しは嫌だ」とか「一息つきたい」とか。終わるときまで落ち着くことがないのが人生なのかもしれないが。不可能な目標設定をしては不幸になるな。
- こうしてみると私は、確かに自分を客観的に見ることは、比較的、弱いかもしれないな。
  - まあ、そういう自覚は既にあって、だからこのように日記を書いているのだけれど。
- 6/12/2008: 大衆化で失われるもの。働くこと（by 森）は大衆化なのかな？
- 6/13/2008: 大衆芸術家の道、長瀬弘樹という作曲家のこと（kenn）。
- 9/10/2008: 某業界の場合は、本当に売り手市場なんだって。
[blog] 遊学「ルネ・デカルト」
- 引用：
  こうして私の愛読書は、一巡してふたたび遺稿『精神指導の規則』にまいもどってきた。ここにはデカルトの生涯にわたる「精神と事物との格闘の軌跡」が読める。この格闘というのは、自然物や事物に犯されそうになった風前の精神が、いっそ犯されるままのほうがどれほど気楽であろうかとあきらめた矢先、相手の一点を凝視したとたんにこれをはねかえしていくという、その格闘のドラマである。これはほとんど「数学の倫理主義」ともいうべき緊張力である。
  
  規則第八に曰く、「探求すべき事物の系列に於いて、我等の悟性の充分に直感しえぬ何ものかがあらわれたならば、そこに停まるべきである」。
  
  まったくその通りだ！デカルトは「街をうろつく前に街をうろつく方法」を知っていたのではなく、「街をうろつけなくて街をうろつく方法をつくった哲学者」であったのだ。
- これまで出会ってきた人の中に、自分で作ったルールにバカみたいに従っている人間が、少数だが複数いた。その人たちはどうもその自分の作ったルール（に従っている自分）に、ある種、陶酔しているようだった。またある人物は、私からは同じ種類に見える別のそういう人物のことを小バカにしているように見えた。その人たちが何でそういうルールを自分に課すに至ったか、という（青臭く、ちょっと恥ずかしいと想像される）話は、その人たちから聞いたことはなかった（し、こっちも聞いたりしなかった）。そういう人は、ちょっと苦手なのだけど、それは何でだろうな？
- ふと不安になって、自分にはこんな子供じみたところはないかな、と思い巡らせてみたりする。だけど多分、ないと思う。元々、そういうタイプの人間じゃないからな…… むしろ、行き当たりばったり的だから、という意味。全然 secure な人生ではない。
- 何で唐突に「デカルト」なんですかって、奥さん？それはね、禅僧に限らず人間とは日々座るものなんですよ。

6/7 (Sat)

[FreeBSD][coding] ATLAS 3.8.2 is just out.
- I found it just after sending the patch for 3.8.1 to the port maintainer... what a timing.
- Release Notes say there are several bug fixes, so let's try again.
- It looks OK for 3.8.2, so I made a patch (from 3.8.0 to 3.8.2) and sent it again.
- ところで（と急に日本語）、件の maintainer の方（は、たまたま日本のひとなのだけど）、実は先日までここ、つまりカナダの、アルバータ州の、エドモントンに、これまた偶然に来ていて、 conference で私の同僚と出会って喋り、その同僚に職場にまで連れてこられた（？）のだそうだ（という返信をいただいた）。しかし世界は狭すぎるね。いやいや、まったく、悪いことはできないもんだな（って、何も悪いことなどしてないですけれど、もちろん）。
- 6/12/2008: committed.

6/6 (Fri)

[phys]ryuon
- code の手入れ：
  - stokes/src/stokes3.c の BD-nl-solver の parsing に間違いがあった。
  - atlas が upgrade できたのを受けて、やっと octave も upgrade できたので、 stokes/octave も通るようにしておく。
  - linux (bue9) で stokes 系の update. しかし swig が軒並みダメ（python 以外は全滅）。何で？
大台が、とうとう、やってきた。
- まあ、何の努力をしないでも、自然に、自動的に、やってくるものである。
- さして自慢できるものもないが、「回り道」に関してはかなり独走状態で独壇場ではなかろうか、とは昨日も思ったばかりだな。
  - あ、「回り道」はきちんと戻ってきてはじめて「回り道」になるんであって、そうでなければ「迷い道」かな。まあいいけど。
- 6/6/2006: これまでの誕生日
- 6/8/2008: バースデー・ランチ。
- 6/29/2008: プレゼント。
- 7/17/2008: 野茂、引退。

6/5 (Thu)

[phys]nanoプロジェクト
- dWLC の simulations:
  - CF-cylinder:
    - 結果の回収。
    - 解析。 transient を除き、初期配置より random number の方が重要なんだな。
    - こういう当たり前なことばかりで、自分のアホさを目の当たりにする日々だが、日々賢くなっているのだと思うことにしよう。
[phys]ryuon
- nonlinear solver の組み込み、続き。
  - NITSOL の C-wrapper (nitsol_c.c) を書く。
  - bd-imp で GSL / NITSOL 両方を呼べるようにコードの統一化。
  - チェック：
    - test/check-bd-imp-nitsol.c をしかるべく update. 結果は確認した。
  - performance:
    - GSL vs NITSOL で速度の比較。非周期の N=20 の場合で４倍から７倍くらいでている。 N が大きいともっと開くかな？
    - 今回の構造変更による performance up は、しかし見られなかった。残念。
- bead-rod model の勉強：
  - 先日の結果から、（今更ながら）その重要性に気づいた、ということ。
  - 文献を（再）収集。いろいろと妄想中。しかしまずは既存の方法を実装してから、かなあ……
  - 『世間では「急がば回れ」と言いますがねえ奥さん、最近ちょっと回り過ぎではないですか？』などと自問してしまう４０目前の今日この頃。（cf. 藤谷さんの日記。）
    - 6/13/2008: さて、回って来ました。
    - 7/15/2008: 知識と知恵とか、学習の効率とか、モヤモヤと。
    - 8/29/2008: 小島直記「回り道を選んだ男たち」は読めなかったよ……。
  - 6/8/2008: 古典論文に目を通す。
[blog] 藤谷さんの日記 (lc)
- どこの誰とも知らない横丁の奥さんに呼びかけたくなるような日に読むと、うんうん、そうだそうだ、と頷いてしまうよね：
  人間としても、自分の頭が自覚以上に出来が悪いことを、僕はもっと認める必要があるのかもしれない。
  
  そこを励ましてくれるのが川上弘美氏の言葉で、この小説は「小説ってのはこういうことを書くとばかにされるものかな」と思っていない、そこがいい、という意味のことを書いてくれている。 …… 川上氏はまた、集中しなさいとも書いてくれた。
  
  辻原登氏の選評は大半を費やして『棺桶』を激賞してくれている。過褒です……というのは表向きの謙遜で、ほんとうはうるっと来るくらい有難かった。 …… これを読むと真面目な気持ちになる。
  （みなさんご承知のことだとは思いますが、この「引用」は自分用に切り刻んでいるので、正しい文脈とかは本文に当たってくださいね。）
  - 6/8/2010: 柳田さんの本音。
- そこに唐突に出てきた平野さんも、いつからかは忘れたし、何でまだ残しているのか分からないけど、アンテナに入っているな。ちょっと手入れをしておかないといかんな、アンテナも。
- 川上さんの方は名前は見たことあるなと思ったら、少なくともここでは見ていたようだ。
- 6/8/2008: 国語の勉強。
- 付記：川上さんの激励を定理にしておく。
[FreeBSD][coding] upgrading ATLAS ports to 3.8.1.
- For a quite long time, I've been stuck at ATLAS 3.6.0, because upgrading to 3.8.0 (the current version in ports) failed during the build. it seems the same symptom here -- PR - ports/119845.
  - cf. freebsd-ports ML;
    FreeBSD Port:atlas-3.8.0,1(Dec 30 21:49:14 PST 2007)
    FreeBSD Port:atlas-3.8.0,1(Jan 6 18:05:04 PST 2008)
- I thought I had tried before to compile the newer version by myself and found it failed. but today i tried again with 3.8.1, and somehow it is done! So I decided to upgrade the port.
- It looks that some patches had already been applied at the main code, and the corrections for other patch looks straightforward.
- It is now working (i guess). My system is FreeBSD 6.2-PRERELEASE (how lazy i am...) with intel core duo T2400 (32bit). i don't know anything about other systems.
- Some notes:
  - Patch should be made by
```
diff -ruN superedit.bak superedit
   
```
    I made a patch file with "-rcN", the context output rather than the unified one (following my own old tips memo). I should update it.
- Ref: freebsd.org, freshports.
- 6/6/2008: update octave code in ryuon-stokes.
- 6/7/2008: the patch has been sent to the maintainer, and...
- 6/12/2008: committed.

6/4 (Wed)

[music] 今日の気になる曲。
- 9:30 切れの直前、（なよなよ系）男性ボーカル、特徴的なフレーズのリピート。
- 付記： Iron & Wine の White Tooth Man という曲らしい。（youtube にあったので、確認済み。）
[phys]ryuon
- nonlinear solver の組み込み。
  - ふと思い立って、以前からやっておかなければと思っていたことに再びアタックする。
    - Newton-GMRES に関して以前調べたとき（11/23/2007）に見つけた NITSOL という library を使ってみようという話。
    - NITSOL の reference は以下のとおり：
      M. Pernice and H.F. Walker, SIAM J.Sci.Comput. (1998) 19, pp.302-318:
      "NITSOL: A Newton Iterative Solver for Nonlinear Systems"
      (http://www.siam.org/journals/sisc/19-1/30384.html)
  - FORTRAN のコードを C から使う、という毎度ある事例で手こずるのだが、何とか出来た。
    - 以前もはまって、解決したような気がしていたが、その時の記述などが見当たらないので、また１から調べ直した。
    - ポイントを「FORTRAN のライブラリを C から使う事例の１つとして」にまとめておく。
  - libstokes に登録しておく。
    - まだ top level （例えば stokes3 の level）から使えるようにはなっていない。

[coding]FORTRAN のライブラリを C から使う事例の１つとして。

FORTRAN の common block の C からの使い方：

FORTRAN の common block

      common /nitprint/ iplvl, ipunit

へのアクセスは、構造体

extern struct {
  int iplvl;
  int ipunit;
} nitprint_;

を介して行える。

使い方の例 in C：

  nitprint_.iplvl = 1;
  nitprint_.ipunit = 6;

引き数に関数がある FORTRAN のルーチンの呼び方、および FORTRAN で使われる関数の C での実装：

難しいと思い込んでいたが、普通に書いたら普通に出来た。

FORTRAN の subroutine

      subroutine nitsol(n, x, f, jacv, ftol, stptol, 
     $     input, info, rwork, rpar, ipar, iterm, dinpr, dnorm)

      implicit none  

      integer n, input(10), info(6), ipar(*), iterm 
      double precision x(n), ftol, stptol, rwork(*), rpar(*) 
      double precision dinpr, dnorm
      external f, jacv, dinpr, dnorm

というものがあった場合、つまり void 型の関数 f と jacb と double を返す関数 dinpr と dnorm があって、それぞれ FORTRAN では

      subroutine f(n, xcur, fcur, rpar, ipar, itrmf)
      implicit none
      integer itrmf, n, ipar(*)
      double precision xcur(n), fcur(n), rpar(*)

      subroutine jacv(n,xcur,fcur,ijob,v,z,rpar,ipar,itrmjv)
      implicit none
      integer n, ijob, ipar(*), itrmjv
      double precision xcur(n), fcur(n), v(n), z(n), rpar(*)

      subroutine dinpr( n, x, sx, y, sy )
      integer n, sx, sy
      double precision x, y

      subroutine dnorm( n, x, sx )
      integer n, sx
      double precision x

と定義されている場合を考える。

（例１） C から nitsol() だけを呼びたい場合（f と jacv は FORTRAN のものを使う）

FORTRAN の関数の C での function prototype declarations:

void nitsol_ (int *n, double *x,
	      void (*f)(int *n, double *xcur, double *fcur,
			double *rpar, int *ipar, int *itrmf),
	      void (*jacv)(int *n, double *xcur, double *fcur,
			   int *ijob, double *v, double *z,
			   double *rpar, int *ipar, int *itrmjv),
	      double *ftol, double *stptol,
	      int *input, int *info, double *rwork, double *rpar,
	      int *ipar, int *iterm,
	      double (*dinpr)(int* N, 
			      double* X, int* incX, 
			      double* Y, int* incY),
	      double (*dnorm)(int* N, double* X, int* incX));

void fbratu_ (int *n, double *xcur, double *fcur,
	      double *rpar, int *ipar, int *itrmf);

void jacvbratu_ (int *n, double *xcur, double *fcur,
		 int *ijob, double *v, double *z,
		 double *rpar, int *ipar, int *itrmjv);

// BLAS functions
double
ddot_(int* N, 
      double* X, int* incX, 
      double* Y, int* incY);
double
dnrm2_(int* N, 
       double* X, int* incX);

C からの呼び出し：

  int n = MAXN;
  double *x = (double *)malloc (sizeof (double) * MAXN);
  double *rwork = (double *)malloc (sizeof (double) * LRWORK);
  double *rpar = (double *)malloc (sizeof (double) * LRPAR);
  double ftol = 1.0e-6;
  double stptol = 1.0e-6;
  int input[10];
  int info[6];
  int ipar[2];
  int iterm;

  nitsol_(&n, x, f_, jacv_, &ftol, &stptol, 
          input, info, rwork, rpar, ipar, &iterm, ddot_, dnrm2_);

（例２） C で f を定義したい場合

void f_c (int *n, double *xcur, double *fcur,
	  double *rpar, int *ipar, int *itrmf);

という関数を書いて、以下のように普通に呼び出す。

  nitsol_(&n, x, f_c, jacv_, &ftol, &stptol, 
          input, info, rwork, rpar, ipar, &iterm, ddot_, dnrm2_);

すべて GCC と GNU Fortran compiler の話。
- 他の環境は知らない。といっても、実際上何があるのかな？ Intel の compiler くらいかな、思いつくのは。
references:
- 6/4/2008: [ryuon] NITSOL の組み込み。
- 5/8/2003: 行列に関する(個人的な)混乱。

帰宅したらオランダから手紙が届いていた。
- あちらが送ったという知らせを受けたのが 5/15 だったから、結構かかったことになるな。やっぱり遠い（のか、カナダの郵便屋さんの問題なのか……）。

6/3 (Tue)

[phys]ryuon
- dWLC の simulations は Nano に移行する。
- confinement の debug:
  - CF-dumbbell:
    - まず過去の結果を確認。
    - いろいろとパラメータを振ってみるが、それぞれ停止してしまうな。
    - 対策をいろいろと思案中。
  - 10/2/2008: 再始動。
- RYUON-manual:
  - Nano から、 imposed flow parameters の指定方法について、その度に手で計算して確認していたことを、具体的に記述しておく。（System Parametersの部分。）
[phys]nanoプロジェクト
- dWLC の simulations:
  - 応用計算を RYUON からこっちに移行する。
  - CF-cylinder:
    - 既存の結果を整理
    - 足りない部分を（再）計算させるための準備。新 stokes3 （正しく continue ができる）でやり直したいということもある。
  - HI-Shear:
    - 解析。統計平均をとらなければいけないな。
    - 既存の論文を見て extensional flow もやることにする。
    - 毎度のこと、どうパラメータを与えたらいいか迷うので、今度こそきちんと documentation しておく(RYUON-manual)。
    - shear をベースに planar を一通り（20個）作り bue9 に回して計算開始。が、しばらくすると悲鳴を上げている。
    - 見てみると bond が伸びすぎ。よく考えてみたら、 dWLC は単なる Hookean なので、当たり前だな……。
    - モデルにまた手を入れなければならないか。
    - 6/5/2008: 勉強、開始。
    - 6/12/2008: できたてほやほやの bead-rod でテスト、予定通り（ダメだった）。
久しぶりの第一火曜日の買い物。
- やっぱり 10% はでかいよな、と。
- 見返してみると、去年の１１月以来かな。

6/2 (Mon)

朝オフィスにやってくると、各人の cubicle にネーム・プレートが付いていた。
- ここの私たちの cubicle は、実験グループなどが labo に使うようような場所に（理論グループなので）机を並べている。かつここの建物の security の堅さから、この部屋に入るにはカード・キーによる承認が必要で、つまり、ここへは部外者は入ってこれないことになっている。そこにネーム・プレートをつけてもねえ（部屋の扉に付けるならまだしも）、と隣人と苦笑する。他に金の使い方はないものか。
[phys]ryuon
- BD の continuation:
  - stokes-nc への実装、つづき。出来た。
  - チェック： CF-cylinder で 100 step を回した結果と 50 回で中断し、更に 50 step 回した物を比較。問題無しってことで、よしとする。
  - 10/6/2008: stnc の設計思想。
[phys]月刊レポート書き。
- 先日のポスターの内容をベースにしたノートをベースに、省略したり落とされた結果も含めて書き上げる。
- back-numbers:
  - 7/2/2008: ７月号（およびこれ以降）
  - 6/2/2008: ６月号
  - 5/1/2008: ５月号
  - 4/1/2008: ４月号
  - 3/3/2008: ３月号
  - 2/1/2008: ２月号
  - 1/7/2008: １月号
  - これ以前は１２月号のリンクから。

6/1 (Sun)

６月になってしまった。
- いや、気が重い。何しろ大台だからなあ……。（もちろん、半分は冗談。）
[phys]ryuon
- BD の continuation:
  - （ことわざ集で有名な松本さん（cf.6/12/2008）の） MT-FAQを見ながら random number generator の状態をキープする（やっと）。実装完了。
  - stokes-nc に interface を実装中。
今日の夕食は知り合いの家に呼ばれる。
- 基本的にホストのカップルともう一組のカップルと私たち。プラス、ホストの元（および現）テナントの女の子たち。カップルはかおの知り合いで、片方が日本人な国際カップル×２。あと、ここカナダで現役バリバリのトライアスリートの人と、モンよりも重そうな、毛沢東という名の老オス猫が居た。
- 最近は、気がつくと一番年長だったりすることが多かったりするのだけど、今回は年齢レベルでも、あと風来坊レベルでも、我々と同等くらいだったこともあり、社交下手な私でも少しのんびりできた、居心地のよい集まりだったかな。

2008年6月