Mercurial > repos > public > sbplib

comparison +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
comparison
equal deleted inserted replaced
1328:3a286d2d1939 1329:7df63b17e078
21 end 21 end
22 assert(matchingParams,'grid:boundaryOptimized:NonMatchingParameters','The number of parameters per dimensions do not match.'); 22 assert(matchingParams,'grid:boundaryOptimized:NonMatchingParameters','The number of parameters per dimensions do not match.');
23 23
24 % Check that stencil options are passed correctly (if supplied) 24 % Check that stencil options are passed correctly (if supplied)
25 if length(varargin) == n+2 % Stencil options supplied 25 if length(varargin) == n+2 % Stencil options supplied
26 availabe_opts = ["Accurate","accurate","A","Minimal","minimal","M"]; 26 availabe_opts = ["Accurate","accurate","A","acc","Minimal","minimal","M","min"];
27 assert(any(varargin{n+2} == availabe_opts), ... 27 assert(any(varargin{n+2} == availabe_opts), ...
28 'grid:boundaryOptimized:InvalidOption',"The operator option must be 'accurate' or 'minimal.'"); 28 'grid:boundaryOptimized:InvalidOption',"The operator option must be 'accurate' or 'minimal.'");
29 else %If not passed, populate varargin with default option 'accurate' 29 else %If not passed, populate varargin with default option 'accurate'
30 varargin(n+2) = {'accurate'}; 30 varargin(n+2) = {'accurate'};
31 end 31 end
32 32
33 % Specify generating function 33 % Specify generating function
34 switch varargin{n+2} 34 switch varargin{n+2}
35 case {'Accurate','accurate','A'} 35 case {'Accurate','accurate','A','acc'}
36 gridgenerator = @sbp.grid.accurateBoundaryOptimizedGrid; 36 gridgenerator = @sbp.grid.accurateBoundaryOptimizedGrid;
37 case {'Minimal','minimal','M'} 37 case {'Minimal','minimal','M','min'}
38 gridgenerator = @sbp.grid.minimalBoundaryOptimizedGrid; 38 gridgenerator = @sbp.grid.minimalBoundaryOptimizedGrid;
39 end 39 end
40 40
41 X = {}; 41 X = {};
42 h = []; 42 h = [];