Шутка QLineF::angleTo()

[Гурман]

Гоняю приложение в Релизе, и только сейчас заметил, что у меня кое-где не сработала оптимизация ортогональных маршрутов. Запускаю Креатор, в нём отладочную версию - а там всё работает, как надо. Недавно похожее было, но там с доступами к файлам библиотек и диагностикой были проблемы, но сейчас никаких библиотек нет и в помине, используется только QtCore. Проблему быстро обнаружил, она оказалась в функции вычисления угла между двумя линиями QLineF::angleTo(QLineF). В первой колонке таблицы значения, которые эта функция мне выдала при запуске отладочной версии, а во второй колонке - при запуске релиза:

0 180 0 0 180 0 0 180 0 0 180 0 0 180 0 0 180 0

360 180 360 360 180 360 360 180 360 3.50414e-15 180 3.50414e-15 3.50414e-15 180 3.50414e-15 3.50414e-15 180 3.50414e-15

Ну 3.5*10^-15 вместо 0, это ещё ладно, я этот момент заранее предусмотрел fuzzy-сравнением. Но вот разворот на 360 градусов (именно в этих случаях оптимизация не сработала), причём блин ровно  - это уже издевательство какое-то...  

[Fregloin]

что мешает сделать свою реализацию? лично я так и делал

[Гурман]

Да ничего не мешает, только не нужна она. Ну добавил проверку ещё и на 360 градусов, тоже fuzzy на всякий случай, и проблема ушла. Просто как-то это было ну ооочень неожиданно, и нормальной логике не поддаётся. Особенно смачно это смотрится по сравнению с дебаггером, где тихо всё, шито-крыто... А тема тут - на всякий случай, для тех, кто дальше упрётся фейсом об тейбл в похожих ситуациях. Я быстро в Сети ничего такого не нашёл.

[Авварон]

Зарепортьте багу с примером компилябельным. Я косяков не вижу, похоже оптимизатор шалит:

Код:

qreal QLineF::angleTo(const QLineF &l) const
{
    if (isNull() || l.isNull())
        return 0;

    const qreal a1 = angle();
    const qreal a2 = l.angle();

    const qreal delta = a2 - a1;
    const qreal delta_normalized = delta < 0 ? delta + 360 : delta;

    if (qFuzzyCompare(delta, qreal(360)))
        return 0;
    else
        return delta_normalized;
}

[Гурман]

Зарепортьте багу с примером компилябельным. Я косяков не вижу, похоже оптимизатор шалит:

Код:

delta < 0 ? delta + 360 : delta

Вот. Если delta условно "-0", например -3.50414e-15, то будет ответ 360, а при таком же положительном значении - это значение, то есть, "+0". Зависит уже от вычисления angle() каждого отрезка. Фишка в том, что отрезки могут быть разнонаправленными и с четырьмя разными углами к горизонтали. angle() получается немного разным в разных случаях, это же вещественное двойной точности, с неизбежным "шумом" младших разрядов. Чтобы пример получить, надо потратить время на охоту за нужными положениями отрезков. Сейчас у меня его нет, может быть на выходных.

Почему в отладчике всё ровно 0 - это другой вопрос... Это действительно похоже на баг.

[Авварон]

А, собственно бага-то вот она:

Код:

- if (qFuzzyCompare(delta, qreal(360)))
+ if (qFuzzyCompare(delta_normalized, qreal(360)))

[Авварон]

Можете привести корды линий, на к-ых бага происходит? Не могу подобрать:)

[Авварон]

https://codereview.qt-project.org/#/c/121039/

[Гурман]

А, собственно бага-то вот она:

Код:

- if (qFuzzyCompare(delta, qreal(360)))
+ if (qFuzzyCompare(delta_normalized, qreal(360)))

Не отвечает на вопрос - почему в отладчике работает правильно. Условная транслация не видна. Разве только если... отладочная библиотека собрана из другого исходника. 

[Гурман]

Можете привести корды линий, на к-ых бага происходит? Не могу подобрать:)

Чуть попозже, другие дела сейчас. Тоже ловить надо, где именно это вылазит. На сцене под сотню отрезков...

[Авварон]

Пока тестил, нашел еще одну "проблему" - qDebug печатает 359+1-0 (0 - очень малое) как 360, т.е. 359+1-0 < 360, но в консоли именно 360. Причем это "очень малое" ни фига не малое, на самом деле (оно больше чем e-5). То есть фаззи компар не всегда отрабатывает (см тест).

[Igors]

Код:

qreal QLineF::angleTo(const QLineF &l) const
{
    if (isNull() || l.isNull())
        return 0;

    const qreal a1 = angle();
    const qreal a2 = l.angle();

    const qreal delta = a2 - a1;
    const qreal delta_normalized = delta < 0 ? delta + 360 : delta;

    if (qFuzzyCompare(delta, qreal(360)))
        return 0;
    else
        return delta_normalized;
}

А, собственно бага-то вот она:

Код:

- if (qFuzzyCompare(delta, qreal(360)))
+ if (qFuzzyCompare(delta_normalized, qreal(360)))

Да, но это еще не все. В одну сторону округляем, в др нет. Поэтому возможна такая ситуация

Код

C++ (Qt)
a.angleTo(b)  // возвращает ноль  
b.angleTo(a)  // возвращает НЕ ноль 

Народная примета: если угол в градусах - реализация слабенькая. Fuzzy здесь вообще не к месту, вызывающему все равно самому надо округлять.

[Авварон]

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

Код:

if (a.angleTo(b)) == 0) ...

Вы напишете

Код:

if (qFuzzyIsNull(a.angleTo(b)))) ...

Зачем делать дважды одно и то же?

[Igors]

Вы напишете

Код:

if (qFuzzyIsNull(a.angleTo(b)))) ...

Из этого ничего хорошего не выйдет. Вообще в прикладной задаче есть свой epsilon, который может намного отличаться от используемого в Qt. Поэтому вероятно я напишу примерно

Код

C++ (Qt)
const qreal angleFudgeMin = 1.0e-5;   // может скопирую значение из qFuzzyIsNull, но может и нет  
const qreal angleFudgeMax = 360 - angleFudgeMin;
...
qreal angle = a.angleTo(b);
if (angle < angleFudgeMin || angle > angleFudgeMax) ...  // считаем что вектора совпадают 
 

   

[Авварон]

Ну вот, уже qAbs забыли:)