diff +multiblock/DefCurvilinear.m @ 1331:60c875c18de3 feature/D2_boundary_opt

Merge with feature/poroelastic for Elastic schemes
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Thu, 10 Mar 2022 16:54:26 +0100
parents 7df63b17e078 8aa0909125a4
children 0666629aa183
line wrap: on
line diff
--- a/+multiblock/DefCurvilinear.m	Thu Feb 17 18:55:11 2022 +0100
+++ b/+multiblock/DefCurvilinear.m	Thu Mar 10 16:54:26 2022 +0100
@@ -74,6 +74,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);