Levenberg-Marquardt


cliffj
06-19-2003, 10:53 AM
Hi all,

I am using MRQMIN (GNU-FORTRAN-77) to try to fit a Gaussian-Exponential convolution function (EMG) to mass spectrometry data. I have been able to modify FUNCS to fit 4 gaussians at defined centers etc, however, the EMG function and derivatives are fairly complex. Gaussj is failing because IPIV(k).gt.1. Any advice on possible causes for this would be greatly appreciated.

Thanks in advance,

John

John.cliff@oregonstate.edu

kibitzer
06-20-2003, 12:40 PM
Hi,

Have you tried using double precision?