Method commitSensitivity

SRC/element/zeroLength/CoupledZeroLength.cpp:958–978

}
 
int
CoupledZeroLength::commitSensitivity(int gradIndex, int numGrads)
{
  // Get nodal displacement sensitivity
  Vector diff(numDOF/2);
  for (int i = 0; i < numDOF/2; i++) {
    diff(i) = theNodes[1]->getDispSensitivity(i+1,gradIndex) - theNodes[0]->getDispSensitivity(i+1,gradIndex);
  }

  dX = diff(dirn1);
  dY = diff(dirn2);
  double depsdh = sqrt(dX*dX + dY*dY);

  // strain neg if to left of X+Y = 0 line
  if (dX < 0.0 || dY < 0.0) {
    if (dX + dY < 0.0)
      depsdh *= -1.0;
  }

  return theMaterial->commitSensitivity(depsdh, gradIndex, numGrads);
}