Сообщество программистов Autodesk в СНГ

ADN Club => AutoLisp / VisualLISP и DCL => Тема начата: Macondo от 30-08-2017, 23:49:43

Название: Exception occurred: 0xC0000005 (Access Violation): alloc, expand или что делать?
Отправлено: Macondo от 30-08-2017, 23:49:43
На чертеже имеется 3000 отрезков (может быть и больше). Нужно найти все точки их пересечений.
Использую примерно такой код с немного преобразованными функциями от Lee Mac (inters вместо intersectwith). Для 3000 отрезков получается около 5 миллионов запусков inters, в результате всё висит, пока не появляется ошибка
; ошибка: Возникло исключение: 0xC0000005 (Нарушение доступа)

Код - Auto/Visual Lisp [Выбрать]
  1. (setq ss nil ss (ssget "_X" '((0 . "LINE")))); набор всех отрезков на чертеже
  2. (if (not ss) (setq ss (ssadd)))
  3. (setq all_lines_list (vl-remove-if-not (function (lambda (x) (= (type x) 'ename))) (mapcar 'cadr (ssnamex ss)))); список всех отрезков
  4. (setq ss nil vla_line_list nil)
  5. (foreach item all_lines_list (setq vla_line_list (append vla_line_list (list (vlax-ename->vla-object item))))); список всех vla-объектов-отрезков
  6. (setq all_inters_points nil all_inters_points (LM:intersectionsinobjlist vla_line_list)); поиск всех точек пересечения
  7.  
  8. (defun LM:intersectionsinobjlist (lst / ob1 rtn )
  9.     (while (setq ob1 (car lst))
  10.         (foreach ob2 (setq lst (cdr lst))
  11.             (setq rtn (cons (LM:intersections ob1 ob2) rtn))
  12.         )
  13.     )
  14.     (apply 'append (reverse rtn))
  15. );defun
  16.  
  17. (defun LM:intersections (ob1 ob2 / lst rtn )
  18.     (setq lst (inters (vlax-curve-getStartPoint ob1) (vlax-curve-getEndPoint ob1) (vlax-curve-getStartPoint ob2) (vlax-curve-getEndPoint ob2)))
  19.     (repeat (/ (length lst) 3)
  20.         (setq rtn (cons (list (car lst) (cadr lst) (caddr lst)) rtn) lst (cdddr lst))
  21.     )
  22.     (reverse rtn)
  23. );defun

На одном из форумов видел сообщение от Александра Ривилиса об использовании функции expand. Не очень понял, за что именно отвечает expand и alloc (и как их подбирать), но значения вроде 1000 для alloc и 10000 для expand проблему не решили. Что тут можно ещё сделать?
Название: Re: Exception occurred: 0xC0000005 (Access Violation): alloc, expand или что делать?
Отправлено: Александр Ривилис от 30-08-2017, 23:56:06
Цитата: Macondo от 30-08-2017, 23:49:43
Что тут можно ещё сделать?
Радикально изменить алгоритм.
Название: Re: Exception occurred: 0xC0000005 (Access Violation): alloc, expand или что делать?
Отправлено: Macondo от 31-08-2017, 00:01:28
При поиске точек пересечения всех отрезков друг с другом в любом случае потребуется перебрать все комбинации отрезков. Сам код несложный, но объём операций огромный. Как его упростить, в какую сторону "копать"?
Название: Re: Exception occurred: 0xC0000005 (Access Violation): alloc, expand или что делать?
Отправлено: Александр Ривилис от 31-08-2017, 00:07:06
1. Не нужны здесь vla-/vlax-методы. Для отрезков достаточно получать dxf-группы 10 и 11
2. Что это за заумное выражение:
Код - Auto/Visual Lisp [Выбрать]
  1.     (repeat (/ (length lst) 3)
  2.         (setq rtn (cons (list (car lst) (cadr lst) (caddr lst)) rtn) lst (cdddr lst))
  3.     )
???
3. Можно перед проверкой на пересечение отрезков проверять факт пересечения BoundingBox обоих отрезков. Если они не пересекаются, то и отрезки не могут пересекаться.
Название: Re: Exception occurred: 0xC0000005 (Access Violation): alloc, expand или что делать?
Отправлено: Macondo от 31-08-2017, 01:53:48
1. Изменил vla- на (cdr (assoc (entget...
2. Заумное выражение сделано для полилиний со множественными пересечениями, тут действительно лишнее.
После этих двух исправлений ошибка исчезла, хотя время ожидания составляет около 5 минут.
А как BoundingBox может помочь? Ведь это дополнительный поиск пересечений, а сейчас один только inters?
Название: Re: Exception occurred: 0xC0000005 (Access Violation): alloc, expand или что делать?
Отправлено: Владимир Шу от 31-08-2017, 07:27:16
Цитата: Macondo от 31-08-2017, 01:53:48
А как BoundingBox может помочь? Ведь это дополнительный поиск пересечений, а сейчас один только inters?

Попробуйте разделить весь массив линий на несколько зон и потом искать пересечения только в линиях внутри зоны, а вообще 
Цитата: Александр Ривилис от 30-09-2016, 15:36:25
тут на форуме обсуждалось использование R-tree для быстрого поиска.
Название: Re: Exception occurred: 0xC0000005 (Access Violation): alloc, expand или что делать?
Отправлено: Александр Ривилис от 31-08-2017, 10:51:07
Цитата: Macondo от 31-08-2017, 01:53:48
А как BoundingBox может помочь? Ведь это дополнительный поиск пересечений, а сейчас один только inters?
Пересечении в данном случае это проверки пересечения диапазонов [X1min,X1max] с [X2min,X2max] и [Y1min,Y1max] с [Y2min,Y2max]. Если пересечений диапазонов нет, то нечего и проверять пересечение самих отрезков при помощи inters.
Кроме того, я подозреваю что ты вызвал каждый раз в цикле (entget) для каждого из отрезков. А нужно один раз пройтись по всем отрезкам и создать список вида '((P1начала P1конца) ... (PNначала PNконца)) и тогда находить пересечения будет намного проще и быстрее.