On Sep 5, 10:14 am, Martin Albrecht <[EMAIL PROTECTED]> wrote: > On Friday 05 September 2008, phil wrote: > > > R.<x,y> = QQ[] > > C = random_matrix(R,10,10) > > Cdet = C.determinant() > > Here's a workaround: > sage: %time d2 = R(C._singular_().det())
Thanks for the tip. After making that change, Sage no longer crashes with an out of memory error on 64 bit Debian. However, I may need to make additional changes. The computation is still going after 3 days. Memory usage is slowly increasing and is now at 6.4 GB. The entries of the matrix are composed of coefficients extracted from some matrix operations on 4 3x3 matrices so there are 38 variables in the polynomial ring. Specifically, if anyone is wondering, I am trying to compute left hand side of equation 7 in "Five point motion estimation made easy" by Li and Hartley (http://users.rsise.anu.edu.au/ ~hongdong/new5pt_cameraREady_ver_1.pdf). The solver given on Li's webpage using Maple within Matlab to compute the determinant at runtime after the coefficients are given as numbers in the problem. I want to pre-compute the determinant with the coefficients specified by constants. That way at run time all you need to do is evaluate the expression replacing the constants with the numerical values. --~--~---------~--~----~------------~-------~--~----~ To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to [EMAIL PROTECTED] For more options, visit this group at http://groups.google.com/group/sage-support URLs: http://www.sagemath.org -~----------~----~----~----~------~----~------~--~---
