Parallel rational matrix inversion in sage

asked 2021-09-14 16:47:34 +0200

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)])

Is something along those lines already implemented in sage?

edit retag flag offensive close merge delete


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.

rburing gravatar imagerburing ( 2021-09-15 08:53:21 +0200 )edit