差分方程式に変換して数値計算を行う
CCS HPCサマーセミナー 並列数値計算アルゴリズム
... まとめ（1/2） • 物質科学の実アプリケーションにおいて使われることが 多い，高速フーリエ変換（FFT）について紹介した． • これまで並列FFTで行われてきた自動チューニングで は，基数の選択や組み合わせ，そしてメモリアクセスの 最適化など，主にノード内の演算性能だけが考慮され てきた． ...
51
数値計算：フーリエ変換
... 8 点高速フーリエ変換： 二回の 4 点高速フーリエ変換 二回の 4 次の対角行列とベクトルの積 4 次の対角行列とベクトルの積：4 回の乗算 (8 点 FFT の乗算回数) = 2 × (4 点 FFT の乗算回数) + 8 4 点高速フーリエ変換： 二回の 2 点高速フーリエ変換 ...
84
CCS HPCサマーセミナー 並列数値計算アルゴリズム
... まとめ（1/2） • 物質科学の実アプリケーションにおいて使われることが 多い，高速フーリエ変換（FFT）について紹介した． • これまで並列FFTで行われてきた自動チューニングで は，基数の選択や組み合わせ，そしてメモリアクセスの 最適化など，主にノード内の演算性能だけが考慮され てきた． ...
51
カーネルベンチマークコード 開発の目的 エクサスケール規模のシミュレーションの核となる数値計算アルゴリズムの中で 特に重要なものについて 数値計算ライブラリ等を用いてそのコストを推定するためにカーネルベンチマークを作成し 評価に使用する 対象計算アルゴリズム 固有値計算 ( 実数密行列 標準固有値計
... 行列の固有値計算では行列を簡単な形式（形状）に変換し、それを中間形式と して取り扱います。理研の研究チームは、帯行列（ゼロでない要素が対角線上 に帯状に分布する行列）を中間形式に採用することによって、前処理の時間の ...
26
数値計算で学ぶ物理学 4 放物運動と惑星運動 地上のように下向きに重力がはたらいているような場においては 物体を投げると放物運動をする 一方 中心星のまわりの重力場中では 惑星は 円 だ円 放物線または双曲線を描きながら運動する ここでは 放物運動と惑星運動を 運動方程式を導出したうえで 数値シミュ
... ②①で作成したプログラムを、時刻、 座標、 座標、水平速度、鉛直速度に加えて、これ らの値を用いて評価した単位質量あたりの運動エネルギー 、ポテンシャルエネルギー 、 力学的エネルギー も出力するように改良して実行せよ（プログラムや計算結果の提出は ...
11
アウトライン ( 古典 ) 分 動 学計算とは ( 古典 ) 分 動 学計算の概要 運動方程式 数値解法 場 境界条件 実例 タンパク質 リガンド間の結合自由エネルギー予測 (MP-CAFEE 法 )
... ff03 溶媒中での環境を直接 QM で計算し、電荷と二面角を求める アラニンテトラペプチドの構造式 Comparison of Multiple Amber Force Fields and Development of Improved Protein Backbone Parameters, Viktor Hornak, Robert Abel, Asim Okur, ...
39
Microsoft Word - 03-数値計算の基礎.docx
... 直接法で解くと時間がかかりすぎることがある。このような場合には，計算を繰り返し て解を十分精度の良いものに近づけていく方法が用いられる。このような方法は反復法 と呼ばれる。ここでは，２つの代表的な方法を採り上げる。 ...
17
数値計算：常微分方程式
... 常微分方程式の標準型 MATLAB ファイル pendulum.m pendulumConstants; % 単振り子のパラメータを定義する timeinterval=0:0.1:10; % 固定ステップ (ステップ幅 0.1s) q0=[pi/3;0]; % 状態変数ベクトルの初期値 ...
86
複雑な立体形状を捉える [3 次元計測 ] カラーカメラでの画像処理では 侵入者を抽出するために 画像の変化を検出する差分処理 ( 背景差分 フレーム間差分 ) を行う しかし 侵入者が 背景と類似した服を着用していると 輝度の違い ( コントラスト ) がなくなり 分離が難しくなる 3 次元情報が
... ダムパターンを照射することで、物体にはテクス チャがなくても距離計測ができるアクティブ型が 研究されている。 光切断法は、精度の高い立体データの取得方法 として知られている。レーザーなどで立体を切断 するように光を照射し、カメラから見た位置から ...
6
としてもよいし,* を省略し, その代わりにスペースを空けてもよい. In[5]:= 2 3 Out[5]= 6 数値計算 厳密な代数計算 整数および有理数については, 厳密な計算がなされます. In[6]:= Out[6]= In[7]:= Out[7]= 2^
... : 方程式を意味する等号は = ではなく = = （イコールが２つ）を用います． 解が代数的に求まらないような方程式に対しても，関数 FindRoot を用いて数値的な近似解を求めることができま す．例えば方程式 Sin[x] + 1 = x は（下のグラフからも分かるように） x = ...
11
光電気変換型コヒーレント光信号再生の研究 Study on Opto-Electronic Coherent Signal Regeneration (1) 光信号を電気信号に変換せず 光のままで雑音除去や波形整形を行う全光型の再生器と (2) 光信号を一旦電気信号に変換して再生処理を行い その後に
... Hu, K. M. Røge, M. Galili, T. Morioka, and L. K. Oxenløwe, "Scalable WDM phase regeneration in a single phase-sensitive amplifier through optical time lenses", Nat. Commun., vol. 9, article #1049 (2018). では毎秒 10 ...
5
数値計算ライブラリの使用方法 「実習編」
... - “N”,”N”は転置するかしないかを指定できる - 1,1を指定している部分は部分行列のとき開始インデックスを指定できる - 今回行列の形状は同じなので DESCAをすべてで使用（それぞれ変更可能) call MPI_Barrier ( MPI_COMM_WORLD, ierr ) z1 = MPI_Wtime() ...
31
ST 計算 1: このST 計算は 放射トラハ ー計算や 開放トラハ ー計算で新たにトラハ ー杭を設置した後 逆放射トラハ ー計算や後方交会計算を行う前の杭間の距離を確認する時に使用して下さい 2:8の時 方向角も出る様にして有りますが 参考値程度でお考え下さい PROGRAM-NO.1 番 No
... 継続事業の場合、内容（精度）に準じて併合ﾄﾗﾊﾞｰや結合ﾄﾗﾊﾞｰ計算で、 補正をして誤差を拡散してある事が、あるからです。 又、事業の内容によっては、球面補正（地球は丸い）をして精度を上げてある現場もあります ...
7
再分割アルゴリズムを用いた力学系の定常分布の数値計算
... この例では，アトラクタの計算時間を約 13% 短縮でき た（表 1 ）．また，定常分布の精度を保持したまま計算 時間を約 70% 短縮できた（表 2 ）． 4. 4 計算例： Duffing 方程式のポアンカレ写像 非自律系の微分方程式の例として，次式で与えられ る Duffing ...
8
数値計算による瀬戸内海の潮汐・潮流計算と海水交換について
... 潮汐残差流の計算結果についてもう少し詳細に見 るために，大阪湾の Case - 1 とCase-2 の潮汐残差流 の結果を拡大して Fig.8 に示す.時計回りの循環流 が卓越するという傾向はどちらのケースの計算結果 にも見られるが，合成潮による Case - 2 の潮汐残差 流の方が幾分大きな値となっている.この傾向は， 他の海域(例えば健灘)の潮汐残差流の結果でも同[r] ...
11
JAXA Repository AIREX: スカラー保存則の保存型差分近似の数値粘性とエントロピー条件適合性
... 22 ただし、 " の値を非常に小さくとった場合、現実に数値計算が行われる有限の格子刻みで得られる数値解では、膨張波 における膨張の最初の過程が厳密解に比べて遅くその形も不自然になることが知られており、そこには定性的な解析と ...
29
連立1次方程式Ax=bの解法：公式にしたがって解くのは，計算量大
... 1. Lisp のしくみ Lisp で中心となるのは，S 式(Symbolic Expression)と呼ばれる記号の構造的な表現である．Lisp ユーザはインタープリタを使って， S 式を作り出したり，変形したり，探索したりできる． Lisp の特色の一つに，リフレクションがある．S ...
13
圧縮性LESを用いたエアリード楽器の発音機構の数値解析 (数値解析と数値計算アルゴリズムの最近の展開)
... $T_{ij}$ を計算し，その偏微分を取 ればよい．したがって，形式的にではあるが， Lighthill の方程式を使って流体の運動がどの様に 音源になるかを調べる事ができる．このような理由で， Lighthill ...
11
ランダムウォークの境界条件・偏微分方程式の数値計算
... ランダムウォークの境界条件・偏微分方程式の数値計算 ランダムウォークの境界条件 ランダムウォークの端スペシャルルール=境界条件 1 ≤ x ≤ m − 2 は普通の場所 , 端 x = 0, x = m − 1 は特別 ( 壁 ) と思おう . x = 0 でのありうるルール = 境界条件 . 壁はランダムウォークの時の言葉 . ...
38
2.2 支配方程式および離散化手法本研究では, キャビティ流れにおける流れと音の連成を再現するため, 流れと音の直接数値解析を行う. 支配方程式は式 (1) に示す三次元圧縮性 Navier-Stockes 方程式であり, 有限差分法による直接計算を行った. Qt x k F k F 0 νk ここ
... 0.5 にならないため， 外挿により推定される値である． Tu0.0Pの流れ場の可視化においてみられたスパン方向に均一な渦 構造に関して， 図 12 からスパン方向全体に速度変動が同期してい ることがわかる．図 12 の ...の比較に関 して，流れ場の可視化からは明らかにならなかった三次元的な流 ...
8