Method allocateMemory

src/engine/ConstraintMatrixAnalyzer.cpp:123–156

}

void ConstraintMatrixAnalyzer::allocateMemory()
{
    // Initialize the row and column headers
    _rowHeaders = new unsigned[_m];
    _columnHeaders = new unsigned[_n];
    _rowHeadersInverse = new unsigned[_m];
    _columnHeadersInverse = new unsigned[_n];

    for ( unsigned i = 0; i < _m; ++i )
    {
        _rowHeaders[i] = i;
        _rowHeadersInverse[i] = i;
    }

    for ( unsigned i = 0; i < _n; ++i )
    {
        _columnHeaders[i] = i;
        _columnHeadersInverse[i] = i;
    }

    // Initialize the row and column counters
    _numRowElements = new unsigned[_m];
    _numColumnElements = new unsigned[_n];

    for ( unsigned i = 0; i < _m; ++i )
        _numRowElements[i] = _A.getRow( i )->getNnz();

    for ( unsigned i = 0; i < _n; ++i )
        _numColumnElements[i] = _At.getRow( i )->getNnz();

    // Work memory
    _workRow = new double[_n];
    _workRow2 = new double[_n];
}

void ConstraintMatrixAnalyzer::gaussianElimination()
{