- 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 ?