jlapeyre / PDL-Fit-Levmar

PDL interface to constrained Levenberg-Marquardt library levmar
Other
3 stars 3 forks source link

Fixups, upgrade to levmar-2.6, make work with PDL >2.057 #6

Closed mohawk2 closed 2 years ago

mohawk2 commented 2 years ago

The last commit fixes a long-standing problem of t/chkjac.t failing for float. The fix was to just cast the return value of exp in the lpp function to (FLOAT). Otherwise a floating-point difference happens, which only for the first value of the function gets amplified into a difference of about 1e-1, far outside of the epsilon which is 1e-4.