delaunaytriangulation function: How to associate triangles with edges?

Question

Burak Pehlivan 2017-10-16

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/361530-delaunaytriangulation-function-how-to-associate-triangles-with-edges

评论： Precise Simulation 2017-10-21

Hi, I'm trying to perform finite volume method (CFD) calculations based on Matlab's built in "delaunaytriangulation" function. The function lists each triangle as a set of vertices or points. Similarly, the function lists each edge/face a set of vertices or points.

But I need each triangle information as a set of faces/edges to be able to loop over all edges for flux.

I get from dt.ConnectivityList command:

dt.ConnectivityList output

Point id_1 Point id_2 Point id_3 (Triangle no 1)

Point id_4 Point id_5 Point id_6 (Triangle no 2)

.... (And so on)

dt.Edge output

Point id_1 Point id_2 (Edge no 1)

POint id_3 Point id_4 (Edge no 2)

... (And so on)

What I need

Edge id_1 Edge id_2 Edge id_3 (Triangle no 1)

Edge id_4 Edge id_5 Edge id_6 (Triangle no 2)

.... (And so on)

Thank you very much.

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

请先登录，再进行评论。

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

Answer 1

KSSV 2017-10-16

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/361530-delaunaytriangulation-function-how-to-associate-triangles-with-edges#answer_285963

DelaunayTriangulation supports lot of functions....check this documentation for the functions it supports. I think edges should work for you.

https://in.mathworks.com/help/matlab/ref/delaunaytriangulation.html

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

KSSV 2017-10-17

Burak Pehlivan commented: Hi, edges command gives the vertex (id) sets of edges. I want the edge (id) sets of triangles.

请先登录，再进行评论。

Answer 2

Precise Simulation 2017-10-16

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/361530-delaunaytriangulation-function-how-to-associate-triangles-with-edges#answer_286016

编辑：Precise Simulation 2017-10-17

在 MATLAB Online 中打开

With the gridedge function included with the Matlab FEA Toolbox you can reconstruct and recover the grid edge numbering like this:

   p = rand(10,2);
   t = delaunay(p);
   [e,ev] = gridedge(t',2);
   help gridedge, e'

where e contains the edge numbering for each cell, and ev the vertices for each edge. Alternativley, the following code snippet should work with a different edge numbering:

   p = rand(10,2);
   t = delaunay(p);
                       
   e = []; e2c = [];
   for i=1:3
     e = [e; t(:,[i mod(i,3)+1])];
     e2c = [e2c; [i*ones(size(t,1),1) [1:size(t,1)]']];
   end
   [~,~,ei] = unique(sort(e,2),'rows');
   e = zeros(3,size(t,1));
   e(e2c(:,1)+3*(e2c(:,2)-1)) = ei;
   e'