linbox
Bug List
Global BlasMatrix< _Field >::apply (Vector1 &y, const Vector2 &x) const
this supposes &x[0]++ == &x[1]
Global BlasMatrixDomain< Field >::areEqual (const Matrix1 &A, const Matrix2 &B)
use refs
Global BlasMatrixDomain< Field >::isIdentity (const Matrix1 &A)

use refs

use refs

Global BlasMatrixDomain< Field >::isIdentityGeneralized (const Matrix1 &A)

use refs

use refs

Global BlasMatrixDomain< Field >::isZero (const Matrix1 &A)
use refs
Global BlasSubmatrix< _Field >::apply (Vector1 &y, const Vector2 &x) const
this supposes &x[0]++ == &x[1]
Global BlasSubmatrix< _Field >::read (std::istream &file)
reading a submatrix should not be allowed !!
globalScope> Global CheckMulAdd (const Integer &alpha, const BlasMatrix< PID_integer > &A, const BlasMatrix< PID_integer > &B, const Integer &beta, const BlasMatrix< PID_integer > &C)
si p n'est pas premier, fgemm fait n'importe quoi (division par alpha)
Global EarlyMultipCRA< Domain_Type >::dot (Integer &z, const Integer &D, const Vect1< Integer > &v1, const Vect2 &v2)
why a dot product here ?
Global EarlyMultipCRA< Domain_Type >::dot (DomainElement &z, const Domain &D, const Vect1< DomainElement, Alloc< DomainElement > > &v1, const Vect2 &v2)
why a dot product here ?
File field/modular.h
move Element& init(const Element&) to FFPACK. use using more..
Global IMLTraits::IMLTraits (const Specifier &S)
not complete
Global LinBox::checkBlasPrime (const Prime p)
why is this hard coded ?
Global LinBox::Identity (const Field &F, typename Field::Element *Id, const size_t ldI, const size_t lig1, const size_t col1, const size_t lig2, const size_t col2)
long et size_t ne cohabitent pas bien.
Global LinBox::Protected::random_rankupdate (Field &F,const Randiter &R, BlasMatrix< Field > &A, int &rank, const RingCategories::IntegerTag &tag)

const !

do perms ?

do perms ?

Global LinBox::solve (Vector &x, typename Ring::Element &d, const BlasMatrix< Ring > &A, const Vector &b, const RingCategories::IntegerTag &tag, const Method::BlasElimination &m)
check we don't copy
Global LQUPMatrix< Field >::LQUPMatrix (BlasMatrix< Field > &A, BlasPermutation< size_t > &P, BlasPermutation< size_t > &Q)
in place ?
globalScope> Global main (int ac, char **av)
: this one completely fails :
globalScope> Global main (int argc, char **argv)

testRandomIterator fails

testRandomIterator fails

the arguments are meaningless

globalScope> Global main (int argc, char *argv[])
q or not q ?
Class MatrixDomain< GF2 >
this is half done and makes MatrixDomain on GF2 hardly usable.
Global MatrixDomain< GF2 >::mulColSpecialized (Vector1 &w, const Matrix &A, const Vector2 &v, VectorCategories::DenseZeroOneVectorTag, VectorCategories::DenseZeroOneVectorTag) const
what's happening here ?
Class MatrixRank< _Ring, _Field, _RandomPrime >
there is no generic rankIn method.
Global MatrixRank< _Ring, _Field, _RandomPrime >::rank (const SparseMatrix< Ring, Row > &A) const
we suppose we can map IRing to Field...
Global MatrixRank< _Ring, _Field, _RandomPrime >::rank (const BlasMatrix< IRing > &A) const
we suppose we can map IRing to Field...
Global Modular< double >::operator= (const Modular< double > &F)
use FFPACK operator
Global ModularBalanced< double >::isMinusOne (const Element &x) const
faux si modulus==2
File rational-reconstruction.h
there seems to be a confusion between Ring _r, [iI]nteger... Besides we have to check that Ring/Integer/NTL_ZZ are the 'same'.
Global subtests::testRingTrivia (const Field &F, const char *name)
BlockRing does not know about 0 and 1 !
File test-charpoly.C

occasionnnaly there is a "SIGFPE, Arithmetic exception." in CRA

testRandomCharpoly is not always tested !!

File test-commentator.C
this test does not even test if anything is outputted (and it passes even though)
File test-ffpack.C
this file does not compile or the test fails for some fields.
File test-ftrmm.C

ftrmm fails on non-double transposed versions.

is fgemm correct on null A and B ?

File test-matrix-stream.C
this tests reveals memory leaks
globalScope> Global test_fgemm (std::ostream &report, const Field &F)

if beta != 0 but A and B are 0 (or don't exist) then C != beta C

needs p prime.

globalScope> Global test_full_multip_matrix (std::ostream &report, size_t PrimeSize, size_t Size, std::pair< size_t, size_t > dims)
it is not possible to allocate some memory and use submatrices ?