Mercurial > repos > public > sbplib

annotate +grid/boundaryOptimized.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 e53b1e25970a
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
1296
2853b655c172 Fix typos in comments
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 1289
diff changeset
1 % Creates a Cartesian grid of dimension length(m)
2853b655c172 Fix typos in comments
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 1289
diff changeset
2 % over the domain xlim, ylim, ...
1289
2fd2e2337b77 Add utility function for constructing a (possibly multidimensional) grid based on the grid points used by the boundary optimized SBP operators
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
3 % The grid is non-equidistant in the boundary regions,
1296
2853b655c172 Fix typos in comments
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 1289
diff changeset
4 % with node placement based on boundary-optimized SBP operators.
1289
2fd2e2337b77 Add utility function for constructing a (possibly multidimensional) grid based on the grid points used by the boundary optimized SBP operators
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
5 % Examples:
1297
e53b1e25970a Change +sbp/+util/ to +sbp/+grid and change function names to camel case
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 1296
diff changeset
6 % g = grid.boundaryOptimized([mx, my], xlim, ylim, order, opt)
e53b1e25970a Change +sbp/+util/ to +sbp/+grid and change function names to camel case
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 1296
diff changeset
7 % g = grid.boundaryOptimized([10, 15], {0,1}, {0,2}, 4) - defaults to 'accurate' stencils
e53b1e25970a Change +sbp/+util/ to +sbp/+grid and change function names to camel case
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 1296
diff changeset
8 % g = grid.boundaryOptimized([10, 15], {0,1}, {0,2}, 4, 'minimal')
e53b1e25970a Change +sbp/+util/ to +sbp/+grid and change function names to camel case
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 1296
diff changeset
9 function g = boundaryOptimized(m, varargin)
1289
2fd2e2337b77 Add utility function for constructing a (possibly multidimensional) grid based on the grid points used by the boundary optimized SBP operators
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
10 n = length(m);
1296
2853b655c172 Fix typos in comments
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 1289
diff changeset
11
1289
2fd2e2337b77 Add utility function for constructing a (possibly multidimensional) grid based on the grid points used by the boundary optimized SBP operators
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
12 % Check that parameters matches dimensions
2fd2e2337b77 Add utility function for constructing a (possibly multidimensional) grid based on the grid points used by the boundary optimized SBP operators
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
13 matchingParams = false;
2fd2e2337b77 Add utility function for constructing a (possibly multidimensional) grid based on the grid points used by the boundary optimized SBP operators
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
14 if length(varargin) == n+1 % Minimal number of arguments
2fd2e2337b77 Add utility function for constructing a (possibly multidimensional) grid based on the grid points used by the boundary optimized SBP operators
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
15 matchingParams = iscell([varargin{1:n}]) && ...
2fd2e2337b77 Add utility function for constructing a (possibly multidimensional) grid based on the grid points used by the boundary optimized SBP operators
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
16 isfloat([varargin{n+1}]);
2fd2e2337b77 Add utility function for constructing a (possibly multidimensional) grid based on the grid points used by the boundary optimized SBP operators
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
17 elseif length(varargin) == n+2 % Stencil options supplied
2fd2e2337b77 Add utility function for constructing a (possibly multidimensional) grid based on the grid points used by the boundary optimized SBP operators
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
18 matchingParams = iscell([varargin{1:n}]) && ...
2fd2e2337b77 Add utility function for constructing a (possibly multidimensional) grid based on the grid points used by the boundary optimized SBP operators
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
19 isfloat([varargin{n+1}]) && ...
2fd2e2337b77 Add utility function for constructing a (possibly multidimensional) grid based on the grid points used by the boundary optimized SBP operators
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
20 ischar([varargin{n+2}]);
1296
2853b655c172 Fix typos in comments
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 1289
diff changeset
21 end
1297
e53b1e25970a Change +sbp/+util/ to +sbp/+grid and change function names to camel case
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 1296
diff changeset
22 assert(matchingParams,'grid:boundaryOptimized:NonMatchingParameters','The number of parameters per dimensions do not match.');
1289
2fd2e2337b77 Add utility function for constructing a (possibly multidimensional) grid based on the grid points used by the boundary optimized SBP operators
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
23
2fd2e2337b77 Add utility function for constructing a (possibly multidimensional) grid based on the grid points used by the boundary optimized SBP operators
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
24 % Check that stencil options are passed correctly (if supplied)
2fd2e2337b77 Add utility function for constructing a (possibly multidimensional) grid based on the grid points used by the boundary optimized SBP operators
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
25 if length(varargin) == n+2 % Stencil options supplied
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: 1297
diff changeset
26 availabe_opts = ["Accurate","accurate","A","acc","Minimal","minimal","M","min"];
1289
2fd2e2337b77 Add utility function for constructing a (possibly multidimensional) grid based on the grid points used by the boundary optimized SBP operators
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
27 assert(any(varargin{n+2} == availabe_opts), ...
1297
e53b1e25970a Change +sbp/+util/ to +sbp/+grid and change function names to camel case
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 1296
diff changeset
28 'grid:boundaryOptimized:InvalidOption',"The operator option must be 'accurate' or 'minimal.'");
1289
2fd2e2337b77 Add utility function for constructing a (possibly multidimensional) grid based on the grid points used by the boundary optimized SBP operators
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
29 else %If not passed, populate varargin with default option 'accurate'
2fd2e2337b77 Add utility function for constructing a (possibly multidimensional) grid based on the grid points used by the boundary optimized SBP operators
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
30 varargin(n+2) = {'accurate'};
2fd2e2337b77 Add utility function for constructing a (possibly multidimensional) grid based on the grid points used by the boundary optimized SBP operators
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
31 end
1296
2853b655c172 Fix typos in comments
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 1289
diff changeset
32
1289
2fd2e2337b77 Add utility function for constructing a (possibly multidimensional) grid based on the grid points used by the boundary optimized SBP operators
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
33 % Specify generating function
2fd2e2337b77 Add utility function for constructing a (possibly multidimensional) grid based on the grid points used by the boundary optimized SBP operators
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
34 switch varargin{n+2}
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: 1297
diff changeset
35 case {'Accurate','accurate','A','acc'}
1297
e53b1e25970a Change +sbp/+util/ to +sbp/+grid and change function names to camel case
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 1296
diff changeset
36 gridgenerator = @sbp.grid.accurateBoundaryOptimizedGrid;
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: 1297
diff changeset
37 case {'Minimal','minimal','M','min'}
1297
e53b1e25970a Change +sbp/+util/ to +sbp/+grid and change function names to camel case
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 1296
diff changeset
38 gridgenerator = @sbp.grid.minimalBoundaryOptimizedGrid;
1289
2fd2e2337b77 Add utility function for constructing a (possibly multidimensional) grid based on the grid points used by the boundary optimized SBP operators
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
39 end
2fd2e2337b77 Add utility function for constructing a (possibly multidimensional) grid based on the grid points used by the boundary optimized SBP operators
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
40
2fd2e2337b77 Add utility function for constructing a (possibly multidimensional) grid based on the grid points used by the boundary optimized SBP operators
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
41 X = {};
2fd2e2337b77 Add utility function for constructing a (possibly multidimensional) grid based on the grid points used by the boundary optimized SBP operators
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
42 h = [];
2fd2e2337b77 Add utility function for constructing a (possibly multidimensional) grid based on the grid points used by the boundary optimized SBP operators
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
43 for i = 1:n
2fd2e2337b77 Add utility function for constructing a (possibly multidimensional) grid based on the grid points used by the boundary optimized SBP operators
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
44 try
2fd2e2337b77 Add utility function for constructing a (possibly multidimensional) grid based on the grid points used by the boundary optimized SBP operators
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
45 [X{i},h(i)] = gridgenerator(varargin{i},m(i),varargin{n+1});
2fd2e2337b77 Add utility function for constructing a (possibly multidimensional) grid based on the grid points used by the boundary optimized SBP operators
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
46 catch exception % Propagate any errors in the grid generation functions.
2fd2e2337b77 Add utility function for constructing a (possibly multidimensional) grid based on the grid points used by the boundary optimized SBP operators
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
47 msgText = getReport(exception);
1297
e53b1e25970a Change +sbp/+util/ to +sbp/+grid and change function names to camel case
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 1296
diff changeset
48 error('grid:boundaryOptimized:InvalidParameter',msgText)
1289
2fd2e2337b77 Add utility function for constructing a (possibly multidimensional) grid based on the grid points used by the boundary optimized SBP operators
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
49 end
2fd2e2337b77 Add utility function for constructing a (possibly multidimensional) grid based on the grid points used by the boundary optimized SBP operators
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
50 end
2fd2e2337b77 Add utility function for constructing a (possibly multidimensional) grid based on the grid points used by the boundary optimized SBP operators
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
51
2fd2e2337b77 Add utility function for constructing a (possibly multidimensional) grid based on the grid points used by the boundary optimized SBP operators
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
52 g = grid.Cartesian(X{:});
2fd2e2337b77 Add utility function for constructing a (possibly multidimensional) grid based on the grid points used by the boundary optimized SBP operators
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
53 g.h = h;
2fd2e2337b77 Add utility function for constructing a (possibly multidimensional) grid based on the grid points used by the boundary optimized SBP operators
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
54 end