Method formSensitivityRHS

SRC/analysis/integrator/DisplacementControl.cpp:861–990

   
int
DisplacementControl::formSensitivityRHS(int passedGradNumber)
{
  sensitivityFlag = 1;
  //this->Activatesensitivity();//Abbas
  // Set a couple of data members
  gradNumber = passedGradNumber;

  // get model
  AnalysisModel* theAnalysisModel = this->getAnalysisModel();
  LinearSOE* theSOE = this->getLinearSOE();

  // Loop through elements
  FE_Element *elePtr;
  FE_EleIter &theEles = theAnalysisModel->getFEs(); 

  while((elePtr = theEles()) != 0) {
    theSOE->addB(elePtr->getResidual(this) ,elePtr->getID()  );
  }

  //  opserr<<" dKdh*deltaUbar after while loop is "<<*dUIJdh<<endln;
  (*Residual)=theSOE->getB();
 
  // opserr<<"RHS():getB()............"<<theSOE->getB()<<endln;
  // (*Residual) += (*dUIJdh);

  //   (*Residual2)=theSOE->getB();

  //  opserr<<" the new RHS is  "<<*Residual<<endln;
  //  opserr<<"where dUIJ is "<<*dUIJdh<<endln;

  
  //    if(CallParam==1) {
  //  opserr<<"inside if statement of gradIndex # "<<gradNumber<<endln;
  //  this->formTangDispSensitivity(dUhatdh,gradNumber);
  //  this->formdLambdaDh(gradNumber);
  // }
  

  //   opserr<<"RHS: dLambdaStepdh is "<<(*dLAMBDAdh)(gradNumber)<<endln;

  double CallDlambda1dh=(*dLAMBDAdh)(gradNumber);
  // opserr<<"RHS:: dphatdh is ....................................//"<<*dphatdh<<endln;
  Residual->addVector(1.0,*phat, CallDlambda1dh ); //needed to calculate dLambdadh
  //opserr<<"The residual sensitivity... dRdh= "<<*Residual<<endln;
  //    opserr<<"Pref=  "<<*phat<<endln;
  //    opserr<<"dLambdadh= "<<CallDlambda1dh<<endln;
  //opserr<<" dPrdh after adding dLamdh= "<<*Residual<<endln;
  Residual->addVector(1.0,*dphatdh, currentLambda ); //needed to calculate dLambdadh
  
  //  Matrix dKdh(size,size);
  //	 dKdh.Zero();
  //	 dKdh=this->getdKdh();
  //	 Residual->addMatrixVector(1.0,dKdh,*deltaUbar,-1.0);
  
  
  /*
    Matrix Kt(size,size);