diff -r 51dbb2cc25d9 -r 8d1c959432c8 src/lib.rs
--- a/src/lib.rs	Tue Aug 04 13:46:29 2026 -0500
+++ b/src/lib.rs	Tue Aug 04 14:12:18 2026 -0500
@@ -21,6 +21,7 @@
 //use pyo3::pyclass;
 
 pub mod dataterm;
+pub mod domain;
 pub mod fb;
 pub mod forward_model;
 pub mod forward_pdps;
@@ -242,9 +243,14 @@
     pub max_scaling_iter: Option<usize>,
 
     #[arg(long)]
-    /// Extra finite dimensional steps to take
+    /// Extra finite dimensional steps to take after merge and random insertion.
+    /// This all methods support this.
     pub extra_weight_optimisation_steps: Option<usize>,
 
+    #[arg(long)]
+    /// Whether to perform random insertions after merge. Not all methods support this.
+    pub random_insertions: Option<bool>,
+
     #[arg(long, value_names = &["ε", "θ", "p"])]
     /// Set the tolerance to ε_k = ε/(1+θk)^p
     pub tolerance: Option<Vec<F>>,