fix isBdNode not assigned a value bug (#12)

* fix `isBdNode` not being assigned a value in gradmatrix.m
lyc102 · Nov 10, 2023 · 884a92a · 884a92a
1 parent 29f31c8
commit 884a92a
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/solver/gradmatrix.m b/solver/gradmatrix.m
@@ -21,9 +21,9 @@
 i = repmat((1:NE)',2,1);
 j = double(edge(:));
 s = [-ones(NE,1),ones(NE,1)];
+isBdNode = false(N,1);
 if any(isBdEdge) % no grad on boundary edges
     bdEdge = edge(isBdEdge,:);
-    isBdNode = false(N,1);
     isBdNode(bdEdge(:)) = true;
     idx = ~(isBdEdge(i) | isBdNode(j)); 
 else