数値計算によればこの
C による数値計算法入門 ( 第 2 版 ) 新装版 サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. このサンプルページの内容は, 新装版 1 刷発行時のものです.
... なお，数値計算のプログラムを学生に一つひとつ作成させるのは，授業時数や学生 がパソコンに費やすことのできる時間の関係から，そう多くを期待できないのが実状 である． それで，指導者のほうでいくつかを取捨選択して，学生に全部または一部を作成さ せるというのが現実的であり，その他のプログラムは指導者が学生に与えて，パソコ ...
33
エクセルによる数値計算と化学への応用
... エクセルにはマクロ（エクセル中の一連の操作をまとめてひつつのコマンドで行うことのできる機 能）としてＶＢＡ(Visual Basic for Application)を持っています。これはＢＡＳＩＣ（プログラム言語の ひとつ）をエクセルの操作や処理の自動化ができるようにしたものです。その意味でＶＢＡはプログラ ムの作り方として Visual Basic ...
12
178 宇宙航空研究開発機構特別資料 JAXA-SP 図 2 ステ ジング型燃料ノズル [2] 図 3 計算格子概略図 3. 数値計算数値解析には NuFD/ Front Flow Red を用いた. LES の支配方程式を以下に示す. 以下の式において, は化学反応による化学種の生成
... ん断層によって微粒化と混合を促進する機構となっている． この燃料ノズルが燃焼器ライナーに取り付けられている． ライナーの上流側 89.5mm の部分は矩形断面を有している． 燃焼室下面には，ステンレス製板が用いられ，燃焼室入口 から 17mm の位置に圧力センサーが取り付けられている． 燃焼器ライナーの下流側 105mm は縮流ノズルとなってい ...
6
数値計算による瀬戸内海の潮汐・潮流計算と海水交換について
... 潮汐残差流の計算結果についてもう少し詳細に見 るために，大阪湾の Case - 1 とCase-2 の潮汐残差流 の結果を拡大して Fig.8 に示す.時計回りの循環流 が卓越するという傾向はどちらのケースの計算結果 にも見られるが，合成潮による Case - 2 の潮汐残差 流の方が幾分大きな値となっている.この傾向は， 他の海域(例えば健灘)の潮汐残差流の結果でも同[r] ...
11
としてもよいし,* を省略し, その代わりにスペースを空けてもよい. In[5]:= 2 3 Out[5]= 6 数値計算 厳密な代数計算 整数および有理数については, 厳密な計算がなされます. In[6]:= Out[6]= In[7]:= Out[7]= 2^
... &=2.71828... に収束する様子を視覚化したい場合には，次のようにします． In[39]:= napier ' Table "$1 $ 1 ! n%^n, &n, 1, 50'#; Show "ListPlot"napier, PlotStyle ( Red#, Plot"Exp"1#, &x, 0, 50', ...
11
ウォッシュロードによる貯水池堆砂の1次元および2次元数値計算法
... due to wash ioad is presented as weH as a l‐ dimensional numerical method 「Γ he 2‐ dilnensional rnethod is based on 3‐ dilnensional rnodels for the flow and the sedilnent concentration`[r] ...
12
CCS HPCサマーセミナー 並列数値計算アルゴリズム
... • ブロックサイズNBについても探索を行うことによって， 最適なブロックサイズを調べることができる． • 今回の実装では，データサイズ およ びMPIプロセス数 が2のべき乗であると仮定してい るため，ブロックサイズNBも2のべき乗に限定して2， 4，8，16，32，64のように変化させている． ...
51
圧縮性LESを用いたエアリード楽器の発音機構の数値解析 (数値解析と数値計算アルゴリズムの最近の展開)
... 発生する渦のスケールはそれよりも遥かに小さい．したがって，流体の細部の構造を正確に再現 するには，音波の解析に使われるメッシュよりもより細かなメッシュが要求される． そこで，メッシュ数と計算時間を節約するために，楽器のサイズをできるだけ小さくし，図 2 に示すような二次元の数値モデルを考える．三次元の等価なものを考えると， ...
11
カーネルベンチマークコード 開発の目的 エクサスケール規模のシミュレーションの核となる数値計算アルゴリズムの中で 特に重要なものについて 数値計算ライブラリ等を用いてそのコストを推定するためにカーネルベンチマークを作成し 評価に使用する 対象計算アルゴリズム 固有値計算 ( 実数密行列 標準固有値計
... EigenExaについて • 理化学研究所ホームページ 2013年12月5日 60 秒でわかるプレスリリース 」より • 「京」を使い世界最高速の固有値計算に成功 • 行列の固有値計算では行列を簡単な形式（形状）に変換し、それを中間形式と して取り扱います。理研の研究チームは、帯行列（ゼロでない要素が対角線上 ...
26
産業動学に関する研究ノート（数値計算編）（1）
... と定まる。プログラムにおいて，この確率は関数 chkentry を呼び出すことで計 算される 16) 。 ・価値関数および投資関数の推測値 プログラムでは，関数 contract を１回呼び出すことで，反復が１回実行され る。関数 contract は，各々の産業構造 z { 1 , , w } に対して関数 optimize を その都度呼び出す。関数 optimize は，投資関数行列の第 z 行[x(z, ...
25
再分割アルゴリズムを用いた力学系の定常分布の数値計算
... Yu NUREKI †a) and Sunao MURASHIGE †b) あらまし 本研究では，カオス的な離散力学系の定常分布に対する有効な計算手法を提案する．定常分布は， 力学系の定常状態において状態空間の各点に軌道が訪れる確率分布を表し，系の複雑な挙動を特徴づける．定常 分布の計算方法として，状態空間を有限個の小領域に分割し，系の時間発展を有限マルコフ連鎖で近似し，定常 ...
8
ランダムウォークの境界条件・偏微分方程式の数値計算
... ランダムウォークの境界条件・偏微分方程式の数値計算 状態数が大きく規則的なマルコフ連鎖の時間発展の数値計算 プチテスト (筆記) やります! 2018-06-05 火 5, 10 分 ( 外部記憶ペーパー作成 )+80 分 ( 筆記 ), 20 ピーナッツ . 配 点も範囲も昨年度と異なります . 介護等体験などで欠席する人は事後に届を出せ ...
38
数値計算法
... 最尤法の直接利用と最小二乗法 最小二乗法を使えないとき＝分布が正規分布でないとき ビンまとめし、ヒストグラムをつくると、１ビンあたりに含ま れるデータ数が十分大きい場合、正規分布で近似できる。 この場合最小二乗法が使えるようになる。 ...
28
Microsoft Word - 03-数値計算の基礎.docx
... 近接するいくつかの区間を１まとめにして，より高次の次数をもつ多項式で関数を近似 する方法である。精度は補間する関数の次数で決まる。例えば，区間[a,b]を偶数個ｎに 等分し，２つまとめにすると，３点を用いて２次多項式を当てはめることができる。こ のとき，小区間を２つ合わせた区間の面積は、 ...
17
2011 年ハイパフォーマンスコンピューティングと計算科学シンポジウム High Performance Computing Symposium 2011 HPCS /1/19 数値計算ポリシーインターフェース付 行列計算ライブラリの開発と評価櫻井隆雄 直野健 片桐孝洋 中島研吾 黒
... 3 にライブラリを呼び出す呼出部の動 作を示した．呼出部は POLICY と EPS の値により 再直交化の実装方式とライブラリへの再投入の有 無といった動作の違いがある． 最直交化方式は EPS が ...際 に POLICY が ACCURACY の場合は改めて検算∥Ax'–b∥（A:入 力行列，x':ライブラリが出力した解，b:右辺項）を 行い，解の信頼性を確認する．このとき，解が要求 ...
9
数値計算法
... 35.5,25.0,34.2,34.6,22.8,27.7,30.6,26.8,23.0,39.3 出力：（標本）平均値、標準偏差 ソースプログラムと出力結果をメイルの本文にして [email protected] までメイルせよ。 実行形 式は添付しないこと。 メイルのタイトルは report1_学籍番 号とすること。 ...
20
数値計算ライブラリの使用方法 「実習編」
... 基本データは２次元配列として、プロセスグリッド上に分散配置 ブロックサイズ (NB)を基本単位として、大きな行列の対応箇所のみ保有 FORTRANのC-Major方式に則って、同一列のブロックは１次元目が連続メモリ アクセスとなるようにメモリ上に格納される ...
31
マルコフ連鎖の時間発展の数値計算
... 7 ( 3 4 ) のみであ り , これが唯一の定常分布 . 2 ⃗ u 1 は確率ベクトルとなるように s = 1 7 ととる . ⃗ u 2 は確率ベクトル でないので適当に s = 1 ととる . 他の取り方をしても最終的には同 じ ⃗ p(t) が求まる . ...
20
CCS HPCサマーセミナー 並列数値計算アルゴリズム
... が一次元分割よりも性能が高くなっていることが分かる． • 二次元分割を行った場合でも，4,096コアにおいては 96%以上が通信時間に費やされている． – 全対全通信において各プロセッサが一度に送る通信量が わずか1KBとなるため，通信時間においてレイテンシが 支配的になるためであると考えられる． ...
51
数値計算：常微分方程式
... 常微分方程式の数値解法 長所 解析的に解けない微分方程式を解くことができる シミュレーション（力学，回路，CG） 短所 方法によって解が異なる ステップ幅によって解が異なる ステップ幅が大きい 計算した解が誤っている ステップ幅が小さい 計算時間が要する ステップ幅を小さく (たとえば半分) にして解が一致 ⇓ もとのステップ幅で計算してよい． 平井 慎一 (立命館大学[r] ...
86