Function specifyAlgorithm

SRC/tcl/commands.cpp:3822–4264

command invoked to allow the SolnAlgorithm object to be built

// command invoked to allow the SolnAlgorithm object to be built
//
int
specifyAlgorithm(ClientData clientData, Tcl_Interp *interp, int argc, 
		 TCL_Char **argv)
{
  // make sure at least one other argument to contain numberer
  if (argc < 2) {
      opserr << "WARNING need to specify an Algorithm type \n";
      return TCL_ERROR;
  }    
  EquiSolnAlgo *theNewAlgo = 0;
  OPS_ResetInputNoBuilder(clientData, interp, 2, argc, argv, &theDomain);

  // check argv[1] for type of Algorithm and create the object
  if (strcmp(argv[1],"Linear") == 0) {
    int formTangent = CURRENT_TANGENT;
    int factorOnce = 0;
    int count = 2;
    while (count < argc) {
      if ((strcmp(argv[count],"-secant") == 0) || (strcmp(argv[count],"-Secant") == 0)) {
	formTangent = CURRENT_SECANT;
      } else if ((strcmp(argv[count],"-initial") == 0) || (strcmp(argv[count],"-Initial") == 0)) {
	formTangent = INITIAL_TANGENT;
      } else if ((strcmp(argv[count],"-factorOnce") == 0) || (strcmp(argv[count],"-FactorOnce") ==0 )) {
	factorOnce = 1;
      }
      count++;
    }
    theNewAlgo = new Linear(formTangent, factorOnce);
  }

  else if (strcmp(argv[1],"Newton") == 0) {
    void *theNewtonAlgo = OPS_NewtonRaphsonAlgorithm();
    if (theNewtonAlgo == 0)
      return TCL_ERROR;

    theNewAlgo = (EquiSolnAlgo *)theNewtonAlgo;
    if (theTest != 0)
      theNewAlgo->setConvergenceTest(theTest);
  }

  else if ((strcmp(argv[1],"NewtonHallM") == 0) || (strcmp(argv[1],"NewtonHall") == 0)) {
    void *theNewtonAlgo = OPS_NewtonHallM();
    if (theNewtonAlgo == 0)
      return TCL_ERROR;

    theNewAlgo = (EquiSolnAlgo *)theNewtonAlgo;
    if (theTest != 0)
      theNewAlgo->setConvergenceTest(theTest);
  }

  else if (strcmp(argv[1],"ModifiedNewton") == 0) {
    void *theNewtonAlgo = OPS_ModifiedNewton();
    if (theNewtonAlgo == 0)
      return TCL_ERROR;

    theNewAlgo = (EquiSolnAlgo *)theNewtonAlgo;
    if (theTest != 0)
      theNewAlgo->setConvergenceTest(theTest);