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 this system on my PC and I need to repeat this calculations in loop so there's a need to speed up the calculations.
Is there a way to somehow parallelize the calculations using Python, e.g. with multiprocessing? Or maybe there is another way to shorten calculations time or split problem into many smaller ones? I've read about PyTrilinos but can't find any examples..