diff --git a/src/ics_bcs.jl b/src/ics_bcs.jl
index 6a95397..9a422d0 100644
--- a/src/ics_bcs.jl
+++ b/src/ics_bcs.jl
@@ -249,10 +249,10 @@ end
 
 function apply_sm_schwarz_coupling_neumann(model::SolidMechanics, bc::CouplingSchwarzBoundaryCondition)
     schwarz_tractions = get_dst_traction(bc)
-    local_to_global_map = get_side_set_local_to_global_map(model.mesh, bc.side_set_id)
-    num_local_nodes = length(local_to_global_map)
+    global_from_local_map = get_side_set_global_from_local_map(model.mesh, bc.side_set_id)
+    num_local_nodes = length(global_from_local_map)
     for local_node ∈ 1:num_local_nodes
-        global_node = local_to_global_map[local_node]
+        global_node = global_from_local_map[local_node]
         node_tractions = schwarz_tractions[:, local_node]
         @debug "Applying Schwarz NBC as $node_tractions"
         model.boundary_force[3*global_node-2:3*global_node] += node_tractions
@@ -397,10 +397,10 @@ end
 function apply_sm_schwarz_contact_neumann(model::SolidMechanics, bc::SMContactSchwarzBC)
     schwarz_tractions = get_dst_traction(bc)
     normals = compute_normal(model.mesh, bc.side_set_id, model)
-    local_to_global_map = get_side_set_local_to_global_map(model.mesh, bc.side_set_id)
-    num_local_nodes = length(local_to_global_map)
+    global_from_local_map = get_side_set_global_from_local_map(model.mesh, bc.side_set_id)
+    num_local_nodes = length(global_from_local_map)
     for local_node ∈ 1:num_local_nodes
-        global_node = local_to_global_map[local_node]
+        global_node = global_from_local_map[local_node]
         node_tractions = schwarz_tractions[:, local_node]
         normal = normals[:, local_node]
         model.boundary_force[3*global_node-2:3*global_node] += transfer_normal_component(
@@ -417,11 +417,11 @@ function local_traction_from_global_force(
     side_set_id::Integer,
     global_force::Vector{Float64},
 )
-    local_to_global_map = get_side_set_local_to_global_map(mesh, side_set_id)
-    num_local_nodes = length(local_to_global_map)
+    global_from_local_map = get_side_set_global_from_local_map(mesh, side_set_id)
+    num_local_nodes = length(global_from_local_map)
     local_traction = zeros(3, num_local_nodes)
     for local_node ∈ 1:num_local_nodes
-        global_node = local_to_global_map[local_node]
+        global_node = global_from_local_map[local_node]
         local_traction[:, local_node] =
             global_force[3*global_node-2:3*global_node]
     end
diff --git a/src/interpolation.jl b/src/interpolation.jl
index b77d52e..bb6722c 100644
--- a/src/interpolation.jl
+++ b/src/interpolation.jl
@@ -684,15 +684,15 @@ function get_side_set_global_to_local_map(mesh::ExodusDatabase, side_set_id::Int
     return global_to_local_map, num_nodes_per_sides, side_set_node_indices
 end
 
-function get_side_set_local_to_global_map(mesh::ExodusDatabase, side_set_id::Integer)
+function get_side_set_global_from_local_map(mesh::ExodusDatabase, side_set_id::Integer)
     side_set_node_indices = Exodus.read_side_set_node_list(mesh, side_set_id)[2]
     unique_node_indices = unique(side_set_node_indices)
     num_nodes = length(unique_node_indices)
-    local_to_global_map = zeros(Int64, num_nodes)
+    global_from_local_map = zeros(Int64, num_nodes)
     for i ∈ 1:num_nodes
-        local_to_global_map[i] = Int64(unique_node_indices[i])
+        global_from_local_map[i] = Int64(unique_node_indices[i])
     end
-    return local_to_global_map
+    return global_from_local_map
 end
 
 function get_square_projection_matrix(
diff --git a/src/schwarz.jl b/src/schwarz.jl
index f084397..b3ebe44 100644
--- a/src/schwarz.jl
+++ b/src/schwarz.jl
@@ -362,8 +362,8 @@ function check_compression(
     compression = false
     nodal_reactions = get_dst_traction(bc)
     normals = compute_normal(mesh, bc.side_set_id, model)
-    local_to_global_map = get_side_set_local_to_global_map(mesh, bc.side_set_id)
-    num_local_nodes = length(local_to_global_map)
+    global_from_local_map = get_side_set_global_from_local_map(mesh, bc.side_set_id)
+    num_local_nodes = length(global_from_local_map)
     for local_node ∈ 1:num_local_nodes
         nodal_reaction = nodal_reactions[:, local_node]
         normal = normals[:, local_node]