Merge pull request #98 from fastalgorithms/dev/chunkgraph-routs

Dev/chunkgraph routs
fastalgorithms · Oct 17, 2024 · 65ed58f · 65ed58f
2 parents fec212b + 7511883
commit 65ed58f
Show file tree

Hide file tree

Showing 25 changed files with 792 additions and 271 deletions.
diff --git a/chunkie/@chunker/onesmat.m b/chunkie/@chunker/onesmat.m
@@ -14,12 +14,12 @@
 %
 % Output:
 %   mat - the matrix discretization of the operator W above. mat is 
-%       (2*chnkr.npt) x (2*chnkr.npt)
+%       (chnkr.npt) x (chnkr.npt)
 %
 % Examples:
-%   mat = normonesmat(chnkr)
+%   mat = onesmat(chnkr)
 %
-% see also ONESMAT
+% see also NORMONESMAT
 
 % author: Travis Askham ([email protected])
 

diff --git a/chunkie/@chunkgraph/build_v2emat.m b/chunkie/@chunkgraph/build_v2emat.m
@@ -21,20 +21,22 @@
 for i = 1:nedges
     j1 = obj.edgesendverts(1,i);
     j2 = obj.edgesendverts(2,i);
-    if j1 == j2
-        nf = nf + 1;
-        ii(nf) = i;
-        jj(nf) = j1;
-        vv(nf) = 2;
-    else
-        nf = nf + 1;
-        ii(nf) = i;
-        jj(nf) = j1;
-        vv(nf) = -1;
-        nf = nf + 1;
-        ii(nf) = i;
-        jj(nf) = j2;
-        vv(nf) = 1;
+    if ~isnan(j1) && ~isnan(j2)
+        if j1 == j2
+            nf = nf + 1;
+            ii(nf) = i;
+            jj(nf) = j1;
+            vv(nf) = 2;
+        else
+            nf = nf + 1;
+            ii(nf) = i;
+            jj(nf) = j1;
+            vv(nf) = -1;
+            nf = nf + 1;
+            ii(nf) = i;
+            jj(nf) = j2;
+            vv(nf) = 1;
+        end
     end
 end