牧野の公開用日誌

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

Copyright 1999-2002. Jun Makino

当面の予定

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

2002/10/31

• というわけで「特定計画」サポート。どうもありがとうございます＞関 係者の皆様。
• 早口言葉。特定領域研究計画研究研究計画調書。

• 今日は駒場で打ち合わせであった。
• N=60k, 4 node で 10,000 particles/sec ですか。それならまあ普通かな。16 ノードの数字は astro-ph/0003162 の図11の右のほうから持ってきてみました。
• とはいえ、意外に弱気という噂なのであんまりいじめないほうがいいかも。
• おお、S様真に偉大。感謝感激。

• こんなものが。 Cohn 先生プンプン。
• それは風邪です。早くおうちに帰って十分休養をとりましょう。＞北のほう。
• GRAPE-6 の実効性能の資料を作成を頼まれていたのを思い出す。折角なので公開。並列化した時の加速率がいまいち。
• それはいいけど Gadget って今でもこんなに遅いの？レフェリーした時に、これ遅い、、、と思ったわけだがあらためてみるとやっぱり遅い。
• なんてことを書いてるとそのうちにこんなことがおきたりして。

• 最近気になること。このChandrasekher氏、カタカナだと「チャンドラシーカ」(例)となって ることがあるんだけど本当にそんな発音なのか？
• KT400 というか MSI KT4V/L。なんかまともに動くぞ。まあ、266A も MSI のはメモリが Memtest でこけるだけで(これはこれで問題なんだが) PCI はまともだったっけ。 あれ、そうじゃなくて DMA はおかしかったよな。さす が KT400 は直ってるのか。

• ほぼ復活。
• nForce2 マザーもでてきますか。KT400 と一緒にテストしてみるか。今年は 2700+ でいくしかなさそうだし。
• 10/27 停電
  計算機群は特に問題なく復活。 GbE スイッチ交換もトラブルなし。4ノード以上使っても問題起きなくなったような気がする。どうもです＞稲葉様、玉造様、平木様。交換作業どうもです＞台坂様
• 最近あんまりどこも壊れないな。とか書くとディスクが死んだりするかな。
• とても赤新聞(英語では yellow journalism かもしれないけど、ページ自体が赤いんだもの)とは思えない格調高い記事。

• 風邪でねこむ。
• 電源落としに夕方大学に。

• 風邪をひく。
• 別に、難しい原稿を読んだので知恵熱がでたとか、だから晩飯くらいおごれとか、そういうことは(思ってないとはいわないが)いってないからね＞原稿書いた人。
• 昨日本屋で見つけた「知の歴史」、大学に送られてきてた。私も少し(本当にほんの少しだけ)お手伝いした本なわけだが、何故かちゃんと名前がのってたりする。しかし作業が速いね。問い合わせがきたのが記録によると9月中旬なんだが。
• 10/25 16:00 矢沢さん
  McMaster が GRAPE-6 導入?本当に? あ、Alison Sills がいるのか、一瞬 Couchman かと思ったがさすがにそんなはずはないな。
• 10/25 特定〆切？
  これは出してないので関係ない。
• 何故人は同じ間違いを何度でも繰り返すのか？うーん、何故？

• 難しい原稿を読む。

• おたんじょうびおめでとう。
• Opteron 1万台 ですか。あれ、 ASCI Q は出来たんだっけ？ Opteron の SSE2 はどれく らいまともなのかなあ？しかし、日本語訳でまで ES が 35.9 TF になってる のはちょっとねえ。ピークと実測をくらべて議論しないでほしい。まあ、2004 年になって Opteron が 2GHz でしか動かないということもないような気もす るけど。でも AMD ならありえるかも。
• しかし、1万台で1億ドルか。クレイ(つーか、 Tera だな)はボロ儲けでないか。
• あ、そういえば SC に申し込みしておかないといかん。
• というわけでせこせこレジストレーションとホテル予約。ホテルは既に conference rate はありませんという冷たい表示がでたので、ホテルリストのなかで適当なのを google すると こんなものがでてきてそこで予約。これ大丈夫かしら？
• SCI2003 来ましたね。私のところにはまた3通来ましたです。ったく。
• 道具ならあなたのほうが、、、
• multi-mass King model...... うーみゅ。世の中適当である。
• 特定なんか出してもなかなか当たらんと思うけど、出さないと絶対当たらないしねえ。
• 展示用ボードをでっち上げる。何故か SC2002 で展示されるといふ。どもです＞関係者の皆様。

