Mercurial > repos > public > sbplib

diff +multiblock/+domain/Annulus.m @ 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 4d472d020ccf
children
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