Ignore:
Timestamp:
Dec 12, 2013 10:40:41 AM (6 years ago)
Author:
forrest
Message:

changes for parallel and idiot

File:
1 edited

Legend:

Unmodified
Added
Removed
  • stable/1.15/Clp/src/CoinAbcBaseFactorization1.cpp

    r1910 r2006  
    14831483                 eArea,saveColumn,startColumnU,tempColumnCount,
    14841484                 elementU,numberInColumn,indexRowU);
     1485    assert (tempColumnCount[0]>=0);
    14851486    for (int iTry=numberZeroOther-1;iTry>=0;iTry--) {
    14861487      pivotWhile(0,numberInPivotRow,numberInPivotColumn,lengthArea,giveUp,
    14871488                 eArea,multipliersL);
     1489    assert (tempColumnCount[0]>=0);
    14881490      // find new pivot row
    14891491      int jPivotColumn=others[iTry];
     
    15061508        // put away last
    15071509        afterPivot(pivotRow,pivotColumn);
     1510        if (reinterpret_cast<unsigned int *>(elementUColumnPlusAddress_+lastEntryByColumnUPlus_)>aBits) {
     1511#ifndef NDEBUG
     1512          printf("? dense\n");
     1513#endif
     1514          return -99;
     1515        }
    15081516        // new ones
    15091517        pivotRow=indexL[jPivotRow];
     
    15271535        lastRow[pivotRow] = -2;
    15281536        numberInRow[pivotRow] = 0;
     1537    assert (tempColumnCount[0]>=0);
    15291538        //store column in L, compress in U and take column out
    15301539        CoinBigIndex l = lengthL_;
     
    15391548        // move to L and move last
    15401549        numberInPivotColumn--;
     1550    assert (tempColumnCount[0]>=0);
    15411551        if (jPivotRow!=numberInPivotColumn) {
    15421552          // swap
     
    15561566          jPivotRow=numberInPivotColumn;
    15571567        }
     1568    assert (tempColumnCount[0]>=0);
    15581569        numberInPivotRow--;
    15591570        CoinFactorizationDouble pivotElement = area[numberInPivotColumn];
     
    15681579        totalElements_-=numberInColumn[pivotColumn];
    15691580        mark[pivotColumn]=0;
     1581    assert (tempColumnCount[0]>=0);
    15701582        for (int jRow=0;jRow<numberInPivotColumn+1;jRow++) {
    15711583          CoinFactorizationDouble temp=areaLast[jRow];
     
    15781590          CoinFactorizationDouble * COIN_RESTRICT area =eArea+jColumn*lengthArea;
    15791591          CoinFactorizationDouble thisPivotValue = area[numberInPivotColumn];
     1592    assert (tempColumnCount[0]>=0);
    15801593          area[numberInPivotColumn]=thisPivotValue;
    15811594          if (fabs(thisPivotValue)>tolerance) {
     
    16011614          }
    16021615        }
     1616    assert (tempColumnCount[0]>=0);
    16031617        // put away last (but remember count is one out)
    16041618        CoinAbcMemcpy(elementL+lSave,multipliersL,numberInPivotColumn+1);
     
    16311645      }
    16321646    }
     1647    assert (tempColumnCount[0]>=0);
    16331648#if 0
    16341649    // zero out extra bits
Note: See TracChangeset for help on using the changeset viewer.