Method eigen

SRC/interpreter/OpenSeesCommands.cpp:264–389

}

int
OpenSeesCommands::eigen(int typeSolver, double shift,
			bool generalizedAlgo, bool findSmallest,
            EigenSOE *providedEigenSOE)
{
    //
    // create a transient analysis if no analysis exists
    //
    bool newanalysis = false;
    if (theStaticAnalysis == 0 && theTransientAnalysis == 0) {
	if (theAnalysisModel == 0)
	    theAnalysisModel = new AnalysisModel();
	if (theTest == 0)
	    theTest = new CTestNormUnbalance(1.0e-6,25,0);
	if (theAlgorithm == 0) {
	    theAlgorithm = new NewtonRaphson(*theTest);
	}
	if (theHandler == 0) {
	    theHandler = new TransformationConstraintHandler();
	}
	if (theNumberer == 0) {
	    RCM *theRCM = new RCM(false);
	    theNumberer = new DOF_Numberer(*theRCM);
	}
	if (theTransientIntegrator == 0) {
        setIntegrator(new Newmark(0.5,0.25), true);
	}
	if (theSOE == 0) {
	    ProfileSPDLinSolver *theSolver;
	    theSolver = new ProfileSPDLinDirectSolver();
	    theSOE = new ProfileSPDLinSOE(*theSolver);
	}

	theTransientAnalysis = new DirectIntegrationAnalysis(*theDomain,
							     *theHandler,
							     *theNumberer,
							     *theAnalysisModel,
							     *theAlgorithm,
							     *theSOE,
							     *theTransientIntegrator,
							     theTest);
	newanalysis = true;
    }

    //
    // create or assign an eigen system and solver
    //
    bool eigenSOEUpdated = false;

    if (providedEigenSOE != 0) {
        theEigenSOE = providedEigenSOE;
        eigenSOEUpdated = true;
    } else {
        if (theEigenSOE != 0) {
            if (theEigenSOE->getClassTag() != typeSolver) {
                //	delete theEigenSOE;
                theEigenSOE = 0;
            }