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法で <literal>M.X = b</literal> を解く．

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

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

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

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

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

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

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

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

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

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

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

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

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

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