Method setSize

DEVELOPER/system/SkylineSPD.cpp:112–261

}

int 
SkylineSPD::setSize(Graph &theGraph)
{
    int oldSize = size;
    int result = 0;
    size = theGraph.getNumVertex();

    // check we have enough space in iDiagLoc and iLastCol
    // if not delete old and create new
    if (size > Bsize) { 
      if (iDiagLoc != 0) delete [] iDiagLoc;
      if (RowTop != 0) delete [] RowTop;
      if (topRowPtr != 0) delete [] topRowPtr;
      if (invD != 0) delete [] invD;
      // if (topRowPtr != 0) free((void *)topRowPtr);
      
      iDiagLoc = new int[size];
      RowTop = new int[size];
      invD = new double[size]; 
      //      topRowPtr = new double *[size] ;
      topRowPtr = (double **)malloc(size *sizeof(double *));


      if (iDiagLoc == 0 || RowTop == 0 || topRowPtr == 0 || invD == 0) {
	opserr << "WARNING SkylineSPD::setSize() : ";
	opserr << " - ran out of memory for iDiagLoc\n";
	size = 0; Asize = 0;
	result = -1;
      }
    }
    
    // zero out iDiagLoc 
    for (int i=0; i<size; i++) {
      iDiagLoc[i] = 0;
    }
    
    // now we go through the vertices to find the height of each col and
    // width of each row from the connectivity information.
    
    Vertex *vertexPtr;
    VertexIter &theVertices = theGraph.getVertices();
    
    while ((vertexPtr = theVertices()) != 0) {
      int vertexNum = vertexPtr->getTag();
      const ID &theAdjacency = vertexPtr->getAdjacency();
      int iiDiagLoc = iDiagLoc[vertexNum];
      int *iiDiagLocPtr = &(iDiagLoc[vertexNum]);
      
      for (int i=0; i<theAdjacency.Size(); i++) {
	int otherNum = theAdjacency(i);
	int diff = vertexNum-otherNum;
	if (diff > 0) {
	  if (iiDiagLoc < diff) {
	    iiDiagLoc = diff;
	    *iiDiagLocPtr = diff;
	  }
	} 
      }