import implicit

# Use the implicit cython library.
user_vecs, item_vecs = implicit.alternating_least_squares(data_sparse,
        factors=20, regularization = 0.1, iterations = 20)

# Convert to csr matrix
user_vecs = sparse.csr_matrix(user_vecs)
item_vecs = sparse.csr_matrix(item_vecs)