fflas-ffpack/checkers/checker_charpoly.inl

96             typename Field::Element_ptr Ac = FFLAS::fflas_new(F,n,n);  local
97             FFLAS::finit(F,n*n,Ac,1);
99             // Ac <- A - lambda.I
100             FFLAS::fassign(F,n,n,A,lda,Ac,n);
102                 F.subin(*(Ac+i*n+i),lambda);
105             // P,Ac,Q <- PLUQ(Ac)
112                 CheckerImplem_PLUQ<Field> checker (F,n,n,Ac,n);
122                 FFPACK::PLUQ(F, FFLAS::FflasNonUnit, n, n, Ac, n, P, Q);
129                 checker.check(Ac,n,FFLAS::FflasNonUnit,R,P,Q);
139             F.assign(det, *Ac);
141                 F.mulin(det,*(Ac+i*n+i));
142             if (n%2 == 1) F.negin(det); // Ac is A -lambda I
152             FFLAS::fflas_delete(Ac,P,Q);