Changeset 4a7a45 in git for factory/fac_multihensel.cc
- Timestamp:
- Jul 29, 2020, 5:33:12 PM (4 years ago)
- Branches:
- (u'spielwiese', '17f1d200f27c5bd38f5dfc6e8a0879242279d1d8')
- Children:
- 700b89d0131c4b3d214f49af1eeb67501875ec89
- Parents:
- 4772b1b4621f2ab68acc9d2b5a568604d4d5e790
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
factory/fac_multihensel.cc
r4772b1 r4a7a45 1 1 /* emacs edit mode for this file is -*- C++ -*- */ 2 3 4 #include "config.h" 5 6 7 #include "cf_assert.h" 2 /* $Id: fac_multihensel.cc 12231 2009-11-02 10:12:22Z hannes $ */ 3 4 #include <config.h> 5 6 #include "assert.h" 8 7 #include "debug.h" 9 8 #include "timing.h" … … 14 13 #include "fac_util.h" 15 14 #include "fac_iterfor.h" 16 #include "fac_multihensel.h"17 15 #include "cf_iter.h" 18 16 19 17 #ifndef HAVE_NTL 20 18 21 TIMING_DEFINE_PRINT(fac_solve) 22 TIMING_DEFINE_PRINT(fac_modpk) 23 TIMING_DEFINE_PRINT(fac_corrcoeff) 24 TIMING_DEFINE_PRINT(fac_extgcd) 19 TIMING_DEFINE_PRINT(fac_solve); 20 TIMING_DEFINE_PRINT(fac_modpk); 21 TIMING_DEFINE_PRINT(fac_corrcoeff); 22 TIMING_DEFINE_PRINT(fac_extgcd); 25 23 26 24 static void 27 extgcdrest ( const CanonicalForm & a, const CanonicalForm & b, const CanonicalForm & s, const CanonicalForm & t, const CanonicalForm & c, CanonicalForm & S, CanonicalForm & T, const modpk & /*pk*/)25 extgcdrest ( const CanonicalForm & a, const CanonicalForm & b, const CanonicalForm & s, const CanonicalForm & t, const CanonicalForm & c, CanonicalForm & S, CanonicalForm & T, const modpk & pk ) 28 26 { 29 27 CanonicalForm sigma = s * c, tau = t * c; … … 169 167 A[i] = remainder( pk( a[i] * C0 ), P0[i], pk ); 170 168 DEBOUTLN( cerr, "the first approximation of the correction coefficients is " << A ); 171 /*#ifdef DEBUGOUTPUT169 #ifdef DEBUGOUTPUT 172 170 if ( check_dummy( A, P, Q ) - C != 0 ) 173 171 { … … 177 175 DEBOUTLN( cerr, " Q " << Q ); 178 176 } 179 #endif */177 #endif 180 178 for ( m = 0; m <= h && ( m == 0 || Dm != 0 ); m++ ) 181 179 { … … 219 217 DEBOUTLN( cerr, "the correction coefficients at step " << m ); 220 218 DEBOUTLN( cerr, "are now " << A ); 221 /*#ifdef DEBUGOUTPUT219 #ifdef DEBUGOUTPUT 222 220 if ( check_dummy( A, P, Q ) - C != 0 ) { 223 221 DEBOUTLN( cerr, "there is an error detected, the correction coefficients do not" ); … … 226 224 DEBOUTLN( cerr, " Q " << Q ); 227 225 } 228 #endif */226 #endif 229 227 } 230 228 DEBDECLEVEL( cerr, "findCorrCoeffs" ); … … 326 324 327 325 bool 328 Hensel ( const CanonicalForm & U, CFArray & G, const CFArray & lcG, const Evaluation & A, const modpk & bound, const Variable & /*x*/)326 Hensel ( const CanonicalForm & U, CFArray & G, const CFArray & lcG, const Evaluation & A, const modpk & bound, const Variable & x ) 329 327 { 330 328 DEBINCLEVEL( cerr, "Hensel" );
Note: See TracChangeset
for help on using the changeset viewer.