diff --git a/Project.toml b/Project.toml
index 486a805..de4dae9 100644
--- 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
index 5c0c28d..8fe9d49 100644
--- 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))