Method analyze

SRC/reliability/analysis/analysis/ImportanceSamplingAnalysis.cpp:95–613

int 
ImportanceSamplingAnalysis::analyze(void)
{

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

	// Declaration of some of the data used in the algorithm
	double gFunctionValue;
	int result, I, seed = 1;
    long int k = 1;
	double det_covariance, phi, h, q;
	int numRV = theReliabilityDomain->getNumberOfRandomVariables();
    int numParam = theOpenSeesDomain->getNumParameters();
	int numLsf = theReliabilityDomain->getNumberOfLimitStateFunctions();
    
	Vector x(numRV);
	Vector z(numRV);
	Vector u(numRV);
	Vector randomArray(numRV);
	static NormalRV aStdNormRV(1,0.0,1.0);
	bool failureHasOccured = false;

	det_covariance = pow(samplingStdv, numRV);


	// Pre-compute some factors to minimize computations inside simulation loop
	static const double twopi = 2.0*acos(-1.0);
	double factor1 = 1.0 / ( pow(twopi,0.5*numRV));
	//double factor2 = 1.0 / ( pow(twopi,0.5*numRV) * sqrt(det_covariance) );
	double factor2 = factor1 / sqrt(det_covariance);


	Vector sum_q(numLsf);
	Vector sum_q_squared(numLsf);
	double var_qbar;
	double pfIn;
	double CovIn;
	char restartFileName[256];
	sprintf(restartFileName,"%s_%s","restart",fileName);
    

	if (analysisTypeTag == 1) {
		// Possible read data from file if this is a restart simulation
		if (printFlag == 2) {
			ifstream inputFile( restartFileName, ios::in );
			inputFile >> k;
			inputFile >> seed;
			if (k == 1 && seed == 1) {
			}
			else { 
				for (int i=1; i<=numLsf; i++) {
					inputFile >> pfIn;
					if (pfIn > 0.0) {
						failureHasOccured = true;
					}
					inputFile >> CovIn;
					sum_q(i-1) = pfIn*k;