Method buildMap

DEVELOPER/core/Subdomain.cpp:882–930

int
Subdomain::buildMap(void)
{
  if (mapBuilt == false) {
	// determine the mapping between local dof and subdomain ana dof
        int numDOF = this->getNumDOF();
	if (map == 0) 
	  map = new ID(numDOF);
	if (map->Size() != numDOF) {
	  delete map;
	  map = new ID(numDOF);
	}

	//	int numExt = theAnalysis->getNumExternalEqn();
	int numInt = theAnalysis->getNumInternalEqn();

	const ID &theExtNodes = this->getExternalNodes();
	int numExtNodes = theExtNodes.Size();
	int locInMap =0;
	for (int i=0; i<numExtNodes; i++) {
	  Node *nodePtr = this->getNode(theExtNodes(i));
	  int numNodeDOF = nodePtr->getNumberDOF();
	  DOF_Group *theDOF = nodePtr->getDOF_GroupPtr();
	  const ID &theLocalID = theDOF->getID();
	  for (int j=0; j<numNodeDOF; j++){
	    int locInSubdomainExt = theLocalID(j)-numInt;
	    (*map)(locInMap)=locInSubdomainExt;
	    locInMap++;
	  }
	}
	mapBuilt = true;

	if (mappedVect == 0) 
	  mappedVect = new Vector(numDOF);
	if (mappedVect->Size() != numDOF) {
	  delete mappedVect;
	  mappedVect = new Vector(numDOF);
	}

	if (mappedMatrix == 0) 
	  mappedMatrix = new Matrix(numDOF,numDOF);
	if (mappedMatrix->noRows() != numDOF) {
	  delete mappedMatrix;
	  mappedMatrix = new Matrix(numDOF,numDOF);
	}
  }
  
  return 0;
}


DomainDecompositionAnalysis *