Mercurial > hg > octave-nkf
diff scripts/linear-algebra/krylov.m @ 5775:ace8d8d26933
[project @ 2006-04-24 19:13:06 by jwe]
author | jwe |
---|---|
date | Mon, 24 Apr 2006 19:13:11 +0000 |
parents | dbdba35033a6 |
children | b3478d7a0486 |
line wrap: on
line diff
--- a/scripts/linear-algebra/krylov.m +++ b/scripts/linear-algebra/krylov.m @@ -122,7 +122,7 @@ nv = columns (V); if (jj != nv) [V(:,jj), V(:,nv)] = swap (V(:,jj), V(:,nv)); - ## XXX FIXME XXX -- H columns should be swapped too. Not done + ## FIXME -- H columns should be swapped too. Not done ## since Block Hessenberg structure is lost anyway. endif V = V(:,1:(nv-1)); @@ -156,7 +156,7 @@ # reduce V per the reflection V(idx,:) = V(idx,:) - av*hv*(hv' * V(idx,:)); if(iter > 1) - ## XXX FIXME XXX -- not done correctly for block case + ## FIXME -- not done correctly for block case H(nu,nu-1) = V(pivot_vec(nu),jj); endif