Mercurial > repos > public > sbplib

comparison surfSym.m @ 550:e860670e72f1 feature/grids

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Merge default
author Jonatan Werpers <jonatan@werpers.com>
date Mon, 28 Aug 2017 14:29:46 +0200
parents bcddbc2beef4
children
comparison
equal deleted inserted replaced
549:ae905a11e32c 550:e860670e72f1
1 function h = surfSym(X,Y,Z)
2 if isvector(X) && isvector(Y)
3 [X,Y] = meshgrid(X,Y);
4 end
5
6 V_nodes = [X(:), Y(:), Z(:)];
7
8 X_centers = neighbourMean(X);
9 Y_centers = neighbourMean(Y);
10 Z_centers = neighbourMean(Z);
11 V_centers = [X_centers(:), Y_centers(:), Z_centers(:)];
12
13
14 N = prod(size(X));
15 nodeIndecies = reshape(1:N, size(X));
16 centerIndecies = reshape(N+(1:prod(size(X)-[1,1])), size(X) - [1,1]);
17
18 % figure()
19 % h = line(V_nodes(:,1),V_nodes(:,2),V_nodes(:,3));
20 % h.LineStyle = 'none';
21 % h.Marker = '.';
22 % h = line(V_centers(:,1),V_centers(:,2),V_centers(:,3));
23 % h.LineStyle = 'none';
24 % h.Marker = '.';
25 % h.Color = Color.red;
26 % axis equal
27
28
29 I_0 = nodeIndecies(1:end-1, 1:end-1);
30 I_1 = nodeIndecies(2:end, 1:end-1);
31 I_2 = nodeIndecies(2:end, 2:end);
32 I_3 = nodeIndecies(1:end-1, 2:end);
33 I_C = centerIndecies;
34
35 S.Vertices = [
36 V_nodes;
37 V_centers;
38 ];
39
40 S.Faces = [
41 I_0(:), I_1(:), I_C(:);
42 I_1(:), I_2(:), I_C(:);
43 I_2(:), I_3(:), I_C(:);
44 I_3(:), I_0(:), I_C(:);
45 ];
46
47 % figure()
48 h = patch(S, 'FaceVertexCData', S.Vertices(:,3),'FaceColor','flat');
49 end
50
51 % Calculate the mean of four neighbours around a patch
52 function M = neighbourMean(A)
53 M = (A(1:end-1, 1:end-1) + A(2:end, 1:end-1) + A(1:end-1, 2:end) + A(2:end, 2:end))/4;
54 end