Method initialize

SRC/reliability/analysis/analysis/SystemAnalysis.cpp:103–219

int
SystemAnalysis::initialize()
{
	// Initial declarations
	double beta;
    NormalRV uRV(1, 0.0, 1.0);
	
	if (strlen(betaFile) > 1 && strlen(rhoFile) > 1) {
		// read data for beta and rho from file, do not take from reliability domain
		ifstream inputBeta( betaFile, ios::in );
		if (!inputBeta) {
			opserr << "SystemAnalysis::initialize ERROR - could not open beta file " << betaFile << endln;
			return -1;
		}
		ifstream inputRho( rhoFile, ios::in );
		if (!inputRho) {
			opserr << "SystemAnalysis::initialize ERROR - could not open rho file " << rhoFile << endln;
			return -1;
		}
		
		numLsf = 0;
		while (inputBeta >> beta)
			numLsf++;
		inputBeta.clear();
		inputBeta.seekg(0);
		
		// allocate arrays
		allBetas = new Vector(numLsf);
		allPf1s = new Vector(numLsf);
		rhos = new Matrix(numLsf,numLsf);
		
		// read data from file into arrays
		for (int i=0; i < numLsf; i++ ) {
			inputBeta >> (*allBetas)(i);
			(*allPf1s)(i) = 1.0 - uRV.getCDFvalue( (*allBetas)(i) );
			
			for (int j=0; j < numLsf; j++ ) {
				if (!inputRho.eof() )
					inputRho >> (*rhos)(i,j);
				else {
					opserr << "SystemAnalysis::initialize ERROR - rho file does not contain the same data size as the beta file"
						   << endln;
					return -1;
				}
			}
		}
		
		inputBeta.close();
		inputRho.close();
		
	} else {
		// proceed with reliability domain data, requires previous FORM analysis objects

		// Number of limit-state functions
		numLsf = theReliabilityDomain->getNumberOfLimitStateFunctions();

		// Number of random variables
		int numRV = theReliabilityDomain->getNumberOfRandomVariables();