ぎょーむ日誌 2005-01-25

2005 年 01 月 25 日 (火)

• 0820 起床． 朝飯． コーヒー． 今朝はなぜかしら三次元光分布問題にトラップされてしまう． 朝からアタマの中で三次元ヴェクトルがぐりぐりする． うう． 0940 自宅発． 曇． 0950 研究室着．
• まだ提出命令でてないんで， ぱいぷ樹木原稿なおしのつづき …… しかし， こんなもの読んでると， また今朝がた考えてた三次元光分布計算問題のうずに 巻き込まれていく ……
• で午前をぜんぶ思考迷宮の探索に費してしまって， またしても 新しい計算方法を考えついてしまった …… あえて名づけるなら 「光格子平面蚕食法，周期境界条件対応版」 とでもいうようなもので， 「見おろし法」の variation のひとつ． さして独創的ではない …… というか， ごくふつーのアイデアなわけで． いやはや， 私の場合，これまでのやり口が奇をてらうものばかり， と言えるのかも． たぶん， こっちの方式は誰かがすでに考えだしているだろう …… あるいはすでにふつーに使われているのかもしれない． ともあれ， 周期境界条件への対応が容易でそれなりに高速， メモリもむちゃくちゃに食いつぶすわけではない， と期待できそうな気もする．
• 昨年 3 月ごろには， 遮光物のほうを「回転」させる方法 なんぞを一生懸命かんがえてたわけだが ……

[去年 3 月ごろ考えてたやり口]

これって ずいぶんと手ぬきしてるんだけど， それでもけっこう大変なんだよね． いちいち座標変換して並べ換えてやらんといかんから．

• こういうかんぢではなく， また光線だの視線だのといった概念も とりあえず持ち出さずに， 「光子」を水平な格子平面上にならべておいて， その平面の水平をたもちつつある直線 (まあ，これは光線の入射角なんだが) にそって移動していけばよい． で， その光子平面が次第に「蚕食」されていく， と． どんなもんでしょね．
• ここに書いたら脳内から以上のアイデアをひとまずは消去できる， と思ったんだけど， なかなか消せん． 暴走してんのか? どこを押せば私は再起動するんだ? …… それはともかく， この方式はけっこうメモリを食いつぶすような気がしてきた． 三次元格子モデル になっちまうのか? これは 4-5 年前まで使っていて， いまは放棄したやりくちだ …… 調べてみたら kubolib20001014 が最後か． それの限界は このあたり のぎょーむ日誌に． ぱいぷ樹木には使えるかもしれんけどねえ．
• ぱいぷ樹木原稿なおしのつづき．
• 昼飯． 昼飯後に毎週の講座セミナーならぬ宣教師殿の 宣教せみなりお …… まぁ， とりあえず 3 週間後の修論発表会で， だれが不信心ものなのかすべて明らかになるんでしょう．
• ぱいぷ樹木原稿なおしがひと区切りついたんで， 印刷出力して共著者・甲山さんにわたす． Springer の batch 処理は遅々としてるよーで， 「はよだせ」督促はまだとのこと．
• で， インドネシアみやげをいただく． 下の図に示されてるよーな， 熱帯林 (カリマンタン島) 樹木の成長データ． うーむ， なんか面白いというかみょーなパターンが ……

• まあ， これは後日にやることにして， 1850 研究室発． 1900 帰宅． 体重 74.4kg． うーむ． 晩飯．
• なンとなく， また三次元光分布問題を考える． 三次元格子モデルは論外． STL の bitset 使っても無理だとわかった …… 1024 × 1024 × 1024 あたりが限界か． これで数百メガバイト使ってしまう (ま， データベイスとかランダムアクセスファイル， って方法もあるけどね)． とゆーことで …… 「光格子平面蚕食法，周期境界条件対応版」 なる方法つかうんなら， すべての遮光おぶぢぇくとは球， どの方向に射影しても円形に「蚕食」される， とするしかありませんなぁ． まあ， 「なにもかも球」 にしたところで近似的になってしまうけれど．
• しかし， 近似でもなんでもいいから， このあたりの計算を 100 倍 ぐらい速くせんことには， ぱいぷ樹木 は改良できん．
• [今日の運動]
  • エアロバイク 50 分間．
  • 腹筋運動 30 × 3 回． 腕立ふせ 5 × 3 回．
• [今日の食卓]
  • 朝 (0850): 米麦 0.7 合． ダイコン・ブナシメジ・ワカメ・コンブ・豆腐の味噌汁．
  • 昼 (1450): 研究室お茶部屋． 米麦 0.7 合． ダイコン・ブナシメジ・ワカメ・コンブ・豆腐の味噌汁．
  • 晩 (1930): 米麦 0.7 合． ダイコン・ブナシメジ・ワカメ・コンブ・豆腐の味噌汁．