At Wed, 21 Mar 2007 12:02:06 -0600, Patrick Alken wrote: > Is this a known issue? I'm trying to use gsl_finite() in testgen > to better sort the eigenvalues so I can compare them to lapack > if there are infinities present. Basically I'm getting segfaults > all over the place whenever I try to access nans/infs in > double-precision mode.
You can use GSL_IEEE_MODE=double-precision,mask-all to avoid that, or just call finite() directly. It does look like the order of the #ifdefs in sys/infnan.c could be improved there to pick up the right functions though. -- Brian Gough
