source: git/Singular/ssiLink.cc @ f5a3a23

spielwiese
Last change on this file since f5a3a23 was 0fb34ba, checked in by Hans Schoenemann <hannes@…>, 13 years ago
start with dir Singular: fix includes and makefile
  • Property mode set to 100644
File size: 40.7 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 <netdb.h>
21#include <sys/wait.h>
22
23
24
25#include <kernel/mod2.h>
26// #include "mod2.h"
27
28#include <Singular/tok.h>
29#include <Singular/ipid.h>
30#include <omalloc/omalloc.h>
31#include <polys/monomials/ring.h>
32#include <polys/matpol.h>
33#include <kernel/ideals.h>
34#include <polys/polys.h>
35#include <kernel/longrat.h>
36#include <kernel/ideals.h>
37#include <misc/intvec.h>
38#include <misc/options.h>
39#include <kernel/timer.h>
40#include <Singular/subexpr.h>
41#include <Singular/silink.h>
42#include <Singular/cntrlc.h>
43#include <Singular/lists.h>
44#include <Singular/blackbox.h>
45#include <Singular/ssiLink.h>
46
47#ifdef HAVE_MPSR
48#include <Singular/mpsr.h>
49#endif
50
51//#if (_POSIX_C_SOURCE >= 200112L) || (_XOPEN_SOURCE >= 600)
52//#define HAVE_PSELECT
53//#endif
54
55#define SSI_VERSION 2
56
57typedef struct
58{
59  FILE *f_read;
60  FILE *f_write;
61  ring r;
62  pid_t pid; /* only valid for fork/tcp mode*/
63  int fd_read,fd_write; /* only valid for fork/tcp mode*/
64  char level;
65  char ungetc_buf; /* status sets to !=0, if ungetc was used, ssiRead* set to 0*/
66  char send_quit_at_exit;
67
68} ssiInfo;
69
70
71link_list ssiToBeClosed=NULL;
72sigset_t ssi_sigmask;
73sigset_t ssi_oldmask;
74#define SSI_BLOCK_CHLD sigprocmask(SIG_SETMASK, &ssi_sigmask, &ssi_oldmask)
75#define SSI_UNBLOCK_CHLD sigprocmask(SIG_SETMASK, &ssi_oldmask, NULL)
76
77// the helper functions:
78void ssiSetCurrRing(const ring r)
79{
80  if (!rEqual(r,currRing,1))
81  {
82    char name[20];
83    int nr=0;
84    do
85    { sprintf(name,"ssiRing%d",nr); nr++; }
86    while(IDROOT->get(name, 0)!=NULL);
87    idhdl h=enterid(omStrDup(name),0,RING_CMD,&IDROOT,FALSE);
88    IDRING(h)=r;
89    r->ref++;
90    rSetHdl(h);
91  }
92}
93// the implementation of the functions:
94void ssiWriteInt(ssiInfo *d,const int i)
95{
96  SSI_BLOCK_CHLD;
97  fprintf(d->f_write,"%d ",i);
98  SSI_UNBLOCK_CHLD;
99  //if (d->f_debug!=NULL) fprintf(d->f_debug,"int: %d ",i);
100}
101
102void ssiWriteString(ssiInfo *d,const char *s)
103{
104  SSI_BLOCK_CHLD;
105  fprintf(d->f_write,"%d %s ",strlen(s),s);
106  SSI_UNBLOCK_CHLD;
107  //if (d->f_debug!=NULL) fprintf(d->f_debug,"stringi: %d \"%s\" ",strlen(s),s);
108}
109
110
111void ssiWriteBigInt(const ssiInfo *d, const number n)
112{
113  // syntax is as follows:
114  // case 2 Q:     3 4 <int>
115  //        or     3 3 <mpz_t nominator>
116  SSI_BLOCK_CHLD;
117  if(SR_HDL(n) & SR_INT)
118  {
119    fprintf(d->f_write,"4 %ld ",SR_TO_INT(n));
120    //if (d->f_debug!=NULL) fprintf(d->f_debug,"bigint: short \"%ld\" ",SR_TO_INT(n));
121  }
122  else if (n->s==3)
123  {
124    fputs("3 ",d->f_write);
125    mpz_out_str(d->f_write,10,n->z);
126    fputc(' ',d->f_write);
127    //gmp_fprintf(d->f_write,"3 %Zd ",n->z);
128    //if (d->f_debug!=NULL) gmp_fprintf(d->f_debug,"bigint: gmp \"%Zd\" ",n->z);
129  }
130  else WerrorS("illiegal bigint");
131  SSI_UNBLOCK_CHLD;
132}
133
134void ssiWriteNumber(const ssiInfo *d, const number n)
135{
136  // syntax is as follows:
137  // case 1 Z/p:   3 <int>
138  // case 2 Q:     3 4 <int>
139  //        or     3 0 <mpz_t nominator> <mpz_t denominator>
140  //        or     3 1  dto.
141  //        or     3 3 <mpz_t nominator>
142  //        or     3 5 <mpz_t raw nom.> <mpz_t raw denom.>
143  //        or     3 6 <mpz_t raw nom.> <mpz_t raw denom.>
144  //        or     3 7 <mpz_t raw nom.>
145  SSI_BLOCK_CHLD;
146  if(rField_is_Zp(d->r))
147  {
148    fprintf(d->f_write,"%d ",(int)(long)n);
149    //if (d->f_debug!=NULL) fprintf(d->f_debug,"number: \"%ld\" ",(int)(long)n);
150  }
151  else if (rField_is_Q(d->r))
152  {
153    if(SR_HDL(n) & SR_INT)
154    {
155      fprintf(d->f_write,"4 %ld ",SR_TO_INT(n));
156      //if (d->f_debug!=NULL) fprintf(d->f_debug,"number: short \"%ld\" ",SR_TO_INT(n));
157    }
158    else if (n->s<2)
159    {
160      //gmp_fprintf(d->f_write,"%d %Zd %Zd ",n->s,n->z,n->n);
161      fprintf(d->f_write,"%d ",n->s+5);
162      mpz_out_str (d->f_write,32, n->z);
163      fputc(' ',d->f_write);
164      mpz_out_str (d->f_write,32, n->n);
165      fputc(' ',d->f_write);
166
167      //if (d->f_debug!=NULL) gmp_fprintf(d->f_debug,"number: s=%d gmp/gmp \"%Zd %Zd\" ",n->s,n->z,n->n);
168    }
169    else /*n->s==3*/
170    {
171      //gmp_fprintf(d->f_write,"3 %Zd ",n->z);
172      fputs("8 ",d->f_write);
173      mpz_out_str (d->f_write,32, n->z);
174      fputc(' ',d->f_write);
175
176      //if (d->f_debug!=NULL) gmp_fprintf(d->f_debug,"number: gmp \"%Zd\" ",n->z);
177    }
178  }
179  else WerrorS("coeff field not implemented");
180  SSI_UNBLOCK_CHLD;
181}
182
183void ssiWriteRing(ssiInfo *d,const ring r)
184{
185  /* 5 <ch> <N> <l1> <v1> ...<lN> <vN> <number of orderings> <ord1> <block0_1> <block1_1> .... */
186  if (d->r!=NULL) rKill(d->r);
187  d->r=r;
188  d->r->ref++;
189  SSI_BLOCK_CHLD;
190  fprintf(d->f_write,"%d %d ",r->ch,r->N);
191
192  int i;
193  for(i=0;i<r->N;i++)
194  {
195    fprintf(d->f_write,"%d %s ",strlen(r->names[i]),r->names[i]);
196  }
197  /* number of orderings:*/
198  i=0;
199  while (r->order[i]!=0) i++;
200  fprintf(d->f_write,"%d ",i);
201  /* each ordering block: */
202  i=0;
203  while(r->order[i]!=0)
204  {
205    fprintf(d->f_write,"%d %d %d ",r->order[i],r->block0[i], r->block1[i]);
206    i++;
207  }
208  SSI_UNBLOCK_CHLD;
209}
210
211void ssiWritePoly(ssiInfo *d, int typ, poly p)
212{
213  SSI_BLOCK_CHLD;
214  fprintf(d->f_write,"%d ",pLength(p));//number of terms
215  SSI_UNBLOCK_CHLD;
216  int i;
217
218  while(p!=NULL)
219  {
220    ssiWriteNumber(d,pGetCoeff(p));
221    //nWrite(fich,pGetCoeff(p));
222    SSI_BLOCK_CHLD;
223    fprintf(d->f_write,"%ld ",p_GetComp(p,d->r));//component
224
225    for(int j=1;j<=rVar(d->r);j++)
226    {
227      fprintf(d->f_write,"%ld ",p_GetExp(p,j,d->r ));//x^j
228    }
229    pIter(p);
230    SSI_UNBLOCK_CHLD;
231  }
232}
233
234void ssiWriteIdeal(ssiInfo *d, int typ,ideal I)
235{
236   // syntax: 7 # of elements <poly 1> <poly2>.....
237   // syntax: 8 <rows> <cols> <poly 1> <poly2>.....
238   matrix M=(matrix)I;
239   SSI_BLOCK_CHLD;
240   if (typ==MATRIX_CMD)
241        fprintf(d->f_write,"%d %d ", MATROWS(M),MATCOLS(M));
242   else
243     fprintf(d->f_write,"%d ",IDELEMS(I));
244    SSI_UNBLOCK_CHLD;
245
246   int i;
247   int tt;
248   if (typ==MODUL_CMD) tt=VECTOR_CMD;
249   else                tt=POLY_CMD;
250
251   for(i=0;i<IDELEMS(I);i++)
252   {
253     ssiWritePoly(d,tt,I->m[i]);
254   }
255}
256void ssiWriteCommand(si_link l, command D)
257{
258  ssiInfo *d=(ssiInfo*)l->data;
259  // syntax: <num ops> <operation> <op1> <op2> ....
260  SSI_BLOCK_CHLD;
261  fprintf(d->f_write,"%d %d ",D->argc,D->op);
262  SSI_UNBLOCK_CHLD;
263  if (D->argc >0) ssiWrite(l, &(D->arg1));
264  if (D->argc < 4)
265  {
266    if (D->argc >1) ssiWrite(l, &(D->arg2));
267    if (D->argc >2) ssiWrite(l, &(D->arg3));
268  }
269}
270
271void ssiWriteProc(ssiInfo *d,procinfov p)
272{
273  ssiWriteString(d,p->data.s.body);
274}
275
276void ssiWriteList(si_link l,lists dd)
277{
278  ssiInfo *d=(ssiInfo*)l->data;
279  SSI_BLOCK_CHLD;
280  fprintf(d->f_write,"%d ",dd->nr+1);
281  SSI_UNBLOCK_CHLD;
282  int i;
283  for(i=0;i<=dd->nr;i++)
284  {
285    ssiWrite(l,&(dd->m[i]));
286  }
287}
288void ssiWriteIntvec(ssiInfo *d,intvec * v)
289{
290  SSI_BLOCK_CHLD;
291  fprintf(d->f_write,"%d ",v->length());
292  int i;
293  for(i=0;i<v->length();i++)
294  {
295    fprintf(d->f_write,"%d ",(*v)[i]);
296  }
297  SSI_UNBLOCK_CHLD;
298}
299
300char *ssiReadString(ssiInfo *d)
301{
302  char *buf;
303  int l;
304  SSI_BLOCK_CHLD;
305  fscanf(d->f_read,"%d ",&l);
306  SSI_UNBLOCK_CHLD;
307  buf=(char*)omAlloc(l+1);
308  SSI_BLOCK_CHLD;
309  fread(buf,1,l,d->f_read);
310  SSI_UNBLOCK_CHLD;
311  buf[l]='\0';
312  return buf;
313}
314
315int ssiReadInt(FILE *fich)
316{
317  int d;
318  SSI_BLOCK_CHLD;
319  fscanf(fich,"%d",&d);
320  SSI_UNBLOCK_CHLD;
321  return d;
322}
323
324number ssiReadBigInt(ssiInfo *d)
325{
326  int sub_type=-1;
327  SSI_BLOCK_CHLD;
328  fscanf(d->f_read,"%d",&sub_type);
329  SSI_UNBLOCK_CHLD;
330  switch(sub_type)
331  {
332   case 3:
333     {// read int or mpz_t or mpz_t, mpz_t
334       number n=nlRInit(0);
335       SSI_BLOCK_CHLD;
336       mpz_inp_str(n->z,d->f_read,0);
337       SSI_UNBLOCK_CHLD;
338       n->s=sub_type;
339       return n;
340     }
341   case 4:
342     {
343       int dd;
344       SSI_BLOCK_CHLD;
345       fscanf(d->f_read,"%d",&dd);
346       SSI_UNBLOCK_CHLD;
347       return INT_TO_SR(dd);
348     }
349   default:
350       Werror("error in reading bigint: invalid subtype %d",sub_type);
351       return NULL;
352   }
353}
354
355number ssiReadNumber(ssiInfo *d)
356{
357  if (rField_is_Q(d->r))
358  {
359     int sub_type=-1;
360     SSI_BLOCK_CHLD;
361     fscanf(d->f_read,"%d",&sub_type);
362     SSI_UNBLOCK_CHLD;
363     switch(sub_type)
364     {
365     case 0:
366     case 1:
367       {// read mpz_t, mpz_t
368         number n=nlRInit(0);
369         mpz_init(n->n);
370         SSI_BLOCK_CHLD;
371         gmp_fscanf(d->f_read,"%Zd %Zd",n->z,n->n);
372         SSI_UNBLOCK_CHLD;
373         n->s=sub_type;
374         return n;
375       }
376
377     case 3:
378       {// read mpz_t
379         number n=nlRInit(0);
380         SSI_BLOCK_CHLD;
381         gmp_fscanf(d->f_read,"%Zd",n->z);
382         SSI_UNBLOCK_CHLD;
383         n->s=3; /*sub_type*/
384         return n;
385       }
386     case 4:
387       {
388         int dd;
389         SSI_BLOCK_CHLD;
390         fscanf(d->f_read,"%d",&dd);
391         SSI_UNBLOCK_CHLD;
392         return INT_TO_SR(dd);
393       }
394     case 5:
395     case 6:
396       {// read raw mpz_t, mpz_t
397         number n=nlRInit(0);
398         mpz_init(n->n);
399         SSI_BLOCK_CHLD;
400         mpz_inp_str (n->z, d->f_read, 32);
401         mpz_inp_str (n->n, d->f_read, 32);
402         SSI_UNBLOCK_CHLD;
403         n->s=sub_type-5;
404         return n;
405       }
406     case 8:
407       {// read raw mpz_t
408         number n=nlRInit(0);
409         SSI_BLOCK_CHLD;
410         mpz_inp_str (n->z, d->f_read, 32);
411         SSI_UNBLOCK_CHLD;
412         n->s=sub_type=3; /*subtype-5*/
413         return n;
414       }
415
416     default: Werror("error in reading number: invalid subtype %d",sub_type);
417              return NULL;
418     }
419  }
420  else if (rField_is_Zp(d->r))
421  {
422    // read int
423    int dd;
424    SSI_BLOCK_CHLD;
425    fscanf(d->f_read,"%d",&dd);
426    SSI_UNBLOCK_CHLD;
427    return (number)dd;
428  }
429  else Werror("coeffs not implemented");
430  return NULL;
431}
432
433ring ssiReadRing(ssiInfo *d)
434{
435/* syntax is <ch> <N> <l1> <v1> ...<lN> <vN> <number of orderings> <ord1> <block0_1> <block1_1> .... */
436  int ch, N,i,l;
437  char **names;
438  SSI_BLOCK_CHLD;
439  fscanf(d->f_read,"%d %d ",&ch,&N);
440  SSI_UNBLOCK_CHLD;
441  names=(char**)omAlloc(N*sizeof(char*));
442  for(i=0;i<N;i++)
443  {
444    names[i]=ssiReadString(d);
445  }
446  // read the orderings:
447  int num_ord; // number of orderings
448  SSI_BLOCK_CHLD;
449  fscanf(d->f_read,"%d",&num_ord);
450  SSI_UNBLOCK_CHLD;
451  int *ord=(int *)omAlloc0((num_ord+1)*sizeof(int));
452  int *block0=(int *)omAlloc0((num_ord+1)*sizeof(int));
453  int *block1=(int *)omAlloc0((num_ord+1)*sizeof(int));
454  SSI_BLOCK_CHLD;
455  for(i=0;i<num_ord;i++)
456  {
457     fscanf(d->f_read,"%d %d %d",&ord[i],&block0[i],&block1[i]);
458  }
459  SSI_UNBLOCK_CHLD;
460  return rDefault(ch,N,names,num_ord,ord,block0,block1);
461}
462
463poly ssiReadPoly(ssiInfo *D)
464{
465// < # of terms> < term1> < .....
466  int n,i,l;
467  n=ssiReadInt(D->f_read);
468  //Print("poly: terms:%d\n",n);
469  poly p;
470  int j;
471  j=0;
472  poly ret=NULL;
473  poly prev=NULL;
474  for(l=0;l<n;l++) // read n terms
475  {
476// coef,comp.exp1,..exp N
477    p=p_Init(D->r);
478    pGetCoeff(p)=ssiReadNumber(D);
479    int d;
480    SSI_BLOCK_CHLD;
481    fscanf(D->f_read,"%d",&d);
482    SSI_UNBLOCK_CHLD;
483    p_SetComp(p,d,D->r);
484    for(i=1;i<=rVar(D->r);i++)
485    {
486      SSI_BLOCK_CHLD;
487      fscanf(D->f_read,"%d",&d);
488      SSI_UNBLOCK_CHLD;
489      p_SetExp(p,i,d,D->r);
490    }
491    p_Setm(p,D->r);
492    p_Test(p,D->r);
493    if (ret==NULL) ret=p;
494    else           pNext(prev)=p;
495    prev=p;
496 }
497 return ret;
498}
499
500ideal ssiReadIdeal(ssiInfo *d)
501{
502  int n,i;
503  ideal I;
504  SSI_BLOCK_CHLD;
505  fscanf(d->f_read,"%d",&n);
506  SSI_UNBLOCK_CHLD;
507  I=idInit(n,1);
508  for(i=0;i<IDELEMS(I);i++) // read n terms
509  {
510    I->m [i]=ssiReadPoly(d);
511  }
512  return I;
513}
514
515matrix ssiReadMatrix(ssiInfo *d)
516{
517  int n,m,i,j;
518  SSI_BLOCK_CHLD;
519  fscanf(d->f_read,"%d %d",&m,&n);
520  SSI_UNBLOCK_CHLD;
521  matrix M=mpNew(m,n);
522  poly p;
523  for(int i=1;i<=MATROWS(M);i++)
524    for(int j=1;j<=MATCOLS(M);j++)
525    {
526      p=ssiReadPoly(d);
527      MATELEM(M,i,j)=p;
528    }
529  return M;
530}
531
532command ssiReadCommand(si_link l)
533{
534  ssiInfo *d=(ssiInfo*)l->data;
535  // syntax: <num ops> <operation> <op1> <op2> ....
536  command D=(command)omAlloc0(sizeof(*D));
537  int argc,op;
538  SSI_BLOCK_CHLD;
539  fscanf(d->f_read,"%d %d",&argc,&op);
540  SSI_UNBLOCK_CHLD;
541  D->argc=argc; D->op=op;
542  leftv v;
543  if (argc >0)
544  {
545    v=ssiRead1(l);
546    memcpy(&(D->arg1),v,sizeof(*v));
547    omFreeBin(v,sleftv_bin);
548  }
549  if (argc <4)
550  {
551    if (D->argc >1)
552    {
553      v=ssiRead1(l);
554      memcpy(&(D->arg2),v,sizeof(*v));
555      omFreeBin(v,sleftv_bin);
556    }
557    if (D->argc >2)
558    {
559      v=ssiRead1(l);
560      memcpy(&(D->arg3),v,sizeof(*v));
561      omFreeBin(v,sleftv_bin);
562    }
563  }
564  else
565  {
566    leftv prev=&(D->arg1);
567    argc--;
568    while(argc >0)
569    {
570      v=ssiRead1(l);
571      prev->next=v;
572      prev=v;
573      argc--;
574    }
575  }
576  return D;
577}
578
579procinfov ssiReadProc(ssiInfo *d)
580{
581  char *s=ssiReadString(d);
582  procinfov p=(procinfov)omAlloc0Bin(procinfo_bin);
583  p->language=LANG_SINGULAR;
584  p->libname=omStrDup("");
585  p->procname=omStrDup("");
586  p->data.s.body=s;
587  return p;
588}
589lists ssiReadList(si_link l)
590{
591  ssiInfo *d=(ssiInfo*)l->data;
592  int nr;
593  SSI_BLOCK_CHLD;
594  fscanf(d->f_read,"%d",&nr);
595  SSI_UNBLOCK_CHLD;
596  lists L=(lists)omAlloc(sizeof(*L));
597  L->Init(nr);
598
599  int i;
600  leftv v;
601  for(i=0;i<nr;i++)
602  {
603    v=ssiRead1(l);
604    memcpy(&(L->m[i]),v,sizeof(*v));
605    omFreeBin(v,sleftv_bin);
606  }
607  return L;
608}
609intvec* ssiReadIntvec(ssiInfo *d)
610{
611  int nr;
612  SSI_BLOCK_CHLD;
613  fscanf(d->f_read,"%d",&nr);
614  SSI_UNBLOCK_CHLD;
615  intvec *v=new intvec(nr);
616  SSI_BLOCK_CHLD;
617  for(int i=0;i<nr;i++)
618  {
619    fscanf(d->f_read,"%d",&((*v)[i]));
620  }
621  SSI_UNBLOCK_CHLD;
622  return v;
623}
624
625void ssiReadBlackbox(leftv res, si_link l)
626{
627  ssiInfo *d=(ssiInfo*)l->data;
628  int throwaway;
629  SSI_BLOCK_CHLD;
630  fscanf(d->f_read,"%d ",&throwaway);
631  SSI_UNBLOCK_CHLD;
632  char *name=ssiReadString(d);
633  int tok;
634  blackboxIsCmd(name,tok);
635  if (tok>MAX_TOK)
636  {
637    blackbox *b=getBlackboxStuff(tok);
638    res->rtyp=tok;
639    b->blackbox_deserialize(&b,&(res->data),l);
640  }
641  else
642  {
643    Werror("blackbox %s not found",name);
644  }
645}
646
647//**************************************************************************/
648
649BOOLEAN ssiOpen(si_link l, short flag, leftv u)
650{
651  if (l!=NULL)
652  {
653    const char *mode;
654    ssiInfo *d=(ssiInfo*)omAlloc0(sizeof(ssiInfo));
655    sigprocmask(SIG_SETMASK, NULL, &ssi_sigmask);
656    sigaddset(&ssi_sigmask, SIGCHLD);
657    if (flag & SI_LINK_OPEN)
658    {
659      if (l->mode[0] != '\0' && (strcmp(l->mode, "r") == 0))
660        flag = SI_LINK_READ;
661      else flag = SI_LINK_WRITE;
662    }
663
664    if (flag == SI_LINK_READ) mode = "r";
665    else if (strcmp(l->mode, "w") == 0) mode = "w";
666    else if (strcmp(l->mode, "fork") == 0) mode = "fork";
667    else if (strcmp(l->mode, "tcp") == 0) mode = "tcp";
668    else if (strcmp(l->mode, "connect") == 0) mode = "connect";
669    else mode = "a";
670
671
672    SI_LINK_SET_OPEN_P(l, flag);
673    l->data=d;
674    omFree(l->mode);
675    l->mode = omStrDup(mode);
676
677    if (l->name[0] == '\0')
678    {
679      if (strcmp(mode,"fork")==0)
680      {
681        link_list n=(link_list)omAlloc(sizeof(link_struct));
682        n->u=u;
683        n->l=l;
684        n->next=(void *)ssiToBeClosed;
685        ssiToBeClosed=n;
686
687        int pc[2];
688        int cp[2];
689        pipe(pc);
690        pipe(cp);
691        pid_t pid=fork();
692        if (pid==0) /*fork: child*/
693        {
694          link_list hh=(link_list)ssiToBeClosed->next;
695          /* we know: l is the first entry in ssiToBeClosed-list */
696          while(hh!=NULL)
697          {
698            ssiInfo *dd=(ssiInfo*)hh->l->data;
699            fclose(dd->f_read);
700            fclose(dd->f_write);
701            if (dd->r!=NULL) rKill(dd->r);
702            omFreeSize((ADDRESS)dd,(sizeof *dd));
703            hh->l->data=NULL;
704            SI_LINK_SET_CLOSE_P(hh->l);
705            link_list nn=(link_list)hh->next;
706            omFree(hh);
707            hh=nn;
708          }
709          ssiToBeClosed->next=NULL;
710          close(pc[1]); close(cp[0]);
711          d->f_read=fdopen(pc[0],"r");
712          d->fd_read=pc[0];
713          d->f_write=fdopen(cp[1],"w");
714          d->fd_write=cp[1];
715          l->data=d;
716          omFree(l->mode);
717          l->mode = omStrDup(mode);
718          singular_in_batchmode=TRUE;
719          SI_LINK_SET_RW_OPEN_P(l);
720          //myynest=0;
721          fe_fgets_stdin=fe_fgets_dummy;
722          if ((u!=NULL)&&(u->rtyp==IDHDL))
723          {
724            idhdl h=(idhdl)u->data;
725            h->lev=0;
726          }
727          loop
728          {
729            leftv h=ssiRead1(l); /*contains an exit.... */
730            if (feErrors != NULL && *feErrors != '\0')
731            {
732              // handle errors:
733              PrintS(feErrors); /* currently quite simple */
734              *feErrors = '\0';
735            }
736            ssiWrite(l,h);
737            h->CleanUp();
738            omFreeBin(h, sleftv_bin);
739          }
740          /* never reached*/
741        }
742        else if (pid>0) /*fork: parent*/
743        {
744          d->pid=pid;
745          close(pc[0]); close(cp[1]);
746          d->f_read=fdopen(cp[0],"r");
747          d->fd_read=cp[0];
748          d->f_write=fdopen(pc[1],"w");
749          d->fd_write=pc[1];
750          SI_LINK_SET_RW_OPEN_P(l);
751          d->send_quit_at_exit=1;
752        }
753        else
754        {
755          Werror("fork failed (%d)",errno);
756          l->data=NULL;
757          omFree(d);
758          return TRUE;
759        }
760      }
761      // ---------------------------------------------------------------------
762      else if (strcmp(mode,"tcp")==0)
763      {
764        int sockfd, newsockfd, portno, clilen;
765        struct sockaddr_in serv_addr, cli_addr;
766        int n;
767        sockfd = socket(AF_INET, SOCK_STREAM, 0);
768        if(sockfd < 0)
769        {
770          WerrorS("ERROR opening socket");
771          l->data=NULL;
772          omFree(d);
773          return TRUE;
774        }
775        memset((char *) &serv_addr,0, sizeof(serv_addr));
776        portno = 1025;
777        serv_addr.sin_family = AF_INET;
778        serv_addr.sin_addr.s_addr = INADDR_ANY;
779        do
780        {
781          portno++;
782          serv_addr.sin_port = htons(portno);
783          if(portno > 50000)
784          {
785            WerrorS("ERROR on binding (no free port available?)");
786            l->data=NULL;
787            omFree(d);
788            return TRUE;
789          }
790        }
791        while(bind(sockfd, (struct sockaddr *) &serv_addr, sizeof(serv_addr)) < 0);
792        Print("waiting on port %d\n", portno);mflush();
793        listen(sockfd,5);
794        newsockfd = accept(sockfd, (struct sockaddr *) &cli_addr, (socklen_t *)&clilen);
795        if(newsockfd < 0)
796        {
797          WerrorS("ERROR on accept");
798          l->data=NULL;
799          omFree(d);
800          return TRUE;
801        }
802        PrintS("client accepted\n");
803        d->fd_read = newsockfd;
804        d->fd_write = newsockfd;
805        d->f_read = fdopen(newsockfd, "r");
806        d->f_write = fdopen(newsockfd, "w");
807        SI_LINK_SET_RW_OPEN_P(l);
808        close(sockfd);
809      }
810      // no ssi-Link on stdin or stdout
811      else
812      {
813        Werror("invalid mode >>%s<< for ssi",mode);
814        l->data=NULL;
815        omFree(d);
816        return TRUE;
817      }
818    }
819    // =========================================================================
820    else /*l->name=NULL*/
821    {
822      // tcp mode
823      if(strcmp(mode,"tcp")==0)
824      {
825        int sockfd, newsockfd, portno, clilen;
826        struct sockaddr_in serv_addr, cli_addr;
827        int n;
828        sockfd = socket(AF_INET, SOCK_STREAM, 0);
829        if(sockfd < 0)
830        {
831          WerrorS("ERROR opening socket");
832          l->data=NULL;
833          omFree(d);
834          return TRUE;
835        }
836        memset((char *) &serv_addr,0, sizeof(serv_addr));
837        portno = 1025;
838        serv_addr.sin_family = AF_INET;
839        serv_addr.sin_addr.s_addr = INADDR_ANY;
840        do
841        {
842          portno++;
843          serv_addr.sin_port = htons(portno);
844          if(portno > 50000)
845          {
846            WerrorS("ERROR on binding (no free port available?)");
847            l->data=NULL;
848            return TRUE;
849          }
850        }
851        while(bind(sockfd, (struct sockaddr *) &serv_addr, sizeof(serv_addr)) < 0);
852        //Print("waiting on port %d\n", portno);mflush();
853        listen(sockfd,5);
854        char* cli_host = (char*)omAlloc(256);
855        char* path = (char*)omAlloc(1024);
856        int r = sscanf(l->name,"%255[^:]:%s",cli_host,path);
857        if(r == 0)
858        {
859          WerrorS("ERROR: no host specified");
860          l->data=NULL;
861          omFree(d);
862          omFree(path);
863          omFree(cli_host);
864          return TRUE;
865        }
866        else if(r == 1)
867        {
868          WarnS("program not specified, using /usr/local/bin/Singular");
869          strcpy(path,"/usr/local/bin/Singular");
870        }
871        char* ssh_command = (char*)omAlloc(256);
872        char* ser_host = (char*)omAlloc(64);
873        gethostname(ser_host,64);
874        sprintf(ssh_command,"ssh %s %s -q --batch --link=ssi --MPhost=%s --MPport=%d &",cli_host,path,ser_host,portno);
875        //Print("client on %s started:%s\n",cli_host,path);
876        omFree(path);
877        omFree(cli_host);
878        if (TEST_OPT_PROT) { Print("running >>%s<<\n",ssh_command); }
879        system(ssh_command);
880        omFree(ssh_command);
881        omFree(ser_host);
882        clilen = sizeof(cli_addr);
883        newsockfd = accept(sockfd, (struct sockaddr *) &cli_addr, (socklen_t *)&clilen);
884        if(newsockfd < 0)
885        {
886          WerrorS("ERROR on accept");
887          l->data=NULL;
888          omFree(d);
889          return TRUE;
890        }
891        //PrintS("client accepted\n");
892        d->fd_read = newsockfd;
893        d->fd_write = newsockfd;
894        d->f_read = fdopen(newsockfd, "r");
895        d->f_write = fdopen(newsockfd, "w");
896        SI_LINK_SET_RW_OPEN_P(l);
897        d->send_quit_at_exit=1;
898        close(sockfd);
899        fprintf(d->f_write,"98 %d %d %u %u\n",SSI_VERSION,MAX_TOK,test,verbose);
900      }
901      // ----------------------------------------------------------------------
902      else if(strcmp(mode,"connect")==0)
903      {
904        char* host = (char*)omAlloc(256);
905        int sockfd, portno, n;
906        struct sockaddr_in serv_addr;
907        struct hostent *server;
908 
909        sscanf(l->name,"%255[^:]:%d",host,&portno);
910        //Print("connect to host %s, port %d\n",host,portno);mflush();
911        if (portno!=0)
912        {
913          sockfd = socket(AF_INET, SOCK_STREAM, 0);
914          if (sockfd < 0) { WerrorS("ERROR opening socket"); return TRUE; }
915          server = gethostbyname(host);
916          if (server == NULL) {  WerrorS("ERROR, no such host");  return TRUE; }
917          memset((char *) &serv_addr, 0, sizeof(serv_addr));
918          serv_addr.sin_family = AF_INET;
919          memcpy((char *)&serv_addr.sin_addr.s_addr,
920                (char *)server->h_addr,
921                server->h_length);
922          serv_addr.sin_port = htons(portno);
923          if (connect(sockfd,(sockaddr*)&serv_addr,sizeof(serv_addr)) < 0)
924          { Werror("ERROR connecting(errno=%d)",errno); return TRUE; }
925          //PrintS("connected\n");mflush();
926          d->f_read=fdopen(sockfd,"r");
927          d->fd_read=sockfd;
928          d->f_write=fdopen(sockfd,"w");
929          d->fd_write=sockfd;
930          SI_LINK_SET_RW_OPEN_P(l);
931          d->send_quit_at_exit=1;
932          omFree(host);
933        }
934        else
935        {
936          l->data=NULL;
937          omFree(d);
938          return TRUE;
939        }
940      }
941      // ======================================================================
942      else
943      {
944        // normal link to a file
945        FILE *outfile;
946        char *filename=l->name;
947 
948        if(filename[0]=='>')
949        {
950          if (filename[1]=='>')
951          {
952            filename+=2;
953            mode = "a";
954          }
955          else
956          {
957            filename++;
958            mode="w";
959          }
960        }
961        outfile=myfopen(filename,mode);
962        if (outfile!=NULL)
963        {
964          if (strcmp(l->mode,"r")==0) d->f_read = outfile;
965          else
966          {
967            d->f_write = outfile;
968            fprintf(d->f_write,"98 %d %d %u %u\n",SSI_VERSION,MAX_TOK,test,verbose);
969          }
970        }
971        else
972        {
973          omFree(d);
974          l->data=NULL;
975          return TRUE;
976        }
977      }
978    }
979  }
980
981  return FALSE;
982}
983
984//**************************************************************************/
985BOOLEAN ssiClose(si_link l)
986{
987  if (l!=NULL)
988  {
989    ssiInfo *d = (ssiInfo *)l->data;
990    if (d!=NULL)
991    {
992      if (d->r!=NULL) rKill(d->r);
993      if ((strcmp(l->mode,"tcp")==0)
994      || (strcmp(l->mode,"fork")==0))
995      {
996        link_list hh=ssiToBeClosed;
997        if (hh!=NULL)
998        {
999          if (hh->l==l)
1000          {
1001             ssiToBeClosed=(link_list)hh->next;
1002             omFreeSize(hh,sizeof(link_struct));
1003          }
1004          else while(hh!=NULL)
1005          {
1006            link_list hhh=(link_list)hh->next;
1007            if (hhh->l==l)
1008            {
1009              hh->next=hhh->next;
1010              omFreeSize(hhh,sizeof(link_struct));
1011              break;
1012            }
1013            hh=(link_list)hh->next;
1014          }
1015        }
1016        if (d->send_quit_at_exit)
1017        { 
1018          fputs("99\n",d->f_write);fflush(d->f_write);
1019        }
1020      }
1021      if (d->f_read!=NULL) fclose(d->f_read);
1022      if (d->f_write!=NULL) fclose(d->f_write);
1023      if (d->pid!=0)
1024      {
1025        kill(d->pid,15);
1026        waitpid(d->pid,NULL,WNOHANG);
1027        kill(d->pid,9); // just to be sure
1028        waitpid(d->pid,NULL,0);
1029      }
1030      omFreeSize((ADDRESS)d,(sizeof *d));
1031    }
1032    l->data=NULL;
1033    SI_LINK_SET_CLOSE_P(l);
1034  }
1035  return FALSE;
1036}
1037
1038//**************************************************************************/
1039leftv ssiRead1(si_link l)
1040{
1041  ssiInfo *d = (ssiInfo *)l->data;
1042  d->ungetc_buf='\0';
1043  leftv res=(leftv)omAlloc0(sizeof(sleftv));
1044  int t=0;
1045  SSI_BLOCK_CHLD;
1046  fscanf(d->f_read,"%d",&t);
1047  SSI_UNBLOCK_CHLD;
1048  //Print("got type %d\n",t);
1049  switch(t)
1050  {
1051    case 1:res->rtyp=INT_CMD;
1052           res->data=(char *)ssiReadInt(d->f_read);
1053           break;
1054    case 2:res->rtyp=STRING_CMD;
1055           res->data=(char *)ssiReadString(d);
1056           break;
1057    case 3:res->rtyp=NUMBER_CMD;
1058           res->data=(char *)ssiReadNumber(d);
1059           break;
1060    case 4:res->rtyp=BIGINT_CMD;
1061           res->data=(char *)ssiReadBigInt(d);
1062           break;
1063    case 15:
1064    case 5:{
1065             d->r=ssiReadRing(d);
1066             d->r->ref++;
1067             res->rtyp=RING_CMD;
1068             res->data=(char*)d->r;
1069             // we are in the top-level, so set the basering to d->r:
1070             ssiSetCurrRing(d->r);
1071             if (t==15) return ssiRead1(l);
1072           }
1073           break;
1074    case 6:res->rtyp=POLY_CMD;
1075           if (d->r==NULL) goto no_ring;
1076           res->data=(char*)ssiReadPoly(d);
1077           break;
1078    case 7:res->rtyp=IDEAL_CMD;
1079           if (d->r==NULL) goto no_ring;
1080           res->data=(char*)ssiReadIdeal(d);
1081           break;
1082    case 8:res->rtyp=MATRIX_CMD;
1083           if (d->r==NULL) goto no_ring;
1084           res->data=(char*)ssiReadMatrix(d);
1085           break;
1086    case 9:res->rtyp=VECTOR_CMD;
1087           if (d->r==NULL) goto no_ring;
1088           res->data=(char*)ssiReadPoly(d);
1089           break;
1090    case 10:res->rtyp=MODUL_CMD;
1091           if (d->r==NULL) goto no_ring;
1092           res->data=(char*)ssiReadIdeal(d);
1093           break;
1094    case 11:
1095           {
1096             res->rtyp=COMMAND;
1097             res->data=ssiReadCommand(l);
1098             int nok=res->Eval();
1099             if (nok) WerrorS("error in eval");
1100             break;
1101           }
1102    case 12: /*DEF_CMD*/
1103           {
1104             res->rtyp=0;
1105             res->name=(char *)ssiReadString(d);
1106             int nok=res->Eval();
1107             if (nok) WerrorS("error in name lookup");
1108             break;
1109           }
1110    case 13: res->rtyp=PROC_CMD;
1111             res->data=ssiReadProc(d);
1112             break;
1113    case 14: res->rtyp=LIST_CMD;
1114             res->data=ssiReadList(l);
1115             break;
1116    case 16: res->rtyp=NONE; res->data=NULL;
1117             break;
1118    case 17: res->rtyp=INTVEC_CMD;
1119             res->data=ssiReadIntvec(d);
1120             break;
1121    case 20: ssiReadBlackbox(res,l);
1122             break;
1123    // ------------
1124    case 98: // version
1125             {
1126                int n98_v,n98_m;
1127                BITSET n98_o1,n98_o2;
1128                SSI_BLOCK_CHLD;
1129                fscanf(d->f_read,"%d %d %u %u\n",&n98_v,&n98_m,&n98_o1,&n98_o2);
1130                SSI_UNBLOCK_CHLD;
1131                if ((n98_v!=SSI_VERSION) ||(n98_m!=MAX_TOK))
1132                {
1133                  Print("incompatible versions of ssi: %d/%d vs %d/%d",
1134                                  SSI_VERSION,MAX_TOK,n98_v,n98_m);
1135                }
1136                #ifndef NDEBUG
1137                if (TEST_OPT_DEBUG)
1138                  Print("// opening ssi-%d, MAX_TOK=%d\n",n98_v,n98_m);
1139                #endif
1140                test=n98_o1;
1141                verbose=n98_o2;
1142                return ssiRead1(l);
1143             }
1144    case 99: ssiClose(l); exit(0);
1145    case 0: if (feof(d->f_read))
1146            {
1147              ssiClose(l);
1148              res->rtyp=DEF_CMD;
1149              break;
1150            }
1151    default: Werror("not implemented (t:%d)",t);
1152             omFreeSize(res,sizeof(sleftv));
1153             res=NULL;
1154             break;
1155  }
1156  return res;
1157no_ring: WerrorS("no ring");
1158  omFreeSize(res,sizeof(sleftv));
1159  return NULL;
1160}
1161//**************************************************************************/
1162BOOLEAN ssiWrite(si_link l, leftv data)
1163{
1164  if(SI_LINK_W_OPEN_P(l)==0)
1165     if (slOpen(l,SI_LINK_OPEN|SI_LINK_WRITE,NULL)) return TRUE;
1166  ssiInfo *d = (ssiInfo *)l->data;
1167  d->level++;
1168  //FILE *fich=d->f;
1169  while (data!=NULL)
1170  {
1171    int tt=data->Typ();
1172    void *dd=data->Data();
1173
1174    switch(tt /*data->Typ()*/)
1175    {
1176          case NONE/* nothing*/:fputs("16 ",d->f_write);
1177                          break;
1178          case STRING_CMD: fputs("2 ",d->f_write);
1179                           ssiWriteString(d,(char *)dd);
1180                           break;
1181          case INT_CMD: fputs("1 ",d->f_write);
1182                        ssiWriteInt(d,(int)(long)dd);
1183                        break;
1184          case BIGINT_CMD:fputs("4 ",d->f_write);
1185                        ssiWriteBigInt(d,(number)dd);
1186                        break;
1187          case NUMBER_CMD:
1188                          if (d->r!=currRing)
1189                          {
1190                            fputs("15 ",d->f_write);
1191                            ssiWriteRing(d,currRing);
1192                            if (d->level<=1) fputc('\n',d->f_write);
1193                          }
1194                          fputs("3 ",d->f_write);
1195                          ssiWriteNumber(d,(number)dd);
1196                        break;
1197          case RING_CMD:fputs("5 ",d->f_write);
1198                        ssiWriteRing(d,(ring)dd);
1199                        break;
1200          case POLY_CMD:
1201          case VECTOR_CMD:
1202                        if (d->r!=currRing)
1203                        {
1204                          fputs("15 ",d->f_write);
1205                          ssiWriteRing(d,currRing);
1206                          if (d->level<=1) fputc('\n',d->f_write);
1207                        }
1208                        if(tt==POLY_CMD) fputs("6 ",d->f_write);
1209                        else             fputs("9 ",d->f_write);
1210                        ssiWritePoly(d,tt,(poly)dd);
1211                        break;
1212          case IDEAL_CMD:
1213          case MODUL_CMD:
1214          case MATRIX_CMD:
1215                        if (d->r!=currRing)
1216                        {
1217                          fputs("15 ",d->f_write);
1218                          ssiWriteRing(d,currRing);
1219                          if (d->level<=1) fputc('\n',d->f_write);
1220                        }
1221                        if(tt==IDEAL_CMD)       fputs("7 ",d->f_write);
1222                        else if(tt==MATRIX_CMD) fputs("8 ",d->f_write);
1223                        else                    fputs("10 ",d->f_write);
1224                        ssiWriteIdeal(d,tt,(ideal)dd);
1225                        break;
1226          case COMMAND:
1227                   fputs("11 ",d->f_write);
1228                   ssiWriteCommand(l,(command)dd);
1229                   break;
1230          case DEF_CMD: /* not evaluated stuff in quotes */
1231                   fputs("12 ",d->f_write);
1232                   ssiWriteString(d,data->Name());
1233                   break;
1234          case PROC_CMD:
1235                   fputs("13 ",d->f_write);
1236                   ssiWriteProc(d,(procinfov)dd);
1237                   break;
1238          case LIST_CMD:
1239                   fputs("14 ",d->f_write);
1240                   ssiWriteList(l,(lists)dd);
1241                   break;
1242          case INTVEC_CMD:
1243                   fputs("17 ",d->f_write);
1244                   ssiWriteIntvec(d,(intvec *)dd);
1245                   break;
1246          default:
1247            if (tt>MAX_TOK)
1248            {
1249              blackbox *b=getBlackboxStuff(tt);
1250              fputs("20 ",d->f_write);
1251              b->blackbox_serialize(b,dd,l);
1252            }
1253            else
1254            {
1255              Werror("not implemented (t:%d, rtyp:%d)",tt, data->rtyp);
1256              d->level=0;
1257              return TRUE;
1258            }
1259            break;
1260    }
1261    if (d->level<=1) { fputc('\n',d->f_write); fflush(d->f_write); }
1262    data=data->next;
1263  }
1264  d->level--;
1265  return FALSE;
1266}
1267
1268si_link_extension slInitSsiExtension(si_link_extension s)
1269{
1270  s->Open=ssiOpen;
1271  s->Close=ssiClose;
1272  s->Kill=ssiClose;
1273  s->Read=ssiRead1;
1274  s->Read2=(slRead2Proc)NULL;
1275  s->Write=ssiWrite;
1276
1277  s->Status=slStatusSsi;
1278  s->type="ssi";
1279  return s;
1280}
1281
1282const char* slStatusSsi(si_link l, const char* request)
1283{
1284  ssiInfo *d=(ssiInfo*)l->data;
1285  if (d==NULL) return "not open";
1286  if (((strcmp(l->mode,"fork")==0)||(strcmp(l->mode,"tcp")==0))
1287  && (strcmp(request, "read") == 0))
1288  {
1289    fd_set  mask, fdmask;
1290    struct timeval wt;
1291    if (d->ungetc_buf) return "ready";
1292    loop
1293    {
1294      /* Don't block. Return socket status immediately. */
1295      wt.tv_sec  = 0;
1296      wt.tv_usec = 0;
1297
1298      FD_ZERO(&mask);
1299      FD_SET(d->fd_read, &mask);
1300      //Print("test fd %d\n",d->fd_read);
1301    /* check with select: chars waiting: no -> not ready */
1302      switch (select(d->fd_read+1, &mask, NULL, NULL, &wt))
1303      {
1304        case 0: /* not ready */ return "not ready";
1305        case -1: /*error*/      return "error";
1306        case 1: /*ready ? */    break;
1307      }
1308    /* yes: read 1 char*/
1309    /* if \n, check again with select else ungetc(c), ready*/
1310      int c=fgetc(d->f_read);
1311      //Print("try c=%d\n",c);
1312      if (c== -1) return "eof"; /* eof or error */
1313      else if (isdigit(c))
1314      { ungetc(c,d->f_read); d->ungetc_buf='\1'; return "ready"; }
1315      else if (c>' ')
1316      {
1317        Werror("unknown char in ssiLink(%d)",c);
1318        return "error";
1319      }
1320      /* else: next char */
1321    }
1322  }
1323  else if (strcmp(request, "read") == 0)
1324  {
1325    if (SI_LINK_R_OPEN_P(l) && (!feof(d->f_read))) return "ready";
1326    else return "not ready";
1327  }
1328  else if (strcmp(request, "write") == 0)
1329  {
1330    if (SI_LINK_W_OPEN_P(l)) return "ready";
1331    else return "not ready";
1332  }
1333  else return "unknown status request";
1334}
1335
1336int slStatusSsiL(lists L, int timeout)
1337{
1338// input: L: a list with links of type
1339//           ssi-fork, ssi-tcp, MPtcp-fork or MPtcp-launch.
1340//           Note: Not every entry in L must be set.
1341//        timeout: timeout for select in micro-seconds
1342//           or -1 for infinity
1343//           or 0 for polling
1344// returns: ERROR (via Werror): L has wrong elements or link not open
1345//           -2: select returns an error
1346//           -1: the read state of all links is eof
1347//           0:  timeout (or polling): none ready,
1348//           i>0: (at least) L[i] is ready
1349  si_link l;
1350  ssiInfo *d;
1351  #ifdef HAVE_MPSR
1352  MP_Link_pt dd;
1353  #endif
1354  int d_fd;
1355  fd_set  mask, fdmask;
1356  FD_ZERO(&fdmask);
1357  FD_ZERO(&mask);
1358  int max_fd=0; /* 1 + max fd in fd_set */
1359
1360  /* timeout */
1361  #ifdef HAVE_PSELECT
1362  struct timespec wt;
1363  struct timespec *wt_ptr=&wt;
1364  #else
1365  struct timeval wt;
1366  struct timeval *wt_ptr=&wt;
1367  #endif
1368  int startingtime = getRTimer()/TIMER_RESOLUTION;  // in seconds
1369  if (timeout== -1)
1370  {
1371    wt_ptr=NULL;
1372  }
1373  else
1374  {
1375    wt.tv_sec  = timeout / 1000000;
1376  #ifdef HAVE_PSELECT
1377    wt.tv_nsec = 1000 * (timeout % 1000000);
1378  #else
1379    wt.tv_usec = timeout % 1000000;
1380  #endif
1381  }
1382
1383  /* signal mask for pselect() */
1384  sigset_t sigmask;
1385  if(sigprocmask(SIG_SETMASK, NULL, &sigmask) < 0)
1386  {
1387    WerrorS("error in sigprocmask()");
1388    return -2;
1389  }
1390  if(sigaddset(&sigmask, SIGCHLD) < 0)
1391  {
1392    WerrorS("error in sigaddset()");
1393    return -2;
1394  }
1395
1396  /* auxiliary variables */
1397  int i;
1398  int j;
1399  int k;
1400  int s;
1401  char fdmaskempty;
1402
1403  /* check the links and fill in fdmask */
1404  for(i=L->nr; i>=0; i--)
1405  {
1406    if (L->m[i].Typ()!=DEF_CMD)
1407    {
1408      if (L->m[i].Typ()!=LINK_CMD)
1409      { WerrorS("all elements must be of type link"); return -2;}
1410      l=(si_link)L->m[i].Data();
1411      if(SI_LINK_OPEN_P(l)==0)
1412      { WerrorS("all links must be open"); return -2;}
1413      if (((strcmp(l->m->type,"ssi")!=0) && (strcmp(l->m->type,"MPtcp")!=0))
1414      || ((strcmp(l->mode,"fork")!=0) && (strcmp(l->mode,"tcp")!=0)
1415        && (strcmp(l->mode,"launch")!=0)))
1416      {
1417        WerrorS("all links must be of type ssi:fork, ssi:tcp, MPtcp:fork\n");
1418        WerrorS("or MPtcp:launch");
1419        return -2;
1420      }
1421    #ifdef HAVE_MPSR
1422      if (strcmp(l->m->type,"ssi")==0)
1423      {
1424        d=(ssiInfo*)l->data;
1425        d_fd=d->fd_read;
1426      }
1427      else
1428      {
1429        dd=(MP_Link_pt)l->data;
1430        d_fd=((MP_TCP_t *)dd->transp.private1)->sock;
1431      }
1432    #else
1433      d=(ssiInfo*)l->data;
1434      d_fd=d->fd_read;
1435    #endif
1436      FD_SET(d_fd, &fdmask);
1437      if (d_fd > max_fd) max_fd=d_fd;
1438    }
1439  }
1440  max_fd++;
1441
1442do_select:
1443  /* copy fdmask to mask */
1444  FD_ZERO(&mask);
1445  for(k = 0; k < max_fd; k++)
1446  {
1447    if(FD_ISSET(k, &fdmask))
1448    {
1449      FD_SET(k, &mask);
1450    }
1451  }
1452
1453  /* check with select: chars waiting: no -> not ready */
1454  #ifdef HAVE_PSELECT
1455  s = pselect(max_fd, &mask, NULL, NULL, wt_ptr, &sigmask);
1456  #else
1457  SSI_BLOCK_CHLD;
1458  s = select(max_fd, &mask, NULL, NULL, wt_ptr);
1459  SSI_UNBLOCK_CHLD;
1460  #endif
1461
1462  if (s==-1)
1463  {
1464    WerrorS("error in select call");
1465    return -2; /*error*/
1466  }
1467  if (s==0)
1468  {
1469    return 0; /*poll: not ready */
1470  }
1471  else /* s>0, at least one ready  (the number of fd which are ready is s)*/
1472  {
1473    j=0;
1474    while (j<=max_fd) { if (FD_ISSET(j,&mask)) break; j++; }
1475    for(i=L->nr; i>=0; i--)
1476    {
1477      if (L->m[i].rtyp==LINK_CMD)
1478      {
1479        l=(si_link)L->m[i].Data();
1480        #ifdef HAVE_MPSR
1481        if (strcmp(l->m->type,"ssi")!=0)
1482        {
1483          // for MP links, return here:
1484          dd=(MP_Link_pt)l->data;
1485          d_fd=((MP_TCP_t *)dd->transp.private1)->sock;
1486          if(j==d_fd) return i+1;
1487        }
1488        else
1489        {
1490          d=(ssiInfo*)l->data;
1491          d_fd=d->fd_read;
1492          if(j==d_fd) break;
1493        }
1494        #else
1495        d=(ssiInfo*)l->data;
1496        d_fd=d->fd_read;
1497        if(j==d_fd) break;
1498        #endif
1499      }
1500    }
1501    // only ssi links:
1502    if (d->ungetc_buf) return i+1;
1503    loop
1504    {
1505      /* yes: read 1 char*/
1506      /* if \n, check again with select else ungetc(c), ready*/
1507      /* setting: d: current ssiInfo, j current fd, i current entry in L*/
1508      int c=fgetc(d->f_read);
1509      //Print("try c=%d\n",c);
1510      if (c== -1) /* eof */
1511      {
1512        FD_CLR(j,&fdmask);
1513        fdmaskempty = 1;
1514        for(k = 0; k < max_fd; k++)
1515        {
1516          if(FD_ISSET(k, &fdmask))
1517          {
1518            fdmaskempty = 0;
1519            break;
1520          }
1521        }
1522        if(fdmaskempty)
1523        {
1524          return -1;
1525        }
1526        if(timeout != -1)
1527        {
1528          timeout = si_max(0,
1529               timeout - 1000000*(getRTimer()/TIMER_RESOLUTION - startingtime));
1530          wt.tv_sec  = timeout / 1000000;
1531          #ifdef HAVE_PSELECT
1532          wt.tv_nsec = 1000 * (timeout % 1000000);
1533          #else
1534          wt.tv_usec = (timeout % 1000000);
1535          #endif
1536        }
1537        goto do_select;
1538      }
1539
1540      else if (isdigit(c))
1541      { ungetc(c,d->f_read); d->ungetc_buf='\1'; return i+1; }
1542      else if (c>' ')
1543      {
1544        Werror("unknown char in ssiLink(%d)",c);
1545        return -2;
1546      }
1547      /* else: next char */
1548    }
1549  }
1550}
1551
1552int ssiBatch(const char *host, const char * port)
1553/* return 0 on success, >0 else*/
1554{
1555  si_link l=(si_link) omAlloc0Bin(sip_link_bin);
1556  char *buf=(char*)omAlloc(256);
1557  sprintf(buf,"ssi:connect %s:%s",host,port);
1558  slInit(l, buf);
1559  if (slOpen(l,SI_LINK_OPEN,NULL)) return 1;
1560  SI_LINK_SET_RW_OPEN_P(l);
1561
1562  idhdl id = enterid(omStrDup("link_ll"), 0, LINK_CMD, &IDROOT, FALSE);
1563  IDLINK(id) = l;
1564
1565  loop
1566  {
1567    leftv h=ssiRead1(l); /*contains an exit.... */
1568    if (feErrors != NULL && *feErrors != '\0')
1569    {
1570      // handle errors:
1571      PrintS(feErrors); /* currently quite simple */
1572      *feErrors = '\0';
1573    }
1574    ssiWrite(l,h);
1575    h->CleanUp();
1576    omFreeBin(h, sleftv_bin);
1577  }
1578  /* never reached*/
1579  exit(0);
1580}
1581
1582// ----------------------------------------------------------------
1583// format
1584// 1 int %d
1585// 2 string <len> %s
1586// 3 number
1587// 4 bigint 4 %d or 3 <mpz_t>
1588// 5 ring
1589// 6 poly
1590// 7 ideal
1591// 8 matrix
1592// 9 vector
1593// 10 module
1594// 11 command
1595// 12 def <len> %s
1596// 13 proc <len> %s
1597// 14 list %d <elem1> ....
1598// 15 setring .......
1599// 16 nothing
1600// 17 intvec <len> ...
1601//
1602// 20 blackbox <name> ...
1603//
1604// 98: verify version: <ssi-version> <MAX_TOK> <OPT1> <OPT2>
1605// 99: quit Singular
Note: See TracBrowser for help on using the repository browser.