dimension_reduction failed for method='se'

[xyang2uchicago]

Hi, I am running dimension_reduction function for my 200k+ cells.
It works well when I call

st.dimension_reduction(adata, method='pca', n_neighbors=n_neighbors, feature='top_pcs', n_components=5, n_jobs=32)

However, it takes a while and was killed when I didn't specify the method, thus using default method='se':

st.dimension_reduction(adata, n_neighbors=n_neighbors, feature='top_pcs', n_components=5, n_jobs=32)

There is error I don't know how to debug:
Can't expand MemType 3: jcol 160965
Traceback (most recent call last):
File "", line 1, in
File "/project/xyang2/software-packages/env/stream/lib/python3.7/site-packages/stream/core.py", line 1179, in dimension_reduction
trans = reducer.fit(input_data)
File "/project/xyang2/software-packages/env/stream/lib/python3.7/site-packages/sklearn/manifold/_spectral_embedding.py", line 651, in fit
random_state=random_state,
File "/project/xyang2/software-packages/env/stream/lib/python3.7/site-packages/sklearn/manifold/_spectral_embedding.py", line 300, in spectral_embedding
laplacian, k=n_components, sigma=1.0, which="LM", tol=eigen_tol, v0=v0
File "/project/xyang2/software-packages/env/stream/lib/python3.7/site-packages/scipy/sparse/linalg/eigen/arpack/arpack.py", line 1645, in eigsh
hermitian=True, tol=tol)
File "/project/xyang2/software-packages/env/stream/lib/python3.7/site-packages/scipy/sparse/linalg/eigen/arpack/arpack.py", line 1073, in get_OPinv_matvec
return SpLuInv(A).matvec
File "/project/xyang2/software-packages/env/stream/lib/python3.7/site-packages/scipy/sparse/linalg/eigen/arpack/arpack.py", line 913, in init
self.M_lu = splu(M)
File "/project/xyang2/software-packages/env/stream/lib/python3.7/site-packages/scipy/sparse/linalg/dsolve/linsolve.py", line 339, in splu
ilu=False, options=_options)
SystemError: gstrf was called with invalid arguments

Anyone can help? Thank you.