Changeset 4a7a45 in git for factory/fac_multihensel.cc


Ignore:
Timestamp:
Jul 29, 2020, 5:33:12 PM (4 years ago)
Author:
Hans Schoenemann <hannes@…>
Branches:
(u'spielwiese', '17f1d200f27c5bd38f5dfc6e8a0879242279d1d8')
Children:
700b89d0131c4b3d214f49af1eeb67501875ec89
Parents:
4772b1b4621f2ab68acc9d2b5a568604d4d5e790
Message:
fix: factorize in Z[x,..] w/o NTL
File:
1 edited

Legend:

Unmodified
Added
Removed
  • factory/fac_multihensel.cc

    r4772b1 r4a7a45  
    11/* 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"
    87#include "debug.h"
    98#include "timing.h"
     
    1413#include "fac_util.h"
    1514#include "fac_iterfor.h"
    16 #include "fac_multihensel.h"
    1715#include "cf_iter.h"
    1816
    1917#ifndef HAVE_NTL
    2018
    21 TIMING_DEFINE_PRINT(fac_solve)
    22 TIMING_DEFINE_PRINT(fac_modpk)
    23 TIMING_DEFINE_PRINT(fac_corrcoeff)
    24 TIMING_DEFINE_PRINT(fac_extgcd)
     19TIMING_DEFINE_PRINT(fac_solve);
     20TIMING_DEFINE_PRINT(fac_modpk);
     21TIMING_DEFINE_PRINT(fac_corrcoeff);
     22TIMING_DEFINE_PRINT(fac_extgcd);
    2523
    2624static 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*/ )
     25extgcdrest ( const CanonicalForm & a, const CanonicalForm & b, const CanonicalForm & s, const CanonicalForm & t, const CanonicalForm & c, CanonicalForm & S, CanonicalForm & T, const modpk & pk )
    2826{
    2927    CanonicalForm sigma = s * c, tau = t * c;
     
    169167        A[i] = remainder( pk( a[i] * C0 ), P0[i], pk );
    170168    DEBOUTLN( cerr, "the first approximation of the correction coefficients is " << A );
    171 /*#ifdef DEBUGOUTPUT
     169#ifdef DEBUGOUTPUT
    172170    if ( check_dummy( A, P, Q ) - C != 0 )
    173171    {
     
    177175        DEBOUTLN( cerr, "              Q " << Q );
    178176    }
    179 #endif*/
     177#endif
    180178    for ( m = 0; m <= h && ( m == 0 || Dm != 0 ); m++ )
    181179    {
     
    219217        DEBOUTLN( cerr, "the correction coefficients at step " << m );
    220218        DEBOUTLN( cerr, "are now " << A );
    221 /*#ifdef DEBUGOUTPUT
     219#ifdef DEBUGOUTPUT
    222220    if ( check_dummy( A, P, Q ) - C != 0 ) {
    223221        DEBOUTLN( cerr, "there is an error detected, the correction coefficients do not" );
     
    226224        DEBOUTLN( cerr, "              Q " << Q );
    227225    }
    228 #endif*/
     226#endif
    229227    }
    230228    DEBDECLEVEL( cerr, "findCorrCoeffs" );
     
    326324
    327325bool
    328 Hensel ( const CanonicalForm & U, CFArray & G, const CFArray & lcG, const Evaluation & A, const modpk & bound, const Variable & /*x*/ )
     326Hensel ( const CanonicalForm & U, CFArray & G, const CFArray & lcG, const Evaluation & A, const modpk & bound, const Variable & x )
    329327{
    330328    DEBINCLEVEL( cerr, "Hensel" );
Note: See TracChangeset for help on using the changeset viewer.