Архив форумов ЦИТФорума

[r_o_o_t]

Hi everybody.

Unfortunately I cannot type in Russian but I can read.

I have changed all doubles to float in my c++ code which perform some scientific computations to increase the spead and memory requirements. And I've faced a very strange behaviour. In fact I cannot compare one float with another with != operator.

Here is the sample code which illustrates this situation:

#include <stdlib.h>

void main()
{
float a=0.0089f;
float b=0.0087f;
float zero=b-a+0.0002f; // do you think it would be zero?
// no, it is 3.34694e-010 )
if(b!=(a+0.0002f)) abort();
}

The compiler is MSVC.
Does anybody have any idea about this trouble?

Thanks in advance!

[r_o_o_t]

It's just a mistype, the question is the same

[CrazyWizard]

Насколько я понял, то не (b!=(a+0.0002f)), а (a!=(b+0.0002f)).
Или ошибка в другом?

[Kefir]

Да, здесь видать одни системщики. Привет от численников! В любом учебнике по программирования сказано что следует избегать сранений типа =, != для чисел с плавающей запятой. Из-за ошибок вычисления, т.к. в абсолютном большинстве случаев эти числа будут неравны. Более того от перестановки мест слагаемых сумма изменяется и еще много чего интересного. А если действительно нужно узнать равны ли числа, то нужно сделать примерно так:

[r_o_o_t]

Спасибо за ответ. Это всё, конечно, да. Но вот с даблом всё ж ведь работало как часы - там только сложение, вычитание и умножение на целые числа. И на флоат я перехожу для увеличения скорости. А если везде fabs(a-b)<... написать, то быстродействие - наоборот - упадёт.
И потом, если бы я, например, корень извлёк, а потом - в квадрат, то понятно, тут же - почти целые числа, больше чем четвёртого знака после запятой нет и быть не может.

[DmitryShm]

Используй длинную арифметику. На олимпиадах когда задачки на большие вычисления возникали, то все подготовленные школьники использовали длинную арифметику. Сейчас реализаций этих библиотек с длинными числами (очень длинными) немерянно. Или используй ядро вычислительной системы наконец. Например, Mathematica.
_________________
love IT

[Kefir]

С double работало потому как там точность выше. Вот и все объяснение. При переходе на float скорость скорее упадет, т.к. в FPU все вычисления один ...(в смысле, все равно) выполняются над double'ами. Если только компилятор не будет задействовать SSE. А будет он их задействовать или нет, кто его знает. Единственная выгода от использования float это уменьшение расхода памяти, что на мой взгляд, сегодня неактуально.

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

Увеличения скорости надо добиваться другими методами, например вынесение операций за пределы цикла...