牧野の公開用日誌

つっても、非公開のを別につけているわけではない。

Copyright 1999-2002. Jun Makino

当面の予定

  冬学期日程: (計算天文学:月曜3限)
　授業開始　１０月　１日（火）
　授業終了　１２月２０日（金）
　冬季休業　１２月２１日（土）〜１月　５日（日）
　授業開始　　１月　６日（月）
　授業終了　　１月２７日（月）
　補講・試験　１月２８日（火）〜２月１７日（月）

2002/6/30

• お、江守君のインタビューの掲載が始まってる。

• AMNH にもどってくる。何故か今日は DHCP が動く。全く謎である。

• FreeBSD のファイルサーバを Linux の機械からマウントした時に、 Linux の NIC が GbE だとなんか死ぬという問題があったんだが、

  fs:/exportdisk /inportdisk nfs rw,bg,mountvers=3,rsize=8192,wsize=8192 0 0

  みたいに version 3 を指定したら動く、、、と思ってよく man page を読む と mountvers は RPC version で、 nfs version は nfsvers で指定するとい うことなんだが、、、なぜこれで動くようになるのかしら？

• オランダ人のものはオランダ人へ、、、というわけでもなく単に先送り にすることで同意。
• Vine で /etc/printcap をいじりまわす。なぜか lpd を kill してもう 一度動かすだけではだめで、 init.d/lpd start すると使えるようになる。な んだかねえ。
• Slackware の機械をリモートからいじりまわして電源いれたらいきなり なんか走るようにしてみる。こちらは rc.local にちょこちょこ書くだけなの で簡単。 root が見る .xinitrc が /root/.xinitrc ではなく /.xinitrc な のはそういうものなのか？
• お、ちゃんと更新情報がでてる。
• 今年の Gordon Bell の選考委員会委員長は Thomas Sterling か、、、

• 今日からプリンストン。うーん、やはり DHCPD に嫌われている。 Vine なんか嫌いだ、、、隣の Fiva+Slack8 は動いてるし、、、(現在の設定は、多 分 DHCP がくれるはずのアドレスにこちらで手で設定している)
• 銀行にいく。こういう時に共同研究者に全面的にお世話になるようでは いかんような気もするが、まあねえ。
• IAU のアブストラクトは送った(富坂様すみません。忘れてました)。次 は校正刷を直して、日本人論文を、、、ええと、 61 ページで図が 40枚？で、 オランダ人論文は、、、 double column で 21 ページで図が 20 枚。これ読 めってか？
• 中里君とこで更新情報が拾えてないみたいだけど、 ここ はちゃんと更新されていたのね。
• そうか、元ボス氏もサッカー気にしたりするのか。US はサッカーとはほ とんど無関係な世界であるような気がする。一昨日マンハッタンをうろうろし ていたら、ニュースの電光掲示板で

  USA eliminated from the World Cup

  とかいうのがでてて、 eliminate というのはこういうときに使うのかと学習 したくらい。ホテルで ESPN とかみてもほとんど野球だし。
• 6/24-28 New Horizons in Globular Cluster Astronomy ?
  すみません。これはさぼってます。

• 6/中? インターフェース〆切
  6/20 Computer Today 〆切
  送った。

• 日曜からニューヨークにきている。 AMNH の LAN はへろへろなので、 結局今は Starbucks で書いてたりして、、、 2時間 20$ は高いか安いか？ とはいえ WLan card さえあればつながるのは結構というか素晴らしく便利。
• 6/17-21 Integrating Stellar Evolution and Stellar Dynamics?
  終わった。ふう。意外にうまくいったような。私としては Open source な population systhesis program とか binary evolution model とか stellar evolution program とかがあるとありがたいわけだが、意外にそういうのを作 りたいという人もいるみたい。

• 文部科学省って Ministry of Education, Culture, Sports, Science and Technology だったんだということを改めて納得するような ページ。
• 事務処理というか、人に書類を押しつける。すみません＞関係者各位。

• これは経産省ので文科省の 700億とは別なんだな。
• 致命的な間違いのうちの 2 つをあげただけなので、それだけしか直してくれないようでは困るのでございます。

• うーむ、なんというか、人というか人の学生と議論するのは難しい。と いうか、、、
• つーか、前に来たメイルには isothermal と書いてあったべさ。それが いつのまに Hernquist に、、、あ、そうか、そうではなくて実は isothermal という可能性は、、、さすがにないと思うが、、、

