Method debug

SRC/reliability/analysis/analysis/BivariateDecomposition.cpp:563–644

double BivariateDecomposition::debug(Vector * point, Vector * dp2prime, Vector * gradG2){


/////////////////////////////// ------------------- for debug only by Quan ----

	static int iii =1; 
	if (iii==1){
		iii++;
		ofstream file_coeff_2("Bivariate_surface_fitting_debug.out");

     	int numOfGridPlanes = numAxes*(numAxes-1)/2;

		for (int ii=0; ii<numOfGridPlanes; ii++) {

			file_coeff_2 <<"\n the "<<ii<<"  th plane: "<<endln;
			PrincipalAxis * axis_1 = theGridPlanes[ii]->getAxisPtr(1); // first axis
			PrincipalAxis * axis_2 = theGridPlanes[ii]->getAxisPtr(2);  // second axis

			 int m = axis_1->getExperimentalPointRule()->getNumberOfPoints();
			 int n = axis_2->getExperimentalPointRule()->getNumberOfPoints();




			Matrix * valueOnGrid = theGridPlanes[ii]->getGridValuesPtr();
			Matrix * valueG2OnGrid=0;
			
			if(isTimeVariant){
				valueG2OnGrid = theGridPlanes[ii]->getGridValuesG2Ptr();
			}

			double x_begin = axis_1->getExperimentalPointRule()->getBeginOfGrid();
			double x_end =   axis_1->getExperimentalPointRule()->getEndOfGrid();
			double x_interval = (x_end-x_begin)/(m-1);

			double y_begin=  axis_2->getExperimentalPointRule()->getBeginOfGrid();
			double y_end =   axis_2->getExperimentalPointRule()->getEndOfGrid();
			double y_interval = (y_end-y_begin)/(n-1);

			Vector * vector1 = axis_1->getAxisDirection();
			Vector * vector2 = axis_2->getAxisDirection();

			Vector tmpp(50); tmpp.Zero();
			file_coeff_2.precision(16);
			for( int pointX = 0; pointX < m; pointX++){
   

				for( int pointY = 0; pointY < n; pointY++){
					
					double x = x_begin + pointX*x_interval;
					double y = y_begin + pointY*y_interval;
					

					tmpp = x*(*vector1)/vector1->Norm()+y*(*vector2)/vector2->Norm();



					file_coeff_2<<"x:"<<x<<"y:"<<y<<endln;