Function sensNodePressure

SRC/runtime/commands/domain/sensitivity.cpp:193–242

}

int
sensNodePressure(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 < 3) {
    opserr << OpenSees::PromptValueError << "want - sensNodePressure nodeTag? paramTag?\n";
    return TCL_ERROR;
  }

  int tag, paramTag;

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

  double dp = 0.0;
  Pressure_Constraint *thePC = theDomain->getPressure_Constraint(tag);
  if (thePC != nullptr) {
    // int ptag = thePC->getPressureNode();
    // Node* pNode = theDomain->getNode(ptag);
    Node *pNode = thePC->getPressureNode();
    if (pNode != nullptr) {

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

      int gradIndex = theParam->getGradIndex();
      dp = pNode->getVelSensitivity(1, gradIndex);
    }
  }

  Tcl_SetObjResult(interp, Tcl_NewDoubleObj(dp));

  return TCL_OK;
}

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