• 10/22 基盤〆切
  出した。
• 誰もこれに言及しないのは、みんな5分で解けたからでは多分なくて、みんな解けていないからだろうか？(私は解けてないことはいうまでもない)。実はみんな5分でできてたりするとちょっといや。永久機関のほうはまあ計算だけだと思うけど。

• これはさすがに疑いようがないよな。
• 10/21 10:00 3-404 打ち合わせ
  無事終了。
• だめすか？そかも。タイムステップ、というかタイムスケールの分布によると思うけど、、、大半のガスが集まってしまうならそれはそれでそこを分割しちゃえばいいわけなので。

  ちなみに三鷹の AMR は、階層毎に違う領域分割をするというなかなか素敵な方法でした。並列化効率については、、、えーと、、、時間刻み一定の時はよかったと思う。

• 今日の Astro-ph。第3積分。
• 田崎さんのOHP の書き方。重要なことがまとめて書いてある。

• 領域分割と独立時間刻みは原理的には整合性が悪い、特に single halo でうまくいくような気がしないというのはその通りなんだけど、まあ、問題は single halo を何台くらいで分割するかと、 single halo のなかに substructure がどれくらいあるかではないかと。 本当にのっぺりした NFW プロファイルみたいなのが一つだけあると、これを細かくわけてしまうと多分 どうしようもなくて、 8 台とかそれくらいがせいぜい。 まあ、このときは独 立時間刻みがどれくらい必要か？という問題はあるわけで、最近の Powers et al でも独立時間刻みの効能は意外にうすいという結果にはなってたような。

  やはり、2nd order で独立時間刻みにすると異常に短い時間刻みにしないと積 分精度がたもてないので、、、例えば、 mass resolution が 10³ くらいになったとしても、 DM だけ考えると dynamical time は 1Myr とかそ んなものなので時間刻みは(一定刻みなら) 0.1Myr くらいで多分 OK。ところ が可変刻みにすると 10⁴軌道積分してもエネルギーに精度がのこっ てないといけないので、 dynamical time あたり数百ステップ必要になるので。 4 次の公式とか使っても 50 くらいはいるんじゃないかなあ(このへん適当。)

  とはいえ、流体入れるとそれではすまないわけだけど、この時はどうせそこら じゅうで高密度領域ができるわけで、領域分割でもいいというわけにはいかな いのかしら?技術的にはむしろ、領域分割して tree を作ったとして、短いス テップで一部の粒子がアップデートされた時にそれを他のプロセッサの tree にどうやって反映させるかと、反映させるべきかどうかをどうやって判定する かが問題な気が。

  まあ、それ以前に GRAPE 使うと single node でも独立時間刻みで性能がでる ような方法は(現状のハードでは)知られていないという問題があるのですが。 プロセッサモジュール作りなおして(チップはいじらない)なんとかするという のは昔考えたんだけど、、、人手が、、、

• 講義準備。去年は初回に破綻したので、今年はそうならないように、、、
• で、科研費作文。〆切が明後日であった。間に合うのか？

• あ、なるほど、独立時間刻みの話なのね。Gadget の方法は原理的に scalable ではないので、並列化するならなんか違う方法を考えないと幸せになれないような 気もするんですが、、、 GRAPE 使うとどうしても相対的に通信が遅いし。
• あれ、下の、 and と or が違う。一応直しておく。

• 10/18 放送大学10:00-
  杉本先生、ロケ隊の皆様、お疲れさまでした。しかしカメラの前で話すのって なんかいつまでたっても苦手というかうまくできない、、、確実に練習不足で はあるんだけど、、、
