Mercurial > repos > public > sbplib
diff +sbp/+util/minimalBoundaryOptimizedGrid.m @ 1287:38653d26225c feature/boundary_optimized_grids
Make accurate/minimalBoundaryOptimizedGrid take the domain limits as input
| author | Vidar Stiernström <vidar.stiernstrom@it.uu.se> |
|---|---|
| date | Wed, 01 Jul 2020 14:54:21 +0200 |
| parents | 6b68f939d023 |
| children | e059a43bb675 |
line wrap: on
line diff
--- a/+sbp/+util/minimalBoundaryOptimizedGrid.m Wed Jul 01 13:43:32 2020 +0200 +++ b/+sbp/+util/minimalBoundaryOptimizedGrid.m Wed Jul 01 14:54:21 2020 +0200 @@ -1,4 +1,6 @@ -function [x,h] = minimalBoundaryOptimizedGrid(L,N,order) +function [x,h] = minimalBoundaryOptimizedGrid(lim,N,order) + L = lim{2} - lim{1}; + assert(L>0,'Limits must be given in increasing order.') %%%% Non-equidistant grid points %%%%% xb = boundaryPoints(order); m = length(xb)-1; % Number of non-equidistant points @@ -11,6 +13,7 @@ %%%% Define grid %%%%%%%% x = h*[xb; linspace(xb(end)+1,L/h-xb(end)-1,N-2*(m+1))'; L/h-flip(xb) ]; + x = x + lim{1}; %%%%%%%%%%%%%%%%%%%%%%%%% end