annotate +multiblock/+domain/Rectangle.m @ 1031:2ef20d00b386 feature/advectionRV

For easier comparison, return both the first order and residual viscosity when evaluating the residual. Add the first order and residual viscosity to the state of the RungekuttaRV time steppers
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Thu, 17 Jan 2019 10:25:06 +0100
parents 4bb298faa8dc
children 375f73edbbd4 b0208b130880
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
590
7789a844ab65 Add stub for Rectangle
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
1 classdef Rectangle < multiblock.Definition
7789a844ab65 Add stub for Rectangle
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
2 properties
7789a844ab65 Add stub for Rectangle
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
3
593
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
4 blockTi % Transfinite interpolation objects used for plotting
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
5 xlims
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
6 ylims
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
7 blockNames % Cell array of block labels
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
8 nBlocks
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
9 connections % Cell array specifying connections between blocks
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
10 boundaryGroups % Structure of boundaryGroups
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
11
590
7789a844ab65 Add stub for Rectangle
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
12 end
7789a844ab65 Add stub for Rectangle
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
13
7789a844ab65 Add stub for Rectangle
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
14
7789a844ab65 Add stub for Rectangle
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
15 methods
593
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
16 % Creates a divided rectangle
590
7789a844ab65 Add stub for Rectangle
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
17 % x and y are vectors of boundary and interface positions.
593
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
18 % blockNames: cell array of labels. The id is default.
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
19 function obj = Rectangle(x,y,blockNames)
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
20 default_arg('blockNames',[]);
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
21
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
22 n = length(y)-1; % number of blocks in the y direction.
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
23 m = length(x)-1; % number of blocks in the x direction.
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
24 N = n*m; % number of blocks
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
25
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
26 if ~issorted(x)
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
27 error('The elements of x seem to be in the wrong order');
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
28 end
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
29 if ~issorted(flip(y))
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
30 error('The elements of y seem to be in the wrong order');
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
31 end
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
32
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
33 % Dimensions of blocks and number of points
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
34 blockTi = cell(N,1);
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
35 xlims = cell(N,1);
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
36 ylims = cell(N,1);
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
37 for i = 1:n
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
38 for j = 1:m
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
39 p1 = [x(j), y(i+1)];
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
40 p2 = [x(j+1), y(i)];
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
41 I = flat_index(m,j,i);
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
42 blockTi{I} = parametrization.Ti.rectangle(p1,p2);
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
43 xlims{I} = {x(j), x(j+1)};
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
44 ylims{I} = {y(i+1), y(i)};
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
45 end
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
46 end
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
47
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
48 % Interface couplings
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
49 conn = cell(N,N);
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
50 for i = 1:n
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
51 for j = 1:m
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
52 I = flat_index(m,j,i);
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
53 if i < n
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
54 J = flat_index(m,j,i+1);
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
55 conn{I,J} = {'s','n'};
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
56 end
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
57
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
58 if j < m
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
59 J = flat_index(m,j+1,i);
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
60 conn{I,J} = {'e','w'};
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
61 end
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
62 end
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
63 end
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
64
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
65 % Block names (id number as default)
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
66 if isempty(blockNames)
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
67 obj.blockNames = cell(1, N);
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
68 for i = 1:N
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
69 obj.blockNames{i} = sprintf('%d', i);
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
70 end
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
71 else
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
72 assert(length(blockNames) == N);
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
73 obj.blockNames = blockNames;
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
74 end
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
75 nBlocks = N;
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
76
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
77 % Boundary groups
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
78 boundaryGroups = struct();
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
79 nx = m;
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
80 ny = n;
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
81 E = cell(1,ny);
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
82 W = cell(1,ny);
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
83 S = cell(1,nx);
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
84 N = cell(1,nx);
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
85 for i = 1:ny
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
86 E_id = flat_index(m,nx,i);
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
87 W_id = flat_index(m,1,i);
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
88 E{i} = {E_id,'e'};
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
89 W{i} = {W_id,'w'};
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
90 end
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
91 for j = 1:nx
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
92 S_id = flat_index(m,j,ny);
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
93 N_id = flat_index(m,j,1);
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
94 S{j} = {S_id,'s'};
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
95 N{j} = {N_id,'n'};
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
96 end
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
97 boundaryGroups.E = multiblock.BoundaryGroup(E);
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
98 boundaryGroups.W = multiblock.BoundaryGroup(W);
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
99 boundaryGroups.S = multiblock.BoundaryGroup(S);
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
100 boundaryGroups.N = multiblock.BoundaryGroup(N);
596
4bb298faa8dc Add more boundary groups in Rectangle. Bugfix in boundaryGroups.all.
Martin Almquist <malmquist@stanford.edu>
parents: 593
diff changeset
101 boundaryGroups.all = multiblock.BoundaryGroup([E,W,S,N]);
4bb298faa8dc Add more boundary groups in Rectangle. Bugfix in boundaryGroups.all.
Martin Almquist <malmquist@stanford.edu>
parents: 593
diff changeset
102 boundaryGroups.WS = multiblock.BoundaryGroup([W,S]);
4bb298faa8dc Add more boundary groups in Rectangle. Bugfix in boundaryGroups.all.
Martin Almquist <malmquist@stanford.edu>
parents: 593
diff changeset
103 boundaryGroups.WN = multiblock.BoundaryGroup([W,N]);
4bb298faa8dc Add more boundary groups in Rectangle. Bugfix in boundaryGroups.all.
Martin Almquist <malmquist@stanford.edu>
parents: 593
diff changeset
104 boundaryGroups.ES = multiblock.BoundaryGroup([E,S]);
4bb298faa8dc Add more boundary groups in Rectangle. Bugfix in boundaryGroups.all.
Martin Almquist <malmquist@stanford.edu>
parents: 593
diff changeset
105 boundaryGroups.EN = multiblock.BoundaryGroup([E,N]);
593
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
106
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
107 obj.connections = conn;
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
108 obj.nBlocks = nBlocks;
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
109 obj.boundaryGroups = boundaryGroups;
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
110 obj.blockTi = blockTi;
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
111 obj.xlims = xlims;
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
112 obj.ylims = ylims;
590
7789a844ab65 Add stub for Rectangle
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
113
7789a844ab65 Add stub for Rectangle
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
114 end
7789a844ab65 Add stub for Rectangle
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
115
7789a844ab65 Add stub for Rectangle
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
116
7789a844ab65 Add stub for Rectangle
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
117 % Returns a multiblock.Grid given some parameters
593
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
118 % ms: cell array of [mx, my] vectors
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
119 % For same [mx, my] in every block, just input one vector.
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
120 function g = getGrid(obj, ms, varargin)
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
121
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
122 default_arg('ms',[21,21])
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
123
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
124 % Extend ms if input is a single vector
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
125 if (numel(ms) == 2) && ~iscell(ms)
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
126 m = ms;
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
127 ms = cell(1,obj.nBlocks);
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
128 for i = 1:obj.nBlocks
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
129 ms{i} = m;
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
130 end
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
131 end
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
132
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
133 grids = cell(1, obj.nBlocks);
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
134 for i = 1:obj.nBlocks
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
135 grids{i} = grid.equidistant(ms{i}, obj.xlims{i}, obj.ylims{i});
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
136 end
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
137
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
138 g = multiblock.Grid(grids, obj.connections, obj.boundaryGroups);
590
7789a844ab65 Add stub for Rectangle
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
139 end
7789a844ab65 Add stub for Rectangle
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
140
7789a844ab65 Add stub for Rectangle
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
141 % label is the type of label used for plotting,
7789a844ab65 Add stub for Rectangle
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
142 % default is block name, 'id' show the index for each block.
7789a844ab65 Add stub for Rectangle
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
143 function show(obj, label, gridLines, varargin)
593
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
144 default_arg('label', 'name')
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
145 default_arg('gridLines', false);
590
7789a844ab65 Add stub for Rectangle
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
146
593
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
147 if isempty('label') && ~gridLines
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
148 for i = 1:obj.nBlocks
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
149 obj.blockTi{i}.show(2,2);
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
150 end
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
151 axis equal
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
152 return
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
153 end
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
154
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
155 if gridLines
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
156 m = 10;
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
157 for i = 1:obj.nBlocks
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
158 obj.blockTi{i}.show(m,m);
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
159 end
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
160 end
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
161
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
162
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
163 switch label
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
164 case 'name'
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
165 labels = obj.blockNames;
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
166 case 'id'
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
167 labels = {};
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
168 for i = 1:obj.nBlocks
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
169 labels{i} = num2str(i);
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
170 end
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
171 otherwise
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
172 axis equal
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
173 return
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
174 end
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
175
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
176 for i = 1:obj.nBlocks
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
177 parametrization.Ti.label(obj.blockTi{i}, labels{i});
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
178 end
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
179
37948bfe9d79 Complete Rectangle definition
Martin Almquist <malmquist@stanford.edu>
parents: 590
diff changeset
180 axis equal
590
7789a844ab65 Add stub for Rectangle
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
181 end
7789a844ab65 Add stub for Rectangle
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
182
7789a844ab65 Add stub for Rectangle
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
183 % Returns the grid size of each block in a cell array
7789a844ab65 Add stub for Rectangle
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
184 % The input parameters are determined by the subclass
7789a844ab65 Add stub for Rectangle
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
185 function ms = getGridSizes(obj, varargin)
7789a844ab65 Add stub for Rectangle
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
186 end
7789a844ab65 Add stub for Rectangle
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
187 end
7789a844ab65 Add stub for Rectangle
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
188 end