Changeset 26 for include/fflas-ffpack/fflas_fgemm.inl

Timestamp:

07/06/07 18:58:19 (2 years ago)

Author:

pernet

Message:

New row echelon form and reduced row echelon form

Files:

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

include/fflas-ffpack/fflas_fgemm.inl

@@ -1142 +1142 @@
                 if (w <= 0) 
                         callClassicMatmul<double> () (F, ta, tb, m, n, k, 
-                                                 alpha, A, lda, B, ldb, beta, C, ldc, kmax,base);
+                                                      alpha, A, lda, B, ldb, beta, C, ldc, kmax,base);
                 else {
-                        if (k < kmax) { // switch on double
-                                // Temporary double matrices
+                        if (k < kmax) { // switch on delayed modulus
                                 DoubleDomain::Element _alpha, _beta;
                         
@@ -1169 +1168 @@
                                 WinoMain (DoubleDomain(), ta, tb, m, n, k, 
                                           _alpha, A, lda, B, ldb, _beta, C, ldc, kmax, w,base);
-                                // Conversion double = >  GFq
+                                // Modular reduction
                                 for (double * Ci = C; Ci != C+m*ldc; Ci+=ldc)
                                         for (size_t j = 0; j < n; ++j)
@@ -1180 +1179 @@
                                                         F.mulin (* (Ci + j), alpha);
                                 }
-                                // Temporary double matrices destruction
                         }
                         else{