Thank you!
On Saturday, November 5, 2022 at 12:15:26 PM UTC+1 dim...@gmail.com wrote:
> I've opened https://trac.sagemath.org/ticket/34724 to deal with this
>
> On Sat, Nov 5, 2022 at 11:08 AM Dima Pasechnik wrote:
> >
> > Well, applying a naive exact linear algebra routine to inexact data,
> >
I've opened https://trac.sagemath.org/ticket/34724 to deal with this
On Sat, Nov 5, 2022 at 11:08 AM Dima Pasechnik wrote:
>
> Well, applying a naive exact linear algebra routine to inexact data,
> and that's what Sage is doing here, is prone to errors.
> sage: A=m.augment(identity_matrix(RR,2))
Well, applying a naive exact linear algebra routine to inexact data,
and that's what Sage is doing here, is prone to errors.
sage: A=m.augment(identity_matrix(RR,2))
sage: A
[-6.12323399573677e-17 -1.72508242466029 1.00
0.000]
[0.579682446302195 6.1232339957367
something is *definitely* unhinged here : On 9.8.beta3 running on Debian
testing on core i7 + 16 GB RAM, after running :
a = RR(-4967757600021511 / 2**106)
b = RR(-7769080564883485 / 2**52)
c = RR( 5221315298319565 / 2**53)
m = matrix([[a, b], [c, -a]])
M = matrix([[var("p%d%d"%(u, v), latex_na
