Function printB

SRC/runtime/commands/analysis/analysis.cpp:641–703

}

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

  int res = 0;

  FileStream outputFile;
  OPS_Stream *output = &opserr;

  bool ret = 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 << "print <filename> .. - failed to open file: "
               << argv[currentArg] << endln;
        return TCL_ERROR;
      }
      output = &outputFile;
    } else if ((strcmp(argv[currentArg], "ret") == 0) ||
               (strcmp(argv[currentArg], "-ret") == 0)) {
      ret = true;
    }
    currentArg++;
  }

  LinearSOE  *theSOE = builder->getLinearSOE();
  if (theSOE != nullptr) {

    // TODO
    builder->formUnbalance();

    if (theSOE->getNumEqn() == 0) {
      opserr << OpenSees::PromptValueError << "System of equations is empty\n";
      return TCL_ERROR;
    }

    const Vector &b = theSOE->getB();

    if (ret) {
      const int size = b.Size();
      Tcl_Obj* list = Tcl_NewListObj(size, nullptr);
      for (int i = 0; i < size; ++i)
        Tcl_ListObjAppendElement(interp, list, Tcl_NewDoubleObj(b[i]));

      Tcl_SetObjResult(interp, list);
    } else {
      *output << b;
      outputFile.close();