Mercurial > repos > public > sbplib

comparison +parametrization/Ti.m @ 272:9fe423849525 feature/beams

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Added evalfunction to Ti.
author Jonatan Werpers <jonatan@werpers.com>
date Tue, 06 Sep 2016 13:22:28 +0200
parents 81e0ead29431
children 875386d0b2ff
comparison
equal deleted inserted replaced
250:9628617bdf14 272:9fe423849525
34 end 34 end
35 35
36 obj.S = @S_fun; 36 obj.S = @S_fun;
37 end 37 end
38 38
39 % Does this funciton make sense?
40 % Should it always be eval?
39 function [X,Y] = map(obj,u,v) 41 function [X,Y] = map(obj,u,v)
40 default_arg('v',u); 42 default_arg('v',u);
41 43
42 if isscalar(u) 44 if isscalar(u)
43 u = linspace(0,1,u); 45 u = linspace(0,1,u);
60 62
61 for i = 1:nv 63 for i = 1:nv
62 p = S(u,v(i)); 64 p = S(u,v(i));
63 X(i,:) = p(1,:); 65 X(i,:) = p(1,:);
64 Y(i,:) = p(2,:); 66 Y(i,:) = p(2,:);
67 end
68 end
69
70 % Evaluate S for each pair of u and v,
71 % Return same shape as u
72 function [x, y] = eval(obj, u, v)
73 x = zeros(size(u));
74 y = zeros(size(u));
75
76 for i = 1:numel(u)
77 p = obj.S(u(i), v(i));
78 x(i) = p(1,:);
79 y(i) = p(2,:);
65 end 80 end
66 end 81 end
67 82
68 function h = plot(obj,nu,nv) 83 function h = plot(obj,nu,nv)
69 S = obj.S; 84 S = obj.S;