Method formTangDispSensitivity

SRC/analysis/integrator/ArcLength.cpp:569–698

/////////////////////////////// Sensitivity Begin///////////////////////////// Added by Abbas obtain the derivative of the tangent displacement (dUhatdh)

//Added by Abbas
//obtain the derivative of the tangent displacement (dUhatdh)
   void
ArcLength::formTangDispSensitivity(int gradNumber)
{
   AnalysisModel *theModel=this->getAnalysisModel();
   int size=theModel->getNumEqn();
   LinearSOE *theLinSOE = this->getLinearSOE(); 
// To get the structural stiffness Matrix
//...............................................................
dphatdh->Zero();
//static Matrix K(size,size);
//K.Zero();
//this->formTangent();
//opserr<<" before getK"<<endln;
//K=this->getK();
//opserr<<"after getK"<<endln;

/*
opserr<<"the tangent printed from the DisplacementControl.cpp"<<endln;
for(int i=0;i<size;i++)
{
for(int j=0;j<size;j++)
{
opserr<<K(i,j)<<"      ";

}
opserr<<endln;

}
*/
// ................................................................


//this->formTangentSensitivity(CURRENT_TANGENT);
//static Matrix dKdh(size,size);
//dKdh.Zero();
//dKdh=this->getdKdh(gradNumber);
// To print dKdh to the Screen 
//................................................................

//opserr<<"dKdh from the DisplacementControl.cpp"<<endln;
//for(int i=0;i<size;i++)
//{
//for(int j=0;j<size;j++)
//{
//opserr<<dKdh(i,j)<<"      ";

//}
//opserr<<endln;
//}

//form dKdh*Uft
//...............................................................

//dphatdh->addMatrixVector(1.0,dKdh,*deltaUhat,-1.0);
//call the tangent (K)
this->formTangent();
theLinSOE->setB(*dphatdh);
if(theLinSOE->solve()<0) {