ぎょーむ日誌 2001-08-(11-20)

2001 年 08 月 11 日 (土)

• 1100 起床． うー， 乱れた生活．
• あ， 雨が降ってる． 東京都の秘密兵器 ヨウ化銀アセトン溶液「雨乞い」装置 のご利益かしらん? いずれにせよ気化熱で地熱が多少は下がってくれるだろう．
• 夜まで …… あまり建設的ではない計算機設定パズルで遊んでしまった． えーと， 昨日書いてた Galeon というのを自宅の ADSL 闇ルーター兼外づけ計算機 Mate で走らせようとしたんだけど， この Mate に入ってる Mozilla のヴァージョン (0.9.3) と VineSeedPlus に入ってる Galeon のヴァージョンが 対応してなかったのである．
• Galeon は Mozilla の高速描画エンジン gecko に 寄生して Web 画面を表示するんで， ここらはきちんと対応させなくては動いてくれない． ということで， ひさびさにrpm --rebuild を連発する作業になった．
• Galeon のような Gnome なソフトウェアを コンパイルするのは存外面倒である． というのも， GUI その他が Gnome というシステム環境に依存しており， しかもその Gnome が無数の部品に分割されていて， それぞれ改良されているからである． たとえば galeon-0.12pre1-1 の場合， こういうヴァージョンの部品が必要とされている．

	GConf >= 1.0.3
	ORBit >= 0.4.0
	gdk-pixbuf >= 0.10.0
	gnome-libs >= 1.2.0
	gnome-vfs >= 0.6.0
	libglade >= 0.13
	libxml >= 1.8.11
	oaf >= 0.6.2

ま， いちいち ftp でばらばらに rpm package とってたころよりは， 今はラクにはなっているんだが． 言うまでもなく apt-get install と apt-cache search が使えるためである．
• 上で一番面倒だったのは Gconf-1.0.3． というのも， まだ VineSeed にも入っていなかったためである． しょうがないんで， GNOME ご本尊 までお参りにいってみたんだけど srpm package なかった． しょうがないんでネット上をふらふらとさまよって 1.0.3 の srpm package をとってくる． で， こいつを rpm --rebuild するためにまた別のpackage が必要で …… という寄り道の連鎖を楽しんでいたわけだ．
• あたらしい Galeon もちゃきちゃき動くようになったんで， (Galeon はひじょーに快適 …… ってんでまたあちこちのサイトに遊びにでかけてしまったんだが) …… さーて， 仕事でもしようか …… と思ったのに， ふと魔がさして latex2html という変換ツールも新しくしたところ …… なんと論文書きに必要不可欠な pLaTeX 一式が 自動削除されちまって代わりに tetex とかいうものが インストールされてしまった． なんぢゃ， これは?
• えーと， この teTeX とゆーのは Tohmas Esser 氏なる人物が保守してる Unix 用 TeX package か． ハノーヴァー大学のヒトみたいね． で， Vine 用の tetex はアスキーから配布されてる 日本語 pTeX よりちと新しい， ということのようだ (なぜならアスキー pTeX も teTeX を日本語化してるから)．
• ということで， この新しい tetex が動くようにいろいろとインストールした 後になって sudo apt-get install task-tetex とすればまとめて何もかもインストールしてくれるようになってる， と気づく． どうも Vine ではこの sudo apt-get install task-X で X 関連をまとめていれてくれるみたいね． apt-cache search task にて調べてみるべし．
• 昨日てきとーに書いたところを直す． たとえば， 昨日のぎょーむ日誌に紹介してた 「親より小さいモノのうちから最大のモノを子供として選ぶ」 なる操作は， 結局こう書き直さなければならなかった．

いままでこんな表現方法を「発明」しなければならなかった 生物学の論文があっただろうか? 「独自性」以外は売り文句のない 「ニセ」アカマツを開発せしめた狂気の一端が うかがわれる一節である …… え， 当方の敵性言語能力がヘボいだけ? うう
• まあ， そもそも上のような書き方だって 「発明」したんじゃなくて借用してきたにすぎんのですが． いやいや， 借用というと貧乏ったらしいので， 極東ものまね国のでたらめモデル屋らしく 「伝統にのっとった表記法を採用している」 とでも書いてやろうかな．
• あ， いつまでも「ニセ」アカマツやっていてはいかん ……
• 今日の食卓
  • 朝 (1130): そうめん．
  • 昼 (1730): 蕎麦．
  • 晩 (2130): なんか食欲わかんのでパン． ゴーヤ・タマネギ・エリンギのオムレツ． 意外とうまくできた．

2001 年 08 月 12 日 (日)

• 0900 起床． 日曜日なんだけど目覚し時計で無理矢理起きる． 生活周期をずらすための努力である． 眠い． 4 時間ぐらいは寝たんだけど．
• お， また雨が降ってる． よしよし．
• StarWars の EpisodeII は ``Attack of the Clones'' ということで， いよいよクローン戦争に突入． 人体複製はこの世界の重要技術らしくて， Episode IV/V/VI のやられ役だった ストームトルーパー たちはことごとく クローンであったとか (連中のオリジナルは誰なのかは そのうち明らかにされるらしい)． クローンなのでいくらやっつけてもよい， というひどく難解な正当化が あの勧善懲悪な活劇を非難されないために 必要とされてるみたいね．
• クローンといえば， 先週発表されたはずの 敵国棟梁の Bush の演説はどういうものだったのかな …… ES 細胞研究への補助を一部容認 か． まぁ， 双方に妥協してみせたわけだ．
• …… なぜかすでに真夜中である． 「ニセ」アカマツもオオシラビソ PipeTree も 微動だにしなかった． どーもこういう二正面作戦は難しいね． もっぱら計算環境とプログラミング言語のことばかり調べていた． わかりやすく言えば， 逃避と呼ばれる行為である．
• 「ニセ」アカマツは最大の難所のモデル説明に関して， ほぼケリがついたので明日からはオオシラビソ中心でいくか． 気持ちとしてはさっさと「ニセ」アカマツ終らせたいんだが．
• 今日の食卓
  • 朝 (0930): 昨晩の残り． パン． ゴーヤ・タマネギ・エリンギのオムレツ． コーヒー．
  • 昼 (1330): そうめん．
  • 晩 (2100): 米 0.7 合． 栄養あるものでも食うかと鶏レバ買って来る． ピーマン・ニラ・ニンニク茎・エノキダケと 炒めて食う．

2001 年 08 月 13 日 (月)

