Mercurial > hg > octave-nkf

comparison liboctave/dbleQR.cc @ 8806:c7864bb74914

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
avoid some GCC warnings
author John W. Eaton <jwe@octave.org>
date Wed, 18 Feb 2009 13:45:20 -0500
parents e9cb742df9eb
children 20dfb885f877
comparison
equal deleted inserted replaced
8805:065a05eb148a 8806:c7864bb74914
237 octave_idx_type k = q.columns (); 237 octave_idx_type k = q.columns ();
238 238
239 if (u.rows () == m && v.rows () == n && u.cols () == v.cols ()) 239 if (u.rows () == m && v.rows () == n && u.cols () == v.cols ())
240 { 240 {
241 OCTAVE_LOCAL_BUFFER (double, w, 2*k); 241 OCTAVE_LOCAL_BUFFER (double, w, 2*k);
242 for (octave_idx_type i = 0; i < u.cols (); i++) 242 for (volatile octave_idx_type i = 0; i < u.cols (); i++)
243 { 243 {
244 ColumnVector utmp = u.column (i), vtmp = v.column (i); 244 ColumnVector utmp = u.column (i), vtmp = v.column (i);
245 F77_XFCN (dqr1up, DQR1UP, (m, n, k, q.fortran_vec (), m, r.fortran_vec (), k, 245 F77_XFCN (dqr1up, DQR1UP, (m, n, k, q.fortran_vec (), m, r.fortran_vec (), k,
246 utmp.fortran_vec (), vtmp.fortran_vec (), w)); 246 utmp.fortran_vec (), vtmp.fortran_vec (), w));
247 } 247 }
313 { 313 {
314 r.resize (k, n + nj); 314 r.resize (k, n + nj);
315 } 315 }
316 316
317 OCTAVE_LOCAL_BUFFER (double, w, kmax); 317 OCTAVE_LOCAL_BUFFER (double, w, kmax);
318 for (octave_idx_type i = 0; i < js.length (); i++) 318 for (volatile octave_idx_type i = 0; i < js.length (); i++)
319 { 319 {
320 ColumnVector utmp = u.column (jsi(i)); 320 ColumnVector utmp = u.column (jsi(i));
321 F77_XFCN (dqrinc, DQRINC, (m, n + i, std::min (kmax, k + i), 321 F77_XFCN (dqrinc, DQRINC, (m, n + i, std::min (kmax, k + i),
322 q.fortran_vec (), q.rows (), 322 q.fortran_vec (), q.rows (),
323 r.fortran_vec (), r.rows (), js(i) + 1, 323 r.fortran_vec (), r.rows (), js(i) + 1,
372 else if (nj > 0 && (js(0) > n-1 || js(nj-1) < 0)) 372 else if (nj > 0 && (js(0) > n-1 || js(nj-1) < 0))
373 (*current_liboctave_error_handler) ("qrinsert: index out of range"); 373 (*current_liboctave_error_handler) ("qrinsert: index out of range");
374 else if (nj > 0) 374 else if (nj > 0)
375 { 375 {
376 OCTAVE_LOCAL_BUFFER (double, w, k); 376 OCTAVE_LOCAL_BUFFER (double, w, k);
377 for (octave_idx_type i = 0; i < js.length (); i++) 377 for (volatile octave_idx_type i = 0; i < js.length (); i++)
378 { 378 {
379 F77_XFCN (dqrdec, DQRDEC, (m, n - i, k == m ? k : k - i, 379 F77_XFCN (dqrdec, DQRDEC, (m, n - i, k == m ? k : k - i,
380 q.fortran_vec (), q.rows (), 380 q.fortran_vec (), q.rows (),
381 r.fortran_vec (), r.rows (), js(i) + 1, w)); 381 r.fortran_vec (), r.rows (), js(i) + 1, w));
382 } 382 }