• 0210158 が改訂されてて驚く。諦めが悪いというだけでは表現しきれま せんね。
• あれ、同学年じゃなかったっけか。なお、下のリンク先では写真より絵 のほうがもっと似てるという話も。(ネタにしてすみません＞本人)
• あああ、人の論文いじってるばっかりでなくて自分のもかかねば。(ちょっと自己主張すると 0210133 の first draft は私がでっち上げたんだけど、計算したのは Holger だし仕上げたのはヨーロッパ人連合だから、、、)
• 私は並列版は sort から insert 方式に切換えたです。だいぶ遅いのがいやなんだけど。 SSE(2) でも使うように書き変えるかな。
• おせっかい。

  unsigned int high = ((1<<(N-p2-1))-1)<<(p2-1);
  unsigned int mid  = ((1<<(p2-p1-1))-1)<<p1;
  unsigned int low  = (1<<p1)-1;
  unsigned int i = (1<<N)-1;
  
  unsigned int il = i & low;
  unsigned int im = i & mid;
  unsigned int ih = i & high;
  im = im << 1;
  ih = ih << 2;
  unsigned int ix = il | im | ih;
  	  

  がなにをしたいか？ということを考える。上のを文字通りに取ると、

  unsigned int ix = ((((((((1<<(N-p2-1))-1)<<2) + 1)<<(p2-p1-1)-1))<<2)+1)<<p1) -1;
  	  

  (括弧があってるような気がしないが、、、)ということがしたいように見えるし、それなら

  unsigned int i = ((1<<(N+2))-1) ^ (1<<(p2+1)) ^ (1<<p1);
  	  

  でよさげにみえるが、多分したいのはそういうことではなくて、

  #define	MASK(n) ((1<<n) -1)
  #define INSERT0(i,n)  (((i & (~MASK(n)))<<1)| (i&MASK(n)))
  i = ...
  unsigned int ix = INSERT0(INSERT0(i,p2),p1);
  	  

  といったことなのかもしれない。なお、良い子はマクロは使わないで inline 関数を使いましょう。

• ゼミの時間を間違える。 10:30 からじゃなかったっけ？まあ、メイルを見ると確かに 10:00 からと書いてはあるわね。
• 10/17 駒場 打ち合わせ 13:30-
  無事済む。
• Netmeeting でなんとかという花輪さんの野望はもろくも崩れ去る。まあ、 MS の作るものなんてそんなものであろう(偏見)。
• 読まないといけないもの一つ処理。依然処理キューは伸びるばかりではある。つーか、本一冊が、、、

• 10/16 10:30 NHK 高間さん
  科学者版プロジェクト×、、、、あ、ちょっと字が違った、、、
• 昨日はケンブリッジは一日中雨だったそうな。
• なんで僕だけ会議に出てるの？爆笑。
• いかん、私の処理能力をはるかに超える速度でみないといけない論文が発生している、、、
• おーい、、、というかこれ。

• ノーベル賞のおかげで公開講演会では話さないで済むかも。ありがたい。
• 375 日雨が降るというのはどこの国？私がいた 3週間の間、ケンブリッ ジは(夜にちょっと降った他)まったく雨が降らなかった。
• みー、陰謀は既に露見、というか露見してもなくならないのか。うみゅ。

• 時差ぼけのまま。明日も休みで助かった気分。もともとスケジュールはそのことを考慮しないで組んであったので。

• というわけで無事日本に帰ってきました。帰りの便 (JL402) はとっても すいていた。
• で、3席確保して寝てたわけですが、途中でぼーっとみていた映画のなかで、ピンポンが凄かった。いや、話がどうとかいうわけではなくて、「スマイル」役の ARATA、某君の生き別れのお兄さんかと思った。

