牧野の公開用日誌

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

2002/7 2002/6 2002/5 2002/4 2002/3 2002/2 2002/1
もっと昔

当面の予定

9/2-9/28 Gerhard Lienhart
9/3 広報委員会 15:00-17:00 1-205
9/4 3号館理研 15:00
9/4 会議 13:00-
9/5 天文大学院2次 10:00-
9/6 浜松 10:30-
9/4 会議 13:00-
9/13 会議 10:00-12:00
9/22 3号館停電 8-12
9/23-10/11 ケンブリッジの予定

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

2002/8/30

8/28-9/21 Piet Hut
bmobile なカードを使うというので RH な機械に入れる。 Ethertap が入ってるらしいので ppxp を入れるが動かない。で、 pppd を動かしてもうまくいかなかったんだけど、 RH の pppd setup tool を動かしたあと自前の ppp-on を動かすと何故か使えるようになる。 RH とか Vine とかってなんかよくわからない。カーネルも /usr/src/linux で make したら動くものができるというわけではないみたいだし。
8/30 会議 13:30-17:00
意外に早く終わった。

2002/8/29

8/29 会議 10:00-
意外に早く終わった。が、疲れた。明日も会議だし。

2002/8/28

そうか、「自分で見易い」ためなのね。なるほど。

      subroutine symchC( m,         C,      tp    )
      integer       i,j, m; real(8) C(0:*), tp(0:*), T, S; S = 0.0d0
      do i = 0, m-1
      do j = 0, m-1 ; tp(i*m+j) =( C(i*m+j) + C(j*m+i) )/2
                   if( S .lt. ABS( C(i*m+j) - C(j*m+i) ) )
     &                 S   =  ABS( C(i*m+j) - C(j*m+i) )
      end do; end do
c
      call maxel( m*m, C, T); C(0:m*m-1)= tp(0:m*m-1)
c
      if( S/T .gt. 1.0E-12 ) write(6,*) 'sym C = ', S/T, m
      end

とか書いてあると、思わず

      subroutine symchC(m,C,tp)
      integer i,j, m
      real(8) C(m,m),tp(m,m)  T, S 
      call maxel( m*m, C, T)
      S = 0.0d0
      do i = 1, m
         do j = 1, m
            tp(i,j) =(C(i,j) + C(j,i) )*0.5D0
            S =MAX(S,ABS(C(i,j)-C(j,i)))
         end do
      end do
c
      C(1:m*m)= tp(1:m*m)
c
      if( S/T .gt. 1.0E-12 ) write(6,*) 'sym C = ', S/T, m
      end

とか、さらには

       
      subroutine symchC(m,C)
      integer i,j, m
      real(8) C(m,m),tp T, S 
      call maxel( m*m, C, T)
      S = 0.0d0
      do i = 1, m-1
         do j = i+1, m
            S =MAX(S,ABS(C(i,j)-C(j,i)))
            tp =(C(i,j) + C(j,i) )*0.5D0
            C(i,j)=tp
            C(j,i)=tp
         end do
      end do
c
      if( S/T .gt. 1.0E-12 ) write(6,*) 'sym C = ', S/T, m
      end

とかしたくなるけど、まあ趣味の問題ではありましょう。(下の2つがあってるかどうかは無保証)

/2 を *0.5D0 と書きたくなるあたり老化現象かも。
一番下が「分かりやすい」と思うのは既に脳が溶けているからではないか？と言われると否定しがたいものがあるし。

            tp =(C(i,j) + C(j,i) )*0.5D0
            C(i,j)=tp
            C(j,i)=tp

は

            C(i,j) = (C(i,j) + C(j,i))*0.5D0
            C(j,i) = C(i,j)

のほうがいいかな。

2002/8/27

