Hi, there seems to be a problem with inverses of matrices with elements in RR. It only occurs very sporadically for me, but here is an example:
a = RR(-4967757600021511 / 2**106) b = RR(-7769080564883485 / 2**52) c = RR( 5221315298319565 / 2**53) m = matrix([[a, b], [c, -a]]) print(m) print() print(~m) On my machines it produces the output [-6.12323399573677e-17 -1.72508242466029] [ 0.579682446302195 6.12323399573677e-17] [ 4.00000000000000 1.72508242466029] [ -0.579682446302195 -6.12323399573676e-17] Clearly the element 4 is wrong (the correct inverse is -m). Is this a known bug? Some system information: SageMath version 9.7, using Python 3.10.5 OS: Ubuntu 20.04.5 LTS CPU: Intel(R) Core(TM) i7-7700 CPU @ 3.60GHz Best regards, Håkan Granath -- You received this message because you are subscribed to the Google Groups "sage-devel" group. To unsubscribe from this group and stop receiving emails from it, send an email to sage-devel+unsubscr...@googlegroups.com. To view this discussion on the web visit https://groups.google.com/d/msgid/sage-devel/49610a42-ee51-4830-8785-79a86c36c621n%40googlegroups.com.
