Function G3Parse_newLinearSOE

SRC/runtime/commands/analysis/solver.cpp:109–167

LinearSOE*
G3Parse_newLinearSOE(ClientData clientData, Tcl_Interp* interp, Tcl_Size argc,
                     G3_Char ** const argv)
{
  G3_Runtime* rt = G3_getRuntime(interp); 

  // argc is checked by calling function;

  LinearSOE *theSOE = nullptr;
  std::string sys_name{argv[1]};
  transform(sys_name.begin(), sys_name.end(), sys_name.begin(), ::tolower);
  auto ctor = soe_table.find(sys_name);

  if (ctor != soe_table.end()) {
    return ctor->second.ss(rt, argc, argv);
  }

#if defined(XARA_USE_MUMPS)
  else if (strcasecmp(argv[1], "mumps") == 0) {
    return TclDispatch_newMumpsLinearSOE(clientData, interp, argc, argv);
  }
#endif

  else if (strcasecmp(argv[1], "Umfpack")==0) {
    // TODO: if "umfpack" is in solver.hpp, this wont be reached
    return TclDispatch_newUmfpackLinearSOE(clientData, interp, argc, argv);
  }

#if defined(OPS_PETSC)
  else if (strcmp(argv[1], "petsc")==0 ||
           strcmp(argv[1], "Petsc")==0) {
    theSOE = TclDispatch_newPetscSOE(clientData, interp, argc, argv);
  }
#endif

  // Diagonal SOE & SOLVER
#ifdef _PARALLEL_INTERPRETERS
  else if (strcmp(argv[1], "ParallelProfileSPD") == 0) {
    ProfileSPDLinSolver *theSolver = new ProfileSPDLinDirectSolver();
    DistributedProfileSPDLinSOE *theParallelSOE =
        new DistributedProfileSPDLinSOE(*theSolver);
    theSOE = theParallelSOE;
    theParallelSOE->setProcessID(OPS_rank);
    theParallelSOE->setChannels(numChannels, theChannels);
  }
#endif

#ifdef _PARALLEL_INTERPRETERS
  else if (strcmp(argv[1], "MPIDiagonal") == 0) {
    setMPIDSOEFlag = true;
  }
#endif

  else {
    opserr << OpenSees::PromptValueError 
           << " system '" 
           << argv[1] << "' is unknown or not installed\n";
    return nullptr;