diff --git a/Project.toml b/Project.toml
index 1333294..057d204 100644
--- a/Project.toml
+++ b/Project.toml
@@ -1,7 +1,7 @@
 name = "DoseCalculations"
 uuid = "b0d14063-c48a-4081-83d5-5a5ab48cb495"
 authors = ["lmejn <33491793+lmejn@users.noreply.github.com>"]
-version = "0.7.0"
+version = "0.7.1"
 
 [deps]
 CoordinateTransformations = "150eb455-5306-5404-9cee-2592286d6298"
diff --git a/src/Structures.jl b/src/Structures.jl
index 52d98b5..6f6c3da 100644
--- a/src/Structures.jl
+++ b/src/Structures.jl
@@ -11,25 +11,13 @@ mesh into a Meshes.SimpleMesh.
 From https://github.com/JuliaIO/MeshIO.jl/issues/67#issuecomment-913353708
 """
 function load_structure_from_ply(filename)
-    mesh = FileIO.load(filename)
-
-    vertexToIdx = Dict()
-    for i = 1:length(mesh.position)
-        vertex = mesh.position[i]
-        if !haskey(vertexToIdx, vertex) # duplicated vertices
-            vertexToIdx[vertex] = i
-        end
-    end
-    faces = []
-    for triangle in mesh
-        i1 = vertexToIdx[triangle[1]]
-        i2 = vertexToIdx[triangle[2]]
-        i3 = vertexToIdx[triangle[3]]
-        push!(faces, (i1, i2, i3))
-    end
-    topology = FullTopology(connect.(faces))
-
-    SimpleMesh([convert(Point{3, Float64}, Point([x[1], x[2], x[3]])) for x in mesh.position], topology)
+	mesh = FileIO.load(filename)
+	points = [Float64.(Tuple(p)) for p in Set(mesh.position)]
+	indices = Dict(p => i for (i, p) in enumerate(points))
+	connectivities = map(mesh) do el
+		Meshes.connect(Tuple(indices[Tuple(p)] for p in el))
+	end
+	Meshes.SimpleMesh(points, connectivities)
 end
 
 #--- Mesh Transformations -----------------------------------------------------