波数の部分空間での数値計算
ଗȨɍɫȮĘർǻ 図 : a)3 次元自由粒子の波数空間におけるエネルギー固有値の分布の様子 b) マクロなサイズの系 L ) における W E) と ΩE) の対応 として与えられる 周期境界条件を満たす波数 kn は kn = πn, L n = 0, ±, ±, 7) となる 長さ L の有限
... 3 古典力学から熱力学へ これまでは，量子力学から出発して，熱力学量を計算する方法を紹介してきた．ここでは，古典力学から 出発して熱力学量を計算する方法を考えよう．量子力学では粒子のエネルギーを求めるためにシュレディン ガー方程式をわざわざ計算する必要があった．しかも，シュレディンガー方程式を手で解くことができのは ...
10
演習 II 2 つの講義の演習 奇数回 : 連続系アルゴリズム 部分 偶数回 : 計算量理論 部分 連続系アルゴリズム部分は全 8 回を予定 前半 2 回 高性能計算 後半 6 回 数値計算 4 回以上の課題提出 ( プログラム + 考察レポート ) で単位
... int source, int tag, MPI_Comm comm, MPI_Status *status) ● buf: 送受信バッファ, count: 要素数, datatype:型, dest/source: 送信先/受信先, tag: 識別子, status: 受信結果 ● 以下の例ではデッドロックします ...
18
178 宇宙航空研究開発機構特別資料 JAXA-SP 図 2 ステ ジング型燃料ノズル [2] 図 3 計算格子概略図 3. 数値計算数値解析には NuFD/ Front Flow Red を用いた. LES の支配方程式を以下に示す. 以下の式において, は化学反応による化学種の生成
... m の円筒に含まれる領域をパイロ ット領域，そこに含まれない外側の領域をメイン領域と定 義した．パイロット領域・メイン領域に分割して算出した 温度分布（断面平均）を図 9 ，速度分布（断面平均）を図 10 にそれぞれ示す．ここで，本研究で対象とした燃焼器 では，メイン流路・パイロットカップは ...
6
2011 年ハイパフォーマンスコンピューティングと計算科学シンポジウム High Performance Computing Symposium 2011 HPCS /1/19 数値計算ポリシーインターフェース付 行列計算ライブラリの開発と評価櫻井隆雄 直野健 片桐孝洋 中島研吾 黒
... EPS の値により 再直交化の実装方式とライブラリへの再投入の有 無といった動作の違いがある． 最直交化方式は EPS が ...TIME でな い場合は速度と正確さのバランスのとれたMGS を 用い，それ以外の場合は速度重視の BCGS を用い ...
9
図 宇宙論解析の流れ 次元データの CMB の例 宇宙論ゆらぎ場 F (θ) の測定 左上図 ゆらぎ場のフーリエ波数分解 右上図 右下図は パ ワースペクトル推定の結果 灰色点は各波数ビンでの測定値 エラーバーを伴う青点は 複数の波数ビンで測定値を平均した結果 エラーバーとして 有限数のフーリエモー
... パワースペクトル推定の統計誤差が与えられれば、 理論モデルと比較することができます。 Pˆ F (l) ←→ P F model [l; P ζ (k), Ω m0 h 2 , Ω b0 h 2 , Ω de , ...] (7) 右辺は、理論モデルによる F 場のパワースペクトル が、原始パワースペクトル、その他宇宙膨張を記述す るための宇宙論パラメータの関数として与えられるこ ...
8
しかし わたしたちの日常生活では 大地は平面であると思って特に差し支えはありませ ん わたしたちの生活空間に対して地球は充分に大きいので 球面の一部であることを地 上で肉眼的に認識するのは非常に難しいのです 地球の層構造 地球は中心部分を金属鉄の核 その周りを岩石でできたマントルが囲み 表層部分は地
... ・武蔵野台地の３層構造 等々力渓谷に見られる地層は、大きく分けて、下から海成の泥岩層、武蔵野礫層、関東 ローム層の３つがあります。地層がたまった順番もこの順で、一番下の海成泥岩は百数十 万年前の上総層群の地層と、12 万年ほど前の東京層（渋谷粘土層）の２つがありますが、 ...
12
カーネルベンチマークコード 開発の目的 エクサスケール規模のシミュレーションの核となる数値計算アルゴリズムの中で 特に重要なものについて 数値計算ライブラリ等を用いてそのコストを推定するためにカーネルベンチマークを作成し 評価に使用する 対象計算アルゴリズム 固有値計算 ( 実数密行列 標準固有値計
... EigenExaについて • 理化学研究所ホームページ 2013年12月5日 60 秒でわかるプレスリリース 」より • 「京」を使い世界最高速の固有値計算に成功 • 行列の固有値計算では行列を簡単な形式（形状）に変換し、それを中間形式と して取り扱います。理研の研究チームは、帯行列（ゼロでない要素が対角線上 ...
26
ランダムウォークの境界条件・偏微分方程式の数値計算
... ランダムウォークの境界条件・偏微分方程式の数値計算 状態数が大きく規則的なマルコフ連鎖の時間発展の数値計算 状態数が大きく規則的なマルコフ連鎖の時間発展の数値計算 例 : ランダムウォークや偏微分方程式 ...
38
計算での注意 : 文字式で計算し 数値計算は最後にする 文字にはそれぞれ意味がある 質量 ss 力 foc 加速度 cclion 速度 loci など 質点 : 大きさの無視できる物体 質量をもつ 自身の周りの回転運動は考えない 例えばコマは その位置を変える事なく回転運動し その運動エネルギーを持
... M の斜面 ABC に沿った運動であり、物体 m が台 M を押す力や台 M が物体 m を押す垂直抗力（作用・反作用の関係にある 2 つの力）とは、常に直交する。よって、これら 2 つの力が行う仕事の量は、いずれにせよゼロである。 ...
204
再分割アルゴリズムを用いた力学系の定常分布の数値計算
... Yu NUREKI †a) and Sunao MURASHIGE †b) あらまし 本研究では，カオス的な離散力学系の定常分布に対する有効な計算手法を提案する．定常分布は， 力学系の定常状態において状態空間の各点に軌道が訪れる確率分布を表し，系の複雑な挙動を特徴づける．定常 ...
8
Ver.10( 住宅 住戸の外皮性能の計算プログラム Ver.02.01~) 第三節熱貫流率及び線熱貫流率 1. 適用範囲 この計算は 用途が住宅である建築物又は建築物の部分における 外皮の部位の熱貫流率及び熱橋等の線熱貫流率の計算について適用する 2. 引用規格 JIS A1420:
... ② 評価品に二重ガラス構造の建具（複数枚の単板ガラス等を用いた一重建具）を使用している場合は、 評価品の中空層の層数が試験体のそれ以上であり、かつ、中空層の厚さが試験体のそれ以上であるこ と。 試験体にガラスを使用し、評価品に腰パネル等を使用している場合は、腰パネル中央部の熱貫流率を一次 ...
43
アウトライン ( 古典 ) 分 動 学計算とは ( 古典 ) 分 動 学計算の概要 運動方程式 数値解法 場 境界条件 実例 タンパク質 リガンド間の結合自由エネルギー予測 (MP-CAFEE 法 )
... バイオグリッド HPCI プロジェクト「新薬開発を加速する「京」インシリコ創薬基盤の構築」 プロジェクト「新薬開発を加速する「京」インシリコ創薬基盤の構築」 プロジェクト「新薬開発を加速する「京」インシリコ創薬基盤の構築」 プロジェクト「新薬開発を加速する「京」インシリコ創薬基盤の構築」 KBDD (K supercomputer-based drug discovery ...
39
密集市街地における換気・通風性能簡易評価ツールの開発 (その2 流体計算部分の開発)」
... 得られればその他のパラメータを算出することが出来る。 本ツールでは計算時間短縮のため建物内に完全に含まれるセルと建物境界面を含むセルにそ れぞれ解析領域全体での平均値を一律に与える簡易的な方法を用いている。 建物境界面を含むセル 建物内に完全に含まれるセル ...
20
CCS HPCサマーセミナー 並列数値計算アルゴリズム
... • これまでに，MPIの集合通信を自動チューニングする 研究が行われている[Faraj and Yuan 05]． • InfiniBandで接続されたマルチコアクラスタにおいて， 全対全通信をノード内とノード間の2段階に分けて行う ことで，性能を向上させる手法も知られている[Kumar et al. 08]． ...
51
CCS HPCサマーセミナー 並列数値計算アルゴリズム
... • これまでに，MPIの集合通信を自動チューニングする 研究が行われている[Faraj and Yuan 05]． • InfiniBandで接続されたマルチコアクラスタにおいて， 全対全通信をノード内とノード間の2段階に分けて行う ことで，性能を向上させる手法も知られている[Kumar et al. 08]． ...
51
エクセルによる数値計算と化学への応用
... 表示の部分の入力は簡素化しようがありませんが、四則計算の部分ではたとえば最初の加算の 部分に数値を２つ入力すれば加減乗除算を自動で行うことはできます。また、数列の和も一番右 下のセルに公式を入力しておけば自動的に計算してくれます。 ...
12
フジタ技術研究報告第 48 号 2012 年 CUDA を用いた並列数値解析手法の一考察 仲沢武志 概 要 従来 大規模な問題や複雑な形状が計算対象の場合では スーパーコンピュータのような高価な計算リソースを使用する必要があった これに対し コンピュータのハード的な進化によって 現在では市販のパーソ
... 例えば実際に物質移行拡散解析を地下水の汚染拡散 の予測に適用する場合などでは、空間メッシュは工場が 立地されるような規模の土地面積であり、時間的には戦 後から 60 年オーダを対象とするようなこともある。このよう な場合には、単一ステップの計算時間は短くても、時系 ...
6
マルコフ連鎖の時間発展の数値計算
... マルコフ連鎖の時間発展の数値計算 状態数が大きく規則的なマルコフ連鎖の時間発展の数値計算 L07-Q1 Quiz( ランダムウォークの時間発展 ) 次の推移確率行列を持つ, 状態空間 {x} = {0, 1, 2, . . . , 99} 上のマルコフ連鎖を考 えよう. ...
20
数値計算ライブラリの使用方法 「実習編」
... 簡単な行列積 これまで行列データは形状と開始インデックスで指定 (A(1,1) or A(1,2)) 配列 Aと分散情報を保持するデスクリプタの組で 注意: プログラム上の配列Aはローカルに確保された配列を代表しているので A(1,1)は現在の ...
31
000/6/3 目次 地球シミュレータの紹介 大気大循環モデル 基礎方程式 スペクトル法 コアとなる計算部分 並列化手法 実行性能 問題点と今後の予定 ワークショップ 計算科学におけるアルゴリズム
... 「地球シミュレータ」計画 ■ 地球変動現象の解明・予測の研究開発を推進 ■ 大規模，高精度シミュレーションにより「仮想地球」を実現 数値シミュレーション 数値シミュレーション 数値シミュレーション 数値シミュレーション ...
32