On 24 juin, 07:50, Jason Grout <[email protected]> wrote: > And 5. Compare (maybe wrap??) to the nice LU decomposition routine that > already exists in sympy! I just found this. It looks like very nice > work. I should probably look at sympy more often when there are missing > bits in Sage!
For symbolic matrices, comparisons seem to give a clear advantage to the sympy code. Think I'll go for that. (as for brokenness of your patch for symbolic matrices, the problem isn't as much inexactness as undecidability of comparisons. Because of that, partial pivoting sticks with a fake zero pivot, which obviously fails). E. --~--~---------~--~----~------------~-------~--~----~ To post to this group, send email to [email protected] 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 -~----------~----~----~----~------~----~------~--~---
