Function specifyAnalysis

SRC/tcl/commands.cpp:2452–2962

command invoked to allow the Analysis object to be built

// command invoked to allow the Analysis object to be built
//
int 
specifyAnalysis(ClientData clientData, Tcl_Interp *interp, int argc, 
		TCL_Char **argv)
{
    // make sure at least one other argument to contain type of system
    if (argc < 2) {
	opserr << "WARNING need to specify an analysis type (Static, Transient)\n";
	return TCL_ERROR;
    }    

    //
    // do nothing if request is for the same analysis type!
    //

    if ((strcmp(argv[1],"Static") == 0) && (theStaticAnalysis != 0))
      return TCL_OK;

    if (((strcmp(argv[1],"VariableTimeStepTransient") == 0) ||
	(strcmp(argv[1],"TransientWithVariableTimeStep") == 0) ||
	 (strcmp(argv[1],"VariableTransient") == 0)) && 
	(theVariableTimeStepTransientAnalysis != 0))
      return TCL_OK;

    if ((strcmp(argv[1],"Transient") == 0) && (theTransientAnalysis != 0))
      return TCL_OK;

    //
    // analysis changing .. delete the old analysis
    //

    if (theStaticAnalysis != 0) {
	delete theStaticAnalysis;
	theStaticAnalysis = 0;
	opserr << "WARNING: analysis .. StaticAnalysis already exists => wipeAnalysis not invoked, problems may arise\n";
    }

    if (theTransientAnalysis != 0) {
	delete theTransientAnalysis;
	theTransientAnalysis = 0;
	theVariableTimeStepTransientAnalysis = 0;
	opserr << "WARNING: analysis .. TransientAnalysis already exists => wipeAnalysis not invoked, problems may arise\n";
    }
    
    // check argv[1] for type of SOE and create it
    if (strcmp(argv[1],"Static") == 0) {
	// make sure all the components have been built,
	// otherwise print a warning and use some defaults
	if (theAnalysisModel == 0) 
	    theAnalysisModel = new AnalysisModel();

	if (theTest == 0) 
	  theTest = new CTestNormUnbalance(1.0e-6,25,0);       
	
	if (theAlgorithm == 0) {
	    opserr << "WARNING analysis Static - no Algorithm yet specified, \n";
	    opserr << " NewtonRaphson default will be used\n";	    

	    theAlgorithm = new NewtonRaphson(*theTest);