• あ、復活してる。S 様感謝。というか大量のメイルが、、、
• あー、おっさんになってしまいましたか。まあ、あの発言ではしょうがないか。 えーと、実は学部の卒業研究の指導教官なので、、、
• 今晩の飛行機で日本に戻ります。
• しかし次々にいろんなことが起こるなあ、、、やはりこれはなにかの陰 謀に違いない。
• Fowler は、、、 "for his theoretical and experimental studies of the nuclear reactions of importance in the formation of the chemical elements in the universe" となってますね。そういう意味ではこれは理論家と実験家の共同受賞なの？

• あー、緊急のメイルも電話もない。平和だ、、、
• レポート書いた。ちょっと寝かしておこうかな。
• あれ、メイルサーバーが、、、死んではいないけど sendmail が御返事 くれない。 S 様見てます？
• 確かに、この 発言はかなり不快。島津の研究者の成果が科学技術政策の成果の証明だと か

  田中さんに『こんな実験をしてみたい』と独創的な提案をした研究者が、だれかいたのではないか

  なんてことを本当にいったのだとしたら、、、うーむ。
• なるほど。まあ、天文、理論というと今のところ Chandrasekhar だけだっけ？ Hoyle は、、、とかいう話はもちろんあったわ けだ。
• 今から天文でノーベル賞をと思ったらやっぱり重力波なんだろうね。 LIGO より後でもいいから space で低周波を LISA より先になんとかするというのがお買い得。

• というわけで昨日の晩は high table dinner というものでございました。 ガウンをはおった Aarseth と Chris。なんというか、極めて謎なものであっ た。私の斜め向かいに座っていたのは Bondi であったという話を後から聞い たような。 隣に座ってたのは現在の College の Master であるところのも と駐日大使の奥さんとのことでした。ヨシカワがどうとかアリマがこうとかい われても私は面識ありません。
• トーク 終わった! かなり受 け狙いでいったわけだが、それなりになんとか。最初の質問がなんであるかの 予測が的中した。私には珍しいことである。
• direct-x? shockwave? なにそれ？食べられるの？みたいな気分になる。
• こないだ A&A をなんとかしたと思ったら次は ApJ ですか、ちょっと多 すぎないか？年間6本以下とかそういう制限ってしちゃいけないのかしら。あ、 ラングミュアだ。おお、それに某氏のもひかれている。そっちにまわそうかな あ。
• とりあえず騒ぎは先送りになった。先送りになっただけというあたりが 激しく問題なんだけど。

• というわけで、 M15 にブラッ クホールなんかあるはずがないという我々の論文とそれに対応して出たGerssen, van der Marel その他による Addendum なるもの。諦めが悪いなあ。1σで BH が あるほうがいいっていったってねえ。
• 宣伝:この結果については、明日の R71a (A会場 14:42) で発表がありま す。
• そういえば、日本語でも「平服」っていったらスーツにネクタイだっけか。
• GRAPE-6 のホストに icc, ifc, gcc-3.2, MPICH 1.2.4 をいれてみる。 別に GbE で性能がでないのが直ったりはしていない。
• あ、なるほど、Bahcall 先生も可能性はあったのね。Bahcall といえば、彼が初めて太陽ニュート リノの計算をして、なんだか妙な結果がでたので悩んでいたら、 Feynmann に 「妙な結果だって、それはとりあえず論文にしちゃえ」というどこかの国のだれかの ボスのようなことをいわれて、本当に論文にしちゃったというような話がのっ てたのは Collins and Pinch の The Golem だったかな。
• ひえー、急速な展開が、、、

• Eggleton を目撃。
• 科研費は、、、特定をどうするか決めておかないといろいろ面倒。うみゅ。
• あれ、公募と計画代表は重複申請できるのか？なんか日本語が曖昧だ。
• 明日のディナーはジャケットとネクタイでね、とホストからお達しがくる。 ケンブリッジに 40 年いる人からは「ネクタイはいらないから」みたいなメイルがきていたわけだが、、、(9/21 の項参照のこと)。 "In Cambridge, informal means tie, it is one of secrets of Cambridge" とはケンブリッジにきて 1 年くらいのオランダ人の お言葉。うーん、奥が深い。
• ちなみに、40年在住氏によれば、 Trinity College とかならともかくChurchill College にそんなものはいらんだろうと思ったんだとのこと。ちなみに上の写真をとってる位置の右後方に IoA がある。
• しかし、 M/L が間違ってるのはともかく、そもそも M/L が一定の解析しかしてないもう片方はどうしたものかしら。
• 某所にあった怪しい積分法の論文を眺める。こんなのを通すレフェリーは阿呆である。

