reduced randomness for breaking geometrical divide-by-zero cases

TomGeorge1234 · TomGeorge1234 · commit b2d7ef245e6c · 2025-02-24T19:45:15.000-08:00
diff --git a/ratinabox/utils.py b/ratinabox/utils.py
@@ -62,10 +62,10 @@ def vector_intercepts(vector_list_a, vector_list_b, return_collisions=False):
     vector_list_a = vector_list_a.reshape(-1, 2, 2)
     vector_list_b = vector_list_b.reshape(-1, 2, 2)
     vector_list_a = vector_list_a + np.random.normal(
-        scale=1e-6, size=vector_list_a.shape
+        scale=1e-9, size=vector_list_a.shape
     )
     vector_list_b = vector_list_b + np.random.normal(
-        scale=1e-6, size=vector_list_b.shape
+        scale=1e-9, size=vector_list_b.shape
     )
 
     N_a = vector_list_a.shape[0]

Original file line number	Diff line number	Diff line change
`@@ -62,10 +62,10 @@ def vector_intercepts(vector_list_a, vector_list_b, return_collisions=False):`
`62`	`62`	`vector_list_a = vector_list_a.reshape(-1, 2, 2)`
`63`	`63`	`vector_list_b = vector_list_b.reshape(-1, 2, 2)`
`64`	`64`	`vector_list_a = vector_list_a + np.random.normal(`
`65`		`- scale=1e-6, size=vector_list_a.shape`
	`65`	`+ scale=1e-9, size=vector_list_a.shape`
`66`	`66`	`)`
`67`	`67`	`vector_list_b = vector_list_b + np.random.normal(`
`68`		`- scale=1e-6, size=vector_list_b.shape`
	`68`	`+ scale=1e-9, size=vector_list_b.shape`
`69`	`69`	`)`
`70`	`70`
`71`	`71`	`N_a = vector_list_a.shape[0]`