2007年01月22日
カテゴリ： カテゴリ未分類
プログラミングのための線形代数という本を読んでいます。
http://wiki.fdiary.net/lacs/
3.8 LU分割の手順(2)例外が生じる場合
3.8.1 並べ替えが必要になる状況
LU分解の途中で1/0ができてしまい計算できないケースを扱います。
◆理屈
今までも扱った pivoting という手法を使います。行を入れ替え何事もなかったかのようにLU分解の手順を続けていきます。
◆実装
実装の上では pivoting は間接参照で扱い読みかえるだけにします。後は0かどうかの判定をする際に気を使うかあるいは一番よさそうな行を毎回選ぶという手がよさそうです。
http://wiki.fdiary.net/lacs/
3.8 LU分割の手順(2)例外が生じる場合
3.8.1 並べ替えが必要になる状況
LU分解の途中で1/0ができてしまい計算できないケースを扱います。
◆理屈
今までも扱った pivoting という手法を使います。行を入れ替え何事もなかったかのようにLU分解の手順を続けていきます。
◆実装
実装の上では pivoting は間接参照で扱い読みかえるだけにします。後は0かどうかの判定をする際に気を使うかあるいは一番よさそうな行を毎回選ぶという手がよさそうです。
お気に入りの記事を「いいね！」で応援しよう
【毎日開催】
15記事にいいね！で1ポイント
PR
キーワードサーチ
▼キーワード検索
大鳥居.
人生ゲーム
コメント新着
© Rakuten Group, Inc.
