Function OPS_MumpsSolver

SRC/interpreter/OpenSeesCommands.cpp:3722–3775

}

void* OPS_MumpsSolver() {
    int icntl14 = 20;
    int icntl7 = 7;
    int matType = 0; // 0: unsymmetric, 1: symmetric positive definite, 2: symmetric general
    while (OPS_GetNumRemainingInputArgs() > 2) {
        const char* opt = OPS_GetString();
        int num = 1;
        if (strcmp(opt, "-ICNTL14") == 0) {
            if (OPS_GetIntInput(&num, &icntl14) < 0) {
                opserr << "WARNING: failed to get icntl14\n";
                return 0;
            }
        } else if (strcmp(opt, "-ICNTL7") == 0) {
            if (OPS_GetIntInput(&num, &icntl7) < 0) {
                opserr << "WARNING: failed to get icntl7\n";
                return 0;
            }
        } else if (strcmp(opt, "-matrixType") == 0) {
            if (OPS_GetIntInput(&num, &matType) < 0) {
                opserr << "WARNING: failed to get -matrixType. Unsymmetric matrix assumed\n";
                return 0;
            }
            if (matType < 0 || matType > 2) {
                opserr << "Mumps Warning: wrong -matrixType value (" << matType << "). Unsymmetric matrix assumed\n";
                matType = 0;
            }
        }
    }

#ifdef _PARALLEL_INTERPRETERS
#ifdef _MUMPS
    MumpsParallelSOE* soe = 0;

    MumpsParallelSolver *solver= new MumpsParallelSolver(icntl7, icntl14);
    soe = new MumpsParallelSOE(*solver, matType);

    MachineBroker* machine = cmds->getMachineBroker();
    Channel** channels = cmds->getChannels();
    int numChannels = cmds->getNumChannels();

    int rank = machine->getPID();
    soe->setProcessID(rank);
    soe->setChannels(numChannels, channels);
    return soe;
#endif
#else
#ifdef _MUMPS
    MumpsSolver *theSolver = new MumpsSolver(icntl7, icntl14);
    theSOE = new MumpsSOE(*theSolver, matType);
    return theSOE;
#endif
#endif
    return 0;
}

// Sensitivity:BEGIN /////////////////////////////////////////////