• 0700 起床． 寝たのは 0430 すぎ． もっと早くに寝ようとしたんだけど 「PipeTree 改造したらまたメモリを多量に浪費するんじゃ なかろーか」 といったつまらん問題が脳裏に去来して， なかなか眠れなかったのである． いつものことながら， 眠い． くそう．
• 朝飯・弁当の準備． 朝飯． コーヒー． シャワー．
• 0750 自宅発． 晴れ． 特急に乗って 0801 京急平和島発． 普通に乗り換えて 0828 京急上大岡発． 0845 研究所着．
• アタマの準備運動と称して， BookPC では不調ぎみの Galeon の設定を再点検してみる． 試行錯誤で調べてみると， どうやら galeon-0.11.5-1 では ウリである tab ブラウザーとしての設定まわりが ヘンらしい， とわかった． なんか設定すると起動できなくなる． tab つきペイジを生成することはできるけど．
• ぼちぼちと PipeTree に取り組んでみる． 「ニセ」アカマツはまたあとから．
• どうでもいいが毎日のように SirCam なる virus に汚染されたメイルが 毎日のよーに地球ふろんちあに来冦してくる． たぶん関係者の中に M$-OutLook 使いがいて こうなってるんだけど， どうもその莫迦の Private アドレスから 送られてるらしくて当人は 200 名以上に迷惑をかけてるのに気づいていないようだ． そんなメイルを全員に転送するここのメイルサーヴァーも いいかげん腐りきった設定だと思うが． つくづく間抜けな管理者だ．
• 近ごろ BookPC が不調ぎみ …… 冷静に考えると， どーも一月ほど前に 64MB の SDRAM (PC100) を 256MB に取り換えてから動作が不安定になってるような気がする …… 256MB を外すのは悔しいので， とりあえずもうひとつのスロットにささってる 128MB を取り除く． 引き出し内で休眠状態にあった 64MB を装着して様子を見る． 合計 320MB．
• さて， とりあえず去年 11 月ぐらいから放置している PipeTree を …… と思ったところで甲山さんから 計算ぎょーむ発注メイル． Pasoh ぢゃない別の熱帯林 (Serimbu)． うう． 「帰無仮説棄却的」というべき AIC の計算． めんどくさい． ああっ， またまた ぱらめーたー比較? 毎度のことながらこの方面の群集生態学は奇怪きわめる．
• この程度の disturbance でぱにっくに陥っていては 生態学計算下請けぎょー者はつとまらん． おもむろに archive から Serimbu 関係のファイルを取り出す． 「こんなこともあろうかと」 Makefile で全系統の動きを自動化しておいたから， さっさと片づきそうだ．
• make all で以前の計算をゼロから完全に再現してみる． 各プログラムの挙動を思い出す． 改造のポイントを整理．
• とりあえず 27 樹種が個別のパラメーターをもつか， 共通のパラメーターをもつかを調べる最尤推定を 処理の流れに組み込んでみる．
• まだ作業途中だけど， 昼飯の弁当を食う．
• 総理大臣が 15 日ではなく今日靖国神社に参拝するとか． さーて， ネットを経由した東京大空襲はまた発動するのか? 単なるデマだった， という可能性もある． しかしもしホントに攻撃準備してるやつらがいたら， 振りあげた拳は叩きつけるだろうなあ． おもしろ半分で． そして何度やられても抜本的な対策をとれないのが， 本邦お役所のめでたさである． 今も昔も．
• 昼過ぎても Serimbu 関係の作業． 1520 計算ぎょーむ終了． 全樹種共通値か個別値かを判別する尤度計算と， 趣旨のよくわからぬ二重フィルタリングした 推定結果を出す． うーん， この計算結果たちは どういうふうに扱われるんだろう …… とはいえ， どうせここから先はよくわからぬ世界なので 発注元でやってもらうことにする． あれこれ考えるのはヤメて， 必要なファイルをそろえて送る．
• 自分の仕事に復帰することにする． その前にちょいと一休み． 考えてみればオオシラビソ PipeTree というのも 甲山さんの下請け & 代理発表だな．
• まずは PipeTree 計算環境を更新する． C++ なプロジェクトとしては一番新しい Pasoh シミュレイターを参考にしながら Makefile などを修正していく．
• 去年 11 月の時点では三次元虚業ライブラリと 三次元実存 Vit の組み合わせはまだ 影もカタチもなかったんで， そのあたりも Pasoh から PipeTree (2001 年 9 月版) にも移動する．
• BookPC 安定して動いてやがるな． よくわからぬ理由で Buffalo (128MB) と プリンストンテクノロジー (256MB) の 組み合せに問題が …… いやいや， あと何日か様子みてみよう
• 1800 研究所発． 普通に乗って 1810 京急杉田発． 特急に乗り換えて 1820 京急上大岡発． 1850 帰宅．
• 三谷幸喜という脚本家 (有名人らしいんだけど知らなかった) の 「気まずい二人」 (ああ，なんて題名) という女優・女性歌手たちとの対談集における 三谷氏の挙動が私と似ている， と指摘を某所よりうける． なんだか気になったので， 独房群からの帰路， 京急平和島の駅の本屋で買う． 角川文庫． 税込 560 円．
• 「途切れる会話，最悪の出だし， 重苦しい沈黙，激しい動揺． 緊張の嵐，もやしと枝豆の話， 呆然とするゲスト，焦る三谷」 と身の引き締まるあおり文句の 本書を読んでみて …… うーん， それほど似てないと思う． だって， 話が続かなくなったり， よそ見したりとか誰にでもありますよねぇ．
• しかしサスペンスな内容でした． 「こう言われたらどう対話を続けよう?」 と自問自答しようとして， 頭の中真っ白に． 研究のハナシだめ， 計算機バナシだめ …… そんな極限状態で女性相手に何を話せばよいのか? 西田ひかるとの恐ろしいやりとり．

  三谷「結局， 何の話を …… あの，本当に怒ってないですか」
  西田「怒ってないですよ」
        にっこりと頬笑む西田

  私はかなりおびえました．
• 十数人の女性相手に対談を続けて行くうちに 著者のしゃべり技術は次第に向上していく． しかしあとがきによると， その後また人とあわない生活が続くと 元の「人見知り」な状態に逆戻りしてしまった， とも …… ああ， 独房群で毎日沈黙を守ってるとホントに 口頭のやりとり技術がダメになるな， とは私も実感してる． もともと得意でもなかったしね ……
• 今日の食卓
  • 朝 (0710): 米 0.5 合． 昨晩の残り． 鶏レバ・ピーマン・ニラ・ニンニク茎・エノキダケの 炒めもの．
  • 昼 (1220): 弁当． 米 0.8 合． 朝と同じ．
  • 晩 (1920): 米 0.7 合． ナス・タマネギ・ニンニク茎・ニラの野菜炒め．

2001 年 08 月 14 日 (火)

• 0720 起床． 途中で一回ごそごそ起きたけど 合計 5 時間ぐらいは寝たかな．
• 朝飯・弁当の準備． 朝飯． シャワー．
• 0755 自宅発． 晴れ． ついに曇りの日々が終ってしまった． これからまた暑くなるんだろうか． 特急に乗って 0802 京急平和島発． 普通に乗り換えて 0830 京急上大岡発． 0845 研究所着． さすが 8 月． 日が照ると朝からぢりぢり暑い．
• 昨日の熱帯林で少しヘンな Kolmogorv-Smirnov 検定やってる甲山さんからメイル． そりゃあ， 生起確率ゼロの標本が得られれば 帰無仮説は必ず棄却されますよ． どんな検定法使っても．
• 確率論的モデルに誤差が混じってる場合は， なンか確率の複合分布 (compound distribution) 考えんといかんなぁ …… と前から考えてたんだけど， 具体的にはどうなるのかよくわからなかった． 今回つらつらと考えてみるに， あたかも誤差分布が中心的であるかのように 考えてみて， その平均値 m が別の分布 F(m) に従う， とすればいいんでないかなぁ …… などと思いついた．
• しかし， これって単に分布の平均が確率分布である， というモデルは想像しやすいだけのことで …… 逆に確率論的モデルの分布の全部のパラメーター を平均を含んだものに無理矢理翻訳して (というのも， たいていの確率分布はパラメーターを組み合わせると 平均になるんで …… しかし陰関数になることが多いか?)， そのあたりが何か誤差分布にしたがう， としてもかまわないわけだ．
• けっきょく， これって超パラメーターを考えるベイズ推定だよね． きっと． 超パラメーターの分布もデータから推定してみせる， というアレ． えーと empirical Baysian estimation， かな． ベイズな推定やったことないんで あまり確信もてないんだが ……
• 値段は数百万円するけど， いったい誰が使うのだろう …… という生態系変動領域 Sun Workstation で bash が使えるようになったので， ちょっとログインして， あちこちごそごそと覗いてみる …… ああ， Solaris ってやっぱ Linux とはだいぶ違うね． う? uname -a とすると SunOS 5.8 と出るな． いいのかな? よくわからん．
• ぢりぢりと危機的状況においこまれつつある PipeTree 改造に着手する． 何が危機的といって， これから大々的な改造しなければならないのである． なぜならいまの設計はかなりぼろぼろだから．
• そもそもパイプモデルというよくわからぬ 「原理」というか信仰にしたがって 樹木を作るようにしてるわけなんだけど …… だからといって 全体をこういう階層構造にする必要はないのである．

	森林 ┬ 樹木 ┬ パイプ ┬ 葉群	     ├ 樹木 ├ パイプ ├ パイプ断片	     │ 樹木 │ パイプ │ パイプ断片	     └ 樹木 └ パイプ └ パイプ断片

