Method performEigenAnalysis

SRC/reliability/analysis/misc/MatrixOperations.cpp:618–769

#endif

int MatrixOperations::performEigenAnalysis(int pBeginMod, int pEndMod)
{

  beginMod=pBeginMod;
  endMod = pEndMod;
 
  Matrix * A =  theMatrix;

//  opserr<< "hessian:"<< * theMatrix<<endln;

  if (A ==0) {opserr<<"error, Hessian does not exist in MatrixOperator::performEigenAnalysis !"<<endln; exit(-1);}

  // Number of equations
  int n = A->noRows();
  if (n != A->noCols()) {
	opserr<<"MatrixOperations::performEigenAnalysis wrong. m!=n"<<endln;
	exit(-1);
  }

  numModes = endMod - beginMod +1; //Q
  // Check for quick return
  if (numModes < 1) {
    numModes = 0;
    return 0;
  }

  // Simple check
  if (numModes > n)
    numModes = n;

  // Allocate storage for eigenvalues
  if (eigenvalue != 0)
    delete [] eigenvalue;
  eigenvalue = new double [n];

  // Real work array (see LAPACK dsbevx subroutine documentation)
  double *work = new double [7*n];

  // Integer work array (see LAPACK dsbevx subroutine documentation)
  int *iwork = new int [5*n];

  // Leading dimension of eigenvectors
  int ldz = n;
  setSizeOfEigenVector(ldz); // alloc memory here.

  // Allocate storage for eigenvectors
  if (eigenvector != 0)
    delete [] eigenvector;
  eigenvector = new double [ldz*numModes];

  // Number of superdiagonals
  int kd = n-1;

  // Matrix data

  
  //double *ab =  theSOE->A;  //
  double *ab = new double [n*n];