Function sensNodeVel

SRC/runtime/commands/domain/sensitivity.cpp:87–137

}

int
sensNodeVel(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 - sensNodeVel nodeTag? dof? paramTag?\n";
    return TCL_ERROR;
  }

  int tag, dof, paramTag;

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

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

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

  int gradIndex = theParam->getGradIndex();

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

  Tcl_SetObjResult(interp, Tcl_NewDoubleObj(value));

  return TCL_OK;
}

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