Method analyze

SRC/reliability/analysis/analysis/FOSMAnalysis.cpp:82–363

int 
FOSMAnalysis::analyze(void)
{

	// Alert the user that the FORM analysis has started
	opserr << "FOSM Analysis is running ... " << endln;

	// Open output file
	ofstream outputFile( fileName, ios::out );
    int result;

	// Get number of random variables and limit-state function
    int nrv = theReliabilityDomain->getNumberOfRandomVariables();
    int numberOfParameters = theOpenSeesDomain->getNumParameters();
    int numLsf = theReliabilityDomain->getNumberOfLimitStateFunctions();
	LimitStateFunction *theLSF;
	
	// Get mean point
	Vector meanVector(nrv);
	// Establish vector of standard deviations
	Vector stdvVector(nrv);
    

    // set start point to be the mean for FOSM
    for (int j = 0; j < nrv; j++) {
        RandomVariable *theRV = theReliabilityDomain->getRandomVariablePtrFromIndex(j);
        int param_indx = theReliabilityDomain->getParameterIndexFromRandomVariableIndex(j);
        Parameter *theParam = theOpenSeesDomain->getParameterFromIndex(param_indx);
        
        meanVector(j) = theRV->getMean();
        stdvVector(j) = theRV->getStdv();
        
        // now we should update the parameter value
        theParam->update( meanVector(j) );
    }
        
        
	// Perform analysis using mean vector
	Vector meanEstimates(numLsf);
    
    // set values in the variable namespace
    if (theGFunEvaluator->setVariables() < 0) {
        opserr << "FOSMAnalysis::analyze() - " << endln
            << " could not set variables in namespace. " << endln;
        return -1;
    }
    
    // run analysis
	if (theGFunEvaluator->runAnalysis() < 0) {
		opserr << "FOSMAnalysis::analyze() - " << endln
			<< " could not run analysis to evaluate limit-state function. " << endln;
		return -1;
	}
	
	// Loop over number of limit-state functions and evaluate
	//LimitStateFunctionIter lsfIter = theReliabilityDomain->getLimitStateFunctions();
	//while ((theLSF = lsfIter()) != 0) {
	for (int lsf = 0; lsf < numLsf; lsf++ ) {