===================================================================
RCS file: /home/cvs/OpenXM/src/kan96xx/Kan/output.c,v
retrieving revision 1.3
retrieving revision 1.7
diff -u -p -r1.3 -r1.7
--- OpenXM/src/kan96xx/Kan/output.c	2001/05/04 01:06:24	1.3
+++ OpenXM/src/kan96xx/Kan/output.c	2011/10/05 05:46:14	1.7
@@ -1,5 +1,8 @@
-/* $OpenXM: OpenXM/src/kan96xx/Kan/output.c,v 1.2 2000/01/16 07:55:40 takayama Exp $ */
+/* $OpenXM: OpenXM/src/kan96xx/Kan/output.c,v 1.6 2006/12/21 05:29:49 takayama Exp $ */
 #include <stdio.h>
+#include <stdlib.h>
+#include <string.h>
+#include <ctype.h>
 #include "datatype.h"
 #include "stackm.h"
 #include "extern.h"
@@ -56,6 +59,8 @@ char *POLYToString(f,multSym,brace)
      int brace;
 {
   extern int Wrap;
+  extern int UseDsmall;
+  extern int COutput;
   int i,j,jj,fi;
   int printed = 0;
   int vi; /* index for variables */
@@ -84,7 +89,9 @@ char *POLYToString(f,multSym,brace)
     printed = 1;
     xnames = dnames = (char **)NULL;
   }else{
-    ringp = f->m->ringp; xnames = f->m->ringp->x; dnames = f->m->ringp->D;
+    ringp = f->m->ringp; xnames = f->m->ringp->x; 
+    if (UseDsmall) dnames = f->m->ringp->Dsmall;
+    else dnames = f->m->ringp->D;
     n = ringp->n;
     xout = ringp->outputOrder;
   }
@@ -125,34 +132,67 @@ char *POLYToString(f,multSym,brace)
     for (jj=0; jj<n*2; jj++) {
       j = xout[jj];
       if (j <n) {
-        if (f->m->e[j].x) {
-          vi++;
-          if (vi != 1) putstr(multStr(multSym),&oq);
-          putstr(xnames[j],&oq);
-          if (f->m->e[j].x >= 2) {
-            putstr("^",&oq);
-            putstr(intToString(f->m->e[j].x),&oq);
-          }else if (f->m->e[j].x < 0) {
-            putstr("^(",&oq);
-            putstr(intToString(f->m->e[j].x),&oq);
-            putstr(")",&oq);
-          }
-        }
+		if (COutput) {
+		  if (f->m->e[j].x) {
+			vi++;
+			if (vi != 1) putstr(multStr(multSym),&oq);
+			if (f->m->e[j].x == 1) {
+			  putstr(xnames[j],&oq);
+			} else {
+			  putstr("sm1power(",&oq);
+			  putstr(xnames[j],&oq);
+			  putstr(",",&oq);
+			  putstr(intToString(f->m->e[j].x),&oq);
+			  putstr(")",&oq);
+			}
+		  }
+		}else{
+		  if (f->m->e[j].x) {
+			vi++;
+			if (vi != 1) putstr(multStr(multSym),&oq);
+			putstr(xnames[j],&oq);
+			if (f->m->e[j].x >= 2) {
+			  putstr("^",&oq);
+			  putstr(intToString(f->m->e[j].x),&oq);
+			}else if (f->m->e[j].x < 0) {
+			  putstr("^(",&oq);
+			  putstr(intToString(f->m->e[j].x),&oq);
+			  putstr(")",&oq);
+			}
+		  }
+		}
       }else {
-        j = j-n;  
-        if (f->m->e[j].D) {
-          vi++;
-          if (vi != 1) putstr(multStr(multSym),&oq);
-          putstr(dnames[j],&oq);
-          if (f->m->e[j].D >= 2) {
-            putstr("^",&oq);
-            putstr(intToString(f->m->e[j].D),&oq);
-          }else if (f->m->e[j].D < 0) {
-            putstr("^(",&oq);
-            putstr(intToString(f->m->e[j].D),&oq);
-            putstr(")",&oq);
-          }
-        }
+		if (COutput) {
+		  j = j-n;  
+		  if (f->m->e[j].D) {
+			vi++;
+			if (vi != 1) putstr(multStr(multSym),&oq);
+			if (f->m->e[j].D == 1) {
+			  putstr(dnames[j],&oq);
+			}else {
+			  if (f->m->e[j].D) putstr("sm1power(",&oq);
+			  putstr(dnames[j],&oq);
+			  putstr(",",&oq);
+			  putstr(intToString(f->m->e[j].D),&oq);
+			  putstr(")",&oq);
+			}
+		  }
+		}else{
+		  j = j-n;  
+		  if (f->m->e[j].D) {
+			vi++;
+			if (vi != 1) putstr(multStr(multSym),&oq);
+			putstr(dnames[j],&oq);
+			if (f->m->e[j].D >= 2) {
+			  putstr("^",&oq);
+			  putstr(intToString(f->m->e[j].D),&oq);
+			}else if (f->m->e[j].D < 0) {
+			  putstr("^(",&oq);
+			  putstr(intToString(f->m->e[j].D),&oq);
+			  putstr(")",&oq);
+			}
+		  }
+		}
       }
     }
     fi++;
@@ -257,9 +297,34 @@ isConstant(f)
   return(1);
 }
 
+int isConstantAll(POLY f)
+{
+  int i;
+  int n;
+  if (f == POLYNULL) return(1);
+  while (f != POLYNULL) {
+	if (!isConstant(f)) return 0;
+	f = f->next;
+  }
+  return 1;
+}
+
 void errorOutput(s)
      char *s;
 {
   fprintf(stderr,"Error(output.c):%s\n",s);
   exit(15);
+}
+
+char **makeDsmall(char **dvars,int n) {
+  char **ans;
+  int i;
+  ans = (char **) sGC_malloc(sizeof(char *)*(n+1));
+  for (i=0; i<n; i++) {
+	ans[i] =  (char *) sGC_malloc(sizeof(char)*(strlen(dvars[i])+1));
+	strcpy(ans[i],dvars[i]);
+	ans[i][0] = tolower(ans[i][0]);
+  }
+  ans[n] = NULL;
+  return ans;
 }