Method computew

SRC/element/forceBeamColumn/ForceBeamColumnCBDI3d.cpp:2787–2824

}

void
ForceBeamColumnCBDI3d::computew(Vector &w, Vector &wp, double xi[],
				const Vector &kappa, const Vector &gamma)
{
  double L = crdTransf->getInitialLength();

  Matrix ls(numSections, numSections);

  Matrix Ginv(numSections, numSections);
  this->getGinv(numSections, xi, Ginv);

  Matrix H(numSections, numSections);

  bool isGamma = false;
  for (int i = 0; i < numSections; i++) {
    if (gamma(i) != 0.0)
      isGamma = true;
  }
  isGamma = CSBDI && isGamma;

  this->getHk(numSections, xi, H);
  ls.addMatrixProduct(0.0, H, Ginv, 1.0);
  w.addMatrixVector(0.0, ls, kappa, L*L);

  if (isGamma) {
    this->getHg(numSections, xi, H);
    ls.addMatrixProduct(0.0, H, Ginv, 1.0);
    w.addMatrixVector(1.0, ls, gamma, L);
    
    this->getHkp(numSections, xi, H);
    ls.addMatrixProduct(0.0, H, Ginv, 1.0);
    wp.addMatrixVector(0.0, ls, kappa, L);
    
    this->getHgp(numSections, xi, H);
    ls.addMatrixProduct(0.0, H, Ginv, 1.0);
    wp.addMatrixVector(1.0, ls, gamma, 1.0);
  }
}

void
ForceBeamColumnCBDI3d::computedwdq(Matrix &dwidq, const Vector &q,