fpc-3.0.0 - запилили баг в арифметику

[Сквозняк]

Раньше числа с плавающей запятой обрабатывались по борландовски, правильно - ячейки памяти с мусором усекались и в дальнейшем не использовались. Теперь же глючит, зато по стандартам - как в питонах и прочих рубях.

Код: Выделить всё

var
Q: double;

begin
Q:=0.1+0.2;
writeln(Q);
writeln(0.1+0.2);
end.

Результат тестовой программы скомпилированной компилятором 2.6.4 версии:

Код: Выделить всё

3.00000000000000E-001
3.0000000000000000E-0001


Результат тестовой программы скомпилированной компилятором 3.0.0 версии:

Код: Выделить всё

2.9999999999999999E-001
3.00000000000000000011E-0001


В новом компиляторе произошла потеря точности при операциях с числами с плавающей запятой. Кто-нибудь в курсе, это случайно произошло или кто-то сознательно портирует баги из других ЯП чтобы в паскале лучше чем там не было?

[Vadim]

Зато если сделать вот так:

Код: Выделить всё

var
  Q: extended;

begin
  Q:=0.1+0.2;
  writeln(Q);
  writeln(0.1+0.2);
end.

то результат получается просто замечательный:

3.00000000000000000011E-0001
3.00000000000000000011E-0001

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

[Сквозняк]

3.00000000000000000011E-0001

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

[Vadim]

А Вы посчитайте, сколько ноликов после запятой было в версии 2.6.4 (там, где у Вас константы, естественно) и сколько этих же ноликов в моём примере. Если Вы после этого не поменяете мнение на счет мусора, то я съем Вашу шляпу. В 2.6.4 Вы мусор просто не видели, а с типом ext - увидели, так что пока нет ничего сногсшибательного...

[Azazaz]

Числа с плавающей точкой - игрушка дьявола, ежжи.
Количество значащих цифр для real типов ограничено, за пределами этого ограничения находится мусор. Для double это – 15 цифр приблизительно, дальше уже идет что угодно, как эта четверка например. В бытовых ситуациях 15 цифр – хватает с головой. Для вызова сатаны можно использовать Extended 10 байтовый.

Добавлено спустя 18 минут 30 секунд:
Тип / размер(byte) / точность(значащих цифр)
float 4 7
double 8 15
extended - 10 19

[zub]

http://www.delphikingdom.com/asp/viewit ... alogid=374

[Сквозняк]

Раз Double "улучшили", то у нас для вещественных чисел осталось только 2 нормальных типа: Extended и Currency?

[zub]

Как раз Double и стоит использовать.

[Сквозняк]

Визуально содержимое нового Double вызывает желание его не использовать. Одно дело когда мусор где-то там за горизонтом и вылезает только в циклах, а другое, когда всё им засыпано.

[zub]

Незнаю про какое "нового Double" вы говорите. это дело давно стандартизовано, и новшеств в нем быть не может.

[Дож]

Больше похоже на то, что Writeln по-другому печатает Double'ы.

Чтобы детальней обсуждать именно арифметику, лучше сделать дамп содержимого переменной Q. То, что там выводит Writeln, да ещё и в десятичном формате, не так интересно.

[Сквозняк]

Незнаю про какое "нового Double" вы говорите. это дело давно стандартизовано, и новшеств в нем быть не может.

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

Добавлено спустя 14 минут 3 секунды:

Больше похоже на то, что Writeln по-другому печатает Double'ы.

Чтобы детальней обсуждать именно арифметику, лучше сделать дамп содержимого переменной Q. То, что там выводит Writeln, да ещё и в десятичном формате, не так интересно.

Именно десятичный формат и интересен - неинтересны не отображаемые последние разряды, неточность из которых выползет не скоро, если вообще успеет. Легко убедиться что поменялась общая логика обработки Double, а не только не во Writeln

Код: Выделить всё

var
S: string;
Q: double;

begin
Q:=0.1+0.2;
writeln(Q);
str(Q,S);
writeln(S);
end.

[zub]

Просветите, что за стандарт такой - борландовский?
Я вот в первом посте вижу повышение точности.
Попробуйте на бумажке написать как выглядят 0.1 0.2 0.3 с двоичной плавающей точкой?

Добавлено спустя 12 минут 40 секунд:
На лицо непонимание плавающей точки. Попробуйте такой вариант

Код: Выделить всё

var
S: string;
Q: double;
begin
Q:=0.1+0.2;
writeln(Q:1:14);
str(Q:1:14,S);
writeln(S);
readln;
end.

Добавлено спустя 56 секунд:
Сейчас на форуме неделя багоискателей. все дружно находят баги на ровном месте))

[Сквозняк]

Просветите, что за стандарт такой - борландовский?

То как Борланд реализовал в своих продуктах.

Попробуйте на бумажке написать как выглядят 0.1 0.2 0.3 с двоичной плавающей точкой?

Там будет не один вариант. Самый простой, это если не прятать, локализовать и т.д. баги. А другие - если этим заниматься.

Код: Выделить всё

var
Q: double;
E: extended;
L: int64;

begin
Q:=0.1+0.2;
L:=trunc(frac(Q)*100000000);
writeln(L);
E:=0.1+0.2;
L:=trunc(frac(E)*100000000);
writeln(L);
readln;
end.

Результат этого теста одинаковый для компиляторов 2.6.4 и 3.0.0 Из чего делается вывод что в extended баг двоичной записи числа отлично локализуется в независимости от того, как этот тип теперь реализован. К сожалению double показал худший результат, и в новом компиляторе его даже не пытаются подправить.

[zub]

Какието далекоидущие выводы. На основе чего? мне непонятно
Double одинаков в транке и в 2.6.4. Различие в процедуре округления к десятичному виду, т.е. только в том что вы видете на экране. в памяти и в прцессе вычислений всё идентично.

Код: Выделить всё

uses sysutils;
procedure DumpBytes(const name:string;const a;const len:integer);
var i:integer;
    pb:pbyte;
begin
  pb:=@a;
  write(format('%10s:',[name]));
  for i:=0 to len-1 do
  begin
    write(format('%3x',[pb^]));
    inc(pb);
  end;
  write(format(' general number:%g',[double(a)]));
  writeln;
end;
var
  q:double;
begin
  q:=0.1;
  DumpBytes('0.1',q,sizeof(q));
  q:=0.2;
  DumpBytes('0.2',q,sizeof(q));
  q:=0.3;
  DumpBytes('0.3',q,sizeof(q));
  q:=0.1+0.2;
  DumpBytes('0.1+0.2',q,sizeof(q));
  readln;
end.

транк:

Код: Выделить всё

0.1: 9A 99 99 99 99 99 B9 3F general number:0,10000000000000001
0.2: 9A 99 99 99 99 99 C9 3F general number:0,20000000000000001
0.3: 33 33 33 33 33 33 D3 3F general number:0,29999999999999999
0.1+0.2: 33 33 33 33 33 33 D3 3F general number:0,29999999999999999

2.6.4

Код: Выделить всё

0.1: 9A 99 99 99 99 99 B9 3F general number:0,10000000000000001
0.2: 9A 99 99 99 99 99 C9 3F general number:0,20000000000000001
0.3: 33 33 33 33 33 33 D3 3F general number:0,29999999999999999
0.1+0.2: 33 33 33 33 33 33 D3 3F general number:0,29999999999999999

привычные нам "короткие" дроби на подобии 0.1, 0.2, 0.3 не имеют в двоичном формате точного конечного представления. это вас и вводит в заблуждение