Changeset 5

Timestamp:

03/13/07 12:01:38 (2 years ago)

Author:

pernet

Message:

- ffpack_frobenius: Add a missing Las-Vegas Check in preconditionning

Files:

• include/fflas-ffpack/ffpack_frobenius.inl (modified) (3 diffs)

include/fflas-ffpack/ffpack_frobenius.inl

@@ -83 +83 @@
                              << " degree sequence is not monotonically not increasing"
                              << std::endl;
-                        //exit (-1);
                         throw CharpolyFailed();
                 }
@@ -126 +125 @@
         for (size_t i=Mk-1; i>=nb_full_blocks+1;  --i)
                 if (dK[i] >= 1){ 
+                        for (size_t j = offset+1; j<R; ++j)
+                                if (!F.isZero(*(K2b + i*ldk + j))){
+                                        std::cerr<<"FAIL C != 0 in preconditionning"<<std::endl;
+                                        throw CharpolyFailed();
+                                }
                         Polynomial P (dK [i]+1);
                         F.assign(P[dK[i]], one);
+                        cerr<<"Recovery in preconditionning P =";
                         for (size_t j=0; j < dK [i]; ++j){
                                 F.neg (P [dK [i]-j-1], *(K2b + i*ldk + (offset-j)));
+                                cerr<<" "<<P[dK[i]-j-1];
                         }
+                        cerr<<endl;
                         frobeniusForm.push_front(P);
                         offset -= dK [i];
@@ -138 +145 @@
         Mk = Ma;
 
-//      for (size_t i= offset+1; i<oldNcurr; ++i)
-//              for (size_t j=0; j<nb_full_blocks+1; ++j){
-//                      //              cerr<<"K + "<<i<<"*ldk + "<<j<<" = "<<(*(K+i*ldk+j))<<endl;
-//                      if (!F.isZero( *(K2b+i*ldk+j) )){
-//                              std::cerr<<"FAIL C != 0 in preconditionning"<<std::endl;
-//                              exit(-1);
-//                      }
-//              }
-        
         if (R<N){
-                
 //              std::cerr<<"Preconditionning failed; missing rank = "<<N-R
 //                       <<" completing the Krylov matrix"