Changeset 34ab5f6 in git for kernel


Ignore:
Timestamp:
Dec 2, 2010, 4:06:01 PM (13 years ago)
Author:
Hans Schoenemann <hannes@…>
Branches:
(u'fieker-DuVal', '117eb8c30fc9e991c4decca4832b1d19036c4c65')(u'spielwiese', 'c5facdfddea2addfd91babd8b9019161dea4b695')
Children:
be1e92ee0c87c9a2bdae2aa22601cba0ac97c857
Parents:
e81572c65614d207af363ed0da23be28de8043dd
Message:
code cleanup: smBareiss

git-svn-id: file:///usr/local/Singular/svn/trunk@13696 2c84dea3-7e68-4137-9b89-c4e89433aadc
File:
1 edited

Legend:

Unmodified
Added
Removed
  • kernel/sparsmat.cc

    re81572c r34ab5f6  
    587587
    588588/*
    589 * the Bareiss elimination:
    590 *   - with x unreduced last rows, pivots from here are not allowed
    591 *   - the method will finish for number of unreduced columns < y
    592 */
    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   loop
    624   {
    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 /*
    647589* the new Bareiss elimination:
    648590*   - with x unreduced last rows, pivots from here are not allowed
Note: See TracChangeset for help on using the changeset viewer.