ralna / RALFit

A non-linear least squares solver
Other
25 stars 6 forks source link

Lanczos example (LAPACK) #23

Open talassio opened 5 years ago

talassio commented 5 years ago

Example program Lanczos.f90 causes exception Warning: Floating divide by zero occurred in call to

call dgesvd( JOBU, JOBVT, n, m, w%Jcopy, n, w%S, w%S, 1, w%S, 1, &

See https://github.com/talassio/RALFit/blob/38b3dc820570ae7b66294229dd7182a3782a513d/libRALFit/src/ral_nlls_internal.f90#L3456

Note: Compiler flag -ieee=stop will break example.