Le 08/05/2015 10:23, Vincent Delecroix a écrit :
Bonjour Thierry,Could you post your matrices Prec, M, mfeEast? I am not able to reproduce the problem with change_ring... Vincent
Hello Vincent. Yes, I join a part of my code, which reproduces the problem... Yours t. -- You received this message because you are subscribed to the Google Groups "sage-support" group. To unsubscribe from this group and stop receiving emails from it, send an email to sage-support+unsubscr...@googlegroups.com. To post to this group, send email to sage-support@googlegroups.com. Visit this group at http://groups.google.com/group/sage-support. For more options, visit https://groups.google.com/d/optout.
lprec=[(22502471856513972497017661440000/940135272187668275673825470054308951, 0, -7044387942611760194091358546315/940135272187668275673825470054308951, 0), (0, 628091429911656478463873925/115769812361316527145683860922368, 0, -12849448754121499186194957/7235613272582282946605241307648), (-7044387942611760194091358546315/940135272187668275673825470054308951, 0, 4500494371302794499403532288000/940135272187668275673825470054308951, 0), (0, -12849448754121499186194957/7235613272582282946605241307648, 0, 525746174591230144444035/226112914768196342081413790864)] mlist=[(-40960, 1423743620867035193/22298143521175, -483387868093604552/7539677263923, 34549818091021825/797033922567), (-1423743620867035193/22298143521175, 98304, -311185622448809473/3269359628156, 3297567640129996/58560241971), (-483387868093604552/7539677263923, 311185622448809473/3269359628156, -3752999689502719/45812984491, 8587337044602103/354376193361), (-34549818091021825/797033922567, 3297567640129996/58560241971, -8587337044602103/354376193361, -4621642813466987640258501351/80595054643327775067341)] Prec=matrix(QQ,lprec) M=matrix(QQ,mlist) print "Prec:",Prec.parent() print "M:",M.parent() LC=Prec.cholesky() print "LC: ",LC.parent() print "LC: ",LC.list() LCt=LC.transpose() print "LCt: ",LCt.parent() # Z=LC*M*LCt print "LC*M*LCt: ",Z.parent() # print "Now, we play..." M1=Z.change_ring(RDF)
<<attachment: tdumont.vcf>>