• お、反応が。え、使った M/L のプロットが間違ってた？どれどれ、、、 あ、本当だ。なんか倍くらい違う。嘘のグラフを論文にのせた人とそれを信用 して使った人のどっちが悪いかってのは難しいけど、でも、その論文の他のグ ラフとつじつまがあってないんだから、どちらもかなりはずかしいかも。
• Rees は土曜も日曜もきている。
• あ、こないだ遊びにきてた Aarseth 先生のお友達の Anna ってのは、そ うか、 TZ object の Zytkow なのか。記念写真とってもらうのであった。それはともかく、それで彼が TZ object を発音する時のアクセントが思いっきり Z のほうにあるのね。これが TZ であって ZT でないのには深いわけがあるというのをなんかで読んだ気がするんだけど、思い出せない。
• Aarseth 先生の本の原稿を拝見させていただく。なんかもうすぐにも出版できそう。 Heggie のいうところの Bible before its publication である。 Old testament か New かとか、そういう話はあるんだが。
• 某 X 線の偉い先生からのメイルも狂った astro-ph がひかれてたので、 京都のせいではないらしい。明日のはちゃんと配送されるのかしら。
• City Centre にいって買物。で、戻ってきたところから email 爆発 が、、、3時間に 70 通くらいが大西洋を飛び交う。まあ、火曜日の astro-ph で一部が明らかになるわけだ。ふう。

• 昨日 Ostriker にそっくりなおじさんが危なっかしげに自転車にのって いるところをみたんだが、あれはやはり本人だったのだろうか？昨日研究所に 顔を出したという話だし。
• 自動販売機でポテトチップを買っているところに Thorsten Naab から声 をかけられる。彼もここのポスドクか。
• 私が寝ているうちに極東と東海岸の間で往復した膨大なメイルをフォロー。 Jeans equation がなんだかわからなくなったので自分の講義ノートをみるとちゃんと書いてある。講義した時の私は今よりはるかに賢かったに違いない。BT を丸写ししたのではという気もちょっとするが。
• あれ、 astro-ph が発狂した？

  Subject: astro-ph daily 0209346 -- 0209379 received 1351
  From: e-print distributor 
  To: astro-ph@yukawa.kyoto-u.ac.jp
  Date: Sun, 6 Oct 2002 00:50:25 +0900 (JST)
  

• Jeans equation が良くわかってないのは私だけではないということを確 認して安心する。
• よーし、投稿したぞ。後は野となれ山となれ、、、あ、いや、そうじゃ なくて、これからが問題なんだけど、まあ、politics は外人任せ。記録に よ ると 9/25 の夕方から 9/27 の午前中までで第 0 次稿が上がってるから、そ れから 10日たってない。なかなか早く投稿できたといえよう。共著者の皆様、 特に Holger, thank you for your hard work in the last 8 days!

• メモ tennet 2520
• 図面どうもです＞久賀様、鈴木様。でも、確かにちょっと意味が薄かっ たかも。肝心の容量がどこにも書いてない。コンセントなんかついてればいいっ てばさ。ないと困るけど。
• だぼらのつもりだったけど、 そうでもない。 まあ、 local group だから星が分解できるんだよね。
• 論文の調整で 100 近いメイルが3大陸、、、でもないか、島ばっかりだったりするけど、、、アムステルダム大学はちょっと前まで海の中だったところにたってるし、、、を飛び交う。なんか疲れた。まあ、明日か明後日には投稿できるかな？

