Function OPS_System

SRC/interpreter/OpenSeesCommands.cpp:1402–1534

}

int OPS_System()
{
    if (OPS_GetNumRemainingInputArgs() < 1) {
    	opserr << "WARNING insufficient args: system type ...\n";
    	return -1;
    }

    const char* type = OPS_GetString();

    // create soe
    LinearSOE* theSOE = 0;

    if ((strcmp(type,"BandGeneral") == 0) || (strcmp(type,"BandGEN") == 0)
    	|| (strcmp(type,"BandGen") == 0)){
	// BAND GENERAL SOE & SOLVER
    	theSOE = (LinearSOE*)OPS_BandGenLinLapack();

    } else if (strcmp(type,"BandSPD") == 0) {
	// BAND SPD SOE & SOLVER
    	theSOE = (LinearSOE*)OPS_BandSPDLinLapack();

    } else if (strcmp(type,"Diagonal") == 0) {
	// Diagonal SOE & SOLVER
	theSOE = (LinearSOE*)OPS_DiagonalDirectSolver();


    } else if (strcmp(type,"MPIDiagonal") == 0) {
#ifdef _PARALLEL_INTERPRETERS
        MPIDiagonalSolver* theSolver = new MPIDiagonalSolver();
        theSOE = new MPIDiagonalSOE(*theSolver);
        setMPIDSOEFlag = true;
#else
	// Diagonal SOE & SOLVER
	theSOE = (LinearSOE*)OPS_DiagonalDirectSolver();
#endif

    } else if (strcmp(type,"SProfileSPD") == 0) {
	// PROFILE SPD SOE * SOLVER
	// now must determine the type of solver to create from rest of args
	theSOE = (LinearSOE*)OPS_SProfileSPDLinSolver();

    } else if (strcmp(type, "ProfileSPD") == 0) {

	theSOE = (LinearSOE*)OPS_ProfileSPDLinDirectSolver();

#ifdef _PARALLEL_INTERPRETERS
    } else if (strcmp(type, "ParallelProfileSPD") == 0) {
        ProfileSPDLinSolver* theSolver = new ProfileSPDLinDirectSolver();
        DistributedProfileSPDLinSOE* theParallelSOE = new DistributedProfileSPDLinSOE(*theSolver);
        theSOE = theParallelSOE;
        auto theMachineBroker = cmds->getMachineBroker();
        auto rank = theMachineBroker->getPID();
        auto numChannels = cmds->getNumChannels();
        auto theChannels = cmds->getChannels();
        theParallelSOE->setProcessID(rank);
        theParallelSOE->setChannels(numChannels, theChannels);
    
#endif