typos. I swear my tests were working great at one point ¯\_(ツ)_/¯

khanlab · Jan 28, 2025 · b51000b · b51000b
1 parent e2c860b
commit b51000b
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 5 deletions.
diff --git a/hippunfold/workflow/scripts/gen_isosurface.py b/hippunfold/workflow/scripts/gen_isosurface.py
@@ -213,7 +213,7 @@ def compute_geodesic_distances(vertices, faces, source_indices):
     ]  # really hope there's no x-y switching fuckery here!
 
 # keep vertices that are in a nice coordinate range
-epsilon - snakemake.params.coords_epsilon
+epsilon = snakemake.params.coords_epsilon
 good_v = np.where(np.logical_and(coord_at_V < (1 - epsilon), coord_at_V > epsilon))[0]
 
 # morphological open

diff --git a/hippunfold/workflow/scripts/laplace_beltrami.py b/hippunfold/workflow/scripts/laplace_beltrami.py
@@ -147,14 +147,12 @@ def solve_laplace_beltrami_open_mesh(vertices, faces, boundary_conditions=None):
 
 APinds = np.array(np.where(boundary_values < 12)[0]).astype(int)
 boundary_conditions = dict(zip(boundary_vertices[APinds], boundary_values[APinds] - 10))
-APcoords = np.ones((len(vertices_orig))) * np.nan
-APcoords[i_concomp] = solve_laplace_beltrami_open_mesh(
+APcoords = solve_laplace_beltrami_open_mesh(
     vertices, faces, boundary_conditions
 )
 PDinds = np.array(np.where(boundary_values > 12)[0]).astype(int)
 boundary_conditions = dict(zip(boundary_vertices[PDinds], boundary_values[PDinds] - 20))
-PDcoords = np.ones((len(vertices_orig))) * np.nan
-PDcoords[i_concomp] = solve_laplace_beltrami_open_mesh(
+PDcoords = solve_laplace_beltrami_open_mesh(
     vertices, faces, boundary_conditions
 )