今日きたメイルにこんなのが
Subject: ３階の足跡について
理学部もなかなか楽しいところである。
ありゃ、会議の日を間違えてた。明日ではなくて明後日ではないか。
論文一本みて、レフェリーレポート書いて、それからちょっとゴミ捨てを、、、と、これくらいは今日中に片付けないといかん。
論文処理すんだ。
富士通のニュースリリース、良く読むと結構素敵。
- 最大構成でのピーク性能は VPP-5000 の 17 倍
- 最大構成で実アプリケーション性能は同 7-8 倍
しかも、 128 ノードで 128 本独立な計算を流した時の性能比較だったりして。まあ、それだけ出れば偉いという考え方もなくはないが。
レポート書く(どうでもいいけど天文関係じゃないからね)。「却下」の2文字ですましてもいいような気もするけど 70行ばかりだらだら書く。良い子のみなさんは
- 図の軸には名前と目盛をつけましょう
- 図の曲線には説明をつけましょう
- 本文中では a, b, c, d なものが図では I, II... だったりしないようにしましょう
ゴミもだいぶ処理した。出田様感謝! とりあえず 404 は 1.5 倍くらい広くなったような。 402 は膨大なゴミを出したわりにはなんにも変わってみえないけど。
8/27 天文大学院1次

2002/8/26

8/26 会議 13:00-
大変であった。
ところで、このFortran プログラム、多次元配列を使わないのはなにか深遠な理由があるのかなあ。 Fortran の(生の C に比べた)利点の一つは、多次元配列が数値計算で使いたいようなやりかたで使えるということである。

2002/8/25

再起動。
並列ツリーコードのベンチマークデータを取ろうと思ったら、あれは MPI を普通に使っていたのであった。邪道なライブラリを使うように書き直す。ついでにライブラリを公開。
富士通のスカラ型スーパーコンピューター。最大構成で 16,384 CPU の機械が 8 CPU で1.7億、、、比例よりずっと安いとしてもなかなか。最大構成って 1000 億くらいにならないか？
K7 は2600+ が発表されたと思ったらP4 は 2.8GHz ですか。

2002/8/23

正しく休暇モード。

2002/8/22

休暇モード、、、のはずだが午後は大学に。ちょっと議論。

2002/8/21

Lunch talk って何かと思ったらこれ(20.08.02)か。まあ、そーゆー人は Haehnelt and Kauffmann でも読んでりゃいいのよ。
とはいえ来月くらいにはこれの話も新しくやった分をまとめないといかん。
宣伝不足はそうかもしれないというかそうなんだけど、宣伝するよか新しい仕事したほうが楽しいような。大体 Haehnelt はこっちの仕事はよくよく知ってるはずだから、なかなか確信犯だわね。まあ、そういうふうに仕事をしたいならそれも一つの研究者としての生き方ではありましょう。なんか大変そうな気はするけど。
8/21 会議 13:00-
大変である。

2002/8/19

防犯カメラヒアリング。去年は建物を間違えたので今年は間違えないように注意。
どうかな？必要性はそれなりに御理解いただけたと思うのですが、、、
host14 が御機嫌斜め。
一度も論文を投稿したこともないしそもそもみたこともない雑誌からレフェリー依頼が来る。何故だ?
今週は私は一応休暇モードなんですが、、、前半はちょっといろいろそうはいかない。後半は本当に休暇とります。
期間限定 Computer Today 原稿一部。なんかえらそうな文章だ、、、大体、計算「機」科学かどうか？
16 ノード占有する計算を投入する。問題は1ノードで待つと全ノードがブロックしちゃうことだな。すみません＞皆様。
どうも計算が遅いと思ったら、粒子数が一桁違ってた、、、

というわけでこういう感じの悪いことを、、、

Mon Aug 19 20:04:44 JST 2002
28741 ?        R      1:23 /usr2/makino/src/pit_planet5/pit_grape6_MPI_lx g6host
13620 ?        R      1:38 /usr2/makino/src/pit_planet5/pit_grape6_MPI_lx g6host
 6670 ?        R      4:49 /usr2/makino/src/pit_planet5/pit_grape6_MPI_lx g6host
11230 ?        R      1:29 /usr2/makino/src/pit_planet5/pit_grape6_MPI_lx g6host
 3414 pts/1    R      5:05 /usr2/makino/src/pit_planet5/pit_grape6_MPI_lx -j 4 -
 6803 ?        R      1:32 /usr2/makino/src/pit_planet5/pit_grape6_MPI_lx g6host
11615 ?        R      1:30 /usr2/makino/src/pit_planet5/pit_grape6_MPI_lx g6host
 9036 ?        R      1:38 /usr2/makino/src/pit_planet5/pit_grape6_MPI_lx g6host
