Глюк в BSD-шной библиотеке

[vak]

Обнаружил ошибку в математической библиотеке под Chimera Linux на процессоре arm64.

• log10l(1e444L) выдаёт inf, а должно быть 444
• log10l(1e555L) выдаёт inf, а должно быть 555
• log10l(1e-666L) выдаёт -inf, а должно быть -666
• log10l(1e-777L) выдаёт -inf, а должно быть -777
• log10l(-1e-999L) выдаёт -inf, а должно быть nan

На интеловском процессоре под той же Chimera Linux работает правильно.

Comments

[archaicos]

А на 10 и 0.1 работает?

А, это ещё и long double... Может быть и ABI проблема.

[vak]

На небольших экспонентах работает.

[archaicos]

А прочие функции с long double? Типа fabsl(), hypotl(), sqrtl(), truncl(), fmodl()?

[vak]

Об эти функции я пока не спотыкался. Разыщу гляну исходники этой библиотеки, может сразу всё станет понятно.

[ircicq]

Это вроде особенность "Apple Silion":
The long double type is a double precision
https://developer.apple.com/documentation/xcode/writing-arm64-code-for-apple-platforms

С этими ARM каждый как хочет так и коверкает архитектуру

[vak]

Это только в компиляторе clang под MacOS тип long double имеет размер 64 бита. На том же arm64 под Ubuntu оно 128 бит, а под Chimera Linux - 80 бит.

Перефразируя Венечку Ерофеева, чтобы не сумел загордиться программист, чтобы человек был грустен и растерян. 😀

[vak]

Нет, вру, под Химерой на arm64 тоже 128 бит. То есть совпадает с Ubuntu.

[chabapok]

у arm64 sizeof(long double) что выдает? 8, 10 или 16?

[vak]

Под MacOS выдаёт 8, под Линуксом 16.

[chabapok]

глюк же в макоси? или под линуксом на arm64 этот баг тоже имеет место?

Я не спец по ieee754, но помоему 64 бита это ничто для такого числа.

[vak]

Глюк под Химера Линуксом. На макоси компилятор не пропускает такие литералы. Экспонента выходит за допустимые пределы double.

[chabapok]

может старая и/или битая/неверно_собранная libm

ubuntu18.04 - нормально считает.

[vak]

На убунте libm от GNU. На химере - musl.

[juan_gandhi]

Оппа, вот это да.

[vak]

Это я полез разбираться, почему некоторые тесты X1 Algol не проходят.