Changeset 1903


Ignore:
Timestamp:
Jan 3, 2013 5:49:30 PM (7 years ago)
Author:
stefan
Message:

code from John to check that fillBasis works correctly (need to be enabled via define)

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/Clp/src/ClpFactorization.cpp

    r1732 r1903  
    15281528               numberElements += matrix->countBasis(pivotTemp + numberRowBasic,
    15291529                                                    numberColumnBasic);
     1530#ifndef NDEBUG
     1531//#define CHECK_CLEAN_BASIS
     1532#ifdef CHECK_CLEAN_BASIS
     1533               int saveNumberElements = numberElements;
     1534#endif
     1535#endif
    15301536               // Not needed for dense
    15311537               numberElements = 3 * numberBasic + 3 * numberElements + 20000;
     
    15451551               indexRowU = coinFactorizationB_->indices();
    15461552               startColumnU = coinFactorizationB_->starts();
     1553#ifdef CHECK_CLEAN_BASIS
     1554               for (int i=0;i<saveNumberElements;i++) {
     1555                 elementU[i]=0.0;
     1556                 indexRowU[i]=-1;
     1557               }
     1558               for (int i=0;i<numberRows;i++)
     1559                 startColumnU[i]=-1;
     1560#endif
    15471561#ifndef COIN_FAST_CODE
    15481562               double slackValue;
     
    15851599               else
    15861600                    numberElements = 0;
     1601#ifdef CHECK_CLEAN_BASIS
     1602               assert (!startColumnU[0]);
     1603               int lastStart=0;
     1604               for (int i=0;i<numberRows;i++) {
     1605                 assert (startColumnU[i+1]>lastStart);
     1606                 lastStart=startColumnU[i+1];
     1607               }
     1608               assert (lastStart==saveNumberElements);
     1609               for (int i=0;i<saveNumberElements;i++) {
     1610                 assert(elementU[i]);
     1611                 assert(indexRowU[i]>=0&&indexRowU[i]<numberRows);
     1612               }
     1613#endif
    15871614               coinFactorizationB_->preProcess ( );
    15881615               coinFactorizationB_->factor (  );
Note: See TracChangeset for help on using the changeset viewer.