Method analyze

SRC/reliability/analysis/analysis/MultiDimVisPrincPlane.cpp:144–652

}
int MultiDimVisPrincPlane::analyze()
{

// --- type =0:   Limit state surface -----
// --- type =1:   Limit function -------
	if (type ==0){
		opserr<<"not implemented yet."<<endln;
		exit(-1);
	
	}
    //	double tolFun =1.0e-6;

// --1. create principalplanes, set grid info ----

	
	if (thePrincipalPlanes[0]==0) {
		for(int i =0; i<numOfPrinPlane; i++){

			thePrincipalPlanes[i]=new PrincipalPlane(i+1, theDesignPoint, 
								   0, 0, theProbabilityTransformation,
								   theGFunEvaluator, 0.0);
		//	thePrincipalPlanes[i]->setDesignPoint( theDesignPoint);

		}  /// 
	}


	ofstream debug("debug.out"  );

     
	//2.  --- compute and form Hessian by FDM; compute A, Rotation R, A=R*H*R'/|dG|, but one dim less.

	int numRV = theReliabilityDomain->getNumberOfRandomVariables();
  	
	//if (theHessian ==0){
		//theHessian = new Hessian(numRV,	theReliabilityDomain,theProbabilityTransformation,
        //                         theGFunEvaluator,theGradGEvaluator, 1.0e-5);
		//theHessian->formReducedHessian(theDesignPtXSpace);
	//}

	if (HessianFileName !=0){  //recorder Hessian

		Matrix tmpHessian(numRV, numRV);
		//tmpHessian.addMatrix(0.0, theHessian->getHessianApproximation(), 1.0);

		ofstream output(HessianFileName);
		
		int ii=0;
		int jj=0;
	
		for(ii=0; ii<numRV; ii++){ 
			for (jj=0; jj<numRV;jj++){
				output << tmpHessian(jj,ii)<<endln ;
			}	
		
		}
		output.close();