GCC странно считает float

>В чём проблема?

В том, что ты используешь float, не зная об особенностях работы с данным типом данных на уровне процессора.

Попробуй double поюзать, чисто по приколу.

Выдаёт 2.00500. Примечательно, что если в цикле знак с < изменить на <= то нарезультат это не повлияет

Тогда, когда число должно быть равно 2.00, оно на самом деле чуточку меньше. Поэтому условие выполняется лишний раз.

А почему оно меньше? Как добиться более точных вычислений?

Если я не ошибаюсь число 5/1000 в двоичной системе представляет бесконечную периодическую дробь.

Есть такой стандарт - IEEE 754. Рекомендую ознакомиться.

Ах вот оно как. Спасибо, смысл понятен

>А почему оно меньше? Как добиться более точных вычислений?

Никогда не делай строгих сравнений с числами с плавающей запятой если для тебя чувствительно отклонение даже на бесконечно малую величину.

Всем спасибо, всё понял из этой статьи http://www.yur.ru/science/computer/IEEE754.htm

ты ошибаешься. КО намекает: двоичная система тут не причем, выше пруф есть от чего оно так (способ хранения в памяти).

и да, мне интересно как ты собрался делить 101 на 1111101000 (на самом деле интересно, я двоичную арифметику давно изучал, может забыл что).

тебе выше AptGet ответил

Есть такой стандарт - IEEE 754. Рекомендую ознакомиться.

>Всем спасибо, всё понял из этой статьи http://www.yur.ru/science/computer/IEEE754.htm

По ссылке фигня какая-то. Вот статья получше http://www.rsdn.ru/article/alg/float.xml#ELB

>ты ошибаешься. КО намекает: двоичная система тут не причем, выше пруф есть от чего оно так (способ хранения в памяти).

Я с капитаном не дружу. Давай развернутый ответ почему я не прав.

ты ошибаешься. КО намекает: двоичная система тут не причем, выше пруф есть от чего оно так (способ хранения в памяти).

Ты ошибаешься. Двоичная система тут тоже причём. Любую конечную двоичную дробь можно записать конечной десятичной дробью, но не всякую конечную десятичную дробь можно записать конечной двоичной дробью.

Кстати про IEEE 754. Если кому-нибудь нужно, могу отправить pdf'ку с IEEE 754-2008. Сам около месяца назад безуспешно искал на торрентах и рапидшарах, но нашёл только на каком-то китайском форуме в виде вложения, доступного только после регистрации =).

>IEEE 754-2008

Хммм.... Редкое коллекционное издание.

Если не трудно, скинь на ящик pathfinder.public@gmail.com

ок. поехали.

объясни мне что такое двоичная дробь. и как она храниться в памяти (таки на всех типах архитектур одинаково?).

насколько я понимаю - все зависит от конкретной реализации. что мешает тебе хранить конечную десятичную дробь в виде числа с фиксированной точкой (самописной реализацией к примеру). в чем я неправ?

возможно я и ошибаюсь, так ты просвети. возможно я чего-то не знаю.

ну я же спросил

и да, мне интересно как ты собрался делить 101 на 1111101000 (на самом деле интересно, я двоичную арифметику давно изучал, может забыл что).[/qoute]

http://kfe.fjfi.cvut.cz/~klimo/nm/ieee754.pdf
1985г

у 101 ведущие нули добавляются? блин, я что-то забыл просто, вот и вопрос задал. хотел попробовать поделить, в столбик )

По ссылке которую я дал есть пример с двоичной дробью.

>и да, мне интересно как ты собрался делить 101 на 1111101000 (на самом деле интересно, я двоичную арифметику давно изучал, может забыл что).

Арифметические действия во всех системах счисления выполняются одинаково. Делить столбиком, как в начальной школе учили.

спасибо, мое упущения. был неправ.

я поэтому и сказал что столбиком собираюсь делить, знаю я что такое двоичное деление. я про двоичные дроби то ли не знал, то ли забыл (не помню, может забыл просто).

*мое упущение (в какой-то момент хотел написать мои упущения похоже, что-то плохо соображаю сегодня)

101.00101 = 2^2 + 2^0 + 2^-3 + 2^-5.

