inverse matrices

[frederic han]

Avec des matrices assez grandes il y a un probleme dans le calcul de l'inverse il me semble:

Code : Tout sélectionner

A:=ranm(20,20):;
(A*A^(-1))

ne donne pas souvent l'indentité quand la taille de A augmente

[parisse]

C'etait un bug stupide, un des entiers de Bezout utilise pour les restes chinois tenait sur un int, alors qu'un entier long etait suppose...
diff vecteur.cc vecteur.cc~ 4335,4336d ... coerce();