Mercurial > hg > octave-lyh
comparison liboctave/CmplxQR.cc @ 7554:40574114c514
implement Cholesky factorization updating
| author | Jaroslav Hajek <highegg@gmail.com> |
|---|---|
| date | Tue, 04 Mar 2008 22:25:50 -0500 |
| parents | 56be6f31dd4e |
| children | 07522d7dcdf8 |
comparison
equal
deleted
inserted
replaced
| 7553:56be6f31dd4e | 7554:40574114c514 |
|---|---|
| 273 } | 273 } |
| 274 | 274 |
| 275 void | 275 void |
| 276 ComplexQR::economize (void) | 276 ComplexQR::economize (void) |
| 277 { | 277 { |
| 278 idx_vector c (':'), i (Range (1, r.columns ())); | |
| 279 q = ComplexMatrix (q.index (c, i)); | |
| 280 r = ComplexMatrix (r.index (i, c)); | |
| 281 #if 0 | |
| 282 octave_idx_type r_nc = r.columns (); | 278 octave_idx_type r_nc = r.columns (); |
| 283 | 279 |
| 284 if (r.rows () > r_nc) | 280 if (r.rows () > r_nc) |
| 285 { | 281 { |
| 286 q.resize (q.rows (), r_nc); | 282 q.resize (q.rows (), r_nc); |
| 287 r.resize (r_nc, r_nc); | 283 r.resize (r_nc, r_nc); |
| 288 } | 284 } |
| 289 #endif | |
| 290 } | 285 } |
| 291 | 286 |
| 292 /* | 287 /* |
| 293 ;;; Local Variables: *** | 288 ;;; Local Variables: *** |
| 294 ;;; mode: C++ *** | 289 ;;; mode: C++ *** |