Mercurial > hg > octave-nkf

comparison liboctave/sparse-dmsolve.cc @ 5775:ace8d8d26933

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
[project @ 2006-04-24 19:13:06 by jwe]
author jwe
date Mon, 24 Apr 2006 19:13:11 +0000
parents 2fe20065a545
children 6b9cec830d72
comparison
equal deleted inserted replaced
5774:e7af222e98b0 5775:ace8d8d26933
418 btmp.insert (ctmp - m * mtmp, 0, 0); 418 btmp.insert (ctmp - m * mtmp, 0, 0);
419 } 419 }
420 } 420 }
421 421
422 // Structurally non-singular blocks 422 // Structurally non-singular blocks
423 // XXX FIXME XXX Should use fine Dulmange-Mendelsohn decomposition here. 423 // FIXME Should use fine Dulmange-Mendelsohn decomposition here.
424 if (dm->rr [1] < dm->rr [2] && dm->cc [2] < dm->cc [3] && 424 if (dm->rr [1] < dm->rr [2] && dm->cc [2] < dm->cc [3] &&
425 !info && !error_state) 425 !info && !error_state)
426 { 426 {
427 ST m = dmsolve_extract (a, pinv, q, dm->rr [1], dm->rr [2], 427 ST m = dmsolve_extract (a, pinv, q, dm->rr [1], dm->rr [2],
428 dm->cc [2], dm->cc [3], nnz_remaining, false); 428 dm->cc [2], dm->cc [3], nnz_remaining, false);