Function printA

SRC/runtime/commands/analysis/analysis.cpp:480–639

}

static int
printA(ClientData clientData, Tcl_Interp *interp, Tcl_Size argc, TCL_Char ** const argv)
{
  // printA <filename> - m <double> -c <double> -k <double>
  assert(clientData != nullptr);
  BasicAnalysisBuilder *builder = (BasicAnalysisBuilder*)clientData;

  int res = 0;

  enum class Format {
    None
  } format = Format::None;

  FileStream outputFile;
  OPS_Stream *output = &opserr;
  LinearSOE  *oldSOE = builder->getLinearSOE();


  // Cant allocate theSolver on stack because theSOE is going to 
  // delete it
  FullGenLinLapackSolver *theSolver = new FullGenLinLapackSolver();
  FullGenLinSOE theSOE(*theSolver);

  builder->set(&theSOE, false);
  // invoke domainChange which constructs a graph and passes
  // it to the SOE. Otherwise, getA() returns null
  builder->domainChanged();


  bool ret = false;
  double m = 0.0, c = 0.0, k = 0.0;
  bool do_mck = false;
  int currentArg = 1;
  while (currentArg < argc) {
    if ((strcmp(argv[currentArg], "file") == 0) ||
        (strcmp(argv[currentArg], "-file") == 0)) {
      currentArg++;
      if (currentArg == argc) {
        opserr << G3_WARN_PROMPT << "-file missing argument\n";
        return TCL_ERROR;
      }

      if (outputFile.setFile(argv[currentArg]) != 0) {
        opserr << "failed to open file: "
               << argv[currentArg] << "\n";
        return TCL_ERROR;
      }
      output = &outputFile;
    }
    else if ((strcmp(argv[currentArg], "ret") == 0) ||
               (strcmp(argv[currentArg], "-ret") == 0)) {
      ret = true;

    } else if ((strcmp(argv[currentArg], "-m") == 0)) {
      currentArg++;
      if (Tcl_GetDouble(interp, argv[currentArg], &m) != TCL_OK) {
        opserr << OpenSees::PromptValueError << "failed to read float following flag -m\n";
        return TCL_ERROR;