Function sensitivityAlgorithm

SRC/tcl/commands.cpp:1324–1417

}

int 
sensitivityAlgorithm(ClientData clientData, Tcl_Interp *interp, int argc, TCL_Char **argv)
{
	bool withRespectToRVs = true;
	bool newalgorithm = false;
	int analysisTypeTag = 1;
	if(theStaticIntegrator != 0) {
	    theSensitivityIntegrator = theStaticIntegrator;
	} else if(theTransientIntegrator != 0) {
	    theSensitivityIntegrator = theTransientIntegrator;
	}
	// 1: compute at each step (default); 2: compute by command

	if (argc < 2) {
		opserr << "ERROR: Wrong number of parameters to sensitivity algorithm." << endln;
		return TCL_ERROR;
	}
	if (theReliabilityBuilder == 0) {
		opserr << "The command 'reliability' needs to be issued before " << endln
		<< " the sensitivity algorithm can be created." << endln;
		return TCL_ERROR;
	}
	else if (theSensitivityIntegrator == 0) {
		opserr << "The sensitivity integrator needs to be instantiated before " << endln
		<< " the sensitivity algorithm can be created." << endln;
		return TCL_ERROR;
	}

	if (strcmp(argv[1],"-computeAtEachStep") == 0)
	  analysisTypeTag = 1;
	else if (strcmp(argv[1],"-computeByCommand") == 0)
	  analysisTypeTag = 2;
	else {
	  opserr << "Unknown sensitivity algorithm option: " << argv[1] << endln;
	  return TCL_ERROR;
	}
	
	ReliabilityDomain *theReliabilityDomain;
	theReliabilityDomain = theReliabilityBuilder->getReliabilityDomain();
	if(newalgorithm){
	  // theSensitivityAlgorithm = new 
	  //   NewSensitivityAlgorithm(theReliabilityDomain, 
	  // 			    &theDomain,
	  // 			    theAlgorithm,
	  // 			    theSensitivityIntegrator,
	  // 			    analysisTypeTag);
	} else {
	//  theSensitivityAlgorithm = new 
	 //   SensitivityAlgorithm(&theDomain,
	//			 theAlgorithm,
	//			 theSensitivityIntegrator,
	//			 analysisTypeTag);
	
	
    IncrementalIntegrator *theIntegrator = 0;
      
	   if (theStaticAnalysis != 0 && theStaticIntegrator != 0) {
               theIntegrator = theStaticIntegrator;