forum.boolean.name - библиотека пересечений

forum.boolean.name (http://forum.boolean.name/index.php)

- Библиотеки (http://forum.boolean.name/forumdisplay.php?f=28)

- - библиотека пересечений (http://forum.boolean.name/showthread.php?t=16494)

библиотека пересечений

Код:

initmath3d%(key$) - key$="boolean"  -в принципе можно

 и без этой функции



distance3d#(x#,y#,z#,x2#,y2#,z2#) - дистанция между двумя

 точками в 3d пространстве



lineoverlap2d%(ax1#,ay1#,ax2#,ay2#,bx1#,by1#,bx2#,by2#) - пересечение

 отрезков



trioverlap%(x#,y#,x1#,y1#,x2#,y2#,x3#,y3#) - принадлежность точки 

к треугольнику



quadoverlap%(x#,y#,x1#,y1#,x2#,y2#,x3#,y3#,x4#,y4#) - принадлежность 

точки к четырехугольнику



ovalsoverlap%(x1#,y1#,r1#,x2#,y2#,r2#)-пересечение двух

 окружностей

X#,Y# - координаты точки.
Переменные имеющие индекс(х1...хn) - координаты фигур

протестируйте кому не лень, или кому пригодится, кому что непонятно или баг, глюк или прочая нечисть - пишите

Ответ: библиотека пересечений

Я конечно не великий програмёр, но мне кажется что нет смысла в некоторых функциях

Например
---------------------
lineoverlap2d%(ax1#,ay1#,ax2#,ay2#,bx1#,by1#,bx2#, by2#)
вместо него использовать старый
RectsOverlap (X1,Y1,1,1,X2,Y2,rectWidth2,rectHeight2) - сделает ту же проверку точки к прямоугольнику
----------------
distance3d#(x#,y#,z#,x2#,y2#,z2#) - мне ни разу не пришлось брать "никкому" не относящиеся точки для проверки, и не представляю для чего это
Лучше старый велосипед entityDistance(ent1,ent2)
На крайняк пивота создать )
-------------------
А вот проверка отрезков, треугольника и окружностей вполне пригодятся
:)

Ответ: библиотека пересечений

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

Цитата:

Randomize

и что ты этим хочешь сказать?

конечно это как прикол

Код:

initmath3d%("boolean")



st=MilliSecs()

  For t=0 To 100000

      a=distance3db#(10,10,10,100,200,150)

  Next 

en=MilliSecs()

sum=en-st

DebugLog sum



st=MilliSecs()

  For t=0 To 100000

      a=distance3d#(10,10,10,100,200,150)

  Next 

en=MilliSecs()

sum=en-st

DebugLog sum



Stop 



Function distance3db#(x1#,y1#,z1#,x2#,y2#,z2#)

  Return Sqr((x2-x1)^2+(y2-y1)^2+(z2-z1)^2)

End Function

Результат
47
29

я конечно ни чего не хочу доказывать, но когда один раз нужно вычислить то это одно дело проще юзать функции,
а когда за цикл это вызывается много*много раз, то библиотека будет быстрее

Ответ: библиотека пересечений

Дистанция 3D:

Код:

Function Distance3D#(x1#,y1#,z1#,x2#,y2#,z2#)

  Return Sqr((x2-x1)^2+(y2-y1)^2+(z2-z1)^2)

End Function

Пересечение линий:

Код:

Global Intersection_X#

Global Intersection_Y#

Global Intersection_AB#

Global Intersection_CD#



Function LinesIntersect(Ax#, Ay#, Bx#, By#, Cx#, Cy#, Dx#, Dy#)

        Rn# = (Ay#-Cy#)*(Dx#-Cx#) - (Ax#-Cx#)*(Dy#-Cy#)

        Rd# = (Bx#-Ax#)*(Dy#-Cy#) - (By#-Ay#)*(Dx#-Cx#)



        If Rd# = 0        

                Return False ; Lines are parralel.        

        

        Else

                Sn# = (Ay#-Cy#)*(Bx#-Ax#) - (Ax#-Cx#)*(By#-Ay#)



                Intersection_AB# = Rn# / Rd#

                Intersection_CD# = Sn# / Rd#

                Intersection_X# = Ax# + Intersection_AB#*(Bx#-Ax#)

                Intersection_Y# = Ay# + Intersection_AB#*(By#-Ay#)

                        

                Return True

        EndIf

End Function

2Д точка об 2д треугольник

Код:

Function PointInTri(Px#, Py#, V0x#, V0y#, V1x#, V1y#, V2x#, V2y#)

        ; vector(e1,v1,v0)

        E1x# = V1x# - V0x#

        E1y# = V1y# - V0y#

        ; vector(e2,v2,v0)

        E2x# = V2x# - V0x#

        E2y# = V2y# - V0y#

        ; crossproduct(h,d,e2)

        Hx# = -E2y#

        Hy# =  E2x# 

        ; a = dotproduct(e1,h)

        A# = (E1x# * Hx#) + (E1y# * Hy#)        

        F# = 1.0 / A#        

        ; vector(s,p,v0)

        Sx# = Px# - V0x#

        Sy# = Py# - V0y#                

        ;u = f * (dotProduct(s,h))

        U# = F# * ((Sx# * Hx#) + (Sy# * Hy#))

        

        ; If the value of the U coordinate is outside the range of values inside the triangle,

        ; then the ray has intersected the plane outside the triangle.

        If (U# < 0) Or (U# > 1)

                Return False

        EndIf

        

        ; crossProduct(q,s,e1)

        Qz# = (Sx# * E1y#) - (E1x# * Sy#)



        ; v = f * dotProduct(d,q)

        V# = F# * Qz#

        

        ; If the value of the V coordinate is outside the range of values inside the triangle,

        ; then the ray has intersected the plane outside the triangle.

        If (V# < 0) Or (V# > 1) Then Return False



        ; U + V together cannot exceed 1.0 or the point is not in the triangle. 

        ; If you imagine the triangle as half a square this makes sense.  U=1 V=1 would be in the 

        ; lower left hand corner which would be in the second triangle making up the square.

        If (U# + V#) > 1 Then Return False



        ; The point was in the triangle. Yay!                

        Return True

                

        ; Note that you could also return the U and V coordinates calculated in this function

        ; if you need those values!



End Function

Столкновение с неповёрнутым прямоугольником

Код:

Function PointInRect(iPointX,iPointY,iXPos1,iYPos1,iXPos2,iYPos2)

        iXPos2=iXPos2 + iXPos1 

        iYPos2=iYPos2 + iYPos1

        Return ((((iPointX-iXPos1) Xor (iPointX-iXPos2)) And ((iPointY-iYPos1) Xor (iPointY-iYPos2))) And $80000000)

End Function

Твои овалы это не овалы, а обычная Distance2D:

Код:

Function Distance2D#(x1#,y1#,x2#,y2#)

Local nx#=x1-x2 ;Длина горизонтального Катета

Local ny#=y1-y2 ;Длина Вертикального Катета

Return Sqr((nx*nx)+(ny*ny)) ;Длина Гипотенузы

End Function

Ответ: библиотека пересечений

Ещё в середине 2007 писал подобные функции..
http://moka.ucoz.ru/load/0-0-0-27-20

Ответ: библиотека пересечений

Цитата:

Сообщение от Randomize (Сообщение 227931)

Дистанция 3D:

Код:

Function Distance3D#(x1#,y1#,z1#,x2#,y2#,z2#)

  Return Sqr((x2-x1)^2+(y2-y1)^2+(z2-z1)^2)

End Function

Цитата:

Сообщение от Spy4433 (Сообщение 227898)

Результат
47
29

А вы в курсе что в Блице оператор возведения в степень не оптимизируется (было бы логично оптимизировать т.к. степень в данном случае константная) и не инлайнится (т.е. всегда вызывается как функция), и раза в 3-4 медленее чем аналог умножением?
Spy4433
Разверни степени в умножение и потом сравни со своей библиотечной функцией.
Чтобы ты понял к чему я веду - нет смысла выносить в длл элементарный код.

ЗЫ
Ты оптимизации включал когда собирал библиотеку? Код ужасный, да и весит 386 кило - слишком толсто :)

Ответ: библиотека пересечений

Цитата:

Сообщение от Платон Александрович (Сообщение 227950)

Нет не в курсе. Даже после твоего поста. Даже больше - не хочу быть в курсе блица. Моей задачей было тупо вбросить функции которые уже давно существуют. Оптимизировать не так то и сложно.

Цитата:

initmath3d%(key$) - key$="boolean" -в принципе можно
и без этой функции

Меня это заставило найти и запостить аналоги.

Ответ: библиотека пересечений

Цитата:

Меня это заставило найти и запостить аналоги.

)))))))))))))))))))))))
эту функцию убрал почти сразу тк самого достала )

кстати насчет веса...
в библиотеке были подключены модули для работы с диалоговыми функциями и реализована работа, я их просто не описал в decls, а так библиотека весит около 90 кб