Я наткнулся на эту тему в поисках того же кода. Вот что я написал в итоге:
public static byte [] Double2Real48(double d)
{
byte [] r48 = new byte[6];
byte [] da = BitConverter.GetBytes(d);
for (int i = 0; i < r48.Length; i++)
r48[i] = 0;
//Copy the negative flag
r48[5] |= (byte)(da[7] & 0x80);
//Get the expoent
byte b1 = (byte)(da[7] & 0x7f);
ushort n = (ushort)(b1 << 4);
byte b2 = (byte)(da[6] & 0xf0);
b2 >>= 4;
n |= b2;
if (n == 0)
return r48;
byte ex = (byte)(n - 1023);
r48[0] = (byte)(ex + 129);
//Copy the Mantissa
r48[5] |= (byte)((da[6] & 0x0f) << 3);//Get the last four bits
r48[5] |= (byte)((da[5] & 0xe0) >> 5);//Get the first three bits
r48[4] = (byte)((da[5] & 0x1f) << 3);//Get the last 5 bits
r48[4] |= (byte)((da[4] & 0xe0) >> 5);//Get the first three bits
r48[3] = (byte)((da[4] & 0x1f) << 3);//Get the last 5 bits
r48[3] |= (byte)((da[3] & 0xe0) >> 5);//Get the first three bits
r48[2] = (byte)((da[3] & 0x1f) << 3);//Get the last 5 bits
r48[2] |= (byte)((da[2] & 0xe0) >> 5);//Get the first three bits
r48[1] = (byte)((da[2] & 0x1f) << 3);//Get the last 5 bits
r48[1] |= (byte)((da[1] & 0xe0) >> 5);//Get the first three bits
return r48;
}
Real48 похож на IEEE 754 в том, что Mantissa будет такой же. Сдвиг битов необходим для того, чтобы Мантисса оказалась в нужном месте.
Показатель Real48 имеет смещение 129, а двойное - 1023.
Отрицательный флаг сохраняется в первом бите последнего байта.
Примечания:
Я не думаю, что этот код будет работать на машине с прямым порядком байтов. Он не проверяет NAN или INF.
Вот код, который конвертирует real48 в double. Он был портирован из компилятора Free Pascal:
static double real2double(byte [] r)
{
byte [] res = new byte[8];
int exponent;
//Return zero if the exponent is zero
if (r[0] == 0)
return (double)0;
//Copy Mantissa
res[0] = 0;
res[1] = (byte)(r[1] << 5);
res[2] = (byte)((r[1] >> 3) | (r[2] << 5));
res[3] = (byte)((r[2] >> 3) | (r[3] << 5));
res[4] = (byte)((r[3] >> 3) | (r[4] << 5));
res[5] = (byte)((r[4] >> 3) | ((r[5] & 0x7f) << 5));
res[6] = (byte)((r[5] & 0x7f) >> 3);
//Copy exponent
//correct exponent
exponent = (r[0] + (1023-129));
res[6] = (byte)(res[6] | ((exponent & 0xf) << 4));
res[7] = (byte)(exponent >> 4);
//Set Sign
res[7] = (byte)(res[7] | (r[5] & 0x80));
return BitConverter.ToDouble(res, 0);
}