Mercurial > octave-nkf

comparison scripts/polynomial/polyfit.m @ 17307:4448cc742880

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Use '!' not operator rather than '~' in core m-files. * scripts/general/interp1.m, scripts/image/rgb2ind.m, scripts/plot/findobj.m, scripts/polynomial/polyfit.m, test/build-sparse-tests.sh, test/classes/@Snork/end.m, test/classes/@Snork/mpower.m, test/classes/@Snork/ne.m, test/classes/@Snork/power.m, test/classes/@Snork/subsasgn.m, test/classes/@Snork/subsref.m, test/classes/@Spork/loadobj.m, test/classes/classes.tst: Use '!' not operator rather than '~' in core m-files.
author Rik <rik@octave.org>
date Wed, 21 Aug 2013 17:04:44 -0700
parents 7babcdb9bc13
children b327c19b6b16
comparison
equal deleted inserted replaced
17306:09543e9c8f40 17307:4448cc742880
113 ## Solve by QR decomposition. 113 ## Solve by QR decomposition.
114 [q, r, k] = qr (v(:, polymask), 0); 114 [q, r, k] = qr (v(:, polymask), 0);
115 p = r \ (q' * y); 115 p = r \ (q' * y);
116 p(k) = p; 116 p(k) = p;
117 117
118 if (n ~= m) 118 if (n != m)
119 q = p; p = zeros (n+1, 1); 119 q = p; p = zeros (n+1, 1);
120 p(polymask) = q; 120 p(polymask) = q;
121 endif 121 endif
122 122
123 if (nargout > 1) 123 if (nargout > 1)
138 C = cholinv (r.'*r)(k, k); 138 C = cholinv (r.'*r)(k, k);
139 catch 139 catch
140 C = NaN * ones (m+1, m+1); 140 C = NaN * ones (m+1, m+1);
141 end_try_catch 141 end_try_catch
142 142
143 if (n ~= m) 143 if (n != m)
144 ## fill matrices if required 144 ## fill matrices if required
145 s.X(:, ~polymask) = 0; 145 s.X(:, !polymask) = 0;
146 s.R = zeros (n+1, n+1); s.R(polymask, polymask) = r; 146 s.R = zeros (n+1, n+1); s.R(polymask, polymask) = r;
147 s.C = zeros (n+1, n+1); s.C(polymask, polymask) = C; 147 s.C = zeros (n+1, n+1); s.C(polymask, polymask) = C;
148 else 148 else
149 s.R = r; 149 s.R = r;
150 s.C = C; 150 s.C = C;