Method analyze

SRC/reliability/analysis/analysis/OutCrossingAnalysis.cpp:93–510

}

int 
OutCrossingAnalysis::analyze(void)
{

	// Alert the user that the analysis has started
	opserr << "Out-Crossing Analysis is running ... " << endln;

	// Declare variables used in this method
	int numRV = theReliabilityDomain->getNumberOfRandomVariables();
	static NormalRV aStdNormRV(1,0.0,1.0);
	int numVel, i, j, k, kk, nodeNumber, dofNumber; 
	double dgduValue, accuSum;
	Vector uStar2(numRV);
	Vector alpha(numRV), alpha2(numRV), alpha_k(numRV), alpha_kk(numRV);
	double beta1, beta2;
	double pf1, pf2;
//	int n_2;
//	double a, b, integral, h, fa, fb, sum_fx2j, sum_fx2j_1, Pmn1;
//	double a, b, c, integral, fa, fb, fc;
	double a, Pmn1;
	char string[500];


	// Determine number of points
	double nsteps = stepsToEnd-stepsToStart;       //total integration steps
	int numPoints = (int)floor(nsteps/sampleFreq);  //  how many time points of computation of outcrossingRate
	numPoints++;   
	double dt = theGFunEvaluator->getDt();  // integration time step
	double Dt = dt*impulseFreq;              // time between impuses.
	double T = dt*sampleFreq*numPoints;     // total time of last analysis


	// Allocate reult vectors
	Vector nu(numPoints);
	Vector ED(numPoints);
	Vector pf(numPoints);
	Vector beta(numPoints);
	Matrix Pmn2(numPoints,numPoints);
	Matrix allAlphas(numRV,numPoints);


	// Open output file and start writing to it
	ofstream outputFile( fileName, ios::out );


	// Loop over number of limit-state functions and perform analysis
	LimitStateFunctionIter lsfIter = theReliabilityDomain->getLimitStateFunctions();
	LimitStateFunction *theLimitStateFunction;
	//for (lsf=1; lsf<=numLsf; lsf++ ) {
	while ((theLimitStateFunction = lsfIter()) != 0) {

	  int lsf = theLimitStateFunction->getTag();

		// Inform the user which limit-state function is being evaluated
		opserr << "Limit-state function number: " << lsf << endln;


		// Start printing results to the output file