source: git/Singular/ipprint.cc @ c7ae4d

fieker-DuValspielwiese
Last change on this file since c7ae4d was d2fc5c5, checked in by Hans Schoenemann <hannes@…>, 8 years ago
removed unused test for QRING_CMD
  • Property mode set to 100644
File size: 10.1 KB
Line 
1/****************************************
2*  Computer Algebra System SINGULAR     *
3****************************************/
4/*
5* ABSTRACT: interpreter: printing
6*/
7
8#include <kernel/mod2.h>
9#include <omalloc/omalloc.h>
10
11#include <misc/intvec.h>
12
13#include <polys/matpol.h>
14
15#include <kernel/polys.h>
16#include <kernel/ideals.h>
17
18#include "tok.h"
19#include "ipid.h"
20#include "subexpr.h"
21#include "ipshell.h"
22#include "ipprint.h"
23#include "attrib.h"
24
25/*2
26* print for: int, string, poly, vector, ideal
27*/
28/*2
29* print for: intvec
30*/
31static BOOLEAN ipPrint_INTVEC(leftv u)
32{
33  intvec *v=(intvec *)u->Data();
34  v->show();
35  PrintLn();
36  return FALSE;
37}
38
39/*2
40* print for: intmat
41*/
42static BOOLEAN ipPrint_INTMAT(leftv u)
43{
44  intvec *v=(intvec *)u->Data();
45  int i,j;
46  for(i=0;i<v->rows();i++)
47  {
48    for(j=0;j<v->cols();j++)
49    {
50      Print(" %5d",IMATELEM(*v,i+1,j+1));
51    }
52    PrintLn();
53  }
54  return FALSE;
55}
56
57/*2
58* internal print for: matrix
59*/
60static void ipPrint_MA0(matrix m, const char *name)
61{
62  if ((MATCOLS(m)>0)&&(MATROWS(m)>0))
63  {
64    char **s=(char **)omAlloc(MATCOLS(m)*MATROWS(m)*sizeof(char*));
65    char *ss;
66    int *l=(int *)omAlloc0(MATCOLS(m)*sizeof(int));
67    int i,j,k;
68    int vl=si_max(colmax/MATCOLS(m),8);
69
70    /* make enough space for the "largest" name*/
71    ss=(char *)omAlloc(14+strlen(name));
72    sprintf(ss,"%s[%d,%d]",name,MATCOLS(m),MATROWS(m));
73    vl=si_max(vl,(int)strlen(ss));
74    omFree(ss);
75
76    /* convert all polys to string */
77    i=MATCOLS(m)*MATROWS(m)-1;
78    ss=pString(m->m[i]);
79    if ((int)strlen(ss)>colmax) { s[i]=NULL; omFree(ss); }
80    else                        s[i]=ss;
81    for(i--;i>=0;i--)
82    {
83      StringSetS("");
84      pString0(m->m[i]);
85      StringAppendS(",");
86      ss=StringEndS();
87      if ((int)strlen(ss)>colmax) s[i]=NULL;
88      else                        s[i]=ss;
89    }
90    /* look up the width of all columns, put it in l[col_nr] */
91    /* insert names for very long entries */
92    for(i=MATROWS(m)-1;i>=0;i--)
93    {
94      for(j=MATCOLS(m)-1;j>=0;j--)
95      {
96        if (s[i*MATCOLS(m)+j]==NULL)
97        {
98          ss=(char *)omAlloc(14+strlen(name));
99          s[i*MATCOLS(m)+j]=ss;
100          ss[0]='\0';
101          sprintf(ss,"%s[%d,%d]",name,i+1,j+1);
102          if ((i!=MATROWS(m)-1) || (j!=MATCOLS(m)-1))
103          {
104            strcat(ss,",");
105            vl=si_max(vl,(int)strlen(ss));
106          }
107        }
108        k=strlen(s[i*MATCOLS(m)+j]);
109        if (k>l[j]) l[j]=k;
110      }
111    }
112    /* does it fit on a line ? */
113    int maxlen=0;
114    for(j=MATCOLS(m)-1;j>=0;j--)
115    {
116      maxlen+=l[j];
117    }
118    if (maxlen>colmax)
119    {
120      /* NO, it does not fit, so retry: */
121      /* look up the width of all columns, clear very long entriess */
122      /* put length in l[col_nr] */
123      /* insert names for cleared entries */
124      for(j=MATCOLS(m)-1;j>=0;j--)
125      {
126        for(i=MATROWS(m)-1;i>=0;i--)
127        {
128          k=strlen(s[i*MATCOLS(m)+j]);
129          if (/*strlen(s[i*MATCOLS(m)+j])*/ k > vl)
130          {
131            omFree((ADDRESS)s[i*MATCOLS(m)+j]);
132            ss=(char *)omAlloc(14+strlen(name));
133            s[i*MATCOLS(m)+j]=ss;
134            ss[0]='\0';
135            sprintf(ss,"%s[%d,%d]",name,i+1,j+1);
136            if ((i!=MATROWS(m)-1) || (j!=MATCOLS(m)-1))
137            {
138              strcat(ss,",");
139            }
140            l[j]=strlen(s[i*MATCOLS(m)+j]);
141            if (l[j]>vl)
142            {
143//#ifdef TEST
144//              PrintS("pagewidth too small in print(matrix)\n");
145//#endif
146              vl=l[j]; /* make large names fit*/
147            }
148            i=MATROWS(m);
149          }
150          else
151          {
152            if (k>l[j]) l[j]=k;
153          }
154        }
155      }
156    }
157    /*output of the matrix*/
158    for(i=0;i<MATROWS(m);i++)
159    {
160      k=l[0];
161      Print("%-*.*s",l[0],l[0],s[i*MATCOLS(m)]);
162      omFree(s[i*MATCOLS(m)]);
163      for(j=1;j<MATCOLS(m);j++)
164      {
165        if (k+l[j]>colmax)
166        {
167          PrintS("\n  ");
168          k=2;
169        }
170        k+=l[j];
171        Print("%-*.*s",l[j],l[j],s[i*MATCOLS(m)+j]);
172        omFree(s[i*MATCOLS(m)+j]);
173      }
174      PrintLn();
175    }
176    /* clean up */
177    omFreeSize((ADDRESS)s,MATCOLS(m)*MATROWS(m)*sizeof(char*));
178    omFreeSize((ADDRESS)l,MATCOLS(m)*sizeof(int));
179  }
180  else Print("%d x %d zero matrix\n",MATROWS(m),MATCOLS(m));
181}
182
183/*2
184* print for: matrix
185*/
186static BOOLEAN ipPrint_MA(leftv u)
187{
188  matrix m=(matrix)u->Data();
189  ipPrint_MA0(m,u->Name());
190  return FALSE;
191}
192
193/*2
194* print for: ring
195*/
196static BOOLEAN ipPrint_RING(leftv u)
197{
198  ring r=(ring)u->Data();
199  PrintS("polynomial ring, over a ");
200  if (rField_is_Ring(r))
201  {
202    if (rField_is_Domain(r)) PrintS("domain");
203    else                     PrintS("ring (with zero-divisors)");
204  }
205  else PrintS("field");
206  if (r->OrdSgn==1) PrintS(", global"); else PrintS(", local/mixed");
207  PrintS(" ordering\n");
208  rWrite(r, TRUE);
209  return FALSE;
210}
211
212#ifdef SINGULAR_4_1
213static BOOLEAN ipPrint_CRING(leftv u)
214{
215  coeffs r=(coeffs)u->Data();
216  if (nCoeff_is_Ring(r))
217  {
218    if (nCoeff_is_Domain(r)) PrintS("domain: ");
219    else                     PrintS("ring (with zero-divisors): ");
220  }
221  else PrintS("field: ");
222  PrintS(nCoeffName(r));
223  return FALSE;
224}
225#endif
226/*2
227* print for: vector
228*/
229static BOOLEAN ipPrint_V(leftv u)
230{
231  polyset m=NULL;
232  int l,j;
233  /*convert into an array of the components*/
234  p_Vec2Polys((poly)u->Data(), &m, &l, currRing);
235  /*output*/
236  PrintS("[");
237  j=0;
238  loop
239  {
240    PrintS(pString(m[j]));
241    j++;
242    if (j<l) PrintS(",");
243    else
244    {
245      PrintS("]\n");
246      break;
247    }
248  }
249  /* clean up */
250  for(j=l-1;j>=0;j--) pDelete(&m[j]);
251  omFreeSize((ADDRESS)m,l*sizeof(poly));
252  return FALSE;
253}
254
255BOOLEAN jjPRINT(leftv res, leftv u)
256{
257  SPrintStart();
258  BOOLEAN bo=FALSE;
259  switch(u->Typ())
260  {
261      case INTVEC_CMD:
262        bo=ipPrint_INTVEC(u);
263        break;
264
265      case INTMAT_CMD:
266        bo=ipPrint_INTMAT(u);
267        break;
268
269      case MATRIX_CMD:
270        bo=ipPrint_MA(u);
271        break;
272
273      case IDEAL_CMD:
274      {
275        char* s = u->String(NULL, FALSE, 2);
276        PrintS(s);
277        PrintLn();
278        omFree(s);
279        break;
280      }
281
282      case MODUL_CMD:
283      {
284        matrix m = id_Module2Matrix(id_Copy((ideal) u->Data(),currRing),currRing);
285        ipPrint_MA0(m, u->Name());
286        id_Delete((ideal *) &m,currRing);
287        break;
288      }
289
290      case VECTOR_CMD:
291        bo=ipPrint_V(u);
292        break;
293
294      case RING_CMD:
295        bo=ipPrint_RING(u);
296        break;
297
298      #ifdef SINGULAR_4_1
299      case CRING_CMD:
300        bo=ipPrint_CRING(u);
301        break;
302      #endif
303
304      default:
305        u->Print();
306        break;
307  }
308  char *s=SPrintEnd();
309  if (u->next==NULL)
310  {
311    int l=strlen(s);
312    if (s[l-1]=='\n') s[l-1]='\0';
313  }
314  res->data=(void*)s;
315  return bo;
316}
317
318
319/*2
320* dbprint
321*/
322BOOLEAN jjDBPRINT(leftv res, leftv u)
323{
324  BOOLEAN print=(printlevel>myynest);
325  if ((u->next!=NULL)&&(u->Typ()==INT_CMD))
326  {
327    print=  (((int)((long)(u->Data()))) > 0);
328    u=u->next;
329  }
330  if (print)
331  {
332    // BOOLEAN r=FALSE;
333    leftv h=u;
334    leftv hh;
335    while (h!=NULL)
336    {
337      hh=h->next;
338      h->next=NULL;
339      if (jjPRINT(res, h)) return TRUE;
340      PrintS((char*)res->data);
341      omFree(res->data);
342      PrintLn();
343      h->next=hh;
344      h=hh;
345    }
346  }
347  return FALSE;
348}
349
350static void ipPrintBetti(leftv u)
351{
352  int i,j;
353  int row_shift=(int)((long)(atGet(u,"rowShift",INT_CMD)));
354  intvec * betti=(intvec *)u->Data();
355  // head line --------------------------------------------------------
356  PrintS("      "); // 6 spaces for no. and :
357  for(j=0;j<betti->cols();j++) Print(" %5d",j); // 6 spaces pro column
358  PrintS("\n------"); // 6 spaces for no. and :
359  for(j=0;j<betti->cols();j++) PrintS("------"); // 6 spaces pro column
360  PrintLn();
361  // the table --------------------------------------------------------
362  for(i=0;i<betti->rows();i++)
363  {
364    Print("%5d:",i+row_shift);
365    for(j=1;j<=betti->cols();j++)
366    {
367      int m=IMATELEM(*betti,i+1,j);
368      if (m==0)
369        PrintS("     -");
370      else
371        Print(" %5d",m);
372    }
373    PrintLn();
374  }
375  // sum --------------------------------------------------------------
376  PrintS("------"); // 6 spaces for no. and :
377  for(j=0;j<betti->cols();j++) PrintS("------"); // 6 spaces pro column
378  PrintS("\ntotal:");
379  for(j=0;j<betti->cols();j++)
380  {
381    int s=0;
382    for(i=0;i<betti->rows();i++)
383    {
384      s+=IMATELEM(*betti,i+1,j+1);
385    }
386    Print(" %5d",s); // 6 spaces pro column
387  }
388  PrintLn();
389}
390
391
392/*2
393* print(...,"format")
394*/
395BOOLEAN jjPRINT_FORMAT(leftv res, leftv u, leftv v)
396{
397/* ==================== betti ======================================== */
398  if ((u->Typ()==INTMAT_CMD)&&(strcmp((char *)v->Data(),"betti")==0))
399  {
400    SPrintStart();
401    ipPrintBetti(u);
402    char *s = SPrintEnd();
403    s[strlen(s)]='\0';
404    res->data=s;
405  }
406  else
407/* ======================== end betti ================================= */
408  {
409    char* ns = omStrDup((char*) v->Data());
410    int dim = 1;
411    if (strlen(ns) == 3 && ns[1] == '2')
412    {
413      dim = 2;
414      ns[1] = ns[2];
415      ns[2] = '\0';
416    }
417    if (strcmp(ns,"%l") == 0)
418    {
419      res->data = (char*) u->String(NULL, TRUE, dim);
420      if (dim == 2)
421      {
422        char* ns = (char*) omAlloc(strlen((char*) res->data) + 2);
423        strcpy(ns, (char*) res->data);
424        omFree(res->data);
425        strcat(ns, "\n");
426        res->data = ns;
427      }
428    }
429    else if (strcmp(ns,"%t") == 0)
430    {
431      SPrintStart();
432      type_cmd(u);
433      res->data = SPrintEnd();
434      if (dim != 2)
435        ((char*)res->data)[strlen((char*)res->data) -1] = '\0';
436    }
437    else if (strcmp(ns,"%;") == 0)
438    {
439      SPrintStart();
440      u->Print();
441      if (dim == 2) PrintLn();
442      res->data = SPrintEnd();
443    }
444    else if  (strcmp(ns,"%p") == 0)
445    {
446      iiExprArith1(res, u, PRINT_CMD);
447    }
448    else if (strcmp(ns,"%b") == 0 && (u->Typ()==INTMAT_CMD))
449    {
450      SPrintStart();
451      ipPrintBetti(u);
452      if (dim == 2) PrintLn();
453      res->data = SPrintEnd();
454    }
455    else
456    {
457      res->data = u->String(NULL, FALSE, dim);
458      if (dim == 2)
459      {
460        char* ns = (char*) omAlloc(strlen((char*) res->data) + 2);
461        strcpy(ns, (char*) res->data);
462        omFree(res->data);
463        strcat(ns, "\n");
464        res->data = ns;
465      }
466    }
467    omFree(ns);
468  }
469  return FALSE;
470}
Note: See TracBrowser for help on using the repository browser.