• ふむふむ、なんとそうであろう(日本語おかしい。気が狂ったかも)。つーか、本当はシミュレーション結果の選択を「自分たちの観測の解釈に都合がいいかどうか」だけではなくて「中身がまともかどうか」という基準も少しでいいから考えてやって欲しいんですが。しかしどうでもいいが著者の面子がすごいというか、なんかあらゆる人が入ってるような。
• しかし、こないだ van der Marel から「手頃なポスドクおらんかね」というのが回ってきた時にも思ったけど、やっぱりもったいなかったよなあ、、、この業界は世界的に人材不足ではある。銀河形成と宇宙論だけが天文学ってわけじゃないのよ。(まあ、じゃあ日本にポストはあるのか？っていわれたらないんだけどさ)。
• しかし、まあ、すでにいまさら文句をいう問題ではないのかもしれないけど、なぜただのテキストを Word 文書で送ってくるのかな？まあ、もちろん答は「Word で書いたから」ということなんだろうけど。あと、 HancomWord5 で読めるのと読めないのの違いはなんだろう？といっても Word の場合バージョン違いとかかな。
• うーん、軍事優先ということ？なんというか、、、

• ひたすら原稿書き。
• 先月分へのリンクを忘れてた。追加。
• 10個再び、みたいな。
• 27T か、、、これがでてくるのかな？

• 3だと思いますが、、、違うのかな？ちなみにこっちはもちろん違うわけですが、説明は違うような。
• 6/10 教務委員会 16- 205
  とりあえず毎週ではなくなったかな？
• うーん、タイムステップに上限を付けるとうまくいきますか、、、これはまあそうすればいいんだけど、破綻しないかな？
• 次の原稿処理開始。
• え、エリック・コタニって近藤先生だったの？それは知らなかった。

• どうも原稿が進まない。何故進まないかというと、理由ははっきりしていて読者層を設定できていないからである。何の原稿かというと「シミュレーション」に関する原稿なんだそうで、まあ常微分方程式の初期値問題の数値積分についてなんか書けばいいんだけど、これが実は難しい。

  なぜ難しいかというと、どうもいろんな本とかレビュー論文とかマニュアルとか的なものから推測すると、割合いろんな分野で常微分方程式の数値解法については迷信とか妄想とかがいろいろ存在しているからである。

  迷信が存在しているのが何故問題であるかというと、まあ、それが本当に迷信である場合の話だけど、つまり、その迷信を信じている人はそれが正しいかどうかをまともに考えたことは多分ないものと思われるわけで、そういう人に何をどういうべきかが良くわからないから。

  というか、もうちょっと考えてみると、仮にも科学者であろう人の間に迷信が蔓延するのはどうしてか？ということが本来問題である。楽観的な解釈としては、「そこんとこの知識が間違っていても大した問題ではないところに迷信がはびこる」というのが考えられる。この立場をとるとすれば、世界に特に問題というべきものはないはずである。以下、実例で考えてみよう。

  ここに「粉体シミュレーション入門」という本がある(今実は手元にないので書名はあやしい)。このなかで離散要素法が紹介されている。これは要するに砂山とかそんなのの動きを、砂つぶ1つ1つの運動を時間積分することでシミュレーションしようという方法である。粒子間相互作用は、もっとも簡単には接触した、というかめりこんだ時にめりこんだ深さに比例した反発力と速度に比例した抵抗が働くとする。

  抵抗の大きさは、問題によるが例えば臨界減衰になるように決める。 問題は時間積分をどうするかということだが、紹介されているのは速度を前進オイラー、位置には2次のアダムス・モールトン修正子を使う方法である。時間刻みについては、Δt<2k^-1/2になるようにとれと書いてある。ここで kは調和振動の項の係数であり、この時間刻みは回転角が2ラジアンになる場合に相当する。 参考文献もある。

  この回転角はやたら大きいと思うかもしれないが、リープフロッグの周期性区間の上端は2なので、もしも使っている積分公式がリープフロッグであり、 速度に比例した抵抗もない系で、さらに計算精度が問題でないのなら間違っていない(リファレンスはそういう場合について議論した文献であった)。が、上の条件はどれもあてはまっていないので、これではだめである。具体的に臨界減衰の時に安定性を評価すれば、真の解とはほど遠い解がでてくることがわかる。上のは安定性限界になってはいて、もしも本当にこの値に時間刻みをとると差分方程式の固有値は -1 になってる。もちろん、これより小さければ収束はする。臨界減衰でなければ安定性限界はもっと小さいのはもちろんではある。

  しかし、まあ、だからといってこれで計算した答が間違っているかというと必ずしもそうとは限らない。というのは、もともと臨界減衰にするかどうかには自由度があるからである。適当な散逸が入っていさえすれば詳細はあまり問題ではないような系を扱っているのであれば、モデル自体が多少変わってもかまわない。この場合には、数値解がモデルの真の解からかなりずれていても問題ではない「かもしれない」。

  が、問題でないかどうかは本当は調べてみないとわからないわけで、だから自分はそういうことをしているとわかってやる分には問題ないけど、気が付いてないといかんということはあるであろう。実態はどうなんでしょうね。

  というわけでなんの考察もなく終わる。付録:ありがちな迷信

  1. 低次の積分公式のほうが安定性が高い。
     線型の問題に関する限り全然そういうことはない(陽解法の場合)。オイラー法よりも2次のルンゲクッタや ABM (PECでの)のほうが安定領域が広い。ルンゲクッタの場合は4次までは次数が高いほど安定領域が広い。ABM はさすがに3次以上は狭くなるけど。 もちろん、加速度が不連続に変化する場合は話が別で、特に多段階法は BDF のような特殊なものを使っていなければ不安定になる。
  2. ルンゲクッタのほうが線型多段階法より精度がよい。
     同じ次数で non-stiff な問題ならば線型多段階法のほうがつねによい。(シンプレクティック法は議論が別)普通の陽的ルンゲクッタの利点はプログラムが簡単なことだけである。(これはとても大きな利点だけど)
  3. ハミルトン系の長時間積分には線型多段階法よりリープフロッグのほうがよい。
     上の逆ですね。これは要求する計算精度に依存するので簡単には言えない。
  4. Multivalue method は multistep method よりよい。
     大抵の場合、「同じ」。数学的には対応する(丸めを別にすれば同じ答をだす公式が存在する)し、丸め誤差についても Multivalue がよいということはない。時間刻み可変でなければ multistep のほうがプログラムが簡単。
  5. 高い精度の結果を短い時間で得るには線型多段階法より外挿法である。
     よほど高い精度でない限りそういうことはまずない。
