Method formTangDispSensitivity

SRC/analysis/integrator/DisplacementControl.cpp:619–744

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

//Added by Abbas
//obtain the derivative of the tangent displacement (dUhatdh)
   Vector *
DisplacementControl::formTangDispSensitivity(Vector *dUhatdh,int gradNumber)
{
  // opserr<<"DisplacementControl:; formTangDispsensitivity:::::::::::Start"<<endln;
   LinearSOE *theLinSOE = this->getLinearSOE(); 
   dUhatdh->Zero();
   dphatdh->Zero();
//   opserr<<"DisplacementControl::dUhatdh is  "<<*dUhatdh<<endln;//

   // To get the structural stiffness Matrix using the full general system of equations
   //...............................................................
 //  this->formTangent();
 //  Matrix K(size,size);
 //   K.Zero();
  // K=this->ActivateSensitivity();
  // // to print K to the screen
   //..........................
  /* 
      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;

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

   // form dKdh

 // Matrix dKdh(size,size);
 //  dKdh.Zero();
//  dKdh=this->getdKdh();
//
   // 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*deltaUhat
  // dUhatdh->addMatrixVector(0.0,dKdh,*deltaUhat,-1.0);