Method analyze

SRC/reliability/analysis/analysis/DP_RSM_Sim_TimeVariant.cpp:201–1194

}

int DP_RSM_Sim_TimeVariant::analyze()
{
  //1  form Hessian;
  //2  get  eigenvectors & eigenvalues
  //3  order eigenVectors
  //4. grid values computation
  //5. fit curve
  //6. simulation


	// --1.1. create grid planes, set grid info ----

	int numOfGridPlanes = (numOfPrincipalAxes+1)*numOfPrincipalAxes/2;
 
    			
	ofstream output( outputFileName, ios::out );



	//1.  --- 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-3);
		//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();
		
	}

	//theRespSurfaceSimulation->setNormOfGradientU(theHessian->getNormOfGradientU());



    if (HessianMatrix !=0) delete HessianMatrix;
	if (rotation !=0) delete rotation;
	 
	//HessianMatrix = new Matrix(theHessian->getReducedHessian()); //A matrix, one dim less the nrv
    HessianMatrix = new Matrix(1,1);