• 計算時間等の求めかた。
  
  % time a.out
  
  計算中に別の処理をしたいときは、
  
  % time a.out &
  
• 計算時間等をfileに書き込むには
  
  % time a.out > cputime
  
  input_file, output_file がある場合は
  
  % (time a.out < input_file > output_file ) >& cputime
  
  計算中に別の処理をしたいときは
  
  % (time a.out < input_file > output_file ) >& cputime &
  
  こうしておかないと処理が済むまで何も出来ません。
• 計算中にlogoutしてその場を離れたい時は、
  例えば、jobgo というスクリプトファイルを
  
  #! /bin/csh
  (time a.out < input_file > output_file ) >& cputime
  
  のように作っておき、
  
  % nohup jobgo &
  
  としてから logout する。但しスクリプトファイル jobgo は
  
  % chmod u+x jobgo
  として実行可のファイルにしておく。
• サーバで動いている全プロセスを見るには(UID, PID, STIME, TIME, CMD 等を確認する)
  % ps -efl
  自分のプロセスだけを見るには
  % ps -efl | grep user-id
• 自分の計算中のプロセス(CMD=a.out, CMD=jobgo等)をとめるには、
  % ps -efl| grep user-id
  で表示された自分のプロセスから、CMD, TIME, PID を確認して、
  % kill PID番号
  自分のプロセスしか止められません。

• 項目
  • 乱数
    1. 一様乱数ξの分布
       ０から１までの分割数を10,20,50,100などとして度数分布を見る。
       サンプリング数を 100, 1000, 10000, 100000,..... と変化させ る。
       度数分布図はグラフにする。
    2. 一様乱数ξの最大値と最小値
       0.0 ≦ drand48() 〈 1.0
    3. 一様乱数ξ1、ξ2から円周率πを求める。
       (ξ1)^2 ＋ (ξ2)^2 ≦ 1.0
       3.14259265....にちかずくか？
    4. 一様乱数ξからサイコロをつくる
       整数部分( 1.0 + ξ )
       
  • 任意の確率密度関数φ(Ｘ)に従う乱数発生(１)
    ξ = φ(ＸminからＸまで積分) / φ(ＸminからＸmaxまで積分)
• 注意
  • cpu time に注意。
• 連絡
  • レポート課題はありませんが、上の項目を自分で試してみること。

• 項目
  • 任意の確率密度関数φ(Ｘ)に従う乱数発生(２)
    1. 逆関数法
       x ＝ Ψ(ξ)
    2. 棄却法
       ξ1, ξ2
    3. 数表による方法
       ξと x の対応表を予めつくっておき、この対応表から ξに対する x を補間法により求める。
       1. XminからXmaxを細かく分割し、 XminからXiまでを数値積分する。
       2. Xiと ξi との対応表を作成する。
       3. 任意のξに対する x を lagrange補間により求める。
       ★ 数値積分(シンプソン法、「台形近似、二次曲線近似」)
       ★ 補間(一次補間、二次補間、ラグランジェ補間)
• 注意
  • Nov.11,'99に配布したプリントの飛来方 向を表す方向余弦(a,b,c)の求めか方に誤りがありました。訂正プリン トを配布しました。このプリントにしたがって、飛来方向を sin(θ), sin(Φ)を求めてみて下さい。今度は上手くいくはずです。(θ,Φ)プロッ トを期待しています.
• 連絡
  • φ(x)＝ 1/x の場合に付いて 逆関数法 及び 棄却法 による乱数 発生を試みること。余力のある学生は 数表による方法 をここみること。 Ｘminはいろいろ(0.1, 0.01, 0.001,...) 変えてみる。

第３回目のレポート課題を出題しました。 締め切りは12月17日(金)

• 項目
  1. モンテカルロ法(1)
  2. 正規乱数を用いたモンテカルロシミュレーション
     • Nov.18,1999 配布のプリント
  3. 指数乱数を用いたモンテカルロシミュレーション
     f(t) ： x まで変化せずにいる確率
     σ ： 単位長(時間、長さ)あたりに変化する確率
     • f ( t + dt ) = f ( t )×( 1 - σ×dt )
     • df / dt = - σf(t)
     • f ( t ) = σ× exp ( -σ×t )
     • t の平均値 λ = 1 / σ
     • f ( t ) = ( 1 / λ ) × Exp ( - t / λ )
     • t = - λ ln (ξ)
• 注意
• 連絡
  • レポート課題(３)
    Nov.18,1999 配布のプリントの課題です。 締め切りは12月17日(金)。
    遅れる場合は12月22日までとします。

• 項目
  1. 方位角Φ(Nov.25配布プリント)の補足説明
     • sin(Φ)の逆関数からは -π/2 〜 π/2
     • t1,t2,t3,t4 の値から、Φ と π-Φ を判定する。
     • 上の判定をせずに、(Θ,Φ)プロットを描いて良いことにします が、
       余力のある学生は挑戦してみること。
  2. モンテカルロ法(2)
     • 増殖過程及び減衰過程
     • 2倍増殖のフローチャート例のプリントを配布
       変数 I の使い方を理解すること。
       Ep = 1000
       Emin = 10
       tp = 0
       tmax = 100
       λ = 1.0
       Nmax = 100
       t = -λ×ln(ξ1)
       dE = E × ξ2
       E = E - dE
       
• 注意
  • 増殖及び減衰がアルゴリズムどおりに計算 されていることを必ず確認すること！
  • Ep÷Emin や Nmaxを大きくすると計算時間が 膨大になるので、少しずつ大きくすること。
• 連絡
  • 今週もレポートはありません。

• 項目
  1. 天頂角Θ
  2. 方位角Φ
• 連絡
  • レポート課題(４)
    課題内容は第３回(Oct.21,1999)の項目４を参照せよ。
    
  • 30000 events のデータを後日e-mailで送ります。
    
  • taka@si.hirosaki-u.ac.jp 宛に "Send me 30000.dat"
    という Subject のメールを送信すること。
    このメールを受信した学生に 30000.dat 送信します。
    
  • このデータでは、TDC 1カウントは 100ピコ秒です。
    
  • 11月25日に配布したプリントを参照すること。
    天頂角分布(x軸:cos(Θ)　対数目盛 、y軸:度数　対数目盛)
    方位角分布(x軸:Φ 、y軸:度数)
    時間差分布(x軸:dt　対数目盛 、y軸:度数　対数目盛) -->レポート課題(1)
    飛来度数分布(x軸:日付 、y軸:度数)
    など
  • 締め切りは2月10日(木)。 遅れる場合は2月17日(木)までとします。

• 項目
• 注意
• 連絡

• 項目
• 注意
• 連絡