Mercurial > repos > public > sbplib

comparison +multiblock/DefCurvilinear.m @ 1274:8aa0909125a4 feature/poroelastic

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add support for plotting and evaluating on curvilinear multiblock staggered grids
author Martin Almquist <malmquist@stanford.edu>
date Tue, 02 Jun 2020 11:44:38 -0700
parents e7a6744499fa
children 60c875c18de3
comparison
equal deleted inserted replaced
1273:cf542444f022 1274:8aa0909125a4
41 ms = obj.getGridSizes(varargin{:}); 41 ms = obj.getGridSizes(varargin{:});
42 42
43 grids = cell(1, obj.nBlocks); 43 grids = cell(1, obj.nBlocks);
44 for i = 1:obj.nBlocks 44 for i = 1:obj.nBlocks
45 grids{i} = grid.equidistantCurvilinear(obj.blockMaps{i}.S, ms{i}); 45 grids{i} = grid.equidistantCurvilinear(obj.blockMaps{i}.S, ms{i});
46 end
47
48 g = multiblock.Grid(grids, obj.connections, obj.boundaryGroups);
49 end
50
51 function g = getLebedevGrid(obj, varargin)
52 ms = obj.getGridSizes(varargin{:});
53
54 grids = cell(1, obj.nBlocks);
55 for i = 1:obj.nBlocks
56 % grids{i} = grid.equidistantCurvilinear(obj.blockMaps{i}.S, ms{i});
57 grids{i} = grid.lebedev2dCurvilinear(obj.blockMaps{i}.S, ms{i});
46 end 58 end
47 59
48 g = multiblock.Grid(grids, obj.connections, obj.boundaryGroups); 60 g = multiblock.Grid(grids, obj.connections, obj.boundaryGroups);
49 end 61 end
50 62