B_T[][] を、ローカルな行列-行列
Title 拡張クロスデータ行列法と共分散行列関数の不偏推定 Author(s) 矢田, 和善 ; 青嶋, 誠 Citation 数理解析研究所講究録 (2015), 1954: Issue Date URL
... いとき， $tr(S_{n})^{2}$ は大きなバイアスをもつ．簡単なシミュレーション実験で検証する． 母集団分布は $N_{p}(0, \Sigma)$ とし， $\Sigma=0.5I_{p}+0.51_{p}1_{p}^{T}$ なる級内相関モデルを考える． ただし， $1_{p}$ は $1_{p}=$ $(1, 1)^{T}$ なる $P$ 次ベクトルである．このとき，最大固有値は ...
11
書き換え可能なゲート素子を持つデバイスを用いた行列計算専用集積回路の設計
... によって記述し、シミュレーションによって、この計算アルゴリズムをハードウェア レベルで動作させるためのモデルを作成した。 '98 年度は、山岡 [5] と私により、先に製作された基板を利用してハウスホルダ法の アルゴリズムを使い、実際に行列の固有値と固有ベクトルの計算をハードウェア上で 動作させた。まず、基板と PC ...
175
群論モデルにもとづく正方行列枠中の要素パターンに対する良さと複雑さ
... にもとづき ，ホログラフィックモデルは1次元の文字列 （string）の形式化にもとづいていて，それらを2次元の 問題に適用している。 パターン認知については変換群のみならず枠組によっ て 影 響 を 受 け る こ と が 知 ら れ て い る（Hamada et al., 2016）。たとえば濱田他（2013）によると，鏡映対称軸 の本数が同じ場合で正六角形枠の鏡映軸2本のパターン の良さ（D 2 : ...
10
超高速行列演算チップの作成
... したら固有値固有ベクトルを求めることができるかを考えるのが本研究の目的であ る。そしてその過程が正しいか確かめるために C 言語を用いプログラムでシミュレー ションすることにした。これにより固有値固有ベクトルを求める際の過程がわかり、 どの部分が重要であるか、またこの過程は時間がかかるため専用チップ上で行うよ うにして速く計算させるなど専用チップの制作へスムーズに移行しやすくなる。 ...
81
ハードウェア記述言語を用いた行列演算専用回路設計
... Project” を選択． 図 B.2: “メニュー” の “ファイル” を選択 22 • 下のウインドウが表われたら，右クリックを押して”Add Module” を選択，ダイアログ ボックスが表われたら，”The Project has not been ...で，”OK” を押す，するとダイアログボックスが表われるので，ACC(*.acc) ファイルに名 ...
73
FPGAを用いた行列計算専用プロセッサの設計
... [1] 、また、計算アルゴリズムの並列化の一つのモデルを作成した [2] 。 専用計算機の計算対象、計算アルゴリズムが決定したところで、次は計算アルゴリ ズムを実際にハードウェア化することが課題である。ハードウェア化に際して、まず 検討しなければならないことは設計手段と使用するテクノロジーであるが、前節で述 ...
139
3061 チュートリアル Fortran FORmula TRANslation Fortran90/95 Fortran90/95 1 HP 有限要素計算における全体剛性行列の作成法 疎行列データ構造の視点から 永井学志橋本一輝 1 はじめに FEM SIMD FEM PDE FEM FEM FE
... Fortran （FORmula TRANslation） は半世紀以上の歴史を持ち、一部では時代遅れと言われなが らも、今なお数値計算に利用する研究者が多いプログラミング言語です。Fortran90/95の機 能や使用例を改めて理解したいという研究者のために、Fortran90/95による近年の有限要素 ...
5
ハードウェア記述言語を用いた行列計算専用プロセッサの設計
... 例するため，次数が数千にも及ぶ実際の数値計算では実用に耐えないほど非常に多 くの時間を必要とする． この計算時間を短縮しようと，これまで並列化コンピュータを用いたソフトウェ ア的な手法が試みられてきたが，やはり限界があり，それほどの改善は見込まれな い．その原因として，現在使われているワークステーションなどは，どの問題にも ...
89
インテル MKL を使用した小行列乗算の高速化 インテル MKL チーム
... インテル® Xeon® Platinum プロセッサー上での MKL_DIRECT_CALL のパフォーマンス システム構成: ハードウェア: インテル® Xeon® Platinum 8180 プロセッサー、2x28 コア、2.50GHz、192GB RAM。オペレーティング・システム: Red Hat* Enterprise Linux* 7.2 LTS。ソフトウェア: インテル® MKL 2018。性能に関するテストに使用され ...
30
表計算による待ち行列のシミュレーション テーマパークアトラクション待ち時間の再現 ( 坂本 ) 表計算による待ち行列のシミュレーション テーマパークアトラクション待ち時間の再現 坂本 紫 Simulation of the queue by the spreadsheet Reproduction
... 文系学部における情報教育として，プログラミング を取り上げている大学は少ないなか[1][2]，プログ ラミング教育を実施するため，次の2つの理由から Excelマクロ（Excel VBA）を取り上げている教育 が紹介されている[1][3]。(1) 文系学部生は初年次に Excelを学ぶ機会が多く，Excelの画面操作，関数 の使い方や種類に慣れていること，(2) 文系学部生 ...
17
動関数が反対称の性質を持つことになる.(2.4.4) の行列式を Slater 行列式とよぶ.(2.4.1) の波動方程式を解くためには,(2.4.4) の Slater 行列式を用いた波動関数 Φ e に関するエネルギー汎関数の最小値を求めればよい. その際に, 一電子軌道に関する極値 ( つまり
... 第一原理計算による XANES 計算の基礎 東京大学・生産技術研究所・溝口照康 本発表では電子状態理論に基づく XANES の理論計算法を説明する．XANES のスペクトル形状 は伝導帯の電子構造を反映しており， 伝導帯の電子構造を適切に計算することができれば XANES に対応するスペクトルも計算できるはずである．近年ではグラフィカルユーザーインターフェー ...
34
行列計算専用大規模集積回路の開発
... ルータを通じてそれぞれのプロセッサにデータを配ることになるが、プロセッサ はデータを受け取るまでの間処理を待たされる。行列を分けた行ベクトルがそれぞ れのプロセッサのメモリに記憶されているとしても、列ベクトルを送るにはプロセッ サの数だけ繰り返さなくてはならない。この場合プロセッサ数を増して計算が速く ...
82
1M5-3 行列因子分解による遺伝子データからの潜在的因子の抽出
... 1. はじめに 医学生物系のデータを扱う場合，遺伝子（タンパク質）とその 役割（機能）の関係は重要で，従来から広く収集活用されている． 遺伝子の機能などの情報は日々発見され，データベースが更 新される．記述に用いる用語については，オントロジーの研究 から次第に整備がすすんでいる．例えば遺伝子機能について は Gene Ontology (GO) [GO 00, GO13] が整備されており，著 ...
3
Excelを用いた行列演算
... ① 単一方程式モデル（一般に最小二乗法） 経済のある特定部分のメカニズムをそれ以外の部分から切り離して研究する場合 例） ある特定財についての生産関数の推定や需要予測など 変数の名称：「独立変数」⇒「従属変数」、「説明変数」⇒「目的変数」、「制御変数」⇒「応答変数」etc. ...
30
空間的自己相関モデルにおける距離行列のべき数の推定と地価決定要因分析への応用
... 空間要因は ρ , α ともに全期間で有意であり、住宅地では バブル期をはさんで空間的な影響が一貫して存在し、FIT の上昇や寄与率からも地価の決定に大きく影響しているこ とがわかる。したがって、住宅地では常に周囲の地点の地 価を参照して地価が決定されるというメカニズムが働いて いるといえる。また、α の推定値は 1.5～2.3 の範囲で安定 ...
6
疎な転置推移確率行列
... サブチーム課題のやり方 (計算科学☆演習 B) III Moodle への課題提出方法パターン 1: 完成したファイルは , 同一の ファイルを , サブチームメンバーごとに Moodle に提出します . Windows をログオフしてなくても , Moodle だけログイン / ログアウ ...
11
今回の内容 CUDA 付属のライブラリ cublas 行列 ベクトル積, 行列 行列積 cusperse 行列格納形式 cufft 余弦波の FFT curand モンテカルロ法による円周率計算 Thrust 913
... 複素数から複素数へのFFTを実行 cufftExecC2C(plan, dev, dev,CUFFT_FORWARD); 結果の表示（実部と虚部の取り出し） cuCrealf(re_host[i]) 複素数の実部を取り出す cuCimagf(re_host[i]) 複素数の虚部を取り出す ...
58
本日の講義内容 固有値 ( 線形代数 ) と応用問題 振動問題 ネットワーク定常問題 固有値計算アルゴリズム 密行列 べき乗法 ヤコビ法 ハウスホルダー三重対角 + 分割統治法 + 逆変換 疎行列 ランチョス法 ヤコビ デビッドソン法 その他 固有値計算ソフトウェア ScaLAPACK EigenE
... ハウスホルダー逆変換 • ハウスホルダー変換自身はそれ自身が逆変換でもあるので、 作用させた逆順に作用していけばよい。 • 近代的なアルゴリズムでは複数のハウスホルダー変換をまとめ てブロック化する。 ...
64
スパース行列のガウス消去における最適ピボッテング・アルゴリズムとそのフォトランプログラム: University of the Ryukyus Repository
... Title スパース行列のガウス消去における最適ピボッテング・ アルゴリズムとそのフォトランプログラム Author(s) 喜屋武, 盛基; 白川, 功; 尾崎, 弘 Citation 琉球大学理工学部紀要. 工学篇 = Bulletin of Science & Engineering Division, University of the Ryukyus. Engineering(7): 137-150 ...
15
に対して 例 2: に対して 逆行列は常に存在するとは限らない 逆行列が存在する行列を正則行列 (regular matrix) という 正則である 逆行列が存在する 一般に 正則行列 A の逆行列 A -1 も正則であり (A -1 ) -1 =A が成り立つ また 2 つの正則行列 A B の積
... 2) 逆行列を求めるプログラムは INVMTX という名前のサブルーチン（C の場合は invmtx と いう名前の関数）として作成せよ。サブルーチン（関数）の中では、別に作成したサブ ルーチン（関数）を参照してもよい。 3) 主プログラム中ではｎの値は固定でよいが、サブルーチン（関数）は任意のｎに対して ...
14