I am currently working for RDF and CDF mainly. I was planning to overload LU function of *sage/matrix/matrix2.pyx* by defining new LU in *sage/matrix/matrix_double_sparce.pyx* .
I couldn't understand what is multiprecision version but I guess you are asking for RR or other inexact field. Please correct me if I am wrong. On Wednesday, February 28, 2024 at 9:51:32 PM UTC+5:30 Nils Bruin wrote: > As Dima points out, SPLU decomposition is probably more useful in > practice, so it's great to expose it. It should probably be exposed through > a different method than PLU decomposition, though. > > For the return format: permutations are much more compact and efficient to > describe permutations than permutation matrices. Also because of the > parent: a permutation matrix is a {0,1}-valued matrix so it can be well > represented over ZZ, rather than floats. So primarily the routine should > expose that scipy exposes. Whether other bells and whistles are required is > another matter. > > Is this only for system floats/complex numbers? or will you also be > implementing multiprecision versions? > -- You received this message because you are subscribed to the Google Groups "sage-devel" group. To unsubscribe from this group and stop receiving emails from it, send an email to sage-devel+unsubscr...@googlegroups.com. To view this discussion on the web visit https://groups.google.com/d/msgid/sage-devel/32c0253b-6456-4522-8638-8dd958687bf6n%40googlegroups.com.