Mercurial > repos > public > sbplib

annotate +grid/boundaryOptimizedCurvilinear.m @ 1329:7df63b17e078 feature/D2_boundary_opt

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add support for boundary optimized grids in DefCurvilinear, and add boundaryOptimizedCurvilinear for constructing a curvilinear grid with boundary optimized grid point placement.
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Mon, 14 Feb 2022 14:55:29 +0100
parents
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
1329
7df63b17e078 Add support for boundary optimized grids in DefCurvilinear, and add boundaryOptimizedCurvilinear for constructing a curvilinear grid with boundary optimized grid point placement.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
1 % Creates a curvilinear grid of dimension length(m).
7df63b17e078 Add support for boundary optimized grids in DefCurvilinear, and add boundaryOptimizedCurvilinear for constructing a curvilinear grid with boundary optimized grid point placement.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
2 % over the logical domain xi_lim, eta_lim, ...
7df63b17e078 Add support for boundary optimized grids in DefCurvilinear, and add boundaryOptimizedCurvilinear for constructing a curvilinear grid with boundary optimized grid point placement.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
3 % The grid point distribution is based on the boundary
7df63b17e078 Add support for boundary optimized grids in DefCurvilinear, and add boundaryOptimizedCurvilinear for constructing a curvilinear grid with boundary optimized grid point placement.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
4 % optimized SBP operators, passed as arguments, order
7df63b17e078 Add support for boundary optimized grids in DefCurvilinear, and add boundaryOptimizedCurvilinear for constructing a curvilinear grid with boundary optimized grid point placement.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
5 % and opt.
7df63b17e078 Add support for boundary optimized grids in DefCurvilinear, and add boundaryOptimizedCurvilinear for constructing a curvilinear grid with boundary optimized grid point placement.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
6 % Examples:
7df63b17e078 Add support for boundary optimized grids in DefCurvilinear, and add boundaryOptimizedCurvilinear for constructing a curvilinear grid with boundary optimized grid point placement.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
7 % g = grid.boundaryOptimizedCurvilinear(mapping, [mx, my], xlim, ylim, order, opt)
7df63b17e078 Add support for boundary optimized grids in DefCurvilinear, and add boundaryOptimizedCurvilinear for constructing a curvilinear grid with boundary optimized grid point placement.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
8 % g = grid.boundaryOptimizedCurvilinear(mapping, [10, 15], {0,1}, {0,2}, 4) - defaults to 'accurate' stencils
7df63b17e078 Add support for boundary optimized grids in DefCurvilinear, and add boundaryOptimizedCurvilinear for constructing a curvilinear grid with boundary optimized grid point placement.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
9 % g = grid.boundaryOptimizedCurvilinear(mapping, [10, 15], {0,1}, {0,2}, 4, 'minimal')
7df63b17e078 Add support for boundary optimized grids in DefCurvilinear, and add boundaryOptimizedCurvilinear for constructing a curvilinear grid with boundary optimized grid point placement.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
10 function g = boundaryOptimizedCurvilinear(mapping, m, varargin)
7df63b17e078 Add support for boundary optimized grids in DefCurvilinear, and add boundaryOptimizedCurvilinear for constructing a curvilinear grid with boundary optimized grid point placement.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
11 n = length(m);
7df63b17e078 Add support for boundary optimized grids in DefCurvilinear, and add boundaryOptimizedCurvilinear for constructing a curvilinear grid with boundary optimized grid point placement.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
12 % TODO: The input parameter check is also present in boundaryOptimizedGrid.m,
7df63b17e078 Add support for boundary optimized grids in DefCurvilinear, and add boundaryOptimizedCurvilinear for constructing a curvilinear grid with boundary optimized grid point placement.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
13 % Consider refactoring it.
7df63b17e078 Add support for boundary optimized grids in DefCurvilinear, and add boundaryOptimizedCurvilinear for constructing a curvilinear grid with boundary optimized grid point placement.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
14
7df63b17e078 Add support for boundary optimized grids in DefCurvilinear, and add boundaryOptimizedCurvilinear for constructing a curvilinear grid with boundary optimized grid point placement.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
15 % Check that parameters matches dimensions
7df63b17e078 Add support for boundary optimized grids in DefCurvilinear, and add boundaryOptimizedCurvilinear for constructing a curvilinear grid with boundary optimized grid point placement.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
16 matchingParams = false;
7df63b17e078 Add support for boundary optimized grids in DefCurvilinear, and add boundaryOptimizedCurvilinear for constructing a curvilinear grid with boundary optimized grid point placement.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
17 if length(varargin) == n+1 % Minimal number of arguments
7df63b17e078 Add support for boundary optimized grids in DefCurvilinear, and add boundaryOptimizedCurvilinear for constructing a curvilinear grid with boundary optimized grid point placement.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
18 matchingParams = iscell([varargin{1:n}]) && ...
7df63b17e078 Add support for boundary optimized grids in DefCurvilinear, and add boundaryOptimizedCurvilinear for constructing a curvilinear grid with boundary optimized grid point placement.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
19 isfloat([varargin{n+1}]);
7df63b17e078 Add support for boundary optimized grids in DefCurvilinear, and add boundaryOptimizedCurvilinear for constructing a curvilinear grid with boundary optimized grid point placement.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
20 elseif length(varargin) == n+2 % Stencil options supplied
7df63b17e078 Add support for boundary optimized grids in DefCurvilinear, and add boundaryOptimizedCurvilinear for constructing a curvilinear grid with boundary optimized grid point placement.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
21 matchingParams = iscell([varargin{1:n}]) && ...
7df63b17e078 Add support for boundary optimized grids in DefCurvilinear, and add boundaryOptimizedCurvilinear for constructing a curvilinear grid with boundary optimized grid point placement.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
22 isfloat([varargin{n+1}]) && ...
7df63b17e078 Add support for boundary optimized grids in DefCurvilinear, and add boundaryOptimizedCurvilinear for constructing a curvilinear grid with boundary optimized grid point placement.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
23 ischar([varargin{n+2}]);
7df63b17e078 Add support for boundary optimized grids in DefCurvilinear, and add boundaryOptimizedCurvilinear for constructing a curvilinear grid with boundary optimized grid point placement.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
24 end
7df63b17e078 Add support for boundary optimized grids in DefCurvilinear, and add boundaryOptimizedCurvilinear for constructing a curvilinear grid with boundary optimized grid point placement.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
25 assert(matchingParams,'grid:boundaryOptimizedCurvilinear:NonMatchingParameters','The number of parameters per dimensions do not match.');
7df63b17e078 Add support for boundary optimized grids in DefCurvilinear, and add boundaryOptimizedCurvilinear for constructing a curvilinear grid with boundary optimized grid point placement.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
26
7df63b17e078 Add support for boundary optimized grids in DefCurvilinear, and add boundaryOptimizedCurvilinear for constructing a curvilinear grid with boundary optimized grid point placement.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
27 % Check that stencil options are passed correctly (if supplied)
7df63b17e078 Add support for boundary optimized grids in DefCurvilinear, and add boundaryOptimizedCurvilinear for constructing a curvilinear grid with boundary optimized grid point placement.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
28 if length(varargin) == n+2 % Stencil options supplied
7df63b17e078 Add support for boundary optimized grids in DefCurvilinear, and add boundaryOptimizedCurvilinear for constructing a curvilinear grid with boundary optimized grid point placement.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
29 availabe_opts = ["Accurate","accurate","A","acc","Minimal","minimal","M","min"];
7df63b17e078 Add support for boundary optimized grids in DefCurvilinear, and add boundaryOptimizedCurvilinear for constructing a curvilinear grid with boundary optimized grid point placement.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
30 assert(any(varargin{n+2} == availabe_opts), ...
7df63b17e078 Add support for boundary optimized grids in DefCurvilinear, and add boundaryOptimizedCurvilinear for constructing a curvilinear grid with boundary optimized grid point placement.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
31 'grid:boundaryOptimizedCurvilinear:InvalidOption',"The operator option must be 'accurate' or 'minimal.'");
7df63b17e078 Add support for boundary optimized grids in DefCurvilinear, and add boundaryOptimizedCurvilinear for constructing a curvilinear grid with boundary optimized grid point placement.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
32 else %If not passed, populate varargin with default option 'accurate'
7df63b17e078 Add support for boundary optimized grids in DefCurvilinear, and add boundaryOptimizedCurvilinear for constructing a curvilinear grid with boundary optimized grid point placement.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
33 varargin(n+2) = {'accurate'};
7df63b17e078 Add support for boundary optimized grids in DefCurvilinear, and add boundaryOptimizedCurvilinear for constructing a curvilinear grid with boundary optimized grid point placement.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
34 end
7df63b17e078 Add support for boundary optimized grids in DefCurvilinear, and add boundaryOptimizedCurvilinear for constructing a curvilinear grid with boundary optimized grid point placement.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
35
7df63b17e078 Add support for boundary optimized grids in DefCurvilinear, and add boundaryOptimizedCurvilinear for constructing a curvilinear grid with boundary optimized grid point placement.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
36 % Specify generating function
7df63b17e078 Add support for boundary optimized grids in DefCurvilinear, and add boundaryOptimizedCurvilinear for constructing a curvilinear grid with boundary optimized grid point placement.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
37 switch varargin{n+2}
7df63b17e078 Add support for boundary optimized grids in DefCurvilinear, and add boundaryOptimizedCurvilinear for constructing a curvilinear grid with boundary optimized grid point placement.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
38 case {'Accurate','accurate','A','acc'}
7df63b17e078 Add support for boundary optimized grids in DefCurvilinear, and add boundaryOptimizedCurvilinear for constructing a curvilinear grid with boundary optimized grid point placement.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
39 gridgenerator = @sbp.grid.accurateBoundaryOptimizedGrid;
7df63b17e078 Add support for boundary optimized grids in DefCurvilinear, and add boundaryOptimizedCurvilinear for constructing a curvilinear grid with boundary optimized grid point placement.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
40 case {'Minimal','minimal','M','min'}
7df63b17e078 Add support for boundary optimized grids in DefCurvilinear, and add boundaryOptimizedCurvilinear for constructing a curvilinear grid with boundary optimized grid point placement.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
41 gridgenerator = @sbp.grid.minimalBoundaryOptimizedGrid;
7df63b17e078 Add support for boundary optimized grids in DefCurvilinear, and add boundaryOptimizedCurvilinear for constructing a curvilinear grid with boundary optimized grid point placement.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
42 end
7df63b17e078 Add support for boundary optimized grids in DefCurvilinear, and add boundaryOptimizedCurvilinear for constructing a curvilinear grid with boundary optimized grid point placement.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
43
7df63b17e078 Add support for boundary optimized grids in DefCurvilinear, and add boundaryOptimizedCurvilinear for constructing a curvilinear grid with boundary optimized grid point placement.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
44 X = {};
7df63b17e078 Add support for boundary optimized grids in DefCurvilinear, and add boundaryOptimizedCurvilinear for constructing a curvilinear grid with boundary optimized grid point placement.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
45 h = [];
7df63b17e078 Add support for boundary optimized grids in DefCurvilinear, and add boundaryOptimizedCurvilinear for constructing a curvilinear grid with boundary optimized grid point placement.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
46 for i = 1:n
7df63b17e078 Add support for boundary optimized grids in DefCurvilinear, and add boundaryOptimizedCurvilinear for constructing a curvilinear grid with boundary optimized grid point placement.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
47 try
7df63b17e078 Add support for boundary optimized grids in DefCurvilinear, and add boundaryOptimizedCurvilinear for constructing a curvilinear grid with boundary optimized grid point placement.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
48 [X{i},h(i)] = gridgenerator(varargin{i},m(i),varargin{n+1});
7df63b17e078 Add support for boundary optimized grids in DefCurvilinear, and add boundaryOptimizedCurvilinear for constructing a curvilinear grid with boundary optimized grid point placement.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
49 catch exception % Propagate any errors in the grid generation functions.
7df63b17e078 Add support for boundary optimized grids in DefCurvilinear, and add boundaryOptimizedCurvilinear for constructing a curvilinear grid with boundary optimized grid point placement.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
50 msgText = getReport(exception);
7df63b17e078 Add support for boundary optimized grids in DefCurvilinear, and add boundaryOptimizedCurvilinear for constructing a curvilinear grid with boundary optimized grid point placement.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
51 error('grid:boundaryOptimizedCurvilinear:InvalidParameter',msgText)
7df63b17e078 Add support for boundary optimized grids in DefCurvilinear, and add boundaryOptimizedCurvilinear for constructing a curvilinear grid with boundary optimized grid point placement.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
52 end
7df63b17e078 Add support for boundary optimized grids in DefCurvilinear, and add boundaryOptimizedCurvilinear for constructing a curvilinear grid with boundary optimized grid point placement.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
53 end
7df63b17e078 Add support for boundary optimized grids in DefCurvilinear, and add boundaryOptimizedCurvilinear for constructing a curvilinear grid with boundary optimized grid point placement.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
54
7df63b17e078 Add support for boundary optimized grids in DefCurvilinear, and add boundaryOptimizedCurvilinear for constructing a curvilinear grid with boundary optimized grid point placement.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
55 g = grid.Curvilinear(mapping, X{:});
7df63b17e078 Add support for boundary optimized grids in DefCurvilinear, and add boundaryOptimizedCurvilinear for constructing a curvilinear grid with boundary optimized grid point placement.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
56 g.logic.h = h;
7df63b17e078 Add support for boundary optimized grids in DefCurvilinear, and add boundaryOptimizedCurvilinear for constructing a curvilinear grid with boundary optimized grid point placement.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
57 end