cholmod_solve, cholmode_l_solveは結構大変 どう大変か、というと こちらからもわかるように cholmod_dense * cholmod_solve (int sys, cholmod_factor *L, cholmod_dense *B, cholmod_common *Common) solveする式のパターンを決めるint sys (たとえばAx =…

前の記事の疎行列用クラスのラッパークラス(spinxform)のヘッダー定義と関数コードを眺めて、疎行列クラスがどのようになっているかを見てみる 疎行列クラスの基本 nxm行列の要素を全部並べると大きすぎる。疎であることを利用して、「番地」と「値」の組だ…

疎行列情報ファイル 例から見よう。３ｘ３行列に５つの非ゼロ要素がある複素数を要素とする行列。./Matrix/c.matというファイルだとする %%MatrixMarket matrix coordinate complex hermitian 3 3 5 1 1 1. 0. 3 1 2. -1. 2 2 1. 0. 3 2 3. 0. 3 3 42. 0. 以…

疎行列情報を記録したファイルを読み込んで、それをハンドリングする、というスタイルを基本にしたい ディレクトリ hogeDirを作り、その下に 以下に示す Makefileを置く。このMakefileは以下に示すディレクトリ構造・ファイル配置で動くように作ってある デ…

Ubuntuのパッケージサイトはこちら Ubuntuのパッケージ管理システムについてはこちら コマンドラインのそれはapt-get Ubuntuのバージョンごとに提供ライブラリをわけて置いてある(こちらのパッケージ一覧に目を通す) trustyを使ったので、trustyのリンクに行…

Suitesparseは疎行列を扱う色々なアルゴリズムの集合体で、テキサス大の先生のサイトが提供元(こちら) cholmodはその一部 疎行列の扱いにあたって、グラフ理論を用いた疎行列の密化をし、そのうえと行列計算を組み合わせる仕組み グラフ理論・グラフの密化の…