• いや， もしホントにそういう信仰だけに生きるのであるなら， これで教義どーりまことにけっこうなモデルである． ぢつは PipeTree プロジェクトは， まずはそこから始まった． しかしそういう狂気に染まってしまうのもイヤなんで (とはいえ，世の中には信仰への誓いによって 平和と思索の静かなる日々を実現してるヒトたちもいるらしい)， ケガれた観測データとの対応を考えねばならない．
• ということで， 1999 年秋 → 2000 年 11 月と PipeTree を 修繕を重ねていくうちに， 上の階層構造をみょーなカタチでまたがるような 構造が追加されていったのである． たとえば幹を形成する Pipe 断片だけをたばねた Trunk とか， Pipe 断片の位置情報だけを記憶している Skelton とか …… つまり上の図で示している世界の秩序は とっくに破綻しちゃってるわけだ．
• ではどうすればよいか? あとからいかなる宗教的解釈が付与されるにせよ， 樹木の地上より上に出てる 構造を生成するのに必要なクラスは， ふつーに考えるなら， たぶん以下のふたつだけであろう．
  • 樹状構造の伸長方向と分岐を決定する たくさんの「先端頂点」
  • 先端頂点が通過した場所を確保する 「断片」
  …… むろん， このほかにも葉群みたいなのも必要だろうし， あるいは LIGNUM の構造 に見られるような「分岐点」あるいは「節点」 (これらは断片と断片を接続してるもの) が必要になってくる場合もあるかもしれない (たとえば枝みたいな構造が時間とともに変形するような場合)． ともあれ， 樹状構造の基本的なホネぐみだけを生成 したいのであれば， 先端頂点と断片， もしくは先端頂点と節点の組み合わせでよい．
• パイプモデルでいう「パイプ」ごとき構造 (?) は 先端から根本の方向にむかって断片をたどるだけで (断片には極性がある)， パイプなる架空概念の形状は一意に決まってしまう． そしてこれにはとくに unique な機能があるわけでもない． ならば独立したクラスとして定義してやる価値はない．
• これまでの PipeTree の Pipe クラスは先端頂点と断片 (の一部) の 混じったものであった． これを「先端頂点」「断片」といった機能ごとに 独立したクラスとして分割しなければならない．
• と， 考えつつコードをどんどん削除していく． またあとでいろいろと書き足さんといかんのだろうけど．
• りすとらもあまり進んでないけど昼飯． とりあえず先端頂点だけの樹木でも作ってみるかなぁ．
• 午後も同じ作業． なんか撤退を繰り返してるようでいらいらするなぁ．
• べつにキレたわけではないけど， 削除に削除を重ねていたらこれだけが残った．

	森林 ┬	     ├
	     │
	     └

とほほ． さてさて， どうなるんだろうね．
• 新しいコードを書く気も失せたので， 甲山さん要求するところのオオシラビソシミュレイターの 満たすべき用件を整理してみる． いままでいただいたメモなど見直してみる．
• よくわからんのは先端部がどういうふうに死ぬか， 先端部に形成される断片の形状はどうなるか， というところだ． えーと， 生産したらすぐに針葉の呼吸をさっぴけ， か． もし， そういうことやったらどうなるんだろうか?
• 私がまだ九州大学の 4 年生やってたころに 書かれたメモの数値を信頼したふりをして そのあたりを計算してみる． 絶体絶命開空度 5.3%． ホントかね．

[葉っぱでの生産速度マイナス呼吸速度] 開空度が 5% あたりより低くなると (つまり暗くなると)， 稼ぎより消費のほうが多くなって死ぬ， という設定のようです． 甲山メモ 1991 年 (10 年モノですな) にあれこれと列挙されてる 信頼性よくわからぬ数値を組み合わせて概算． gnuplot script

(2001.09.05 後記: 上の計算と図はまったくの間違い )
• その他「当年生パイプの太さ-長さ関係はどうなってるんですか?」 うんぬんといった， 夏休み子供電話相談室ふう質問を列挙した メイルを北大に送る．
• BookPC 今日も安定してやがるな． やはり Buffalo の 128MB SDRAM の問題だったのかな …… でも 256MB 増設前は順調だったんだよね． 奇怪な．
• 甲山さんのアイデアについて 質問した以外の部分については理解できたんで (お望みどーりに動くかどうかはわからんけど)， ふたたびプログラムのほうにもどって考える． 何もかもきれいさっぱりと削除したあとは， 自由にアイデアが出るなぁ． さーて， どうしてくれよう．
• 新しいコード書きはやめといて， 帰宅することにする． 1740 研究所発． 普通に乗って 1755 京急杉田発． 特急に乗り換えて 1800 京急上大岡発． 1840 帰宅．
• 早くに帰ったので掃除機をかける． ダニ対策である．
• 旧科技庁系の健康保健組合に強制的に入らされてるんだけど …… 30 周年記念品とやらで 体脂肪計つき体重計なんぞ頼みもしないのに送って来やがった． うちにはもう体重計あるんだよね (最近乗ってないけど)． 未開封のまま Yahoo! オークションで売り飛ばそうかな． 希望小売価格を調べると …… なンと 16000 円． 保険料とりすぎてんぢゃないの?
• あ， 気象庁のサーヴァー， やられましたか ． 他もやられてるみたいね． お役所系はどこもまぬけだねぇ． 本格的な襲撃は 15 日か?
• 早くに帰ったので「ニセ」アカマツの続きでもやろう， と思ってるんだけど， どうにも眠い． 眠らずに起きてるだけ状態．
• とゆーことで， さっさと寝ますかね．
• 今日の食卓
  • 朝 (0730): 米 0.7 合． 鶏レバ・ピーマン・ニラ・ニンニク茎・エノキダケの 炒めものの残り．
  • 昼 (1220): 弁当． 米 0.7 合． ナス・タマネギ・ニンニク茎・ニラの野菜炒めの残り．
  • 晩 (2010): 米 0.7 合． 朝昼の残り． 残りものばかり食ってた一日であった．

2001 年 08 月 15 日 (水)

