Function TclCommand_newNewtonHallM

SRC/runtime/commands/analysis/algorithm.cpp:278–336

int
TclCommand_newNewtonHallM(ClientData clientData, Tcl_Interp* interp, Tcl_Size argc, TCL_Char**const argv)
{
  BasicAnalysisBuilder *builder = (BasicAnalysisBuilder *)clientData;
  assert(builder != nullptr);

  int method = 0;
  double iFactor = .1;
  double alpha = .01;
  double c = 100;
  double data[2];
  
  int numData = 1;
  if (OPS_GetDoubleInput(&numData,&data[0]) < 0) {
    opserr << "WARNING invalid data reading 2 hall factors\n";
    return 0;
  }
  iFactor = data[0];

  for (int i=2; i<argc; i++) {

      if (strcmp(argv[i], "-exp")==0 || 
          strcmp(argv[i], "-Exp")==0) {
        numData = 1;
        if (OPS_GetDoubleInput(&numData,&data[0]) < 0) {
          opserr << "WARNING invalid data reading 2 hall factors\n";
          return 0;
        } else 
          alpha = data[0];

      } else if (strcmp(argv[i], "-sigmoid")==0 || 
                 strcmp(argv[i], "-Sigmoid")==0) {
        method = 1;
        int numData = 2;
        if (OPS_GetDoubleInput(&numData,&data[0]) < 0) {
          opserr << "WARNING invalid data reading 2 hall factors\n";
          return 0;
        } else {
          alpha = data[0];
          c = data[1];
        }

      } else if (strcmp(argv[i], "-constant")==0 || 
                 strcmp(argv[i], "-Constant")==0) {
        method = 2;
        int numData = 1;
        if (OPS_GetDoubleInput(&numData,&data[0]) < 0) {
          opserr << "WARNING invalid data reading 2 hall factors\n";
          return 0;
        } else {
          c = data[0];
        }
      }
    }

  auto algorithm = new NewtonHallM(iFactor, method, alpha, c);
  builder->set(algorithm);
  return TCL_OK;