Parallelize Scipy iterative methods for linear equation systems(bicgstab) in Python

Question

I need to solve linear equations system Ax = b, where A is a sparse CSR matrix with size 500 000 x 500 000. I\'am using scipy.bicgstab and it takes almost 10min to solve thi