ぎょーむ日誌 2001-05-01

2001 年 05 月 01 日 (火)

• 0650 起床． 寝たのは 0500． 眠い． 結局， 週末時差が解消されなかった． いやはや．
• 朝飯・弁当の準備． 朝飯． コーヒー． シャワー． ひげそり．
• 0720 自宅発． 曇． 気温低い． 0728 京急平和島発． 普通に乗り換えて 0751 京急上大岡発． 0810 研究所着．
• 先日 openssh のパッケイジを入れ換えて以来， どーにも不調な 「パスワードなしで ssh を通す」 設定に苦労する． 朝の 40 分以上費してみたが努力が実らず． $HOME/.shosts や /etc/ssh/*_config あたりはちゃんと設定してるつもりなんだがな ……
• さらに 30 分粘る． ssh-keygen コマンドとか? …… 違うな． うーむ． いったんあきらめる． どうも表ネットと闇ネットが錯綜してるあたりが 問題なような気もする．
• あきらめて Pasoh 自動処理 Makefile 書きに戻る． なぜかこれがやる気がおこらん …… 原始的とはいえインターフェイス作りだから， だろうか．
• 書類書きでスサんでしまったらしい (*BSD と Ruby を布教してる) 長谷川宣教師 から爆弾テロ． Linux は /usr/share/misc が無いからけしからんとか わけのわからぬケチをつけたあげくに， BSD の当該ディレクトリを圧縮して郵送するという …… ココロはタリバン， 手口はゆなぼまーのごとき無法な御仁である．
• 白濁しがちな脳内状態に耐えつつ， 光計算のところまで終わらせる． なんかバカというか 洗練されてないというか 柔軟性の欠如した 自動処理機構を作ってしまったんで， 修正していてもうんざりする． どうにかならんものかね． 午前中はこれで終わり． 昼飯にする．
• 天気も悪くないので， 外に 30 分ほどふらふらと出かける． まず富岡総合公園に直行する裏道探索の続き． 車の多い横須賀街道を避けて南東に進むと， すぐに石川島建機・日本飛行機といった 工場にはばまれる． 住宅地のある尾根にいたる断崖を登るのは やめたほうがよさそう． 技術的に困難というより 登攀中に警察に通報される危険性が高いためである． 残り時間も少なくなったんで独房群の周りを ふらふらする． コンクリートで固めたり蓋をしている 幅 1m ぐらいのドブ川の横の道にまで わざわざ「杉田川散歩道」なんぞという 看板を 50m おきに設置する 神奈川県担当部署の正気を疑う． 1255 独房帰着． 眠い．
• とりあえず make all とやると Pasoh 解析を最初からやりなおすような Makefile を書く． だからなんだというんだ， というかんぢなんだが． 時刻は 1315．
• いちど「通し」でやり直してみるか． おっとその前にフィルターを もうひとつ入れないといけないんだった． えーと ……

  第二回観測日と第三観測日の間が 4.5 × 365 日以上 6.5 × 365 日未満のものだけを選ぶ

  …… か． 何かのまじないなんだろうな．
• Makefile の help など書き直す． 台本はこんなところか． 1410 ごろから 全体を「通し」で実行してみる実験開始． たぶん 40 分ぐらいかかるだろうなぁ．
• 50ha (弱) の開空度計算は例によって 25 分ぐらい． パラメーター推定にやけに時間かかってるうえに 樹種によって収束しないでコケてしまう． ああ， やっぱりパラメーターはできるだけリニアーに つないだほうが良さそう， という教訓が確認された．
• つまりこういうことである． 甲山さんのごようぼーによって 胸高直径 (dbh) 肥大成長式のカタチを

	Δdbh = p0 * dbh * exp( -p1 * dbh - p2 * darkness )

…… てなふうにしてたんだけど， こういう非線形なのは難しいんだよね． パラメーター (p0, p1, p2) の値を決めるのが． p0 がつねに正であることに着目して， 上の式を

	Δdbh = dbh * exp( p0' -p1 * dbh - p2 * darkness )

と変形してやると， すーっと推定値が計算できるのである． うまく理屈で説明できぬケド， これは有効な経験則かもしれぬ ―― 「最尤推定するパラメーターは できるだけ線形にくっつけろ」．
• 471 樹種のパラメーター推定と簡単な作図で だいたい 17 分 20 秒か …… とゆーことで，

  前処理 (約 2 分)

  38 万本の第一段選抜， 50ha を 800 個の 25m プロットに 分割．

  シミュレイション (約 27 分)

  Pasoh シミュレイターなどのコンパイルと 開空度ごとき指標の計算．

  推定 (約 17 分)

  成長モデルのパラメーター推定と 簡単な作図．

  …… ま， ホントに全部やり直したら 45 分かかるか． しかしモデルの変更が成長式の部分だけなら， 最後の 17 分をやり直すだけでよい． 樹冠の形状まで買えてしまうと その前の 27 分の過程を通過しなければならぬ．
• 問題は， この推定されたパラメーターの使いみちなんだが …… 当方にはアイデアがまるで無いんで， とりあえず発注者の甲山さんにおまかせしますか．
• メイルなど書くと時刻は 1600． さーて， つぎはどうしようか． 三次元計算ライブラリー改良の続きか． 最尤推定法まわりの再検討か． それともそろそろ Pasoh 解説の web page でも作るか．
• とりあえずストレッチでもして一休みしよう． 屋上に上がる． 高くて広々としていて 誰もいないので気分転換になる．
• fj.comp.lang.c なんぞを見てると， 前置・後置演算子の ++ を使うときに i ++ とゆーふーに変数と演算子の間に 空白を入れるような書き方は 「問題外」 とばっさり切り捨てられてるのを見つけて戦慄する …… ああ， 私は Perl でも C++ でもその流儀で 書いてしまってる (そういえば C++ とくっつけて書くなぁ)． とりあえず Pasoh 関係のコードは sed で一括置換して証拠隠滅する． くわばらくわばら．
• うーん， Linux カーネルなんかのソースコード見ると， 単項・二項演算子は前後に空白置いても良さそう． 前置・後置だけが 「変数間空白」 御法度なのかな …… どうもそのようだ． あと， 複数行に式を分けるときには 演算子を行末に置かない， とかいう決まりもあるみたいだ． これも知らなかった．
• とりあえず三次元ライブラリー改善作業に取り組む． 二週間ほど放置しておいたので， 以前の作業の痕跡を調べるところから ……
• バテたので撤退． 1800 研究所発． ちょっと古本屋に寄り道． 1824 京急杉田発． 快特に乗り換えて 1829 京急上大岡発． 普通に乗り換えて 1850 京急川崎発． 1900 京急大森町着． 商店街で買物． いつもの 100 円ショップで折りたたみ式 ノコギリが 100 円で売られていたので， なんとなく買ってしまう． ノコギリなんか何に使うんだ． 雪洞でも掘るのか． 第一京浜は国道 16 号ではなく 15 号であることに気づく． いやはや． 1920 帰宅．
• 今月 24 日の中間評価とやらに関するメイルが来ている． 中間評価といってもこちらが審査されるんじゃなくて， 地球フロンティアという やくざな組織が外部評価されるのである． ところが間抜け評価委員の誰かが 下っぱにも報告させろとか言ったせいで， こちらも準備しなければならぬ仕儀に． なに? 密室内で 22 分間もかけて リモセンだか海洋だかの専門家相手に 森林モデルのハナシしろって? いやはや． そういう細かいところまで評価したいわけかね．
• 晩飯食ったあとガス代の振込みに出る． 1535 円．
• クレジットカードの支払請求来る． 今月は本を買ったので出費多い． 東京めたりっく通信は月額 5971 円か． さーて， まだ請求が来てないけど， NTT への支払いはこの金額以上に減じているはず， と期待してるんだが．
• 眠気に耐えられず 2 時間ほど寝てしまった． 起きたら真夜中． 晩飯食った後に眠くなったら 朝まで寝てしまうようにすべきか． いったん起きると なかなか眠れないんだよね．
• 今日の食卓
  • 朝 (0655): 米 0.6 合． 昨晩の味噌汁残り． ニラ・マイタケ炒め．
  • 昼 (1215): 弁当． 米 0.8 合． ニラ・マイタケ炒め．
  • 晩 (1950): 米 0.8 合． チンゲンサイ・タマネギ・マイタケ炒め． 味噌汁の残り．