diff --git a/Orb_AutFN.jl b/Orb_AutFN.jl
index 01a2bd7..533b012 100644
--- a/Orb_AutFN.jl
+++ b/Orb_AutFN.jl
@@ -53,6 +53,7 @@ function sparsify!{T}(U::SparseMatrixCSC{T}, eps=eps(T))
     return U
 end
 
+sparsify{T}(U::SparseMatrixCSC{T}, eps=eps(T)) = sparsify!(deepcopy(U), eps)
 small_to_zero!{T}(A::AbstractArray{T}, eps=eps(T)) = A[abs(A) .< eps] = zero(T)
 
 function init_model(Uπs)