On Wednesday, December 19, 2012 11:07:03 PM UTC, William wrote: > > > sage: D=Matrix(IntegerModRing(), > [[-1,1,0,1,1,0],[1,-3,1,0,0,0],[0,1,-2,0,0,0],[1,0,0,-3,0,0],[1,0,0,0,-4,1],[0,0,0,0,1,-5]]);D > > > [-1 1 0 1 1 0] > > [ 1 -3 1 0 0 0] > > [ 0 1 -2 0 0 0] > > [ 1 0 0 -3 0 0] > > [ 1 0 0 0 -4 1] > > [ 0 0 0 0 1 -5] > > sage: X = D.LLL_gram(); >
This uses Pari lllgramint(), which assumes that the matrix is positive definite. If the matrix is not positive definite, Pari may not return. sage: D.is_positive_definite() False -- You received this message because you are subscribed to the Google Groups "sage-support" group. To post to this group, send email to sage-support@googlegroups.com. To unsubscribe from this group, send email to sage-support+unsubscr...@googlegroups.com. Visit this group at http://groups.google.com/group/sage-support?hl=en.
