source: git/Singular/ssiLink.cc @ e8c21d

spielwiese
Last change on this file since e8c21d was e8c21d, checked in by Hans Schoenemann <hannes@…>, 14 years ago
better name git-svn-id: file:///usr/local/Singular/svn/trunk@13137 2c84dea3-7e68-4137-9b89-c4e89433aadc
  • Property mode set to 100644
File size: 26.5 KB
Line 
1/****************************************
2 * Computer Algebra System SINGULAR     *
3 ****************************************/
4/***************************************************************
5 * File:    ssiLink.h
6 *  Purpose: declaration of sl_link routines for ssi
7 *  Version: $Id$
8 ***************************************************************/
9#include <stdio.h>
10#include <fcntl.h>
11#include <errno.h>
12#include <unistd.h>
13#include <sys/types.h>
14#include <signal.h>
15#include <sys/types.h>          /* for portability */
16#include <sys/select.h>
17#include <sys/socket.h>
18#include <netinet/in.h>
19#include <ctype.h>   /*for isdigit*/
20#include <unistd.h>
21#include <netdb.h>
22
23
24#include "mod2.h"
25
26#include "tok.h"
27#include "ipid.h"
28#include <omalloc.h>
29#include <kernel/ring.h>
30#include <kernel/matpol.h>
31#include <kernel/ideals.h>
32#include <kernel/polys.h>
33#include <kernel/longrat.h>
34#include <kernel/ideals.h>
35#include "subexpr.h"
36#include "silink.h"
37#include "lists.h"
38#include "ssiLink.h"
39
40typedef struct
41{
42  FILE *f_read;
43  FILE *f_write;
44  ring r;
45  pid_t pid; /* only valid for fork/tcp mode*/
46  int fd_read,fd_write; /* only valid for fork/tcp mode*/
47  char level;
48} ssiInfo;
49
50// the helper functions:
51void ssiSetCurrRing(const ring r)
52{
53  if (!rEqual(r,currRing,1))
54  {
55    char name[20];
56    int nr=0;
57    do
58    { sprintf(name,"ssiRing%d",nr); nr++; }
59    while(IDROOT->get(name, 0)!=NULL);
60    idhdl h=enterid(omStrDup(name),0,RING_CMD,&IDROOT,FALSE);
61    IDRING(h)=r;
62    r->ref++;
63    rSetHdl(h);
64  }
65}
66// the implementation of the functions:
67void ssiWriteInt(ssiInfo *d,const int i)
68{
69  fprintf(d->f_write,"%d ",i);
70  //if (d->f_debug!=NULL) fprintf(d->f_debug,"int: %d ",i);
71}
72
73void ssiWriteString(ssiInfo *d,const char *s)
74{
75  fprintf(d->f_write,"%d %s ",strlen(s),s);
76  //if (d->f_debug!=NULL) fprintf(d->f_debug,"stringi: %d \"%s\" ",strlen(s),s);
77}
78
79
80void ssiWriteBigInt(const ssiInfo *d, const number n)
81{
82  // syntax is as follows:
83  // case 2 Q:     3 4 <int>
84  //        or     3 3 <mpz_t nominator>
85  if(SR_HDL(n) & SR_INT)
86  {
87    fprintf(d->f_write,"4 %ld ",SR_TO_INT(n));
88    //if (d->f_debug!=NULL) fprintf(d->f_debug,"bigint: short \"%ld\" ",SR_TO_INT(n));
89  }
90  else if (n->s==3)
91  {
92    gmp_fprintf(d->f_write,"3 %Zd ",n->z);
93    //if (d->f_debug!=NULL) gmp_fprintf(d->f_debug,"bigint: gmp \"%Zd\" ",n->z);
94  }
95  else Werror("illiegal bigint");
96}
97
98void ssiWriteNumber(const ssiInfo *d, const number n)
99{
100  // syntax is as follows:
101  // case 1 Z/p:   3 <int>
102  // case 2 Q:     3 4 <int>
103  //        or     3 0 <mpz_t nominator> <mpz_t denominator>
104  //        or     3 1  dto.
105  //        or     3 3 <mpz_t nominator>
106  if(rField_is_Zp(d->r))
107  {
108    fprintf(d->f_write,"%d ",(int)(long)n);
109    //if (d->f_debug!=NULL) fprintf(d->f_debug,"number: \"%ld\" ",(int)(long)n);
110  }
111  else if (rField_is_Q(d->r))
112  {
113    if(SR_HDL(n) & SR_INT)
114    {
115      fprintf(d->f_write,"4 %ld ",SR_TO_INT(n));
116      //if (d->f_debug!=NULL) fprintf(d->f_debug,"number: short \"%ld\" ",SR_TO_INT(n));
117    }
118    else if (n->s<2)
119    {
120      gmp_fprintf(d->f_write,"%d %Zd %Zd ",n->s,n->z,n->n);
121      //if (d->f_debug!=NULL) gmp_fprintf(d->f_debug,"number: s=%d gmp/gmp \"%Zd %Zd\" ",n->s,n->z,n->n);
122    }
123    else /*n->s==3*/
124    {
125      gmp_fprintf(d->f_write,"3 %Zd ",n->z);
126      //if (d->f_debug!=NULL) gmp_fprintf(d->f_debug,"number: gmp \"%Zd\" ",n->z);
127    }
128  }
129  else WerrorS("coeff field not implemented");
130}
131
132void ssiWriteRing(ssiInfo *d,const ring r)
133{
134  /* 5 <ch> <N> <l1> <s1>....<ln> <sN> */
135  if (d->r!=NULL) rKill(d->r);
136  d->r=r;
137  d->r->ref++;
138  fprintf(d->f_write,"%d %d ",r->ch,r->N);
139
140  int i;
141  for(i=0;i<r->N;i++)
142  {
143    fprintf(d->f_write,"%d %s ",strlen(r->names[i]),r->names[i]);
144  }
145  /* number of orderings:*/
146  i=0;
147  while (r->order[i]!=0) i++;
148  fprintf(d->f_write,"%d ",i);
149  /* each ordering block: */
150  i=0;
151  while(r->order[i]!=0)
152  {
153    fprintf(d->f_write,"%d %d %d ",r->order[i],r->block0[i], r->block1[i]);
154    i++;
155  }
156}
157
158void ssiWritePoly(ssiInfo *d, int typ, poly p)
159{
160  fprintf(d->f_write,"%d ",pLength(p));//number of terms
161  int i;
162
163  while(p!=NULL)
164  {
165    ssiWriteNumber(d,pGetCoeff(p));
166    //nWrite(fich,pGetCoeff(p));
167    fprintf(d->f_write,"%ld ",p_GetComp(p,d->r));//component
168
169    for(int j=1;j<=rVar(d->r);j++)
170    {
171      fprintf(d->f_write,"%ld ",p_GetExp(p,j,d->r ));//x^j
172    }
173    pIter(p);
174  }
175}
176
177void ssiWriteIdeal(ssiInfo *d, int typ,ideal I)
178{
179   // syntax: 7 # of elements <poly 1> <poly2>.....
180   // syntax: 8 <rows> <cols> <poly 1> <poly2>.....
181   matrix M=(matrix)I;
182   if (typ==MATRIX_CMD)
183        fprintf(d->f_write,"%d %d ", MATROWS(M),MATCOLS(M));
184   else
185     fprintf(d->f_write,"%d ",IDELEMS(I));
186
187   int i;
188   int tt;
189   if (typ==MODUL_CMD) tt=VECTOR_CMD;
190   else                tt=POLY_CMD;
191
192   for(i=0;i<IDELEMS(I);i++)
193   {
194     ssiWritePoly(d,tt,I->m[i]);
195   }
196}
197void ssiWriteCommand(si_link l, command D)
198{
199  ssiInfo *d=(ssiInfo*)l->data;
200  // syntax: <num ops> <operation> <op1> <op2> ....
201  fprintf(d->f_write,"%d %d ",D->argc,D->op);
202  if (D->argc >0) ssiWrite(l, &(D->arg1));
203  if (D->argc < 4)
204  {
205    if (D->argc >1) ssiWrite(l, &(D->arg2));
206    if (D->argc >2) ssiWrite(l, &(D->arg3));
207  }
208}
209
210void ssiWriteProc(ssiInfo *d,procinfov p)
211{
212  ssiWriteString(d,p->data.s.body);
213}
214
215void ssiWriteList(si_link l,lists dd)
216{
217  ssiInfo *d=(ssiInfo*)l->data;
218  fprintf(d->f_write,"%d ",dd->nr+1);
219  int i;
220  for(i=0;i<=dd->nr;i++)
221  {
222    ssiWrite(l,&(dd->m[i]));
223  }
224}
225
226char *ssiReadString(ssiInfo *d)
227{
228  char *buf;
229  int l;
230  fscanf(d->f_read,"%d ",&l);
231  buf=(char*)omAlloc(l+1);
232  fread(buf,1,l,d->f_read);
233  buf[l]='\0';
234  return buf;
235}
236
237int ssiReadInt(FILE *fich)
238{
239  int d;
240  fscanf(fich,"%d",&d);
241  return d;
242}
243
244number ssiReadBigInt(ssiInfo *d)
245{
246   int sub_type=-1;
247   fscanf(d->f_read,"%d",&sub_type);
248   switch(sub_type)
249   {
250   case 3:
251     {// read int or mpz_t or mpz_t, mpz_t
252       number n=nlRInit(0);
253       gmp_fscanf(d->f_read,"%Zd",n->z);
254       n->s=sub_type;
255       return n;
256     }
257   case 4: { int dd; fscanf(d->f_read,"%d",&dd); return INT_TO_SR(dd); }
258   default: Werror("error in reading number: invalid subtype %d",sub_type);
259            return NULL;
260   }
261}
262
263number ssiReadNumber(ssiInfo *d)
264{
265  if (rField_is_Q(d->r))
266  {
267     int sub_type=-1;
268     fscanf(d->f_read,"%d",&sub_type);
269     switch(sub_type)
270     {
271     case 0:
272     case 1:
273       {// read int or mpz_t or mpz_t, mpz_t
274        number n=nlRInit(0);
275        mpz_init(n->n);
276        gmp_fscanf(d->f_read,"%Zd %Zd",n->z,n->n);
277        n->s=sub_type;
278        return n;
279       }
280
281     case 3:
282       {// read int or mpz_t or mpz_t, mpz_t
283         number n=nlRInit(0);
284         gmp_fscanf(d->f_read,"%Zd",n->z);
285         n->s=sub_type;
286         return n;
287       }
288     case 4: { int dd; fscanf(d->f_read,"%d",&dd); return INT_TO_SR(dd); }
289     default: Werror("error in reading number: invalid subtype %d",sub_type);
290              return NULL;
291     }
292  }
293  else if (rField_is_Zp(d->r))
294  {
295    // read int
296    int dd;
297    fscanf(d->f_read,"%d",&dd);
298    return (number)dd;
299  }
300  else Werror("coeffs not implemented");
301  return NULL;
302}
303
304ring ssiReadRing(ssiInfo *d)
305{
306/* syntax is <ch> <N> <l1> <s1> ....<lN> <sN> */
307  int ch, N,i,l;
308  char **names;
309  fscanf(d->f_read,"%d %d ",&ch,&N);
310  names=(char**)omAlloc(N*sizeof(char*));
311  for(i=0;i<N;i++)
312  {
313    names[i]=ssiReadString(d);
314  }
315  // read the orderings:
316  int num_ord; // number of orderings
317  fscanf(d->f_read,"%d",&num_ord);
318  int *ord=(int *)omAlloc0((num_ord+1)*sizeof(int));
319  int *block0=(int *)omAlloc0((num_ord+1)*sizeof(int));
320  int *block1=(int *)omAlloc0((num_ord+1)*sizeof(int));
321  for(i=0;i<num_ord;i++)
322  {
323     fscanf(d->f_read,"%d %d %d",&ord[i],&block0[i],&block1[i]);
324  }
325  return rDefault(ch,N,names,num_ord,ord,block0,block1);
326}
327
328poly ssiReadPoly(ssiInfo *D)
329{
330// < # of terms> < term1> < .....
331  int n,i,l;
332  n=ssiReadInt(D->f_read);
333  //Print("poly: terms:%d\n",n);
334  poly p;
335  int j;
336  j=0;
337  poly ret=NULL;
338  poly prev=NULL;
339  for(l=0;l<n;l++) // read n terms
340  {
341// coef,comp.exp1,..exp N
342    p=p_Init(D->r);
343    pGetCoeff(p)=ssiReadNumber(D);
344    int d;
345    fscanf(D->f_read,"%d",&d);
346    p_SetComp(p,d,D->r);
347    for(i=1;i<=rVar(D->r);i++)
348    {
349      fscanf(D->f_read,"%d",&d);
350      p_SetExp(p,i,d,D->r);
351    }
352    p_Setm(p,D->r);
353    p_Test(p,D->r);
354    if (ret==NULL) ret=p;
355    else           pNext(prev)=p;
356    prev=p;
357 }
358 return ret;
359}
360
361ideal ssiReadIdeal(ssiInfo *d)
362{
363  int n,i;
364  ideal I;
365  fscanf(d->f_read,"%d",&n);
366  I=idInit(n,1);
367  for(i=0;i<IDELEMS(I);i++) // read n terms
368  {
369    I->m [i]=ssiReadPoly(d);
370  }
371  return I;
372}
373
374matrix ssiReadMatrix(ssiInfo *d)
375{
376  int n,m,i,j;
377  fscanf(d->f_read,"%d %d",&m,&n);
378  matrix M=mpNew(m,n);
379  poly p;
380  for(int i=1;i<=MATROWS(M);i++)
381    for(int j=1;j<=MATCOLS(M);j++)
382    {
383      p=ssiReadPoly(d);
384      MATELEM(M,i,j)=p;
385    }
386  return M;
387}
388
389command ssiReadCommand(si_link l)
390{
391  ssiInfo *d=(ssiInfo*)l->data;
392  // syntax: <num ops> <operation> <op1> <op2> ....
393  command D=(command)omAlloc0(sizeof(*D));
394  int argc,op;
395  fscanf(d->f_read,"%d %d",&argc,&op);
396  D->argc=argc; D->op=op;
397  leftv v;
398  if (argc >0)
399  {
400    v=ssiRead1(l);
401    memcpy(&(D->arg1),v,sizeof(*v));
402    omFreeBin(v,sleftv_bin);
403  }
404  if (argc <4)
405  {
406    if (D->argc >1)
407    {
408      v=ssiRead1(l);
409      memcpy(&(D->arg2),v,sizeof(*v));
410      omFreeBin(v,sleftv_bin);
411    }
412    if (D->argc >2)
413    {
414      v=ssiRead1(l);
415      memcpy(&(D->arg3),v,sizeof(*v));
416      omFreeBin(v,sleftv_bin);
417    }
418  }
419  else
420  {
421    leftv prev=&(D->arg1);
422    argc--;
423    while(argc >0)
424    {
425      v=ssiRead1(l);
426      prev->next=v;
427      prev=v;
428      argc--;
429    }
430  }
431  return D;
432}
433
434procinfov ssiReadProc(ssiInfo *d)
435{
436  char *s=ssiReadString(d);
437  procinfov p=(procinfov)omAlloc0Bin(procinfo_bin);
438  p->language=LANG_SINGULAR;
439  p->libname=omStrDup("");
440  p->procname=omStrDup("");
441  p->data.s.body=s;
442  return p;
443}
444lists ssiReadList(si_link l)
445{
446  ssiInfo *d=(ssiInfo*)l->data;
447  int nr;
448  fscanf(d->f_read,"%d",&nr);
449  lists L=(lists)omAlloc(sizeof(*L));
450  L->Init(nr);
451
452  int i;
453  leftv v;
454  for(i=0;i<nr;i++)
455  {
456    v=ssiRead1(l);
457    memcpy(&(L->m[i]),v,sizeof(*v));
458    omFreeBin(v,sleftv_bin);
459  }
460  return L;
461}
462
463//**************************************************************************/
464BOOLEAN ssiOpen(si_link l, short flag)
465{
466  const char *mode;
467  ssiInfo *d=(ssiInfo*)omAlloc0(sizeof(ssiInfo));
468  if (flag & SI_LINK_OPEN)
469  {
470    if (l->mode[0] != '\0' && (strcmp(l->mode, "r") == 0))
471      flag = SI_LINK_READ;
472    else flag = SI_LINK_WRITE;
473  }
474
475  if (flag == SI_LINK_READ) mode = "r";
476  else if (strcmp(l->mode, "w") == 0) mode = "w";
477  else if (strcmp(l->mode, "fork") == 0) mode = "fork";
478  else if (strcmp(l->mode, "tcp") == 0) mode = "tcp";
479  else if (strcmp(l->mode, "connect") == 0) mode = "connect";
480  else mode = "a";
481
482
483  if (l->name[0] == '\0')
484  {
485    if (strcmp(mode,"fork")==0)
486    {
487      int pc[2];
488      int cp[2];
489      pipe(pc);
490      pipe(cp);
491      pid_t pid=fork();
492      if (pid==0) /*child*/
493      {
494        close(pc[1]); close(cp[0]);
495        d->f_read=fdopen(pc[0],"r");
496        d->fd_read=pc[0];
497        d->f_write=fdopen(cp[1],"w");
498        d->fd_write=cp[1];
499        l->data=d;
500        omFree(l->mode);
501        l->mode = omStrDup(mode);
502        SI_LINK_SET_RW_OPEN_P(l);
503        myynest=0;
504        fe_fgets_stdin=fe_fgets_dummy;
505        loop
506        {
507          leftv h=ssiRead1(l); /*contains an exit.... */
508          if (feErrors != NULL && *feErrors != '\0')
509          {
510            // handle errors:
511            PrintS(feErrors); /* currently quite simple */
512            *feErrors = '\0';
513          }
514          ssiWrite(l,h);
515          h->CleanUp();
516          omFreeBin(h, sleftv_bin);
517        }
518        /* never reached*/
519      }
520      else if (pid>0)
521      {
522        d->pid=pid;
523        close(pc[0]); close(cp[1]);
524        d->f_read=fdopen(cp[0],"r");
525        d->fd_read=cp[0];
526        d->f_write=fdopen(pc[1],"w");
527        d->fd_write=pc[1];
528        SI_LINK_SET_RW_OPEN_P(l);
529      }
530      else
531      {
532        Werror("fork failed (%d)",errno);
533      }
534    }
535    else if (strcmp(mode,"tcp")==0)
536    {
537      int sockfd, newsockfd, portno, clilen;
538      struct sockaddr_in serv_addr, cli_addr;
539      int n;
540      sockfd = socket(AF_INET, SOCK_STREAM, 0);
541      if(sockfd < 0)
542      {
543        WerrorS("ERROR opening socket");
544        return TRUE;
545      }
546      memset((char *) &serv_addr,0, sizeof(serv_addr));
547      portno = 1024;
548      serv_addr.sin_family = AF_INET;
549      serv_addr.sin_addr.s_addr = INADDR_ANY;
550      do
551      {
552        portno++;
553        serv_addr.sin_port = htons(portno);
554        if(portno > 50000)
555        {
556          WerrorS("ERROR on binding (no free port available?)");
557          return TRUE;
558        }
559      }
560      while(bind(sockfd, (struct sockaddr *) &serv_addr, sizeof(serv_addr)) < 0);
561      Print("waiting on port %d\n", portno);mflush();
562      listen(sockfd,5);
563      newsockfd = accept(sockfd, (struct sockaddr *) &cli_addr, (socklen_t *)&clilen);
564      if(newsockfd < 0)
565      {
566        WerrorS("ERROR on accept");
567        return TRUE;
568      }
569      PrintS("client accepted\n");
570      d->fd_read = newsockfd;
571      d->fd_write = newsockfd;
572      d->f_read = fdopen(newsockfd, "r");
573      d->f_write = fdopen(newsockfd, "w");
574      SI_LINK_SET_RW_OPEN_P(l);
575      close(sockfd);
576    }
577    // stdin or stdout
578    else if (flag == SI_LINK_READ)
579    {
580      d->f_read = stdin;
581      mode = "r";
582    }
583    else
584    {
585      d->f_write = stdout;
586      mode = "a";
587    }
588  }
589  else /*l->name=NULL*/
590  {
591    // tcp mode
592    if(strcmp(mode,"tcp")==0)
593    {
594      int sockfd, newsockfd, portno, clilen;
595      struct sockaddr_in serv_addr, cli_addr;
596      int n;
597      sockfd = socket(AF_INET, SOCK_STREAM, 0);
598      if(sockfd < 0)
599      {
600        WerrorS("ERROR opening socket");
601        return TRUE;
602      }
603      memset((char *) &serv_addr,0, sizeof(serv_addr));
604      portno = 1024;
605      serv_addr.sin_family = AF_INET;
606      serv_addr.sin_addr.s_addr = INADDR_ANY;
607      do
608      {
609        portno++;
610        serv_addr.sin_port = htons(portno);
611        if(portno > 50000)
612        {
613          WerrorS("ERROR on binding (no free port available?)");
614          return TRUE;
615        }
616      }
617      while(bind(sockfd, (struct sockaddr *) &serv_addr, sizeof(serv_addr)) < 0);
618      //Print("waiting on port %d\n", portno);mflush();
619      listen(sockfd,5);
620      char* cli_host = (char*)omAlloc(256);
621      char* path = (char*)omAlloc(1024);
622      int r = sscanf(l->name,"%255[^:]:%s",cli_host,path);
623      if(r == 0)
624      {
625        WerrorS("ERROR: no host specified");
626        return TRUE;
627      }
628      else if(r == 1)
629      {
630        path = "Singular";
631      }
632      char* ssh_command = (char*)omAlloc(256);
633      char* ser_host = (char*)omAlloc(64);
634      gethostname(ser_host,64);
635      sprintf(ssh_command,"ssh %s %s -q --batch --link=ssi --MPhost=%s --MPport=%d &",cli_host,path,ser_host,portno);
636      system(ssh_command);
637      //Print("client on %s started:%s\n",cli_host,path);
638      clilen = sizeof(cli_addr);
639      newsockfd = accept(sockfd, (struct sockaddr *) &cli_addr, (socklen_t *)&clilen);
640      if(newsockfd < 0)
641      {
642        WerrorS("ERROR on accept");
643        return TRUE;
644      }
645      //PrintS("client accepted\n");
646      d->fd_read = newsockfd;
647      d->fd_write = newsockfd;
648      d->f_read = fdopen(newsockfd, "r");
649      d->f_write = fdopen(newsockfd, "w");
650      SI_LINK_SET_RW_OPEN_P(l);
651      close(sockfd);
652    }
653    else if(strcmp(mode,"connect")==0)
654    {
655      char* host = (char*)omAlloc(256);
656      int sockfd, portno, n;
657      struct sockaddr_in serv_addr;
658      struct hostent *server;
659 
660      sscanf(l->name,"%255[^:]:%d",host,&portno);
661      //Print("connect to host %s, port %d\n",host,portno);mflush();
662      if (portno!=0)
663      {
664        sockfd = socket(AF_INET, SOCK_STREAM, 0);
665        if (sockfd < 0) { WerrorS("ERROR opening socket"); return TRUE; }
666        server = gethostbyname(host);
667        if (server == NULL) {  WerrorS("ERROR, no such host");  return TRUE; }
668        memset((char *) &serv_addr, 0, sizeof(serv_addr));
669        serv_addr.sin_family = AF_INET;
670        memcpy((char *)&serv_addr.sin_addr.s_addr,
671              (char *)server->h_addr,
672              server->h_length);
673        serv_addr.sin_port = htons(portno);
674        if (connect(sockfd,(sockaddr*)&serv_addr,sizeof(serv_addr)) < 0)
675        { WerrorS("ERROR connecting"); return TRUE; }
676        //PrintS("connected\n");mflush();
677        d->f_read=fdopen(sockfd,"r");
678        d->fd_read=sockfd;
679        d->f_write=fdopen(sockfd,"w");
680        d->fd_write=sockfd;
681        l->data=d;
682        SI_LINK_SET_RW_OPEN_P(l);
683        omFree(host);
684      }
685      else return TRUE;
686    }
687    else
688    {
689      // normal ascii link to a file
690      FILE *outfile;
691      char *filename=l->name;
692
693      if(filename[0]=='>')
694      {
695        if (filename[1]=='>')
696        {
697          filename+=2;
698          mode = "a";
699        }
700        else
701        {
702          filename++;
703          mode="w";
704        }
705      }
706      outfile=myfopen(filename,mode);
707      if (outfile!=NULL)
708      {
709        if (strcmp(l->mode,"r")==0) d->f_read = outfile;
710        else d->f_write = outfile;
711      }
712      else
713      {
714        omFree(d);
715        return TRUE;
716      }
717    }
718  }
719  l->data=d;
720
721  omFree(l->mode);
722  l->mode = omStrDup(mode);
723  SI_LINK_SET_OPEN_P(l, flag);
724  return FALSE;
725}
726
727//**************************************************************************/
728LINKAGE BOOLEAN ssiClose(si_link l)
729{
730  ssiInfo *d = (ssiInfo *)l->data;
731  if (d!=NULL)
732  {
733    if ((strcmp(l->mode,"tcp")==0)
734    || (strcmp(l->mode,"fork")==0))
735    {
736      fprintf(d->f_write,"99\n");fflush(d->f_write);
737    }
738    if (d->f_read!=NULL) fclose(d->f_read);
739    if (d->f_write!=NULL) fclose(d->f_write);
740    if (d->r!=NULL) rKill(d->r);
741    if (d->pid!=0) { kill(d->pid,15); kill(d->pid,9); }
742    omFreeSize((ADDRESS)d,(sizeof *d));
743  }
744  l->data=NULL;
745  SI_LINK_SET_CLOSE_P(l);
746  return FALSE;
747}
748
749//**************************************************************************/
750LINKAGE leftv ssiRead1(si_link l)
751{
752  ssiInfo *d = (ssiInfo *)l->data;
753  leftv res=(leftv)omAlloc0(sizeof(sleftv));
754  int t=0;
755  fscanf(d->f_read,"%d",&t);
756  //Print("got type %d\n",t);
757  switch(t)
758  {
759    case 1:res->rtyp=INT_CMD;
760           res->data=(char *)ssiReadInt(d->f_read);
761           break;
762    case 2:res->rtyp=STRING_CMD;
763           res->data=(char *)ssiReadString(d);
764           break;
765    case 3:res->rtyp=NUMBER_CMD;
766           res->data=(char *)ssiReadNumber(d);
767           break;
768    case 4:res->rtyp=BIGINT_CMD;
769           res->data=(char *)ssiReadBigInt(d);
770           break;
771    case 15:
772    case 5:{
773             d->r=ssiReadRing(d);
774             d->r->ref++;
775             res->rtyp=RING_CMD;
776             res->data=(char*)d->r;
777             // we are in the top-level, so set the basering to d->r:
778             ssiSetCurrRing(d->r);
779             if (t==15) return ssiRead1(l);
780           }
781           break;
782    case 6:res->rtyp=POLY_CMD;
783           if (d->r==NULL) goto no_ring;
784           res->data=(char*)ssiReadPoly(d);
785           break;
786    case 7:res->rtyp=IDEAL_CMD;
787           if (d->r==NULL) goto no_ring;
788           res->data=(char*)ssiReadIdeal(d);
789           break;
790    case 8:res->rtyp=MATRIX_CMD;
791           if (d->r==NULL) goto no_ring;
792           res->data=(char*)ssiReadMatrix(d);
793           break;
794    case 9:res->rtyp=VECTOR_CMD;
795           if (d->r==NULL) goto no_ring;
796           res->data=(char*)ssiReadPoly(d);
797           break;
798    case 10:res->rtyp=MODUL_CMD;
799           if (d->r==NULL) goto no_ring;
800           res->data=(char*)ssiReadIdeal(d);
801           break;
802    case 11:
803           {
804             res->rtyp=COMMAND;
805             res->data=ssiReadCommand(l);
806             int nok=res->Eval();
807             if (nok) WerrorS("error in eval");
808             break;
809           }
810    case 12: /*DEF_CMD*/
811           {
812             res->rtyp=0;
813             res->name=(char *)ssiReadString(d);
814             int nok=res->Eval();
815             if (nok) WerrorS("error in name lookup");
816             break;
817           }
818    case 13: res->rtyp=PROC_CMD;
819             res->data=ssiReadProc(d);
820             break;
821    case 14: res->rtyp=LIST_CMD;
822             res->data=ssiReadList(l);
823             break;
824    case 16: res->rtyp=NONE; res->data=NULL;
825             break;
826    case 99: ssiClose(l); exit(0);
827    case 0: if (feof(d->f_read))
828            {
829              ssiClose(l);
830              res->rtyp=DEF_CMD;
831              break;
832            }
833    default: Werror("not implemented (t:%d)",t);
834             omFreeSize(res,sizeof(sleftv));
835             res=NULL;
836             break;
837  }
838  return res;
839no_ring: WerrorS("no ring");
840  return NULL;
841}
842//**************************************************************************/
843LINKAGE BOOLEAN ssiWrite(si_link l, leftv data)
844{
845  if(!SI_LINK_W_OPEN_P(l)) slOpen(l,SI_LINK_OPEN|SI_LINK_WRITE);
846  ssiInfo *d = (ssiInfo *)l->data;
847  d->level++;
848  //FILE *fich=d->f;
849  while (data!=NULL)
850  {
851    int tt=data->Typ();
852    void *dd=data->Data();
853
854    switch(tt /*data->Typ()*/)
855    {
856          case NONE/* nothing*/:fprintf(d->f_write,"16 ");
857                          break;
858          case STRING_CMD: fprintf(d->f_write,"2 ");
859                           ssiWriteString(d,(char *)dd);
860                           break;
861          case INT_CMD: fprintf(d->f_write,"1 ");
862                        ssiWriteInt(d,(int)(long)dd);
863                        break;
864          case BIGINT_CMD:fprintf(d->f_write,"4 ");
865                        ssiWriteBigInt(d,(number)dd);
866                        break;
867          case NUMBER_CMD:
868                          if (d->r!=currRing)
869                          {
870                            fprintf(d->f_write,"15 ");
871                            ssiWriteRing(d,currRing);
872                            if (d->level<=1) fprintf(d->f_write,"\n");
873                          }
874                          fprintf(d->f_write,"3 ");
875                          ssiWriteNumber(d,(number)dd);
876                        break;
877          case RING_CMD:fprintf(d->f_write,"5 ");
878                        ssiWriteRing(d,(ring)dd);
879                        break;
880          case POLY_CMD:
881          case VECTOR_CMD:
882                        if (d->r!=currRing)
883                        {
884                          fprintf(d->f_write,"15 ");
885                          ssiWriteRing(d,currRing);
886                          if (d->level<=1) fprintf(d->f_write,"\n");
887                        }
888                        if(tt==POLY_CMD) fprintf(d->f_write,"6 ");
889                        else             fprintf(d->f_write,"9 ");
890                        ssiWritePoly(d,tt,(poly)dd);
891                        break;
892          case IDEAL_CMD:
893          case MODUL_CMD:
894          case MATRIX_CMD:
895                        if (d->r!=currRing)
896                        {
897                          fprintf(d->f_write,"15 ");
898                          ssiWriteRing(d,currRing);
899                          if (d->level<=1) fprintf(d->f_write,"\n");
900                        }
901                        if(tt==IDEAL_CMD)       fprintf(d->f_write,"7 ");
902                        else if(tt==MATRIX_CMD) fprintf(d->f_write,"8 ");
903                        else                    fprintf(d->f_write,"10 ");
904                        ssiWriteIdeal(d,tt,(ideal)dd);
905                        break;
906          case COMMAND:
907                   fprintf(d->f_write,"11 ");
908                   ssiWriteCommand(l,(command)dd);
909                   break;
910          case DEF_CMD: /* not evaluated stuff in quotes */
911                   fprintf(d->f_write,"12 ");
912                   ssiWriteString(d,data->Name());
913                   break;
914          case PROC_CMD:
915                   fprintf(d->f_write,"13 ");
916                   ssiWriteProc(d,(procinfov)dd);
917                   break;
918          case LIST_CMD:
919                   fprintf(d->f_write,"14 ");
920                   ssiWriteList(l,(lists)dd);
921                   break;
922          default: Werror("not implemented (t:%d, rtyp:%d)",tt, data->rtyp);
923                   d->level=0;
924                   return TRUE;
925    }
926    if (d->level<=1) { fprintf(d->f_write,"\n"); fflush(d->f_write); }
927    data=data->next;
928  }
929  d->level--;
930  return FALSE;
931}
932
933si_link_extension slInitSsiExtension(si_link_extension s)
934{
935  s->Open=(slOpenProc)ssiOpen;
936  s->Close=(slCloseProc)ssiClose;
937  s->Kill=(slKillProc)ssiClose;
938  s->Read=(slReadProc)ssiRead1;
939  s->Read2=(slRead2Proc)NULL;
940  s->Write=(slWriteProc)ssiWrite;
941
942  s->Status=slStatusSsi;
943  s->type="ssi";
944  return s;
945}
946const char* slStatusSsi(si_link l, const char* request)
947{
948  ssiInfo *d=(ssiInfo*)l->data;
949  if (d==NULL) return "not open";
950  if (((strcmp(l->mode,"fork")==0)||(strcmp(l->mode,"tcp")==0))
951  && (strcmp(request, "read") == 0))
952  {
953    fd_set  mask, fdmask;
954    struct timeval wt;
955    loop
956    {
957      /* Don't block. Return socket status immediately. */
958      wt.tv_sec  = 0;
959      wt.tv_usec = 0;
960
961      FD_ZERO(&mask);
962      FD_SET(d->fd_read, &mask);
963      //Print("test fd %d\n",d->fd_read);
964    /* check with select: chars waiting: no -> not ready */
965      switch (select(d->fd_read+1, &mask, NULL, NULL, &wt))
966      {
967        case 0: /* not ready */ return "not ready";
968        case -1: /*error*/      return "error";
969        case 1: /*ready ? */    break;
970      }
971    /* yes: read 1 char*/
972    /* if \n, check again with select else ungetc(c), ready*/
973      int c=fgetc(d->f_read);
974      //Print("try c=%d\n",c);
975      if (c== -1) return "eof";
976      else if (isdigit(c))
977      { ungetc(c,d->f_read); return "ready"; }
978      else if (c>' ')
979      {
980        Werror("unknown char in ssiLink(%d)",c);
981        return "error";
982      }
983      /* else: next char */
984    }
985  }
986  else if (strcmp(request, "read") == 0)
987  {
988    if (SI_LINK_R_OPEN_P(l) && (!feof(d->f_read))) return "ready";
989    else return "not ready";
990  }
991  else if (strcmp(request, "write") == 0)
992  {
993    if (SI_LINK_W_OPEN_P(l)) return "ready";
994    else return "not ready";
995  }
996  else return "unknown status request";
997}
998
999int ssiBatch(const char *host, const char * port)
1000/* return 0 on success, >0 else*/
1001{
1002  si_link l=(si_link) omAlloc0Bin(sip_link_bin);
1003  char *buf=(char*)omAlloc(256);
1004  sprintf(buf,"ssi:connect %s:%s",host,port);
1005  slInit(l, buf);
1006  slOpen(l,SI_LINK_OPEN);
1007  SI_LINK_SET_RW_OPEN_P(l);
1008  loop
1009  {
1010    leftv h=ssiRead1(l); /*contains an exit.... */
1011    if (feErrors != NULL && *feErrors != '\0')
1012    {
1013      // handle errors:
1014      PrintS(feErrors); /* currently quite simple */
1015      *feErrors = '\0';
1016    }
1017    ssiWrite(l,h);
1018    h->CleanUp();
1019    omFreeBin(h, sleftv_bin);
1020  }
1021  /* never reached*/
1022  exit(0);
1023}
1024
1025// ----------------------------------------------------------------
1026// format
1027// 1 int %d
1028// 2 string <len> %s
1029// 3 number
1030// 4 bigint 4 %d or 3 <mpz_t>
1031// 5 ring
1032// 6 poly
1033// 7 ideal
1034// 8 matrix
1035// 9 vector
1036// 10 module
1037// 11 command
1038// 12 def <len> %s
1039// 13 proc <len> %s
1040// 14 list %d <elem1> ....
1041// 15 setring .......
1042// 16 nothing
1043//
1044// 99: quit Singular
Note: See TracBrowser for help on using the repository browser.