Fix boundary_edge_matrix when there is no boundary

SciML · Apr 25, 2023 · b934738 · b934738 · DanielVandH · Apr 25, 2023
1 parent 71a3a69
commit b934738
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 2 deletions.
diff --git a/Project.toml b/Project.toml
@@ -1,7 +1,7 @@
 name = "FiniteVolumeMethod"
 uuid = "d4f04ab7-4f65-4d72-8a28-7087bc7f46f4"
 authors = ["Daniel VandenHeuvel <danj.vandenheuvel@gmail.com>"]
-version = "0.4.4"
+version = "0.4.5"
 
 [deps]
 ChunkSplitters = "ae650224-84b6-46f8-82ea-d812ca08434e"

diff --git a/src/geometry.jl b/src/geometry.jl
@@ -526,7 +526,7 @@ function boundary_edge_matrix(tri::Triangulation)
         bn_nodes = get_convex_hull_indices(tri)
         nedges = num_boundary_edges(bn_nodes)
         for edge_idx in 1:nedges
-            left_node = get_boundary_node(bn_nodes, edge_idx)
+            left_node = get_boundary_nodes(bn_nodes, edge_idx)
             right_node = get_right_boundary_node(tri, edge_idx, boundary_index)
             adjacent_node = get_adjacent(tri, left_node, right_node)
             push!(E, (left_node, right_node, adjacent_node, type))