gf_linsolve

概要

X = gf_linsolve('gmres', spmat M, vec b[, int restart][, precond P][,'noisy'][,'res', r][,'maxiter', n])
X = gf_linsolve('cg', spmat M, vec b [, precond P][,'noisy'][,'res', r][,'maxiter', n])
X = gf_linsolve('bicgstab', spmat M, vec b [, precond P][,'noisy'][,'res', r][,'maxiter', n])
{U, cond} = gf_linsolve('lu', spmat M, vec b)
{U, cond} = gf_linsolve('superlu', spmat M, vec b)
{U, cond} = gf_linsolve('mumps', spmat M, vec b, ... ['sym'])

説明 :

各種線形システムソルバー．

コマンドリスト

X = gf_linsolve('gmres', spmat M, vec b[, int restart][, precond P][,'noisy'][,'res', r][,'maxiter', n])

GMRES法で M.X = b を解く．

任意で前処理として P を使用します．restartパラメーターのデフォルト値は50です．

X = gf_linsolve('cg', spmat M, vec b [, precond P][,'noisy'][,'res', r][,'maxiter', n])

共役勾配法で M.X = b を解きます．

オプションで前処理 P を使用します．

X = gf_linsolve('bicgstab', spmat M, vec b [, precond P][,'noisy'][,'res', r][,'maxiter', n])

双共役勾配安定化法で M.X = b を解く．

オプションで前処理 P を使用します．

{U, cond} = gf_linsolve('lu', spmat M, vec b)

gf_linsolve('superlu',...) のエイリアス

{U, cond} = gf_linsolve('superlu', spmat M, vec b)

M.U = b を解くには，SuperLUソルバ(疎LU分解)を適用します．

条件数推定値 cond は，解 U とともに返されます．

{U, cond} = gf_linsolve('mumps', spmat M, vec b, ... ['sym'])

MUMPSソルバーを使用して， M.U = b を解きます．

右辺 b は、複数の右辺を一度に解くために、複数の列を持つ行列にすることもできます。

sym オプションを指定すると、対称バージョンのMUMPSソルバーが使用されます。