I noticed that the two smooth_linear_interpolation functions defined in the interpolation.cxx file are slightly different. In the second declaration of the smooth_linear_interpolation: smooth_linear_interpolation(adouble* y, adouble& x, adouble* Xdata, adouble* Ydata, int n), the following lines are commented:
// Si = smooth_heaviside( (x-(Xdata[0]-100*a)), a ) - smooth_heaviside( (x-Xdata[0]), a );
// *y+=Si*Ydata[0];
// Si = smooth_heaviside( (x-Xdata[n-1]), a ) - smooth_heaviside( (x-(Xdata[n-1]+100*a)), a );
// *y+=Si*Ydata[n-1];
whereas these lines are not commented in the first declaration of the smooth_linear_interpolation (smooth_linear_interpolation(adouble* y, adouble& x, DMatrix& Xdata, DMatrix& Ydata, int n))
I noticed that the two
smooth_linear_interpolationfunctions defined in theinterpolation.cxxfile are slightly different. In the second declaration of thesmooth_linear_interpolation:smooth_linear_interpolation(adouble* y, adouble& x, adouble* Xdata, adouble* Ydata, int n), the following lines are commented:whereas these lines are not commented in the first declaration of the
smooth_linear_interpolation(smooth_linear_interpolation(adouble* y, adouble& x, DMatrix& Xdata, DMatrix& Ydata, int n))Which one is the intended one?