В C89, как я могу обрезать и разделить плавучую точку с двойной точностью на два 32-битных слова?

Question

При ориентации C99, я мог бы усечь и разделить значение double на два 32-битовых целых чисел со следующим кодом:

#include <stdint.h> 

void split(double d, unsigned long *most_significant_word, unsigned long *least_significant_word) 
{ 
    uint64_t i = (uint64_t)d; 
    *most_significant_word = i >> 32; 
    *least_significant_word = i & 0xffffffff; 
} 

C89 однако, не кажется, чтобы определить 64-битное целое типа, так что я могу Не используйте компилятор для выполнения усечения. Даже если усечение не было необходимым (значение уже представляло целое число), я также не мог использовать операторы битов, такие как & или >>, поскольку они не работают с значениями double.

Итак, как можно реализовать вышеприведенную функцию split() в чистом C89 (и, таким образом, не полагаясь на 64-битные целые числа), возвращая 21/32-битные слова, которые составляют 53-битовое целое, сохраненное в значении double?

Answer 1

Игнорирование знака (как это делает исходный код) он должен быть простым, как это:

void split(double d, unsigned long *most_significant_word, unsigned long *least_significant_word) 
{ 
    *most_significant_word = d/4294967296.; // d>>32 in double 
    *least_significant_word = fmod(d, 4294967296.); 
} 

Для того, чтобы отрицательные числа во внимание, действуют на абсолютное значение, то вычислить двоичное дополнение:

void split(double d, unsigned long *most_significant_word, unsigned long *least_significant_word) 
{ 
    double dabs = d < 0 ? -d : d; 
    *most_significant_word = dabs/4294967296.; 
    *least_significant_word = fmod(dabs, 4294967296.); 
    if (d < 0) { 
     *most_significant_word = ~*most_significant_word; 
     *least_significant_word = ~*least_significant_word + 1; 
     if (!*least_significant_word) *most_significant_word += 1; 
    } 
}