• Cruithne 。あ、 Mikkola が Nature に論文のせたやつのことか。そういえば 彼は Symplectic Integrator の本を書いてるとか。谷川さんではない日本人 と共著という話だったけど、吉田さんなのかしら。
• とはいえ、 Seppo といえばアルコール、世界中のだれから聞く話でもそ うなのがすごい。フィンランドはそーゆー国だという話は Linus Torvalds 自伝にもあったので、まあ、そーゆーものなのかもしれないけど。
• SSE を触ってみる。確かにそれでできることなら x87 の 3-4倍速い。
• ifc について日本語でもっとも詳しい情報があるのはここかな。
• 朝の5時(現地時刻)からメイルが来る。いつ寝てるんだ？これから寝るのか？
• あれ、 SSE だと 3 倍でるのに、 SSE2 だと x87 と変わらないの？なんかどっか壊れてるのかしら？
• ひえー 。ちなみに、解析的な結果は最初のアイデアも含めて私ではないほうの著 者からでてるのでございます。
• しかし、天文で必要な数学ってとても ラフマニノフまでいってないような気が。「ガラパゴス」とかいう話も昔 あったし。
• 吉田さんかどうかは未確認なんだけど、今は Saha と書いてるとのこと でありました。 しかし Acting Director? 古在さんが台長みたいなものかしら。

• あら、カウンタの桁が変わってる。
• アイルランド人はアメリカに帰っていたらしい。ので昨日の記述はちょっ と嘘。
• (主に先週末に)とった写真 。
• ふーん。

  千葉大学先進科学教育センター教官　　
  　　　　　　　　　　　　　　　　　　　
  募集人員　　　教授または助教授 1名
  所属　　　　　先進科学教育センター
  専門分野　　　宇宙物理学（理論シミュレーション）
  

• 三鷹用に GRAPE-6 ライブラリをいじりまわす。これで X(仮名) も Y(仮名) も怖くない。
• 来週のセミナーは Blackholes in star clusters? --- M82 and M15. でいきます。 "I like this question mark." とはタイトルをみた Aarseth 先生のお言葉。
• 極めて重要なメモ。 ifc -Vaxlib で Vax Fortran 化する。さすがにも とは DEC Fortran だっただけのことはある。
• それはいいんだけど dvec.h を読み込むとコンパイルエラーがでるんだ が、、、

• g77 (3.2) で私の書いたプログラムは確かに遅いんだが、Aarseth 先 生のは速い、、、何故だ？アセンブラみても違いがわからん。
• あ、新しい図がきた。なかなか素晴らしいではないの。現在の瞬間だと日本にいるドイツ人、イギリスにいる日本人、オランダにいるアイルランド人、アメリカにいるオランダ人とオランダにいるオランダ人で論文書いているのである。
• 99611
• 723 は7号館2階ですね。2階は階段教室はなかったような気がするから中 くらいの部屋？駒場の教室番号ってなんか変な風についてたような。今資料が ないので記憶で書くと、1号館は1階が 101 から、2階が 151 から、5号館は 551-(1F), 571-(2F), 591-(3F), 7号館は 701-, 721-, 741-, 761, 11号館は 階とは無関係に1101 からの通し番号、 13 号館は 1311 が1階,1321 が2階に あったような(これちょっと怪しい)。 12号館にも教室があるような気もする けど使ったことがない。
• なんだか次々と原稿が送られてくる。
• いんちきでもなんでも、8日に1本はすごい。 refereed journal で年に40本というのは、論文の中身を自分でみてたら不可能だと思うが、、、半分くらいしかみてないとすればできなくはないかな。
• セミナーは GAIA の分光器。 R=10000 くらいの分光器をつむんだそうで。 JASMINE はそんなのないよね。
• "So, this is a very provocative paper." そだよん。
• Intel Compiler はどうも振舞いが分からん。 P4 で -tpp7 -xW とする とかえって遅くなったりするし。昔のベクトル化コンパイラみたく、ベクトル 化したら遅くなるのでもめったやたらベクトル化するせいかな。でも、 -tpp6 より -tpp7 のが遅いのはなぜだ？
• メモ

• 牧野のホームに戻る