Вот двоичная дробь, дальше сам поймёшь?

объясни мне что такое двоичная дробь.

То же, что и десятичная, например 1001100,010010011101001101.

и как она храниться в памяти (таки на всех типах архитектур одинаково?).

Архитектура и компьютеры вообще тут не причём, мы говорим про представление числа в разных системах счисления.

что мешает тебе хранить конечную десятичную дробь в виде числа с фиксированной точкой (самописной реализацией к примеру).

Потеря точности, причём при любой фиксированной длине дробной части.

В качестве эксперимента попробуй перевести число 0,3 в двоичную систему счисления.

слушай, а ты тут причем вообще? я вроде бы людям вопрос задавал, с которыми разговаривал. спасибо конечно, но мне ответили. к чему этот пост? *налетели коршуны мля*

>Кстати про IEEE 754. Если кому-нибудь нужно, могу отправить pdf'ку с IEEE 754-2008

М.б. вы выложите ее в открытый доступ, на slil.ru, например?

спс за ответ. был неправ.

> я вроде бы людям вопрос задавал, с которыми разговаривал

Это форум, да ещё и древовидный, да к тому же с показом комментариев в плоском виде сортированном по времени. Привыкай.

>Потеря точности, причём при любой фиксированной длине дробной части.

В качестве эксперимента попробуй перевести число 0,3 в двоичную систему счисления.

Что мешает хранить число в виде обыкновенной дроби? :)

М.б. вы выложите ее в открытый доступ, на slil.ru, например?

http://zalil.ru/29741926

Что мешает хранить число в виде обыкновенной дроби? :)

Это уже другой вопрос. =)

>http://zalil.ru/29741926

Спасибо.

да, наверное это особенность данного форума. и на самом деле непривычно просто. на будущее учту.

>http://zalil.ru/29741926

Спасибо

А почему оно меньше? Как добиться более точных вычислений?

Использовать рациональные вместо вещественных, например.

Или нужно сказать самому себе в первую очередь, что хотим получить таким алгоритмом и потом уже компилятору: [code=c] #include<stdio.h> #include<math.h> int main() { float a = 1; float h = 0.005; while ((fabs(a - 2.0) > fabs(a + h - 2.0))) a = a + h; printf(«%f\n», a); return 0; } [\code] Таким алгоритмом получают не 2, а наибольшее число X, которое ближе к 2, чем X+0.005.

[code=c] #include<stdio.h> #include<math.h> int main() { float a = 1; float h = 0.005; while ((fabs(a - 2.0) > fabs(a + h - 2.0))) a = a + h; printf(«%f\n», a); return 0; } [/code] Топикстартер предложил алгоритм, которым получают не 2, а наибольшее число X, которое ближе к 2, чем X+0.005. Если это то, что он хотел, то вот улучшенный вариант его алгоритма.

Сорри, [code] что-то заглючил. Но идея думаю понятна. В таком цикле лишний раз прибавление не произойдет. Мы с минимальной ошибкой вычислений окажемся около 2: ответом будет 1.9999999 или 2.0000001 - тут уже могут спорить знатоки двоичных представлений и архитектур.

юзай printf («%.3f\n», a);

Ну не в этом же трабла у топикстартера.

На самом деле, предлагаемые варианты типа:

а) использовать точное представление конечных десятичных дробей б) использовать точное представление рациональных чисел

правильные, но излишне масштабно решают проблему топикстартера, которая заключалась собственно в следующем:

В предложенном алгоритме какое условие остановки цикла использовать, чтоб полученное a было наибольшим числом, ближайшим к числу 2?

Условие остановки, предложенное топикстартером просто нужно заменить на предложенное анонимусом

fabs(a - 2) > fabs(a + 0.005 - 2)

что словами можно описать как:

«увеличивать a по 0.005 до тех пор, пока текущее значение не окажется ближе к 2, чем следующее».

При таком подходе a будет максимально близко после окончания вычислений к 2, и лишнего шага не произойдет.

Это вообще-то проходят на самых первых уроках курса численных методов или типа того.

http://ru.wikipedia.org/wiki/%D0%9C%D0%B0%D1%88%D0%B8%D0%BD%D0%BD%D1%8B%D0%B9...