エントリーの編集
![loading...](https://b.st-hatena.com/bdefb8944296a0957e54cebcfefc25c4dcff9f5f/images/v4/public/common/loading@2x.gif)
エントリーの編集は全ユーザーに共通の機能です。
必ずガイドラインを一読の上ご利用ください。
記事へのコメント0件
- 注目コメント
- 新着コメント
このエントリーにコメントしてみましょう。
注目コメント算出アルゴリズムの一部にLINEヤフー株式会社の「建設的コメント順位付けモデルAPI」を使用しています
![アプリのスクリーンショット](https://b.st-hatena.com/bdefb8944296a0957e54cebcfefc25c4dcff9f5f/images/v4/public/entry/app-screenshot.png)
- バナー広告なし
- ミュート機能あり
- ダークモード搭載
関連記事
数値計算とかの備忘録(仮)
とあるシミュレーションの問題では,シフト方程式と呼ばれる線形方程式 という,係数行列Aの対角成分が... とあるシミュレーションの問題では,シフト方程式と呼ばれる線形方程式 という,係数行列Aの対角成分がシフト量σだけ変化した方程式を解く必要が出てくるそうだ.Aは大規模で疎なものを対象としており,σは複数個現れるそうだ.これを解く方法としては,σだけ変化した係数行列をつくり,それをKrylov部分空間法なんかで解いてあげれば良いのだが,σの数だけ係数行列を新たに生成するのはあまり効率がよろしくない. そこで,シフト方程式の係数行列を直接生成せずに,(シード方程式と呼ばれている)を解く際に,ついでにシフト方程式も解いてあげて,全体の計算量を節約しよう,という解法が生み出された.さまざまな解法が生み出されているが,その1つに,BiCG法を拡張することで得られるShifted BiCG法がある. 論文に擬似コードが載っていたので,matlabで実装することにした.コード例はこんな感じ.シフト量σは複