+36
- 01
- 02
- 03
- 04
- 05
- 06
- 07
- 08
- 09
- 10
- 11
- 12
- 13
- 14
- 15
- 16
- 17
- 18
- 19
- 20
- 21
- 22
- 23
- 24
- 25
- 26
- 27
- 28
- 29
- 30
- 31
- 32
- 33
- 34
- 35
- 36
- 37
- 38
- 39
- 40
- 41
- 42
- 43
- 44
- 45
- 46
- 47
- 48
- 49
- 50
- 51
- 52
- 53
- 54
- 55
- 56
- 57
- 58
- 59
- 60
- 61
- 62
- 63
- 64
- 65
- 66
- 67
- 68
- 69
- 70
- 71
- 72
- 73
- 74
- 75
- 76
- 77
- 78
- 79
- 80
- 81
- 82
- 83
- 84
- 85
- 86
- 87
- 88
- 89
- 90
- 91
- 92
- 93
#include <iostream>
#include <cmath>
struct Point3D
{
float x,y,z;
Point3D () {}
Point3D (float x, float y, float z) : x(x), y(y), z(z) {}
Point3D& operator -= (const Point3D& p) { x-=p.x; y-=p.y; z-=p.z; return *this; }
Point3D operator - (const Point3D& p) const { Point3D p2(*this); return (p2-=p); }
Point3D& operator *= (const float f) { x*=f; y*=f; z*=f; return *this; }
Point3D operator * (const float f) const { Point3D p2(*this); return (p2*=f); }
};
float Dot (const Point3D& p1, const Point3D& p2) { return p1.x*p2.x + p1.y*p2.y + p1.z*p2.z; }
struct Face
{
Point3D n;
float nc;
float Dist (const Point3D& p) const { return Dot(p,n)-nc; }
};
int show_float(float src)
{
union
{
int i;
float f;
} u;
u.f = src;
return u.i;
}
float from_int(int src)
{
union
{
int i;
float f;
} u;
u.i = src;
return u.f;
}
template<typename T>
T& operator<<(T& str, const Point3D& p)
{
str << std::hex << "Point3D(from_int(0x" << show_float(p.x) << "), from_int(0x" << show_float(p.y) << "), from_int(0x" << show_float(p.z) << "))";
return str;
}
struct SPoint
{
Point3D p;
bool DoCorrectFace(const Face& face)
{
bool correct = true;
float j=1.0f;
Point3D np=p;
for (;;)
{
float ad = face.Dist(np);
if (ad<=0.0f)
break;
correct=false;
np = p - (face.n*(ad*j));
j += 1.0f;
}
p=np;
return correct;
}
};
using namespace std;
int main()
{
cout << "Hello World!" << endl;
SPoint spoint;
spoint.p = Point3D(from_int(0x41c6940e), from_int(0x427352a6), from_int(0xc166e2d0));
cout << "Initial p:" << endl;
cout << spoint.p << endl;
cout << "Corrected:" << endl;
Face f;
f.n = Point3D(from_int(0x3d6cc83b), from_int(0x3f0e8841), from_int(0x3f5422bd));
f.nc = from_int(0x41bac3dc);
bool result = spoint.DoCorrectFace(f);
cout << spoint.p << endl;
cout << "Done: " << result << endl;
return 0;
}
говно в gcc
g++ (rev5, Built by MinGW-W64 project) 4.8.1
вывод в -O2 -DNDEBUG :
Hello World!
Initial p:
Point3D(from_int(0x41c6940e), from_int(0x427352a6), from_int(0xc166e2d0))
Corrected:
Point3D(from_int(0x41c6940e), from_int(0x427352a6), from_int(0xc166e2d1))
Done: 0
вывод в -O3 -DNDEBUG:
Hello World!
Initial p:
Point3D(from_int(0x41c6940e), from_int(0x427352a6), from_int(0xc166e2d0))
Corrected:
Point3D(from_int(0x41c6940e), from_int(0x427352a6), from_int(0xc166e2d0))
Done: 0
внимание вопрос: может ли быть такое, что DoCorrectFace не изменил точку ни на бит, но вернул false? В gcc может!
Запостил: 
TarasB,
23 Февраля 2014
Follow us!
-O2 -DNDEBUG и -O3 -DNDEBUG одинаково выдают хуйню (точка не изменилась после вызова), правильно точка обрабатывается вообще без -O
Это другой чел, который нарвался на баг в моём крестопроекте, которые не проявляется в студии.
Неделю пытались понять, почему такая хуйня.
Теперь поняли. Плавающийпитухпроблемы.
Сдается мне что там пахнет UB и gcc сделал тебе spoint.DoTrollFace(f).
>>Это другой чел, который нарвался на баг в моём крестопроекте, которые не проявляется в студии.
А-а. То есть по-прежнему непортируемый мс-онли код.
Мы вдвоём неделю искали, почему он оказался мс-онли.
http://artdesign21.narod.ru/image/1/ko_034.jpg
Все единички, баг повторить не удалось ;(
P.S. gcc любит оптимизировать операции с флоатами оставляя их в регистрах сопроцессора (вспомни мой пост о sin(x) != sin(x). Возможно проблема в этом. Как экстендеды они не равны, а когда меняешь и суешь во флоат - ниче не меняется.
UPD: А что делает DoCorrectFace?
Двигает точку p так, чтобы face.dist(p) стало <= 0.0f
Проблема в том, что на бумаге написать, как надо двигать, очень просто. А чтоб на компе было ок - тяжело.
Hello World!
Initial p:
Point3D(from_int(0x41c6940e), from_int(0x427352a6), from_int(0xc166e2d0))
Corrected:
Point3D(from_int(0x41c6940e), from_int(0x427352a6), from_int(0xc166e2d0))
Done: 0
Может это винды. На луникс влом ребутить.
> На луникс влом ребутить.
Вот на луниксе баг почему-то не получилось повторить.
P.S. Есть вероятность, что от вывода баг пройдет, как и тогда с синусами. Что очень хреново ;(
именно поэтому мы неделю и не могли его поймать
Тогда опция -ffloat-store по идее уберет баг.
Но все-таки зря ты сравниваешь с 0. Имхо, стоит взять какой-то маленький эпсилон и сравнивать с ним. Если один объект воткнется в другой на тысячную часть пикселя - никто не заметит.
Я-то думал, что уж на таких трюизмах-то можно забить на эпсилоны, а вот хуй, нельзя.
А вот в SSE флоаты и даблы настоящие, поэтому у меня баг не проявился ;(
А есть ли портабельный какой-нибидь
__pitux_gcc_builtin_set_fpy_precision?
Ну нету у FPU float'а, нету. Его можно только сэмулировать через выгрузку-загрузку в память. Что и делает -ffloat-store.
По идее портабельным решением будет volatile переменная. Компилятор не сможет убрать запись и чтение, и все будет норм.
P.S. Кстати баг, подобный твоему, когда-то был в пыхе. И он позволял специально подобранным числом повешать к хуям то ли парсер флоатов то ли конвертер флоата в строку.
А, ты вон про че... Я уже и забыл, что у штеудов такой флаг был...
И я его ставил в дельфийском хулионе. Для ускорения типа.
Кстати, вот очередной пример, когда глобальное состояние это очень хуёвая идея. Была бы эта точность в самой команде (например как в sse - разные инструкции для разной точности), а не в control word - таких проблем бы не было.
С другой стороны, создавать десяток почти одинаковых команд - усложнять декодер, жизнь кодерам, лишние байты в коде команды, особенно когда надо выполнить пачку - алгоритм. Оно ж для программеров на асме - байтоёбов по-опеределению.
И меня честно говоря уже раздражает насколько сильно распух ассемблер для последних интелов.
Смутно вспоминаю срачи, когда у интела была избыточная точность и она кому-то мешала, что они запилили сиё.
На сколько процентов ускоряет? Есть статистика?
Через год представляю картину: Тарас портирует хулион на ARM.
Джве недели постигали где баг армоблядском проце...
О сколько нам открытий чудных
Готовит плавальный питух.
И опыт, сын багов коварных,
И гцц - кодеров друг.
Интересно. А тут, на гкоде такое было?
Убирает.
> стоит взять какой-то маленький эпсилон и сравнивать с ним
+1
Зато int ad - да!
O_o. Но почему?
А ну хотя. Оно же не влияет на вычисления, которые приводят к кривому ad. Там тогда всю цепочку вычислений надо через volatile прогонять...
volatile float x,y,z; тоже. Мне кажется если ставить то там, что поставит КРЕСТ на скорости.
Последнее приводит только к изменению from_int(0xc166e2d0)) на from_int(0xc166e2d1)) - то есть как без оптимизаций, но done у меня и без оптимизаций false.
Не, оно не рекомендация. Просто точность теряется на минусе в dist(). Вот туда если добавляешь volatile переменные - баг проходит. Но это не решение проблемы ;(
А проблема тут в том, что у интеловского FPU нет флоатов.
Ну как это нет? Внутри нет, а чтобы сохранять/загружать из памяти есть.
Может printfы и конвертят таким образом в single precision? Мы же надеялись на то что volatile прогонит через память.
gcc 4.82 без O
На всех итерациях (штук 20): 32CFDDE8. На последней: B54DA3CE.
g++ -Wall -O3
Initial p:
Point3D(from_int(0x41c6940e), from_int(0x427352a6), from_int(0xc166e2d0))
Corrected:
B5551390
Point3D(from_int(0x41c6940e), from_int(0x427352a6), from_int(0xc166e2d0))
Done: 1
> На всех итерациях (штук 20): 32CFDDE8. На последней: B54DA3CE.
Именно так.
точка НЕ изменилась, но возвращаемый результат ноль говорит, что точка была некорректной, то есть её таки меняли
http://ideone.com/C5koEO
зацените возвращаемый результат
Я понял. -O1 и выше там нормально.
У меня Done: 0 всегда с любыми O (без принтов).
Единственное отличие from_int(0xc166e2d0)) и from_int(0xc166e2d1)).
Видать действительно в регистрах кешится, а принты каким-то образом очищают это.
Ну вспомни код, который я кидал давным-давно. Там два синуса давали разный результат. Если выводишь один из них - все равно багует. Выводишь оба - оп, и нету бага. Настоящий баг шрёдингера ;)
Можно ли считать, что мой крестоскилл за эти недели вырос так, как он не вырос за год до этого?
Сорри, его не я постил. Я сраный плагиатор ;(
Комментарии с анализом получили больше плюсов чем сам ГК.
("На луникс влом ребутить.": ставь виртуалбокс, ставь какой дебьян без гуя (или с xfce). 512М памяти и 20Г диска хватает за глаза. если хочешь совсем свежак ставь дебьян sid (aptosid), но его надо апдейтить часто (как минимум раз в месяц). и он ломается иногда когда миграции происходят - https://release.debian.org/transitions/ - но меня как вторичная/третичная система устраивает.)
надо epsilon добавить на строке 64 попробовать.
Ну вот у меня 64 битная ось. И, как я понимаю, наличие SSE на такой платформе подразумевается. Поэтому gcc собрал код через sse, и точность выше флоата не подскочила - бага нет.
У Тараса SSE нет (целерон 600), все считается на FPU, и проявляется баг.
да и вроде же gcc SSE без пинка не подключает. по крайней мере раньше надо было -march/-mcpu задавать - просто -О не хватало.
На 32 же по дефолту стоят весьма консервативные march/mcpu. Что-то типа 386 или 686. Поэтому там и надо указывать...
к сожалению у меня для 32бит нет окружения - с -m32 ничего не компилится вообще.
Аналогично, надо бы установить.
unknown software exception (0xc0000096) в приложении по адресу 0x00405025.
ВНЕЗАПНО
God made natural numbers; all else is the work of man
Leopold Kronecker
http://cache.desktopnexus.com/thumbnails/695809-bigthumbnail.jpg