DelaunayTri still hangs or crashes MATLAB, any workaround?

Question

Illya 2014-12-29

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/168452-delaunaytri-still-hangs-or-crashes-matlab-any-workaround

评论： Illya 2014-12-29

Here is a simple test to reproduce the situation. The following example contains only 3 intersecting triangles.

    x = [ 15.5; 15.872042489999444; 14.0; 15.5;   
          16.35853698188107;  15.385547998117818; 14.5; 16.35853698188107;
          17.0; 15.388922001685851; 16.355162978313039; 17.0];
    y = [ 6.5; 5.0949970000008307; 4.5; 6.5;   
          4.9795712257177893; 5.210422774283872;  4.5; 4.9795712257177893;
          5.5; 5.2238170053246717; 4.9661769946769896; 5.5];
   con = [1 2; 2 3; 3 4; 5 6; 6 7; 7 8; 9 10; 10 11; 11 12];
   tri=DelaunayTri(x,y,con);  % after this call MATLAB never returns back 
                              % slowly eating-up memory until it eventually crashes.

Given the data-sets I'm dealing with such cases are not rare, ant it is almost impossible to detect/isolate them.

I tested it against R2010 and R2011 (win64), no success. Could anybody, please, run this test on newer versions?

Any clues are appreciated.

Thanks

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Sean de Wolski 2014-12-29

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/168452-delaunaytri-still-hangs-or-crashes-matlab-any-workaround#answer_163593

在 MATLAB Online 中打开

On my Win7x64 laptop, R2014b:

tic
x = [ 15.5; 15.872042489999444; 14.0; 15.5;
  16.35853698188107;  15.385547998117818; 14.5; 16.35853698188107;
  17.0; 15.388922001685851; 16.355162978313039; 17.0];
y = [ 6.5; 5.0949970000008307; 4.5; 6.5;
  4.9795712257177893; 5.210422774283872;  4.5; 4.9795712257177893;
  5.5; 5.2238170053246717; 4.9661769946769896; 5.5];
con = [1 2; 2 3; 3 4; 5 6; 6 7; 7 8; 9 10; 10 11; 11 12];
tri=DelaunayTri(x,y,con);  % after this call MATLAB never returns back
% slowly eating-up memory until it eventually
toc

Yields

Warning: Duplicate data points have been detected
and removed.
 The Triangulation indices and Constraints are
 defined with respect to the unique set of points
 in DelaunayTri property X. 
Warning: Intersecting edge constraints have been
split, this may have added new points into the
triangulation. 
Elapsed time is 0.001876 seconds.