32236 ?        R      4:49 /usr2/makino/src/pit_planet5/pit_grape6_MPI_lx g6host
15885 ?        R      4:52 /usr2/makino/src/pit_planet5/pit_grape6_MPI_lx g6host
25523 ?        R      4:49 /usr2/makino/src/pit_planet5/pit_grape6_MPI_lx g6host
15890 ?        R      2:15 /usr2/makino/src/pit_planet5/pit_grape6_MPI_lx g6host
25513 ?        R      4:50 /usr2/makino/src/pit_planet5/pit_grape6_MPI_lx g6host
19378 ?        R      4:53 /usr2/makino/src/pit_planet5/pit_grape6_MPI_lx g6host
13592 ?        R      4:53 /usr2/makino/src/pit_planet5/pit_grape6_MPI_lx g6host
16542 ?        R      4:50 /usr2/makino/src/pit_planet5/pit_grape6_MPI_lx g6host

2002/8/18

XFree86 でスクリーンセーバーが効かないようにする方法を調べる。
```
       xinit -- :0 -p 100000 -v -s 100000
       
```
とかするとうまくいくのかな？とりあえず 10000 を 1 にすると1分でスクリーンセーバーが動くというのは確認できたんだけど。 Xserver の man page を見てもいま一つよくわからない。

2002/8/17

通った! わーい。

2002/8/16

適当な図と適当なフィッティングフォーミュラを大量生産する。うーん、計算速度モデルはホストのキャッシュメモリの効果をいれないとまるであわない、、、

2002/8/15

某研究会の集録担当から「おまえの作った PDF ファイルは汚い。ただちになんとかしろ」みたいなメイルが来たのでなんとかしようとする。
1. 最近は dvipdfm とか pdftex とか使えば PS を通さないでいきなり pdf ファイルを吐くらしいということがわかった。Note の vine にはとりあえず pdftex はあるのでこれを使ってみる。
2. 問題は pdftex は epsf を扱えないことである。 pdf ファイルなら(jpg,pngも)インクルードする (graphicx パッケージで includegraphics を使う。 \usepackage{graphics} \includegraphics[width=5cm]{test.pdf} みたいな感じ) とのことなので ps2pdf で図のファイルを pdf にすると、、、ありゃ、 boundingbox が破壊されている。
3. epstopdf なるものを使うと boundingbox が再構成されることが判明。
4. 何故か jpeg ファイルを xv で eps にしたものを ps2pdf やら epstopdf で pdf にしたものを includegraphics で読み込むと色が破綻する。しょうがないので jpeg ファイルを指定するとなんかエラーになる。
5. png ファイルを指定すると何故かうまくいくのでそれでいいことにする。
うーむ、、、なんかいきあたりばったりですね。
UA が倒産(するかも)ですか。それはなんかすごいような。

2002/8/14

スイッチもらってくる。あれ、、、まあ、それはそのうちに。スイッチが 1 階層で済むと GAMMA が動くかもしれないという希望があるんだな。これがまともに使えると TCP/IP ソケットよりはるかに速い。とりあえず 8 ポートで実験してみるか。
再投稿されたこと確認。まあ、これは多分通るであろう。

2002/8/13

復讐は確かにそこに何かの真実があるような気が。
ちなみに、日本語で「コンセント」という場合には普通こんなふうに差し込む側と差し込まれる側の両方をさすような気がする(ここにもそう書いてあるし)んだけど、これが「情報コンセント」になると何故か差し込まれる側で、しかも部屋の壁についてる時だけそういう、、、かと思って Google でみてみると、必ずしもそうでもないというか「情報コンセント室」とかいう用例のほうが多いのね。で、用例のほとんど全てが大学である。実は特殊な文部科学省用語とか高等教育局(だっけ？こーゆーのよく知らない)用語だったりするのかしら。
長い論文はとりあえず OK を出す。
短いほうは再投稿されたはず、、、だよね？
自分の論文。だいぶ進んだ気分。
しかしさらにまた別の長いものに関する不吉なメイルが海の向こうから、、、
Computer today見本が届く。えーと、なんで私のが巻頭なの？リンク先(8/19まで限定？)にあるように最初はそういう話ではなかったはずなんだけど、、、

