ぎょーむ日誌 2001-03-04

2001 年 03 月 04 日 (日)

• 1120 起床． いかん， 時差が生じて来た．
• 「ニセ」アカマツのでたらめな三次元計算系の修正． 枝の出かたが変になる原因は PolarXyz::Rotate という回転の関数にある， と特定． 第一象現から第四象現までぐるっとまわってほしいのに， 第一象現の中で車のワイパーみたいな 往復運動やってるんですよね．
• あぁ， ヴェクトルがくるくる回ったり 反転したりして， なかなかややこしーぞ． ♪ たーおせー， 火を吐く大怪獣 …… debugging ハイである． しかし， なぜこの Canna 強化辞書は 「大怪獣」を成句として一発変換するのか ?
• 1330 よーやくにして三次元バグ取れる． PolarXyz::Rotate から呼び出される PolarXyz::MultipleScalar でスカラー値が負のときに 仰角φのヘンな補正をやってることが原因であった．

[「復元」されゆく樹木] ちゃんと表示されているようです． 広島大の院生の皆さんが執念深くばらばらに 解体してしまった樹木個体は， いまここに「ニセ」アカマツたる形態を 与えられて復元しつつある．

これは freeware のレイトレイシングソフトウェア POV-Ray を使った作画である．
• 針葉をつけてみる …… うまくいってるかな …… いや， 変だ． 三次元計算， 「だいたい」うまくいってるんだけど， ときどき変になる場合がある， と不幸にも発見してしまった． やはり回転がらみである． 腰をすえて改善にとりくむ．
• うーん， うまくいくはずだと思うんだけどうまくいかない． ちょっと気分転換にパンケイキを焼いてかじる．
• 1720 今度こそホントに直る． 結局， さっき直したと思った部分が 完全には直っていなかったのである．
• ああ， 最後に気になるのはこのクラスというか Perl モジュールは汎用的に使ってる kubolib 中の三次元座標系クラス (C++) から移植されたモノである， という点だ． 分家は修復されたけど， 本家は大丈夫なんだろうか． あっちのほうが かなりテストを経てるはずなんだが ……
• ま， 今回の事件は Perl への移植ミス， ということにしておこう．

[葉っぱつけてみました] こう遠目に見るとそれっぽく見えてるかも …… ズームすると手抜きがばれます． 写真 と比べると葉の量が足りませんね． ついでに 全 8 個体 (55 kb, jpg) とか．

開空度計算などという， わけのわからぬコトをやらなければ， Perl と POV-Ray だけでも これぐらいはできてしまうのである …… とか言うわりには 苦闘してるけど．
• おおっ， ばかばかしいけど愉快な描画手抜き方法発見． 針葉がついている茎だけ 「樹皮の色」 ではなく 「葉っぱ色」 で塗ってやると， (実際にはすかすかであっても) あたかもそこにたくさん針葉があるように見える． なかなかセコくて楽しい偽装工作です．

[欺瞞工作] ここまで拡大するとバレますが ……

• いろいろと描画して遊んで晩飯時刻まですごす．
• 晩飯食ったあとは遊び疲れたのか， でれでれと本など読んですごす．
• 久しぶりに今日の朗読
  ``Analytical Population Dynamics'' (Royama, 1992)
  Part I Theoretical bases of population dynamics
  Chapter 2. Structures and patterns of population processes
  • 2.4 Density-dependent/independent processes
    • 2.4.1 Density-dependent/independent process as stochastic processes
    • 2.4.2 Order of a density-dependent/independent process
• ということで， いままでは個体群密度だけで決まる決定論的な過程ばかりを 論考していたわけだが， ここからはさらに密度には影響されないとして定義される density independent な要素もつけ加わる． これは確率変数の形でモデル内に取り込まれる． で， 見ると …… 短期データばかり取り扱っている私などにはアタマ痛いことに その density independent なやつも k ステップ前まで 遡って影響されるような定式化である． おそるべし時系列モデリング (私はどちらかというと空間こんしゃすもでらー?)． 密度も h ステップ前まで考慮するので， AR(h,k) などと書かれたりする． うう． 今日はちょっと読んだだけで挫折．
• 風呂などにのんびりとつかる． 過去三日分ぐらいの無精髭をそる． 長時間湯につかっていたので， でると少しふらふら． 爪を切る． 私は身だしなみがいいかげんであるが， 爪ばかりは丹念に切っている． ちょっと伸びてるだけでも， 打鍵するとき気になるからである．
• なんか急に仕事が進まなくなった …… どうして「ニセ」アカマツ育てようとすると エラー出るのか …… あっ !! 思い出した． 新しく生成された茎の三次元座標を 与える関数をまだ書いてなかった． ははは．
• いやいや， その関数 RedpineTree::StemLocator はすでに書いてるんだけど， どこからも「お呼びではない」 状態に置かれてたから 動いてくれなかったワケね． いやはや．
• よくわからぬところの多い動態モデルで 生成した樹木の作画でもやってみるか ―― なーんか， えらく時間くうな． そうか針葉の位置計算にてこずってるな． 中止． 一番新しいのは 7 年生か． こいつを作画してみよう． これまたエラく時間かかる． 個体を構成するオブジェクト数 (部品数) 15 万ちょっと． PentiumIII 650 MHz では苦しいか ? POV-Ray はかなり愚直に光計算やるからなぁ． 5 分以上経過したけどまだ終わらん．
• そういえば POV-Ray に「雑に描画しろ」 というコマンドラインオプションあったな． あれはどう指定するんだっけ．
• 12 分たっても終わらないので中止． 「粗くやってくれ」 x-povray +Q0 と指定してみる …… ダメだ， これでもいつまで待っても終わらん．
• とりあえずシミュレイター側に「針葉出力ナシ」 モードをつけるか． いや針葉なんぞ出さないのを default 設定にしよう．
• 針葉を出さなくても苦闘してるな． 茎数が 9 年目で 44686 本か．
• とりあえず 10000 本弱の 7 年生で描画してみると …… 針葉なければ数秒で描いてくれるな …… ぬう， なんたる錯乱状態．

[7 年生「ニセ」アカマツ] 破綻しつつあるな …… これは対照区の個体なんで 成長が良い． しかも光計算してないんで， こみあいによる抑制効果がない．

ああ， これがさらに 2 年も「育つ」と ……

[9 年生「ニセ」アカマツ] うに

まぁ， いいや． 描画機構を作ったのは， こういう破綻なんかを視覚化するためだし．
• 時刻は 2350． 明日の飯の準備でもやって寝るか．
• 今日の食卓
  • 朝 (1140): 米 0.4 合． 昨晩のスープの残りで雑炊． コマツナの残り．
  • 昼 (1530): 朝の残り．
  • 晩 (1900): 100 円ショップで買った蕎麦を 50 円ぶん食べる．