Ошибка с конечными симпти сегментами - PullRequest
Новая
       47

Ошибка с конечными симпти сегментами

1 голос
/ 28 сентября 2019

Я пересекаю плоскость с 12 сегментами, которые представляют края куба.У меня проблема в том, что все пересечения создают одну точку, даже если точка не должна быть частью отрезка, потому что они должны быть конечными, верно?

Код: 

from sympy import Point3D, Plane, intersection, Segment3D

# Cube's vertices
v = (Point3D(500, 500, 500), Point3D(-500, 500, 500), Point3D(-500, -500, 500),
     Point3D(500, -500, 500), Point3D(500, 500, -500), Point3D(-500, 500, -500),
     Point3D(-500, -500, -500), Point3D(500, -500, -500))
# Cube´s edges
a = (Segment3D(v[0], v[1]), Segment3D(v[1], v[2]),
     Segment3D(v[2], v[3]), Segment3D(v[3], v[0]),
     Segment3D(v[0], v[4]), Segment3D(v[1], v[5]),
     Segment3D(v[2], v[6]), Segment3D(v[3], v[7]),
     Segment3D(v[4], v[5]), Segment3D(v[5], v[6]),
     Segment3D(v[6], v[7]), Segment3D(v[7], v[4]))

# Example plane which should generate 3 points
plano = Plane(Point3D(450, 400, 400), Point3D(400, 450, 400), Point3D(400, 400, 450))

bad = []
good = []
for i in range(12):
    inter = intersection(plano, a[i])
    # This should be executed when the intersection generates anything, but is always executed:
    if inter:
        bad.append(inter[0])
        # This comparation should not be necessary, checks if point is in range desired
        if abs(inter[0][0]) <= 500 and abs(inter[0][1]) <= 500 and abs(inter[0][2]) <= 500: 
            good.append(inter[0])

print(len(bad), bad)
print(len(good), good)

Вывод: 

12 [Point3D(250, 500, 500), Point3D(-500, 1250, 500), Point3D(1250, -500, 500), Point3D(500, 250, 500), Point3D(500, 500, 250), Point3D(-500, 500, 1250), Point3D(-500, -500, 2250), Point3D(500, -500, 1250), Point3D(1250, 500, -500), Point3D(-500, 2250, -500), Point3D(2250, -500, -500), Point3D(500, 1250, -500)]
3 [Point3D(250, 500, 500), Point3D(500, 250, 500), Point3D(500, 500, 250)]

9 из 12 точек не являются частью какого-либо сегмента

Ответы [ 2 ]

1 голос
/ 28 сентября 2019

Это исправлено в мастере SymPy, но исправление еще не выпущено.По мастеру я получаю: 

3 [Point3D(250, 500, 500), Point3D(500, 250, 500), Point3D(500, 500, 250)]
3 [Point3D(250, 500, 500), Point3D(500, 250, 500), Point3D(500, 500, 250)]

(что, как я полагаю, то, что вы ожидаете)

0 голосов
/ 28 сентября 2019

В настоящее время это ошибка (в версии 1.4 или ранее), которая будет исправлена ​​в версии 1.5.https://github.com/sympy/sympy/pull/16637

...