Method computeSensitivities

SRC/analysis/integrator/DisplacementControl.cpp:1049–1133

   int 
DisplacementControl::computeSensitivities(void)
{
  LinearSOE *theSOE = this->getLinearSOE();
  
  // Zero out the old right-hand side of the SOE
  theSOE->zeroB();

  // Form the part of the RHS which are independent of parameter
  this->formIndependentSensitivityRHS();

  AnalysisModel *theModel = this->getAnalysisModel();   
  Domain *theDomain=theModel->getDomainPtr();//Abbas
  ParameterIter &paramIter = theDomain->getParameters();
  Parameter *theParam;

  // De-activate all parameters
  while ((theParam = paramIter()) != 0)
    theParam->activate(false);
  
  // Now, compute sensitivity wrt each parameter
  int  numGrads = theDomain->getNumParameters();
  paramIter = theDomain->getParameters();

 
  while ((theParam = paramIter()) != 0) {
    // Activate this parameter
    theParam->activate(true);
    
    // Zero the RHS vector
    theSOE->zeroB();
    
    // Get the grad index for this parameter
    int gradIndex = theParam->getGradIndex();
                  
    // this->formResidualDispSensitivity( dUIJdh, gradIndex);//./.
 
    // Form the RHS
    
    //  this->formTangDispSensitivity(dUhatdh,gradIndex);
    this->formSensitivityRHS(gradIndex);
     
    this->formTangent(tangFlag);
    theSOE->solve();
    *dUIJdh=theSOE->getX();// sensitivity of the residual displacement
 
    //  opserr<<"computeSensitivities: dUfRdh is  "<<*dUIJdh<<endln;//
    //    opserr<<"deltaUbar is "<<*deltaUbar<<endln;   
    this->formTangDispSensitivity(dUhatdh,gradIndex);
    double dlamdh = this->getLambdaSensitivity(gradIndex);

    // To obtain the response sensitivity 
    theSOE->setB(*Residual);
    theSOE->solve();
    //theSOE->getX();
       
    //Vector *x=new Vector(size);
    //(*x)=theSOE->getX();