• 0600 起床． 昨晩は 2300 すぎに寝た． ひさしぶりに 6 時間以上まとめて眠れた．
• 朝飯・弁当の準備． 朝飯． シャワー． コーヒー． よく寝たあとはコーヒーがうまい．
• 0650 自宅発． おお， 今朝はまた素晴しい曇天． 最高気温は 30 ℃以下だろう． 特急に乗って 0657 京急平和島発． 普通に乗り換えて 0725 京急上大岡発． 0735 研究所着．
• とりあえず， 昨日から考えてる Abies PipeTree の現段階でのアイデアをメモにしてみる． といっても資源配分のやりかたを変える， というだけなんだが． 他の細かい部分については甲山さんからの返事待ち．
• 上の Abies ペイジ見て初めて知ったんだけど， オオシラビソって アオモリトドマツ の別名だったのか! ああ， 無知無知． ま， モデル屋やってるぶんには差しつかえないけど．
• お役所系まぬけサーヴァーやられ続報． 消防庁， 通信途絶． 靖国神社， アクセス異常集中． なぜか本願寺もやられたとか．
• 簡単なアイデアメモはひととーりできあがり． 急に眠くなる． 早起きのむくいか?
• Abies PipeTree はいったんおいて， 「ニセ」アカマツの図の作成に取り組む．
• ぜんぜん関係ないけど， makedepend というヘッダーファイルの依存関係を 調べてくれるコマンドがある． これを Makefile の中に

	depend:
        	makedepend ${INCLUDE_OPT} -Y ${SRCS} 2>/dev/null

と書いておいて (-Y という指定は標準ヘッダーファイル置場 /usr/include を読ませないため， /dev/null はその warning を捨てるため …… こういうやりかたがよいのかどうかわからん)， make depend すると

	unival.o: unival.h
	main.o: ml.h typedef.h plugin.h parameter.h matrix.h
	proc_arg_sample.o: ml.h typedef.h plugin.h parameter.h unival.h
	newton.o: ml.h typedef.h plugin.h parameter.h matrix.h
	parameter.o: parameter.h unival.h
	matrix.o: matrix.h ml.h typedef.h plugin.h parameter.h
	likelihood.o: ml.h typedef.h plugin.h parameter.h
	output_plot.o: ml.h typedef.h plugin.h parameter.h

というような必要にして十分な依存関係を列挙した行を 自動的に生成して， Makefile 自身に追加してくれるのである．
• う， makedepend で遊んでいて午前中終ってしまった …… 反省しつつ昼飯． ネット上の東京大空襲情報を調べてみたり ……
• お役所系まぬけサーヴァーやられ続報． 文部科学省， 沈黙． 旧科学技術庁は昨年に引続き二度目の失態． ま， ここ独房群 (旧科学技術庁系) の腐れたファイアーウォール見てれば 突破されるのも当り前という気もするが． ただたんにハコつけただけで防御したつもりになってるのが お役人的管理者のめでたさである．
• 昼過ぎても「ニセ」アカマツ作図の続き．
• 「ニセ」アカマツに関して， 以前からほうっておいた問題に取り組む． 実験につかわれた (ほんもの) アカマツの齢の定義の部分である． 「ニセ」アカマツの中では， 毎年成長シーズンの終り (秋ですかね) に 一歳年齢が増える， としているので …… と考えつつ， 観測年と観測データーとシミュレイターで 整合性が取れるようにする．
• 途中で， 結果出力に関するヘンなバグを見つけて直す． よくまぁ， いままで問題無く動いていたもんだ．
• 異端審問官も兼任しておられる某宣教師殿から 職務質問メイルが来信 (しかも恐ろしいことに一通目は 意図を隠蔽した誘導尋問になっていたのである!)． 西日本某闇ネットを不調に陥れた犯人を捜査中とのことで， 幸運にも私はアリバイが認められ 容疑者リストからは外れたようである． うかつな返信していれば 今ごろ火あぶりになっていたところだ．
• 5 枚予定している「ニセ」アカマツ図のうち， 4 枚は関連プログラムからデータ受け取って 仕上がりに至る処理を自動化できた． 腹へったし， 撤退しますかね． 1745 研究所発． 普通に乗って 1755 京急杉田発． 特急に乗り換えて 1800 京急上大岡発． 1840 帰宅．
• 晩飯食ってから， 眠気をこらえてお役所サーヴァーのやられ具合を調べていると …… 意外や意外， いったんは完全にやられてた文部科学省， 柔軟にも機動防御に転じている． netcraft の調査によると ……

	        OS, Web Server and Hosting History for www.mext.go.jp         
	┌─────┬────┬──────┬────────┬──────┐
	│    OS    │ Server │Last changed│   IP address   │  Netblock  │
	│          │        │            │                │   Owner    │
	├─────┼────┼──────┼────────┼──────┤
	│          │Simple, │            │                │            │
	│          │Secure  │            │                │      Hitach│
	│Solaris   │Web     │ 14-Aug-2001│ 210.174.162.214│     Capital│
	│          │Server  │            │                │ Corporation│
	│          │1.1     │            │                │            │
	├─────┼────┼──────┼────────┼──────┤
	│          │Apache/ │            │                │      Hitach│
	│Solaris   │1.3.12  │ 11-Apr-2001│ 210.174.162.214│     Capital│
	│          │(Unix)  │            │                │ Corporation│
	├─────┼────┼──────┼────────┼──────┤
	  ...

