Changeset 34ab5f6 in git
- Timestamp:
- Dec 2, 2010, 4:06:01 PM (13 years ago)
- Branches:
- (u'spielwiese', '0d6b7fcd9813a1ca1ed4220cfa2b104b97a0a003')
- Children:
- be1e92ee0c87c9a2bdae2aa22601cba0ac97c857
- Parents:
- e81572c65614d207af363ed0da23be28de8043dd
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
kernel/sparsmat.cc
re81572c r34ab5f6 587 587 588 588 /* 589 * the Bareiss elimination:590 * - with x unreduced last rows, pivots from here are not allowed591 * - the method will finish for number of unreduced columns < y592 */593 void sparse_mat::smBareiss(int x, int y)594 {595 if ((x > 0) && (x < nrows))596 {597 tored -= x;598 this->smToredElim();599 }600 if (y < 1) y = 1;601 if (act <= y)602 {603 this->smCopToRes();604 return;605 }606 normalize = this->smCheckNormalize();607 if (normalize) this->smNormalize();608 this->smPivot();609 this->smSelectPR();610 this->smElim();611 crd++;612 this->smColToRow();613 act--;614 this->smRowToCol();615 this->smZeroElim();616 if (tored != nrows)617 this->smToredElim();618 if (act < y)619 {620 this->smCopToRes();621 return;622 }623 loop624 {625 if (normalize) this->smNormalize();626 this->smPivot();627 oldpiv = piv;628 this->smSelectPR();629 this->smElim();630 crd++;631 this->smColToRow();632 act--;633 this->smRowToCol();634 this->smZeroElim();635 if (tored != nrows)636 this->smToredElim();637 if (act < y)638 {639 if (TEST_OPT_PROT) PrintS(".\n");640 this->smCopToRes();641 return;642 }643 }644 }645 646 /*647 589 * the new Bareiss elimination: 648 590 * - with x unreduced last rows, pivots from here are not allowed
Note: See TracChangeset
for help on using the changeset viewer.