c: некорректное значение в цикле вычислений

[qptr]

Перевод дробного десятичного числа в систему счисления Q. До этого был код с целой переменной и кастингом double в int, результат тот же. Теперь избавился от int вообще и написал так:

#include <stdio.h>
#include <stdlib.h>
#include <math.h>

#define Q 5

int main(int ac, char *av[]) {

  double x = 0.592;
  double y = 0;

  while(x != 0) {
    x = x * Q;
    y = floor(x);
    printf("x = %f, y = %f\n", x, y);
    x = x - y;
    printf("x - y = %f\n", x);
  }

return 0;

Вывод программы:


x = 2.960000, y = 2.000000
x - y = 0.960000
x = 4.800000, y = 4.000000
x - y = 0.800000
x = 4.000000, y = 3.000000 // ??????????? как так получается ?

[Mandrivnyk]

Все дело в том, что floor() возвращает максимальное целое число, которое не больше аргумента.
Возможно, из-за погрешности работы с числами с плавающей точкой, в переменной x хранится не 4.000000, а, например, 3.99999999998. Вот и получается.
Попробуйте прибавлять при вызове функции какую-нибудь мелкую константную величину.
И, кстати... из-за той же погрешности, x вряд ли когда-нибудь будет равен 0...

[qptr]

Ну и ну. Подзабыл я уже программирование за долгий перерыв, хотя догадки почему так происходит, были. Буду дома значит вкуривать в ieee floating point format.