Function sensNodeDisp

SRC/runtime/commands/domain/sensitivity.cpp:34–85

int
sensNodeDisp(ClientData clientData, Tcl_Interp *interp, int argc,
             TCL_Char ** const argv)
{
  assert(clientData != nullptr);
  Domain *theDomain = (Domain *)clientData;

  // make sure at least one other argument to contain type of system
  if (argc < 4) {
    opserr << OpenSees::PromptValueError << "want - sensNodeDisp nodeTag? dof? paramTag?\n";
    return TCL_ERROR;
  }

  int tag, dof, paramTag;

  if (Tcl_GetInt(interp, argv[1], &tag) != TCL_OK) {
    opserr
        << OpenSees::PromptValueError << "nodeDisp nodeTag? dof? paramTag?- could not read nodeTag? ";
    return TCL_ERROR;
  }
  if (Tcl_GetInt(interp, argv[2], &dof) != TCL_OK) {
    opserr << OpenSees::PromptValueError << "nodeDisp nodeTag? dof? paramTag?- could not read dof? ";
    return TCL_ERROR;
  }
  if (Tcl_GetInt(interp, argv[3], &paramTag) != TCL_OK) {
    opserr << OpenSees::PromptValueError << "nodeDisp paramTag? dof? paramTag?- could not read "
              "paramTag? ";
    return TCL_ERROR;
  }

  //
  //
  Node *theNode = theDomain->getNode(tag);
  if (theNode == nullptr) {
    opserr << "sensNodeDisp: node " << tag << " not found" << "\n";
    return TCL_ERROR;
  }

  Parameter *theParam = theDomain->getParameter(paramTag);
  if (theParam == nullptr) {
    opserr << "sensNodeDisp: parameter " << paramTag << " not found" << "\n";
    return TCL_ERROR;
  }

  int gradIndex = theParam->getGradIndex();

  double value = theNode->getDispSensitivity(dof, gradIndex);

  Tcl_SetObjResult(interp, Tcl_NewDoubleObj(value));

  return TCL_OK;
}

int
sensNodeVel(ClientData clientData, Tcl_Interp *interp, int argc,