Mercurial > repos > public > sbplib

changeset 1343:09a5783a3d37 feature/D2_boundary_opt

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Improve on default choice of number of grid points for each block
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Fri, 14 Oct 2022 14:42:32 +0200
parents d1dad4fbfe22
children b4e5e45bd239
files +multiblock/+domain/Annulus.m
diffstat 1 files changed, 8 insertions(+), 1 deletions(-) [+]
line wrap: on
line diff
--- a/+multiblock/+domain/Annulus.m	Fri Aug 26 14:19:29 2022 +0200
+++ b/+multiblock/+domain/Annulus.m	Fri Oct 14 14:42:32 2022 +0200
@@ -68,7 +68,14 @@
         end
 
         function ms = getGridSizes(obj, m)
-            ms = {[m m], [m m]};
+            mx = m;
+            % Use same grid spacing along inner 
+            % half circle as in radial direction
+            ds = pi*(obj.r_inner);
+            dr = (obj.r_outer-obj.r_inner);
+            my = ceil(ds/dr*(mx-1))+1;
+
+            ms = {[mx my], [mx my]};
         end
     end
 end