Method commitSensitivity

SRC/element/zeroLength/ZeroLengthVG_HG.cpp:1419–1437

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

  double depsdh;
  int ret = 0;
  for (int mat=0; mat<numMaterials1d; mat++) {
    // compute strain and rate; set as current trial for material
    depsdh = this->computeCurrentStrain1d(mat,diff);
    ret += theMaterial1d[mat]->commitSensitivity(depsdh, gradIndex, numGrads);
  }

  return ret;  
}