- Timestamp:
- Nov 22, 2007, 10:33:09 AM (16 years ago)
- Branches:
- (u'spielwiese', 'fe61d9c35bf7c61f2b6cbf1b56e25e2f08d536cc')
- Children:
- efdcc51591d42e61b55ae0e5cbb4586828e3a826
- Parents:
- 9e23c68bff0ebd5ecd51b2b03137e3d8e7afd0ee
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
factory/ffreval.cc
r9e23c6 r54cbca 60 60 } 61 61 return *this; 62 } 62 } 63 63 64 64 /* forward declarations:*/ … … 117 117 DEBOUTLN( cerr, "G = " << G ); 118 118 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 } 120 125 121 126 DEBOUTLN( cerr, "found evaluation b = " << b ); … … 135 140 bt = b; 136 141 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 } 138 148 139 149 int dd=degree( Dbt ); … … 179 189 lcDD[2] = lcF; 180 190 B *= lcF; 181 191 B_is_F=true; 182 192 } 183 193 //else if ( gcd( (DD[1] = Gb / Db), Db ) == 1 ) { … … 196 206 lcDD[2] = lcG; 197 207 B *= lcG; 198 208 B_is_F=false; 199 209 } 200 210 else … … 279 289 // '(4*B^7*C^2-50*B^5*C^3+208*B^3*C^4-288*B*C^5)' 280 290 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 } 283 297 return fin_ezgcd( F, G, b, true ); 284 298 }
Note: See TracChangeset
for help on using the changeset viewer.