По следам вопроса о битовом представлении вещественных чисел и моего ответа на него.
Хочу программно для любого вещественного типа определить, сколько бит в нём отводится под мантиссу, а сколько под экспоненту. Для этого написал следующий код (в нём бит под знак считается отдельно от мантиссы, поэтому числа на 1 меньше):
https://ideone.com/YuIWNc - код на Си (float, double, long double)
https://ideone.com/342B4S - код на Си++ (float, double, long double)
https://ideone.com/VURQnw - код на Си++ (float, double, long double, __float128)
#include <cstdio>
template <typename typed> void count(unsigned *result_m, unsigned *result_e)
{
typed x = 1, exp;
unsigned res, e;
for (res=0; x!=0; ++res) x/=2;
for (exp=1,e=0; exp*2<res; ++e) exp*=2;
*result_e = e+1;
*result_m = res-exp+1;
}
int main(void)
{
unsigned f_m, f_e, d_m, d_e, ld_m, ld_e, f128_m, f128_e;
count<float>(&f_m, &f_e);
count<double>(&d_m, &d_e);
count<long double>(&ld_m, &ld_e);
count<__float128>(&f128_m, &f128_e);
printf(" S M E SZ\n");
printf("float: 1 %3u %2u %3u\n", f_m, f_e, 8 * sizeof(float));
printf("double: 1 %3u %2u %3u\n", d_m, d_e, 8 * sizeof(double));
printf("long double: 1 %3u %2u %3u\n", ld_m, ld_e, 8 * sizeof(long double));
printf("__float128: 1 %3u %2u %3u\n", f128_m, f128_e, 8 * sizeof(__float128));
}
Получается так:
S M E SZ
float: 1 23 8 32
double: 1 52 11 64
long double: 1 63 15 128
__float128: 1 112 15 128
Для float
, double
и даже __float128
всё работает (Википедия, IEEE 754-2008).
А вот с long double
возникают проблемы:
1+63+15 = 79
- 79 бит. Вместо 80. Где ещё один бит?long double
представляет собой 10-байтные числа, ноsizeof
вернул 16.
Как-то можно получить 10?
unsigned
использован для счётчиков. 64-битных целочисленных переменных в коде нет вообще, битовых операций - тоже (/=2
и*=2
делаются над вещественными). Почему я считаю, что так можно, я описал в том ответе (старался описать подробно). Если так по какой-то причине нельзя, то хотелось бы эту причину увидеть ответом. – Qwertiy♦ 23 окт '17 в 14:35