Changeset 58 for include/fflas-ffpack/ffpack_frobenius.inl

Timestamp:

03/11/08 17:50:02 (10 months ago)

Author:

pernet

Message:

Apply fixes related to sage ticket #2453

Files:

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

include/fflas-ffpack/ffpack_frobenius.inl

@@ -46 +46 @@
                 for (size_t j = 0; j < N; ++j)
                         g.random( *(K + i*ldk +j) );
-
+        for (size_t i = 0; i < noc; ++i)
+                g.nonzerorandom (*(K + i*ldk +i));
+        
         // Computing the bloc Krylov matrix [U AU .. A^(c-1) U]^T
         for (size_t i = 1; i<c; ++i)
@@ -80 +82 @@
                 while ( (d<c) && (row_idx<R) && (Qk[row_idx] == i)) {i++; row_idx++; d++;}
                 if (d > dold){
-                        std::cerr << "FAIL in preconditionning phase:"
-                             << " degree sequence is not monotonically not increasing"
-                             << std::endl;
+                        // std::cerr << "FAIL in preconditionning phase:"
+                        //           << " degree sequence is not monotonically not increasing"
+                        //           << std::endl;
                         throw CharpolyFailed();
                 }
@@ -127 +129 @@
                         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;
+                                        //std::cerr<<"FAIL C != 0 in preconditionning"<<std::endl;
                                         throw CharpolyFailed();
                                 }
@@ -224 +226 @@
                 size_t R = SpecRankProfile (F, Ma, Ncurr, Arp, ldarp, deg-1, rp);
                 if (R < Ncurr){
-                        std::cerr<<"FAIL R<Ncurr"<<std::endl;
+                        //std::cerr<<"FAIL R<Ncurr"<<std::endl;
                         throw CharpolyFailed();
                 }
@@ -240 +242 @@
                         if ((block_idx)&&(it_idx > dK[block_idx-1])){
                                 throw CharpolyFailed();
-                                std::cerr<<"FAIL d non decroissant"<<std::endl;
+                                //std::cerr<<"FAIL d non decroissant"<<std::endl;
                                 //exit(-1);
                         }
@@ -312 +314 @@
                 if (LUdivine (F, FflasNonUnit, FflasNoTrans, Mk, Mk , K2 + (Ncurr-Mk)*ldk, ldk, P, Q, FfpackLQUP) < Mk){
                         // should never happen (not a LAS VEGAS check)
-                        std::cerr<<"FAIL R2 < MK"<<std::endl;
+                        //std::cerr<<"FAIL R2 < MK"<<std::endl;
                         //                      exit(-1);
                 }
@@ -353 +355 @@
                         for (size_t j=0; j<nb_full_blocks+1; ++j){
                                 if (!F.isZero( *(K+i*ldk+j) )){
-                                        std::cerr<<"FAIL C != 0"<<std::endl;
+                                        //std::cerr<<"FAIL C != 0"<<std::endl;
                                         throw CharpolyFailed();
                                 }