computes the number of negative eigenvalues of (
). This implementation of the Sturm Sequence loop exploits IEEE arithmetic and has no conditionals in the innermost loop. The signbit for real
is assumed to be bit 32. Double-precision
differ in the order of the double precision word storage and, consequently, in the signbit location. For
, the double precision word is stored in the big-endian word order and the signbit is assumed to be bit 32. For
, the double precision word is stored in the little-endian word order and the signbit is assumed to be bit 64.