Parallel rational matrix inversion in sage
Assume a large 2000x2000 rational matrix (matrix in QQ) for which I would like to compute the inverse. I noticed there are some papers on parallelizing matrix inversion - is any method for parallel matrix inversion already implemented in sage? Ideally, I would like to have the following:
q = 8380416
A = matrix(QQ, [[randint(0, q) for _ in range(1024)] for _ in range(1024)])
A.parallel_inverse()
Is something along those lines already implemented in sage?
I don't know, but let me ask why do you need this inverse? To solve what problem? If computing the inverse can be avoided then it is probably a good idea to do so.