Mercurial > repos > public > sbplib

diff +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
line wrap: on
line diff
--- a/+multiblock/DefCurvilinear.m	Mon Jun 01 19:00:31 2020 -0700
+++ b/+multiblock/DefCurvilinear.m	Tue Jun 02 11:44:38 2020 -0700
@@ -48,6 +48,18 @@
             g = multiblock.Grid(grids, obj.connections, obj.boundaryGroups);
         end
 
+        function g = getLebedevGrid(obj, varargin)
+            ms = obj.getGridSizes(varargin{:});
+
+            grids = cell(1, obj.nBlocks);
+            for i = 1:obj.nBlocks
+                % grids{i} = grid.equidistantCurvilinear(obj.blockMaps{i}.S, ms{i});
+                grids{i} = grid.lebedev2dCurvilinear(obj.blockMaps{i}.S, ms{i});
+            end
+
+            g = multiblock.Grid(grids, obj.connections, obj.boundaryGroups);
+        end
+
         function h = show(obj, label, gridLines, varargin)
             default_arg('label', 'name')
             default_arg('gridLines', false);