This is the mail archive of the
gsl-discuss@sourceware.org
mailing list for the GSL project.
Re: generalized eigensystems
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