Method setDomain

DEVELOPER/recorder/SumElementForcesRecorder.cpp:148–221

int 
SumElementForcesRecorder::setDomain(Domain &theDom)
{
  theDomain = &theDom;

  // delete old data
  if (theElements != 0)
    delete [] theElements;
  if (data != 0)
    delete data;

  // set numEle
  numEle = eleID.Size();
  if (numEle == 0) {
    opserr << "WARNING SumElementForcesRecorder::initialize() - no elements tags passed in input!\n";
    return 0; 
  }

  // create theElements, an array of pointers to elements
  theElements = new Element *[numEle];
  if (theElements == 0) {
    opserr << "WARNING SumElementForcesRecorder::initialize() - out of memory\n";
    numEle = 0; // set numEle = 0, in case record() still called
    return -1;
  }

  // zero theElements array
  for (int k=0; k<numEle; k++)
    theElements[k] = 0;

  //
  // loop over the list of elements, 
  //    if element exists add it's pointer o the array
  //    get its resisting force, check size to determine compatible with others
  // 

  int sizeArray = -1;

  for (int i=0; i<numEle; i++) {
    int eleTag = eleID(i);

    Element *theEle = theDomain->getElement(eleTag);
    
    if (theEle != 0) {

      const Vector &force = theEle->getResistingForce();
      int forceSize = force.Size();
      if (sizeArray == -1) {
	sizeArray = forceSize;
	theElements[i] = theEle;
      } else if (sizeArray != forceSize) {
	opserr << "WARNING: forces mismatch - element: " << eleTag << " will not be included\n";
      } else {
	theElements[i] = theEle;	
      }
    }
  }