use ludivine for ludecomp

[jakobkroeker]

No description provided.

[jakobkroeker]

I can't see at a first sight, how to use ludivine for PLU decomposition in general.

How it should work for this example?

loadPackage "FastLinearAlgebra"
debug Core
R =  ZZp(101,Strategy=>"Ffpack")
m = mutableMatrix matrix{{0_R,6_R},{0_R,0_R},{0_R,4_R}}
 LUdecomposition m
(P,Q) = rawLQUPFactorization(raw m)

see also https://groups.google.com/forum/#!topic/ffpack-devel/THx0TmVCSa8

[mikestillman]

Here is what I think it should be:

i9 : LUdecomposition m

o9 = ({0, 1, 2}, | 1   . |, | . 6 |)
                 | .   1 |  | . . |
                 | -33 . |