Changeset 102daa in git
- Timestamp:
- Apr 22, 2013, 4:32:37 PM (11 years ago)
- Branches:
- (u'spielwiese', 'fe61d9c35bf7c61f2b6cbf1b56e25e2f08d536cc')
- Children:
- e1a22119c89dfbc8533ea4d335f2fa9ee5ec9811
- Parents:
- 6fcd65b41235516be599e23554d96064636cb46c
- git-author:
- Martin Lee <martinlee84@web.de>2013-04-22 16:32:37+02:00
- git-committer:
- Martin Lee <martinlee84@web.de>2013-05-02 11:42:41+02:00
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
factory/cfNewtonPolygon.cc
r6fcd65b r102daa 23 23 #include "templates/ftmpl_functions.h" 24 24 #include "algext.h" 25 26 void out_cf(const char *s1,const CanonicalForm &f,const char *s2);27 25 28 26 static … … 923 921 CanonicalForm g= gcd (newtonPolyg[0][0], newtonPolyg[0][1]); //maybe it's better to use plain intgcd 924 922 925 926 printf ("newtonPolyg[0]= %d\n", newtonPolyg[0][0]);927 printf ("newtonPolyg[0]= %d\n", newtonPolyg[0][1]);928 out_cf ("g= ", g, "\n");929 930 923 int i= 1; 931 924 while (!g.isOne() && i < sizeOfNewtonPolygon) … … 933 926 g= gcd (g, newtonPolyg[i][0]); 934 927 g= gcd (g, newtonPolyg[i][1]); 935 printf ("newtonPolyg[i]= %d\n", newtonPolyg[i][0]);936 printf ("newtonPolyg[i]= %d\n", newtonPolyg[i][1]);937 out_cf ("g= ", g, "\n");938 928 i++; 939 929 } … … 1039 1029 ASSERT (factorize (F).length() <= 2, " expected irreducible polynomial"); 1040 1030 1041 out_cf ("F= ", F, "\n");1042 1031 bool isRat= isOn (SW_RATIONAL); 1043 1032 if (isRat) … … 1064 1053 Fp= Fp (y+E[2], y); 1065 1054 1066 out_cf ("Fp0= ", Fp, "\n");1067 1068 1055 if (tdeg == totaldegree (Fp)) 1069 1056 { … … 1073 1060 if (factors.length() == 2 && factors.getLast().exp() == 1) 1074 1061 { 1075 printf ("irred1\n");1076 1062 if (isRat) 1077 1063 On (SW_RATIONAL); … … 1087 1073 Fp= Fp (y+E[2], y); 1088 1074 1089 out_cf ("Fp1= ", Fp, "\n");1090 1091 1075 if (tdeg == totaldegree (Fp)) 1092 1076 { … … 1096 1080 if (factors.length() == 2 && factors.getLast().exp() == 1) 1097 1081 { 1098 printf ("irred2\n");1099 1082 if (isRat) 1100 1083 On (SW_RATIONAL); … … 1126 1109 if (factors.length() == 2 && factors.getLast().exp() == 1) 1127 1110 { 1128 out_cf ("Fp2= ", Fp, "\n");1129 printf ("tdeg= %d\n", tdeg);1130 printf ("totaldegree (Fp)= %d\n", totaldegree (Fp));1131 printf ("degree (Fp, 1)= %d\n", degree (Fp,1));1132 printf ("degree (Fp,2)= %d\n", degree (Fp,2));1133 printf ("degree (F,1)= %d\n", degree (F,1));1134 printf ("degree (F,2)= %d\n", degree (F,2));1135 out_cf ("factors.getFirst()= ", factors.getFirst().factor(), "\n");1136 printf ("char= %d\n", getCharacteristic ());1137 printf ("irred3\n");1138 1111 if (isRat) 1139 1112 On (SW_RATIONAL);
Note: See TracChangeset
for help on using the changeset viewer.