Function ops_eigenAnalysis

SRC/interpreter/PythonAnalysisBuilder.cpp:1149–1255

}

PyObject *ops_eigenAnalysis(PyObject *self, PyObject *args)
{
    OPS_ResetCommandLine(PyTuple_Size(args), 0, args);

    // check inputs
    int numArgs = OPS_GetNumRemainingInputArgs();
    if(numArgs < 1) {
	PyErr_SetString(PyExc_RuntimeError,"ERROR eigen('type',numModes) or eigen('PythonSparse',numModes,dict)");
	return NULL;
    }

    // get parameters
    bool generalizedAlgo = true;// 0 - frequency/generalized (default),1 - standard, 2 - buckling
    int typeSolver = EigenSOE_TAGS_ArpackSOE;
    double shift = 0.0;
    bool findSmallest = true;
    int numEigen = 0;
    EigenSOE *providedEigenSOE = nullptr;

    // Check if first argument is "PythonSparse" - if so, expect numModes and dict
    std::string firstArg = OPS_GetString();
    if (firstArg == "PythonSparse" || firstArg == "PythonCompressedSparseEigen") {
        // PythonSparse eigen solver - requires numModes as second argument, dict as third
        // First get numModes
        int numData = 1;
        if(OPS_GetIntInput(&numData,&numEigen) < 0) {
            PyErr_SetString(PyExc_RuntimeError,"ERROR eigen('PythonSparse',numModes,dict) - invalid numModes");
            return NULL;
        }
        
        // Then get dict argument
        extern void *OPS_SparsePythonEigenSolver();
        void *eigenSOEPtr = OPS_SparsePythonEigenSolver();
        if (eigenSOEPtr == NULL) {
            return NULL;
        }
        providedEigenSOE = static_cast<EigenSOE *>(eigenSOEPtr);
        typeSolver = providedEigenSOE->getClassTag();
    } else {
        // Traditional eigen analysis - parse type and numModes
        // check type of eigenvalue analysis
        if(numArgs > 1) {
            if(firstArg=="frequency"||firstArg=="-frequency"||firstArg=="generalized"||firstArg=="-generalized") {
                generalizedAlgo = true;
            } else if(firstArg=="standard"||firstArg=="-standard") {
                generalizedAlgo = false;
            } else if(firstArg=="-findLargest") {
                findSmallest = false;
            } else if(firstArg=="genBandArpack"||firstArg=="--genBandArpack"||
                      firstArg=="genBandArpackEigen"||firstArg=="-genBandArpackEigen") {
                typeSolver = EigenSOE_TAGS_ArpackSOE;
            } else if(firstArg=="symmBandLapack"||firstArg=="-symmBandLapack"||
                      firstArg=="symmBandLapackEigen"||firstArg=="-symmBandLapackEigen") {
                typeSolver = EigenSOE_TAGS_SymBandEigenSOE;
            } else if(firstArg=="fullGenLapack"||firstArg=="-fullGenLapack"||
                      firstArg=="fullGenLapackEigen"||firstArg=="-fullGenLapackEigen") {
                typeSolver = EigenSOE_TAGS_FullGenEigenSOE;
            } else {