Method fitCurve

SRC/reliability/analysis/analysis/UnivariateDecomposition.cpp:113–194

*/

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


	

	ofstream tt("point_shape_coeff_uni.out");


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

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


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

			if (isTimeVariant){
				coefficients_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("univariate_coefficient_recorder.out");
	for( ii=0; ii<numAxes; ii++){
		opserr<<"\n\ncoefficients["<<ii<<"]"<<endln;
		file_coeff<<"coefficients["<<ii<<"]"<<endln;

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