Mercurial > octave-antonio

comparison scripts/linear-algebra/onenormest.m @ 8507:cadc73247d65

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
style fixes
author John W. Eaton <jwe@octave.org>
date Tue, 13 Jan 2009 14:08:36 -0500
parents bc982528de11
children e07e93c04080
comparison
equal deleted inserted replaced
8506:bc982528de11 8507:cadc73247d65
187 numpar = sum (partest); 187 numpar = sum (partest);
188 replacements = 2*(rand (n,numpar) < 0.5) - 1; 188 replacements = 2*(rand (n,numpar) < 0.5) - 1;
189 S(:,partest) = replacements; 189 S(:,partest) = replacements;
190 endif 190 endif
191 ## Now test for parallel vectors within S. 191 ## Now test for parallel vectors within S.
192 partest = any ( (S' * S - eye (t)) == n ); 192 partest = any ((S' * S - eye (t)) == n);
193 if (any (partest)) 193 if (any (partest))
194 numpar = sum (partest); 194 numpar = sum (partest);
195 replacements = 2*(rand (n,numpar) < 0.5) - 1; 195 replacements = 2*(rand (n,numpar) < 0.5) - 1;
196 S(:,partest) = replacements; 196 S(:,partest) = replacements;
197 endif 197 endif