changeset 6614:c5d871a4f1ef octave-forge

optim/inst/leasqr.m: fix for the one-residue case
author i7tiol
date Tue, 26 Jan 2010 13:34:10 +0000
parents 3f9e79a4297b
children 07e143f5c587
files main/optim/inst/leasqr.m
diffstat 1 files changed, 1 insertions(+), 1 deletions(-) [+]
line wrap: on
line diff
--- a/main/optim/inst/leasqr.m	Mon Jan 25 18:00:11 2010 +0000
+++ b/main/optim/inst/leasqr.m	Tue Jan 26 13:34:10 2010 +0000
@@ -274,7 +274,7 @@
     sgoal=(1-stol)*sprev;
     msk = dp ~= 0;
     prt(:, msk) = prt(:, msk) .* wt(:, ones (1, sum (msk)));
-    nrm(msk) = sumsq (prt(:, msk));
+    nrm(msk) = sumsq (prt(:, msk), 1);
     msk = nrm > 0;
     nrm(msk) = 1 ./ sqrt (nrm(msk));
     prt = prt .* nrm(ones (1, m), :);