Removed python-size upper limit to number of immersed boundaries

Author: danieljvickers

toolchain/mfc/case_validator.py

@@ -18,7 +18,7 @@
 from typing import Any, Dict, List, Set
 
 from .common import MFCException
-from .params.definitions import CONSTRAINTS, NI
+from .params.definitions import CONSTRAINTS
 from .state import CFG
 
 # Physics documentation for check methods.
@@ -559,7 +559,7 @@ def check_ibm(self):
         ib_state_wrt = self.get("ib_state_wrt", "F") == "T"
 
         self.prohibit(ib and n <= 0, "Immersed Boundaries do not work in 1D (requires n > 0)")
-        self.prohibit(ib and (num_ibs <= 0 or num_ibs > NI), f"num_ibs must be between 1 and {NI}")
+        self.prohibit(ib and num_ibs <= 0, "num_ibs must be >= 1 when ib is enabled")
         self.prohibit(not ib and num_ibs > 0, "num_ibs is set, but ib is not enabled")
         self.prohibit(ib_state_wrt and not ib, "ib_state_wrt requires ib to be enabled")
 

toolchain/mfc/params/definitions.py

@@ -639,7 +639,7 @@ def get_value_label(param_name: str, value: int) -> str:
     # Counts (must be positive)
     "num_fluids": {"min": 1, "max": 10},
     "num_patches": {"min": 0, "max": 10},
-    "num_ibs": {"min": 0, "max": NI},
+    "num_ibs": {"min": 0},
     "num_source": {"min": 1},
     "num_probes": {"min": 1},
     "num_integrals": {"min": 1},
@@ -1157,9 +1157,13 @@ def _load():
     for j in range(1, 4):
         _ib_attrs[f"vel({j})"] = (A_REAL, _ib_tags)
         _ib_attrs[f"angular_vel({j})"] = (A_REAL, _ib_tags)
-    REGISTRY.register_family(IndexedFamily(
-        base_name="patch_ib", attrs=_ib_attrs, tags=_ib_tags, max_index=NI,
-    ))
+    REGISTRY.register_family(
+        IndexedFamily(
+            base_name="patch_ib",
+            attrs=_ib_attrs,
+            tags=_ib_tags,
+        )
+    )
 
     # acoustic sources (4 sources)
     for i in range(1, NA + 1):

toolchain/mfc/params/registry.py

@@ -98,9 +98,7 @@ def __init__(
         for fam in families.values():
             for attr_name, (ptype, tags) in fam.attrs.items():
                 key = f"{fam.base_name}(1)%{attr_name}"
-                self._examples[key] = ParamDef(
-                    name=key, param_type=ptype, tags=set(tags)
-                )
+                self._examples[key] = ParamDef(name=key, param_type=ptype, tags=set(tags))
 
     def _resolve_family(self, name: str) -> Optional[ParamDef]:
         """Try to resolve a name against indexed families. Returns ParamDef or None."""
@@ -196,9 +194,7 @@ def freeze(self) -> None:
         """
         if not self._frozen:
             self._frozen = True
-            self._all_params_view = _FamilyAwareMapping(
-                self._params, self._families
-            )
+            self._all_params_view = _FamilyAwareMapping(self._params, self._families)
 
     @property
     def is_frozen(self) -> bool:
@@ -393,9 +389,7 @@ def get_json_schema(self) -> Dict[str, Any]:
             for attr_name, (ptype, _tags) in fam.attrs.items():
                 # Escape the attr name but replace sub-indices with \(\d+\)
                 attr_pattern = re.sub(r"\(\d+\)", "__IDX__", attr_name)
-                attr_pattern = re.escape(attr_pattern).replace(
-                    "__IDX__", r"\(\d+\)"
-                )
+                attr_pattern = re.escape(attr_pattern).replace("__IDX__", r"\(\d+\)")
                 pattern = f"^{base_esc}\\(\\d+\\)%{attr_pattern}$"
                 if pattern not in pattern_props:
                     pattern_props[pattern] = ptype.json_schema

toolchain/mfc/viz/interactive.py

@@ -434,7 +434,6 @@ def _get_cached_3d_mesh(
         return _mesh3_cache.get(key)
 
 
-
 # (PyVista code removed — server-side rendering uses kaleido on Linux)