如何在Numpy矩阵运算中保持gmpy2-mpfr的精度

Question

我在Numpy矩阵中使用多精度有理数（mpfr）对象

matrix([[ mpfr('-366998.93593422093364191959435957721088073331222596080623233278164906447646654043966366647797',300),
     mpfr('-366997.28868432286431885359868309613943011772698563764930700121744888472828510537502286003536',300),
     mpfr('-366997.28868432286431885359868309613943011772698563764930700121744888472828510537502286003536',300),
     mpfr('-366997.28868432286431885359868309613943011772698563764930700121744888472828510537502310955189',300),
     mpfr('-366997.33936304224917822062156336656390364691713762458391131405889211470102834400572590888586',300),
     mpfr('-366997.28868432286431885359868309613943011772698563764930700121744888472828510537502286003536',300)],
    [ mpfr('-40813927.104656436832435886099653290386078894027773129049451436960078610548203287954114434382',300),
     mpfr('-10418349883335.380900703935580692318458974868691020694148304775624032110383967472053357462067',300),
     mpfr('-40813927.104656436832435886099653290386078894027773129049451436960078610548203287954114434382',300),
     mpfr('-40813927.104656436832435886099653290386078894027773129049451436960078610548203287954114434382',300),
     mpfr('-40813927.104656436832435886099653290386078894027773129049451436960078610548203287954114434382',300),
     mpfr('-40813927.104656436832435886099653290386078894027773129049451436960078610548203287954114434382',300)]], dtype=object)

但是当计算矩阵的逆时，我会失去精度。在

^{pr2}$

那么如何保持mpfr的精度呢？ 如有任何建议，我们将不胜感激！在

Answer 1

出于性能原因，numpy使用LAPACK库，并且必须将矩阵元素转换为标准的double类型。您可能想尝试mpmath进行多精度矩阵求逆。在