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)

説明 :

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

コマンドリスト

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)

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