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.
the double precision word is stored in the big-endian word order and
the signbit is assumed to be bit 32. For
double precision word is stored in the little-endian word order and
the signbit is assumed to be bit 64.