…… となっていて， つまり 14 日 (日本時間の今日) 攻撃を受けたあとになって Simple, Secure Web Server なんて聞いたこともないソフトウェアに切り替えているのである (この httpd についてはいくら調べても情報が出て来なかった …… 開発したらしい http://www.raptor.com/ というところには 接続できない状態だし)． つまり攻撃側が対応できないような無名のものに 交換してしまったわけだ． あらかじめこの展開を想定してたのかな? 今回ばかりはちょっとスゴイね． (後記：とホメてやったのに 22 時すぎまた落とされてる …… この程度の対応では焼け石に水か)．
• などと， よけーなことばかり調べてるもんだから， 仕事がススまん．
• 今日の食卓
  • 朝 (0620): 米 0.7 合． ゴーヤ・ナス・タマネギ・シイタケの炒めもの． 夏は沖縄・鹿児島産のゴーヤ (ニガウリ) が安い． 旧東海道商店街の八百屋だと 一本 (30cm くらい) が 100 円 (税込み)．
  • 昼 (1215): 弁当． 米 0.7 合． 朝と同じ． 近ごろの飯炊きで工夫していること． 米をよく研ぎよく洗う． こうすると弁当箱につめても コメのニオイがほとんどしない．
  • 晩 (1900): 米 0.7 合． 炒めものの残りで炒飯．

2001 年 08 月 16 日 (木)

• 0700 起床． 7 時間ぐらい寝たか． 平日なのによく寝た．
• 朝飯・弁当の準備． 朝飯． コーヒー． シャワー．
• 0750 自宅発． オホーツク海に居座る高気圧の 特急に乗って 0801 京急平和島発． 普通に乗り換えて 0825 京急上大岡発． 0840 研究所着．
• 昨日， 帰宅前に仕掛けておいた 「ニセ」アカマツ独立反復試行終ってる． 私が帰宅してから 3 時間以上も計算してたのか? …… あ， 作図するとき以外は 必要ない「ニセ」枝生成の三次元計算をやってる． 無駄な寄り道をやらせてしまった． しかし結果は前と変わらない． 問題ナシ．
• 朝から 2 時間ちかくも北大闇サーヴァー hosho の設定修正に苦闘する． というのも， 以前に他ならぬ私が rpm package の 依存関係を無茶苦茶にしてしまったからである …… ようやく hosho でも apt による 更新ができるようになった．
• 5 枚予定してる「ニセ」アカマツの図のうち， ひとつ残していた図 4 に取り組む． これは「ニセ」アカマツ一試行の 時間発展の例を示すものである．
• 作図プログラムは単にシミュレイターで 描画用データを生成し， その 1 ステップごとに 描画プログラム・画像変換プログラムを を呼び出すだけの簡単なもの． しかし作図視点位置を決めたり， A4 の紙におさめるような画像変換の 試行錯誤で手間取る．
• 1200 すぎにだいたいできたんで 昼飯の弁当など食う．
• なんか今日はもうネット空襲は終ったみたいね．
• 午後もとりあえず作図自動化のつづき． 図 4 はこんなもんか．

[時間発展] 「ニセ」アカマツシミュレイター (作図用データ生成) → POV-Ray (三次元レンダラー) → ImageMagick (cropping，EPS に画像変換) → pLaTeX (配列，label 貼り) という工程で生成． 枝のように見えるのは枝ではなくて 「ニセ」枝である． 「枝は定義されてないからあるはずないんだ …… 枝なんて幻覚だっ」 と読者が錯乱してくれたら楽しいなぁ． ようこそ， 「ニセ」わーるどへ．

• ということで， figure ディレクトリで make とだけ指示すれば fig.1 ... 5 の清書版を順々に生成していくようになった．
• 図の説明をぢりぢりと書く …… けど終らんうちに撤退． 1740 研究所発． 普通に乗って 1755 京急杉田発． 特急に乗り換えて 1800 京急上大岡発． 普通に乗り換えて 1822 京急川崎発． 1850 帰宅．
• ここ二日ほどやや多めに寝てるので， 晩飯のせいで消化器官に血液取られた脳が 機能低下したあげくにいきなり昏睡状態に陥る， ということはないんだけど …… どうにも眠いな． 晩飯多いのかな． 野菜のカサが多いから わが栄養器官は苦闘してるような気もする． よく噛んで食べたんだけどな．
• ま， 図はあとは説明つければいいとして， 問題は表だな． ちょっとめんどいかも．
• Abies PipeTree は止まったまま． 大丈夫なのか．
• 今日の食卓
  • 朝 (0720): 米 0.7 合． ゴーヤ・タマネギ・シイタケ・ニンニク茎の 炒めもの．
  • 昼 (1220): 弁当． 米 0.7 合． 朝と同じ．
  • 晩 (1900): 米 0.7 合． ニラ・ゴーヤ・タマネギ・シイタケ・ニンニク茎の 炒飯．

2001 年 08 月 17 日 (金)

• 0700 起床． 途中でいっぺん起きたけど， それでも合計 7 時間ぐらいは寝たか． 夏の月→金に独房群がよいを続けるとなると， これぐらいは休んだほうがよさそうである．
• 朝飯・弁当の準備． 朝飯． コーヒー． シャワー．
• 0750 自宅発． 今朝は残念ながら晴． 特急に乗って 0701 京急平和島発． 普通に乗り換えて 0728 京急上大岡発． 0745 研究所着．
• 「ニセ」アカマツ続けようとしたんだけど， そのまえにちょっと Abies PipeTree の メモを直す． 枝の先端部における伸長・肥大・分枝が ある資源量に対して一意に決まる， という甲山新構想を真にウケるなら， 旧 PipeTree の Pipelet 形成につけてた 意味不明な制約を外すことができるハズ …… というもの． はてさて， どうなるか． PipeTree の Pipe とは何なのか， ますますわからなくなってきた．
• 「ニセ」アカマツ図の説明つづき．
• 図の説明できた …… と思ったら午前終了． どうしてこんなに時間かかるんだろうね． 誰かたすけてー
• 昼飯を食いつつ， ネット上をふらふら． matpack って何だか面白そうだな …… 独逸なる数値計算・描画の C++ ライブラリです． しかし 28MB もあるんで， ここの腐りかけた回線ではダウンロード難しそう． 自宅の ADSL のほうがよほど速い．
• 昼飯食ってからもしつこく「ニセ」アカマツの続き． 図はだいたいできたとして (ただしもう一枚追加? の懸念アリ)， 表はどうしたもんだろうか． えーと， 個体ごとにパラメーター並べるのと， 生産力を決めるファクターを並べたものと …… 分割するかひとつにするか? …… ちょっと試作してみないとわからんなぁ．
• 問題は作表に必要な数値の出力だ． 計算結果を再加工するのはちょっと複雑化するか? とはいえ， シミュレイターを改造するのはイヤだ …… ふと思いついたのだけど， 「ニセ」アカマツの構造は簡単で，

	 Sim ┬ 「ニセ」アカマツ ┬ 「ニセ」茎	     ├ 「ニセ」アカマツ ├ 「ニセ」茎	     │     :            │     :	     └ 「ニセ」アカマツ └ 「ニセ」茎

となっており， 最上層の Sim だけが Perl スクリプト， それ以外の階層は Perl モジュールとなっている． で， この Sim というのは単なる周旋屋であり， 「ニセ」アカマツを一本ずつ呼び出しては 「つぎはコレやってください」 と仕事を順番にまわしていくだけで (蛇足ながら「ニセ」アカマツも請け負った仕事の何割かは 下の階層に丸投げしてしまう)， まっとうな計算は何もやっていない．
• ということは， Sim を詐称する別の Perl スクリプトが 「ニセ」アカマツたちを乗っ取ってしまえば …… そう， 8 個体の「ニセ」アカマツを自在にあやつり， 作表に必要な結果だけを計算させることができる． うーむ， しかし Redpine::TableLineOutput というような関数の追加は必要かも …… あ， それは新しいファイルにちょいと書けばいいのか． よし， このニセ「ニセ」アカマツ作戦でいこう．
• …… 1715 よーやく表とその説明できた． 「ニセ」アカマツたちは 上記「ニセ」Sim に操作されて 作表のためだけに動作するのである． その結果を受け取って LaTeX でカッコいい組版をするところで苦闘してしまった．
• いつものごとく make すれば清書版が出て来るようにする． この方式のありがたいところは， やりなおしがラク， ということだ． たとえば元データに入力間違いあって 表の中の数値が何もかも変わってしまう …… という厳しい状況になったとしよう． そういうときにも慌てず騒がず make だけでよいのである． 更新されたデータで再構築された 「ニセ」アカマツたちが呼び出され， 必要な作業だけやって …… となるわけ．
• ひとくぎりついたから撤退しますかね． 1740 研究所発． 普通に乗って 1750 京急杉田発． 快特に乗り換えて 1755 京急上大岡発． 1815 京急川崎着． 川崎ヨドバシにふらふらと． よく見るとあまり品揃えのよくない 書籍コーナーで迷ったすえ 「Python 入門」 を買う． 税込 5000 円ほどだけど， 店員のにーちゃんにむけて ポイントカードなる合成樹脂断片を しゅっと飛ばすと無料になった． わーい …… ってすでにヨドバシに 5 万円以上を貢いでるってコトぢゃないか． ああ， ヨドバシの精神的奴隷状態． 普通に乗り換えて 1831 京急川崎発． 1850 帰宅．
• 晩飯を食いつつ NHK ラジオ第一 19 時ニュース を聞く． 今週の担当は凛とした美声の 大迫あゆみ さんである． どきどき
• 同番組は 1855 に終る． しかし， ここから関東在住数百万の大迫ふぁんが最も緊張を強いられる 3 分間が始まる． この地方の気象情報である． 引き続き大迫さんが担当して気象協会のおじさんと やりとりするものである．
• これはよくある天気予報とは一線を画している． 今日の天気の概況ならびに明日とそれ以降の予報を 3 分間に押しこむというだけでもたいへんそうなのに …… なんと， この 1855 天気情報では， 単調になりがちなこういう報道を 少しでも楽しいものにしようと， 大迫さんが毎日あれこれと工夫されるのである．
• 19 時ニュースのほうはあらかじめ渡されている台本を 感情を交えずに淡々と読む (ことを強制されてる?)， というかんぢなのに対して， こちらは現在の天気に関する所感をアドリブのように 自分のコトバを用いつつも， 内容を精選してごくごく短くに要約するところから始まる． また 「それでは関東南部ではどうでしょうか?」 といった先方との受け答えも， なんというか型どおりにならぬように イントネイションをつけるなど努力されている．
• われわれはそれまでの 55 分間とは一転した演出に驚き， 「ああ，今日も無事に終るだろうか」 とはらはらしつつ， 局限された時間内に定型的な内容を報道をするにあたっても， 聴き手に親しみをもってもらおうとする 大迫さんの心意気に深い感銘をうけるのである．
• う， 考えてみれば， ここ 10 日以上， ラジオ以外から人間の声を聞くことがなかったような気がする …… まぁ， 独房群とはそういうところだ．
• 今日の食卓
  • 朝 (0715): 米 0.7 合． ゴーヤ・タマネギ・ニラ・エノキダケの 炒めもの． 昨日， 100 円ショップで ウスターソースを ひさしぶりに 買ってしまったので， それと塩胡椒・トウガラシで味つけ．
  • 昼 (1210): 弁当． 米 0.7 合． 朝と同じ．
  • 晩 (1900): 米 0.7 合． 朝と同じ．

2001 年 08 月 18 日 (土)

• 0900 起床． 週末のわりには終末でない起床時間．
• 昨日は真夜中までプログラムを直していた． めずらしく夜間に仕事， というわけでもなく， 当面のプロジェクトにはまるで関係ない mlfitting (多次元非線型モデル対応の最尤推定プログラム) の手直しである．
• 直したかったのは STL 依存性を削れるところは削りたい， という部分である． とくに整数・倍精度浮動小数点数・文字列・ブール型に 対応している汎用要素 Unival で 必要もないのに map や strstream といったのを使ってる部分． map は enum+配列に， strstream は C++ ではなく C の標準関数である sprintf だの strtod だのに置換してみる．
• 私は C++ にはここ数年呪われつづけているので， 数画面ぶんの呪われエラー警告文が出たりしても 冷静に対応できる． しかしながら， C++ ではない C の世界， というのはよくわかっていなくて， core ファイルどかんと吐いて異常終了されようものなら かなりアセってしまう．
• 今朝もその手直しを続けて， ようやく問題なく動くようになった． STL 依存を少なくすることでコンパイルの速度は かなり向上した． しかし実行速度については …… この mlfitting では全然かわらないね． さて， ここで試験できた置換技法を シミュレイション用久保ライブラリに適用してみると， 何かいいコトあるだろうか?
• 昨日，川崎ヨドバシで買った ヘビの表紙の 「Python 入門」 にひととーり目を通してみる． これもまた当面のぎょーむには何の関係もない．
• 日本語で書かれた Python 本は少なく， これとこれの続刊である 「Python プログラミング」 が主たるところであろうか． どちらも (というか原書は一冊) 1996 年刊行でちと古くて， Python-1.3 に準拠してる． 現在の主流は 1.5 ですでに 2.x も出回っている． なお 「初めての Python」 てのもあるけど， これは，ちょっと という書評もある．
• で， ヘビ本ですごいのは (と言いますか私がこの本を買った理由は) 最初のほうから 「永続オブジェクト」 なんかを解説してるところである． おぶぢぇくと指向プログラミングも結構ですが， やっててイヤなのは面倒な構造のクラスの インスタンスがプログラム終了と同時に消えちまう， という点である．
• 現在の森林や樹木シミュレイターだのは インスタンスを再現するのに必要な情報を Begin-End で囲まれた階層構造の中に押しこんで テキストファイルに出したり読みこんだりしている． そのために必要な Filereader クラスも 自作ライブラリにある． しかし入出力にはえらく時間かかり， ライブラリ使ってなおそのコーディングは面倒である．
• このあたりを解決するためには， 永続オブジェクトか外部データベイスを使う， というのが上策のようである． で， この 5 年前刊行のヘビ本は初心者むけなのに， それが解説されているのである． いや， まぁたいていの OO な言語では 誰かが永続化の仕掛けを作ってネット上で公開してはいますよ (あの Perl にも !)． しかし C++ も Ruby も入門書みてもそんなことは まるっきり言及されてもいない．
• で， Python では標準ライブラリで簡単に永続オブジェクト作れる， と． たいへんに結構． どうも全部テキストに直しちゃうみたいで (最近のやつは違うのかもしれないけど)， 速度の点ではあまり速くはないかも． あと， つーどくしてわかったのは， Python が短期間開発のための プロトタイピング言語としての側面を 強調している点である． つまり必要な処理は後からどんどん C/C++ で置き換えてくれ， という方針がはっきりしている． その機能も充実している． Perl と比較した場合に こういったところが特徴的なのではなかろーか ……
• …… と思ったら， 昨日の「ぎょーむ日誌」の Python 本うんぬんの記述を 摘発した異端審問官もとい 某宣教師殿が Python ではなく Ruby を， との布教活動． ごくすなおに Python ではなく Ruby， とするご利益を伺ってみたのだけど， この質問はある種の不敬罪のようなものに該当するらしく， 何か意味ありげな沈黙がただ返ってくるばかりである．
• Ruby 作者の手になる 「オブジェクト指向スクリプト言語 Ruby」 に掲載されている比較によると， Ruby は Unix 的で 「オブジェクト指向に基づいて設計されたスクリプト言語」， Python は OS からは中立的で 「スクリプト言語にオブジェクト指向機能を追加したもの」 …… うーむ．
• どの言語にせよ， 苦労しながらひととーり使ってみないとわからん， ということかな． たとえば， 全部 Perl で書いてしまった「ニセ」アカマツを 書き直すとしたら， Python が良いのか Ruby がよいのか (はたまたどうなるか全くわからぬ Perl6 か)? これなんかは実際やってみないと判定できそうにないなぁ．
• けっきょく仕事はススまぬ一日でした． ライブラリの細かい手直し少しやっただけ． 日が沈んでからふらふらと外に出る． 奇跡のように涼しい． 大森町商店街のほうに行ってみる． 電気代支払う． 3687 円． 冷房代のせいで高くなっている．
• 今日の食卓
  • 朝 (0930): 蕎麦．
  • 昼 (1330): スパゲッティー． ゴーヤ・ニンニクのソース．
  • 晩 (2130): 米 0.9 合． ゴーヤ・タマネギ・ニンジン・マイタケ・卵 の炒めもの．

2001 年 08 月 19 日 (日)

• 0800 起床． やはり平日からよく寝ておけば 週末だからといってバテたりしないのかも．
• と言いつつ， 午前中は半分寝てるような状態． Harvard Forest で感光紙を貼りつけ作業を終えて 機嫌のよくなった浦口さんが 雑談メイルのやりとりに応じてくださる． なぜか麻布高校山岳部 1959 年 夏山合宿映画のネット上で公開されてる 動画ファイルに関する議論．
• うーむ， このころは北アルプス国立公園内でも焚き火 OK だったのか． すごいすごい． キスリングとか個人装備は この 30 年後と変わらんな． とはいえ， 去年に槍沢から登ったときにはほとんど見掛けなかったけど． あの日本的な登山用ずた袋．
• 今朝の東京の最低気温 19.6 ℃!! 素晴しい． このままコトを荒立てずに秋になってくれ． ああ， しかし私の「夏休みの宿題」が ……
• ちょっとシミュレイション用ライブラリを手直し． 昼すぎてから「ニセ」アカマツつづき． 無謀な気がするんだけど， 結果・考察まで進めてみよう， という試み．
• ちょっと中断して申請していたパスポートを取りにいく． 1630 までに有楽町に行かねばならない．
• 1535 自宅発． 曇り． 涼しい． たいへんよろしい． 1545 京急平和島発． 品川で乗り換えて 1605 JR 有楽町着． 交通会館とやらで 15000 円の収入印紙というのを買い， それを貼りつけた用紙を差し出して 無事にパスポート受け取る． 赤い表紙． 2011 年 8 月まで有効． 用事おしまい．
• 前回来たときの顛末として， ここらにはろくな本屋が無いと書いたら， すぐ近くに 八重洲ブックセンター ありますよと 竹中さん から教えていただいたので， 今日はそちらを威力偵察していくつもり．
• 有楽町から JR 東京八重洲口までどれぐらいの距離なのか? あっ， うかつにも地図を持って来るのを忘れた． しかし線路ぞいに進めば自位置を失することもあるまい． 幸いにもデイパックには方位磁針つけてる． 方位北々東． 確認． 進出予想距離いっせん …… いや， はっぴゃくめーとる? 1610 交通会館出る． 外堀通という道路ぞいにすすむ …… 6 分後， とくに遭難もせずに八重洲ブックセンター着． ふう．

[黄金金次郎] 八重洲ブックセンターなる JR 東京八重洲口にある巨大書店． 店頭に金ぴか像を配したりするあたりが 東京人のセンスのわけのわからなさだ．

• おお， さすがにでかい …… と感心しつつ 3F 科学書おきばへ． 生物学関連ざっと見てから， 数学関連のほうに行ってみる． やはり， 浜松町の最大本屋 Dan や JR 横浜の有隣堂よりもよほど揃っている． もっともいんぱくとあったのは 「五次方程式の代数的一般解法」なる自費出版本． 一般解が存在しないとされていた従来の定説をくつがえした， か． すごすぎー
• 統計学コーナーはかなりじっくり見たけど買いたいものナシ． 計算機本のほうで 3 冊買う (後述)． 13000 円の出費． うう．
• 7F 輸入本置き場に行ってさらに驚愕． 科学書についてはこちらを見るべきだ． 統計学本しか見てないけど， 英語本置き場のほうがよほど良い． しかし値段がみょーに高い． もしやと思って帰宅してから Amazon.co.jp の値段と照合してみたら …… やっぱり高いぞ八重洲値段． たとえば， 前の版に比べると パラメーター推定方法などがずいぶんと実用的になった Matrix Population Models : Construction, Analysis, and Interpretation (Hal Caswell) が Amazon.co.jp 値段 7785 円なのに， 八重洲値段は 11000 円超． さらに八重洲 12000 円ほどの Numerical Methods of Statistics は Amazon 8984 円． 計算機本なんかでも比較してみたんだけど， この傾向は同じである． さらに Amazon.co.jp の場合， 5000 円以上買うと送料無料 (500 円の払い戻しまでアリ)． 大森の自宅から JR 東京まで往復すると 620 円．
• ということで， 私のようなびんぼー人にとっての 八重洲ブックセンターの利用法は: (1) JR 東京による機会があったら覗いてみる． (2) 科学書は輸入本置き場を見る． (3) しかし実際に購入するならネットから買う． ということになりそう．
• そういえば， 独房群からの帰路途中の京急上大岡にも 八重洲の出店があるけど， どういうかんぢなんだろうなあ． 今度しらべてみるか．
• 2 時間ちかく費した八重洲ブックセンター発． 1810． JR 東京地下の古本屋なぞのぞいてみる． 1826 JR 東京発． 品川で乗り換えて 1900 帰宅．
• 晩飯食ってから買って来た三冊の計算機本をぱらぱらとながめる．
• 「MySQL & mSQL」 ． データベイス本 (原書は 1999 年のもの)． 森林・樹木シミュレイターのデータを 動態コードと分離するなら， この RDBM か， ということで． C++ によるオブジェクト指向データベイスうんぬん という解説あり． むろん Perl，Python，Java，Ruby といった他の言語からも 呼び出す手段が準備されている．
• 先日買った 「Python 入門」 が面白かったので， その続刊 「Python プログラミング」 これらの本は 1996 年刊と古くて， 呪われな C++ との結合に関しては あまり詳しく書かれていない． ネットで調べるしかないか．
• まぁ， C の拡張モジュールが (Perl に比べて) 簡単にかけるみたいなんで， Python → C → C++ とつないでしまえばよいけど． どうせ， いまやろーかなぁと考えてるのは 自作の C++ な開空度計算ライブラリを 何かのスクリプト言語から呼び出すことなんで．
• このヘビ本， 永続的情報に一章さいている． 内容は shelve をつかった dbm との連携， か． どうもこれがメインのワザみたいね． 高速化を狙うなら cPickle 使えってことのようだ．
• 最後の一冊 「入門 Scheme」 ． Scheme は Lisp の一方言である． Linux 上では GNU が GUI 環境設定言語として 布教しつつある guile なんかが Scheme 実装の例である．
• なんでこんなのに興味を持ったかというと， やはり永続オブジェクトがらみである． 以前に 「普通のやつらの上を行け ---Beating the Averages---」 なる Lisp 賞讃文書を読んだ． この著者いわく， Lisp ってのはラテン語みたいなもんじゃなくて ネット通販サイトなんか作るときにもすごーく役に立つんだぜ， ひひひ． へえ， なんであんな計算科学者用の (((()))) という見てくれの言語なんかが 商売の役にたつんだろう …… と思って読み進めてみたら， どうも複雑なデータ構造有する永続オブジェクトの 書き出し・読みこみ素早く保守と拡張が簡単， と言いたいんではなかろーかと推測できた．
• こういう機能に関しては …… なるほど Perl,Python,Ruby にも eval なんてのはあるけど， おそらく Lisp 一族はそのあたりがキモなんで 性能がかなりよい， とか． どうなんだろう．
• メイラー Mew を動かしている Emacs の emacs-lisp も含めて， 私は Lisp はほとんどわからない． しかし， 思い返せば Lisp の影響うけたとある言語を かなり駆使していた時期がある． Mathematica である．
• Mathematica はあのふにゃふにゃしたフロントエンドは横にどけて， リスト {...{...{...}}} を駆使したコードで Mathematica カーネルを直接たたく， というのがもっとも効果的な使い方だった (ただし，今はこの懐かしい言語を使っていない …… 少なからぬカネがかかるソフトウェアなんで)． で， たしかに永続的なオブジェクトを書き出すことは容易であり， 読みこんですぐに eval (とわざわざ明示的に指示する必要もなし) もできた …… あれみたいなもんかな．
• …… などなどと， 当面のぎょーむにはまったく関係ないことばかり書いてたら すでに 2300 すぎ! 明日の飯の準備もしてさっさと寝なくては． ああ， 結局この週末も仕事ススまなかったな． 夏休みの宿題， 間に合うのだろうか?
• 今日の食卓
  • 朝 (0830): 米 0.7 合． 昨晩の残り．
  • 昼 (1330): 蕎麦．
  • 晩 (1940): 米 0.7 合． ニンジン・タマネギ・ニラ・卵の炒飯．

2001 年 08 月 20 日 (月)

• 0700 起床． なぜか眠い． 7 時間ぐらい寝たはずなんだが．
• 朝飯・弁当の準備． あれ炊飯終ってない． タイマーの設定間違えたか． シャワー． 朝飯． コーヒー． 朝のラジオ， なんかヘンだ …… 腕時計で時刻を確認． 0540．
• …… これは， つまり …… 昨晩， 半分寝てるアタマで目覚し時計を設定したときに， 間違えて時刻を 2 時間すすめてしまった!! うう． 訂正． 0500 起床． 5 時間しか寝てない．
• とはいえ， 気温高くなる前に移動したいので 独房群にむかう． 0555 自宅発． 曇だけど日がさしている． 特急に乗って 0605 京急平和島発． 普通に乗り換えて 0635 京急上大岡発． 0650 研究所着． いつもより 2 時間早くついてしまった．
• 眠い． とりあえずメイル書いたり， matpack コンパイルしてみたり．
• matpack …… こりゃあ， カンタンだわ． 八木アンテナのローブを描画させる こういうデモプログラムなんかも入っていたんだけど，

        #include "matpack.h"
        int main (void) 
        {
                Scene scene;
                ...
                scene.Transform(...)
                ...
                scene.AddLight( Vector3D(-0.5,-4,0.5),  white, 1);
                ...
                Vector3D camera_position( distance, 1.2*distance, 1.4*distance), 
                         look_direction (-distance,-1.2*distance,-1.4*distance);
                scene.Look(camera_position, look_direction, FieldOfView(45), 0);
                ...
                MpImage image(size,size);
                scene.Open(image);
                ...
                scene.Show();
                ...
        }

…… の最後に

                image.WritePngFile("isosurf-demo3.png");

と追加するだけで， 下のような図を PNG ファイルとして出力してくれる．

出力先は各種ビットマップだけでなく， Postscript なども含まれている． なかなかよさそう．
• おっと， こんなことやって遊んでいたらもう 0900 か． 朝から何をやってんだか．
• とりあえず Abis PipeTree を側面から支援するはずの VitPolygon クラスを定義してみる．
• 途中で 1 時間ほど， これまた仕事と関係ない， Perl モジュールの導入というか perl -MCPAN -e 'install ...' が使いこなせなくてぢたばたする．
• VitPolygon の基底たる Vit で string を使った名前の登録なんかをやってる． なんでこんなモノいれたんだろう? わけわからんけど， 基底にさかのぼって削除． ついでに Location も 純粋仮想関数として定義する必要ナシ． モノを考えずに作ってるというか ……
• いや， Vit::Location は必要だったみたいだ． うう．
• 早起きしたわりには， というか早起きしためにというか …… あまり仕事のはかどらなかった長い午前終了． 昼飯食う．
• 甲山さんからアロメトリーメイル． Abies では (シュート長さ) = 28 * (シュート太さ)^1.1 ? ほとんど線形だな ……
• 昼からも三次元ライブラリ関係の改善． VitPolygon という変な名前はやめて VitTriangle とする． 核心部は三次元空間にひらひらと飛ぶ三角形を ビームで撃ったときに当たったかハズれたか， の判定である． いずれ詳しく説明するとして， やることは 3 × 3 の逆行列の計算である． こんなかんぢ．

	...        // Preparation to generate 3x3 inverse matrix
        double  x0( p0.X() ), y0( p0.Y() ), z0( p0.Z() );
        double  x1( p1.X() ), y1( p1.Y() ), z1( p1.Z() );
        double  x2( p2.X() ), y2( p2.Y() ), z2( p2.Z() );
        double  xy01_10( x0 * y1 - x1 * y0 );
        double  xy12_21( x1 * y2 - x2 * y1 );
        double  xy20_02( x2 * y0 - x0 * y2 );
        double  yz01_10( y0 * z1 - y1 * z0 );
        double  yz12_21( y1 * z2 - y2 * z1 );
        double  yz20_02( y2 * z0 - y0 * z2 );
        double  zx01_10( z0 * x1 - z1 * x0 );
        double  zx12_21( z1 * x2 - z2 * x1 );
        double  zx20_02( z2 * x0 - z0 * x2 );
        // Determinant
        double  det( x0 * yz12_21 + x1 * yz20_02 + x2 * yz01_10 );
	...
        // Vector composition
        Xyz     gaze_vector( gaze->NormalizedDirection() );
        double  gx( gaze_vector.X() );
        double  gy( gaze_vector.Y() );
        double  gz( gaze_vector.Z() );
        double  a0( yz12_21 * gx + zx12_21 * gy + xy12_21 * gz );
        double  a1( yz20_02 * gx + zx20_02 * gy + xy20_02 * gz );
        double  a2( yz01_10 * gx + zx01_10 * gy + xy01_10 * gz );
	...

• 私はふだんはもうすこし読みやすいコード書くんだけど， この部分はあえてこう書いた． そう， いつかの失敗を彷彿とさせるんだけど， つねに愚かしくなりがちな高速化のためである． かかる自己満足的な努力はあれこれあって …… たとえば， みればわかるように一度も割り算を使っていない， とか (しかしあとで小規模な実験してみたら， ぢつはこういう工夫はあまり意味なかった)．
• さてと， あとは{ a0, a1, a2 } の符号がことごとく同じであれば「当たり」， でなければ「ハズれ」となる …… んだけど …… あれれ? うまくいかないなぁ．
• バグとれぬまま無念の撤退． 1800 研究所発． 普通に乗って 1810 京急杉田発． 特急に乗り換えて 1819 京急上大岡発． 1900 帰宅．
• 晩飯終ってからも気合のバグとり．
• 2230 よーやくまともに動いた． すごくつまらん間違いだった． 最後の { a0, a1, a2 } の算出でしくじっていた (上のコードは修整ずみ)．

[VitTriangle] 三角形の障害物 (緑)． 赤線はそれにぶつかってる「視線」． 黄線はわきをすりぬけたもの． これは POV-Ray による作図． (左図クリックすると拡大)

• こんなチンケな三角形のためにまるまる一日つぶしてしまった． 前途遼遠， 先ゆき暗雲． 嗚呼．
• しかも久しぶりに galeon やめて Netscape で見てみると …… 過去 10 日ほどの「ぎょーむ日誌」レイアウトかなりがたがた． これは Netscape が悪いんだい， とべそをかきつつ， ぢりぢり修整する． ああ， 睡眠時間が．
• (一言近況) 「ニセ」アカマツの難所は越えたけど， いまだ終らず． オオシラビソな PipeTree はドロ沼の予感． 大型で強い勢力の台風 11 号は接近中． いったいどうなるんでしょうね． (2001.08.20)
• 今日の食卓
  • 朝 (0530): 米 0.7 合． ゴーヤ・タマネギ・ニンジン・マイタケの 炒めもの．
  • 昼 (1210): 弁当． 米 0.7 合． 朝と同じ．
  • 晩 (1910): 米 0.7 合． 朝と同じ．