Changeset 2982 for trunk/linbox/linbox/algorithms/rational-cra-full-multip.h

Timestamp:

07/02/08 04:21:49 (5 months ago)

Author:

dumas

Message:

Ported to the improved full-multip-cra
(compute the inverse of the product only once)

Files:

• trunk/linbox/linbox/algorithms/rational-cra-full-multip.h (modified) (1 diff)

trunk/linbox/linbox/algorithms/rational-cra-full-multip.h

@@ -75 +75 @@
                     std::vector<Integer>::iterator t0_it = num.begin();
                     std::vector<Integer>::const_iterator t_it = _tab_it->begin();
+                    Integer invprod;
+                    precomputeInvProd(invprod, Product(), _mod_it->operator()() );
                     for( ; t0_it != num.end(); ++t0_it, ++t_it)
-                        this->normalizesmallbigreconstruct(*t0_it, Product(), *t_it, _mod_it->operator()() );
+                        this->smallbigreconstruct(*t0_it, *t_it, invprod );
                     Product.mulin(*_mod_it);
+
+                    // Moding out and normalization
+                    for(t0_it = num.begin();t0_it != num.end(); ++t0_it) {
+                        *t0_it %= Product();
+                        Integer tmp(*t0_it);
+                        normalize(*t0_it, tmp, Product());
+                    }
                 }
             }