Changeset 64a501 in git for factory/cf_gcd.cc
- Timestamp:
- Jul 25, 2011, 1:59:06 PM (12 years ago)
- Branches:
- (u'jengelh-datetime', 'ceac47cbc86fe4a15902392bdbb9bd2ae0ea02c6')(u'spielwiese', '00e2e9c41af3fde1273eb3633f4c0c7c3db2579d')
- Children:
- 21b8f4cad82601dd1e18c8bf2923808fb6827958
- Parents:
- 22071069039f22a2772d1e4d329f934d243f95da
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
factory/cf_gcd.cc
r220710 r64a501 1020 1020 CanonicalForm d; 1021 1021 #if 1 1022 do{ d = gcd_poly( f, g ); } 1023 while ((!fdivides(d,f)) || (!fdivides(d,g))); 1022 d = gcd_poly( f, g ); 1024 1023 #else 1025 while(1) 1026 { 1027 d = gcd_poly( f, g ); 1028 if ((fdivides(d,f)) && (fdivides(d,g))) break; 1029 printf("g"); fflush(stdout); 1030 } 1024 d = gcd_poly( f, g ); 1025 printf("g"); fflush(stdout); 1031 1026 #endif 1032 1027 return abs( d ); … … 1041 1036 CanonicalForm F = f * l, G = g * l; 1042 1037 Off( SW_RATIONAL ); 1043 do { l = gcd_poly( F, G ); } 1044 while ((!fdivides(l,F)) || (!fdivides(l,G))); 1038 l = gcd_poly( F, G ); 1045 1039 On( SW_RATIONAL ); 1046 1040 return abs( l );
Note: See TracChangeset
for help on using the changeset viewer.