Changeset 17a710 in git for factory/cfNewtonPolygon.cc

Timestamp:

Apr 19, 2013, 3:34:15 PM (11 years ago)

Author:

Martin Lee <martinlee84@…>

Branches:

(u'spielwiese', '4a9821a93ffdc22a6696668bd4f6b8c9de3e6c5f')

Children:

f7ed7dc2ce876bee21478ae1bf88e2734992bcd8

Parents:

efd410fc3eac114fbca114eb66d8070b61317f92

git-author:

Martin Lee <martinlee84@web.de>2013-04-19 15:34:15+02:00

git-committer:

Martin Lee <martinlee84@web.de>2013-05-02 11:42:37+02:00

Message:

chg: debug output and experimental code

File:

• factory/cfNewtonPolygon.cc (modified) (9 diffs)

factory/cfNewtonPolygon.cc

@@ -23 +23 @@
 #include "templates/ftmpl_functions.h"
 #include "algext.h"
+
+void out_cf(const char *s1,const CanonicalForm &f,const char *s2);
 
 static
@@ -921 +923 @@
   CanonicalForm g= gcd (newtonPolyg[0][0], newtonPolyg[0][1]); //maybe it's better to use plain intgcd
 
+
+  printf ("newtonPolyg[0]= %d\n", newtonPolyg[0][0]);
+  printf ("newtonPolyg[0]= %d\n", newtonPolyg[0][1]);
+  out_cf ("g= ", g, "\n");
+
   int i= 1;
   while (!g.isOne() && i < sizeOfNewtonPolygon)
@@ -926 +933 @@
     g= gcd (g, newtonPolyg[i][0]);
     g= gcd (g, newtonPolyg[i][1]);
+  printf ("newtonPolyg[i]= %d\n", newtonPolyg[i][0]);
+  printf ("newtonPolyg[i]= %d\n", newtonPolyg[i][1]);
+    out_cf ("g= ", g, "\n");
     i++;
   }
@@ -1029 +1039 @@
   ASSERT (factorize (F).length() <= 2, " expected irreducible polynomial");
 
+  out_cf ("F= ", F, "\n");
   bool isRat= isOn (SW_RATIONAL);
   if (isRat)
@@ -1053 +1064 @@
   Fp= Fp (y+E[2], y);
 
+  out_cf ("Fp0= ", Fp, "\n");
+
   if (tdeg == totaldegree (Fp))
   {
@@ -1060 +1073 @@
       if (factors.length() == 2 && factors.getLast().exp() == 1)
       {
+        printf ("irred1\n");
         if (isRat)
           On (SW_RATIONAL);
@@ -1073 +1087 @@
   Fp= Fp (y+E[2], y);
 
+  out_cf ("Fp1= ", Fp, "\n");
+
   if (tdeg == totaldegree (Fp))
   {
@@ -1080 +1096 @@
       if (factors.length() == 2 && factors.getLast().exp() == 1)
       {
+        printf ("irred2\n");
         if (isRat)
           On (SW_RATIONAL);
@@ -1109 +1126 @@
           if (factors.length() == 2 && factors.getLast().exp() == 1)
           {
+            out_cf ("Fp2= ", Fp, "\n");
+            printf ("tdeg= %d\n", tdeg);
+            printf ("totaldegree (Fp)= %d\n", totaldegree (Fp));
+            printf ("degree (Fp, 1)= %d\n", degree (Fp,1));
+            printf ("degree (Fp,2)= %d\n", degree (Fp,2));
+            printf ("degree (F,1)= %d\n", degree (F,1));
+            printf ("degree (F,2)= %d\n", degree (F,2));
+            out_cf ("factors.getFirst()= ", factors.getFirst().factor(), "\n");
+            printf ("char= %d\n", getCharacteristic ());
+            printf ("irred3\n");
             if (isRat)
               On (SW_RATIONAL);