Method fitCurve

SRC/reliability/analysis/analysis/BivariateDecomposition.cpp:166–322

*/

int BivariateDecomposition::fitCurve() {  // Lagrange interpolation . refer "Decomposition methods for structural reliability analysis" H. Xu, etc. 2005

  int ii;
	for(ii=0; ii<numAxes; ii++){ // for a special axis
		thePrincipalAxes[ii]->computeShapeFuncCoeff();
	}

//* --- deal  with axes --------
	
	ofstream tt("point_shape_coeff.out");


	

	for(   ii =0; ii<numAxes; ii++){
//		tt << "ii="<<ii<<endln;	
		if (axisCoeff[ii] !=0) delete axisCoeff[ii];

		int m = thePrincipalAxes[ii]->getExperimentalPointRule()->getNumberOfPoints();
		axisCoeff[ii] = new Vector(m);
		axisCoeff[ii]->Zero();


		if (isTimeVariant){ //---
			
			if (axisCoeff_2[ii] !=0) delete axisCoeff_2[ii];
			axisCoeff_2[ii] = new Vector(m);
			axisCoeff_2[ii]->Zero();
		
		
		} //---
		
		for ( int point=0; point < m; point++){
			Vector * theCoeff = thePrincipalAxes[ii]->getShapeFuncCoeff(point);
			axisCoeff[ii]->addVector(1.0, *theCoeff, thePrincipalAxes[ii]->getValueOnAxis(point));

			if (isTimeVariant){
				axisCoeff_2[ii]->addVector(1.0, *theCoeff, thePrincipalAxes[ii]->getValueG2OnAxis(point));			
			}
			tt<<"point num="<<point<<", shape_coeff:"<<endln;
			for(int i =0; i<theCoeff->Size(); i++)
				tt<<(*theCoeff)(i)<<endln;

		} // for each point
	
	} // for ii
	tt.close();  
	
	//* -- debug 
	ofstream file_coeff("Bivariate_coefficient_recorder.out");
	for( ii=0; ii<numAxes; ii++){
		opserr<<"\n\axisCoeff["<<ii<<"]"<<endln;
		file_coeff<<"axisCoeff["<<ii<<"]"<<endln;

		for (int m=0; m<(*axisCoeff[ii]).Size(); m++){
			opserr<<(*axisCoeff[ii])(m)<<endln;
			file_coeff<<(*axisCoeff[ii])(m)<<endln;
		}