diff --git a/geometry/3d_beam_nx/beam_3d.jl b/geometry/3d_beam_nx/beam_3d.jl
index efac829..53b92b4 100644
--- a/geometry/3d_beam_nx/beam_3d.jl
+++ b/geometry/3d_beam_nx/beam_3d.jl
@@ -136,16 +136,12 @@ function lsq_fit(elements, field)
     SparseArrays.droptol!(A, 1.0e-6)
     nz = get_nonzero_rows(A)
     dropzeros!(A)
-    #try
-        F = ldltfact(A)
+    F = ldltfact(A)
 
-        x = zeros(size(b)...)
-        x[nz, :] = F \ b[nz, :]
-        nodal_values = Dict(i => vec(x[i,:]) for i=1:size(x,1))
-    #catch
-        err("Problem with zeros or tolerances")
-        dump(nz)
-    #end
+    x = zeros(size(b)...)
+    x[nz, :] = F \ b[nz, :]
+
+    nodal_values = Dict(i => vec(x[i,:]) for i=1:size(x,1))
     return nodal_values
 end