Revert "domain decomposition"

.gitignore

@@ -51,4 +51,3 @@ pytestdebug.log
 docs/build
 docs/source/pythonapi/generated/
 
-*.csv

mcdc/__init__.py

@@ -23,7 +23,6 @@
     uq,
     print_card,
     reset_cards,
-    domain_decomposition,
 )
 from mcdc.main import run, prepare
 from mcdc.visualizer import visualize

mcdc/card.py

@@ -80,12 +80,6 @@ def reset(self):
             "ww": np.ones([1, 1, 1, 1]),
             "ww_width": 2.5,
             "ww_mesh": make_card_mesh(),
-            "domain_decomposition": False,
-            "dd_idx": 0,
-            "dd_mesh": make_card_mesh(),
-            "dd_exchange_rate": 0,
-            "dd_repro": False,
-            "dd_work_ratio": np.array([1]),
             "weight_roulette": False,
             "wr_threshold": 0.0,
             "wr_survive": 1.0,

mcdc/constant.py

@@ -14,7 +14,6 @@
 EVENT_TIME_BOUNDARY = 1 << 7
 EVENT_LATTICE = 1 << 8
 EVENT_SURFACE_MOVE = 1 << 9
-EVENT_DOMAIN = 1 << 10
 
 # Gyration raius type
 GYRATION_RADIUS_ALL = 0
@@ -39,12 +38,6 @@
 PREC = 1.0 + 1e-5  # Precision factor to determine if a distance is smaller
 BANKMAX = 100  # Default maximum active bank
 
-# Domain Decomp mesh crossing flags
-MESH_X = 0
-MESH_Y = 1
-MESH_Z = 2
-MESH_T = 3
-
 # RNG LCG parameters
 RNG_G = nb.uint64(2806196910506780709)
 RNG_C = nb.uint64(1)

mcdc/input_.py

@@ -1173,50 +1173,6 @@ def weight_window(x=None, y=None, z=None, t=None, window=None, width=None):
     return card
 
 
-def domain_decomposition(
-    x=None,
-    y=None,
-    z=None,
-    t=None,
-    exchange_rate=100,
-    bank_size=1e5,
-    work_ratio=None,
-    repro=True,
-):
-    card = mcdc.input_deck.technique
-    card["domain_decomposition"] = True
-    card["domain_bank_size"] = int(1e5)
-    card["dd_exchange_rate"] = int(exchange_rate)
-    card["dd_repro"] = repro
-    dom_num = 1
-    # Set mesh
-    if x is not None:
-        card["dd_mesh"]["x"] = x
-        dom_num *= len(x)
-    if y is not None:
-        card["dd_mesh"]["y"] = y
-        dom_num *= len(y)
-    if z is not None:
-        card["dd_mesh"]["z"] = z
-        dom_num += len(z)
-    if t is not None:
-        card["dd_mesh"]["t"] = t
-        dom_num += len(t)
-    # Set work ratio
-    if work_ratio is None:
-        card["dd_work_ratio"] = None
-    elif work_ratio is not None:
-        card["dd_work_ratio"] = work_ratio
-    card["dd_idx"] = 0
-    card["dd_xp_neigh"] = []
-    card["dd_xn_neigh"] = []
-    card["dd_yp_neigh"] = []
-    card["dd_yn_neigh"] = []
-    card["dd_zp_neigh"] = []
-    card["dd_zn_neigh"] = []
-    return card
-
-
 def iQMC(
     g=None,
     t=None,