Mercurial > repos > public > sbplib

comparison +parametrization/old/curve_discretise.m @ 248:81e0ead29431 feature/beams

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Fixed package names in +parametrization
author Jonatan Werpers <jonatan@werpers.com>
date Mon, 05 Sep 2016 16:34:22 +0200
parents 3a3cf386bb7e
children
comparison
equal deleted inserted replaced
247:c2ca9717db4d 248:81e0ead29431
31 31
32 % fprintf('ni = %d, nj = %d\n',ni,nj); 32 % fprintf('ni = %d, nj = %d\n',ni,nj);
33 33
34 if do_plot 34 if do_plot
35 fprintf('n:%d max: %f min: %f\n', n, max(d),min(d)); 35 fprintf('n:%d max: %f min: %f\n', n, max(d),min(d));
36 p = grid.map_curve(g,t); 36 p = parametrization.map_curve(g,t);
37 figure 37 figure
38 show(g,t,h); 38 show(g,t,h);
39 end 39 end
40 40
41 end 41 end
42 42
43 function [t,p,d] = curve_discretise_n(g,n) 43 function [t,p,d] = curve_discretise_n(g,n)
44 t = linspace(0,1,n); 44 t = linspace(0,1,n);
45 t = equalize_d(g,t); 45 t = equalize_d(g,t);
46 d = D(g,t); 46 d = D(g,t);
47 p = grid.map_curve(g,t); 47 p = parametrization.map_curve(g,t);
48 end 48 end
49 49
50 function d = D(g,t) 50 function d = D(g,t)
51 p = grid.map_curve(g,t); 51 p = parametrization.map_curve(g,t);
52 52
53 d = zeros(1,length(t)-1); 53 d = zeros(1,length(t)-1);
54 for i = 1:length(d) 54 for i = 1:length(d)
55 d(i) = norm(p(:,i) - p(:,i+1)); 55 d(i) = norm(p(:,i) - p(:,i+1));
56 end 56 end
69 end 69 end
70 end 70 end
71 71
72 72
73 function show(g,t,hh) 73 function show(g,t,hh)
74 p = grid.map_curve(g,t); 74 p = parametrization.map_curve(g,t);
75 75
76 76
77 77
78 h = grid.plot_curve(g); 78 h = parametrization.plot_curve(g);
79 h.LineWidth = 2; 79 h.LineWidth = 2;
80 axis equal 80 axis equal
81 hold on 81 hold on
82 h = plot(p(1,:),p(2,:),'.'); 82 h = plot(p(1,:),p(2,:),'.');
83 h.Color = [0.8500 0.3250 0.0980]; 83 h.Color = [0.8500 0.3250 0.0980];
84 h.MarkerSize = 24; 84 h.MarkerSize = 24;
85 hold off 85 hold off
86 86
87 n = grid.curve_normals(g,t); 87 n = parametrization.curve_normals(g,t);
88 hold on 88 hold on
89 for i = 1:length(t) 89 for i = 1:length(t)
90 p0 = p(:,i); 90 p0 = p(:,i);
91 p1 = p0 + hh*n(:,i); 91 p1 = p0 + hh*n(:,i);
92 l = [p0, p1]; 92 l = [p0, p1];