increase check tolerance for reflectance

NanoComp · Apr 12, 2022 · de9c776 · de9c776
1 parent fac8638
commit de9c776
Showing 1 changed file with 24 additions and 19 deletions.
diff --git a/python/tests/test_mode_decomposition.py b/python/tests/test_mode_decomposition.py
@@ -178,7 +178,7 @@ def test_grating_3d(self):
         wvl = 0.5  # wavelength
         fcen = 1/wvl
 
-        dpml = 2.0  # PML thickness
+        dpml = 1.0  # PML thickness
         dsub = 3.0  # substrate thickness
         dair = 3.0  # air padding
         hcyl = 0.5  # cylinder height
@@ -201,27 +201,29 @@ def test_grating_3d(self):
                              center=mp.Vector3(0,0,-0.5*sz+dpml),
                              component=src_cmpt)]
 
+        symmetries = [mp.Mirror(direction=mp.X,phase=-1),
+                      mp.Mirror(direction=mp.Y,phase=+1)]
+
         sim = mp.Simulation(resolution=resolution,
                             cell_size=cell_size,
                             sources=sources,
                             default_material=SiO2,
                             boundary_layers=boundary_layers,
-                            k_point=k_point)
-
-        flux = sim.add_mode_monitor(fcen,
-                                    0,
-                                    1,
-                                    mp.ModeRegion(center=mp.Vector3(0,0,-0.5*sz+dpml+0.5*dsub),
-                                                  size=mp.Vector3(sx,sy,0)))
-
-        stop_cond = mp.stop_when_fields_decayed(20,
-                                                src_cmpt,
-                                                mp.Vector3(0,0,0.5*sz-dpml-0.5*dair),
-                                                1e-6)
+                            k_point=k_point,
+                            symmetries=symmetries)
+
+        refl_pt = mp.Vector3(0,0,-0.5*sz+dpml+0.5*dsub)
+        refl_flux = sim.add_mode_monitor(fcen,
+                                         0,
+                                         1,
+                                         mp.ModeRegion(center=refl_pt,
+                                                       size=mp.Vector3(sx,sy,0)))
+
+        stop_cond = mp.stop_when_energy_decayed(20,1e-6)
         sim.run(until_after_sources=stop_cond)
 
-        input_flux = mp.get_fluxes(flux)
-        input_flux_data = sim.get_flux_data(flux)
+        input_flux = mp.get_fluxes(refl_flux)
+        input_flux_data = sim.get_flux_data(refl_flux)
 
         sim.reset_meep()
 
@@ -238,12 +240,13 @@ def test_grating_3d(self):
                             sources=sources,
                             geometry=geometry,
                             boundary_layers=boundary_layers,
-                            k_point=k_point)
+                            k_point=k_point,
+                            symmetries=symmetries)
 
         refl_flux = sim.add_mode_monitor(fcen,
                                          0,
                                          1,
-                                         mp.ModeRegion(center=mp.Vector3(0,0,-0.5*sz+dpml+0.5*dsub),
+                                         mp.ModeRegion(center=refl_pt,
                                                        size=mp.Vector3(sx,sy,0)))
         sim.load_minus_flux_data(refl_flux,input_flux_data)
 
@@ -314,9 +317,11 @@ def test_grating_3d(self):
         print("tran:, {}, {}".format(Tsum,Tflux))
         print("sum:,  {}, {}".format(Rsum+Tsum,Rflux+Tflux))
 
-        self.assertAlmostEqual(Rsum,Rflux,places=2)
+        ## to obtain agreement for two decimal digits,
+        ## the resolution must be increased to 200
+        self.assertAlmostEqual(Rsum,Rflux,places=1)
         self.assertAlmostEqual(Tsum,Tflux,places=2)
-        self.assertAlmostEqual(Rsum+Tsum,1.00,places=2)
+        self.assertAlmostEqual(Rsum+Tsum,1.00,places=1)
 
 
 if __name__ == '__main__':