Проверка на пересечение отрезков. Оптимизация

[_PROGRAMM_]

Добрый вечер. В моей программе одним из проблемных мест является проверка пересечения отрезков. В Интернете есть описание алгоритмов http://e-maxx.ru/algo/segments_intersection_checking . я использую первый. Только работать приходится с float. Переделал, получилось вот что:

#include <stdio.h>
#include <iostream>

using namespace std;

struct Point
{
   float x,y;
};

bool InOneLine(float a,float b, float c, float d)
{
     if (a>b) swap(a,b);           // Находятся ли наши отрезки на одной прямой 
     if (c>d) swap(c,d);
     return max(a,c)<=min(b,d);
}


float area(Point a, Point b, Point c)
{
     return ((a.x-c.x)*(b.y-c.y))-((b.x-c.x)*(a.y-c.y)); // поиск знаковой площади
}


bool Cross(Point a,Point b, Point c, Point d)
{
     if (InOneLine(a.x,b.x,c.x,d.x)&&InOneLine(a.y,b.y,c.y,d.y))
     {
         float a1 = area(a,b,c);  // находим площади
         float a2 = area(a,b,d);
         float a3 = area(c,d,a);
         float a4 = area(c,d,b);

         if(  (((int)a1|0x80000000)^((int)a2|0x80000000)==0)&& // Сравниваем float'ы
              (((int)a3|0x80000000)^((int)a4|0x80000000)==0)) 
               return true;     
         else
               return false;

     }
     else
     {
         return false;
     }
}

int main(int argc, char * argv[])
{
    Point a,b,c,d;
    a.x=-5;
    a.y=1;

    b.x=1;
    b.y=5;

    c.x=0;
    c.y=0;

    d.x=3;
    d.y=-6;




    if (Cross(a,b,c,d))
    {
       cout<<"Отрезки пересекаются";
    }
    else 
       cout<<"Отрезки не пересекаются";

    
    return 0;
}

Будет ли он оптимальнее(по скорости) алгоритма по ссылке? Выигрыш в один такт тоже важен(только если ущерб коду не огромен)

[ultimatet41]

вот в этом куске кода можно слегка упростить
с:

bool Cross(Point a,Point b, Point c, Point d)
{
     if (InOneLine(a.x,b.x,c.x,d.x)&&InOneLine(a.y,b.y,c.y,d.y))
     {
         float a1 = area(a,b,c);  // находим площади
         float a2 = area(a,b,d);
         float a3 = area(c,d,a);
         float a4 = area(c,d,b);

         if(  (((int)a1|0x80000000)^((int)a2|0x80000000)==0)&& // Сравниваем float'ы
              (((int)a3|0x80000000)^((int)a4|0x80000000)==0)) 
               return true;     
         else
               return false;

     }
     else
     {
         return false;
     }
}

на:

bool Cross(Point a, Point b, Point c, Point d)
{
     if (InOneLine(a.x, b.x, c.x, d.x) && InOneLine(a.y, b.y, c.y, d.y))
     {
         float a1 = area(a,b,c);  // находим площади
         float a2 = area(a,b,d);
         float a3 = area(c,d,a);
         float a4 = area(c,d,b);

         if(  (( (int) a1 | 0x80000000) ^ ( (int) a2 | 0x80000000 ) == 0) && // Сравниваем float'ы
              (( (int) a3 | 0x80000000) ^ ( (int) a4 | 0x80000000) == 0)) 
               return true;     
       }
     
     return false;
     
}

[_PROGRAMM_]

Спасибо.
Насколько это извращение:

 if(  (( (int) a1 | 0x80000000) ^ ( (int) a2 | 0x80000000 ) == 0) && // Сравниваем float'ы
              (( (int) a3 | 0x80000000) ^ ( (int) a4 | 0x80000000) == 0))

будет быстрее простого

if ((a1*a2<=0)&&(a3*a4<=0)

?

[Fenex]

Вряд ли кто сможет сказать наверняка. Проще написать два теста. Будет видно чётко и ясно.

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

[_PROGRAMM_]

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

Разве в моем компьютере числа кодируются по другому?

[DiemonStar]

Разве в моем компьютере числа кодируются по другому?

тут речь шла о скорости вычислений. если один из вариантов будет у вас работать быстрее, то не факт, что на других ПК будет так же...

[_PROGRAMM_]

Ну ведь в первом случае все равно придется использовать сопроцессор.

[DiemonStar]

Сообщение от _PROGRAMM_

Ну ведь в первом случае все равно придется использовать сопроцессор.

сопроцессор тоже не панацея. битовая логика + сдвиги с учетом конвейера отработают быстрее.

[Fenex]

Существуют вообще, казалось бы, парадоксальные вещи. Две команды на одном и том же процессоре в одном порядке могут выполняться быстрее, чем если их поменять местами. Производители процессоров пишут толстые книги, описывающие все тонкости работы с их детящами и возможности оптимизации.

[rrrFer]

bool InOneLine(float a,float b, float c, float d)
{
     if (a>b) swap(a,b);           // Находятся ли наши отрезки на одной прямой 
     if (c>d) swap(c,d);
     return max(a,c)<=min(b,d);
}

Че это? - комментарий левый, что делает функция? - что такое "одна прямая" вообще?

float area(Point a, Point b, Point c)
{
     return ((a.x-c.x)*(b.y-c.y))-((b.x-c.x)*(a.y-c.y)); // поиск знаковой площади
}

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

Насколько это извращение:

 if(  (( (int) a1 | 0x80000000) ^ ( (int) a2 | 0x80000000 ) == 0) && // Сравниваем float'ы
              (( (int) a3 | 0x80000000) ^ ( (int) a4 | 0x80000000) == 0))

будет быстрее простого

if ((a1*a2<=0)&&(a3*a4<=0)

Я не думаю что это одно и тоже, поэтому и сравнивать скорости их выполнения вобще никак нельзя (че сравнивать мопед с чайной чашкой?) И работать этот фрагмент быстро не будет, ИМХО, когда ты выполняешь приведение типов у тебя опять же создается временный объект (в твоем случае типа int)