python - 用scipy.optimize的时候怎么处理大矩阵？

我有一个用户对商品评分的数据集，大概两万多行，里面有两万个用户和四千样商品。因为要用pymc3，我建立了一个稀疏的user-item-matrix，希望能得到一个user profile matrix 和 item profile matrix。我用的model是probabilistic matrix factorization，用的推断方法是MCMC。首先要用Scipy的optimize function 来求出start。 我试了很多里面的function，例如bfgs,lbfgs之类，但是全部都：

 [1]    45535 killed     python pmf2.py

我想问有什么方法可以解决这个问题？？万分感谢！