linbox
|
!!!
RandomPermutation avec P de type [Matrix-Blas]Permutation
: L = [[L1,0],[A,L2]] ;U = [[U1,B],[0,U2]] ; LU = [[ rec(L1,U1), ftrmm(L1,B)],[ftrmm(A,U1),fgemm(A,B)+rec(L2,U2) ]] de même UL
create BMD.applyP(A,P,LinBoxTag::Left) ; avec P : BlasPermutation ou P : MatrixPermutation
BlasPermutation a un ordre p
et une taille r
distinctes !!!
RandomPermutation avec P de type [Matrix-Blas]Permutation
: L = [[L1,0],[A,L2]] ;U = [[U1,B],[0,U2]] ; LU = [[ rec(L1,U1), ftrmm(L1,B)],[ftrmm(A,U1),fgemm(A,B)+rec(L2,U2) ]] de même UL
create BMD.applyP(A,P,LinBoxTag::Left) ; avec P : BlasPermutation ou P : MatrixPermutation
BlasPermutation a un ordre p
et une taille r
distinctes !!!
long
here ? its multipliers? its output form variants?
primary form? Relative primary form of A? Primary form is similar to A and finest which is a direct sum of companion matrices Cfi such that, for all i, j, gcd(fi, fj) = 1 or fi or fj. Relative primary form is coarsest such decomposition.
this may soon be reversed, in fact.
NULL
permutation F.isInvertible()
InvertibleRandomIter
check ftrsm fails nicely with non invertible A !