debug numba wrinkle

Author: ilhamv

mcdc/constant.py

@@ -74,14 +74,13 @@
 # Geometry events
 EVENT_SURFACE_CROSSING = 1 << 1
 EVENT_LATTICE_CROSSING = 1 << 2
-EVENT_DOMAIN_CROSSING = 1 << 3
-EVENT_LOST = 1 << 4
+EVENT_LOST = 1 << 3
 # Collision/reaction events
-EVENT_COLLISION = 1 << 5
+EVENT_COLLISION = 1 << 4
 # Miscellanies
-EVENT_TIME_CENSUS = 1 << 6
-EVENT_TIME_BOUNDARY = 1 << 7
-EVENT_IQMC_MESH = 1 << 8
+EVENT_TIME_CENSUS = 1 << 5
+EVENT_TIME_BOUNDARY = 1 << 6
+EVENT_IQMC_MESH = 1 << 7
 
 # Materials
 MATERIAL = 0

mcdc/main.py

@@ -453,7 +453,10 @@ def visualize(
         plt.ylabel(second_key + " [cm]")
         plt.title(reference_key + " = %.2f cm" % reference + ", time = %.2f s" % t)
         if save_as is not None:
-            plt.savefig(save_as + ".png")
+            if len(time) > 1:
+                plt.savefig(f"{save_as}_{t:03}.png")
+            else:
+                plt.savefig(save_as + ".png")
             plt.clf()
         else:
             plt.show()

mcdc/transport/kernel.py

@@ -504,7 +504,6 @@ def move_to_event(P_arr, mcdc, data):
 
     # Distance to domain
     speed = physics.particle_speed(P_arr, mcdc, data)
-    d_domain = INF
 
     # Distance to time boundary
     d_time_boundary = speed * (settings["time_boundary"] - P["t"])
@@ -525,14 +524,6 @@ def move_to_event(P_arr, mcdc, data):
 
     distance = d_boundary
 
-    # Check distance to domain
-    if d_domain < distance - COINCIDENCE_TOLERANCE:
-        distance = d_domain
-        P["event"] = EVENT_DOMAIN_CROSSING
-        P["surface_ID"] = -1
-    elif geometry.check_coincidence(d_domain, distance):
-        P["event"] += EVENT_DOMAIN_CROSSING
-
     # Check distance to collision
     if d_collision < distance - COINCIDENCE_TOLERANCE:
         distance = d_collision

mcdc/transport/simulation.py

@@ -483,12 +483,6 @@ def step_particle(P_arr, prog, data):
     # Surface and domain crossing
     if P["event"] & EVENT_SURFACE_CROSSING:
         geometry.surface_crossing(P_arr, prog, data)
-        if P["event"] & EVENT_DOMAIN_CROSSING:
-            if mcdc["surfaces"][P["surface_ID"]]["boundary_condition"] == BC_NONE:
-                kernel.domain_crossing(P_arr, prog)
-
-    elif P["event"] & EVENT_DOMAIN_CROSSING:
-        kernel.domain_crossing(P_arr, prog)
 
     # Census time crossing
     if P["event"] & EVENT_TIME_CENSUS: