Method commitSensitivity

SRC/element/zeroLength/ZeroLengthSection.cpp:968–989

}
    
int
ZeroLengthSection::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);
  }

  // Set some references to make the syntax nicer
  Vector &dedh = *v;
  const Matrix &tran = *A;

  dedh.Zero();

  // Compute element basic deformations ... v = A*(u2-u1)
  for (int i = 0; i < order; i++)
    for (int j = 0; j < numDOF/2; j++)
      dedh(i) += -diff(j)*tran(i,j);

  return theSection->commitSensitivity(dedh, gradIndex, numGrads);
}