Method eigen

SRC/analysis/analysis/StaticDomainDecompositionAnalysis.cpp:240–341

int 
StaticDomainDecompositionAnalysis::eigen(int numMode, bool generalized, bool findSmallest)
{
  int result = 0;
  Domain *the_Domain = this->getDomainPtr();

  if (theEigenSOE == 0) {
    opserr << "StaticDomainDecompositionAnalysis::eigen() - no eigen solver has been set\n";
    return -1;
  }

  // check for change in Domain since last step. As a change can
  // occur in a commit() in a domaindecomp with load balancing
  // this must now be inside the loop

  int stamp = the_Domain->hasDomainChanged();
  if (stamp != domainStamp) {
    domainStamp = stamp;
    result = this->domainChanged();
    if (result < 0) {
      opserr << "StaticDomainDecompositionAnalysis::eigen() - domainChanged failed";
      return -1;
    }	
  }

    //
    // zero A and M
    //

    theEigenSOE->zeroA();
    theEigenSOE->zeroM();

    //
    // form K
    //

    FE_EleIter &theEles = theAnalysisModel->getFEs();    
    FE_Element *elePtr;

    while((elePtr = theEles()) != 0) {
      elePtr->zeroTangent();
      elePtr->addKtToTang(1.0);
      if (theEigenSOE->addA(elePtr->getTangent(0), elePtr->getID()) < 0) {
	opserr << "WARNING StaticAnalysis::eigen() -";
	opserr << " failed in addA for ID " << elePtr->getID();	    
	result = -2;
      }
    }

    //
    // if generalized is true, form M
    //

    if (generalized == true) {
      int result = 0;
      FE_EleIter &theEles2 = theAnalysisModel->getFEs();    
      while((elePtr = theEles2()) != 0) {     
	elePtr->zeroTangent();