KuboLog 2005-08-15

2005 年 08 月 15 日 (月)

0800 起床．朝飯．コーヒー． 0920 自宅発．晴． 0935 研究室着．なかなか暑い．

いきなり関係ない問題を検討しはじめる (じつは登校前からひたすら考えつづけてる)．ちょうど三年ほど前に三次元空間内での「円筒と直線の衝突」とゆー計算問題を検討したわけだが …… これの計算量もうちょい減らせんかと考えてみた．とりあえずの別解としては

二直線間の最小距離 L と最近接点 X を求める (2×2 行列計算)
1. の最小距離 L が円筒直径以下なら 3. にすすむ
円筒両端中心点と点 X の距離がどちらも円筒軸長より短ければ「あたり」

ただし，これはじつは円筒ではなく両端が半球になってる円筒みたいなカタチ，ということになってしまう (3. の計算において)．これを改良するには ……

[円筒と最近接点]

最近接点が円筒内にあるかどうかの判定には少なくとも二段階の判定が必要．過程 3. の改良案として，円筒端中心点 - 点 X - 円筒端中心点のなす角度 (の余弦) を余弦定理で計算し，これが円筒形状から決まる下限より大きい，とする方法もありえるな．

アタマの準備運動がおわったことにして，母子里問題にとりくむ．

とりあえず Canopy クラスにパラメーターファイルとハコデータを読みこませて， 1300 個ばかりの Fbox を生成 (new) させるところまでできた．昼飯．

昼飯後もプログラミングの続き．三次元空間内に生成したハコたちを reference で結合していく．つなぎかたにも三通りあって， top_fbox (直上), effectors (光計算に影響するハコ), neighbors (葉が「移動」しうるハコ) ……

観測してる林分の一番上のハコの top_fbox はどうなるかとゆーと，これは dummy_top_box とゆーのを生成してくっつけている．このあたり欠測点のとりあつかいに注意が必要で ……

あ，ここで dummy_top_box とやらでの光強度を 100% ではなく 99.9% に設定しとくと，局所明るさの尤度評価のときにあれこれ都合よろしいな，という二項分布の性質を悪用した巧妙なごまかしワザを発見する．

お茶部屋で一休み．院生密度いよいよ低い．絶滅危惧なのか? しかし試みに北大生協で「誘引エサ」買ってきてお茶部屋においてみると，稀少野生院生たちが次々ととらっぷされていく．ということでお茶部屋雑談 …… 蛇足ながら，小菅せんせーその他数名は夜行性なので昼間には捕獲できない．

いんちき光計算を実現する世界の構築のつづき．結合した fbox の重みの定義．

[方向と重み]

光量のもっとも単純なモデルを考えた場合，面積あたりのエネルギーは sinφ に比例する (φ は「視線」の仰角)．かかる格子状配列の世界だと

sinφ = 
		sqrt(dz² /
		(dx² + dy² + dz²))

と簡単化される (d* は視点座標と視線上点の座標の差)．

にせデータを生成して，このへんの各種試験運転．三次元空間ばぐがぼろぼろと見つかる．

修正作業してるあいだにいろいろと．統計こんさるとか．近ごろ卒業生たちのあいだで流行してる， ssh & りもーとぶらうざわざ，とか．

見つかったばぐは全部修正．とりあえず，空間構造をあらわす「配管配線工事」，半分ぐらいは片づいたか． 2035 研究室発． 2050 帰宅． 2205 自宅発北大構内走． 2250 帰宅．体重 73.4kg．晩飯．

また一個 Gmail 発送作業．どんどん Gmail 使いが増えてくね．ところで本日，私の web browser から Gmail にうまくアクセスできんかった．ログインすると「読み込み中」とかで止まる．一日中原因がよくわからなかったんだけど …… さっきふと思いついて browser のキャッシュを全部クリア したら，問題なく Gmail にアクセスできるようになった．なかなかナゾ．

北大構内走中に上記「配管配線工事」の意外なめんどくささに気づく．てきとーにやっても問題なさそうに思ってたんだけど，じつはてきとーすぎると境界でヘンなことになったりするかも，という可能性ありそう．たとえば「はしっこ」で葉っぱがすかすかになるといった artifact だ．

これを避けるためには葉っぱ「移動」なんかのときに慎重な加重計算ともなう乱数発生が必要になるか．いやはや …… と書いたところで少しマシな対策を思いつく．移動先の「加重」を数値ではなく reference の個数で表現する，という方式はどうかな? たとえば fbox1 の加重が 1 で fbox2 の加重が 2 なら「行き先候補」配列を (fbox1, fbox2, fbox2) としておく，という方式．境界めんどうも reference つけかえできりぬける，と．

ちょっとそのあたり実装してみる．うまくいきそうだ．

[今日の運動]

北大構内走 45 分間．

[今日の食卓]

朝 (0850): 米麦 0.7 合．タマネギ・キャベツ・ニラ・ブナシメジ・ホタテの炒めもの．
昼 (1240): 研究室お茶部屋．米麦 0.5 合．タマネギ・キャベツ・ニラ・ブナシメジ・ホタテの炒めもの．
晩 (2320): 米麦 0.7 合．ミニトマト．タマネギ・キャベツ・ニラ・ブナシメジ・ホタテの炒めもの．麻婆豆腐．

ぎょーむ日誌 2005-08-15

2005 年 08 月 15 日 (月)