-void ConvertToIeeeExtended(double num, unsigned char *bytes)
-{
- int sign;
- int expon;
- double fMant, fsMant;
- unsigned long hiMant, loMant;
-
- if (num < 0) {
- sign = 0x8000;
- num *= -1;
- } else {
- sign = 0;
- }
-
- if (num == 0) {
- expon = 0; hiMant = 0; loMant = 0;
- }
- else {
- fMant = frexp(num, &expon);
- if ((expon > 16384) || !(fMant < 1)) { /* Infinity or NaN */
- expon = sign|0x7FFF; hiMant = 0; loMant = 0; /* infinity */
- }
- else { /* Finite */
- expon += 16382;
- if (expon < 0) { /* denormalized */
- fMant = ldexp(fMant, expon);
- expon = 0;
- }
- expon |= sign;
- fMant = ldexp(fMant, 32);
- fsMant = floor(fMant);
- hiMant = FloatToUnsigned(fsMant);
- fMant = ldexp(fMant - fsMant, 32);
- fsMant = floor(fMant);
- loMant = FloatToUnsigned(fsMant);
- }
- }
-
- /* disable the warning about 'possible loss of data' */
- #ifdef _MSC_VER
- #pragma warning(disable: 4244)
- #endif /* Visual C++ */
-
- bytes[0] = (expon >> 8) & 0xff;
- bytes[1] = expon & 0xff;
- bytes[2] = (unsigned char) ((hiMant >> 24) & 0xff);
- bytes[3] = (unsigned char) ((hiMant >> 16) & 0xff);
- bytes[4] = (unsigned char) ((hiMant >> 8) & 0xff);
- bytes[5] = (unsigned char) (hiMant & 0xff);
- bytes[6] = (unsigned char) ((loMant >> 24) & 0xff);
- bytes[7] = (unsigned char) ((loMant >> 16) & 0xff);
- bytes[8] = (unsigned char) ((loMant >> 8) & 0xff);
- bytes[9] = (unsigned char) (loMant & 0xff);
-
- #ifdef _MSC_VER
- #pragma warning(default: 4244)
- #endif /* Visual C++ */
-}
+/****************************************************************
+ * The following two routines make up for deficiencies in many
+ * compilers to convert properly between unsigned integers and
+ * floating-point. Some compilers which have this bug are the
+ * THINK_C compiler for the Macintosh and the C compiler for the
+ * Silicon Graphics MIPS-based Iris.
+ ****************************************************************/