2002/8/12

みー、だから人には向き不向きというのがあるんだってば。
両側コンセントだが、その目的であろうということは想像されるわけだけど別にその目的にはそれは必要ではないのではないかという重大な疑問がわきあがって、、、きません？
PPxP と BMH10-Jメモ

2002/8/11

何故日曜なのにややこしいメイルが、、、(西のほうからね)

2002/8/10

お休み。

2002/8/9

うーむ、人の論文を見ているわけでもないのに自分のが進まない。

2002/8/8

やはり本郷にきてどうも感覚が違うのは、事務が書類を書いたり見積りをとったりしてくれることがあるということだったりして。これは別に駒場の事務が悪いとか仕事をしないとかいう問題ではなくて、負荷の違いなんだけど。本郷は暇というわけではないけど、駒場は教官も事務官も忙し過ぎるような。僕がいたのはもう 3 年前になるけど、それからより忙しくなってるみたいだし。
論文改訂稿一つ処理。次は短いのいってみよう。
うーん、なんだかくれくれ君モードな私、、、
というわけで人の論文の処理はとりあえずつけたと。自分のが進まない。
8/8 荷物配送 3 号館 14:00-
白川郷から回収。

2002/8/7

頭の中で防犯カメラと TNO 探索が、、、
でもって Klez が、、、 Eudora でもかかりますか。まあ、実行されちゃえば同じか。うっとうしいことこの上ないですね。
8/5 長坂さん午後
すみません。今日の午前中は(珍しく) 1号館のほうにおりました。
なんか知らんけど下旬は会議ばっかりだなあ、、、
今度は 20 ページになってきた、、、でもってほぼ通ったかな論文のほうも改訂版が。ふう。
イギリスでは計算機に蛾が飛び込むらしい。うーむ、単なるファンガードでなくて細かいメッシュとかにしておかないとだめかしら。
さらにまたもうひとつ別の何か、、、
今度は CELL ですか。 TF-1, Blue Gene と来た延長にあるものに PS3 を託すのは無謀だと思うけど、バックアッププランはどうせ準備されてるんだろうね。

2002/8/6

メモ (157.82.60.213) at 00-e0-29-0a-78-c9
ns 上には存在してないではないの。SMC EtherPower II 10/100 であると。
8/6 三鷹
KBO は最近は TNO というらしい。なかなか難しい。それはそれとしてかなり驚くべき結果だよね。
別件で小久保君と相談。時間刻みの多項式で係数を調節する方向で。

2002/8/5

うーむ、やはりというかなんというか、あそこがああなってるのには深い理由があったのね（これでは全く意味不明ですね）。

2002/8/2

Rutgers G6 が動かないとのメイル。理由を適当に推測したら大当たりであった。モジュールに正しくカーネルのオプションを与えればよくて、そのためにはちゃんと kernel の source tree にいれればいいんだろうけど、そうするとまた distribution 依存になったりするし。
しかし、 PIII 700 MHz よかもうちょっと速い機械につけて欲しいです。
ひたすら論文書き。一応今月中に 2 つ形をつけないといかん。細かい仕上げはイギリスいってる時にやるとしても。

2002/8/1

ファイル更新。
長大ドイツ人論文を見る。
コンパックの営業がくる。
```
       HP ZX2000 (Itanic2, 900MHz): 748,000円
       HP ZX6000 (Itanic2, 1GHz): 2,333,000円
       ES45 (EV6 1.25GHz x 4): 14,800,000 円
       
```
0 の数が 1-2 個違うとしか思えない。 zx1 って DDR266 x 2 なのね。で、 900MHz だと SPECfp2000 で 1050 か。 2.4GHz P4 (400MHz FSB) で 825 だから、3GHz がでる秋にはまた P4 に追いつかれるのね。
みーん、論文見るだけで一日終わったぞ、、、 double column で 19 ページもあるし、、、ということは普通に single column double space で出したら字だけで 50 ページ。こないだ見た長いのよりさらにもっと長い。
パソコン物理実地指導のページちょっと更新。見ても怒らないように。

th visit since Apr 1st, 1999.

牧野のホームに戻る