ぎょーむ日誌 2004-03-09

2004 年 03 月 09 日 (火)

• 0850 起床． ばて． コーヒー． 0920 自宅発． 曇． 0930 研究室着．
• さーて， 明日から関東流浪の旅だし， さっさと林冠射影問題をやらねば， と思っていたのに …… Matt の熱帯林計算こんさるたんと業で午前中がツブれてしまった． 計算状況はひどく呪われきっていたんだけど， 数時間かけてどうにか脱出への糸口が解明された． 原因は甲山学派などでよく用いられる「上側 BA」 なる奇妙な指標の奇妙な性質によるもので …… これはふつーの確率分布で表現するのは 原理的に不可能 とわかった． 確率分布をあてはめたければ切断分布など 条件つき 確率分布を適用せんといかん． なんでいままでコレに気づかなかったんだろうか． 私が上側 BA とか使わんから， かな．
• 雪ふってるんで北大構内走やめ． そもそも時間がない． 昼飯．
• 熱帯林からまだ開放されづうだうだと．
• 1400 すぎによーやく林冠射影問題に復帰． まずは天球上の光源の強度の設定だが …… 天球分割三角形の面積で規格化してみるか． いまだに釈然としないところもあるけど． まあ， 森林全体を「十分に」多いつくす半球に放射される エネルギー量を制御してる (近似的に)， という状況のつもりなんだし．
• つぎに球状葉をヤメて円盤状葉にしてみる． どうも葉っぱが平面的でないといまいち面白くない． ということで， プログラム変更してみるんだが， 急いでるときにかぎってバグがなかなかとれん． うう．

[円盤状葉群]

林冠がホントに傾いてるわけではなく， こういう射影林冠によって 「ナナメ上に光源があるときの光減衰」 の様子を計算してる状況を図示したもの． 葉群にぶちあたると光柱が削りとられる， という電算機図による比喩表現．

• えーい， 計算結果がヘンだと思ってたら， 先週の 木曜日 にやった計算が間違ってた． 林冠回転において， 「林冠を持ち上げる」 角度は正しくはこうだ (光源が天頂にあれば回転しない; 地平線方向にあれば 90°回転):
  • X 軸まわりの回転方向を決定する; 時計まわりに φ すなわち
    • cosφ = v • (0, 0, 1) = vz
    • sinφ = (1 - cos2φ)0.5 = (vx2 + vy2)0.5 = 1 / c
• さて， 円盤状葉群を採用した理由は， 下のようなふたつのタイプの林冠で差が出てほしいなぁ， という希望によるもの． 左側の林冠では円盤葉が全部「真上」をむいてるのに対して， 右側ではでたらめな方向をむいてる， という状況． よーし， 左側をコヒーレント林冠 (coherent canopy)， 右側をランダム林冠 (random canopy) とでも命名してみよう． どちらも林冠の極限形態になっている．

• 計算に要する時間は， と …… 葉群数 1000， 回廊数 256， 光源数 256 という設定で dellpn (Pentium4 2.4GHz) に計算させると 61 秒． Perl 数値計算だからこんなもんかねぇ． 呪われ言語 C++ に置き換えたところで， せいぜい数倍速くなるぐらいだろう．
• ちなみに明日からの出張で持ち運ぶ ThinkPad 240Z (Celeron 0.5GHz) だと同じ条件のもとで 250 秒． まあ， クロック数をそのまま反映してるよーで． せめて ThinkPad X31 (PentiumM 1.6GHz) が壊れてなかったら．
• 納得のいく計算結果は得られつつあるが， いまいち面白くないのう． 2110 研究室発． 2130 帰宅． 体重 73.6kg． 晩飯．
• 晩飯後もコヒーレント & ランダム林冠をヒネってるうちに， 面白い状況が存在することに気づいた． 天頂光源に全エネルギーを集中する 「門司-佐伯極限」， これは「極限」の名にふさわしく (ってまた自分で勝手に命名してみたんだけど)， ひじょーに偏った分布を生成するとわかった …… よしよし， これを使えば明後日の横浜独房群での Dennis さんへの説明がかなり楽ちんになるぞ．
• げ， もう真夜中すぎてるぢゃん． まあ， 明朝はそんなに早い出発ではないけれど．
• [今日の運動]
  • 腹筋運動 30 × 3 回． 腕立ふせ 3 × 3 回．
• [今日の食卓]
  • 朝 : 食ってない．
  • 昼 (1300): 弁当． 研究室お茶部屋． 米麦 0.7 合． キャベツ・モヤシ・ブナシメジ・煮干の味噌汁．
  • 晩 (2220): 米麦 0.9 合． キャベツ・モヤシ・ブナシメジ・煮干の味噌汁．