This is the mail archive of the gsl-discuss@sourceware.org mailing list for the GSL project.

Index Nav: [Date Index] [Subject Index] [Author Index] [Thread Index]
Message Nav: [Date Prev] [Date Next] [Thread Prev] [Thread Next]
Other format: [Raw text]

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
Index Nav: [Date Index] [Subject Index] [Author Index] [Thread Index]
Message Nav: [Date Prev] [Date Next] [Thread Prev] [Thread Next]