Преобразование 80-битной расширенной точности в java

Я не уверен, можно ли это считать ответом или вообще можно ответить, потому что... Вы уверены, что это стандартный 80-битный формат с плавающей запятой? Потому что в таком виде это выглядит немного странно. Возьмите свой первый пример:

40 00 9e...

Здесь 0x4000−16383=1 — показатель степени. Мантисса начинается со старшего бита 1, что соответствует целой части, поэтому она должна иметь форму 1,... x 2^1, что больше 2. Но вы говорите, что она должна быть 1235. Не имеет смысла.

Я создал этот уродливый хак на Java, который как бы собирает двойника:

long high = 0x40_00L;
long low = 0x9e_06_52_14_1e_f0_db_f6L;
long e = (((high & 0x7FFFL) - 16383) + 1023) & 0x7FFL;
long ld = ((high & 0x8000L) << 48)
          | (e << 52)
          | ((low >>> 11) & 0xF_FFFF_FFFF_FFFFL);
System.out.printf("%16X\n", ld);
double d = Double.longBitsToDouble(ld);
System.out.println(d);

Предполагая, что ввод является правильным 80-битным значением, этот фрагмент должен работать нормально, за исключением того, что он не округляется должным образом, не проверяет переполнение, не обрабатывает ненормализованные значения и, вероятно, не может обрабатывать специальные значения, такие как NaN. Для этого ввода он печатает

4003C0CA4283DE1B
2.46913578

Не 1235! Затем я использовал этот код C++ (скомпилированный с помощью GCC, где long double в 80-битном формате):

#include <stdio.h>

int main() {
    long long high = 0x4000;
    long long low = 0x9e0652141ef0dbf6;
    long double d;
    char *pd = (char*)&d, *ph = (char*)&high, *pl = (char*)&low;
    for (int i = 0; i < 8; ++i) {
        pd[i] = pl[i];
    }
    for (int i = 0; i < 2; ++i) {
        pd[8 + i] = ph[i];
    }
    printf("%lf\n", (double) d);
}

Он также печатает 2,469136.