Disable subpixel smoothing by default in MaterialGrid class objects (…

…#2758) * disable subpixel smoothing by default in MaterialGrid * enable subpixel smoothing in test_adjoint_cyl.py to retain default behavior
NanoComp · Jan 13, 2024 · 4f67cb6 · 4f67cb6
1 parent 0dc44da
commit 4f67cb6
Show file tree

Hide file tree

Showing 3 changed files with 10 additions and 4 deletions.
diff --git a/doc/docs/Python_User_Interface.md b/doc/docs/Python_User_Interface.md
@@ -4691,7 +4691,7 @@ def __init__(self,
              medium2: meep.geom.Medium,
              weights: numpy.ndarray = None,
              grid_type: str = 'U_DEFAULT',
-             do_averaging: bool = True,
+             do_averaging: bool = False,
              beta: float = 0,
              eta: float = 0.5,
              damping: float = 0):

diff --git a/python/geom.py b/python/geom.py
@@ -597,7 +597,7 @@ def __init__(
         medium2: Medium,
         weights: np.ndarray = None,
         grid_type: str = "U_DEFAULT",
-        do_averaging: bool = True,
+        do_averaging: bool = False,
         beta: float = 0,
         eta: float = 0.5,
         damping: float = 0,

diff --git a/python/tests/test_adjoint_cyl.py b/python/tests/test_adjoint_cyl.py
@@ -72,7 +72,11 @@ def get_source(m):
 
 def forward_simulation(design_params, m, far_x):
     matgrid = mp.MaterialGrid(
-        mp.Vector3(Nr, 0, Nz), SiO2, Si, weights=design_params.reshape(Nr, 1, Nz)
+        mp.Vector3(Nr, 0, Nz),
+        SiO2,
+        Si,
+        weights=design_params.reshape(Nr, 1, Nz),
+        do_averaging=True,
     )
 
     geometry = [
@@ -113,7 +117,9 @@ def forward_simulation(design_params, m, far_x):
 
 def adjoint_solver(design_params, m, far_x):
 
-    design_variables = mp.MaterialGrid(mp.Vector3(Nr, 0, Nz), SiO2, Si)
+    design_variables = mp.MaterialGrid(
+        mp.Vector3(Nr, 0, Nz), SiO2, Si, do_averaging=True
+    )
     design_region = mpa.DesignRegion(
         design_variables,
         volume=mp.Volume(