Bug #29954 for Algorithm-CurveFit: Math::MatrixReal->new_from_cols fails when $deriv->value is undef

When I call curve_fit on a parabola formula A*(x - x_0)^2 + y_0 where the initial guess for x-value of the vertex of the parabola (x_0) equals one of the x-data values submitted to the subroutine (x), it causes the following error: Math::MatrixReal::new_from_string(): syntax error in input stringString is [ ] [ 6.6666666666667 ] [ 13.3333333333333 ] [ 20 ] [ 26.6666666666667 ] [ 33.3333333333333 ] [ 40 ] --- at /usr/lib/perl5/site_perl/5.8.8/Algorithm/CurveFit.pm line 149 Please see attached test script which reproduces the error. The following parameter values were used: formula: A*(x-x_0)^2+y_0 params: (['A',-5,'0.05'], ['x_0',2,'0.02'], ['y_0',20,'0.2']) variable: x xdata: 0, 0.666666666666667, 1.33333333333333, 2, 2.66666666666667, 3.33333333333333, 4, 4.66666666666667, 5.33333333333333, 6 ydata: 0, 10, 16.6666666666667, 20, 16.6666666666667, 10, 0, -16.6666666666667, -36.6666666666667, -60 maximum iterations: 100 The cause of this error is an unexpected value of undef returned from $deriv->value: When I ask Math::Symbolic::Derivative to take the partial derivative of A*(x - x_0)^2 + y_0 with respect to x_0, and simplify, it returns the expression (-2A*(x - x_0)^2)/(x - x_0) When I do the same by hand, I get -2A*(x - x_0) The two are equivalent except when x == x_0. Then the Symbolic expression returns undef, while mine by hand returns 0. Algorithm::CurveFit has no handling for a returned value of undef at this point, so it proceeds and the error is generated at Math::MatrixReal::new_from_string(). Running perl v5.8.8 built for i686-linux, using Math::Symbolic version 0.603, and Algorithm::CurveFit version 1.04. I have submitted bug #55842 at CPAN ( https://rt.cpan.org/Ticket/Display.html?id=55842 ) for Math-Symbolic, please check there for updates and more information. - John Refior

#!/usr/bin/perl use strict; use warnings; use Algorithm::CurveFit; my $formula = 'A*(x-x_0)^2+y_0'; my @parameters = ( [ 'A', -5, '0.05' ], [ 'x_0', 2, '0.02' ], [ 'y_0', 20, '0.2' ] ); my @xDataPts = (0, 0.666666666666667, 1.33333333333333, 2, 2.66666666666667, 3.33333333333333, 4, 4.66666666666667, 5.33333333333333, 6); my @yDataPts = (0, 10, 16.6666666666667, 20, 16.6666666666667, 10, 0, -16.6666666666667, -36.6666666666667, -60); my $squareResidual = Algorithm::CurveFit->curve_fit( formula => $formula, params => \@parameters, # values in @parameters will be modified in-place variable => 'x', xdata => \@xDataPts, # data corresponsing to variable 'x' ydata => \@yDataPts, # data corresponsing to 'y' maximum_iterations => 100, # fit ends if an iteration introduces smaller changes than the accuracy ) or warn "Error-parabola"; # "Parabola function could not be fitted to data"; print "Done!\n"; print "Square Residual is $squareResidual.\n";