linbox
|
compare to cblas_dgemm instead.
compare to cblas_dgemm instead.
resize if it>nb !!
make long legends oblique.
ftrmm has an 'alpha' but mul/mulin in BMd don't... That could be useful for ...
benchmark ftrsm too here.
DenseMatrixBase<_element
>
and DenseSubmatrix<_element
>
or DenseMatrix<_field
>
... !!!
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,BlasTag::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,BlasTag::Left) ; avec P : BlasPermutation ou P : MatrixPermutation
BlasPermutation a un ordre p
et une taille r
distinctes !!!
use commentator.
use commentator.
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 !