Changeset 54cbca in git for factory


Ignore:
Timestamp:
Nov 22, 2007, 10:33:09 AM (16 years ago)
Author:
Hans Schönemann <hannes@…>
Branches:
(u'spielwiese', 'fe61d9c35bf7c61f2b6cbf1b56e25e2f08d536cc')
Children:
efdcc51591d42e61b55ae0e5cbb4586828e3a826
Parents:
9e23c68bff0ebd5ecd51b2b03137e3d8e7afd0ee
Message:
* avoid return 0 in gcd


git-svn-id: file:///usr/local/Singular/svn/trunk@10435 2c84dea3-7e68-4137-9b89-c4e89433aadc
File:
1 edited

Legend:

Unmodified
Added
Removed
  • factory/ffreval.cc

    r9e23c6 r54cbca  
    6060  }
    6161  return *this;
    62 }   
     62}
    6363
    6464/* forward declarations:*/
     
    117117        DEBOUTLN( cerr, "G = " << G );
    118118        if( ! fin_findeval( F, G, Fb, Gb, Db, b, delta, degF, degG ) )
    119            return 0;
     119        {
     120           Off(SW_USE_EZGCD_P);
     121           d *= gcd_poly_p(F,G);
     122           On(SW_USE_EZGCD_P);
     123           return d;
     124        }
    120125
    121126        DEBOUTLN( cerr, "found evaluation b = " << b );
     
    135140            bt = b;
    136141            if( ! fin_findeval( F, G, Fbt, Gbt, Dbt, bt, delta + 1, degF, degG ) )
    137                return 0;
     142            {
     143              Off(SW_USE_EZGCD_P);
     144              d *= gcd_poly_p(F,G);
     145              On(SW_USE_EZGCD_P);
     146              return d;
     147            }
    138148
    139149            int dd=degree( Dbt );
     
    179189                lcDD[2] = lcF;
    180190                B *= lcF;
    181                 B_is_F=true;
     191                B_is_F=true;
    182192            }
    183193            //else  if ( gcd( (DD[1] = Gb / Db), Db ) == 1 ) {
     
    196206                lcDD[2] = lcG;
    197207                B *= lcG;
    198                 B_is_F=false;
     208                B_is_F=false;
    199209              }
    200210              else
     
    279289        //     '(4*B^7*C^2-50*B^5*C^3+208*B^3*C^4-288*B*C^5)'
    280290        if( ! b.step() ) // Error: run out of points
    281            return 0;
    282 
     291        {
     292          Off(SW_USE_EZGCD_P);
     293          d = gcd_poly_p(F,G);
     294          On(SW_USE_EZGCD_P);
     295          return d;
     296        }
    283297        return fin_ezgcd( F, G, b, true );
    284298    }
Note: See TracChangeset for help on using the changeset viewer.