• というわけで原稿のほうは適当にけりをつける。上のような話をきちんと議論すべきかもしれないけど、ページ数的に無理。大体そんなのよみたい人がいるかってのが問題。
• でもって次の原稿。何書くんだったかな？

• つーか、偏平率は任意なのか。
• 騒音は何 dB か。
• そうだよね＞岡本様(タタリ氏とは違うほう)。営業時間とか定休日とかいった制約があるにしても、ちょっと考えられない。
• 1年ちょっとまえ、某研究会の SOC Chair と LOC Chair の会話。参加者名簿みながら。

  S: この Takashi O.(全然仮名になってない) って誰だっけ？こんな人京都にいた？
  L: あ、それH大の H さんとこにいた O 君ですよ。
  S: え、彼はだってたた、、、(と言いかけてやっと気が付く)、あ、そか、そうだよね。

  私のようなアホウはそうそういないと思うけど、まあ、そういうことも世の中にはあります。

• 原稿+原稿用プログラム。
• Maximal rotation のマクローリン楕円体ってどんな格好だっけ？

• ひたすら原稿+原稿用プログラム書き。

• 記憶がない。なにしてたっけ？

• 13440058
• なんか会議が増殖、、、
• でもって会議の後物品検査に対応して作文して作文したら談話会の時間、、、なんか原稿書きすらする時間が、、、

• ノートにいれた Vine 2.5 が日本語は使いやすいので、同じような環境をデスクトップにもいれたくなるわけだが、 Vine だとどうするのかな？ grub だし、 fstab がよくわからない代物だし、Slackware のようなわけにはいかなそう。 RH はなんかはいったんだけど、これは LILO だしねえ。
• GCC 3.1がでてる。ついに SSE2 サポート(といっても SIMD にはならないんだそうな)ですか。Northwood のなにか買ってみるか。
• 某所で見つけたGA による TSPの解なんだそうだけど、これどうみても間違ってるよね。
• 6/3 15- 会議
  大先生というのは固有名詞なんである。

• 75674
• quota の使い方を勉強する。
  1. kernel config で quota support を on にする。
  2. fstab のオプションに usrquota を付ける。こんなふう。

     /dev/hdx1       /bar        ext3        defaults,usrquota   1   3
     	      

  3. rc.local あたりに quotacheck, quotaon を -auvg でつけて置く。
  4. reboot!
  5. edquota username で個人設定をする。こんなふう?

     	Disk quotas for user foo (uid xxx):
       Filesystem                   blocks       soft       hard     inodes     soft     hard
       /dev/hdx1                  332852   300000   350000      4285   10000   20000
     	      

     単位は KB らしい。これだと 350 MB が上限のはず。 inodes はファイルの数、 soft/hard は 見ての通り。
• とりあえずリンク だけ。なんとも、、、

• 牧野のホームに戻る