Mercurial > repos > public > sbplib

annotate +sbp/+grid/minimalBoundaryOptimizedGrid.m @ 1297:e53b1e25970a feature/boundary_optimized_grids

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Change +sbp/+util/ to +sbp/+grid and change function names to camel case
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Tue, 07 Jul 2020 16:08:08 +0200
parents +sbp/+util/minimalBoundaryOptimizedGrid.m@e059a43bb675
children 0ffb5bfa65e4
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
1287
38653d26225c Make accurate/minimalBoundaryOptimizedGrid take the domain limits as input
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 1285
diff changeset
1 function [x,h] = minimalBoundaryOptimizedGrid(lim,N,order)
1288
e059a43bb675 Error-check format of limit arguments
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 1287
diff changeset
2 assert(iscell(lim) && numel(lim) == 2,'The limits should be cell arrays with 2 elements.');
1287
38653d26225c Make accurate/minimalBoundaryOptimizedGrid take the domain limits as input
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 1285
diff changeset
3 L = lim{2} - lim{1};
1288
e059a43bb675 Error-check format of limit arguments
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 1287
diff changeset
4 assert(L>0,'Limits must be given in increasing order.');
1285
6b68f939d023 Add utility functions for constructing the grids used by the boundary optimized D1Nonequidistant opsets
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
5 %%%% Non-equidistant grid points %%%%%
6b68f939d023 Add utility functions for constructing the grids used by the boundary optimized D1Nonequidistant opsets
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
6 xb = boundaryPoints(order);
6b68f939d023 Add utility functions for constructing the grids used by the boundary optimized D1Nonequidistant opsets
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
7 m = length(xb)-1; % Number of non-equidistant points
6b68f939d023 Add utility functions for constructing the grids used by the boundary optimized D1Nonequidistant opsets
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
8 assert(N-2*(m+1)>=0,'Not enough grid points to contain the boundary region. Requires at least %d points.',2*(m+1));
6b68f939d023 Add utility functions for constructing the grids used by the boundary optimized D1Nonequidistant opsets
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
9 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
6b68f939d023 Add utility functions for constructing the grids used by the boundary optimized D1Nonequidistant opsets
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
10
6b68f939d023 Add utility functions for constructing the grids used by the boundary optimized D1Nonequidistant opsets
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
11 %%%% Compute h %%%%%%%%%%
6b68f939d023 Add utility functions for constructing the grids used by the boundary optimized D1Nonequidistant opsets
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
12 h = L/(2*xb(end) + N-1-2*m);
6b68f939d023 Add utility functions for constructing the grids used by the boundary optimized D1Nonequidistant opsets
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
13 %%%%%%%%%%%%%%%%%%%%%%%%%
6b68f939d023 Add utility functions for constructing the grids used by the boundary optimized D1Nonequidistant opsets
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
14
6b68f939d023 Add utility functions for constructing the grids used by the boundary optimized D1Nonequidistant opsets
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
15 %%%% Define grid %%%%%%%%
6b68f939d023 Add utility functions for constructing the grids used by the boundary optimized D1Nonequidistant opsets
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
16 x = h*[xb; linspace(xb(end)+1,L/h-xb(end)-1,N-2*(m+1))'; L/h-flip(xb) ];
1287
38653d26225c Make accurate/minimalBoundaryOptimizedGrid take the domain limits as input
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 1285
diff changeset
17 x = x + lim{1};
1285
6b68f939d023 Add utility functions for constructing the grids used by the boundary optimized D1Nonequidistant opsets
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
18 %%%%%%%%%%%%%%%%%%%%%%%%%
6b68f939d023 Add utility functions for constructing the grids used by the boundary optimized D1Nonequidistant opsets
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
19 end
6b68f939d023 Add utility functions for constructing the grids used by the boundary optimized D1Nonequidistant opsets
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
20
6b68f939d023 Add utility functions for constructing the grids used by the boundary optimized D1Nonequidistant opsets
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
21 function xb = boundaryPoints(order)
6b68f939d023 Add utility functions for constructing the grids used by the boundary optimized D1Nonequidistant opsets
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
22 switch order
6b68f939d023 Add utility functions for constructing the grids used by the boundary optimized D1Nonequidistant opsets
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
23 case 4
6b68f939d023 Add utility functions for constructing the grids used by the boundary optimized D1Nonequidistant opsets
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
24 x0 = 0.0000000000000e+00;
6b68f939d023 Add utility functions for constructing the grids used by the boundary optimized D1Nonequidistant opsets
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
25 x1 = 7.7122987842562e-01;
6b68f939d023 Add utility functions for constructing the grids used by the boundary optimized D1Nonequidistant opsets
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
26 xb = [x0 x1]';
6b68f939d023 Add utility functions for constructing the grids used by the boundary optimized D1Nonequidistant opsets
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
27 case 6
6b68f939d023 Add utility functions for constructing the grids used by the boundary optimized D1Nonequidistant opsets
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
28 x0 = 0.0000000000000e+00;
6b68f939d023 Add utility functions for constructing the grids used by the boundary optimized D1Nonequidistant opsets
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
29 x1 = 4.0842950991998e-01;
6b68f939d023 Add utility functions for constructing the grids used by the boundary optimized D1Nonequidistant opsets
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
30 x2 = 1.1968523189207e+00;
6b68f939d023 Add utility functions for constructing the grids used by the boundary optimized D1Nonequidistant opsets
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
31 xb = [x0 x1 x2]';
6b68f939d023 Add utility functions for constructing the grids used by the boundary optimized D1Nonequidistant opsets
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
32 case 8
6b68f939d023 Add utility functions for constructing the grids used by the boundary optimized D1Nonequidistant opsets
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
33 x0 = 0.0000000000000e+00;
6b68f939d023 Add utility functions for constructing the grids used by the boundary optimized D1Nonequidistant opsets
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
34 x1 = 4.9439570885261e-01;
6b68f939d023 Add utility functions for constructing the grids used by the boundary optimized D1Nonequidistant opsets
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
35 x2 = 1.4051531374839e+00;
6b68f939d023 Add utility functions for constructing the grids used by the boundary optimized D1Nonequidistant opsets
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
36 xb = [x0 x1 x2]';
6b68f939d023 Add utility functions for constructing the grids used by the boundary optimized D1Nonequidistant opsets
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
37 case 10
6b68f939d023 Add utility functions for constructing the grids used by the boundary optimized D1Nonequidistant opsets
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
38 x0 = 0.0000000000000e+00;
6b68f939d023 Add utility functions for constructing the grids used by the boundary optimized D1Nonequidistant opsets
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
39 x1 = 5.8556160757529e-01;
6b68f939d023 Add utility functions for constructing the grids used by the boundary optimized D1Nonequidistant opsets
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
40 x2 = 1.7473267488572e+00;
6b68f939d023 Add utility functions for constructing the grids used by the boundary optimized D1Nonequidistant opsets
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
41 x3 = 3.0000000000000e+00;
6b68f939d023 Add utility functions for constructing the grids used by the boundary optimized D1Nonequidistant opsets
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
42 xb = [x0 x1 x2 x3]';
6b68f939d023 Add utility functions for constructing the grids used by the boundary optimized D1Nonequidistant opsets
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
43 case 12
6b68f939d023 Add utility functions for constructing the grids used by the boundary optimized D1Nonequidistant opsets
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
44 x0 = 0.0000000000000e+00;
6b68f939d023 Add utility functions for constructing the grids used by the boundary optimized D1Nonequidistant opsets
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
45 x1 = 4.6552112904489e-01;
6b68f939d023 Add utility functions for constructing the grids used by the boundary optimized D1Nonequidistant opsets
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
46 x2 = 1.4647984306493e+00;
6b68f939d023 Add utility functions for constructing the grids used by the boundary optimized D1Nonequidistant opsets
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
47 x3 = 2.7620429464763e+00;
6b68f939d023 Add utility functions for constructing the grids used by the boundary optimized D1Nonequidistant opsets
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
48 x4 = 4.0000000000000e+00;
6b68f939d023 Add utility functions for constructing the grids used by the boundary optimized D1Nonequidistant opsets
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
49 xb = [x0 x1 x2 x3 x4]';
6b68f939d023 Add utility functions for constructing the grids used by the boundary optimized D1Nonequidistant opsets
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
50 otherwise
6b68f939d023 Add utility functions for constructing the grids used by the boundary optimized D1Nonequidistant opsets
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
51 error('Invalid operator order %d.',order);
6b68f939d023 Add utility functions for constructing the grids used by the boundary optimized D1Nonequidistant opsets
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
52 end
6b68f939d023 Add utility functions for constructing the grids used by the boundary optimized D1Nonequidistant opsets
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
53 end